Changeset 26923

Timestamp:

Feb 11, 2010, 6:32:28 PM (16 years ago)

Author:

Paul Price

Message:

Don't do or record pattern correction for video cells.

Location:

trunk

Files:

• ppStack/src/ppStackPhotometry.c (modified) (2 diffs)
• psModules/src/detrend/pmPatternIO.c (modified) (1 diff)

trunk/ppStack/src/ppStackPhotometry.c

@@ -66 +66 @@
     }
 
+    pmModelClassSetLimits(PM_MODEL_LIMITS_IGNORE);
     if (!psphotReadoutKnownSources(config, photView, inSources)) {
         // This is likely a data quality issue
@@ -91 +92 @@
         pmReadout *photRO = pmFPAviewThisReadout(photView, photFile->fpa); // Readout with the sources
         pmDetections *detections = psMetadataLookupPtr(NULL, photRO->analysis, "PSPHOT.DETECTIONS"); // detections
-        if (detections) {
-            psAssert (detections->allSources, "missing sources?");
-            psMetadataAddS32(options->stats, PS_LIST_TAIL, "NUM_SOURCES", 0, "Number of sources detected", detections->allSources->n);
-        } else {
-            psMetadataAddS32(options->stats, PS_LIST_TAIL, "NUM_SOURCES", 0, "Number of sources detected", 0);
-        }
+        if (detections) {
+            psAssert (detections->allSources, "missing sources?");
+            psMetadataAddS32(options->stats, PS_LIST_TAIL, "NUM_SOURCES", 0, "Number of sources detected", detections->allSources->n);
+        } else {
+            psMetadataAddS32(options->stats, PS_LIST_TAIL, "NUM_SOURCES", 0, "Number of sources detected", 0);
+        }
         psMetadataAddF32(options->stats, PS_LIST_TAIL, "TIME_PHOT", 0, "Time to do photometry", psTimerMark("PPSTACK_PHOT"));
     }

trunk/psModules/src/detrend/pmPatternIO.c

@@ -77 +77 @@
     }
 
-    pmFPAview *thisView = pmFPAviewAlloc(view->nRows); // Copy of input view
-    *thisView = *view;
-
-    for (int i = 0; i < cell->readouts->n; i++) {
-        pmReadout *readout = cell->readouts->data[i];
-        thisView->readout = i;
-        if (!pmReadoutWritePattern(readout, file->fits)) {
-            psError(PS_ERR_IO, false, "Failed to write %dth readout", i);
-            psFree(thisView);
-            return false;
-        }
-    }
-    psFree(thisView);
+    // Only do the FIRST readout --- don't want to write lots of headers
+    pmReadout *readout = cell->readouts->data[0];
+    if (!pmReadoutWritePattern(readout, file->fits)) {
+        psError(PS_ERR_IO, false, "Failed to write readout");
+        return false;
+    }
     return true;
 }