Changeset 14614 for branches/eam_branch_20070817/ppSim/src/ppSimInsertSources.c

Timestamp:

Aug 22, 2007, 3:05:47 PM (19 years ago)

Author:

eugene

Message:

load / apply psf model from external source

File:

• branches/eam_branch_20070817/ppSim/src/ppSimInsertSources.c (modified) (4 diffs)

branches/eam_branch_20070817/ppSim/src/ppSimInsertSources.c

@@ -3 +3 @@
 bool ppSimInsertSources (pmReadout *readout, psImage *expCorr, psArray *stars, pmConfig *config) {
 
+    bool mdok;
 
     assert (stars);
+    assert (readout);
     
     pmCell *cell = readout->parent;
     pmChip *chip = cell->parent;
 
-    // XXX we probably need this to get the size right
-# if (0)
-    bool mdok;
+    // XXX this is an estimate of the sky noise based on the inputs to the image simulation.
+    // XXX update this to allow the estimate based on the measured sky background
+    // XXX this is missing the gain.
     psMetadata *recipe = psMetadataLookupMetadata(NULL, config->recipes, PPSIM_RECIPE); // Recipe
 
@@ -29 +31 @@
     // Rough noise estimate, appropriate for entire cell (use for source radius?)
     float roughNoise = sqrtf(PS_SQR(readnoise) + (darkRate + skyRate) * expTime);
-# endif
 
     int x0Chip = psMetadataLookupS32(NULL, chip->concepts, "CHIP.X0");
@@ -43 +44 @@
     int binning = psMetadataLookupS32(NULL, config->arguments, "BINNING"); // Binning in x and y
 
-    pmPSF *psf = psMetadataLookupPtr (NULL, chip->analysis, "PSF");
+    pmPSF *psf = psMetadataLookupPtr (&mdok, chip->analysis, "PSPHOT.PSF");
+    assert (psf);
 
     int dX = PM_CELL_TO_CHIP (0.0, x0Cell, xParityCell, binning);
@@ -74 +76 @@
         pmModelSetFlux (model, star->flux);
 
+        // XXX let the flux limit be a user-defined number of sky sigmas (not just 1.0)
+        float radius = model->modelRadius (model->params, roughNoise);
+        radius = PS_MAX (radius, 1.0);
+
         // construct a source, with model flux pixels set, based on the model
-        pmSource *source = pmSourceFromModel (model, readout, PM_SOURCE_TYPE_STAR);
+        pmSource *source = pmSourceFromModel (model, readout, radius, PM_SOURCE_TYPE_STAR);
 
         // XXX add the sources to a source array
 
+        // insert the source flux in the image
         pmSourceAddWithOffset (source, PM_MODEL_OP_FULL, 0xff, dX, dY);
     }