Changeset 26629 for trunk/ppImage/src/ppImageDetrendReadout.c

Timestamp:

Jan 19, 2010, 2:30:36 PM (16 years ago)

Author:

Paul Price

Message:

Recipe value GAIN.OVERRIDE specifies whether the gain should be set to unity in the event it is non-finite in ppImage and ppMerge. This seems to be a problem for GPC1, where some gains are not set because they haven't been measured. GAIN.OVERRIDE is turned on for ppImage on GPC1, and off for everything else.

File:

• trunk/ppImage/src/ppImageDetrendReadout.c (modified) (1 diff)

trunk/ppImage/src/ppImageDetrendReadout.c

@@ -89 +89 @@
             psFree (binning);
         }
+
+        // Check that the gain is set
+        float gain = psMetadataLookupF32(NULL, input->parent->concepts, "CELL.GAIN"); // Gain for cell
+        if (!isfinite(gain)) {
+            psMetadata *recipe = psMetadataLookupMetadata(NULL, config->recipes, RECIPE_NAME); // Recipe
+            psAssert(recipe, "Should be there!");
+            bool override = psMetadataLookupBool(NULL, recipe, "GAIN.OVERRIDE"); // Override the bad gain?
+            if (override) {
+                psWarning("CELL.GAIN is not set for readout (%d,%d,%d) --- setting to unity.",
+                          view->chip, view->cell, view->readout);
+                psMetadataItem *item = psMetadataLookup(input->parent->concepts, "CELL.GAIN"); // Gain item
+                psAssert(item, "Should be there!");
+                item->data.F32 = 1.0;
+            } else {
+                psError(PS_ERR_BAD_PARAMETER_VALUE, false, "CELL.GAIN is not set for readout (%d,%d,%d)",
+                        view->chip, view->cell, view->readout);
+                return false;
+            }
+        }
+
         pmReadoutGenerateVariance(input, noiseImage, true);
         psFree (noiseImage);