Changeset 25612

Timestamp:

Sep 27, 2009, 10:53:05 AM (17 years ago)

Author:

eugene

Message:

set the PSF fit radius and aperture based on the window function (identical to psf model generation)

File:

• branches/eam_branches/20090715/psphot/src/psphotChoosePSF.c (modified) (2 diffs)

branches/eam_branches/20090715/psphot/src/psphotChoosePSF.c

@@ -76 +76 @@
     // assert (status);
 
-    // XXX we have calculated a Gaussian window function, use that for both the PSF fit radius
-    // and the aperture radius (scaling SIGMA)
-    // XXX put the scale factors into the recipe
+    // We have calculated a Gaussian window function, use that for both the PSF fit radius and
+    // the aperture radius (scaling SIGMA)
     float gaussSigma = psMetadataLookupF32(&status, recipe, "MOMENTS_GAUSS_SIGMA");
     float fitScale = psMetadataLookupF32(&status, recipe, "PSF_FIT_RADIUS_SCALE");
@@ -84 +83 @@
     options->fitRadius = (int)(fitScale*gaussSigma);
     options->apRadius = (int)(apScale*gaussSigma);
+
+    // XXX use the same radii for standard analysis as for the PSF creation
+    psMetadataAddF32(recipe, PS_LIST_TAIL, "PSF_FIT_RADIUS", PS_META_REPLACE, "fit radius", options->fitRadius);
+    psMetadataAddF32(recipe, PS_LIST_TAIL, "PSF_APERTURE", PS_META_REPLACE, "psf aperture", options->apRadius);
 
     // XXX ROBUST seems to be too agressive given the small numbers.