Changeset 25614

Timestamp:

Sep 27, 2009, 11:11:05 AM (17 years ago)

Author:

eugene

Message:

add test dump stats function

File:

• branches/eam_branches/20090715/psphot/src/psphotOutput.c (modified) (1 diff)

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

@@ -31 +31 @@
     }
     return background;
+}
+
+// dump source stats for psf stars
+bool psphotDumpStats (psArray *sources, char *stage) {
+
+    char filename[64];
+    snprintf (filename, 64, "psf.%s.dat", stage);
+    FILE *f = fopen (filename, "w");
+    for (int i = 0; i < sources->n; i++) {
+        pmSource *source = sources->data[i];
+        if (!(source->mode & PM_SOURCE_MODE_PSFSTAR)) continue;
+
+        pmModel *model = source->modelPSF;
+        if (!model) continue;
+
+        // int xc = source->peak->x - source->pixels->col0;
+        // int yc = source->peak->y - source->pixels->row0;
+        // float mcore = source->modelFlux ? source->modelFlux->data.F32[yc][xc] : NAN;
+        // float mpeak = model ? model->params->data.F32[PM_PAR_I0] : NAN;
+        // bool subtracted = source->tmpFlags & PM_SOURCE_TMPF_SUBTRACTED;
+        // fprintf (stderr, "%d %d : %d : %f %f : %f %f\n", source->peak->x, source->peak->y, subtracted, source->peak->flux, source->pixels->data.F32[yc][xc], mcore, mpeak); 
+
+        fprintf (f, "%6.1f %6.1f : %6.1f %6.1f : %8.3f %8.3f %8.3f : %f : %f %f %f : %f\n",
+                 source->peak->xf, source->peak->yf, 
+                 model->params->data.F32[PM_PAR_XPOS], model->params->data.F32[PM_PAR_YPOS], 
+                 source->psfMag, source->apMag, source->errMag,
+                 model->params->data.F32[PM_PAR_I0], 
+                 model->params->data.F32[PM_PAR_SXX], model->params->data.F32[PM_PAR_SXY], model->params->data.F32[PM_PAR_SYY], 
+                 model->params->data.F32[PM_PAR_7]);
+    }
+    fclose (f);
+    return true;
 }