two errors that the HP-UX compiler could not correct
[u/mrichter/AliRoot.git] / TRD / AliTRDtrack.h
CommitLineData
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
9class AliTRDcluster;
10
11class 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
38public:
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//-----------------------------------------------------------------
89class AliTRDseed : public AliTRDtrack {
90 Float_t fdEdx[200];
91public:
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