]>
Commit | Line | Data |
---|---|---|
46d29e70 | 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 <TObject.h> | |
8 | ||
9 | class AliTRDcluster; | |
10 | ||
11 | class AliTRDtrack : public TObject { | |
12 | ||
13 | // Represents reconstructed TRD track | |
14 | ||
15 | Int_t fLab; // track label | |
16 | Double_t fChi2; // total chi2 value for the track | |
17 | Float_t fdEdx; // dE/dx | |
18 | ||
19 | Double_t fAlpha; // rotation angle | |
20 | Double_t fX; // running local X-coordinate of the track (time bin) | |
21 | ||
22 | Double_t fY; // Y-coordinate of the track | |
23 | Double_t fZ; // Z-coordinate of the track | |
24 | Double_t fC; // track curvature | |
25 | Double_t fE; // C*x0 | |
26 | Double_t fT; // tangent of the track dip angle | |
27 | ||
28 | Double_t fCyy; // covariance | |
29 | Double_t fCzy, fCzz; // matrix | |
30 | Double_t fCcy, fCcz, fCcc; // of the | |
31 | Double_t fCey, fCez, fCec, fCee; // track | |
32 | Double_t fCty, fCtz, fCtc, fCte, fCtt; // parameters | |
33 | ||
34 | Short_t fN; // number of clusters associated with the track | |
35 | UInt_t fIndex[200]; // global indexes of these clusters | |
36 | ||
37 | ||
38 | public: | |
39 | ||
40 | AliTRDtrack() { fN=0;} | |
41 | AliTRDtrack(UInt_t index, const Double_t xx[5], | |
42 | const Double_t cc[15], Double_t xr, Double_t alpha); | |
43 | AliTRDtrack(const AliTRDtrack& t); | |
44 | ||
45 | Int_t Compare(TObject *o); | |
46 | ||
47 | Double_t GetAlpha() const {return fAlpha;} | |
48 | Double_t GetC() const {return fC;} | |
49 | Int_t GetClusterIndex(Int_t i) const {return fIndex[i];} | |
50 | Double_t GetChi2() const {return fChi2;} | |
51 | void GetCovariance(Double_t cov[15]) const; | |
52 | Double_t GetdEdX() const {return fdEdx;} | |
53 | Double_t GetEta() const {return fE;} | |
54 | Int_t GetLabel() const {return fLab;} | |
55 | Int_t GetNclusters() const {return fN;} | |
56 | Double_t GetP() const { | |
57 | return TMath::Abs(GetPt())*sqrt(1.+GetTgl()*GetTgl()); | |
58 | } | |
59 | Double_t GetPredictedChi2(const AliTRDcluster*) const ; | |
60 | Double_t GetPt() const {return 0.3*0.2/GetC()/100;} | |
61 | void GetPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const ; | |
62 | Double_t GetSigmaC2() const {return fCcc;} | |
63 | Double_t GetSigmaTgl2() const {return fCtt;} | |
64 | Double_t GetSigmaY2() const {return fCyy;} | |
65 | Double_t GetSigmaZ2() const {return fCzz;} | |
66 | Double_t GetTgl() const {return fT;} | |
67 | Double_t GetX() const {return fX;} | |
68 | Double_t GetY() const {return fY;} // returns running Y | |
69 | Double_t GetZ() const {return fZ;} | |
70 | ||
71 | Bool_t IsSortable() const {return kTRUE;} | |
72 | ||
73 | Int_t PropagateTo(Double_t xr, | |
74 | Double_t x0=8.72,Double_t rho=5.86e-3,Double_t pm=0.139); | |
75 | void PropagateToVertex( | |
76 | Double_t x0=36.66,Double_t rho=1.2e-3,Double_t pm=0.139); | |
77 | Int_t Rotate(Double_t angle); | |
78 | ||
79 | void SetLabel(Int_t lab=0) {fLab=lab;} | |
80 | void SetdEdx(Float_t dedx) {fdEdx=dedx;} | |
81 | ||
82 | void Update(const AliTRDcluster* c, Double_t chi2, UInt_t i); | |
83 | ||
84 | ClassDef(AliTRDtrack,1) // TRD reconstructed tracks | |
85 | }; | |
86 | ||
87 | ||
88 | //----------------------------------------------------------------- | |
89 | class AliTRDseed : public AliTRDtrack { | |
90 | Float_t fdEdx[200]; | |
91 | public: | |
92 | AliTRDseed():AliTRDtrack(){} | |
93 | AliTRDseed(UInt_t index, const Double_t xx[5], | |
94 | const Double_t cc[15], Double_t xr, Double_t alpha): | |
95 | AliTRDtrack(index, xx, cc, xr, alpha) {} | |
96 | void SetSampledEdx(Float_t q, Int_t i) { | |
97 | Double_t c=GetC(), e=GetEta(), t=GetTgl(), x=GetX(); | |
98 | q *= TMath::Sqrt((1-(c*x-e)*(c*x-e))/(1+t*t)); | |
99 | fdEdx[i]=q; | |
100 | } | |
101 | void CookdEdx(Double_t low=0.05, Double_t up=0.70); | |
102 | }; | |
103 | ||
104 | #endif |