Changeset 31451 for trunk/psModules/src/objects/pmSourceIO_PS1_CAL_0.c
- Timestamp:
- May 5, 2011, 11:02:53 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/src/objects/pmSourceIO_PS1_CAL_0.c
r31153 r31451 71 71 psF32 xPos, yPos; 72 72 psF32 xErr, yErr; 73 psF32 errMag, chisq;73 psF32 psfMagErr, chisq; 74 74 75 75 pmChip *chip = readout->parent->parent; … … 121 121 chisq = model->chisq; 122 122 123 // need to determine the PSF photometry error: source-> errMagis the error on the 'best' model mag.124 errMag= model->dparams->data.F32[PM_PAR_I0] / model->params->data.F32[PM_PAR_I0];123 // need to determine the PSF photometry error: source->psfMagErr is the error on the 'best' model mag. 124 psfMagErr = model->dparams->data.F32[PM_PAR_I0] / model->params->data.F32[PM_PAR_I0]; 125 125 } else { 126 126 xPos = source->peak->xf; … … 132 132 axes.theta = NAN; 133 133 chisq = NAN; 134 errMag= NAN;134 psfMagErr = NAN; 135 135 } 136 136 … … 167 167 psMetadataAdd (row, PS_LIST_TAIL, "Y_PSF_SIG", PS_DATA_F32, "Sigma in PSF y coordinate", yErr); 168 168 psMetadataAdd (row, PS_LIST_TAIL, "CAL_PSF_MAG", PS_DATA_F32, "Calibrated Magnitude from PSF Fit", calMag); 169 psMetadataAdd (row, PS_LIST_TAIL, "CAL_PSF_MAG_SIG", PS_DATA_F32, "Calibrated Magnitude Error", errMag);169 psMetadataAdd (row, PS_LIST_TAIL, "CAL_PSF_MAG_SIG", PS_DATA_F32, "Calibrated Magnitude Error", psfMagErr); 170 170 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG", PS_DATA_F32, "PSF fit instrumental magnitude", source->psfMag); 171 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", errMag);171 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", psfMagErr); 172 172 psMetadataAdd (row, PS_LIST_TAIL, "PEAK_FLUX_AS_MAG", PS_DATA_F32, "Peak flux expressed as magnitude", peakMag); 173 173 psMetadataAdd (row, PS_LIST_TAIL, "SKY", PS_DATA_F32, "Sky level", source->sky); … … 284 284 // XXX use these to determine PAR[PM_PAR_I0]? 285 285 source->psfMag = psMetadataLookupF32 (&status, row, "PSF_INST_MAG"); 286 source-> errMag= psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG");286 source->psfMagErr = psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG"); 287 287 PAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? pow(10.0, -0.4*source->psfMag) : NAN; 288 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source-> errMag: NAN;288 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source->psfMagErr : NAN; 289 289 290 290 pmPSF_AxesToModel (PAR, axes);
Note:
See TracChangeset
for help on using the changeset viewer.
