Changeset 20865

Timestamp:

Dec 1, 2008, 4:02:38 PM (17 years ago)

Author:

Paul Price

Message:

Add softening parameter to weights when calculating crNsigma.

File:

• trunk/psphot/src/psphotSourceSize.c (modified) (2 diffs)

trunk/psphot/src/psphotSourceSize.c

@@ -34 +34 @@
         psError(PS_ERR_BAD_PARAMETER_VALUE, true, "PSPHOT.CR.GROW is not positive.");
         return false;
+    }
+
+    float soft = psMetadataLookupF32(&status, recipe, "PSPHOT.CR.NSIGMA.SOFTEN"); // Softening parameter
+    if (!status || !isfinite(soft) || soft < 0.0) {
+        psWarning("PSPHOT.CR.NSIGMA.SOFTEN not set; defaulting to zero.");
+        soft = 0.0;
     }
 
@@ -113 +119 @@
         float dcR = 4*weight[yPeak][xPeak] + weight[yPeak+1][xPeak-1] + weight[yPeak-1][xPeak+1];
         float nR = cR / sqrtf(dcR);
+
+        // Soften weights (add systematic error)
+        float softening = 6 * soft * source->peak->flux; // Softening for weights
+        dCX += softening;
+        dCY += softening;
+        dCL += softening;
+        dCR += softening;
 
         // P(chisq > chisq_obs; Ndof) = gamma_Q (Ndof/2, chisq/2)