Changeset 32628 for branches/eam_branches/ipp-20110906/psphot/src/psphotRadialApertures.c

Timestamp:

Nov 8, 2011, 2:41:42 PM (15 years ago)

Author:

eugene

Message:

match sources before final linear fit (to enable forced photometry of PSF mags); only a single detection per image allowed for an object; distinguish conditions for Radial Aperture analysis and Petrosian analysis; no Radial Aperture analysis for chisq image

File:

• branches/eam_branches/ipp-20110906/psphot/src/psphotRadialApertures.c (modified) (3 diffs)

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

@@ -26 +26 @@
     int num = psphotFileruleCount(config, filerule);
 
+    // skip the chisq image (optionally?)
+    int chisqNum = psMetadataLookupS32 (&status, config->arguments, "PSPHOT.CHISQ.NUM");
+    if (!status) chisqNum = -1;
+
     // loop over the available readouts
     for (int i = 0; i < num; i++) {
+        if (i == chisqNum) continue; // skip chisq image
+
         if (!psphotRadialAperturesReadout (config, view, filerule, i, recipe, entry)) {
             psError (PSPHOT_ERR_CONFIG, false, "failed on measure extended source aperture-like parameters for %s entry %d", filerule, i);
@@ -218 +224 @@
         pmSource *source = sources->data[i];
 
+        // if we have checked the source validity on the basis of the object set, then 
+        // we either skip these tests below or we skip the source completely
+        if (source->tmpFlags & PM_SOURCE_TMPF_RADIAL_SKIP) continue;
+        if (source->tmpFlags & PM_SOURCE_TMPF_RADIAL_KEEP) goto keepSource;
+
         // skip PSF-like and non-astronomical objects
         if (source->type == PM_SOURCE_TYPE_DEFECT) continue;
@@ -235 +246 @@
         if (source->peak->x > region->x1) continue;
         if (source->peak->y > region->y1) continue;
+
+    keepSource:
 
         // allocate pmSourceExtendedParameters, if not already defined