Update sky plots to show the patch outlines as the default

lsst · Jul 26, 2023 · 4068c48 · 4068c48
1 parent 5022aff
commit 4068c48
Show file tree

Hide file tree

Showing 7 changed files with 18 additions and 4 deletions.
diff --git a/python/lsst/analysis/tools/actions/plot/skyPlot.py b/python/lsst/analysis/tools/actions/plot/skyPlot.py
@@ -34,7 +34,7 @@
 
 from ...interfaces import KeyedData, KeyedDataSchema, PlotAction, Scalar, Vector
 from ...statistics import nansigmaMad
-from .plotUtils import addPlotInfo, mkColormap, plotProjectionWithBinning, sortAllArrays
+from .plotUtils import addPlotInfo, generateSummaryStats, mkColormap, plotProjectionWithBinning, sortAllArrays
 
 
 class SkyPlot(PlotAction):
@@ -92,6 +92,9 @@ def getInputSchema(self, **kwargs) -> KeyedDataSchema:
             base.append(("z", Vector))
             base.append(("statMask", Vector))
 
+        if self.plotOutlines:
+            base.append(("patch", Vector))
+
         return base
 
     def __call__(self, data: KeyedData, **kwargs) -> Mapping[str, Figure] | Figure:
@@ -204,8 +207,8 @@ def makePlot(
         fig = plt.figure(dpi=300)
         ax = fig.add_subplot(111)
 
-        if sumStats is None:
-            sumStats = {}
+        if self.plotOutlines:
+            sumStats = generateSummaryStats(data, skymap, plotInfo)
 
         if plotInfo is None:
             plotInfo = {}

diff --git a/python/lsst/analysis/tools/atools/astrometricRepeatability.py b/python/lsst/analysis/tools/atools/astrometricRepeatability.py
@@ -109,6 +109,7 @@ def setDefaults(self):
         self.process.buildActions.residual.buildAction.buildAction = RAcosDec()
         self.process.buildActions.x = LoadVector(vectorKey="coord_ra")
         self.process.buildActions.y = LoadVector(vectorKey="coord_dec")
+        self.process.buildActions.patch = LoadVector(vectorKey="patch")
 
         self.produce = SkyPlot()
 

diff --git a/python/lsst/analysis/tools/atools/astrometryWithReference.py b/python/lsst/analysis/tools/atools/astrometryWithReference.py
@@ -165,6 +165,8 @@ def setDefaults(self, coordinate):
         self.process.buildActions.xStars.vectorKey = "coord_ra_target"
         self.process.buildActions.yStars = LoadVector()
         self.process.buildActions.yStars.vectorKey = "coord_dec_target"
+        self.process.buildActions.patch = LoadVector()
+        self.process.buildActions.patch.vectorKey = "patch"
 
         self.produce = SkyPlot()
         self.produce.plotTypes = ["stars"]

diff --git a/python/lsst/analysis/tools/atools/magDiff.py b/python/lsst/analysis/tools/atools/magDiff.py
@@ -52,6 +52,8 @@ def setDefaults(self):
         self.process.buildActions.xGalaxies.vectorKey = "coord_ra"
         self.process.buildActions.yGalaxies = LoadVector()
         self.process.buildActions.yGalaxies.vectorKey = "coord_dec"
+        self.process.buildActions.patch = LoadVector()
+        self.process.buildActions.patch.vectorKey = "patch"
         self.process.buildActions.galaxyStatMask = SnSelector()
         self.process.buildActions.galaxyStatMask.fluxType = "{band}_psfFlux"
 

diff --git a/python/lsst/analysis/tools/atools/photometry.py b/python/lsst/analysis/tools/atools/photometry.py
@@ -56,6 +56,9 @@ def setDefaults(self):
         self.process.buildActions.xStars.vectorKey = "coord_ra"
         self.process.buildActions.yStars = LoadVector()
         self.process.buildActions.yStars.vectorKey = "coord_dec"
+        self.process.buildActions.patch = LoadVector()
+        self.process.buildActions.patch.vectorKey = "patch"
+
         self.process.buildActions.starStatMask = SnSelector()
         self.process.buildActions.starStatMask.fluxType = "{band}_psfFlux"
 

diff --git a/python/lsst/analysis/tools/atools/simpleDiaPlot.py b/python/lsst/analysis/tools/atools/simpleDiaPlot.py
@@ -36,6 +36,7 @@ def setDefaults(self):
         self.process.buildActions.ras.vectorKey = "ra"
         self.process.buildActions.decs = LoadVector()
         self.process.buildActions.decs.vectorKey = "decl"
+
         # TODO: update column name to 'dec' once column names are standardized,
         # i.e., RFC-863
 

diff --git a/python/lsst/analysis/tools/atools/skyObject.py b/python/lsst/analysis/tools/atools/skyObject.py
@@ -44,14 +44,16 @@ def setDefaults(self):
         self.process.buildActions.statMask = SnSelector()
         self.process.buildActions.statMask.threshold = -1e12
         self.process.buildActions.statMask.fluxType = "{band}_psfFlux"
+        self.process.buildActions.patch = LoadVector()
+        self.process.buildActions.patch.vectorKey = "patch"
 
         self.produce.plot = SkyPlot()
         self.produce.plot.plotTypes = ["any"]
         self.produce.plot.plotName = "skyObject_{band}"
         self.produce.plot.xAxisLabel = "R.A. (degrees)"
         self.produce.plot.yAxisLabel = "Dec. (degrees)"
         self.produce.plot.zAxisLabel = "Sky Object ap09Flux (nJy)"
-        self.produce.plot.plotOutlines = False
+        self.produce.plot.plotOutlines = True
         self.produce.plot.fixAroundZero = True