Changeset 13864

Timestamp:

Jun 18, 2007, 12:32:49 PM (19 years ago)

Author:

eugene

Message:

adding trace dump output; skip sources with apResid not finite

File:

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

trunk/psphot/src/psphotApResid.c

@@ -62 +62 @@
     Npsf = 0;
 
+    FILE *dumpFile = NULL;
+    if (psTraceGetLevel("psphot") >= 5) {
+        dumpFile = fopen ("apresid.dat", "w");
+    }
+
     // select all good PM_SOURCE_TYPE_STAR entries
     for (int i = 0; i < sources->n; i++) {
@@ -80 +85 @@
         }
 
+        if (!isfinite(source->apMag) || !isfinite(source->apMag)) {
+            Nfail ++;
+            continue;
+        }
+
+        // sanity clipping : if the model is so discrepant, but your expectation in the recipe
         apResid->data.F64[Npsf] = source->apMag - source->psfMag;
 
-        xPos->data.F64[Npsf]    = model->params->data.F32[PM_PAR_XPOS];
-        yPos->data.F64[Npsf]    = model->params->data.F32[PM_PAR_YPOS];
-
-        flux->data.F64[Npsf]    = pow(10.0, -0.4*source->psfMag);
-        r2rflux->data.F64[Npsf] = PS_SQR(model->radiusFit) / flux->data.F64[Npsf];
-
-        mask->data.U8[Npsf] = 0;
-
-        // XXX sanity clip?
-        // XXX need to see if all data were tossed?
-        // XXX need to subtract median?
-        // XXX need to put this in the config data...
-        // if (fabs(apResid->data.F64[Npsf]) > 0.2) continue;
         if ((MAX_AP_OFFSET > 0) && (fabs(apResid->data.F64[Npsf]) > MAX_AP_OFFSET)) {
             Nfail ++;
@@ -100 +98 @@
         }
 
+        xPos->data.F64[Npsf]    = model->params->data.F32[PM_PAR_XPOS];
+        yPos->data.F64[Npsf]    = model->params->data.F32[PM_PAR_YPOS];
+
+        flux->data.F64[Npsf]    = pow(10.0, -0.4*source->psfMag);
+        r2rflux->data.F64[Npsf] = PS_SQR(model->radiusFit) / flux->data.F64[Npsf];
+
+        mask->data.U8[Npsf] = 0;
+
         dMag->data.F64[Npsf] = model->dparams->data.F32[PM_PAR_I0] / model->params->data.F32[PM_PAR_I0];
+
+        if (psTraceGetLevel("psphot") >= 5) {
+            fprintf (dumpFile, "%f %f  %f %f %f  %x\n", 
+                     xPos->data.F64[Npsf], yPos->data.F64[Npsf], 
+                     source->psfMag, dMag->data.F64[Npsf], apResid->data.F64[Npsf], 
+                     mask->data.U8[Npsf]);
+        }
 
         psVectorExtend (mask,    100, 1);
@@ -111 +124 @@
         Npsf ++;
     }
+
+    if (psTraceGetLevel("psphot") >= 5) {
+        fclose (dumpFile);
+    }
+    
     psLogMsg ("psphot.apresid", PS_LOG_DETAIL, "measure aperture residuals for %d objects (%d skipped, %d failed, %ld invalid)\n",
               Npsf, Nskip, Nfail, sources->n - Npsf - Nskip - Nfail);