Changeset 30101 for branches/eam_branches/ipp-20101205/psphot/src/psphotStackPSF.c

Timestamp:

Dec 17, 2010, 9:57:34 AM (16 years ago)

Author:

eugene

Message:

various modifications to psphotStack: add needed metadata to header; enable multiple PSF-matched images; work to speed up radial aperture photometry; radialAper is now an array (one element per PSF size); better feedback on extended fit skips; clarify processing: main analysis is on the input images (raw or convolved) while radial apertures is on the psf-matched images

File:

• branches/eam_branches/ipp-20101205/psphot/src/psphotStackPSF.c (modified) (2 diffs)

branches/eam_branches/ipp-20101205/psphot/src/psphotStackPSF.c

@@ -12 +12 @@
     bool autoPSF = psMetadataLookupBool (&mdok, psphotRecipe, "PSPHOT.STACK.TARGET.PSF.AUTO");
 
+    // Get the recipe values
+    psMetadata *recipe = psMetadataLookupMetadata(NULL, config->recipes, "PPSTACK"); // ppStack recipe
+    psAssert(recipe, "We've thrown an error on this before.");
+        
+    char *psfModel = psMetadataLookupStr(NULL, recipe, "PSF.MODEL"); // Model for PSF
+
     if (autoPSF) {
-        // Get the recipe values
-        psMetadata *recipe = psMetadataLookupMetadata(NULL, config->recipes, "PPSTACK"); // ppStack recipe
-        psAssert(recipe, "We've thrown an error on this before.");
-
         int psfInstances = psMetadataLookupS32(NULL, recipe, "PSF.INSTANCES"); // Number of instances for PSF
         float psfRadius = psMetadataLookupF32(NULL, recipe, "PSF.RADIUS"); // Radius for PSF
-        const char *psfModel = psMetadataLookupStr(NULL, recipe, "PSF.MODEL"); // Model for PSF
         int psfOrder = psMetadataLookupS32(NULL, recipe, "PSF.ORDER"); // Spatial order for PSF
 
@@ -49 +50 @@
 
         float targetFWHM = psMetadataLookupF32 (&mdok, psphotRecipe, "PSPHOT.STACK.TARGET.PSF.FWHM");
-        psAssert (isfinite(targetFWHM), "missing psphot recipe value PSPHOT.STACK.TARGET.PSF.FWHM");
+        if (!mdok) {
+            psVector *fwhmValues = psMetadataLookupVector(&mdok, psphotRecipe, "PSPHOT.STACK.TARGET.PSF.FWHM"); // Magnitude offsets
+            psAssert (mdok, "missing psphot recipe value PSPHOT.STACK.TARGET.PSF.FWHM");
+            targetFWHM = fwhmValues->data.F32[0];
+        }
 
         float Sxx = sqrt(2.0)*targetFWHM / 2.35;
 
         // XXX probably should make the model type (and par 7) optional from recipe
-        psf = pmPSFBuildSimple("PS_MODEL_PS1_V1", Sxx, Sxx, 0.0, 1.0);
+        psf = pmPSFBuildSimple(psfModel, Sxx, Sxx, 0.0, 1.0);
         if (!psf) {
             psError(PSPHOT_ERR_PSF, false, "Unable to build dummy PSF.");