Changeset 30772 for branches/eam_branches/ipp-20110213/psphot/src/psphotStackPSF.c

Timestamp:

Mar 1, 2011, 6:21:28 PM (15 years ago)

Author:

eugene

Message:

fix threading for psphotStack; use standard PS_ARGUMENTS_ macros for ppSub, ppStack, pswarp, ppImage, psphotStack; fix FPA.EXPOSURE in ppImageAddNoise; fix accounting for psf bad and poor weights; added flags for substantial poor masks due to each poor mask bit type; fix the measurement of the fwhm in pmPSFtoFWHM

File:

• branches/eam_branches/ipp-20110213/psphot/src/psphotStackPSF.c (modified) (1 diff)

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

@@ -56 +56 @@
         }
 
-        float Sxx = sqrt(2.0)*targetFWHM / 2.35;
+        // measured scale factors (fwhm = Sxx * 2.35 * scaleFactor / sqrt(2.0))
+        // GAUSS  : 1.000
+        // PGAUSS : 1.006
+        // QGAUSS : 1.151
+        // RGAUSS : 0.883
+        // PS1_V1 : 1.134
+        
+        float scaleFactor = NAN;
+        if (!strcmp(psfModel, "PS_MODEL_GAUSS")) {
+            scaleFactor = 1.000;
+        }
+        if (!strcmp(psfModel, "PS_MODEL_PGAUSS")) {
+            scaleFactor = 1.0006;
+        }
+        if (!strcmp(psfModel, "PS_MODEL_QGAUSS")) {
+            scaleFactor = 1.151;
+        }
+        if (!strcmp(psfModel, "PS_MODEL_RGAUSS")) {
+            scaleFactor = 0.883;
+        }
+        if (!strcmp(psfModel, "PS_MODEL_PS1_V1")) {
+            scaleFactor = 1.134;
+        }
+        psAssert (isfinite(scaleFactor), "invalid model for PSF"); 
+
+        float Sxx = sqrt(2.0)*targetFWHM / 2.35 / scaleFactor;
 
         // XXX probably should make the model type (and par 7) optional from recipe
-        psf = pmPSFBuildSimple(psfModel, Sxx, Sxx, 0.0, 1.0);
+        // psf = pmPSFBuildSimple(psfModel, Sxx, Sxx, 0.0, 1.0);
+        psf = pmPSFBuildSimple(psfModel, Sxx, Sxx, 0.0, 0.2);
         if (!psf) {
             psError(PSPHOT_ERR_PSF, false, "Unable to build dummy PSF.");