Changeset 31451 for trunk/psModules/src/objects/pmSourceIO_PS1_DEV_1.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_DEV_1.c
r31153 r31451 70 70 psF32 xPos, yPos; 71 71 psF32 xErr, yErr; 72 psF32 errMag, chisq;72 psF32 psfMagErr, chisq; 73 73 74 74 // let's write these out in S/N order … … 103 103 chisq = model->chisq; 104 104 105 // need to determine the PSF photometry error: source-> errMagis the error on the 'best' model mag.106 errMag= model->dparams->data.F32[PM_PAR_I0] / model->params->data.F32[PM_PAR_I0];105 // need to determine the PSF photometry error: source->psfMagErr is the error on the 'best' model mag. 106 psfMagErr = model->dparams->data.F32[PM_PAR_I0] / model->params->data.F32[PM_PAR_I0]; 107 107 } else { 108 108 xPos = source->peak->xf; … … 114 114 axes.theta = NAN; 115 115 chisq = NAN; 116 errMag= NAN;116 psfMagErr = NAN; 117 117 } 118 118 … … 128 128 psMetadataAdd (row, PS_LIST_TAIL, "Y_PSF_SIG", PS_DATA_F32, "Sigma in PSF y coordinate", yErr); 129 129 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG", PS_DATA_F32, "PSF fit instrumental magnitude", source->psfMag); 130 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", errMag);130 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", psfMagErr); 131 131 psMetadataAdd (row, PS_LIST_TAIL, "PEAK_FLUX_AS_MAG", PS_DATA_F32, "Peak flux expressed as magnitude", peakMag); 132 132 psMetadataAdd (row, PS_LIST_TAIL, "SKY", PS_DATA_F32, "Sky level", source->sky); … … 253 253 // XXX use these to determine PAR[PM_PAR_I0]? 254 254 source->psfMag = psMetadataLookupF32 (&status, row, "PSF_INST_MAG"); 255 source-> errMag= psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG");255 source->psfMagErr = psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG"); 256 256 PAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? pow(10.0, -0.4*source->psfMag) : NAN; 257 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source-> errMag: NAN;257 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source->psfMagErr : NAN; 258 258 259 259 pmPSF_AxesToModel (PAR, axes);
Note:
See TracChangeset
for help on using the changeset viewer.
