Changeset 29312 for branches/eam_branches/ipp-20100823/psModules/src/objects/pmSourcePhotometry.c

Timestamp:

Oct 4, 2010, 2:48:37 PM (16 years ago)

Author:

eugene

Message:

dipole measurements redefined to match definitions from Armin (do not filter pixels by S/N)

File:

• branches/eam_branches/ipp-20100823/psModules/src/objects/pmSourcePhotometry.c (modified) (2 diffs)

branches/eam_branches/ipp-20100823/psModules/src/objects/pmSourcePhotometry.c

@@ -444 +444 @@
     }
 
+    // NOTE: until 2010.10.01, these measurements included a 3sigma-per-pixel significance
+    // this followed what we understood as the definition given to us
+    // by Armin, but it always seemed a poor idea -- a faint source is unlikely to have any 3sigma pixels.
+    // changed to remove the per-pixel filter.
+
     for (int iy = 0; iy < flux->numRows; iy++) {
         for (int ix = 0; ix < flux->numCols; ix++) {
             // only count up the stats in the unmarked region (ie, the aperture)
+            // skip the marked pixels; these are not relevant
             if (mask->data.PS_TYPE_IMAGE_MASK_DATA[iy][ix] & markVal) {
                 continue;
@@ -455 +461 @@
             }
 
-            float SN = flux->data.F32[iy][ix] / sqrt(variance->data.F32[iy][ix]);
-
-            if (SN > +FLUX_LIMIT) {
+            float flux = flux->data.F32[iy][ix];
+
+            if (flux > 0.0) {
                 nGood ++;
-                fGood += fabs(flux->data.F32[iy][ix]);
-            }
-
-            if (SN < -FLUX_LIMIT) {
+                fGood += fabs(flux);
+            } else {
                 nBad ++;
-                fBad += fabs(flux->data.F32[iy][ix]);
+                fBad += fabs(flux-);
             }
         }