Changeset 33838

Timestamp:

May 1, 2012, 11:43:50 AM (14 years ago)

Author:

bills

Message:

Limit moments->Mrf minimum value minKronRadius.
For sources with Moments->SN < 10 limit the maximum value to 'radius'
which is passed in as an argument taken from the value of PSF_MOMENTS_RADIUS in
readout->analysis

File:

• trunk/psModules/src/objects/pmSourceMoments.c (modified) (1 diff)

trunk/psModules/src/objects/pmSourceMoments.c

@@ -319 +319 @@
     }
 
-    source->moments->Mrf = RF/RS;
     source->moments->Mrh = RH/RS;
 
-    // if Mrf (first radial moment) is very small, we are getting into low-significance
+    // if Mrf = RF/RS (first radial moment) is very small, we are getting into low-significance
     // territory.  saturate at minKronRadius.  conversely, if Mrf is >> radius for faint
     // sources, we are clearly making an error.  saturate at radius.
-    float kronRefRadius = MAX(minKronRadius, source->moments->Mrf);
+    float kronRefRadius = MAX(minKronRadius, RF/RS);
     if (source->moments->SN < 10) {
         kronRefRadius = MIN(radius, kronRefRadius);
     }
+    source->moments->Mrf = kronRefRadius;
 
     // *** now calculate the kron flux values using the 1st radial moment