Changeset 34072 for branches/eam_branches/ipp-20120601/psphot/src

Timestamp:

Jun 25, 2012, 5:15:13 PM (14 years ago)

Author:

eugene

Message:

some tests to check working analysis; add some additional stats to psphotFitSourcesLinear

File:

• branches/eam_branches/ipp-20120601/psphot/src/psphotReadout.c (modified) (1 diff)

branches/eam_branches/ipp-20120601/psphot/src/psphotReadout.c

@@ -9 +9 @@
 }
 
+# if (0)
+// TEST CODE, can be removed
+bool psphotDumpFlux (pmConfig *config, const pmFPAview *view, const char *filerule) {
+
+    bool status = false;
+
+    // find the currently selected readout
+    pmFPAfile *file = pmFPAfileSelectSingle(config->files, filerule, 0); // File of interest
+    psAssert (file, "missing file?");
+
+    pmReadout *readout = pmFPAviewThisReadout(view, file->fpa);
+    psAssert (readout, "missing readout?");
+
+    pmDetections *detections = psMetadataLookupPtr (&status, readout->analysis, "PSPHOT.DETECTIONS");
+    psAssert (detections, "missing detections?");
+
+    psArray *sources = detections->allSources;
+    psAssert (sources, "missing sources?");
+
+    static int npass;
+    char filename[64];
+    snprintf (filename, 64, "mags.%d.dat", npass);
+    FILE *ftest = fopen (filename, "w");
+    for (int j = 0; j < sources->n; j++) {
+        pmSource *source = sources->data[j];
+
+        float psfMag;
+        status = pmSourcePhotometryModel (&psfMag, NULL, source->modelPSF);
+
+        float psfMagNorm;
+        float Io = source->modelPSF->params->data.F32[PM_PAR_I0];
+        source->modelPSF->params->data.F32[PM_PAR_I0] = 1.0;
+        status = pmSourcePhotometryModel (&psfMagNorm, NULL, source->modelPSF);
+        source->modelPSF->params->data.F32[PM_PAR_I0] = Io;
+
+        // double apTrend = pmTrend2DEval (psf->ApTrend, (float)source->peak->x, (float)source->peak->y);
+        fprintf (ftest, "%d %d %d  %f %f %f %f  %f %f\n", j, source->peak->x, source->peak->y, source->modelPSF->params->data.F32[PM_PAR_I0], source->modelPSF->params->data.F32[PM_PAR_SXX], source->modelPSF->params->data.F32[PM_PAR_SYY], source->modelPSF->params->data.F32[PM_PAR_SXY], psfMag, psfMagNorm);
+    }
+    fclose (ftest);
+    npass++;
+
+    return true;
+}
+# endif
+
 bool psphotReadout(pmConfig *config, const pmFPAview *view, const char *filerule) {