4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
11 class AliTRDtrack : public TObject {
13 // Represents reconstructed TRD track
15 Int_t fLab; // track label
16 Double_t fChi2; // total chi2 value for the track
17 Float_t fdEdx; // dE/dx
19 Double_t fAlpha; // rotation angle
20 Double_t fX; // running local X-coordinate of the track (time bin)
22 Double_t fY; // Y-coordinate of the track
23 Double_t fZ; // Z-coordinate of the track
24 Double_t fC; // track curvature
26 Double_t fT; // tangent of the track dip angle
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
34 Short_t fN; // number of clusters associated with the track
35 UInt_t fIndex[200]; // global indexes of these clusters
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);
45 Int_t Compare(TObject *o);
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());
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;}
71 Bool_t IsSortable() const {return kTRUE;}
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);
79 void SetLabel(Int_t lab=0) {fLab=lab;}
80 void SetdEdx(Float_t dedx) {fdEdx=dedx;}
82 void Update(const AliTRDcluster* c, Double_t chi2, UInt_t i);
84 ClassDef(AliTRDtrack,1) // TRD reconstructed tracks
88 //-----------------------------------------------------------------
89 class AliTRDseed : public AliTRDtrack {
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));
101 void CookdEdx(Double_t low=0.05, Double_t up=0.70);