Changeset 6825 for branches/rel10_ifa/psModules/src/objects/pmSourceFitSet.c

Timestamp:

Apr 10, 2006, 10:21:36 AM (20 years ago)

Author:

magnier

Message:

adding chisq-flux trend corrections; adding option to use non-poisson errors for fits

File:

• branches/rel10_ifa/psModules/src/objects/pmSourceFitSet.c (modified) (3 diffs)

branches/rel10_ifa/psModules/src/objects/pmSourceFitSet.c

@@ -3 +3 @@
  *  @author EAM, IfA
  *
- *  @version $Revision: 1.1.2.3 $ $Name: not supported by cvs2svn $
- *  @date $Date: 2006-04-01 02:47:20 $
+ *  @version $Revision: 1.1.2.4 $ $Name: not supported by cvs2svn $
+ *  @date $Date: 2006-04-10 20:21:36 $
  *
  *  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
@@ -182 +182 @@
     psS32 nPix = source->pixels->numRows * source->pixels->numCols;
 
+    // local sky variance
+    psF32 dSky = source->moments->dSky;
+
     // construct the coordinate and value entries
     psArray *x = psArrayAlloc(nPix);
@@ -206 +209 @@
             y->data.F32[nPix] = source->pixels->data.F32[i][j];
             // psMinimizeLMChi2 takes wt = 1/dY^2
-            yErr->data.F32[nPix] = 1.0 / source->weight->data.F32[i][j];
+            if (PM_PSF_POISSON_WEIGHTS) {
+                yErr->data.F32[nPix] = 1.0 / source->weight->data.F32[i][j];
+            } else {
+                yErr->data.F32[nPix] = 1.0 / dSky;
+            }
             nPix++;
         }