Changeset 17935 for branches/eam_branch_20080511/ppSim/src/ppSimPhotomReadout.c

Timestamp:

Jun 5, 2008, 3:29:09 AM (18 years ago)

Author:

eugene

Message:

added forced photometry code; compiles, not tested

File:

• branches/eam_branch_20080511/ppSim/src/ppSimPhotomReadout.c (modified) (4 diffs)

branches/eam_branch_20080511/ppSim/src/ppSimPhotomReadout.c

@@ -37 +37 @@
     }
 
+# if 0    
     // set the photcode for this image
     if (!psphotAddPhotcode (recipe, config, view, "PPSIM.CHIP")) {
@@ -42 +43 @@
         return false;
     }
+# endif
 
     // find the currently selected readout. 
@@ -119 +121 @@
     psFree (sources); // only frees the merged references
 
-    // linear fit to real + forced sources
-    sources = ppSimMergeSources (realSources, forceSources);
-    psphotFitSourcesLinear (readout, sources, recipe, psf, FALSE); // XXX option to NOT subtract
-    psphotReplaceAllSources (sources, recipe);
-    psFree (sources); // only frees the merged references
+    // linear fit to real + forced sources (if any were selected)
+    if (forceSources->n) {
+      sources = ppSimMergeSources (realSources, forceSources);
+      psphotFitSourcesLinear (readout, sources, recipe, psf, FALSE); // XXX option to NOT subtract
+      psphotReplaceAllSources (sources, recipe);
+      psFree (sources); // only frees the merged references
+    }
 
     // XXX do we need to measure aperture photometry corrections?
@@ -157 +161 @@
     psMetadataAdd (fakeReadout->analysis, PS_LIST_TAIL, "PSPHOT.SOURCES", PS_DATA_ARRAY, "fake photometry ", fakeSources);
 
-    // psMetadataAdd (readout->analysis, PS_LIST_TAIL, "PPSIM.FORCE.SOURCES", PS_DATA_ARRAY, "forced photometry ", forceSources);
+    psMetadataAdd (readout->analysis, PS_LIST_TAIL, "PPSIM.FORCE.SOURCES", PS_DATA_ARRAY, "forced photometry ", forceSources);
     psFree (forceSources);
     return true;