PropagateToDCA in case of track and vertex (M.Ivanov)

[u/mrichter/AliRoot.git] / STEER / AliESDtrack.cxx

diff --git a/STEER/AliESDtrack.cxx b/STEER/AliESDtrack.cxx
index d3629a02f2486694525e507647a8be2f2fa30866..44755fb2f59bc259f8bd170f0f49746a6e361ecc 100644 (file)
--- a/STEER/AliESDtrack.cxx
+++ b/STEER/AliESDtrack.cxx
@@ -409,6 +409,10 @@ Bool_t AliESDtrack::UpdateTrackParams(const AliKalmanTrack *t, ULong_t flags){
     if (!fIp) fIp=new AliExternalTrackParam(*t);
     else fIp->Set(*t);
   case kTPCout:
+    if (flags & kTPCout){
+      if (!fOp) fOp=new AliExternalTrackParam(*t);
+      else fOp->Set(*t);
+    }
     fTPCncls=t->GetNumberOfClusters();    
     fTPCchi2=t->GetChi2();
     
@@ -843,7 +847,7 @@ Bool_t AliESDtrack::RelateToVertex
   //Propagate to the DCA
   Double_t crv=0.299792458e-3*b*GetParameter()[4];
   Double_t tgfv=-(crv*x - snp)/(crv*y + TMath::Sqrt(1.-snp*snp));
-  sn=tgfv/TMath::Sqrt(1.+ tgfv*tgfv); cs=TMath::Sqrt(1.+ sn*sn);
+  sn=tgfv/TMath::Sqrt(1.+ tgfv*tgfv); cs=TMath::Sqrt(1.- sn*sn);
 
   x = xv*cs + yv*sn;
   yv=-xv*sn + yv*cs; xv=x;