Add GaaP 1p0 flux to skyObject[Source] histograms

This also adds vertical solid black lines at the "expected" value (0.0 here) and a shaded PDF(0, 1) distribution on the now density normalized S/N plots for reference to both skyObject and skySource histogram plots.
lsst · Dec 8, 2022 · 232701a · 232701a
1 parent 1e39efc
commit 232701a
Show file tree

Hide file tree

Showing 2 changed files with 31 additions and 8 deletions.
diff --git a/python/lsst/analysis/tools/analysisPlots/skyObject.py b/python/lsst/analysis/tools/analysisPlots/skyObject.py
@@ -61,16 +61,27 @@ def setDefaults(self):
         self.prep.selectors.skyObjectSelector = SkyObjectSelector()
 
         self.process.buildActions.hist_psf_flux = LoadVector(vectorKey="{band}_psfFlux")
-        self.process.buildActions.hist_09_flux = LoadVector(vectorKey="{band}_ap09Flux")
+        self.process.buildActions.hist_ap09_flux = LoadVector(vectorKey="{band}_ap09Flux")
+        self.process.buildActions.hist_gaap1p0_flux = LoadVector(vectorKey="{band}_gaap1p0Flux")
         self.process.buildActions.hist_psf_sn = SNCalculator(fluxType="{band}_psfFlux")
-        self.process.buildActions.hist_09_sn = SNCalculator(fluxType="{band}_ap09Flux")
+        self.process.buildActions.hist_ap09_sn = SNCalculator(fluxType="{band}_ap09Flux")
+        self.process.buildActions.hist_gaap1p0_sn = SNCalculator(fluxType="{band}_gaap1p0Flux")
 
         self.produce = HistPlot()
 
         self.produce.panels["panel_flux"] = HistPanel()
         self.produce.panels["panel_flux"].label = "Flux (nJy)"
-        self.produce.panels["panel_flux"].hists = dict(hist_psf_flux="psfFlux", hist_09_flux="ap09Flux")
+        self.produce.panels["panel_flux"].hists = dict(
+            hist_psf_flux="psfFlux", hist_ap09_flux="ap09Flux", hist_gaap1p0_flux="gaap1p0Flux"
+        )
+        self.produce.panels["panel_flux"].expectedValue = 0.0
+        self.produce.panels["panel_flux"].validate()
 
         self.produce.panels["panel_sn"] = HistPanel()
         self.produce.panels["panel_sn"].label = "S/N"
-        self.produce.panels["panel_sn"].hists = dict(hist_psf_sn="psfFlux S/N", hist_09_sn="ap09Flux S/N")
+        self.produce.panels["panel_sn"].hists = dict(
+            hist_psf_sn="psf S/N", hist_ap09_sn="ap09 S/N", hist_gaap1p0_sn="gaap1p0 S/N"
+        )
+        self.produce.panels["panel_sn"].expectedValue = 0.0
+        self.produce.panels["panel_sn"].histDensity = True
+        self.produce.panels["panel_sn"].validate()
diff --git a/python/lsst/analysis/tools/analysisPlots/skySource.py b/python/lsst/analysis/tools/analysisPlots/skySource.py
@@ -65,16 +65,28 @@ def setDefaults(self):
         self.prep.selectors.skySourceSelector = SkySourceSelector()
 
         self.process.buildActions.hist_psf_flux = LoadVector(vectorKey="psfFlux")
-        self.process.buildActions.hist_09_flux = LoadVector(vectorKey="ap09Flux")
+        self.process.buildActions.hist_ap09_flux = LoadVector(vectorKey="ap09Flux")
+        self.process.buildActions.hist_gaap1p0_flux = LoadVector(vectorKey="gaap1p0Flux")
         self.process.buildActions.hist_psf_sn = SNCalculator(fluxType="psfFlux")
-        self.process.buildActions.hist_09_sn = SNCalculator(fluxType="ap09Flux")
+        self.process.buildActions.hist_ap09_sn = SNCalculator(fluxType="ap09Flux")
+        self.process.buildActions.hist_gaap1p0_sn = SNCalculator(vectorKey="gaap1p0Flux")
 
         self.produce = HistPlot()
 
         self.produce.panels["panel_flux"] = HistPanel()
         self.produce.panels["panel_flux"].label = "Flux (nJy)"
-        self.produce.panels["panel_flux"].hists = dict(hist_psf_flux="psfFlux", hist_09_flux="ap09Flux")
+        self.produce.panels["panel_flux"].hists = dict(
+            hist_psf_flux="psfFlux", hist_ap09_flux="ap09Flux", hist_gaap1p0_flux="gaap1p0Flux"
+        )
+
+        self.produce.panels["panel_flux"].expectedValue = 0.0
+        self.produce.panels["panel_flux"].validate()
 
         self.produce.panels["panel_sn"] = HistPanel()
         self.produce.panels["panel_sn"].label = "S/N"
-        self.produce.panels["panel_sn"].hists = dict(hist_psf_sn="psfFlux S/N", hist_09_sn="ap09Flux S/N")
+        self.produce.panels["panel_sn"].hists = dict(
+            hist_psf_sn="psf S/N", hist_ap09_sn="ap09 S/N", hist_gaap1p0_sn="gaap1p0 S/N"
+        )
+        self.produce.panels["panel_sn"].expectedValue = 0.0
+        self.produce.panels["panel_sn"].histDensity = True
+        self.produce.panels["panel_sn"].validate()