Changeset 8138

Timestamp:

Aug 4, 2006, 2:38:15 AM (20 years ago)

Author:

eugene

Message:

track failures in log message, trace for source results

File:

• trunk/psphot/src/psphotBlendFit.c (modified) (4 diffs)

trunk/psphot/src/psphotBlendFit.c

@@ -5 +5 @@
 
     int Nfit = 0;
+    int Npsf = 0;
+    int Next = 0;
+    int Nfail = 0;
     bool status;
 
@@ -49 +52 @@
         // skip sources which are insignificant flux?
         if (source->modelPSF->params->data.F32[1] < 0.1) {
-            psTrace ("psphotBlendFit", 5, "skipping near zero source: %f, %f : %f\n",
+            psTrace ("psphotBlendFit", 5, "skipping near-zero source: %f, %f : %f\n",
                      source->modelPSF->params->data.F32[1],
                      source->modelPSF->params->data.F32[2],
@@ -60 +63 @@
         Nfit ++;
 
-        psTrace ("psphot.blend", 5, "trying source at %f, %f\n", source->moments->x, source->moments->y);
-        
         // try fitting PSFs, then try extended sources
-        if (psphotFitBlend (readout, source, psf)) continue;
-        if (psphotFitBlob (readout, source, sources, psf)) continue;
+        if (psphotFitBlend (readout, source, psf)) { 
+            psTrace ("psphotBlendFit", 5, "source at %7.1f, %7.1f is psf", source->moments->x, source->moments->y);
+            Npsf ++;
+            continue;
+        }
+        if (psphotFitBlob (readout, source, sources, psf)) {
+            psTrace ("psphotBlendFit", 5, "source at %7.1f, %7.1f is ext", source->moments->x, source->moments->y);
+            Next ++;
+            continue;
+        }
 
-        psTrace ("psphot.blend", 5, "failed fits\n");
+        psTrace ("psphotBlendFit", 5, "source at %7.1f, %7.1f failed", source->moments->x, source->moments->y);
+        Nfail ++;
 
         // re-subtract PSF for object, leave local sky
@@ -74 +84 @@
     }
 
-    psLogMsg ("psphot", 3, "fit PSF models: %f sec for %d objects (%d total)\n", psTimerMark ("psphot"), Nfit, sources->n);
+    psLogMsg ("psphot", 3, "fit models: %f sec for %d objects (%d psf, %d ext, %d failed, %d skipped)\n", psTimerMark ("psphot"), Nfit, Npsf, Next, Nfail, sources->n - Nfit);
     return (true);
 }