Changeset 19203 for trunk/Mars/mimage

Timestamp:

09/01/18 14:10:18 (6 years ago)

Author:

tbretz

Message:

Still requires numerical protection for a few very unlikely numerical uncertainty cases.

File:

• trunk/Mars/mimage/MHillasExt.cc (modified) (1 diff)

trunk/Mars/mimage/MHillasExt.cc

@@ -299 +299 @@
     const Double_t meandtw  = sumdtw/sumw;
 
-    fTimeSpread          = TMath::Sqrt(sumt2 /cnt  - meant *meant);
-    fTimeSpreadWeighted  = TMath::Sqrt(sumt2w/sumw - meantw*meantw);
-
-    fSlopeSpread         = TMath::Sqrt(sumdt2 /cnt  - meandt *meandt);
-    fSlopeSpreadWeighted = TMath::Sqrt(sumdt2w/sumw - meandtw*meandtw);
+    const Double_t rad_ts   = sumt2  /cnt  - meant *meant;
+    const Double_t rad_tsw  = sumt2w /sumw - meantw*meantw;
+    const Double_t rad_ss   = sumdt2 /cnt  - meandt *meandt;
+    const Double_t rad_ssw  = sumdt2w/sumw - meandtw*meandtw;
+
+    fTimeSpread          = rad_ts <0 ? 0 : TMath::Sqrt(rad_ts);
+    fTimeSpreadWeighted  = rad_tsw<0 ? 0 : TMath::Sqrt(rad_tsw);
+
+    fSlopeSpread         = rad_ss <0 ? 0 : TMath::Sqrt(rad_ss);
+    fSlopeSpreadWeighted = rad_ssw<0 ? 0 : TMath::Sqrt(rad_ssw);
 
     //