Details
-
Type:
Story
-
Status: Done
-
Resolution: Done
-
Fix Version/s: None
-
Component/s: jointcal
-
Labels:
-
Story Points:8
-
Epic Link:
-
Sprint:AP F18-5
-
Team:Alert Production
Description
import lsst.daf.persistence as dafPersist
|
tract = 9813
|
visit = 23692
|
ccd = 62
|
root = '/datasets/hsc/repo/rerun/RC/w_2018_41/DM-16011-jointcal'
|
butler = dafPersist.Butler(root)
|
photoCalib = butler.get('jointcal_photoCalib', visit=visit, ccd=ccd, tract=tract)
|
>>> photoCalib.getInstFluxMag0()
|
-278609934862735.6
|
>>> butler.get('jointcal_photoCalib_filename', visit=visit, ccd=ccd, tract=tract.getId())
|
['/datasets/hsc/repo/rerun/RC/w_2018_41/DM-16011-jointcal/jointcal-results/HSC-R/9813/jointcal_photoCalib-0023692-062.fits']
|
>>> str(photoCalib)
|
'TransformBoundedField on Box2I(Point2I(0, 0), Extent2I(2048, 4176)) with mean: -3.59376e-15 err: 2.98536e-17'
|
>>> photoCalib.instFluxToMaggies(100)
|
-3.589247456278528e-13
|
>>> photoCalib.instFluxToMaggies(100, lsst.geom.Point2D(0,0))
|
3.9565953264617065e-12
|
Krzysztof Findeisen: do you mind taking on this review? It's smaller than I'd initially said, because the proper fix is being pushed off to
DM-16305, and this solution is just a stop-gap (but with more robust tests to watch for the triggering condition). It's ~60 lines of afw and ~180 lines of jointcal, so not so big.If jira doesn't show it, here's the afw PR: https://github.com/lsst/afw/pull/409