Changeset 30784 for branches/eam_branches/ipp-20110213/psphot/src/psphotReadout.c

Timestamp:

Mar 3, 2011, 3:11:09 PM (15 years ago)

Author:

eugene

Message:

fix the calculation of chisq values for constant errors (re-calculated, do not attempt to correct raw value)

File:

• branches/eam_branches/ipp-20110213/psphot/src/psphotReadout.c (modified) (5 diffs)

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

@@ -9 +9 @@
 }
 
+// for now, let's store the detections on the readout->analysis for each readout
+bool psphotDumpChisqs (pmConfig *config, const pmFPAview *view, const char *filerule)
+{
+    static int npass = 0;
+    char filename[64];
+
+    bool status = true;
+
+    int num = psphotFileruleCount(config, filerule);
+
+    snprintf (filename, 64, "chisq.%02d.dat", npass);
+    FILE *f = fopen (filename, "w");
+
+    // loop over the available readouts
+    for (int i = 0; i < num; i++) {
+
+        // find the currently selected readout
+        pmFPAfile *file = pmFPAfileSelectSingle(config->files, filerule, i); // 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?");
+
+        for (int i = 0; i < sources->n; i++) {
+            pmSource *source = sources->data[i];
+            if (!source) continue;
+
+            pmModel *model = pmSourceGetModel (NULL, source);
+            if (!model) continue;
+        
+            if (source->mode & PM_SOURCE_MODE_NONLINEAR_FIT) {
+                fprintf (f, "%f %f %f %d %d %f  1 NONLINEAR\n", model->mag, model->params->data.F32[1], model->chisq, model->nDOF, model->nPix, model->chisqNorm);
+            } else {
+                fprintf (f, "%f %f %f %d %d %f  0 LINEAR\n", model->mag, model->params->data.F32[1], model->chisq, model->nDOF, model->nPix, model->chisqNorm);
+            }
+        }
+    }
+    fclose (f);
+    npass ++;
+
+    return true;
+}
+
 bool psphotReadout(pmConfig *config, const pmFPAview *view, const char *filerule) {
 
@@ -147 +196 @@
     // linear PSF fit to source peaks, subtract the models from the image (in PSF mask)
     psphotFitSourcesLinear (config, view, filerule, false); // pass 1 (detections->allSources)
+    psphotDumpChisqs (config, view, filerule);
 
     // identify CRs and extended sources (only unmeasured sources are measured)
@@ -157 +207 @@
     // replace model flux, adjust mask as needed, fit, subtract the models (full stamp)
     psphotBlendFit (config, view, filerule); // pass 1 (detections->allSources)
+    psphotDumpChisqs (config, view, filerule);
 
     // replace all sources
@@ -164 +215 @@
     // NOTE : apply to ALL sources (extended + psf)
     psphotFitSourcesLinear (config, view, filerule, true); // pass 2 (detections->allSources)
+    psphotDumpChisqs (config, view, filerule);
 
     // if we only do one pass, skip to extended source analysis
@@ -209 +261 @@
         // NOTE: apply to ALL sources
         psphotFitSourcesLinear (config, view, filerule, true); // pass 3 (detections->allSources)
+        psphotDumpChisqs (config, view, filerule);
     }