295 | | * Run Joe's script [CalculateHPFfiberpositions.py] with the argument {{{fromscratch}}} and the data file you just edited the ACQ X,Y coordinates of cCFB, sCFB, BIB, LRS2-R, LRS2-B, as shown: |
| 295 | * Run Joe's script {{{CalculateHPFfiberpositions.py}}} (attached to this page, you can download it) with the argument {{{fromscratch}}} and the data file you just edited the ACQ X,Y coordinates of cCFB, sCFB, BIB, LRS2-R, LRS2-B, as shown: |
| 331 | |
| 332 | * to do the analysis, download or copy the Goldilocks reductions of the data |
| 333 | * measure the mean flux in the region [1002:1051,1017:1024] on the 0th extension for all 14 exposures (Order 19, as usual for HPF SNR) |
| 334 | * for each exposure, get the tracker coordinates from header keywords RTFTHETA, RTFPHI, RTFRHO and use hetillum to calculate the illumination at the HPF fiber position for each exposure: |
| 335 | {{{ |
| 336 | hetillum '[ 0.987817, 1.052546, 2.660203]' '[-0.008194,0.0042778]' 256 |
| 337 | ^theta ^phi ^rho ^,^ x,y of HPF sci fiber |
| 338 | }}} |
| 339 | * also get the average guider transparency and IQ from the RA NR for each of the 14 exposures. If there were big variations or changes you can try to use those measurements to compensate but it may be safer to redo the test under more stable conditions. |
| 340 | * consider using the supermongo codes (or perhaps re-write this in python) attached at the bottom of this page to calculate the x,y offset of the science fiber. In principle the first and last positions should be the brightest (after correcting for the changing pupil illumination), and the others should fall off symmetrically around them. I made some supermongo code that assumed a PSF and predicted the flux in each of the 14 positions. It allowed the user to guess an x,y offset and it would give you something like a chi-squared value to indicate the goodness of fit. I tweaked the x,y inputs until I found a reasonable minimum. This approach is ripe for improvement! Below is an example graph after doing some tests to find the best x,y offset: |
| 341 | |
| 342 | [[Image(hpfcentering.png, 300px)]] |
| 343 | |