GetPt() changed to GetSignedPt(). Get1Pt() changed to GetSigned1Pt() (Markus)
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 7 Aug 2007 21:15:43 +0000 (21:15 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 7 Aug 2007 21:15:43 +0000 (21:15 +0000)
24 files changed:
EMCAL/AliEMCALTracker.cxx
HLT/ITS/AliHLTITStrack.cxx
HLT/TRD/AliTRDtrackerHLT.cxx
ITS/AliITStrackMI.cxx
ITS/AliITStrackV2.cxx
ITS/AliITStrackerMI.cxx
ITS/oldmacros/AliITSImpactParametersPP.C
STEER/AliKFParticle.cxx
STEER/AliKFParticleTest.C
STEER/AliTPCComparisonMI.C
TPC/AliTPCComparisonMI.C
TPC/AliTPCseed.cxx
TPC/AliTPCtrack.cxx
TPC/AliTPCtracker.cxx
TPC/AliTPCtrackerParam.cxx
TPC/TPCcalib/AliTPCcalibTracks.cxx
TRD/AliTRDCalibraFillHisto.cxx
TRD/AliTRDtrack.cxx
TRD/AliTRDtracker.cxx
TRD/AliTRDtrackingAnalysis.cxx
TRD/Macros/AliTRDComparison.C
TRD/Macros/AliTRDComparisonV2.C
TRD/Macros/AliTRDbackTrackAnalysis.C
TRD/Macros/AliTRDseedsAnalysis.C

index 127917f..d0b4599 100644 (file)
@@ -895,7 +895,7 @@ Int_t AliEMCALTracker::CreateMatches()
                                candidate->SetIndexT(it);
                                candidate->SetIndexC(ic);
                                candidate->SetDistance(dist);
-                               candidate->SetPt(track->GetPt());
+                               candidate->SetPt(track->GetSignedPt());
                                fMatches->Add(candidate);
                                count++;
                        }
index d484257..ac0b195 100644 (file)
@@ -58,8 +58,8 @@ Int_t AliHLTITStrack::Compare(const TObject *o) const {
   // This function compares tracks according to the their curvature
   //-----------------------------------------------------------------
   AliHLTITStrack *t=(AliHLTITStrack*)o;
-  Double_t co=TMath::Abs(t->Get1Pt());
-  Double_t c =TMath::Abs(Get1Pt());
+  Double_t co=t->OneOverPt();
+  Double_t c =OneOverPt();
   //  Double_t co=t->GetSigmaY2()*t->GetSigmaZ2();
   //  Double_t c =GetSigmaY2()*GetSigmaZ2();
   if (c>co) return 1;
index cf4ddd4..33d2e14 100644 (file)
@@ -793,7 +793,7 @@ Int_t AliTRDtrackerHLT::PropagateBack(AliESDEvent *event)
     /**/
     // inter-tracks competition ???
     if ((TMath::Abs(track->GetC() - p4) / TMath::Abs(p4) < 0.2) || 
-        (TMath::Abs(track->GetPt())                      > 0.8)) {
+        (track->Pt()                                     > 0.8)) {
       
       fHBackfit->Fill(4);
       
@@ -813,7 +813,7 @@ Int_t AliTRDtrackerHLT::PropagateBack(AliESDEvent *event)
         // Sign only gold tracks
        if (track->GetChi2() / track->GetNumberOfClusters() < 4) {   
          if ((seed->GetKinkIndex(0)      ==   0) &&
-              (TMath::Abs(track->GetPt()) <  1.5)) {
+              (track->Pt()                <  1.5)) {
             UseClusters(track);
          }
        }
@@ -3723,7 +3723,7 @@ Int_t AliTRDtrackerHLT::FindClusters(Int_t sector, Int_t t0, Int_t t1
     tany *= -1.0;
   }
 
-  Double_t sy2       = ExpectedSigmaY2(x0,track->GetTgl(),track->GetPt());
+  Double_t sy2       = ExpectedSigmaY2(x0,track->GetTgl(),track->GetSignedPt());
   Double_t sz2       = ExpectedSigmaZ2(x0,track->GetTgl());
   Double_t road      = 15.0 * TMath::Sqrt(track->GetSigmaY2() + sy2);
   if (road > 6.0) {
@@ -3846,7 +3846,7 @@ Int_t AliTRDtrackerHLT::FindClusters(Int_t sector, Int_t t0, Int_t t1
       fHFindCl[iCheckPoint]->Fill(checkPoint[iCheckPoint]);
 
     if (checkPoint[3]) {
-      if (track->GetPt() > 0) fHMinYPos->Fill(minY);
+      if (track->GetSignedPt() > 0) fHMinYPos->Fill(minY);
       else fHMinYNeg->Fill(minY);
 
     fHMinD->Fill(minD[0], minD[1]);
index b4a712f..15f1a78 100644 (file)
@@ -146,7 +146,7 @@ Int_t AliITStrackMI::CorrectForMaterial(Double_t d, Double_t x0) {
   //------------------------------------------------------------------
   //d=x0=0.;
 
-  Double_t p2=(1.+ GetTgl()*GetTgl())/(Get1Pt()*Get1Pt());
+  Double_t p2=(1.+ GetTgl()*GetTgl())/(GetSigned1Pt()*GetSigned1Pt());
   Double_t et   = p2 + GetMass()*GetMass();
   Double_t beta2=p2/et;
   et = sqrt(et);  
@@ -158,8 +158,8 @@ Int_t AliITStrackMI::CorrectForMaterial(Double_t d, Double_t x0) {
      //Double_t theta2=1.0259e-6*14*14/28/(beta2*p2)*TMath::Abs(d)*9.36*2.33;
      Cov(5) += theta2*(1.- GetSnp()*GetSnp())*(1. + GetTgl()*GetTgl());
      Cov(9) += theta2*(1. + GetTgl()*GetTgl())*(1. + GetTgl()*GetTgl());
-     Cov(13)+= theta2*GetTgl()*Get1Pt()*(1. + GetTgl()*GetTgl());
-     Cov(14)+= theta2*GetTgl()*Get1Pt()*GetTgl()*Get1Pt();
+     Cov(13)+= theta2*GetTgl()*GetSigned1Pt()*(1. + GetTgl()*GetTgl());
+     Cov(14)+= theta2*GetTgl()*GetSigned1Pt()*GetTgl()*GetSigned1Pt();
   }
 
   //Energy losses************************
@@ -177,7 +177,7 @@ Int_t AliITStrackMI::CorrectForMaterial(Double_t d, Double_t x0) {
      //}
 
      Par(4)*=(1.- et/p2*dE);
-     Double_t delta44 = (dE*Get1Pt()*et/p2);
+     Double_t delta44 = (dE*GetSigned1Pt()*et/p2);
      delta44*=delta44;
      Cov(14)+= delta44/400.;
   }
index f604e70..f672c78 100644 (file)
@@ -97,8 +97,8 @@ Int_t AliITStrackV2::Compare(const TObject *o) const {
   // This function compares tracks according to the their curvature
   //-----------------------------------------------------------------
   AliITStrackV2 *t=(AliITStrackV2*)o;
-  //Double_t co=TMath::Abs(t->Get1Pt());
-  //Double_t c =TMath::Abs(Get1Pt());
+  //Double_t co=OneOverPt();
+  //Double_t c =OneOverPt();
   Double_t co=t->GetSigmaY2()*t->GetSigmaZ2();
   Double_t c =GetSigmaY2()*GetSigmaZ2();
   if (c>co) return 1;
@@ -332,7 +332,7 @@ Bool_t AliITStrackV2::Improve(Double_t x0,Double_t xyz[3],Double_t ers[3]) {
 
   Double_t dy = Par(0) - yv, dz = Par(1) - zv;
   Double_t r2=GetX()*GetX() + dy*dy;
-  Double_t p2=(1.+ GetTgl()*GetTgl())/(Get1Pt()*Get1Pt());
+  Double_t p2=(1.+ GetTgl()*GetTgl())/(GetSigned1Pt()*GetSigned1Pt());
   Double_t beta2=p2/(p2 + GetMass()*GetMass());
   x0*=TMath::Sqrt((1.+ GetTgl()*GetTgl())/(1.- GetSnp()*GetSnp()));
   Double_t theta2=14.1*14.1/(beta2*p2*1e6)*x0;
index d3b4fc4..2304a27 100644 (file)
@@ -380,7 +380,7 @@ Int_t AliITStrackerMI::Clusters2Tracks(AliESDEvent *event) {
          delete t;
          continue;
        }
-       if (TMath::Abs(1/t->Get1Pt())<AliITSReconstructor::GetRecoParam()->GetMinPtForProlongation()) {
+       if (t->Pt()<AliITSReconstructor::GetRecoParam()->GetMinPtForProlongation()) {
          delete t;
          continue;
        }
@@ -2148,7 +2148,7 @@ Double_t AliITStrackerMI::GetMatchingChi2(AliITStrackMI * track1, AliITStrackMI
   vec(1,0)=track1->GetZ()   - track3.GetZ();
   vec(2,0)=track1->GetSnp() - track3.GetSnp();
   vec(3,0)=track1->GetTgl() - track3.GetTgl();
-  vec(4,0)=track1->Get1Pt() - track3.Get1Pt();
+  vec(4,0)=track1->GetSigned1Pt() - track3.GetSigned1Pt();
   //
   TMatrixD cov(5,5);
   cov(0,0) = track1->GetSigmaY2()+track3.GetSigmaY2();
@@ -2262,7 +2262,7 @@ Double_t AliITStrackerMI::GetInterpolatedChi2(AliITStrackMI * forwardtrack, AliI
     npoints++;
   }
   if (npoints>1) return 
-                  TMath::Max(TMath::Abs(0.3*forwardtrack->Get1Pt())-0.5,0.)+
+                   TMath::Max(0.3*forwardtrack->OneOverPt()-0.5,0.)+
                   //2*forwardtrack->fNUsed+
                   res/TMath::Max(double(npoints-forwardtrack->GetNSkipped()),
                                  1./(1.+forwardtrack->GetNSkipped()));
@@ -2486,7 +2486,7 @@ AliITStrackMI *  AliITStrackerMI::GetBest2Tracks(Int_t trackID1, Int_t trackID2,
   //
   AliITStrackMI * track10=(AliITStrackMI*) arr1->UncheckedAt(0);
   AliITStrackMI * track20=(AliITStrackMI*) arr2->UncheckedAt(0);
-  if (TMath::Abs(1./track10->Get1Pt())>0.5+TMath::Abs(1/track20->Get1Pt())) return track10;
+  if (track10->Pt()>0.5+track20->Pt()) return track10;
 
   for (Int_t itrack=0;itrack<entries1;itrack++){
     AliITStrackMI * track=(AliITStrackMI*) arr1->UncheckedAt(itrack);
@@ -2549,11 +2549,11 @@ AliITStrackMI *  AliITStrackerMI::GetBest2Tracks(Int_t trackID1, Int_t trackID2,
   //
   w1 = (d2/(d1+d2)+ 2*s2/(s1+s2)+
        +s2/(s1+s2)*0.5*(chi22+2.)/(chi21+chi22+4.)
-       +1.*TMath::Abs(1./track10->Get1Pt())/(TMath::Abs(1./track10->Get1Pt())+TMath::Abs(1./track20->Get1Pt()))
+       +1.*track10->Pt()/(track10->Pt()+track20->Pt())
        );
   w2 = (d1/(d1+d2)+ 2*s1/(s1+s2)+
        s1/(s1+s2)*0.5*(chi21+2.)/(chi21+chi22+4.)
-       +1.*TMath::Abs(1./track20->Get1Pt())/(TMath::Abs(1./track10->Get1Pt())+TMath::Abs(1./track20->Get1Pt()))
+       +1.*track20->Pt()/(track10->Pt()+track20->Pt())
        );
 
   Double_t sumw = w1+w2;
index f05c54b..c9a517b 100644 (file)
@@ -110,7 +110,7 @@ void   AliITSImpactParametersPP(Int_t evFirst=0,Int_t evLast=0) {
       itstrack = (AliITStrackV2*)tArray.At(i);
 
       // pt
-      pt = 1./TMath::Abs(itstrack->Get1Pt());
+      pt = itstrack->Pt();
 
       // primary vertex from other tracks in event
       Bool_t goodVtx = kFALSE;
index 6413938..d05973b 100644 (file)
@@ -38,7 +38,7 @@ AliKFParticle::AliKFParticle( const AliExternalTrackParam &track, Int_t PID )
   Double_t energy = TMath::Sqrt( mass*mass + fP[3]*fP[3] + fP[4]*fP[4] + fP[5]*fP[5]);
   fP[6] = energy;
   fP[7] = 0;
-  fQ = (track.Get1Pt() >0 ) ?1 :-1;
+  fQ = (track.GetSigned1Pt() >0 ) ?1 :-1; // fQ = track->GetSign() would do the same thing
   fNDF = 0;
   fChi2 = 0;
   fAtProductionVertex = 0;
index e05775f..f64062b 100644 (file)
@@ -168,7 +168,7 @@ void RunV0(  AliESDEvent *event )
       if ( !( pTrack->GetStatus()&AliESDtrack::kITSrefit ) ) continue;
       Int_t indi[12];
       if( pTrack->GetITSclusters(indi) <5 ) continue;
-      Int_t PDG = ( pTrack->Get1Pt() <0 ) ?321 :211;
+      Int_t PDG = ( pTrack->GetSigned1Pt() <0 ) ?321 :211;
 
       //* take MC PDG  
       { 
index d604b06..01e0454 100644 (file)
@@ -616,7 +616,7 @@ Int_t TPCCmpTr::TreeGenLoop(Int_t eventNr)
        //if multifound find track with biggest pt - returning tracks - loos energy
        //
        if (fIndexRecTracks[fMCInfo->fLabel*4]+1){
-         Double_t pt = TMath::Abs(1/fTPCTrack->Get1Pt());
+         Double_t pt = fTPCTrack-Pt();
          for (Int_t i=1;i<4;i++){
            if (fIndexRecTracks[fMCInfo->fLabel*4+i]>=0){
              AliTPCtrack * track = (AliTPCtrack*)fTracks->UncheckedAt(fIndexRecTracks[fMCInfo->fLabel*4+i]);
index db549fb..4b2ef4e 100644 (file)
@@ -1345,7 +1345,7 @@ void AliTPCComparisonDraw::ResPt()
   p2->cd();
   p2->Draw();
   TH1F* his2 =  new TH1F("Ptresolution","Ptresolution",40,-5.,5.);
-  fTree->Draw("100.*(abs(1./fTPCTrack.Get1Pt())-MC.fTrackRef.Pt())/MC.fTrackRef.Pt()>>Ptresolution","MC.fRowsWithDigits>100&&RC.fTPCTrack.fN>50&&RC.fMultiple==1"+cteta1+cpos1+cprim);
+  fTree->Draw("100.*(fTPCTrack.Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt()>>Ptresolution","MC.fRowsWithDigits>100&&RC.fTPCTrack.fN>50&&RC.fMultiple==1"+cteta1+cpos1+cprim);
   AliLabelAxes(his2, "#Delta p_{t} / p_{t} [%]", "entries");
   his2->Fit("gaus");
   his2->Draw();
@@ -1464,7 +1464,7 @@ TH1F * AliTPCComparisonDraw::ResPtvsPt(const char* selection, const char * quali
   Double_t maxRes = 10.;
   TH2F* hRes2 = new TH2F("hRes2", "residuals", nBins, bins, nBinsRes, -maxRes, maxRes);
   
-  fTree->Draw("100.*(abs(1./fTPCTrack.Get1Pt())-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
+  fTree->Draw("100.*(fTPCTrack.Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
 
   TH1F* hMean=0;
   TH1F* hRes = CreateResHisto(hRes2, &hMean);
@@ -1484,7 +1484,7 @@ TH1F * AliTPCComparisonDraw::MeanPtvsPt(const char* selection, const char * qual
   Double_t maxRes = 10.;
   TH2F* hRes2 = new TH2F("hRes2", "residuals", nBins, bins, nBinsRes, -maxRes, maxRes);
   
-  fTree->Draw("100.*(1./fTPCTrack.Get1Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
+  fTree->Draw("100.*(fTPCTrack.Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
 
   TH1F* hMean=0;
   TH1F* hRes = CreateResHisto(hRes2, &hMean);
index 065e6ca..77b7b39 100644 (file)
@@ -471,11 +471,11 @@ Int_t AliTPCseed::Compare(const TObject *o) const {
   }
   else {
     Float_t f2 =1;
-    f2 = 1-20*TMath::Sqrt(t->GetSigma1Pt2())/(TMath::Abs(t->Get1Pt())+0.0066);
+    f2 = 1-20*TMath::Sqrt(t->GetSigma1Pt2())/(t->OneOverPt()+0.0066);
     if (t->fBConstrain) f2=1.2;
 
     Float_t f1 =1;
-    f1 = 1-20*TMath::Sqrt(GetSigma1Pt2())/(TMath::Abs(Get1Pt())+0.0066);
+    f1 = 1-20*TMath::Sqrt(GetSigma1Pt2())/(OneOverPt()+0.0066);
 
     if (fBConstrain)   f1=1.2;
  
index f73f238..872afa2 100644 (file)
@@ -186,8 +186,8 @@ Int_t AliTPCtrack::Compare(const TObject *o) const {
   // This function compares tracks according to the their curvature
   //-----------------------------------------------------------------
   AliTPCtrack *t=(AliTPCtrack*)o;
-  //Double_t co=TMath::Abs(t->Get1Pt());
-  //Double_t c =TMath::Abs(Get1Pt());
+  //Double_t co=t->OneOverPt();
+  //Double_t c = OneOverPt();
   Double_t co=t->GetSigmaY2()*t->GetSigmaZ2();
   Double_t c =GetSigmaY2()*GetSigmaZ2();
   if (c>co) return 1;
index 1995916..5c1219a 100644 (file)
@@ -232,7 +232,7 @@ Int_t AliTPCtracker::FollowProlongation(AliTPCseed& t, Int_t rf) {
     UInt_t index=0;
     Double_t maxchi2=kMaxCHI2;
     const AliTPCRow &krow=fSectors[s][nr];
-    Double_t pt=1./t.Get1Pt();
+    Double_t pt=t.GetSignedPt();
     Double_t sy2=AliTPCcluster::SigmaY2(t.GetX(),t.GetTgl(),pt);
     Double_t sz2=AliTPCcluster::SigmaZ2(t.GetX(),t.GetTgl());
     Double_t road=4.*sqrt(t.GetSigmaY2() + sy2), y=t.GetY(), z=t.GetZ();
@@ -370,7 +370,7 @@ Int_t AliTPCtracker::FollowBackProlongation
     AliTPCcluster *cl=0;
     Int_t index=0;
     Double_t maxchi2=kMaxCHI2;
-    Double_t pt=1./seed.Get1Pt();
+    Double_t pt=seed.GetSignedPt();
     Double_t sy2=AliTPCcluster::SigmaY2(seed.GetX(),seed.GetTgl(),pt);
     Double_t sz2=AliTPCcluster::SigmaZ2(seed.GetX(),seed.GetTgl());
     Double_t road=4.*sqrt(seed.GetSigmaY2() + sy2), z=seed.GetZ();
index 3f4c088..022fa1e 100644 (file)
@@ -1259,7 +1259,7 @@ void AliTPCtrackerParam::CompareTPCtracks(
       cmptrk.eta = part->Eta();
       cmptrk.r = TMath::Sqrt(part->Vx()*part->Vx()+part->Vy()*part->Vy());
       
-      cmptrk.pt   = 1/TMath::Abs(geatrack->Get1Pt());
+      cmptrk.pt   = geatrack->Pt();
       cmptrk.cosl = TMath::Cos(TMath::ATan(geatrack->GetTgl()));
       cmptrk.p    = cmptrk.pt/cmptrk.cosl;
     
@@ -1299,7 +1299,7 @@ void AliTPCtrackerParam::CompareTPCtracks(
       cmptrk.dP2 = kaltrack->GetSnp()-geatrack->GetSnp();
       cmptrk.dP3 = kaltrack->GetTgl()-geatrack->GetTgl();
       cmptrk.dP4 = kaltrack->GetC()-geatrack->GetC();
-      cmptrk.dpt = 1/kaltrack->Get1Pt()-1/geatrack->Get1Pt();
+      cmptrk.dpt = 1/kaltrack->GetSigned1Pt()-1/geatrack->GetSigned1Pt();
     
       // get covariance matrix
       // beware: lines 3 and 4 in the matrix are inverted!
@@ -1460,7 +1460,7 @@ void AliTPCtrackerParam::CookTrack(Double_t pt,Double_t eta) {
 
   // get P and Cosl from track
   cosl = TMath::Cos(TMath::ATan(fTrack.GetTgl()));
-  p    = 1./TMath::Abs(fTrack.Get1Pt())/cosl;
+  p    = fTrack.Pt()/cosl;
 
   trkKine[0] = p;
 
index 4d32288..60fc571 100644 (file)
@@ -404,7 +404,7 @@ Bool_t AliTPCcalibTracks::AcceptTrack(AliTPCseed * track){
   if (track->GetNumberOfClusters() < kMinClusters) return kFALSE;
   Float_t ratio = track->GetNumberOfClusters() / (track->GetNFoundable() + 1.);
   if (ratio < kMinRatio) return kFALSE;
-  Float_t mpt = track->Get1Pt();
+  Float_t mpt = track->GetSigned1Pt();
   if (TMath::Abs(mpt) > kMax1pt) return kFALSE;
   //if (TMath::Abs(track->GetZ())>240.) return kFALSE;
   //if (TMath::Abs(track->GetZ())<10.) return kFALSE;
index 3784800..558dc5f 100644 (file)
@@ -807,7 +807,7 @@ Bool_t AliTRDCalibraFillHisto::UpdateHistograms(AliTRDcluster *cl, AliTRDtrack *
   fPar1[time] = t->GetZ();
   fPar2[time] = t->GetSnp();
   fPar3[time] = t->GetTgl();
-  fPar4[time] = t->Get1Pt();
+  fPar4[time] = t->GetSigned1Pt();
 
   // Store the info bis of the tracklet
   Int_t *rowcol   = CalculateRowCol(cl);
index 7c06033..3213b82 100644 (file)
@@ -700,7 +700,7 @@ Bool_t AliTRDtrack::PropagateTo(Double_t xk, Double_t xx0, Double_t xrho)
   {
 
     // Energy losses************************
-    Double_t p2    = (1.0 + GetTgl()*GetTgl()) / (Get1Pt()*Get1Pt());
+    Double_t p2    = (1.0 + GetTgl()*GetTgl()) / (GetSigned1Pt()*GetSigned1Pt());
     Double_t beta2 = p2 / (p2 + GetMass()*GetMass());
     if (beta2<1.0e-10 || (5940.0 * beta2/(1.0 - beta2 + 1.0e-10) - beta2) < 0.0) {
       return kFALSE;
@@ -853,6 +853,7 @@ Int_t AliTRDtrack::UpdateMI(AliTRDcluster *c, Double_t chisq, Int_t index, Doubl
 
   return kTRUE;
 
+
 }                     
 
 // //_____________________________________________________________________________
@@ -949,6 +950,7 @@ Int_t AliTRDtrack::UpdateMI(AliTRDcluster *c, Double_t chisq, Int_t index, Doubl
 //   fCtt-=k30*oldty+k31*oldtz;
 //   fCct-=k40*oldty+k41*oldtz;
 //   //
+
 //   fCcc-=k40*oldcy+k41*oldcz;                 
 //   //
   
index c6d1eb8..c258d16 100644 (file)
@@ -706,7 +706,7 @@ Int_t AliTRDtracker::PropagateBack(AliESDEvent *event)
     /**/
     // inter-tracks competition ???
     if ((TMath::Abs(track->GetC() - p4) / TMath::Abs(p4) < 0.2) || 
-        (TMath::Abs(track->GetPt())                      > 0.8)) {
+        (track->Pt()                                     > 0.8)) {
       
       fHBackfit->Fill(4);
       
@@ -726,7 +726,7 @@ Int_t AliTRDtracker::PropagateBack(AliESDEvent *event)
         // Sign only gold tracks
        if (track->GetChi2() / track->GetNumberOfClusters() < 4) {   
          if ((seed->GetKinkIndex(0)      ==   0) &&
-              (TMath::Abs(track->GetPt()) <  1.5)) {
+              (track->Pt()                <  1.5)) {
             UseClusters(track);
          }
        }
@@ -3514,7 +3514,7 @@ Int_t AliTRDtracker::FindClusters(Int_t sector, Int_t t0, Int_t t1
     tany *= -1.0;
   }
 
-  Double_t sy2       = ExpectedSigmaY2(x0,track->GetTgl(),track->GetPt());
+  Double_t sy2       = ExpectedSigmaY2(x0,track->GetTgl(),track->GetSignedPt());
   Double_t sz2       = ExpectedSigmaZ2(x0,track->GetTgl());
   Double_t road      = 15.0 * TMath::Sqrt(track->GetSigmaY2() + sy2);
   if (road > 6.0) {
@@ -3637,7 +3637,7 @@ Int_t AliTRDtracker::FindClusters(Int_t sector, Int_t t0, Int_t t1
       fHFindCl[iCheckPoint]->Fill(checkPoint[iCheckPoint]);
 
     if (checkPoint[3]) {
-      if (track->GetPt() > 0) fHMinYPos->Fill(minY);
+      if (track->GetSignedPt() > 0) fHMinYPos->Fill(minY);
       else fHMinYNeg->Fill(minY);
 
     fHMinD->Fill(minD[0], minD[1]);
index d45064a..356dd08 100644 (file)
@@ -260,7 +260,7 @@ void AliTRDtrackingAnalysis::DrawResolutionPt(int startEvent, int stopEvent)
 
       if (!(status & AliESDtrack::kTRDout)) continue;
       if (!(status & AliESDtrack::kTRDrefit)) continue;
-      if (TMath::Abs(esdTrack->GetOuterParam()->GetPt()) < 1.0) continue;
+      if (esdTrack->GetOuterParam()->Pt() < 1.0) continue;
 
       int ch=0;
       while(param->GetX() > fGeo->GetTime0(ch)+2) ch++;
@@ -287,7 +287,7 @@ void AliTRDtrackingAnalysis::DrawResolutionPt(int startEvent, int stopEvent)
        AliExternalTrackParam out(*param);
        out.PropagateTo(aRef->LocalX(),bz);
        
-       double dp = aRef->Pt() + out.GetPt();
+       double dp = aRef->Pt() + out.GetSignedPt();
        double dy = 10. * (aRef->LocalY() - out.GetY()); // in mm
 
        fDeltaPt->Fill(100. * dp / aRef->Pt());
index dc7030b..c670102 100644 (file)
@@ -651,7 +651,7 @@ void MakeAlias(){
   comp.fTree->SetAlias("dt","(abs(MC.fPdg)==2212)*dtp+(abs(MC.fPdg)==211)*dtpi+(abs(MC.fPdg)==321)*dtk"); //delta time
   comp.fTree->SetAlias("Beta","sqrt(TR.P()**2/(TR.P()**2+MC.fMass**2))");
   comp.fTree->SetAlias("dtsigma","10+(1-Beta)*280");
-  comp.fTree->SetAlias("dptrel","(abs(TRD0.GetPt())-MC.fTOFReferences[0].Pt())/MC.fTOFReferences[0].Pt()");
+  comp.fTree->SetAlias("dptrel","(abs(TRD0.GetSignedPt())-MC.fTOFReferences[0].Pt())/MC.fTOFReferences[0].Pt()");
   //
   comp.fTree->SetAlias("dphi","atan2(TR.fY,TR..fX)-atan2(TR.fPy,TR..fPx)");
   comp.fTree->SetAlias("dphi0","acos((TR.fX*TR.fPx+TR.fY*TR.fPy)/(TR.Pt()*TR.R()))");
index 8ee54af..7b9982d 100644 (file)
@@ -241,7 +241,7 @@ Int_t AliTRDComparisonV2
 
         hc->Fill(esd->GetTRDclusters(0));
 
-        Float_t pt_1=TMath::Abs(out.Get1Pt());
+        Float_t pt_1=out.OneOverPt();
         hpt->Fill((pt_1 - 1/ptg)/(1/ptg)*100.);
 
         Float_t y=out.GetY();
index c1b82d3..6d2ddf4 100644 (file)
@@ -453,7 +453,7 @@ void AliTRDbackTrackAnalysis() {
 
     hfound->Fill(p->Pt());
 
-    Pt = TMath::Abs(rt->GetPt()); 
+    Pt = rt->Pt(); 
     Double_t cc = TMath::Abs(rt->GetSigmaC2()); 
     mct->GetPxPyPzXYZ(mcPx,mcPy,mcPz,x,y,z,-1);
     rt->GetPxPyPz(Px,Py,Pz);      
index 682b3f9..abb2fd4 100644 (file)
@@ -489,7 +489,7 @@ void AliTRDseedsAnalysis() {
                 label, rt->GetSeedLabel());
          return;
        }
-       Pt = TMath::Abs(rt->GetPt()); 
+       Pt = rt->GetPt(); 
        Double_t cc = TMath::Abs(rt->GetSigmaC2()); 
        mct->GetPxPyPzXYZ(mcPx,mcPy,mcPz,x,y,z,-1);
        rt->GetPxPyPz(Px,Py,Pz);