Changeset 34091 for trunk/ppStack/src/ppStackPSF.c

Timestamp:

Jun 26, 2012, 1:09:06 PM (14 years ago)

Author:

mhuber

Message:

adding option to set target FWHM mean of simple PSF for convolved stacks

File:

• trunk/ppStack/src/ppStackPSF.c (modified) (2 diffs)

trunk/ppStack/src/ppStackPSF.c

@@ -4 +4 @@
 
 pmPSF *ppStackPSF(const pmConfig *config, int numCols, int numRows,
-                  const psArray *psfs, const psVector *inputMask)
+                  const psArray *psfs, ppStackOptions *options)
 {
     bool mdok = false;
     pmPSF *psf = NULL;
+
+    const psVector *inputMask = options->inputMask;
 
     // Get the recipe values
@@ -55 +57 @@
         }
 
+        //MEH -- preliminary add option -1 for choosing mean FWHM of inputs. Nsig<0 for scaling +Nsig by mean
+        if (psfFWHM<0.0){
+            float NpsfFWHMsig = psMetadataLookupF32(&mdok, recipe, "PSF.OUTPUT.NFWHMSIG"); // Nsig to include
+            if (!mdok){
+                NpsfFWHMsig = 0.0;
+            }
+            float clippedMean = options->clippedMean;
+            float clippedStdev = options->clippedStdev;
+            psfFWHM = clippedMean + NpsfFWHMsig*clippedStdev;
+            if (NpsfFWHMsig<0.0){
+                psfFWHM = clippedMean + NpsfFWHMsig * -1.0*clippedStdev/clippedMean;
+            }
+        }
+
+
+
         float Mxx = M_SQRT2 * psfFWHM / 2.35;