Changeset 36327 for branches/eam_branches/ipp-20130904/psphot/src/psphotPetrosianStats.c

Timestamp:

Nov 30, 2013, 1:35:54 PM (13 years ago)

Author:

eugene

Message:

add a number of additional flag bits, esp for model fitting

File:

• branches/eam_branches/ipp-20130904/psphot/src/psphotPetrosianStats.c (modified) (4 diffs)

branches/eam_branches/ipp-20130904/psphot/src/psphotPetrosianStats.c

@@ -21 +21 @@
     if (!profile->binSB) {
         psLogMsg ("psphot", PS_LOG_DETAIL, "no petrosian profile, skipping source %f, %f", source->peak->xf, source->peak->yf);
+        source->mode2 |= PM_SOURCE_MODE2_PETRO_NO_PROFILE;
         return true;
     }
@@ -118 +119 @@
                 petRadius    = InterpolateValues     (1.0, 0.0, petRatio->data.F32[nOut], refRadius->data.F32[nOut], PETROSIAN_RATIO);
                 petRadiusErr = InterpolateValuesErrX (1.0, 0.0, petRatio->data.F32[nOut], refRadius->data.F32[nOut], PETROSIAN_RATIO, 0.0, petRatioErr->data.F32[nOut]);
+                source->mode2 |= PM_SOURCE_MODE2_PETRO_RATIO_ZEROBIN;
             } else {
               // petRadius    = InterpolateValues     (petRatio->data.F32[nOut-1], refRadius->data.F32[nOut-1], petRatio->data.F32[nOut], refRadius->data.F32[nOut], PETROSIAN_RATIO);
               if (nOut > 1) {
-                petRadius    = InterpolateValuesQuadratic (&petRatio->data.F32[nOut-2],   &refRadius->data.F32[nOut-2],   PETROSIAN_RATIO);
+                petRadius    = InterpolateValuesQuadratic (&petRatio->data.F32[nOut-2], &refRadius->data.F32[nOut-2],   PETROSIAN_RATIO);
               } else {
                 petRadius    = InterpolateValuesQuadratic (&petRatio->data.F32[nOut-3], &refRadius->data.F32[nOut-3], PETROSIAN_RATIO);
@@ -158 +160 @@
             fprintf (stderr, "nan pet radius\n");
         }
+        source->mode2 |= PM_SOURCE_MODE2_PETRO_INSIG_RATIO;
     }
 
@@ -190 +193 @@
     bool found50 = false;
     bool found90 = false;
+
     // XXX use bisection to do this faster:
     for (int i = 0; !(found50 && found90) && i < refRadius->n; i++) {