Using TMath functions instead of math.h ones
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 14 Mar 2004 17:24:25 +0000 (17:24 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 14 Mar 2004 17:24:25 +0000 (17:24 +0000)
ITS/AliITSNeuralTrack.cxx
ITS/AliITSNeuralTrack.h
ITS/AliITSNeuralTracker.cxx
ITS/AliITSVertexerIons.cxx

index f19f4be..233d454 100644 (file)
@@ -347,7 +347,7 @@ Int_t AliITSNeuralTrack::PropagateTo(Double_t rk)
        // and the multiple scattering effects, which respectively have the effect
        // of changing the curvature and widening the covariance matrix.
 
-       if (rk < fabs(fDt)) {
+       if (rk < TMath::Abs(fDt)) {
                Error("PropagateTo", Form("Impossible propagation to r (=%17.15g) < Dt (=%17.15g)", rk, fDt));
                return 0;
        }
@@ -1101,7 +1101,7 @@ Double_t AliITSNeuralTrack::ArgZ(Double_t r) const
        Double_t arg;
        arg = (r * r - fDt * fDt) / (1. + fC * fDt);
        if (arg < 0.) {
-               if (fabs(arg) < 1.E-6) arg = 0.;
+               if (TMath::Abs(arg) < 1.E-6) arg = 0.;
                else {
                        Error("ArgZ", "Square root argument error: %17.15g < 0", arg);
                        return 10.;
index a00d014..3521e49 100644 (file)
@@ -60,10 +60,10 @@ public:
        Double_t GetZ(Double_t r) const;    // z = dz + (tanl / C) * asin(argz(rho))
 
        Double_t GetP()       {return GetPt() * (1.0 + fTanL * fTanL);}
-       Double_t GetPt()      {return 0.299792658 * 0.2 * fField * fabs(1./fC/100.);}
+       Double_t GetPt()      {return 0.299792658 * 0.2 * fField * TMath::Abs(1./fC/100.);}
        Double_t GetPz()      {return GetPt() * fTanL;}
-       Double_t GetE()       {return sqrt(fMass*fMass + GetPt()*GetPt());}
-       Double_t GetLambda()  {return atan(fTanL);}
+       Double_t GetE()       {return TMath::Sqrt(fMass*fMass + GetPt()*GetPt());}
+       Double_t GetLambda()  {return TMath::ATan(fTanL);}
        Int_t    GetPDGcode() const {return fPDG;}
        Double_t GetdEdX();
 
index 95a0725..28a4754 100644 (file)
@@ -44,7 +44,7 @@
 
 class AliITSNeuralPoint;
 
-using namespace std;
+//using namespace std;
 ClassImp(AliITSNeuralTracker)
 
 //--------------------------------------------------------------------------------------------
@@ -399,7 +399,7 @@ Bool_t AliITSNeuralTracker::PassCurvCut
        */
        // NEW VERSION
        if (den != 0.) {
-               curv = fabs(num / den);
+               curv = TMath::Abs(num / den);
                if (curv > fCurvCut[curvindex]) return kFALSE;
        }
        else
@@ -415,7 +415,7 @@ Bool_t AliITSNeuralTracker::PassCurvCut
        arc1 /= 2.0 * curv;
        arc2 /= 2.0 * curv;
        if (arc1 == 0.0 || arc2 == 0.0) return kFALSE;
-       helmatch = fabs(z1 / arc1 - z2 / arc2);
+       helmatch = TMath::Abs(z1 / arc1 - z2 / arc2);
        return (helmatch >= fHelixMatchCutMin[reflayer] && helmatch <= fHelixMatchCutMax[reflayer]);
        // END NEW VERSION
 }
@@ -481,7 +481,7 @@ Int_t AliITSNeuralTracker::PassAllCuts
        Double_t den = r1*r2*sqrt(dx*dx + dy*dy);
        Double_t curv = 0.;
        if (den != 0.) {
-               curv = fabs(num / den);
+               curv = TMath::Abs(num / den);
                if (curv > fCurvCut[curvindex]) return 3;
        }
        else
@@ -498,7 +498,7 @@ Int_t AliITSNeuralTracker::PassAllCuts
        arc1 /= 2.0 * curv;
        arc2 /= 2.0 * curv;
        if (arc1 == 0.0 || arc2 == 0.0) return kFALSE;
-       helmatch = fabs(z1 / arc1 - z2 / arc2);
+       helmatch = TMath::Abs(z1 / arc1 - z2 / arc2);
        if (helmatch < fHelixMatchCutMin[reflayer] || helmatch > fHelixMatchCutMax[reflayer]) return 5;
        
        return 0;
index 012c334..0eeac2f 100644 (file)
@@ -327,7 +327,7 @@ Double_t AliITSVertexerIons::FindMaxAround(Double_t point, TH1F *h, Double_t dis
   for(Int_t i=0;i<h->GetNbinsX();i++) {
     Int_t content=(Int_t)h->GetBinContent(i);
     Double_t center=(Double_t)h->GetBinCenter(i);
-    if(fabs(center-point)>distance) continue;
+    if(TMath::Abs(center-point)>distance) continue;
     if(content>max_content) {max_content=content;max_bin=i;}
   }
   Double_t max=h->GetBinCenter(max_bin);