Further modifications in OpenOutput and WriteCluster
[u/mrichter/AliRoot.git] / TRD / AliTRDtrack.h
CommitLineData
46d29e70 1#ifndef ALITRDTRACK_H
b3a5a838 2#define ALITRDTRACK_H
46d29e70 3
4/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
b3a5a838 5 * See cxx source for full Copyright notice */
46d29e70 6
5443e65e 7#include <AliKalmanTrack.h>
b3a5a838 8#include <TMath.h>
46d29e70 9
10class AliTRDcluster;
b3a5a838 11class AliTPCtrack;
5443e65e 12
13const unsigned kMAX_CLUSTERS_PER_TRACK=210;
46d29e70 14
5443e65e 15class AliTRDtrack : public AliKalmanTrack {
46d29e70 16
17// Represents reconstructed TRD track
18
46d29e70 19public:
20
5443e65e 21 AliTRDtrack():AliKalmanTrack(){}
a819a5f7 22 AliTRDtrack(const AliTRDcluster *c, UInt_t index, const Double_t xx[5],
46d29e70 23 const Double_t cc[15], Double_t xr, Double_t alpha);
24 AliTRDtrack(const AliTRDtrack& t);
5443e65e 25 AliTRDtrack(const AliKalmanTrack& t, Double_t alpha);
46d29e70 26
a819a5f7 27 Int_t Compare(const TObject *o) const;
28 void CookdEdx(Double_t low=0.05, Double_t up=0.70);
46d29e70 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];}
a819a5f7 33 Float_t GetClusterdQdl(Int_t i) const {return fdQdl[i];}
5443e65e 34
46d29e70 35 void GetCovariance(Double_t cov[15]) const;
5443e65e 36 Float_t GetdEdx() const {return fdEdx;}
46d29e70 37 Double_t GetEta() const {return fE;}
5443e65e 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());}
46d29e70 45 Double_t GetP() const {
46 return TMath::Abs(GetPt())*sqrt(1.+GetTgl()*GetTgl());
47 }
fd621f36 48 Double_t GetPredictedChi2(const AliTRDcluster*, Double_t h01) const ;
5443e65e 49 Double_t GetPt() const {return 1./Get1Pt();}
46d29e70 50 void GetPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const ;
5443e65e 51 void GetGlobalXYZ(Double_t &x, Double_t &y, Double_t &z) const ;
52 Int_t GetSeedLabel() const { return fSeedLab; }
46d29e70 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;}
5443e65e 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;}
46d29e70 62
b3a5a838 63 Int_t PropagateTo(Double_t xr, Double_t x0=8.72, Double_t rho=5.86e-3);
5443e65e 64 void ResetCovariance();
46d29e70 65 Int_t Rotate(Double_t angle);
66
46d29e70 67 void SetdEdx(Float_t dedx) {fdEdx=dedx;}
5443e65e 68 void SetLikelihoodElectron(Float_t l) { fLhElectron = l; };
69
70 void SetSampledEdx(Float_t q, Int_t i) {
71 Double_t s=GetSnp(), t=GetTgl();
72 q*= TMath::Sqrt((1-s*s)/(1+t*t));
73 fdQdl[i]=q;
74 }
75
76 void SetSeedLabel(Int_t lab) { fSeedLab=lab; }
77
fd621f36 78 Int_t Update(const AliTRDcluster* c, Double_t chi2, UInt_t i,
79 Double_t h01);
80
81
46d29e70 82
bbf92647 83protected:
84
5443e65e 85 Int_t fSeedLab; // track label taken from seeding
bbf92647 86 Float_t fdEdx; // dE/dx
87
88 Double_t fAlpha; // rotation angle
89 Double_t fX; // running local X-coordinate of the track (time bin)
90
91 Double_t fY; // Y-coordinate of the track
92 Double_t fZ; // Z-coordinate of the track
bbf92647 93 Double_t fE; // C*x0
94 Double_t fT; // tangent of the track dip angle
b3a5a838 95 Double_t fC; // track curvature
bbf92647 96
97 Double_t fCyy; // covariance
98 Double_t fCzy, fCzz; // matrix
b3a5a838 99 Double_t fCey, fCez, fCee; // of the
100 Double_t fCty, fCtz, fCte, fCtt; // track
101 Double_t fCcy, fCcz, fCce, fCct, fCcc; // parameters
bbf92647 102
5443e65e 103 UInt_t fIndex[kMAX_CLUSTERS_PER_TRACK]; // global indexes of clusters
104 Float_t fdQdl[kMAX_CLUSTERS_PER_TRACK]; // cluster amplitudes corrected
105 // for track angles
a819a5f7 106
107 Float_t fLhElectron; // Likelihood to be an electron
bbf92647 108
a819a5f7 109 ClassDef(AliTRDtrack,2) // TRD reconstructed tracks
a2b90f83 110
46d29e70 111};
112
113
46d29e70 114#endif