-
Type:
Story
-
Status: Done (View Workflow)
-
Priority:
Normal
-
Resolution: Done
-
Affects Version/s: None
-
Fix Version/s: None
-
Component/s: None
-
Labels:None
-
Sprint:2DDRP-2021 A3
The following code with
dataId = {'visit': 46032, 'arm': 'r', 'spectrograph': 1}
returns an empty image; in fact, continuum contains no fibres.
from pfs.drp.stella import SpectrumSet exp = butler.get("calexp", dataId) detMap = butler.get("detectorMap", dataId) spectra = butler.get("pfsArm", dataId) profiles = butler.get("fiberProfiles", dataId) traces = profiles.makeFiberTracesFromDetectorMap(detMap) continua = fitContinuumTask.run(SpectrumSet.fromPfsArm(spectra)) im = continua.makeImage(exp.getBBox(), traces)