Changeset 12978 for branches/eam_02_branch/psphot/src/psphotRadiusChecks.c

Timestamp:

Apr 23, 2007, 4:54:30 PM (19 years ago)

Author:

eugene

Message:

moving to maskObj and maskView, re-wrote psphotEnsemblePSF using modelFlux, renamed to psphotFitSourcesLinear

File:

• branches/eam_02_branch/psphot/src/psphotRadiusChecks.c (modified) (5 diffs)

branches/eam_02_branch/psphot/src/psphotRadiusChecks.c

@@ -25 +25 @@
 bool psphotCheckRadiusPSF (pmReadout *readout, pmSource *source, pmModel *model)
 {
+    psF32 *PAR = model->params->data.F32;
+
+    // XXX do we have a better value for the sky noise level?  not really...
     pmMoments *moments = source->moments;
-    // do we have a better value for the sky noise level?
-    // not really...
 
     // set the fit radius based on the object flux limit and the model
@@ -46 +47 @@
     }
 
-    bool status = pmSourceRedefinePixels (source, readout, model->params->data.F32[2], model->params->data.F32[3], model->radiusFit);
+    bool status = pmSourceRedefinePixels (source, readout, PAR[PM_PAR_XPOS], PAR[PM_PAR_YPOS], model->radiusFit);
+
+    // set the mask to flag the excluded pixels
+    psImageKeepCircle (source->maskObj, PAR[PM_PAR_XPOS], PAR[PM_PAR_YPOS], model->radiusFit, "OR", PM_MASK_MARK);
     return status;
 }
 
 bool psphotCheckRadiusPSFBlend (pmReadout *readout, pmSource *source, pmModel *model, float dR) {
+
+    psF32 *PAR = model->params->data.F32;
 
     pmMoments *moments = source->moments;
@@ -63 +69 @@
     }
 
-    bool status = pmSourceRedefinePixels (source, readout, model->params->data.F32[2], model->params->data.F32[3], model->radiusFit);
+    bool status = pmSourceRedefinePixels (source, readout, PAR[PM_PAR_XPOS], PAR[PM_PAR_YPOS], model->radiusFit);
+
+    // set the mask to flag the excluded pixels
+    psImageKeepCircle (source->maskObj, PAR[PM_PAR_XPOS], PAR[PM_PAR_YPOS], model->radiusFit, "OR", PM_MASK_MARK);
     return status;
 }
@@ -86 +95 @@
 bool psphotCheckRadiusEXT (pmReadout *readout, pmSource *source, pmModel *model) {
 
+    psF32 *PAR = model->params->data.F32;
+
     pmMoments *moments = source->moments;
     if (moments == NULL) return false;
@@ -94 +105 @@
 
     // redefine the pixels if needed
-    bool status = pmSourceRedefinePixels (source, readout, model->params->data.F32[2], model->params->data.F32[3], model->radiusFit);
+    bool status = pmSourceRedefinePixels (source, readout, PAR[PM_PAR_XPOS], PAR[PM_PAR_YPOS], model->radiusFit);
+
+    // set the mask to flag the excluded pixels
+    psImageKeepCircle (source->maskObj, PAR[PM_PAR_XPOS], PAR[PM_PAR_YPOS], model->radiusFit, "OR", PM_MASK_MARK);
     return status;
 }