Changeset 11352 for trunk/ppImage/src/ppImagePhotom.c

Timestamp:

Jan 27, 2007, 5:02:24 PM (19 years ago)

Author:

eugene

Message:

update psphot and psastro in ppImage

File:

• trunk/ppImage/src/ppImagePhotom.c (modified) (6 diffs)

trunk/ppImage/src/ppImagePhotom.c

@@ -7 +7 @@
 // the top portion of this function is equivalent to psphotParseCamera,
 // but different in a few important ways.  
+// XXX move this section to ppImageParseCamera?
 bool ppImagePhotom (pmConfig *config, pmFPAview *view) {
 
@@ -27 +28 @@
         // we set the freeLevel to be FPA so this pmFPAfile may be used as input to psastro
         // XXX make this the resonsibility of a ppImageAstrom call?
-        pmFPAfile *outfile = pmFPAfileDefineOutput (config, input->fpa, "PSPHOT.OUTPUT");
-        outfile->freeLevel = PM_FPA_LEVEL_FPA;
+        // XXX test for psastro requests?
+        // XXX need to make this an INTERNAL file if output is not requested.
+        bool saveOutput = psMetadataLookupBool (NULL, recipe, "SAVE.OUTPUT");
+        if (saveOutput) {
+          pmFPAfile *outfile = pmFPAfileDefineOutput (config, input->fpa, "PSPHOT.OUTPUT");
+          outfile->freeLevel = PM_FPA_LEVEL_FPA;
+        }
     }
 
     // if we have requested PSPHOT.SRC (externally supplied sources), attempt to resolve it 
     if (psMetadataLookupPtr(NULL, config->arguments, "PSPHOT.SRC")) { 
-        pmFPAfileDefineFromArgs (&status, config, "PSPHOT.SRC", "SRC");
+        pmFPAfileDefineFromArgs (&status, config, "PSPHOT.SRC", "PSPHOT.SRC");
         if (!status) {
             psError(PSPHOT_ERR_CONFIG, false, "Failed to find/build PSPHOT.SRC");
@@ -42 +48 @@
     // select recipe options supplied on command line
     psMetadata *recipe  = psMetadataLookupPtr (&status, config->recipes, PSPHOT_RECIPE);
+
+    // NOTE the difference here from psphotParseCamera: the mask and weight
+    // images are supplied by ppImage and need not / should not be loaded here.
+
+    // optionally load the PSF Model and/or fixed stars
+    pmFPAfileBindFromArgs (NULL, input, config, "PSPHOT.PSF", "PSPHOT.PSF");
 
     // set default recipe values here
@@ -60 +72 @@
     // optionally save the residual image 
     if (psMetadataLookupBool(NULL, recipe, "SAVE.RESID")) {
-        pmFPAfileDefineOutput (config, input->fpa, "PSPHOT.RESID");
+        if (!pmFPAfileDefineOutput (config, input->fpa, "PSPHOT.RESID")) {
+            psError(PSPHOT_ERR_CONFIG, false, "Cannot find a rule for PSPHOT.BACKMDL");
+            return false;
+        }
     }
 
@@ -70 +85 @@
     // optionally save the background model (small FITS image)
     if (psMetadataLookupBool(NULL, recipe, "SAVE.BACKMDL")) {
-        pmFPAfileDefineFromFPA (config, input->fpa, DX, DY, "PSPHOT.BACKMDL");
+        if (!pmFPAfileDefineFromFPA (config, input->fpa, DX, DY, "PSPHOT.BACKMDL")) {
+            psError(PSPHOT_ERR_CONFIG, false, "Cannot find a rule for PSPHOT.BACKMDL");
+            return false;
+        }
     }
     // optionally save the full background image
     if (psMetadataLookupBool(NULL, recipe, "SAVE.BACKGND")) {
-        pmFPAfileDefineFromFPA (config, input->fpa,  1,  1, "PSPHOT.BACKGND");
+        if (!pmFPAfileDefineFromFPA (config, input->fpa,  1,  1, "PSPHOT.BACKGND")) {
+            psError(PSPHOT_ERR_CONFIG, false, "Cannot find a rule for PSPHOT.BACKGND");
+            return false;
+        }
     }
     // optionally save the background-subtracted image
     if (psMetadataLookupBool(NULL, recipe, "SAVE.BACKSUB")) {
-        pmFPAfileDefineFromFPA (config, input->fpa,  1,  1, "PSPHOT.BACKSUB");
+        if (!pmFPAfileDefineFromFPA (config, input->fpa,  1,  1, "PSPHOT.BACKSUB")) {
+            psError(PSPHOT_ERR_CONFIG, false, "Cannot find a rule for PSPHOT.");
+            return false;
+        }
     }
     // optionally save the PSF Model
     if (psMetadataLookupBool(NULL, recipe, "SAVE.PSF")) {
-        pmFPAfileDefineOutput (config, input->fpa, "PSPHOT.PSF.SAVE");
+        if (!pmFPAfileDefineOutput (config, input->fpa, "PSPHOT.PSF.SAVE")) {
+            psError(PSPHOT_ERR_CONFIG, false, "Cannot find a rule for PSPHOT.PSF.SAVE");
+            return false;
+        }
     }
-    // optionally load the PSF Model
-    if (psMetadataLookupBool(NULL, recipe, "LOAD.PSF")) {
-        pmFPAfileDefineInput (config, input->fpa, "PSPHOT.PSF.LOAD");
+
+    // optionally save output plots
+    if (psMetadataLookupBool(NULL, recipe, "SAVE.PLOTS")) {
+        if (!pmFPAfileDefineOutput (config, input->fpa, "SOURCE.PLOT.MOMENTS")) {
+            psTrace ("psphot", 3, "Cannot find a rule for SOURCE.PLOT.MOMENTS");
+        }
+        if (!pmFPAfileDefineOutput (config, input->fpa, "SOURCE.PLOT.PSFMODEL")) {
+            psTrace ("psphot", 3, "Cannot find a rule for SOURCE.PLOT.PSFMODEL");
+        }
     }
 
@@ -92 +125 @@
     // pmFPAfileConstruct (config->files, format, config->camera, "PSPHOT.PSF_SAMPLE");
 
+    // XXX probably need to deactivate all files and activate the psphot ones explicitly
+
+    // In this section, we iterate over the cells and readout for this chip
+    // and run the psphot analysis routine here
     // XXX worry about this not having unique PHOTCODES...
     // view->cell = -1;