https://pfspipe.ipmu.jp/jira This file is an XML representation of an issue en-us 8.3.4 803005 13-09-2019 https://pfspipe.ipmu.jp/jira/browse/PIPE2D-616 DRP 2-D Pipeline <p>With respect to determining the alignment for medium resolution exposures, as discussed today during the SpS hardware telecon, please compute using the 2D DRP second moments for the following medium-resolution HgAr visit sets:  </p> <p>3443...3450<br/> 18052...18059</p> <p>and the following lines:</p> <p>696.7352, 727.494 , 795.0362, 912.5471, 966.0435</p> <p>Note that as the above data are based on the HgAr lamp, the continuum needs to be removed.</p> <p>These data would be used to compare with the plots generated by <a href="https://pfspipe.ipmu.jp/jira/secure/ViewProfile.jspa?name=fmadec" class="user-hover" rel="fmadec">fmadec</a> using the SEP python module (<a href="https://sep.readthedocs.io/en/v1.0.x/" class="external-link" rel="nofollow">https://sep.readthedocs.io/en/v1.0.x/</a>) in slide 12 of the following presentation: <a href="https://sumire.pbworks.com/w/file/fetch/140734437/LAM_Subaru_SM1_20200715.pdf" class="external-link" rel="nofollow">https://sumire.pbworks.com/w/file/fetch/140734437/LAM_Subaru_SM1_20200715.pdf</a></p> <p>So a similar set of plots should be generated for this ticket.</p> PIPE2D-616

Compute 2nd moments for Med Res data

Task Normal Done Done hassan hassan Wed, 15 Jul 2020 17:52:15 +0000 Mon, 4 Jan 2021 20:23:19 +0000 Fri, 4 Sep 2020 21:44:33 +0000 0 5 <p>Please start with processing the post-migration visit_set (18052...18059) which should currently be available on tigress.</p> <p>We need to be quite careful here. The background (including continuum) subtraction has to be done carefully, and a non-trivial moment calculation used (HSM would be my first choice)</p> <p>Is it not necessary to use SExtractor in order to compare apples-to-apples with the plots in the presentation?</p> <p>Please let me know if this is still desired, or if it is moot following the acquisition of new data.</p> <p>Yes, this is still desired. The data haven't yet been acquired from SM1. I'm hoping by early next week.</p> <p>Data taken. The command sequence was: </p> <div class="preformatted panel" style="border-width: 1px;"><div class="preformattedContent panelContent"> <pre>start at low res take 1xflat + 3xarc for rep in 1..5: move to med res take 1xflat + 3xarc move to low res take 1xflat + 3xarc </pre> </div></div> <p>And the arcs are neon=50 argon=50 krypton=50</p> <p>The acquisition log is:</p> <div class="preformatted panel" style="border-width: 1px;"><div class="preformattedContent panelContent"> <pre>set_id|sequence_type|name|start_date |min_visit|max_visit|cmd_str | ------|-------------|----|-------------------|---------|---------|--------------------------------------------------------------------------------| 1348|flats | |2020-08-17 17:02:30| 23743| 23743|expose flat halogen=3 comments=test_gratingMotion_0 | 1349|arcs | |2020-08-17 17:03:34| 23744| 23746|expose arc neon=50 argon=50 krypton=50 duplicate=3 comments=test_gratingMotion_0| 1350|flats | |2020-08-17 17:11:45| 23747| 23747|expose flat halogen=3 comments=test_gratingMotion_1 | 1351|arcs | |2020-08-17 17:12:49| 23748| 23750|expose arc neon=50 argon=50 krypton=50 duplicate=3 comments=test_gratingMotion_1| 1352|flats | |2020-08-17 17:20:58| 23751| 23751|expose flat halogen=3 comments=test_gratingMotion_1 | 1353|arcs | |2020-08-17 17:22:03| 23752| 23754|expose arc neon=50 argon=50 krypton=50 duplicate=3 comments=test_gratingMotion_1| 1354|flats | |2020-08-17 17:30:12| 23755| 23755|expose flat halogen=3 comments=test_gratingMotion_2 | 1355|arcs | |2020-08-17 17:31:17| 23756| 23758|expose arc neon=50 argon=50 krypton=50 duplicate=3 comments=test_gratingMotion_2| 1356|flats | |2020-08-17 17:39:26| 23759| 23759|expose flat halogen=3 comments=test_gratingMotion_2 | 1357|arcs | |2020-08-17 17:40:30| 23760| 23762|expose arc neon=50 argon=50 krypton=50 duplicate=3 comments=test_gratingMotion_2| 1358|flats | |2020-08-17 17:48:39| 23763| 23763|expose flat halogen=3 comments=test_gratingMotion_3 | 1359|arcs | |2020-08-17 17:49:43| 23764| 23766|expose arc neon=50 argon=50 krypton=50 duplicate=3 comments=test_gratingMotion_3| 1360|flats | |2020-08-17 17:57:53| 23767| 23767|expose flat halogen=3 comments=test_gratingMotion_3 | 1361|arcs | |2020-08-17 17:58:57| 23768| 23770|expose arc neon=50 argon=50 krypton=50 duplicate=3 comments=test_gratingMotion_3| 1362|flats | |2020-08-17 18:07:06| 23771| 23771|expose flat halogen=3 comments=test_gratingMotion_4 | 1363|arcs | |2020-08-17 18:08:11| 23772| 23774|expose arc neon=50 argon=50 krypton=50 duplicate=3 comments=test_gratingMotion_4| 1364|flats | |2020-08-17 18:16:19| 23775| 23775|expose flat halogen=3 comments=test_gratingMotion_4 | 1365|arcs | |2020-08-17 18:17:23| 23776| 23778|expose arc neon=50 argon=50 krypton=50 duplicate=3 comments=test_gratingMotion_4| 1366|flats | |2020-08-17 18:25:33| 23779| 23779|expose flat halogen=3 comments=test_gratingMotion_5 | 1367|arcs | |2020-08-17 18:26:37| 23780| 23782|expose arc neon=50 argon=50 krypton=50 duplicate=3 comments=test_gratingMotion_5| 1368|flats | |2020-08-17 18:34:49| 23783| 23783|expose flat halogen=3 comments=test_gratingMotion_5 | 1369|arcs | |2020-08-17 18:35:53| 23784| 23786|expose arc neon=50 argon=50 krypton=50 duplicate=3 comments=test_gratingMotion_5| </pre> </div></div> <p>The measurements are on the tiger cluster in <tt>/tigress/pprice/pipe2d-616</tt>, in the same format as used for measurements by the LSST pipeline (FITS tables, first HDU is of the most interest).</p> <p>The below plot shows the measured shape of all lines as a function of x,y on the detector. The user will probably want to threshold the list of lines by the flux, multiple measures of which are in the catalogs.</p> <p><span class="image-wrap" style=""><img src="https://pfspipe.ipmu.jp/jira/secure/attachment/12837/12837_shape.png" style="border: 0px solid black" /></span></p> <div class="code panel" style="border-width: 1px;"><div class="codeContent panelContent"> <pre class="code-java"> sqlite&gt; SELECT visit, GROUP_CONCAT(arm), pfsDesignId, field, lamps, exptime, attenuator, dither, shift, focus FROM raw WHERE visit BETWEEN 23743 AND 23782 GROUP BY visit; 23743|r,b|1099528409104|FLAT|Quartz|3.0|-9998.0|0.0|0.0|-0.0003 23744|r,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23745|r,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23746|r,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23747|b,m|1099528409104|FLAT|Quartz|3.0|-9998.0|0.0|0.0|-0.0003 23748|m,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23749|m,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23750|b,m|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23751|b,r|1099528409104|FLAT|Quartz|3.0|-9998.0|0.0|0.0|-0.0003 23752|b,r|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23753|b,r|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23754|b,r|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23755|m,b|1099528409104|FLAT|Quartz|3.0|-9998.0|0.0|0.0|-0.0003 23756|m,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23757|m,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23758|m,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23759|b,r|1099528409104|FLAT|Quartz|3.0|-9998.0|0.0|0.0|-0.0003 23760|b,r|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23761|b,r|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23762|b,r|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23763|b,m|1099528409104|FLAT|Quartz|3.0|-9998.0|0.0|0.0|-0.0003 23764|m,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23765|m,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23766|m,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23767|b,r|1099528409104|FLAT|Quartz|3.0|-9998.0|0.0|0.0|-0.0003 23768|r,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23769|r,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23770|r,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23771|m,b|1099528409104|FLAT|Quartz|3.0|-9998.0|0.0|0.0|-0.0003 23772|b,m|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23773|m,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23774|b,m|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23775|b,r|1099528409104|FLAT|Quartz|3.0|-9998.0|0.0|0.0|-0.0003 23776|b,r|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23777|b,r|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23778|r,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23779|m,b|1099528409104|FLAT|Quartz|3.0|-9998.0|0.0|0.0|-0.0003 23780|m,b|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23781|b,m|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 23782|b,m|1099528409104|COMPARISON|Ar,Kr,Ne|50.0|-9998.0|0.0|0.0|-0.0003 (lsst-scipipe) pprice@tiger2-sumire:/projects/HSC/PFS/Subaru $ ingestPfsImages.py /projects/HSC/PFS/Subaru --pfsConfigDir /projects/HSC/PFS/Subaru/drp/pfsDesign --config parse.pfsDesignId=1099528409104 --ignore-ingested -- <span class="code-quote">'/projects/HSC/PFS/Subaru/raw/2020-08-*/sps/PFSA*.fits'</span> The last r/m biases were 6 weeks ago (2020-06-30), probably taken as part of the stability test (PIPE2D-618). Will push forward assuming those are good enough. Steps: 1. Bootstrap a valid detectorMap. 2. Generate fiberTrace. 3. <span class="code-object">Process</span> the arc to get a really good detectorMap. 4. Use custom script to run LSST measurement framework on arc lines. (lsst-scipipe) pprice@tiger2-sumire:/tigress/pprice/pipe2d-616 $ ingestPfsCalibs.py /projects/HSC/PFS/Subaru --output=/projects/HSC/PFS/Subaru/CALIB-price --validity=1800 --doraise --mode=copy ~/pfs/drp_pfs_data/detectorMap/detectorMap-2020Apr-m1.fits (lsst-scipipe) pprice@tiger2-sumire:/tigress/pprice/pipe2d-616 $ bootstrapDetectorMap.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/bootstrap --doraise --flatId visit=23747 arm=m --arcId visit=23748 arm=m bootstrap INFO: Median difference from detectorMap: -1.268960,0.469803 pixels bootstrap INFO: Fit 277/336 points, rms: x=0.722616 y=2.441186 total=2.069997 pixels bootstrap INFO: Updating detectorMap... bootstrap INFO: Median difference from detectorMap: 0.698470,0.597844 pixels bootstrap INFO: Fit 191/229 points, rms: x=0.345339 y=2.602193 total=2.158558 pixels (lsst-scipipe) pprice@tiger2-sumire:/tigress/pprice/pipe2d-616 $ ingestPfsCalibs.py /projects/HSC/PFS/Subaru --output=/projects/HSC/PFS/Subaru/CALIB-price --validity=1800 --doraise --mode=copy /projects/HSC/PFS/Subaru/rerun/price/pipe2d-616/bootstrap/DETECTORMAP/pfsDetectorMap-023748-m1.fits (lsst-scipipe) pprice@tiger2-sumire:/tigress/pprice/pipe2d-616 $ constructFiberTrace.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/fiberTrace --doraise --batch-type=smp --cores=1 --id visit=23747 (lsst-scipipe) pprice@tiger2-sumire:/tigress/pprice/pipe2d-616 $ ingestPfsCalibs.py /projects/HSC/PFS/Subaru --output=/projects/HSC/PFS/Subaru/CALIB-price --validity=1800 --doraise --mode=copy /projects/HSC/PFS/Subaru/rerun/price/pipe2d-616/fiberTrace/FIBERTRACE/pfsFiberTrace-2020-08-17-023747-m1.fits (lsst-scipipe) pprice@tiger2-sumire:/tigress/pprice/pipe2d-616 $ reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23748 arm=m reduceArc.fitDetectorMap INFO: Final fit: chi2=1121.654719 xRMS=0.154362 yRMS=1.054392 (0.049562 nm, 18.599770 km/s) from 445/473 lines reduceArc.fitDetectorMap INFO: Fit quality from reserved lines: chi2=258.705408 xRMS=0.185894 yRMS=0.881486 (0.041434 nm, 15.549659 km/s) from 48 lines (9.1%) reduceArc.fitDetectorMap INFO: Softening errors by 0.115597 pixels (0.005434 nm, 2.039167 km/s) to yield chi^2/dof=1 reduceArc.fitDetectorMap INFO: Softened fit: chi2=929.695791 xRMS=0.092336 yRMS=0.802256 (0.037601 nm, 14.109900 km/s) from 445/526 lines reduceArc.fitDetectorMap INFO: Softened fit quality from reserved lines: chi2=105.743892 xRMS=0.108362 yRMS=0.697862 (0.032708 nm, 12.273840 km/s) from 48 lines (9.1%) Having difficulty getting good wavelength fits to the m arm with the <span class="code-keyword">new</span> GlobalDetectorMap (tickets/PIPE2D-358; getting systematic wavelength residuals), so revert to the current version; will work on fitting m later. reduceArc.calibrateWavelengths INFO: FiberId 2, rms 0.003757 nm (0.081 pix) from 54/62 (0.002162 nm = 0.047 pix <span class="code-keyword">for</span> 10 reserved points), 714.90-885.63 nm reduceArc.calibrateWavelengths INFO: Fiber 2: wavelength correction -0.015799 +/- 0.010838 nm reduceArc.calibrateWavelengths INFO: FiberId 63, rms 0.002983 nm (0.064 pix) from 53/63 (0.001688 nm = 0.036 pix <span class="code-keyword">for</span> 10 reserved points), 717.59-885.63 nm reduceArc.calibrateWavelengths INFO: Fiber 63: wavelength correction -0.018139 +/- 0.024752 nm reduceArc.calibrateWavelengths INFO: FiberId 192, rms 0.003791 nm (0.082 pix) from 48/57 (0.004014 nm = 0.087 pix <span class="code-keyword">for</span> 10 reserved points), 714.90-885.63 nm reduceArc.calibrateWavelengths INFO: Fiber 192: wavelength correction -0.017744 +/- 0.044794 nm reduceArc.calibrateWavelengths INFO: FiberId 255, rms 0.002789 nm (0.060 pix) from 50/63 (0.002168 nm = 0.047 pix <span class="code-keyword">for</span> 10 reserved points), 714.90-885.63 nm reduceArc.calibrateWavelengths INFO: Fiber 255: wavelength correction -0.016435 +/- 0.028260 nm reduceArc.calibrateWavelengths INFO: FiberId 339, rms 0.002716 nm (0.058 pix) from 49/61 (0.003449 nm = 0.074 pix <span class="code-keyword">for</span> 10 reserved points), 714.90-885.63 nm reduceArc.calibrateWavelengths INFO: Fiber 339: wavelength correction -0.017634 +/- 0.021243 nm reduceArc.calibrateWavelengths INFO: FiberId 401, rms 0.003812 nm (0.082 pix) from 50/63 (0.002040 nm = 0.044 pix <span class="code-keyword">for</span> 10 reserved points), 714.90-885.63 nm reduceArc.calibrateWavelengths INFO: Fiber 401: wavelength correction -0.011027 +/- 0.027518 nm reduceArc.calibrateWavelengths INFO: FiberId 464, rms 0.004156 nm (0.090 pix) from 56/60 (0.004277 nm = 0.092 pix <span class="code-keyword">for</span> 10 reserved points), 714.90-885.63 nm reduceArc.calibrateWavelengths INFO: Fiber 464: wavelength correction -0.010072 +/- 0.030665 nm reduceArc.calibrateWavelengths INFO: FiberId 525, rms 0.004800 nm (0.103 pix) from 46/55 (0.005414 nm = 0.117 pix <span class="code-keyword">for</span> 10 reserved points), 714.90-885.63 nm reduceArc.calibrateWavelengths INFO: Fiber 525: wavelength correction -0.015462 +/- 0.032322 nm reduceArc.calibrateWavelengths INFO: FiberId 587, rms 0.003302 nm (0.071 pix) from 48/57 (0.003198 nm = 0.069 pix <span class="code-keyword">for</span> 10 reserved points), 714.90-885.63 nm reduceArc.calibrateWavelengths INFO: Fiber 587: wavelength correction -0.010947 +/- 0.030819 nm reduceArc.calibrateWavelengths INFO: FiberId 650, rms 0.003758 nm (0.081 pix) from 49/60 (0.004431 nm = 0.096 pix <span class="code-keyword">for</span> 10 reserved points), 714.90-885.63 nm reduceArc.calibrateWavelengths INFO: Fiber 650: wavelength correction -0.020529 +/- 0.026221 nm (lsst-scipipe) pprice@tiger2-sumire:/tigress/pprice/pipe2d-616 $ constructFiberTrace.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/fiberTrace --doraise --batch-type=smp --cores=1 --id visit=23743 arm=r (lsst-scipipe) pprice@tiger2-sumire:/tigress/pprice/pipe2d-616 $ ingestPfsCalibs.py /projects/HSC/PFS/Subaru --output=/projects/HSC/PFS/Subaru/CALIB-price --validity=1800 --doraise --mode=copy /projects/HSC/PFS/Subaru/rerun/price/pipe2d-616/fiberTrace/FIBERTRACE/pfsFiberTrace-2020-08-17-023743-r1.fits (lsst-scipipe) pprice@tiger2-sumire:/tigress/pprice/pipe2d-616 $ reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23744 arm=r reduceArc.calibrateWavelengths INFO: FiberId 339, rms 0.007783 nm (0.090 pix) from 80/94 (0.005586 nm = 0.065 pix <span class="code-keyword">for</span> 10 reserved points), 640.40-955.00 nm reduceArc.calibrateWavelengths INFO: Fiber 339: wavelength correction 0.045030 +/- 0.006753 nm reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23744..23746 arm=r -j 3 -c reduceExposure.doWriteCalexp=True reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23748..23750 arm=m -j 3 -c reduceExposure.doWriteCalexp=True reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23752..23754 arm=r -j 3 -c reduceExposure.doWriteCalexp=True reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23756..23758 arm=m -j 3 -c reduceExposure.doWriteCalexp=True reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23760..23762 arm=r -j 3 -c reduceExposure.doWriteCalexp=True reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23764..23766 arm=m -j 3 -c reduceExposure.doWriteCalexp=True reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23768..23770 arm=r -j 3 -c reduceExposure.doWriteCalexp=True reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23772..23774 arm=m -j 3 -c reduceExposure.doWriteCalexp=True reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23776..23778 arm=r -j 3 -c reduceExposure.doWriteCalexp=True reduceArc.py /projects/HSC/PFS/Subaru --calib=/projects/HSC/PFS/Subaru/CALIB-price --rerun=price/pipe2d-616/arc --id visit=23780..23782 arm=m -j 3 -c reduceExposure.doWriteCalexp=True (lsst-scipipe) pprice@tiger2-sumire:/tigress/pprice/pipe2d-616 $ python pipe2d-616.py /projects/HSC/PFS/Subaru --calib /projects/HSC/PFS/Subaru/CALIB-price --rerun price/pipe2d-616/arc --id visit=23744..23746^23748..23750^23752..23754^23756..23758^23760..23762^23764..23766^23768..23770^23772..23774^23776..23778^23780..23782 arm=r^m -j 20 </pre> </div></div> <p>Added <tt>ext&#95;simpleShape&#95;SimpleShape</tt>, <tt>ext&#95;shapeHSM&#95;HsmSourceMoments</tt>, and <tt>ext&#95;shapeHSM&#95;HsmSourceMomentsRound</tt> measurements at <a href="https://pfspipe.ipmu.jp/jira/secure/ViewProfile.jspa?name=rhl" class="user-hover" rel="rhl">rhl</a>'s request.</p> <p>This provides a clean sample:</p> <div class="code panel" style="border-width: 1px;"><div class="codeContent panelContent"> <pre class="code-java"> &gt;&gt;&gt; gaussMag = -2.5*np.log10(catalog[<span class="code-quote">"base_GaussianFlux_instFlux"</span>]) &gt;&gt;&gt; psfMag = -2.5*np.log10(catalog[<span class="code-quote">"base_PsfFlux_instFlux"</span>]) &gt;&gt;&gt; saturated = catalog[<span class="code-quote">"base_PixelFlags_flag_saturatedCenter"</span>] &gt;&gt;&gt; select = (psfMag &lt; -9) &amp; (psfMag &gt; -14) &amp; ~saturated &amp; (gaussMag - psfMag &gt; -0.55) </pre> </div></div> <p>I'm a bit concerned that this removes low-wavelength low-<tt>fiberId</tt> lines (bottom right-hand corner). Decreasing the <tt>gaussMag - psfMag</tt> limit gets them back, at the cost of admitting some contamination.</p> <p>I've added a <tt>numNear</tt> column, which is the number of lines within 0.4 nm. I now like the selection:</p> <div class="code panel" style="border-width: 1px;"><div class="codeContent panelContent"> <pre class="code-java"> &gt;&gt;&gt; select = (psfMag &lt; -9) &amp; (psfMag &gt; -14) &amp; ~saturated &amp; (catalog[<span class="code-quote">"numNear"</span>] == 0) </pre> </div></div> <p>Work presented to Fabrice last week. Further work necessary, but need more discussion with Fabrice, and new tickets will be filed.</p> Development Rank 0|02qps2:ji Reviewers hassan Sprint 2DDRP-2021 A Story Points 3.0