Changeset 25433 for branches/eam_branches/20090715/psphot/src/psphotPetrosianAnalysis.c

Timestamp:

Sep 17, 2009, 2:26:32 PM (17 years ago)

Author:

eugene

Message:

change radius for extended sources to use footprint; clean up some of the visualizations; plug some leaks

File:

• branches/eam_branches/20090715/psphot/src/psphotPetrosianAnalysis.c (modified) (3 diffs)

branches/eam_branches/20090715/psphot/src/psphotPetrosianAnalysis.c

@@ -12 +12 @@
     // XXX temporary user-supplied systematic sky noise measurement (derive from background model)
     float skynoise = psMetadataLookupF32 (&status, recipe, "SKY.NOISE");
+
+# if (0)
+    // if backModel or backStdev are missing, the values of sky and/or skyErr will be set to NAN
+    // XXX use this to set skynoise
+    pmReadout *backModel = psphotSelectBackground (config, view);
+    pmReadout *backStdev = psphotSelectBackgroundStdev (config, view);
+# endif
 
     // S/N limit to perform full non-linear fits
@@ -33 +40 @@
         if (source->type == PM_SOURCE_TYPE_DEFECT) continue;
         if (source->type == PM_SOURCE_TYPE_SATURATED) continue;
+        if (source->mode & PM_SOURCE_MODE_DEFECT) continue;
+        if (source->mode & PM_SOURCE_MODE_SATSTAR) continue;
+        if (!(source->mode & PM_SOURCE_MODE_EXT_LIMIT)) continue;
 
         // limit selection to some SN limit
         assert (source->peak); // how can a source not have a peak?
         if (source->peak->SN < SN_LIM) continue;
-        if (source->extNsigma < 10.0) continue; // XXX this should not be hardwired
 
         // limit selection by analysis region
@@ -50 +59 @@
         }
 
-        psphotPetrosianProfile (source, skynoise);
+        psphotPetrosianProfile (readout, source, skynoise);
+
+        pmSourceSub (source, PM_MODEL_OP_FULL, maskVal);
+        source->tmpFlags |= PM_SOURCE_TMPF_SUBTRACTED;
     }
+
+    psphotVisualShowResidualImage (readout);
 
     // pause and wait for user input: