The data used in this KPM was the i-band single visit data from the public HSC engineering run.

To create the plot I took the stars from each visit that were selected for the psf modeling and computed the difference in ellipticities between the model of the psf and the sdss moments. I then computed the average dot product of the residuals for all pairs of stars at different separation distances.

Once this was computed for each visit I then took the median from all of the i-band visits to create the points in the plot.

Reading off the plot, the residual correlations are:

• less than 1 arcmin: 6e-5 requirement: 500e-5
• beween 1-5 arcmin: 2e-5 (see DM-3047) requirement: 500e-5

Both of these are well below the requiremnt

I should add that to find all the pairs and compute the averages I used code written by Mike Jarvis called treecorr. It is designed to quickly calculate correlation functions. It was not strictly necessary to use this software for the HSC data in this KPM because the number of stars in a single visit is relatively small. However, it will be necessary when trying to do calculations on larger datasets. Stile also uses treecorr to calculate the correlation functions.

We've now shipped the S15 release including the KPM measurement required by this epic. Is any further work required before we close it?

Thanks for the review request!

Comments/questions:

• I think it'd be good to add a reference to http://ls.st/srd (the "Galaxy shear measurement accuracy, and PSF ellipticity residuals") section to make it easy to trace what's being measured
• It will be easier if the measurement code is in a git repository
• It would be better from my PoV to put all TE* measurements into one JIRA Issue (or even a "KPM testing definition document") – having them scattered over four is a bit of a pain to review and track. I'll understand if there are management/procedure issues that make this difficult.
• From the code: I wonder if this is correct?

  useFlag = src.schema.find('calib_psfUsed').key

  mask = src.get(useFlag)==1

  Does this snippet select only stars that were used to derive the PSF? If yes: I'd argue we should actually measure the residuals only on stars not used to derive the PSF model, as the goal is to establish that the interpolated PSF residuals are small anywhere in the focal plane (i.e., at locations of galaxies whose shapes we'll measure). For example, if I was evil I could devise an "interpolation scheme" where PSF at the position of each star is exactly the image of the star, while elsewhere it's zero – that would pass your test. If no: ignore.

Noticing that this is an FY15 task, if it's more consistent with the work scheduling methodology, it'd be OK with me to resolve the above as a part of FY16 KPM measurement work. Also, similar to DM-3057 another pair of eyes (David Nidever [X]'s) on the code would be good to have.