Changeset 31384 for branches/eam_branches/ipp-20110404/psphot/src/psphotRadialApertures.c

Timestamp:

Apr 26, 2011, 5:54:39 PM (15 years ago)

Author:

eugene

Message:

pmSourceCopy and pmPeakCopy should copy all appropriate values; do not fit sources with moments failures

File:

• branches/eam_branches/ipp-20110404/psphot/src/psphotRadialApertures.c (modified) (2 diffs)

branches/eam_branches/ipp-20110404/psphot/src/psphotRadialApertures.c

@@ -95 +95 @@
         if (source->type == PM_SOURCE_TYPE_SATURATED) continue;
         if (source->mode & PM_SOURCE_MODE_DEFECT) continue;
-        if (source->mode & PM_SOURCE_MODE_SATSTAR) continue;
+
+        // XXX measure radial apertures even for saturated stars
+        // if (source->mode & PM_SOURCE_MODE_SATSTAR) continue;
 
         // limit selection to some SN limit
@@ -257 +259 @@
         float SBstdv = sqrt((fluxStd->data.F32[i] / nPix) - PS_SQR(SBmean));
 
+        // XXX report the total flux or the mask-corrected flux?
         // flux->data.F32[i]    = SBmean * Area;
+        // fluxErr->data.F32[i] = sqrt(fluxErr->data.F32[i]) * Area / nPix;
+
+        fluxErr->data.F32[i] = sqrt(fluxErr->data.F32[i]);
         fluxStd->data.F32[i] = SBstdv * Area;
-        fluxErr->data.F32[i] = sqrt(fluxErr->data.F32[i]) * Area / nPix;
-
-        // fill->data.F32[i] /= Area;
+        fill->data.F32[i] /= Area;
+
         psTrace ("psphot", 5, "radial bins: %3d  %5.1f : %8.1f +/- %7.2f : %8.1f +/- %8.1f : %4.2f %6.1f\n", 
                  i, aperRadii->data.F32[i], flux->data.F32[i], fluxErr->data.F32[i], SBmean, SBstdv, fill->data.F32[i], Area);