Update of tracking code for tilted pads
[u/mrichter/AliRoot.git] / TRD / AliTRDtrack.h
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 class AliTRDcluster;
11 class AliTPCtrack; 
12
13 const unsigned kMAX_CLUSTERS_PER_TRACK=210; 
14
15 class AliTRDtrack : public AliKalmanTrack {
16
17 // Represents reconstructed TRD track
18
19 public:
20
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); 
26
27    Int_t    Compare(const TObject *o) const;
28    void     CookdEdx(Double_t low=0.05, Double_t up=0.70);   
29
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];}    
34
35    void     GetCovariance(Double_t cov[15]) const;  
36    Float_t  GetdEdx()  const {return fdEdx;}
37    Double_t GetEta()   const {return fE;}
38
39    void     GetExternalCovariance(Double_t cov[15]) const ;   
40    void     GetExternalParameters(Double_t& xr, Double_t x[5]) const ;
41
42    Double_t GetLikelihoodElectron() const { return fLhElectron; };
43
44    Double_t Get1Pt()   const {return TMath::Abs(fC*GetConvConst());} 
45    Double_t GetP()     const {  
46      return TMath::Abs(GetPt())*sqrt(1.+GetTgl()*GetTgl());
47    }
48    Double_t GetPredictedChi2(const AliTRDcluster*, Double_t h01) 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;}
62
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);
67
68    void     SetdEdx(Float_t dedx) {fdEdx=dedx;}  
69    void     SetLikelihoodElectron(Float_t l) { fLhElectron = l; };  
70
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));
74               fdQdl[i]=q;
75             }     
76
77    void     SetSeedLabel(Int_t lab) { fSeedLab=lab; }
78
79    Int_t    Update(const AliTRDcluster* c, Double_t chi2, UInt_t i, 
80                    Double_t h01);
81
82
83
84 protected:
85
86    Int_t    fSeedLab;     // track label taken from seeding  
87    Float_t  fdEdx;        // dE/dx 
88
89    Double_t fAlpha;       // rotation angle
90    Double_t fX;           // running local X-coordinate of the track (time bin)
91
92    Double_t fY;           // Y-coordinate of the track
93    Double_t fZ;           // Z-coordinate of the track
94    Double_t fC;           // track curvature
95    Double_t fE;           // C*x0
96    Double_t fT;           // tangent of the track dip angle   
97
98    Double_t fCyy;                         // covariance
99    Double_t fCzy, fCzz;                   // matrix
100    Double_t fCcy, fCcz, fCcc;             // of the
101    Double_t fCey, fCez, fCec, fCee;       // track
102    Double_t fCty, fCtz, fCtc, fCte, fCtt; // parameters   
103
104    UInt_t  fIndex[kMAX_CLUSTERS_PER_TRACK];  // global indexes of clusters  
105    Float_t fdQdl[kMAX_CLUSTERS_PER_TRACK];   // cluster amplitudes corrected 
106                                              // for track angles    
107                            
108    Float_t fLhElectron;    // Likelihood to be an electron    
109
110    ClassDef(AliTRDtrack,2) // TRD reconstructed tracks
111
112 };                     
113
114
115 #endif