Changeset 17557 for trunk/ppSim/src/ppSimSetPSF.c

Timestamp:

May 7, 2008, 10:55:21 AM (18 years ago)

Author:

eugene

Message:

fix errors with star density normalizations and flux consistency

File:

• trunk/ppSim/src/ppSimSetPSF.c (modified) (3 diffs)

trunk/ppSim/src/ppSimSetPSF.c

@@ -1 +1 @@
 # include "ppSim.h"
-static char *defaultModel = "PS_MODEL_GAUSS";
+static char *defaultModel = "PS_MODEL_QGAUSS";
 
 bool ppSimSetPSF (pmChip *chip, pmConfig *config) {
 
-    bool status;
+    bool status, mdok;
     pmPSF *psf = NULL;
     pmTrend2D *param = NULL;
+
+    psMetadata *recipe = psMetadataLookupMetadata(&mdok, config->recipes, PPSIM_RECIPE); // Recipe
 
     // the pmPSF IO functions stores the PSF on the chip->analysis
@@ -14 +16 @@
     }
 
-    // no supplied PSF, build one using supplied value for seeing
-    // XXX need to correct for the pixel scale
-    float seeing   = psMetadataLookupF32(&status, config->arguments, "SEEING"); // Seeing sigma (pix)
+    // no supplied PSF, build one using supplied value for seeing seeing is already corrected
+    // for the pixel scale, and is converted from FWHM to SIGMA (this is done in
+    // ppSimArguments)
+    float seeing   = psMetadataLookupF32(&status, recipe, "SEEING"); // Seeing SIGMA (pixels)
 
     char *psfModelName = psMetadataLookupStr(&status, config->arguments, "PSF.MODEL"); // Name of PSF model
@@ -50 +53 @@
     psEllipsePol pol;
 
-    // supply the semi-major axis
+    // supply the semi-major axis (these are SIGMA values in PIXELS)
     axes.major = seeing;
     axes.minor = seeing;