Changeset 31449
- Timestamp:
- May 5, 2011, 10:40:58 AM (15 years ago)
- Location:
- trunk
- Files:
-
- 4 edited
-
ppSim/src/ppSimInsertGalaxies.c (modified) (2 diffs)
-
ppSim/src/ppSimInsertStars.c (modified) (1 diff)
-
ppSub/src/ppSubFlagNeighbors.c (modified) (2 diffs)
-
ppSub/src/ppSubReadoutPhotometry.c (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/ppSim/src/ppSimInsertGalaxies.c
r29002 r31449 162 162 source->psfMag = -2.5*log10(galaxy->flux); 163 163 source->extMag = -2.5*log10(galaxy->flux); 164 source-> errMag= sqrt(Area*PS_SQR(roughNoise) + galaxy->flux) / galaxy->flux;164 source->psfMagErr = sqrt(Area*PS_SQR(roughNoise) + galaxy->flux) / galaxy->flux; 165 165 166 166 // insert the source flux in the image … … 173 173 174 174 float par8 = (model->params->n == 8) ? model->params->data.F32[7] : 0.0; 175 fprintf (outfile, "%8.3f %8.3f %10.2f %2d %7.3f %5.3f %5.3f %5.3f %5.3f %5.3f\n", galaxy->x, galaxy->y, flux, 1, source->psfMag, source-> errMag, axes.major, axes.minor, axes.theta, par8);175 fprintf (outfile, "%8.3f %8.3f %10.2f %2d %7.3f %5.3f %5.3f %5.3f %5.3f %5.3f\n", galaxy->x, galaxy->y, flux, 1, source->psfMag, source->psfMagErr, axes.major, axes.minor, axes.theta, par8); 176 176 177 177 psArrayAdd (sources, 100,source); -
trunk/ppSim/src/ppSimInsertStars.c
r29002 r31449 121 121 122 122 source->psfMag = -2.5*log10(star->flux); 123 source-> errMag= sqrt(Area*PS_SQR(roughNoise) + flux) / flux;123 source->psfMagErr = sqrt(Area*PS_SQR(roughNoise) + flux) / flux; 124 124 125 125 // set the expected model errors 126 model->dparams->data.F32[PM_PAR_I0] = source-> errMag* model->params->data.F32[PM_PAR_I0];126 model->dparams->data.F32[PM_PAR_I0] = source->psfMagErr * model->params->data.F32[PM_PAR_I0]; 127 127 128 128 float par8 = (model->params->n == 8) ? model->params->data.F32[7] : 0.0; 129 fprintf (outfile, "%8.3f %8.3f %10.2f %2d %7.3f %5.3f %5.3f %5.3f %5.3f %5.3f\n", star->x, star->y, star->flux, 0, source->psfMag, source-> errMag, axes.major, axes.minor, axes.theta, par8);129 fprintf (outfile, "%8.3f %8.3f %10.2f %2d %7.3f %5.3f %5.3f %5.3f %5.3f %5.3f\n", star->x, star->y, star->flux, 0, source->psfMag, source->psfMagErr, axes.major, axes.minor, axes.theta, par8); 130 130 131 131 // if psfConvolve is TRUE, we will (elsewhere) convolve the image we a PSF -
trunk/ppSub/src/ppSubFlagNeighbors.c
r31156 r31449 112 112 bool positive = !matchRef && (sourceM1->imageID == 1); 113 113 positive |= matchRef && (sourceM1->imageID == 2); 114 float SN1 = isfinite(sourceM1-> errMag) ? 1.0 / sourceM1->errMag: NAN;114 float SN1 = isfinite(sourceM1->psfMagErr) ? 1.0 / sourceM1->psfMagErr : NAN; 115 115 if (positive) { 116 116 source->diffStats->SNp = SN1; … … 138 138 bool positive = !matchRef && (sourceM1->imageID == 1); 139 139 positive |= matchRef && (sourceM1->imageID == 2); 140 float SN1 = isfinite(sourceM1-> errMag) ? 1.0 / sourceM1->errMag: NAN;141 float SN2 = isfinite(sourceM2-> errMag) ? 1.0 / sourceM2->errMag: NAN;140 float SN1 = isfinite(sourceM1->psfMagErr) ? 1.0 / sourceM1->psfMagErr : NAN; 141 float SN2 = isfinite(sourceM2->psfMagErr) ? 1.0 / sourceM2->psfMagErr : NAN; 142 142 if (positive) { 143 143 source->diffStats->SNp = SN1; -
trunk/ppSub/src/ppSubReadoutPhotometry.c
r31435 r31449 176 176 177 177 fprintf(sourceFile, "%f %f %f %f %f %f %f %f %d %f %f %f %f %f\n", 178 x, y, source->psfMag, source-> errMag, chi2, source->crNsigma, source->extNsigma,178 x, y, source->psfMag, source->psfMagErr, chi2, source->crNsigma, source->extNsigma, 179 179 source->pixWeight, source->mode, xMoment, yMoment, xxMoment, xyMoment, yyMoment); 180 180 }
Note:
See TracChangeset
for help on using the changeset viewer.
