Changeset 20712

Timestamp:

Nov 12, 2008, 5:50:30 PM (17 years ago)

Author:

eugene

Message:

var_soft = var + soft_frac
nsigma_^{2 = var_soft * rej
(not the other way around)}

File:

• trunk/psModules/src/imcombine/pmStack.c (modified) (2 diffs)

trunk/psModules/src/imcombine/pmStack.c

@@ -8 +8 @@
  *  @author GLG, MHPCC
  *
- *  @version $Revision: 1.44 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2008-11-01 02:59:33 $
+ *  @version $Revision: 1.45 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2008-11-13 03:50:30 $
  *  Copyright 2004-2007 Institute for Astronomy, University of Hawaii
  *
@@ -353 +353 @@
                   float rej2 = PS_SQR(rej); // Rejection level squared
                   for (int i = 0; i < num; i++) {
-                      pixelVariances->data.F32[i] *= rej2;
+                      // Systematic error contributes to the rejection level
+                      pixelVariances->data.F32[i] += PS_SQR(sys * pixelData->data.F32[i]);
 #ifdef VARIANCE_FACTORS
                       // Variance factor contributes to the rejection level
                       pixelVariances->data.F32[i] *= varFactors->data.F32[pixelSources->data.U16[i]];
 #endif
-                      // Systematic error contributes to the rejection level
-                      pixelVariances->data.F32[i] += PS_SQR(sys * pixelData->data.F32[i]);
+                      pixelVariances->data.F32[i] *= rej2;
                   }
               }