Changeset 34415 for branches/eam_branches/ipp-20120905/psphot/src/psphotKronIterate.c

Timestamp:

Sep 7, 2012, 10:06:48 AM (14 years ago)

Author:

eugene

Message:

when growing the CR mask, copy the values into the existing readout->mask array, otherwise we lose the association to source->maskView; mask a small region for sources that have already been detected (the exclusion radius is curretnly set to 4 pixels, and is not user-specified; deactivate some test code; saturate linear fit fluxes to min valid flux (probably superfluous); calls to psphotFitSourcesLinear in psphotReadout should skip negative fluxes

File:

• branches/eam_branches/ipp-20120905/psphot/src/psphotKronIterate.c (modified) (2 diffs)

branches/eam_branches/ipp-20120905/psphot/src/psphotKronIterate.c

@@ -321 +321 @@
             if (!source->peak) continue; // XXX how can we have a peak-less source?
 
+# if (0)
+# define TEST_X 653
+# define TEST_Y 466
+        if ((fabs(source->peak->xf - TEST_X) < 5) && (fabs(source->peak->yf - TEST_Y) < 5)) {
+            fprintf (stderr, "test object\n");
+        }
+# undef TEST_X
+# undef TEST_Y
+# endif
+
             // check status of this source's moments
             if (!source->moments) continue;
@@ -541 +551 @@
     if (RF <= 0. || RS <= 0 || !isfinite(MrfTry)) {
         // We did not get a good measurement
-        source->moments->Mrf = NAN;
+        source->moments->Mrf = NAN;
         source->moments->KronFlux  = NAN;
         source->moments->KronFluxErr  = NAN;