Changeset 9813 for trunk/psphot/src/psphotRadiusChecks.c
- Timestamp:
- Oct 31, 2006, 9:40:34 AM (20 years ago)
- File:
-
- 1 edited
-
trunk/psphot/src/psphotRadiusChecks.c (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotRadiusChecks.c
r9270 r9813 24 24 25 25 // set the fit radius based on the object flux limit and the model 26 model->radius TMP= (RADIUS_TYPE) (modelRadiusPSF (model->params, PSF_FIT_NSIGMA*moments->dSky) + PSF_FIT_PADDING);27 if (isnan(model->radius TMP)) psAbort ("apply_psf_model", "error in radius");26 model->radiusFit = (RADIUS_TYPE) (modelRadiusPSF (model->params, PSF_FIT_NSIGMA*moments->dSky) + PSF_FIT_PADDING); 27 if (isnan(model->radiusFit)) psAbort ("apply_psf_model", "error in radius"); 28 28 29 29 if (source->mode & PM_SOURCE_MODE_SATSTAR) { 30 model->radius TMP*= 2;30 model->radiusFit *= 2; 31 31 } 32 32 33 bool status = pmSourceRedefinePixels (source, readout, model->params->data.F32[2], model->params->data.F32[3], model->radius TMP);33 bool status = pmSourceRedefinePixels (source, readout, model->params->data.F32[2], model->params->data.F32[3], model->radiusFit); 34 34 return status; 35 35 } … … 41 41 42 42 // set the fit radius based on the object flux limit and the model 43 model->radius TMP= (RADIUS_TYPE) (modelRadiusPSF (model->params, PSF_FIT_NSIGMA*moments->dSky) + dR + PSF_FIT_PADDING);44 if (isnan(model->radius TMP)) psAbort ("apply_psf_model", "error in radius");43 model->radiusFit = (RADIUS_TYPE) (modelRadiusPSF (model->params, PSF_FIT_NSIGMA*moments->dSky) + dR + PSF_FIT_PADDING); 44 if (isnan(model->radiusFit)) psAbort ("apply_psf_model", "error in radius"); 45 45 46 46 if (source->mode & PM_SOURCE_MODE_SATSTAR) { 47 model->radius TMP*= 2;47 model->radiusFit *= 2; 48 48 } 49 49 50 bool status = pmSourceRedefinePixels (source, readout, model->params->data.F32[2], model->params->data.F32[3], model->radius TMP);50 bool status = pmSourceRedefinePixels (source, readout, model->params->data.F32[2], model->params->data.F32[3], model->radiusFit); 51 51 return status; 52 52 } … … 74 74 75 75 // set the fit radius based on the object flux limit and the model 76 model->radius TMP= (RADIUS_TYPE) (modelRadiusEXT (model->params, EXT_FIT_NSIGMA*moments->dSky) + EXT_FIT_PADDING);77 if (isnan(model->radius TMP)) psAbort (__func__, "error in radius");76 model->radiusFit = (RADIUS_TYPE) (modelRadiusEXT (model->params, EXT_FIT_NSIGMA*moments->dSky) + EXT_FIT_PADDING); 77 if (isnan(model->radiusFit)) psAbort (__func__, "error in radius"); 78 78 79 79 // redefine the pixels if needed 80 bool status = pmSourceRedefinePixels (source, readout, model->params->data.F32[2], model->params->data.F32[3], model->radius TMP);80 bool status = pmSourceRedefinePixels (source, readout, model->params->data.F32[2], model->params->data.F32[3], model->radiusFit); 81 81 return status; 82 82 }
Note:
See TracChangeset
for help on using the changeset viewer.
