Add high flux linearity constraint to avoid craziness.

lsst · Oct 31, 2024 · a0b7bca · a0b7bca
1 parent 65a4a88
commit a0b7bca
Showing 1 changed file with 8 additions and 1 deletion.
diff --git a/python/lsst/cp/pipe/utils.py b/python/lsst/cp/pipe/utils.py
@@ -1339,7 +1339,14 @@ def __call__(self, pars):
             log_w = np.sqrt(-2.*np.log(fact*self._w[self._w > 0]))
             constraint = np.hstack([constraint, log_w])
 
-        return np.hstack([resid, constraint])
+        # Don't let it get to >5% correction.
+        values = pars[self.par_indices["values"]]
+        if np.abs(values[-1])/self._nodes[-1] > 0.25:
+            extra_constraint = 1e10
+        else:
+            extra_constraint = 0
+
+        return np.hstack([resid, constraint, extra_constraint])
 
 
 def getReadNoise(exposure, ampName, taskMetadata=None, log=None):