Changeset 17630 for branches/eam_branch_20080511/ppSim/src/ppSimCreate.c

Timestamp:

May 11, 2008, 10:10:06 AM (18 years ago)

Author:

eugene

Message:

working on injecting fake sources in real images (also forced photometry)

File:

• branches/eam_branch_20080511/ppSim/src/ppSimCreate.c (modified) (3 diffs)

branches/eam_branch_20080511/ppSim/src/ppSimCreate.c

@@ -40 +40 @@
     }
 
-    // define the output image file
+    // define the output image file -- this is the basis for the ppSimLoop
     pmFPAfile *file = pmFPAfileDefineOutput(config, fpa, OUTPUT_FILE);
     if (!file) {
@@ -53 +53 @@
         return NULL;
     }
+    // XXX we should not require the output image to be written
     file->save = true;
 
@@ -58 +59 @@
     config->formatName = psStringCopy (file->formatName);
 
-    // have we supplied a psf model?
-    if (psMetadataLookupPtr(NULL, config->arguments, "PSPHOT.PSF")) {
-        bool status = false;
+    // For photometry, we operate on the chip-mosaicked image.  we create a copy of the mosaicked
+    // image for psphot so we can write out a clean image
+    if (options->doPhotom) {
 
-        // tie the psf file to the chipMosaic 
-        pmFPAfileBindFromArgs(&status, file, config, "PSPHOT.PSF.LOAD", "PSPHOT.PSF");
-        if (!status) {
-            psError(PS_ERR_UNKNOWN, false, "Failed to find/build PSPHOT.PSF.LOAD");
-            psFree(fpa);
-            psFree(file);
-            return NULL;
+        // we need a chip image if we perform photometry (is it necessary to build it if we don't use it?)
+        pmFPAfile *chipImage = pmFPAfileDefineChipMosaic(config, output->fpa, "PPSIM.CHIP");
+        if (!chipImage) {
+            psError(PS_ERR_IO, false, _("Unable to generate new file from PPSIM.CHIP"));
+            psFree(options);
+            return NULL;
+        }
+        if (chipImage->type != PM_FPA_FILE_IMAGE) {
+            psError(PS_ERR_IO, true, "PPSIM.CHIP is not of type IMAGE");
+            psFree(options);
+            return NULL;
+        }
+
+        // this file is just used as a carrier; output files (eg, PSPHOT.RESID) are defined by
+        // psphotDefineFiles
+        pmFPAfile *psphotInput = pmFPAfileDefineFromFile (config, chipImage, 1, 1, "PSPHOT.INPUT");
+        PS_ASSERT (psphotInput, false);
+
+        // define associated psphot input/output files
+        if (!psphotDefineFiles (config, psphotInput)) {
+            psError(PSPHOT_ERR_CONFIG, false, "Trouble defining the additional input/output files for psphot");
+            return false;
         }
+    } else {
+        // have we supplied a psf model?  this happens in psphotDefineFiles if we request a photometry
+        // analysis.  however, even if we do not, a psf model may be used to generate the fake
+        // sources.
+        if (psMetadataLookupPtr(NULL, config->arguments, "PSPHOT.PSF")) {
+            // tie the psf file to the chipMosaic 
+            pmFPAfileBindFromArgs(&status, file, config, "PSPHOT.PSF.LOAD", "PSPHOT.PSF");
+            if (!status) {
+                psError(PS_ERR_UNKNOWN, false, "Failed to find/build PSPHOT.PSF.LOAD");
+                psFree(fpa);
+                psFree(file);
+                return NULL;
+            }
+        }
     }
 
+    // PPSIM.SOURCES carries the constructed, fake sources with their true parameters
     // XXX only invoke this code for OBJECT types of images
-    // PPSIM.SOURCES carries the constructed, fake sources with their true parameters
     pmFPAfile *simSources = pmFPAfileDefineOutput (config, file->fpa, "PPSIM.SOURCES");
     if (!simSources) {