From: hristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Wed, 27 Oct 2004 13:14:02 +0000 (+0000)
Subject: Storing impact parameters in the ESD track (Yu.Belikov)
X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=399fb957e8f8efa9c1863fe99eab1dbf066bc8ea

Storing impact parameters in the ESD track (Yu.Belikov)
---

diff --git a/ITS/AliITStrackerMI.cxx b/ITS/AliITStrackerMI.cxx
index 47af8e6863b..e7f69d5b0b1 100644
--- a/ITS/AliITStrackerMI.cxx
+++ b/ITS/AliITStrackerMI.cxx
@@ -462,6 +462,9 @@ Int_t AliITStrackerMI::RefitInward(AliESD *event) {
 
 	 if (fTrackToFollow.Propagate(fv+a,xv)) {
             fTrackToFollow.UpdateESDtrack(AliESDtrack::kITSrefit);
+            Float_t d=fTrackToFollow.GetD(GetX(),GetY());
+            Float_t z=fTrackToFollow.GetZ()-GetZ();
+            fTrackToFollow.GetESDtrack()->SetImpactParameters(d,z);
             //UseClusters(&fTrackToFollow);
             {
             AliITSclusterV2 c; c.SetY(yv); c.SetZ(GetZ());
diff --git a/ITS/AliITStrackerV2.cxx b/ITS/AliITStrackerV2.cxx
index e57517be4ce..bb19687bc31 100644
--- a/ITS/AliITStrackerV2.cxx
+++ b/ITS/AliITStrackerV2.cxx
@@ -376,6 +376,9 @@ Int_t AliITStrackerV2::RefitInward(AliESD *event) {
 
 	 if (fTrackToFollow.Propagate(fv+a,xv)) {
             fTrackToFollow.UpdateESDtrack(AliESDtrack::kITSrefit);
+            Float_t d=fTrackToFollow.GetD(GetX(),GetY());
+            Float_t z=fTrackToFollow.GetZ()-GetZ();
+            fTrackToFollow.GetESDtrack()->SetImpactParameters(d,z);
             UseClusters(&fTrackToFollow);
             {
             AliITSclusterV2 c; c.SetY(yv); c.SetZ(GetZ());
diff --git a/STEER/AliESDtrack.cxx b/STEER/AliESDtrack.cxx
index eb33fc742bc..8fd19354fcb 100644
--- a/STEER/AliESDtrack.cxx
+++ b/STEER/AliESDtrack.cxx
@@ -33,6 +33,8 @@ AliESDtrack::AliESDtrack() :
 fFlags(0),
 fLabel(0),
 fTrackLength(0),
+fD(0),
+fZ(0),
 fStopVertex(0),
 fRalpha(0),
 fRx(0),
@@ -116,6 +118,7 @@ AliESDtrack::AliESDtrack(const AliESDtrack& track):TObject(track){
   fFlags = track.fFlags;
   fLabel =track.fLabel;
   fTrackLength =track.fTrackLength;
+  fD=track.fD; fZ=track.fZ;
   for (Int_t i=0;i<kSPECIES;i++) fTrackTime[i] =track.fTrackTime[i];
   for (Int_t i=0;i<kSPECIES;i++)  fR[i] =track.fR[i];
   fStopVertex =track.fStopVertex;
diff --git a/STEER/AliESDtrack.h b/STEER/AliESDtrack.h
index 55322e2ceea..81f21c3bb50 100644
--- a/STEER/AliESDtrack.h
+++ b/STEER/AliESDtrack.h
@@ -38,6 +38,7 @@ public:
   void SetStatus(ULong_t flags) {fFlags|=flags;}
   void ResetStatus(ULong_t flags) {fFlags&=~flags;}
   Bool_t UpdateTrackParams(const AliKalmanTrack *t, ULong_t flags);
+  void SetImpactParameters(Float_t xy,Float_t z) {fD=xy; fZ=z;}
   void SetIntegratedLength(Double_t l) {fTrackLength=l;}
   void SetIntegratedTimes(const Double_t *times);
   void SetESDpid(const Double_t *p);
@@ -53,6 +54,7 @@ public:
   Bool_t GetPxPyPzAt(Double_t x, Double_t p[3]) const;
   Bool_t GetXYZAt(Double_t x, Double_t r[3]) const;
 
+  void GetImpactParameters(Float_t &xy,Float_t &z) const {xy=fD; z=fZ;}
   Double_t GetIntegratedLength() const {return fTrackLength;}
   void GetIntegratedTimes(Double_t *times) const;
   Double_t GetMass() const;
@@ -186,7 +188,9 @@ protected:
   ULong_t   fFlags;        // Reconstruction status flags 
   Int_t     fLabel;        // Track label
   Int_t     fID;           // Unique ID of the track
-  Float_t   fTrackLength;         // Track length
+  Float_t   fTrackLength;  // Track length
+  Float_t   fD;            // Impact parameter in XY-plane
+  Float_t   fZ;            // Impact parameter in Z 
   Float_t   fTrackTime[kSPECIES]; // TOFs estimated by the tracking
   Float_t   fR[kSPECIES];         // combined "detector response probability"
 
@@ -274,7 +278,7 @@ protected:
   Float_t fRICHsignal;     // detector's PID signal (beta for RICH)
   Float_t fRICHr[kSPECIES];// "detector response probabilities" (for the PID)
   	
-  ClassDef(AliESDtrack,10)  //ESDtrack 
+  ClassDef(AliESDtrack,11)  //ESDtrack 
 };
 
 #endif