- Timestamp:
- Apr 24, 2011, 9:47:48 PM (15 years ago)
- Location:
- branches/eam_branches/ipp-20110404/psModules/src/objects
- Files:
-
- 16 edited
-
pmPSFtryFitPSF.c (modified) (2 diffs)
-
pmSource.c (modified) (2 diffs)
-
pmSourceIO_CMF_PS1_DV1.c (modified) (2 diffs)
-
pmSourceIO_CMF_PS1_DV2.c (modified) (2 diffs)
-
pmSourceIO_CMF_PS1_SV1.c (modified) (2 diffs)
-
pmSourceIO_CMF_PS1_V1.c (modified) (2 diffs)
-
pmSourceIO_CMF_PS1_V2.c (modified) (2 diffs)
-
pmSourceIO_CMP.c (modified) (3 diffs)
-
pmSourceIO_PS1_CAL_0.c (modified) (5 diffs)
-
pmSourceIO_PS1_DEV_0.c (modified) (2 diffs)
-
pmSourceIO_PS1_DEV_1.c (modified) (5 diffs)
-
pmSourceIO_RAW.c (modified) (2 diffs)
-
pmSourceIO_SMPDATA.c (modified) (2 diffs)
-
pmSourceIO_SX.c (modified) (1 diff)
-
pmSourceMatch.c (modified) (2 diffs)
-
pmSourceOutputs.c (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/ipp-20110404/psModules/src/objects/pmPSFtryFitPSF.c
r31153 r31360 122 122 psfTry->fitMag->data.F32[i] = source->psfMag; 123 123 psfTry->metric->data.F32[i] = source->apMag - source->psfMag; 124 psfTry->metricErr->data.F32[i] = source-> errMag;124 psfTry->metricErr->data.F32[i] = source->psfMagErr; 125 125 126 126 // XXX this did not work: modifies shape of psf too much … … 131 131 source->peak->xf, source->peak->yf, 132 132 source->modelPSF->params->data.F32[PM_PAR_XPOS], source->modelPSF->params->data.F32[PM_PAR_YPOS], 133 source->psfMag, source->apMag, source-> errMag,133 source->psfMag, source->apMag, source->psfMagErr, 134 134 source->modelPSF->params->data.F32[PM_PAR_I0], 135 135 source->modelPSF->params->data.F32[PM_PAR_SXX], source->modelPSF->params->data.F32[PM_PAR_SXY], -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSource.c
r31312 r31360 137 137 source->psfFluxErr = NAN; 138 138 source->extMag = NAN; 139 source-> errMag= NAN;139 source->psfMagErr = NAN; 140 140 source->apMag = NAN; 141 141 source->apMagRaw = NAN; … … 207 207 source->psfFluxErr = in->psfFluxErr; 208 208 source->extMag = in->extMag; 209 source-> errMag = in->errMag;209 source->psfMagErr = in->psfMagErr; 210 210 source->apMag = in->apMag; 211 211 source->apMagRaw = in->apMagRaw; -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_CMF_PS1_DV1.c
r31153 r31360 125 125 psMetadataAdd (row, PS_LIST_TAIL, "PLTSCALE", PS_DATA_F32, "plate scale at source (arcsec/pixel)", outputs.pltScale); 126 126 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG", PS_DATA_F32, "PSF fit instrumental magnitude", source->psfMag); 127 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", source-> errMag);127 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", source->psfMagErr); 128 128 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_FLUX", PS_DATA_F32, "PSF fit instrumental magnitude", source->psfFlux); 129 129 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_FLUX_SIG",PS_DATA_F32, "Sigma of PSF instrumental magnitude", source->psfFluxErr); … … 256 256 257 257 source->psfMag = psMetadataLookupF32 (&status, row, "PSF_INST_MAG"); 258 source-> errMag= psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG");258 source->psfMagErr = psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG"); 259 259 source->apMag = psMetadataLookupF32 (&status, row, "AP_MAG"); 260 260 261 261 // XXX this scaling is incorrect: does not include the 2 \pi AREA factor 262 262 PAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? pow(10.0, -0.4*source->psfMag) : NAN; 263 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source-> errMag: NAN;263 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source->psfMagErr : NAN; 264 264 265 265 pmPSF_AxesToModel (PAR, axes); -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_CMF_PS1_DV2.c
r31153 r31360 124 124 psMetadataAdd (row, PS_LIST_TAIL, "PLTSCALE", PS_DATA_F32, "plate scale at source (arcsec/pixel)", outputs.pltScale); 125 125 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG", PS_DATA_F32, "PSF fit instrumental magnitude", source->psfMag); 126 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", source-> errMag);126 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", source->psfMagErr); 127 127 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_FLUX", PS_DATA_F32, "PSF fit instrumental magnitude", source->psfFlux); 128 128 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_FLUX_SIG",PS_DATA_F32, "Sigma of PSF instrumental magnitude", source->psfFluxErr); … … 274 274 275 275 source->psfMag = psMetadataLookupF32 (&status, row, "PSF_INST_MAG"); 276 source-> errMag= psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG");276 source->psfMagErr = psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG"); 277 277 source->apMag = psMetadataLookupF32 (&status, row, "AP_MAG"); 278 278 279 279 // XXX this scaling is incorrect: does not include the 2 \pi AREA factor 280 280 PAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? pow(10.0, -0.4*source->psfMag) : NAN; 281 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source-> errMag: NAN;281 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source->psfMagErr : NAN; 282 282 283 283 pmPSF_AxesToModel (PAR, axes); -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_CMF_PS1_SV1.c
r31153 r31360 125 125 psMetadataAdd (row, PS_LIST_TAIL, "PLTSCALE", PS_DATA_F32, "plate scale at source (arcsec/pixel)", outputs.pltScale); 126 126 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG", PS_DATA_F32, "PSF fit instrumental magnitude", source->psfMag); 127 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", source-> errMag);127 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", source->psfMagErr); 128 128 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_FLUX", PS_DATA_F32, "PSF fit instrumental flux (counts)", source->psfFlux); 129 129 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_FLUX_SIG",PS_DATA_F32, "Sigma of PSF instrumental flux", source->psfFluxErr); … … 273 273 // XXX use these to determine PAR[PM_PAR_I0]? 274 274 source->psfMag = psMetadataLookupF32 (&status, row, "PSF_INST_MAG"); 275 source-> errMag= psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG");275 source->psfMagErr = psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG"); 276 276 source->apMag = psMetadataLookupF32 (&status, row, "AP_MAG"); 277 277 278 278 // XXX this scaling is incorrect: does not include the 2 \pi AREA factor 279 279 PAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? pow(10.0, -0.4*source->psfMag) : NAN; 280 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source-> errMag: NAN;280 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source->psfMagErr : NAN; 281 281 282 282 pmPSF_AxesToModel (PAR, axes); -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_CMF_PS1_V1.c
r31153 r31360 118 118 psMetadataAdd (row, PS_LIST_TAIL, "PLTSCALE", PS_DATA_F32, "plate scale at source (arcsec/pixel)", outputs.pltScale); 119 119 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG", PS_DATA_F32, "PSF fit instrumental magnitude", source->psfMag); 120 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", source-> errMag);120 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", source->psfMagErr); 121 121 psMetadataAdd (row, PS_LIST_TAIL, "AP_MAG", PS_DATA_F32, "magnitude in standard aperture", source->apMag); 122 122 psMetadataAdd (row, PS_LIST_TAIL, "AP_MAG_RADIUS", PS_DATA_F32, "radius used for aperture mags", outputs.apRadius); … … 241 241 // XXX use these to determine PAR[PM_PAR_I0]? 242 242 source->psfMag = psMetadataLookupF32 (&status, row, "PSF_INST_MAG"); 243 source-> errMag= psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG");243 source->psfMagErr = psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG"); 244 244 source->apMag = psMetadataLookupF32 (&status, row, "AP_MAG"); 245 245 246 246 // XXX this scaling is incorrect: does not include the 2 \pi AREA factor 247 247 PAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? pow(10.0, -0.4*source->psfMag) : NAN; 248 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source-> errMag: NAN;248 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source->psfMagErr : NAN; 249 249 250 250 pmPSF_AxesToModel (PAR, axes); -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_CMF_PS1_V2.c
r31153 r31360 116 116 psMetadataAdd (row, PS_LIST_TAIL, "PLTSCALE", PS_DATA_F32, "plate scale at source (arcsec/pixel)", outputs.pltScale); 117 117 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG", PS_DATA_F32, "PSF fit instrumental magnitude", source->psfMag); 118 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", source-> errMag);118 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", source->psfMagErr); 119 119 psMetadataAdd (row, PS_LIST_TAIL, "AP_MAG", PS_DATA_F32, "magnitude in standard aperture", source->apMag); 120 120 psMetadataAdd (row, PS_LIST_TAIL, "AP_MAG_RADIUS", PS_DATA_F32, "radius used for aperture mags", outputs.apRadius); … … 247 247 // XXX use these to determine PAR[PM_PAR_I0]? 248 248 source->psfMag = psMetadataLookupF32 (&status, row, "PSF_INST_MAG"); 249 source-> errMag= psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG");249 source->psfMagErr = psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG"); 250 250 source->apMag = psMetadataLookupF32 (&status, row, "AP_MAG"); 251 251 252 252 // XXX this scaling is incorrect: does not include the 2 \pi AREA factor 253 253 PAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? pow(10.0, -0.4*source->psfMag) : NAN; 254 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source-> errMag: NAN;254 dPAR[PM_PAR_I0] = (isfinite(source->psfMag)) ? PAR[PM_PAR_I0] * source->psfMagErr : NAN; 255 255 256 256 pmPSF_AxesToModel (PAR, axes); -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_CMP.c
r29004 r31360 137 137 axes = pmPSF_ModelToAxes (PAR, 20.0); 138 138 139 float errMag = isfinite(source->errMag) ? source->errMag: 999;139 float psfMagErr = isfinite(source->psfMagErr) ? source->psfMagErr : 999; 140 140 141 141 psLineInit (line); … … 143 143 psLineAdd (line, "%6.1f ", PAR[PM_PAR_YPOS]); 144 144 psLineAdd (line, "%6.3f ", PS_MIN (99.0, source->psfMag + ZERO_POINT)); 145 psLineAdd (line, "%03d ", PS_MIN (999, (int) errMag));145 psLineAdd (line, "%03d ", PS_MIN (999, (int)psfMagErr)); 146 146 psLineAdd (line, "%2d ", type); 147 147 psLineAdd (line, "%3.1f ", lsky); … … 280 280 source->psfMag = atof (array->data[2]); 281 281 source->extMag = atof (array->data[6]); 282 source-> errMag= atof (array->data[3]) / 1000.0;282 source->psfMagErr = atof (array->data[3]) / 1000.0; 283 283 source->apMag = atof (array->data[7]); 284 284 axes.major = atof (array->data[8]); -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_PS1_CAL_0.c
r31153 r31360 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); -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_PS1_DEV_0.c
r31153 r31360 113 113 psMetadataAdd (row, PS_LIST_TAIL, "Y_PSF_SIG", PS_DATA_F32, "Sigma in PSF y coordinate", yErr); 114 114 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG", PS_DATA_F32, "PSF fit instrumental magnitude", PS_MIN (99.0, source->psfMag)); 115 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", PS_MIN (99.0, source-> errMag));115 psMetadataAdd (row, PS_LIST_TAIL, "PSF_INST_MAG_SIG", PS_DATA_F32, "Sigma of PSF instrumental magnitude", PS_MIN (99.0, source->psfMagErr)); 116 116 psMetadataAdd (row, PS_LIST_TAIL, "PEAK_FLUX_AS_MAG", PS_DATA_F32, "Peak flux expressed as magnitude", PS_MIN (99.0, peakMag)); 117 117 psMetadataAdd (row, PS_LIST_TAIL, "SKY", PS_DATA_F32, "Sky level", source->sky); … … 212 212 // XXX use these to determine PAR[PM_PAR_I0]? 213 213 source->psfMag = psMetadataLookupF32 (&status, row, "PSF_INST_MAG"); 214 source-> errMag= psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG");214 source->psfMagErr = psMetadataLookupF32 (&status, row, "PSF_INST_MAG_SIG"); 215 215 216 216 pmPSF_AxesToModel (PAR, axes); -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_PS1_DEV_1.c
r31153 r31360 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); -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_RAW.c
r31153 r31360 110 110 fprintf (f, "%7.1f %7.1f %7.1f %8.4f %7.4f %7.4f ", 111 111 PAR[PM_PAR_XPOS], PAR[PM_PAR_YPOS], source->sky, 112 source->psfMag, source-> errMag, dPos);112 source->psfMag, source->psfMagErr, dPos); 113 113 114 114 for (j = 4; j < model->params->n; j++) { … … 174 174 fprintf (f, "%7.1f %7.1f %7.1f %8.4f %7.4f %7.4f ", 175 175 PAR[PM_PAR_XPOS], PAR[PM_PAR_YPOS], source->sky, 176 source->extMag, source-> errMag, dPos);176 source->extMag, source->psfMagErr, dPos); 177 177 178 178 for (j = 4; j < model->params->n; j++) { -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_SMPDATA.c
r31153 r31360 106 106 psMetadataAdd (row, PS_LIST_TAIL, "Y_PIX", PS_DATA_F32, "", yPos); 107 107 psMetadataAdd (row, PS_LIST_TAIL, "MAG_RAW", PS_DATA_F32, "", PS_MIN (99.0, source->psfMag + ZERO_POINT)); 108 psMetadataAdd (row, PS_LIST_TAIL, "MAG_ERR", PS_DATA_F32, "", PS_MIN (999, 1000*source-> errMag));108 psMetadataAdd (row, PS_LIST_TAIL, "MAG_ERR", PS_DATA_F32, "", PS_MIN (999, 1000*source->psfMagErr)); 109 109 psMetadataAdd (row, PS_LIST_TAIL, "MAG_GAL", PS_DATA_F32, "", PS_MIN (99.0, source->extMag + ZERO_POINT)); 110 110 psMetadataAdd (row, PS_LIST_TAIL, "MAG_AP", PS_DATA_F32, "", PS_MIN (99.0, source->apMag + ZERO_POINT)); … … 194 194 source->psfMag = psMetadataLookupF32 (&status, row, "MAG_RAW") - ZERO_POINT; 195 195 source->extMag = psMetadataLookupF32 (&status, row, "MAG_GAL") - ZERO_POINT; 196 source-> errMag= psMetadataLookupF32 (&status, row, "MAG_ERR") * 0.001;196 source->psfMagErr = psMetadataLookupF32 (&status, row, "MAG_ERR") * 0.001; 197 197 source->apMag = psMetadataLookupF32 (&status, row, "MAG_AP") - ZERO_POINT; 198 198 -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceIO_SX.c
r29004 r31360 88 88 psLineAdd (line, "%11.3f", PAR[PM_PAR_YPOS]); 89 89 psLineAdd (line, "%9.4f", source->psfMag); 90 psLineAdd (line, "%9.4f", source-> errMag);90 psLineAdd (line, "%9.4f", source->psfMagErr); 91 91 psLineAdd (line, "%13.4f", source->sky); 92 92 psLineAdd (line, "%9.2f", axes.major); -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceMatch.c
r29004 r31360 88 88 pmSource *source = sources->data[i]; // Source of interest 89 89 if (!source || (source->mode & SOURCE_MASK) || !isfinite(source->psfMag) || 90 !isfinite(source-> errMag) || source->psfMag > SOURCE_FAINTEST) {90 !isfinite(source->psfMagErr) || source->psfMag > SOURCE_FAINTEST) { 91 91 continue; 92 92 } … … 102 102 (*y)->data.F32[num] = ySrc; 103 103 (*mag)->data.F32[num] = source->psfMag; 104 (*magErr)->data.F32[num] = source-> errMag;104 (*magErr)->data.F32[num] = source->psfMagErr; 105 105 (*indices)->data.S32[num] = i; 106 106 num++; -
branches/eam_branches/ipp-20110404/psModules/src/objects/pmSourceOutputs.c
r31312 r31360 103 103 outputs->yErr = dPAR[PM_PAR_YPOS]; 104 104 } else { 105 outputs->xErr = fwhmMajor * source-> errMag/ 2.35;106 outputs->yErr = fwhmMinor * source-> errMag/ 2.35;105 outputs->xErr = fwhmMajor * source->psfMagErr / 2.35; 106 outputs->yErr = fwhmMinor * source->psfMagErr / 2.35; 107 107 } 108 108 if (isfinite(PAR[PM_PAR_SXX]) && isfinite(PAR[PM_PAR_SXY]) && isfinite(PAR[PM_PAR_SYY])) { … … 126 126 outputs->xPos = source->moments->Mx; 127 127 outputs->yPos = source->moments->My; 128 outputs->xErr = fwhmMajor * source-> errMag/ 2.35;129 outputs->yErr = fwhmMinor * source-> errMag/ 2.35;128 outputs->xErr = fwhmMajor * source->psfMagErr / 2.35; 129 outputs->yErr = fwhmMinor * source->psfMagErr / 2.35; 130 130 } else { 131 131 outputs->xPos = source->peak->xf;
Note:
See TracChangeset
for help on using the changeset viewer.
