// 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;
}
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.;
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();
class AliITSNeuralPoint;
-using namespace std;
+//using namespace std;
ClassImp(AliITSNeuralTracker)
//--------------------------------------------------------------------------------------------
*/
// NEW VERSION
if (den != 0.) {
- curv = fabs(num / den);
+ curv = TMath::Abs(num / den);
if (curv > fCurvCut[curvindex]) return kFALSE;
}
else
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
}
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
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;
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);