4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 #include <AliKalmanTrack.h>
13 const unsigned kMAX_CLUSTERS_PER_TRACK=210;
15 class AliTRDtrack : public AliKalmanTrack {
17 // Represents reconstructed TRD track
21 AliTRDtrack():AliKalmanTrack(){}
22 AliTRDtrack(const AliTRDcluster *c, UInt_t index, const Double_t xx[5],
23 const Double_t cc[15], Double_t xr, Double_t alpha);
24 AliTRDtrack(const AliTRDtrack& t);
25 AliTRDtrack(const AliKalmanTrack& t, Double_t alpha);
27 Int_t Compare(const TObject *o) const;
28 void CookdEdx(Double_t low=0.05, Double_t up=0.70);
30 Double_t GetAlpha() const {return fAlpha;}
31 Double_t GetC() const {return fC;}
32 Int_t GetClusterIndex(Int_t i) const {return fIndex[i];}
33 Float_t GetClusterdQdl(Int_t i) const {return fdQdl[i];}
35 void GetCovariance(Double_t cov[15]) const;
36 Float_t GetdEdx() const {return fdEdx;}
37 Double_t GetEta() const {return fE;}
39 void GetExternalCovariance(Double_t cov[15]) const ;
40 void GetExternalParameters(Double_t& xr, Double_t x[5]) const ;
42 Double_t GetLikelihoodElectron() const { return fLhElectron; };
44 Double_t Get1Pt() const {return TMath::Abs(fC*GetConvConst());}
45 Double_t GetP() const {
46 return TMath::Abs(GetPt())*sqrt(1.+GetTgl()*GetTgl());
48 Double_t GetPredictedChi2(const AliTRDcluster*) const ;
49 Double_t GetPt() const {return 1./Get1Pt();}
50 void GetPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const ;
51 void GetGlobalXYZ(Double_t &x, Double_t &y, Double_t &z) const ;
52 Int_t GetSeedLabel() const { return fSeedLab; }
53 Double_t GetSigmaC2() const {return fCcc;}
54 Double_t GetSigmaTgl2() const {return fCtt;}
55 Double_t GetSigmaY2() const {return fCyy;}
56 Double_t GetSigmaZ2() const {return fCzz;}
57 Double_t GetSnp() const {return fX*fC - fE;}
58 Double_t GetTgl() const {return fT;}
59 Double_t GetX() const {return fX;}
60 Double_t GetY() const {return fY;}
61 Double_t GetZ() const {return fZ;}
63 Int_t PropagateTo(Double_t xr,
64 Double_t x0=8.72,Double_t rho=5.86e-3,Double_t pm=0.139);
65 void ResetCovariance();
66 Int_t Rotate(Double_t angle);
68 void SetdEdx(Float_t dedx) {fdEdx=dedx;}
69 void SetLikelihoodElectron(Float_t l) { fLhElectron = l; };
71 void SetSampledEdx(Float_t q, Int_t i) {
72 Double_t s=GetSnp(), t=GetTgl();
73 q*= TMath::Sqrt((1-s*s)/(1+t*t));
77 void SetSeedLabel(Int_t lab) { fSeedLab=lab; }
79 Int_t Update(const AliTRDcluster* c, Double_t chi2, UInt_t i);
85 Int_t fSeedLab; // track label taken from seeding
86 Float_t fdEdx; // dE/dx
88 Double_t fAlpha; // rotation angle
89 Double_t fX; // running local X-coordinate of the track (time bin)
91 Double_t fY; // Y-coordinate of the track
92 Double_t fZ; // Z-coordinate of the track
93 Double_t fC; // track curvature
95 Double_t fT; // tangent of the track dip angle
97 Double_t fCyy; // covariance
98 Double_t fCzy, fCzz; // matrix
99 Double_t fCcy, fCcz, fCcc; // of the
100 Double_t fCey, fCez, fCec, fCee; // track
101 Double_t fCty, fCtz, fCtc, fCte, fCtt; // parameters
103 UInt_t fIndex[kMAX_CLUSTERS_PER_TRACK]; // global indexes of clusters
104 Float_t fdQdl[kMAX_CLUSTERS_PER_TRACK]; // cluster amplitudes corrected
107 Float_t fLhElectron; // Likelihood to be an electron
109 ClassDef(AliTRDtrack,2) // TRD reconstructed tracks