-
Type:
Story
-
Status: Done (View Workflow)
-
Priority:
Normal
-
Resolution: Done
-
Affects Version/s: None
-
Fix Version/s: None
-
Component/s: None
-
Labels:None
-
Story Points:1
-
Sprint:2DDRP-2021 A 4
Kazumi Murata reports that the variance for a pfsObject from the integration test is much too large compared to the flux.
>>> import numpy as np >>> import matplotlib.pyplot as plt >>> from lsst.daf.persistence import Butler >>> from pfs.datamodel.utils import calculatePfsVisitHash >>> butler = Butler("INTEGRATION/rerun/integration/pipeline") >>> visits = [24, 25] >>> spectrum = butler.get("pfsObject", catId=1, tract=0, patch="0,0", objId=55, nVisit=len(visits), pfsVisitHash=calculatePfsVisitHash(visits)) >>> np.median(spectrum.flux/np.sqrt(spectrum.variance)) 1.2713986241730218e-07 >>> plt.plot(spectrum.wavelength, spectrum.flux, "k-") >>> plt.plot(spectrum.wavelength, np.sqrt(spectrum.variance), "r-") >>> plt.plot(spectrum.wavelength, spectrum.sky, "b-") >>> plt.show()