Changeset 30865 for branches/eam_branches/ipp-20110213/psModules/src/objects/pmSourceMoments.c

Timestamp:

Mar 10, 2011, 5:23:09 PM (15 years ago)

Author:

eugene

Message:

measure a core flux in the moment analysis - this can be used with the kron flux to flag (and exclude?) detections in the badly-subtracted background regions

File:

• branches/eam_branches/ipp-20110213/psModules/src/objects/pmSourceMoments.c (modified) (5 diffs)

branches/eam_branches/ipp-20110213/psModules/src/objects/pmSourceMoments.c

@@ -82 +82 @@
     psF32 Sum = 0.0;
     psF32 Var = 0.0;
+    psF32 SumCore = 0.0;
+    psF32 VarCore = 0.0;
     psF32 R2 = PS_SQR(radius);
     psF32 minSN2 = PS_SQR(minSN);
@@ -248 +250 @@
 
     int nKronPix = 0;
+    int nCorePix = 0;
+    int nInner = 0;
+    int nOuter = 0;
     Sum = Var = 0.0;
     float SumInner = 0.0;
@@ -255 +260 @@
 
         psF32 yDiff = row - yCM;
-        if (fabs(yDiff) > radKron) continue;
+        if (fabs(yDiff) > radKouter) continue;
 
         psF32 *vPix = source->pixels->data.F32[row];
@@ -272 +277 @@
 
             psF32 xDiff = col - xCM;
-            if (fabs(xDiff) > radKron) continue;
+            if (fabs(xDiff) > radKouter) continue;
 
             // radKron is just a function of (xDiff, yDiff)
@@ -293 +298 @@
             }
 
+            // use sigma (fixed by psf) not a radKron based value
+            if (r < sigma) {
+                SumCore += pDiff;
+                VarCore += wDiff;
+                nCorePix ++;
+            }
+
             if ((r > radKinner) && (r < radKron)) {
                 SumInner += pDiff;
+                nInner ++;
             }
             if ((r > radKron)  && (r < radKouter)) {
                 SumOuter += pDiff;
+                nOuter ++;
             }
         }
     }
-    source->moments->KronFlux = Sum;
-    source->moments->KronFinner = SumInner;
-    source->moments->KronFouter = SumOuter;
-    source->moments->KronFluxErr = sqrt(Var);
+    // *** should I rescale these fluxes by pi R^2 / nNpix?
+    source->moments->KronCore    = SumCore       * M_PI * PS_SQR(sigma) / nCorePix;
+    source->moments->KronCoreErr = sqrt(VarCore) * M_PI * PS_SQR(sigma) / nCorePix;
+    source->moments->KronFlux    = Sum       * M_PI * PS_SQR(radKron) / nKronPix;
+    source->moments->KronFluxErr = sqrt(Var) * M_PI * PS_SQR(radKron) / nKronPix;
+    source->moments->KronFinner = SumInner * M_PI * (PS_SQR(radKron)   - PS_SQR(radKinner)) / nInner;
+    source->moments->KronFouter = SumOuter * M_PI * (PS_SQR(radKouter) -   PS_SQR(radKron)) / nOuter;
 
     psTrace ("psModules.objects", 4, "Mrf: %f  KronFlux: %f  Mxx: %f  Mxy: %f  Myy: %f  Mxxx: %f  Mxxy: %f  Mxyy: %f  Myyy: %f  Mxxxx: %f  Mxxxy: %f  Mxxyy: %f  Mxyyy: %f  Mxyyy: %f\n",