Changeset 32287

Timestamp:

Sep 2, 2011, 9:05:35 AM (15 years ago)

Author:

eugene

Message:

measure inserted flux (to test validity of flux model)

File:

• branches/eam_branches/ipp-20110710/ppSim/src/ppSimInsertStars.c (modified) (5 diffs)

branches/eam_branches/ipp-20110710/ppSim/src/ppSimInsertStars.c

@@ -6 +6 @@
     PTR = NULL;
 
+float imageSum (psImage *image);
 
 bool ppSimInsertStars (pmReadout *readout, psImage *expCorr, psArray *stars, pmConfig *config) {
@@ -111 +112 @@
         float radius = model->modelRadius (model->params, nSigmaLim * roughNoise);
         radius = PS_MAX (radius, 1.0);
+        radius = 20.0;
 
         // construct a source, with model flux pixels set, based on the model
@@ -127 +129 @@
 
         float par8 = (model->params->n == 8) ? model->params->data.F32[7] : 0.0;
-        fprintf (outfile, "%8.3f %8.3f %10.2f  %2d  %7.3f %5.3f  %5.3f %5.3f %5.3f %5.3f\n", star->x, star->y, star->flux, 0, source->psfMag, source->psfMagErr, axes.major, axes.minor, axes.theta, par8);
+        float starFlux = 0.0;
 
         // if psfConvolve is TRUE, we will (elsewhere) convolve the image we a PSF
@@ -135 +137 @@
         } else {
             // insert the source flux in the image
+            float sum1 = imageSum(source->pixels);
             pmSourceAddWithOffset (source, PM_MODEL_OP_FULL, 0xff, dX, dY);
+            float sum2 = imageSum(source->pixels);
+            starFlux = sum2 - sum1;
             
             // insert the source flux in the noise image
@@ -149 +154 @@
             RESET(source->blends);
         }
+        fprintf (outfile, "%8.3f %8.3f %10.2f  %2d  %7.3f %5.3f  %5.3f %5.3f %5.3f %5.3f  :  %f %f\n", star->x, star->y, starFlux, 0, source->psfMag, source->psfMagErr, axes.major, axes.minor, axes.theta, par8, model->params->data.F32[PM_PAR_I0], radius);
 
         // add the sources to the source array