Changeset 34179 for branches/eam_branches/ipp-20120627/psphot/src/psphotSourceSize.c

Timestamp:

Jul 19, 2012, 12:02:03 PM (14 years ago)

Author:

eugene

Message:

subtract background in psphotReadoutMinimal (should be optional); set ApResidPSF to 0.0 if PSF stars not available

File:

• branches/eam_branches/ipp-20120627/psphot/src/psphotSourceSize.c (modified) (2 diffs)

branches/eam_branches/ipp-20120627/psphot/src/psphotSourceSize.c

@@ -456 +456 @@
     psImageMaskType maskVal = options->maskVal | options->markVal;
 
+    // in the ppSub context, do we get sensible values for ApResid?
+    float ApResidPSF = options->ApResid;
+    if (!isfinite(ApResidPSF)) {
+      ApResidPSF = 0.0;
+    }
+    float ApSysErrPSF = options->ApSysErr;
+    if (!isfinite(ApSysErrPSF)) {
+      ApSysErrPSF = 0.0;
+    }
+
     for (psS32 i = 0 ; i < sources->n ; i++) {
 
@@ -521 +531 @@
         }
 
-        // XXX in the ppSub context, do we get sensible values for ApResid?
-
         // set nSigmaMAG to include both systematic and poisson error terms.  we include a hard
         // floor on the Ap Sys Err (to be a bit generous).  XXX put the floor in the recipe...
-        float nSigmaMAG = (dMag - options->ApResid) / hypot(source->psfMagErr, hypot(options->ApSysErr, 0.02));
+        float nSigmaMAG = (dMag - ApResidPSF) / hypot(source->psfMagErr, hypot(ApSysErrPSF, 0.02));
         source->extNsigma = nSigmaMAG;