Changeset 20116

Timestamp:

Oct 13, 2008, 3:21:07 PM (18 years ago)

Author:

Paul Price

Message:

Adding source output: more values than is in a CMF file.

File:

• trunk/ppSub/src/ppSubReadout.c (modified) (3 diffs)

trunk/ppSub/src/ppSubReadout.c

@@ -547 +547 @@
                              psTimerClear("PPSUB_PHOT"));
         }
+
+#ifdef TESTING
+        // Record data about sources: not everything gets into the output CMF files
+        {
+            pmReadout *photRO = pmFPAviewThisReadout(view, photFile->fpa); // Readout with the sources
+            psArray *sources = psMetadataLookupPtr(NULL, photRO->analysis, "PSPHOT.SOURCES"); // Sources
+            FILE *sourceFile = fopen("sources.dat", "w"); // File for sources
+            fprintf(sourceFile, "# x y mag mag_err psf_chisq cr_nsigma ext_nsigma psf_qf flags m_x m_y m_xx m_xy m_yy\n");
+            for (int i = 0; i < sources->n; i++) {
+                pmSource *source = sources->data[i];
+                if (!source) {
+                    continue;
+                }
+
+                float x, y;             // Position of source
+                float chi2;             // chi^2 for source
+                if (source->modelPSF) {
+                    x = source->modelPSF->params->data.F32[PM_PAR_XPOS];
+                    y = source->modelPSF->params->data.F32[PM_PAR_YPOS];
+                    chi2 = source->modelPSF->chisq;
+                } else if (source->peak) {
+                    x = source->peak->xf;
+                    y = source->peak->yf;
+                    chi2 = NAN;
+                } else {
+                    psWarning("No position available for source.");
+                    continue;
+                }
+
+                float xMoment = NAN, yMoment = NAN, xxMoment = NAN, xyMoment = NAN, yyMoment = NAN;
+                if (source->moments) {
+                    xMoment = source->moments->Mx;
+                    yMoment = source->moments->My;
+                    xxMoment = source->moments->Mxx;
+                    xyMoment = source->moments->Mxy;
+                    yyMoment = source->moments->Myy;
+                }
+
+                fprintf(sourceFile, "%f %f %f %f %f %f %f %f %d %f %f %f %f %f",
+                        x, y, source->psfMag, source->errMag, chi2, source->crNsigma, source->extNsigma,
+                        source->pixWeight, source->mode, xMoment, yMoment, xxMoment, xyMoment, yyMoment);
+            }
+            fclose(sourceFile);
+        }
+#endif
+
     }
 
@@ -594 +640 @@
     }
 
+#ifdef TESTING
     // Significance image
     {
@@ -603 +650 @@
         psFree(sig);
     }
+#endif
 
     psFree(outRO);