Changeset 24951 for branches/czw_branch/cleanup/magic/censorObjects/src/censorLoop.c

Timestamp:

Jul 30, 2009, 5:20:29 PM (17 years ago)

Author:

watersc1

Message:

Finished up my edits to the detrend cleanup, and some changes to my
copy of burntool and the pslib astrometry. Detrend cleanup has not
been tested yet. That's up next.

Location:

branches/czw_branch/cleanup

Files:

• . (modified) (1 prop)
• magic/censorObjects/src/censorLoop.c (modified) (3 diffs)

branches/czw_branch/cleanup/magic/censorObjects/src/censorLoop.c

@@ -43 +43 @@
     }
 
+    // select the input data sources
+    pmFPAfile *inMask = psMetadataLookupPtr (NULL, config->files, "CENSOR.INPUT.MASK");
+    if (!inMask) {
+        psError(CENSOR_ERR_CONFIG, true, "Can't find input masks!\n");
+        return false;
+    }
+
     pmFPAview *view = pmFPAviewAlloc (0);
+    bool badMaskFile = false;
+
+    // check on the mask file at this level if possible
+    badMaskFile = censorBadFileByView(badMaskFile, inMask, config, view);
 
     // files associated with the science image
@@ -49 +60 @@
 
     while ((chip = pmFPAviewNextChip (view, input->fpa, 1)) != NULL) {
-        psTrace ("censor", 4, "Chip %d: %x %x\n", view->chip, chip->file_exists, chip->process);
-        if (!chip->process || !chip->file_exists) { continue; }
+        psTrace ("censor", 4, "Chip %d: %x %x\n", view->chip, chip->file_exists, chip->process);
+        if (!chip->process || !chip->file_exists) { continue; }
+
+        // check on the mask file at this level if possible
+        badMaskFile = censorBadFileByView(badMaskFile, inMask, config, view);
+
         if (!pmFPAfileIOChecks (config, view, PM_FPA_BEFORE)) ESCAPE;
 
         while ((cell = pmFPAviewNextCell (view, input->fpa, 1)) != NULL) {
-            psTrace ("censor", 4, "Cell %d: %x %x\n", view->cell, cell->file_exists, cell->process);
-            if (!cell->process || !cell->file_exists) { continue; }
+            psTrace ("censor", 4, "Cell %d: %x %x\n", view->cell, cell->file_exists, cell->process);
+            if (!cell->process || !cell->file_exists) { continue; }
+
+            // check on the mask file at this level if possible
+            badMaskFile = censorBadFileByView(badMaskFile, inMask, config, view);
+
             if (!pmFPAfileIOChecks (config, view, PM_FPA_BEFORE)) ESCAPE;
 
             // process each of the readouts
             while ((readout = pmFPAviewNextReadout (view, input->fpa, 1)) != NULL) {
+
+                // check on the mask file at this level if possible
+                badMaskFile = censorBadFileByView(badMaskFile, inMask, config, view);
+
                 if (!pmFPAfileIOChecks (config, view, PM_FPA_BEFORE)) ESCAPE;
                 if (!readout->data_exists) { continue; }
 
-                if (!censorReadout (readout, config, view)) ESCAPE;
+                if (!censorReadout (readout, config, view, badMaskFile)) ESCAPE;
 
                 if (!pmFPAfileIOChecks (config, view, PM_FPA_AFTER)) ESCAPE;
@@ -78 +101 @@
     return true;
 }
+
+bool censorBadFileByView (bool badFile, pmFPAfile *file, pmConfig *config, pmFPAview *view) {
+
+    // get the current level
+    pmFPALevel level = pmFPAviewLevel (view);
+
+    // can we try to open this file? (if not, leave the existing badView intact)
+    if (level != file->fileLevel) return badFile;
+
+    // free, and if needed, realloc badView below
+    badFile = false;
+
+    char *filename = pmFPAfileName(file, view, config);
+    if (filename) {
+        // activate the file
+        pmFPAfileActivate (config->files, true, file->name);
+        psFree (filename);
+        return false;
+    }
+
+    // de-activate the file
+    pmFPAfileActivate (config->files, false, file->name);
+    return true;
+}