3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
7 //-------------------------------------------------------
10 // Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
12 // The track parameterization is fixed in the following way:
13 // param0: local Y-coordinate of a track (cm)
14 // param1: local Z-coordinate of a track (cm)
15 // param2: local sine of the track momentum azimuth angle
16 // param3: tangent of the track momentum dip angle
17 // param4: 1/pt (1/(GeV/c))
19 //-------------------------------------------------------
21 #include <AliKalmanTrack.h>
24 #include "AliTPCreco.h"
25 #include "AliExternalTrackParam.h"
28 class TTreeSRedirector;
29 //_____________________________________________________________________________
30 class AliTPCtrack : public AliKalmanTrack {
33 AliTPCtrack(Double_t x, Double_t alpha, const Double_t p[5],
34 const Double_t cov[15], Int_t index);
35 AliTPCtrack(const AliESDtrack& t, TTreeSRedirector *pcstream=0);
36 AliTPCtrack(const AliTPCtrack& t);
37 AliTPCtrack& operator=(const AliTPCtrack& o);
38 virtual ~AliTPCtrack() {}
40 Int_t Compare(const TObject *o) const;
42 void SetdEdx(Double_t dedx) {fdEdx=dedx;}
43 Double_t GetdEdx() const {return fdEdx;}
44 Double_t GetPIDsignal() const {return GetdEdx();}
46 Int_t GetClusterIndex(Int_t i) const {return fIndex[i];}
47 void SetClusterIndex(Int_t i, Int_t idx) {fIndex[i]=idx;}
49 Double_t GetC() const {return AliExternalTrackParam::GetC(GetBz());}
50 Double_t GetC(Double_t b) const {return AliExternalTrackParam::GetC(b);}
52 Double_t GetPredictedChi2(const AliCluster *cluster) const;
53 Bool_t PropagateTo(Double_t xr, Double_t rho=0.9e-3, Double_t x0=28.94);
54 Bool_t Update(const AliCluster *c, Double_t chi2, Int_t i);
55 Bool_t Rotate(Double_t alpha) {
56 return AliExternalTrackParam::Rotate(GetAlpha()+alpha);
59 Bool_t PropagateToVertex(const AliESDVertex *v,
60 Double_t rho=1.2e-3, Double_t x0=36.66);
61 void ResetClusters() {SetNumberOfClusters(0); SetChi2(0.);}
62 void UpdatePoints();//update points
63 Float_t* GetPoints() {return fPoints;}
65 Float_t Density(Int_t row0, Int_t row1); //calculate cluster density
66 Float_t Density2(Int_t row0, Int_t row1); //calculate cluster density
67 Double_t GetD(Double_t x=0, Double_t y=0) const {
68 return AliExternalTrackParam::GetD(x,y,GetBz());
70 AliExternalTrackParam & GetReference(){ return fReference;}
71 void UpdateReference(){ new (&fReference) AliExternalTrackParam(*this);}
72 Int_t GetKinkIndex(Int_t i) const{ return fKinkIndexes[i];}
73 Int_t* GetKinkIndexes() { return fKinkIndexes;}
74 Int_t GetV0Index(Int_t i) const{ return fV0Indexes[i];}
75 Int_t* GetV0Indexes() { return fV0Indexes;}
77 void SetFirstPoint(Int_t f) {fFirstPoint=f;}
78 void SetLastPoint(Int_t f) {fLastPoint=f;}
79 void SetRemoval(Int_t f) {fRemoval=f;}
80 void SetLab2(Int_t f) {fLab2=f;}
81 void SetKinkIndex(Int_t i, Int_t idx) {fKinkIndexes[i]=idx;}
82 void SetBConstrain(Bool_t b) {fBConstrain=b;}
83 void SetNFoundable(Int_t n) {fNFoundable=n;}
84 void SetNShared(Int_t s) {fNShared=s;}
86 Int_t GetFirstPoint() const {return fFirstPoint;}
87 Int_t GetLastPoint() const {return fLastPoint;}
88 Int_t GetRemoval() const {return fRemoval;}
89 Int_t GetLab2() const {return fLab2;}
90 Bool_t GetBConstrain() const {return fBConstrain;}
91 Int_t GetNShared() const {return fNShared;}
92 Int_t GetNFoundable() const {return fNFoundable;}
95 Double_t fdEdx; // dE/dx
97 Int_t fIndex[kMaxRow]; // indices of associated clusters
98 Float_t fPoints[4]; //first, max dens row end points of the track and max density
100 Float_t fSdEdx; // sigma of dedx
102 Int_t fNFoundable; //number of foundable clusters - dead zone taken to the account
103 Bool_t fBConstrain; // indicate seeding with vertex constrain
104 Int_t fLastPoint; // last cluster position
105 Int_t fFirstPoint; // first cluster position
106 Int_t fRemoval; // removal factor
107 Int_t fTrackType; // track type - 0 - normal - 1 - kink - 2 -V0 3- double found
108 Int_t fLab2; // index of corresponding track (kink, V0, double)
109 Int_t fNShared; // number of shared points
110 AliExternalTrackParam fReference; // track parameters at the middle of the chamber
111 Float_t fKinkPoint[12]; //radius, of kink, dfi and dtheta
112 Int_t fKinkIndexes[3]; // kink indexes - minus = mother + daughter
113 Int_t fV0Indexes[3]; // kink indexes - minus = mother + daughter
115 ClassDef(AliTPCtrack,4) // Time Projection Chamber reconstructed tracks