Changeset 36115 for trunk/psphot/src

Timestamp:

Sep 13, 2013, 10:37:42 AM (13 years ago)

Author:

eugene

Message:

petrosian ratio error was mis-calculated : partial fix using only skynoise for the outer annulus (needs more careful treatment)

File:

• trunk/psphot/src/psphotPetrosianRadialBins.c (modified) (2 diffs)

trunk/psphot/src/psphotPetrosianRadialBins.c

@@ -36 +36 @@
     pmSourceRadialProfile *profile = source->extpars->petProfile;
 
-    float skyModelErrorSQ = PS_SQR(skynoise);
+    // float skyModelErrorSQ = PS_SQR(skynoise);
 
     int nMax = radiusMax;
@@ -133 +133 @@
         if (radius->data.F32[i] > Rmax) {
             // calculate the value for the nOut bin
-            float value, dvalue;
+            float value; // , dvalue;
             if (values->n > 0) {
                 psVectorStats (stats, values, NULL, NULL, 0);
                 value = stats->robustMedian;
-                dvalue = stats->robustStdev;
+                // dvalue = stats->robustStdev;
             } else {
                 value = NAN;
-                dvalue = NAN;
+                // dvalue = NAN;
             }
 
             binSB->data.F32[nOut] = value;
-            binSBstdev->data.F32[nOut] = sqrt(PS_SQR(dvalue) / values->n + skyModelErrorSQ);
+            // binSBstdev->data.F32[nOut] = sqrt(PS_SQR(dvalue) / values->n + skyModelErrorSQ);
+            binSBstdev->data.F32[nOut] = skynoise / sqrt(values->n);
             binFill->data.F32[nOut] = values->n / binArea->data.F32[nOut];