Changeset 30975 for branches/eam_branches/ipp-20110213/psphot/src/psphotSourceFits.c

Timestamp:

Mar 18, 2011, 2:10:35 PM (15 years ago)

Author:

eugene

Message:

we had 3 items conceptually equivalent to the brightness of a peak: "value", "flux", and "SN"; Ive modified pmPeak to carry explicitly named elements "rawFlux", "smoothFlux", and "detValue". Ive also added the corresponding variance values for rawFlux and smoothFlux. This lets us choose independently of the peak detection process whether flux comparisons are done relative to the smoothed or unsmoothed image. It also simplifies tests of the peak flux compared to something else

File:

• branches/eam_branches/ipp-20110213/psphot/src/psphotSourceFits.c (modified) (3 diffs)

branches/eam_branches/ipp-20110213/psphot/src/psphotSourceFits.c

@@ -51 +51 @@
     psArrayAdd (sourceSet, 16, source);
 
-    psTrace ("psphot", 4, "fitting blended source at %f %f : %f\n", source->peak->xf, source->peak->yf, source->peak->flux);
+    psTrace ("psphot", 4, "fitting blended source at %f %f : %f\n", source->peak->xf, source->peak->yf, source->peak->rawFlux);
 
     // we need to include all blends in the fit (unless primary is saturated?)
@@ -69 +69 @@
 
         // XXX assume local sky is 0.0?
-        model->params->data.F32[PM_PAR_I0] = blend->peak->flux;
+        model->params->data.F32[PM_PAR_I0] = blend->peak->rawFlux;
         model->params->data.F32[PM_PAR_XPOS] = blend->peak->xf;
         model->params->data.F32[PM_PAR_YPOS] = blend->peak->yf;
@@ -83 +83 @@
         psArrayAdd (sourceSet, 16, blend);
 
-        psTrace ("psphot", 5, "adding source at %f %f : %f\n", blend->peak->xf, blend->peak->yf, blend->peak->flux);
+        psTrace ("psphot", 5, "adding source at %f %f : %f\n", blend->peak->xf, blend->peak->yf, blend->peak->rawFlux);
 
         // free to avoid double counting model