]>
Commit | Line | Data |
---|---|---|
1 | #ifndef ALITRDTRACK_H | |
2 | #define ALITRDTRACK_H | |
3 | ||
4 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
5 | * See cxx source for full Copyright notice */ | |
6 | ||
7 | #include <AliKalmanTrack.h> | |
8 | #include <TMath.h> | |
9 | ||
10 | #include "AliBarrelTrack.h" | |
11 | #include "AliTRDgeometry.h" | |
12 | #include "TVector2.h" | |
13 | ||
14 | class AliTRDcluster; | |
15 | class AliTPCtrack; | |
16 | class AliESDtrack; | |
17 | ||
18 | const unsigned kMAX_CLUSTERS_PER_TRACK=210; | |
19 | ||
20 | class AliTRDtrack : public AliKalmanTrack { | |
21 | ||
22 | // Represents reconstructed TRD track | |
23 | ||
24 | public: | |
25 | ||
26 | AliTRDtrack():AliKalmanTrack(){} | |
27 | AliTRDtrack(const AliTRDcluster *c, UInt_t index, const Double_t xx[5], | |
28 | const Double_t cc[15], Double_t xr, Double_t alpha); | |
29 | AliTRDtrack(const AliTRDtrack& t); | |
30 | AliTRDtrack(const AliKalmanTrack& t, Double_t alpha); | |
31 | AliTRDtrack(const AliESDtrack& t); | |
32 | ||
33 | Int_t Compare(const TObject *o) const; | |
34 | void CookdEdx(Double_t low=0.05, Double_t up=0.70); | |
35 | ||
36 | Double_t GetAlpha() const {return fAlpha;} | |
37 | Int_t GetSector() const { | |
38 | //if (fabs(fAlpha) < AliTRDgeometry::GetAlpha()/2) return 0; | |
39 | return Int_t(TVector2::Phi_0_2pi(fAlpha)/AliTRDgeometry::GetAlpha())%AliTRDgeometry::kNsect;} | |
40 | ||
41 | Double_t GetC() const {return fC;} | |
42 | Int_t GetClusterIndex(Int_t i) const {return fIndex[i];} | |
43 | Float_t GetClusterdQdl(Int_t i) const {return fdQdl[i];} | |
44 | ||
45 | void GetCovariance(Double_t cov[15]) const; | |
46 | Double_t GetdEdx() const {return fdEdx;} | |
47 | Double_t GetPIDsignal() const {return GetdEdx();} | |
48 | Double_t GetEta() const {return fE;} | |
49 | ||
50 | void GetExternalCovariance(Double_t cov[15]) const ; | |
51 | void GetExternalParameters(Double_t& xr, Double_t x[5]) const ; | |
52 | ||
53 | Double_t GetLikelihoodElectron() const { return fLhElectron; }; | |
54 | ||
55 | Double_t Get1Pt() const {return (1e-9*TMath::Abs(fC)/fC + fC)*GetConvConst(); } | |
56 | Double_t GetP() const { | |
57 | return TMath::Abs(GetPt())*sqrt(1.+GetTgl()*GetTgl()); | |
58 | } | |
59 | Double_t GetPredictedChi2(const AliTRDcluster*, Double_t h01) const ; | |
60 | Double_t GetPt() const {return 1./Get1Pt();} | |
61 | void GetPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const ; | |
62 | void GetGlobalXYZ(Double_t &x, Double_t &y, Double_t &z) const ; | |
63 | Int_t GetSeedLabel() const { return fSeedLab; } | |
64 | Double_t GetSigmaC2() const {return fCcc;} | |
65 | Double_t GetSigmaTgl2() const {return fCtt;} | |
66 | Double_t GetSigmaY2() const {return fCyy;} | |
67 | Double_t GetSigmaZ2() const {return fCzz;} | |
68 | Double_t GetSnp() const {return fX*fC - fE;} | |
69 | Double_t GetTgl() const {return fT;} | |
70 | Double_t GetX() const {return fX;} | |
71 | Double_t GetY() const {return fY;} | |
72 | Double_t GetZ() const {return fZ;} | |
73 | UInt_t * GetBackupIndexes() {return fIndexBackup;} | |
74 | UInt_t * GetIndexes() {return fIndex;} | |
75 | ||
76 | ||
77 | Double_t GetYat(Double_t xk) const { | |
78 | //----------------------------------------------------------------- | |
79 | // This function calculates the Y-coordinate of a track at the plane x=xk. | |
80 | // Needed for matching with the TOF (I.Belikov) | |
81 | //----------------------------------------------------------------- | |
82 | Double_t c1=fC*fX - fE, r1=TMath::Sqrt(1.- c1*c1); | |
83 | Double_t c2=fC*xk - fE, r2=TMath::Sqrt(1.- c2*c2); | |
84 | return fY + (xk-fX)*(c1+c2)/(r1+r2); | |
85 | } | |
86 | void SetStop(Bool_t stop) {fStopped=stop;} | |
87 | Bool_t GetStop() const {return fStopped;} | |
88 | ||
89 | Int_t PropagateTo(Double_t xr, Double_t x0=8.72, Double_t rho=5.86e-3); | |
90 | void ResetCovariance(); | |
91 | void ResetCovariance(Float_t mult); | |
92 | void ResetClusters() { SetChi2(0.); SetNumberOfClusters(0); } | |
93 | Int_t Rotate(Double_t angle); | |
94 | ||
95 | void SetdEdx(Float_t dedx) {fdEdx=dedx;} | |
96 | void SetLikelihoodElectron(Float_t l) { fLhElectron = l; }; | |
97 | ||
98 | void SetSampledEdx(Float_t q, Int_t i) { | |
99 | Double_t s=GetSnp(), t=GetTgl(); | |
100 | q*= TMath::Sqrt((1-s*s)/(1+t*t)); | |
101 | fdQdl[i]=q; | |
102 | } | |
103 | ||
104 | void SetSeedLabel(Int_t lab) { fSeedLab=lab; } | |
105 | ||
106 | Int_t Update(const AliTRDcluster* c, Double_t chi2, UInt_t i, | |
107 | Double_t h01); | |
108 | Int_t UpdateMI(const AliTRDcluster* c, Double_t chi2, UInt_t i, | |
109 | Double_t h01, Int_t plane); | |
110 | ||
111 | // | |
112 | void GetBarrelTrack(AliBarrelTrack *track); | |
113 | void AddNWrong() {fNWrong++;} | |
114 | ||
115 | Int_t GetNWrong() const {return fNWrong;} | |
116 | Int_t GetNRotate() const {return fNRotate;} | |
117 | // | |
118 | ||
119 | ||
120 | protected: | |
121 | ||
122 | Int_t fSeedLab; // track label taken from seeding | |
123 | Float_t fdEdx; // dE/dx | |
124 | ||
125 | Double_t fAlpha; // rotation angle | |
126 | Double_t fX; // running local X-coordinate of the track (time bin) | |
127 | Bool_t fStopped; // track stop indication | |
128 | ||
129 | Double_t fY; // Y-coordinate of the track | |
130 | Double_t fZ; // Z-coordinate of the track | |
131 | Double_t fE; // C*x0 | |
132 | Double_t fT; // tangent of the track momentum dip angle | |
133 | Double_t fC; // track curvature | |
134 | ||
135 | Double_t fCyy; // covariance | |
136 | Double_t fCzy, fCzz; // matrix | |
137 | Double_t fCey, fCez, fCee; // of the | |
138 | Double_t fCty, fCtz, fCte, fCtt; // track | |
139 | Double_t fCcy, fCcz, fCce, fCct, fCcc; // parameters | |
140 | ||
141 | UInt_t fIndex[kMAX_CLUSTERS_PER_TRACK]; // global indexes of clusters | |
142 | UInt_t fIndexBackup[kMAX_CLUSTERS_PER_TRACK]; //backup indexes of clusters - used in iterations | |
143 | Float_t fdQdl[kMAX_CLUSTERS_PER_TRACK]; // cluster amplitudes corrected | |
144 | // for track angles | |
145 | ||
146 | Float_t fLhElectron; // Likelihood to be an electron | |
147 | Int_t fNWrong; // number of wrong clusters | |
148 | Int_t fNRotate; | |
149 | ||
150 | ClassDef(AliTRDtrack,2) // TRD reconstructed tracks | |
151 | ||
152 | }; | |
153 | ||
154 | ||
155 | #endif |