Fix in GetXYZAt(...) : fMomentum is NOT px,py,pz

[u/mrichter/AliRoot.git] / STEER / AOD / AliAODTrack.cxx

diff --git a/STEER/AOD/AliAODTrack.cxx b/STEER/AOD/AliAODTrack.cxx
index bbda27e1a2761602142470f953b0b81b377801f9..a33d5baf584da800732e65c09f067e25a93ff86e 100644 (file)
--- a/STEER/AOD/AliAODTrack.cxx
+++ b/STEER/AOD/AliAODTrack.cxx
@@ -905,7 +905,7 @@ Bool_t AliAODTrack::GetXYZAt(Double_t x, Double_t b, Double_t *r) const
   Double_t radPos2 = fPosition[0]*fPosition[0]+fPosition[1]*fPosition[1];  
   Double_t radMax  = 45.; // approximately ITS outer radius
   if (radPos2 < radMax*radMax) { // inside the ITS     
-     alpha = TMath::ATan2(fMomentum[1],fMomentum[0]);
+    alpha = fMomentum[1]; //TMath::ATan2(fMomentum[1],fMomentum[0]); // fMom is pt,phi,theta!
   } else { // outside the ITS
      Float_t phiPos = TMath::Pi()+TMath::ATan2(-fPosition[1], -fPosition[0]);
      alpha = 
@@ -927,7 +927,7 @@ Bool_t AliAODTrack::GetXYZAt(Double_t x, Double_t b, Double_t *r) const
   
   // Get the vertex of origin and the momentum
   TVector3 ver(fPosition[0],fPosition[1],fPosition[2]);
-  TVector3 mom(fMomentum[0],fMomentum[1],fMomentum[2]);
+  TVector3 mom(Px(),Py(),Pz());
   //
   // avoid momenta along axis
   if (TMath::Abs(mom[0])<kSafe) mom[0] = TMath::Sign(kSafe*TMath::Abs(mom[1]), mom[0]);
@@ -958,7 +958,6 @@ Bool_t AliAODTrack::GetXYZAt(Double_t x, Double_t b, Double_t *r) const
   r[0] = x;
   r[1] = param0 + dx*(f1+f2)/(r1+r2);
   r[2] = param1 + dx*(r2 + f2*(f1+f2)/(r1+r2))*param3;//Thanks to Andrea & Peter
-
   return Local2GlobalPosition(r,alpha);
 }