[PIPE2D-807] pfsObject variance is much too large

Mon, 19 Apr 2021 18:45:16 +0000

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()

PFS-JIRA

[PIPE2D-807] pfsObject variance is much too large