added protection against floating point exception (Alexander)
authorjthaeder <jthaeder@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 11 Sep 2011 23:01:40 +0000 (23:01 +0000)
committerjthaeder <jthaeder@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 11 Sep 2011 23:01:40 +0000 (23:01 +0000)
TPC/AliTPCseed.cxx

index 1e8857d..be244e8 100644 (file)
@@ -1332,9 +1332,18 @@ Float_t  AliTPCseed::CookdEdxAnalytical(Double_t low, Double_t up, Int_t type, I
       sumDN++;
     }
   }
+  //
+  Float_t mean = 0;
+  Float_t meanL = 0;  
+  Float_t meanD = 0;           // lower half mean
+  if (sumn > 1e-30)   mean =suma/sumn;
+  if (sumLN > 1e-30)  meanL =sumL/sumLN;
+  if (sumDN > 1e-30)  meanD =(sumD/sumDN);
+  /*
   Float_t mean =suma/sumn;
   Float_t meanL = sumL/sumLN;  
   Float_t meanD =(sumD/sumDN);           // lower half mean
+  */
   Float_t rms  =TMath::Sqrt(TMath::Abs(suma2/sumn-mean*mean));
   Float_t mean2=suma2/sumn;
   Float_t mean3=suma3/sumn;