Protections agains Arithmetic Exc. added

diff --git a/HBTAN/AliHBTTrackPoints.cxx b/HBTAN/AliHBTTrackPoints.cxx
index f7ec1a042dc0afda38ad323ea8422e58a6ddd24b..82552b38a09864ca4c6f0612e5c8282506494928 100644 (file)
--- a/HBTAN/AliHBTTrackPoints.cxx
+++ b/HBTAN/AliHBTTrackPoints.cxx
@@ -166,10 +166,40 @@ void AliHBTTrackPoints::MakePoints( Float_t dr, Float_t r0, Double_t x, Double_t
   for(Int_t i = 0; i<fN; i++)
    {
      Double_t rc = r0 + i*dr;
-     Double_t factorPhi = TMath::ASin( (c2*rc + cst1/rc)/cst2 );//factor phi od rc
+     Double_t ftmp = (c2*rc + cst1/rc)/cst2;
+     if (ftmp > 1.0) 
+      {
+        Warning("AliHBTTrackPoints","ASin argument > 1 %f:",ftmp);
+        ftmp=1.0;
+      }
+     else if (ftmp < -1.0) 
+      {
+        Warning("AliHBTTrackPoints","ASin argument < -1 %f:",ftmp);
+        ftmp=-1.0;
+      }
+      
+     Double_t factorPhi = TMath::ASin( ftmp );//factor phi od rc
      Double_t phi = phi0global + factorPhi - factorPhi0;
      
-     Double_t factorZ = TMath::ASin(c2*TMath::Sqrt((rc*rc-dcasq)/cst2))*par[3]/c2;
+     ftmp = (rc*rc-dcasq)/cst2;
+     if (ftmp < 0.0) 
+      {
+        Warning("AliHBTTrackPoints","Sqrt argument < 0: %f",ftmp);
+        ftmp=1.0;
+      }
+     
+     ftmp = c2*TMath::Sqrt(ftmp);
+     if (ftmp > 1.0) 
+      {
+        Warning("AliHBTTrackPoints","ASin argument > 1: %f",ftmp);
+        ftmp=1.0;
+      }
+     else if (ftmp < -1.0) 
+      {
+        Warning("AliHBTTrackPoints","ASin argument < -1: %f",ftmp);
+        ftmp=-1.0;
+      }
+     Double_t factorZ = TMath::ASin(ftmp)*par[3]/c2;
      fZ[i] = z0 + factorZ - factorZ0;
      fX[i] = rc*TMath::Cos(phi);
      fY[i] = rc*TMath::Sin(phi);