Changeset 21338

Timestamp:

Feb 5, 2009, 2:04:35 PM (17 years ago)

Author:

Paul Price

Message:

Correct the correction factor for the covariance produced by the (potentially multiple) smoothing

File:

• branches/pap_branch_20090128/psphot/src/psphotSignificanceImage.c (modified) (1 diff)

branches/pap_branch_20090128/psphot/src/psphotSignificanceImage.c

@@ -100 +100 @@
 
     float factor = guess ? 4.0 * M_PI * PS_SQR(SIGMA_SMTH) : 4.0 * M_PI * PS_SQR(SIGMA_SMTH);
+
+    // Correct the correction factor for the covariance produced by the (potentially multiple) smoothing
+    psKernel *kernel = psImageSmoothKernel(SIGMA_SMTH, NSIGMA_SMTH); // Kernel used for smoothing
+    psKernel *covar = psImageCovarianceCalculate(kernel, readout->covariance); // Covariance matrix
+    psFree(kernel);
+    factor /= psImageCovarianceFactor(covar);
+    psFree(covar);
+
     // record the effective area and significance scaling factor
     float effArea = 8.0 * M_PI * PS_SQR(SIGMA_SMTH);