Changeset 31264

Timestamp:

Apr 11, 2011, 12:07:50 PM (15 years ago)

Author:

eugene

Message:

fix case were threshold to max range is surprisingly small

File:

• tags/ipp-20110406/psModules/src/objects/pmFootprintCullPeaks.c (modified) (1 diff)

tags/ipp-20110406/psModules/src/objects/pmFootprintCullPeaks.c

@@ -102 +102 @@
         float beta = nsigma_delta * ALPHA;
         float beta2 = PS_SQR(beta);
-        int nBins = sqrt(4.0*(maxFlux - min_threshold)/beta2) + 10; // let's be extra generous
+        float value = 4.0*(maxFlux - min_threshold)/beta2;
+        int nBins = (value > 1.0) ? sqrt(value) + 10 : 10; // let's be extra generous
 
         // create a vector to store the threshold bins used for each peak