Changeset 19532 for trunk/ppStack/src/ppStack.h

Timestamp:

Sep 11, 2008, 6:12:42 PM (18 years ago)

Author:

Paul Price

Message:

Correcting stack variances for the imperfect convolution. Changed (again) the way the deviations are calculated; they're now straight chi^{2, which allows us to estimate how the variance factor needs to be changed in order to correct for it. Moving additional values into pmSubtractionKernels so that this can be calculated from a kernel that has been read in (very useful for debugging, also nice to have a record of the quality of the fit).}

File:

• trunk/ppStack/src/ppStack.h (modified) (2 diffs)

trunk/ppStack/src/ppStack.h

@@ -89 +89 @@
                                const psArray *readouts, // Input readouts
                                const psArray *regions, // Array with array of regions used in each PSF match
-                               const psArray *kernels // Array with array of kernels used in each PSF match
+                               const psArray *kernels, // Array with array of kernels used in each PSF match
+                               const psVector *addVariance // Additional variance for rejection
     );
 
@@ -128 +129 @@
 
 /// Convolve image to match specified seeing
-bool ppStackMatch(pmReadout *readout, // Readout to be convolved; replaced with output
-                  psArray **regions, // Array of regions used in each PSF matching, returned
-                  psArray **kernels, // Array of kernels used in each PSF matching, returned
+bool ppStackMatch(pmReadout *readout,   // Readout to be convolved; replaced with output
+                  psArray **regions,    // Array of regions used in each PSF matching, returned
+                  psArray **kernels,    // Array of kernels used in each PSF matching, returned
+                  float *chi2,          // Chi^2 from the stamps
                   const psArray *sources, // Array of sources
                   const pmPSF *psf,     // Target PSF