Add WholeSkyAnalysisTask to the DRP HSC reprocessing pipelines

lsst · May 29, 2024 · 3f9e3cf · 3f9e3cf
1 parent d54c695
commit 3f9e3cf
Show file tree

Hide file tree

Showing 4 changed files with 61 additions and 0 deletions.
diff --git a/pipelines/HSC/DRP-RC2.yaml b/pipelines/HSC/DRP-RC2.yaml
@@ -69,6 +69,17 @@ tasks:
     class: lsst.pipe.tasks.isolatedStarAssociation.IsolatedStarAssociationTask
     config:
       band_order: ['i', 'r', 'g', 'z', 'y', 'N921']
+  objectTableMetricWholeSkyPlot:
+    class: lsst.analysis.tools.tasks.WholeSkyAnalysisTask
+    config:
+      atools.wholeSkyMetric: WholeSkyPlotTool
+      atools.wholeSkyMetric.xLimits: [145.0, 340.0]
+      atools.wholeSkyMetric.yLimits: [-5.0, 10.0]
+      atools.wholeSkyMetric.bands: ["g", "r", "i", "z", "y"]
+      atools.wholeSkyMetric.plotKeys: parameters.wholeSkyPlotKeysWithoutBand
+      atools.wholeSkyMetric.keysWithBand: parameters.wholeSkyPlotKeysWithBand
+      python: |
+        from lsst.analysis.tools.atools import *
 subsets:
   step1:
     subset:
@@ -244,6 +255,8 @@ subsets:
       - analyzeObjectTableCore
       - analyzeObjectTableSurveyCore
       - catalogMatchTract
+      - makeMetricTableObjectTableCore
+      - objectTableMetricWholeSkyPlot
       - photometricCatalogMatch
       - photometricRefCatObjectTract
       - plotPropertyMapTract

diff --git a/pipelines/HSC/DRP-RC2_subset.yaml b/pipelines/HSC/DRP-RC2_subset.yaml
@@ -139,6 +139,17 @@ tasks:
         config.connections.matchedCatalog = 'matchedCatalogTractStarsSNR5to80'
         from lsst.faro.measurement import ModelPhotRepTask
         config.measure.retarget(ModelPhotRepTask)
+  objectTableMetricWholeSkyPlot:
+    class: lsst.analysis.tools.tasks.WholeSkyAnalysisTask
+    config:
+      atools.wholeSkyMetric: WholeSkyPlotTool
+      atools.wholeSkyMetric.xLimits: [125.0, 175.0]
+      atools.wholeSkyMetric.yLimits: [-5.0, 10.0]
+      atools.wholeSkyMetric.bands: ["g", "r", "i", "z", "y"]
+      atools.wholeSkyMetric.plotKeys: parameters.wholeSkyPlotKeysWithoutBand
+      atools.wholeSkyMetric.keysWithBand: parameters.wholeSkyPlotKeysWithBand
+      python: |
+        from lsst.analysis.tools.atools import *
   psfPhotRepStar1:
     class: lsst.faro.measurement.TractMatchedMeasurementTask
     config:
@@ -416,6 +427,8 @@ subsets:
       - analyzeObjectTableCore
       - analyzeObjectTableSurveyCore
       - catalogMatchTract
+      - makeMetricTableObjectTableCore
+      - objectTableMetricWholeSkyPlot
       - photometricCatalogMatch
       - photometricRefCatObjectTract
       - plotPropertyMapTract

diff --git a/pipelines/HSC/DRP-ci_hsc.yaml b/pipelines/HSC/DRP-ci_hsc.yaml
@@ -50,6 +50,17 @@ tasks:
     class: lsst.pipe.tasks.isolatedStarAssociation.IsolatedStarAssociationTask
     config:
       band_order: ['i', 'r']
+  objectTableMetricWholeSkyPlot:
+    class: lsst.analysis.tools.tasks.WholeSkyAnalysisTask
+    config:
+      atools.wholeSkyMetric: WholeSkyPlotTool
+      atools.wholeSkyMetric.xLimits: [295.0, 345.0]
+      atools.wholeSkyMetric.yLimits: [-7.5, 7.5]
+      atools.wholeSkyMetric.bands: ["g", "r", "i"]
+      atools.wholeSkyMetric.plotKeys: parameters.wholeSkyPlotKeysWithoutBand
+      atools.wholeSkyMetric.keysWithBand: parameters.wholeSkyPlotKeysWithBand
+      python: |
+        from lsst.analysis.tools.atools import *
 subsets:
   analysis_tools:
     subset:
@@ -60,6 +71,8 @@ subsets:
       - analyzeObjectTableSurveyCore
       - analyzeSourceTableCore
       - catalogMatchTract
+      - makeMetricTableObjectTableCore
+      - objectTableMetricWholeSkyPlot
       - photometricCatalogMatch
       - photometricRefCatObjectTract
       - plotPropertyMapTract

diff --git a/pipelines/_ingredients/HSC/DRP.yaml b/pipelines/_ingredients/HSC/DRP.yaml
@@ -8,6 +8,28 @@ imports:
   - $ANALYSIS_TOOLS_DIR/pipelines/coaddQualityCore.yaml
   - $ANALYSIS_TOOLS_DIR/pipelines/matchedVisitQualityCore.yaml
   - $ANALYSIS_TOOLS_DIR/pipelines/visitQualityCore.yaml
+  - $ANALYSIS_TOOLS_DIR/pipelines/wholeSkyCore.yaml
+parameters:
+  wholeSkyPlotKeysWithBand: [
+    "shapeSizeFractionalDiff_{band}_highSNStars_median",
+    "shapeSizeFractionalDiff_{band}_highSNStars_sigmaMad",
+    "e1Diff_{band}_highSNStars_median",
+    "e1Diff_{band}_highSNStars_sigmaMad",
+    "e2Diff_{band}_highSNStars_median",
+    "e2Diff_{band}_highSNStars_sigmaMad",
+    "skyFluxStatisticMetric_{band}_medianSky",
+    "skyFluxStatisticMetric_{band}_meanSky",
+    "skyFluxStatisticMetric_{band}_stdevSky",
+  ]
+  wholeSkyPlotKeysWithoutBand: [
+    "wPerpPSFP_wPerp_psfFlux_sigmaMAD",
+    "yPerpPSF_yPerp_psfFlux_sigmaMAD",
+    "skippedDeblenderMetrics_numSkippedPeaks",
+    "skippedDeblenderMetrics_numSkippedBlends",
+    "skippedDeblenderMetrics_numBlendParentTooBig",
+    "skippedDeblenderMetrics_numBlendTooManyPeaks",
+    "skippedDeblenderMetrics_numBlendTooManyMasked",
+  ]
 tasks:
   isolatedStarAssociation:
     class: lsst.pipe.tasks.isolatedStarAssociation.IsolatedStarAssociationTask