Add parameter nIterSigmaClipPtc

lsst · May 13, 2020 · 92a6008 · 92a6008
1 parent c37e90b
commit 92a6008
Showing 1 changed file with 7 additions and 1 deletion.
diff --git a/python/lsst/cp/pipe/ptc.py b/python/lsst/cp/pipe/ptc.py
@@ -147,9 +147,14 @@ class MeasurePhotonTransferCurveTaskConfig(pexConfig.Config):
     )
     nSigmaClipPtc = pexConfig.Field(
         dtype=float,
-        doc="Sigma cut for afwMatch.StatisticsControl()",
+        doc="Sigma cut for afwMath.StatisticsControl()",
         default=5.5,
     )
+    nIterSigmaClipPtc = pexConfig.Field(
+        dtype=int,
+        doc="Number of sigma-clipping iterations for afwMath.StatisticsControl()",
+        default=1,
+    )
     maxIterationsPtcOutliers = pexConfig.Field(
         dtype=int,
         doc="Maximum number of iterations for outlier rejection in PTC.",
@@ -541,6 +546,7 @@ def measureMeanVarPair(self, exposure1, exposure2, region=None):
 
         statsCtrl = afwMath.StatisticsControl()
         statsCtrl.setNumSigmaClip(self.config.nSigmaClipPtc)
+        statsCtrl.setNumIter(self.config.nIterSigmaClipPtc)
         #  Clipped mean of images; then average of mean.
         mu1 = afwMath.makeStatistics(im1Area, afwMath.MEANCLIP, statsCtrl).getValue()
         mu2 = afwMath.makeStatistics(im2Area, afwMath.MEANCLIP, statsCtrl).getValue()