Changeset 9576
- Timestamp:
- Oct 13, 2006, 5:13:31 PM (20 years ago)
- Location:
- trunk/psphot/src
- Files:
-
- 7 edited
-
psphotArguments.c (modified) (1 diff)
-
psphotBlendFit.c (modified) (3 diffs)
-
psphotEnsemblePSF.c (modified) (1 diff)
-
psphotParseCamera.c (modified) (1 diff)
-
psphotReadout.c (modified) (2 diffs)
-
psphotReplaceUnfit.c (modified) (1 diff)
-
psphotSourceFits.c (modified) (10 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psphot/src/psphotArguments.c
r9270 r9576 132 132 psMetadataAddStr (config->arguments, PS_LIST_TAIL, "OUTPUT", 0, "", argv[1]); 133 133 134 psTrace( __func__, 1, "Done with psphotArguments...\n");134 psTrace("psphot", 1, "Done with psphotArguments...\n"); 135 135 return (config); 136 136 } -
trunk/psphot/src/psphotBlendFit.c
r9270 r9576 52 52 // skip sources which are insignificant flux? 53 53 if (source->modelPSF->params->data.F32[1] < 0.1) { 54 psTrace ("psphot BlendFit", 5, "skipping near-zero source: %f, %f : %f\n",54 psTrace ("psphot", 5, "skipping near-zero source: %f, %f : %f\n", 55 55 source->modelPSF->params->data.F32[1], 56 56 source->modelPSF->params->data.F32[2], … … 65 65 // try fitting PSFs, then try extended sources 66 66 if (psphotFitBlend (readout, source, psf)) { 67 psTrace ("psphot BlendFit", 5, "source at %7.1f, %7.1f is psf", source->moments->x, source->moments->y);67 psTrace ("psphot", 5, "source at %7.1f, %7.1f is psf", source->moments->x, source->moments->y); 68 68 Npsf ++; 69 69 continue; 70 70 } 71 71 if (psphotFitBlob (readout, source, sources, psf)) { 72 psTrace ("psphot BlendFit", 5, "source at %7.1f, %7.1f is ext", source->moments->x, source->moments->y);72 psTrace ("psphot", 5, "source at %7.1f, %7.1f is ext", source->moments->x, source->moments->y); 73 73 Next ++; 74 74 continue; 75 75 } 76 76 77 psTrace ("psphot BlendFit", 5, "source at %7.1f, %7.1f failed", source->moments->x, source->moments->y);77 psTrace ("psphot", 5, "source at %7.1f, %7.1f failed", source->moments->x, source->moments->y); 78 78 Nfail ++; 79 79 … … 84 84 } 85 85 86 psLogMsg ("psphot ", 3, "fit models: %f sec for %d objects (%d psf, %d ext, %d failed, %ld skipped)\n", psTimerMark ("psphot"), Nfit, Npsf, Next, Nfail, sources->n - Nfit);86 psLogMsg ("psphot.psphotBlendFit", 3, "fit models: %f sec for %d objects (%d psf, %d ext, %d failed, %ld skipped)\n", psTimerMark ("psphot"), Nfit, Npsf, Next, Nfail, sources->n - Nfit); 87 87 return (true); 88 88 } -
trunk/psphot/src/psphotEnsemblePSF.c
r9529 r9576 86 86 psPlane min = psImageBicubeMin (bicube); 87 87 88 psTrace ("psphot EnsemblePSF", 5, "peak coord: %f %f -> %f %f\n",88 psTrace ("psphot", 5, "peak coord: %f %f -> %f %f\n", 89 89 modelEXT->params->data.F32[PM_PAR_XPOS], modelEXT->params->data.F32[PM_PAR_YPOS], min.x + ix, min.y + iy); 90 90 -
trunk/psphot/src/psphotParseCamera.c
r9270 r9576 83 83 psFree (chips); 84 84 85 psTrace( __func__, 1, "Done with psphotParseCamera...\n");85 psTrace("psphot", 1, "Done with psphotParseCamera...\n"); 86 86 87 87 psErrorClear(); // some metadata lookup may have failed -
trunk/psphot/src/psphotReadout.c
r9565 r9576 24 24 psphotModelTest (readout, recipe); 25 25 26 if (psTraceGetLevel("psphot .images") > 5) {26 if (psTraceGetLevel("psphot") > 5) { 27 27 psphotSaveImage (NULL, readout->image, "image.fits"); 28 28 psphotSaveImage (NULL, readout->mask, "mask.fits"); … … 95 95 // measure aperture photometry corrections 96 96 if (!psphotApResid (readout, sources, recipe, psf)) { 97 psTrace ("psphot Readout", 4, "failure on psphotApResid");97 psTrace ("psphot", 4, "failure on psphotApResid"); 98 98 psError(PSPHOT_ERR_PHOTOM, false, "Measure aperture photometry corrections"); 99 99 return false; -
trunk/psphot/src/psphotReplaceUnfit.c
r9270 r9576 13 13 if (source->modelPSF == NULL) continue; 14 14 15 psTrace ( __func__, 3, "replacing object at %f,%f\n",15 psTrace ("psphot", 3, "replacing object at %f,%f\n", 16 16 source->modelPSF->params->data.F32[2], source->modelPSF->params->data.F32[3]); 17 17 -
trunk/psphot/src/psphotSourceFits.c
r9270 r9576 12 12 return status; 13 13 } 14 psTrace ("psphot .blend", 5, "trying blend...\n");14 psTrace ("psphot", 5, "trying blend...\n"); 15 15 16 16 // save the PSF model from the Ensemble fit … … 79 79 // if this one failed, skip it 80 80 if (!psphotEvalPSF (blend, model)) { 81 psTrace ("psphot .blend", 5, "failed on blend %d of %ld\n", i, modelSet->n);81 psTrace ("psphot", 5, "failed on blend %d of %ld\n", i, modelSet->n); 82 82 continue; 83 83 } … … 86 86 psFree(blend->modelPSF); 87 87 blend->modelPSF = psMemIncrRefCounter (model); 88 psTrace ("psphot .blend", 5, "fitted blend as PSF\n");88 psTrace ("psphot", 5, "fitted blend as PSF\n"); 89 89 pmModelSub (source->pixels, source->mask, model, false, false); 90 90 blend->mode |= PM_SOURCE_MODE_SUBTRACTED; … … 96 96 // evaluate the primary object 97 97 if (!psphotEvalPSF (source, PSF)) { 98 psTrace ("psphot .blend", 5, "failed on blend 0 of %ld\n", modelSet->n);98 psTrace ("psphot", 5, "failed on blend 0 of %ld\n", modelSet->n); 99 99 psFree (PSF); 100 100 return false; 101 101 } 102 102 103 psTrace ("psphot .blend", 5, "fitted primary as PSF\n");103 psTrace ("psphot", 5, "fitted primary as PSF\n"); 104 104 pmModelSub (source->pixels, source->mask, PSF, false, false); 105 105 psFree (source->modelPSF); … … 140 140 } 141 141 142 psTrace ("psphot .blend", 5, "fitted as PSF\n");142 psTrace ("psphot", 5, "fitted as PSF\n"); 143 143 pmModelSub (source->pixels, source->mask, PSF, false, false); 144 144 … … 189 189 if (!pmSourceMoments (source, EXT_MOMENTS_RAD)) return false; 190 190 191 psTrace ("psphot .blend", 5, "trying blob...\n");191 psTrace ("psphot", 5, "trying blob...\n"); 192 192 193 193 // this temporary source is used as a place-holder by the psphotEval functions below … … 219 219 220 220 if (okEXT && okDBL) { 221 psTrace ("psphot .blend", 5, "blob chisq: %f vs %f\n", chiEXT, chiDBL);221 psTrace ("psphot", 5, "blob chisq: %f vs %f\n", chiEXT, chiDBL); 222 222 // XXX EAM : a bogus bias: need to examine this better 223 223 if (3*chiEXT > chiDBL) goto keepDBL; … … 237 237 psFree (DBL); 238 238 pmModelSub (source->pixels, source->mask, EXT, false, false); 239 psTrace ("psphot .blend", 5, "blob as EXT: %f %f\n", EXT->params->data.F32[2], EXT->params->data.F32[3]);239 psTrace ("psphot", 5, "blob as EXT: %f %f\n", EXT->params->data.F32[2], EXT->params->data.F32[3]); 240 240 241 241 // save new model … … 250 250 pmModelSub (source->pixels, source->mask, (pmModel *) DBL->data[0], false, false); 251 251 pmModelSub (source->pixels, source->mask, (pmModel *) DBL->data[1], false, false); 252 psTrace ("psphot .blend", 5, "blob as DBL: %f %f\n", ONE->params->data.F32[2], ONE->params->data.F32[3]);252 psTrace ("psphot", 5, "blob as DBL: %f %f\n", ONE->params->data.F32[2], ONE->params->data.F32[3]); 253 253 254 254 // drop old model, save new second model... … … 339 339 340 340 if ((source->moments->Sx < 1e-3) || (source->moments->Sx < 1e-3)) { 341 psTrace ("psphot SourceFits", 5, "problem source: moments: %f %f\n", source->moments->Sx, source->moments->Sy);341 psTrace ("psphot", 5, "problem source: moments: %f %f\n", source->moments->Sx, source->moments->Sy); 342 342 } 343 343
Note:
See TracChangeset
for help on using the changeset viewer.
