Changeset 31363
- Timestamp:
- Apr 24, 2011, 11:06:23 PM (15 years ago)
- Location:
- branches/eam_branches/ipp-20110404
- Files:
-
- 7 edited
-
ppSim/src/ppSimInsertGalaxies.c (modified) (2 diffs)
-
ppSim/src/ppSimInsertStars.c (modified) (1 diff)
-
ppSub/src/ppSub.c (modified) (1 diff)
-
ppSub/src/ppSubFlagNeighbors.c (modified) (2 diffs)
-
ppSub/src/ppSubReadoutPhotometry.c (modified) (1 diff)
-
psastro/src/psastroConvert.c (modified) (2 diffs)
-
pswarp/src/pswarpTransformSources.c (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/eam_branches/ipp-20110404/ppSim/src/ppSimInsertGalaxies.c
r29002 r31363 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); -
branches/eam_branches/ipp-20110404/ppSim/src/ppSimInsertStars.c
r29002 r31363 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 -
branches/eam_branches/ipp-20110404/ppSub/src/ppSub.c
r30874 r31363 94 94 95 95 if (data) { 96 psString dump_file = psMetadataLookupStr(NULL, data->config->arguments, " -dumpconfig");96 psString dump_file = psMetadataLookupStr(NULL, data->config->arguments, "DUMP_CONFIG"); 97 97 if (dump_file) { 98 98 if (!pmConfigDump(data->config, dump_file)) { -
branches/eam_branches/ipp-20110404/ppSub/src/ppSubFlagNeighbors.c
r31156 r31363 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; -
branches/eam_branches/ipp-20110404/ppSub/src/ppSubReadoutPhotometry.c
r30619 r31363 175 175 176 176 fprintf(sourceFile, "%f %f %f %f %f %f %f %f %d %f %f %f %f %f\n", 177 x, y, source->psfMag, source-> errMag, chi2, source->crNsigma, source->extNsigma,177 x, y, source->psfMag, source->psfMagErr, chi2, source->crNsigma, source->extNsigma, 178 178 source->pixWeight, source->mode, xMoment, yMoment, xxMoment, xyMoment, yyMoment); 179 179 } -
branches/eam_branches/ipp-20110404/psastro/src/psastroConvert.c
r31161 r31363 102 102 pmSource *source = sources->data[n]; 103 103 104 psTrace ("psastro", 6, "mag: %f +/- %f, mode: %x, skip: %x\n", source->psfMag, source-> errMag, source->mode, skip);104 psTrace ("psastro", 6, "mag: %f +/- %f, mode: %x, skip: %x\n", source->psfMag, source->psfMagErr, source->mode, skip); 105 105 106 106 if (source->mode & skip) { … … 166 166 obj->pix->yErr = dPAR[PM_PAR_YPOS]; 167 167 obj->Mag = source->psfMag; 168 obj->dMag = source-> errMag;168 obj->dMag = source->psfMagErr; 169 169 170 170 // XXX do we have the information giving the readout and cell offset? -
branches/eam_branches/ipp-20110404/pswarp/src/pswarpTransformSources.c
r31159 r31363 101 101 new->psfMag = source->psfMag; 102 102 new->extMag = source->extMag; 103 new-> errMag = source->errMag;103 new->psfMagErr = source->psfMagErr; 104 104 new->apMag = source->apMag; 105 105 new->pixWeightNotBad = source->pixWeightNotBad; … … 116 116 (isfinite(new->psfMag) ? pow(10.0, -0.4*new->psfMag) : NAN); 117 117 new->modelPSF->dparams->data.F32[PM_PAR_I0]= 118 (isfinite(new->psfMag) ? new-> errMag*pow(10.0, -0.4*new->psfMag) : NAN);118 (isfinite(new->psfMag) ? new->psfMagErr*pow(10.0, -0.4*new->psfMag) : NAN); 119 119 120 120 #if 0
Note:
See TracChangeset
for help on using the changeset viewer.
