]>
Commit | Line | Data |
---|---|---|
3e87ef69 | 1 | // @(#) $Id$ |
2 | ||
108615fc | 3 | #ifndef ALIL3TRACK_H |
4 | #define ALIL3TRACK_H | |
5 | ||
6d4f1901 | 6 | class AliL3Vertex; |
de3c3890 | 7 | class AliL3SpacePointData; |
6d4f1901 | 8 | |
108615fc | 9 | class AliL3Track { |
10 | ||
11 | private: | |
12 | ||
b1ed0288 | 13 | Int_t fNHits; //Number of hits |
108615fc | 14 | Int_t fMCid; //Assigned id from MC data. |
eeddc64d | 15 | |
16 | Double_t fKappa; // Signed curvature (projected to a circle) | |
17 | Double_t fRadius; // Radius of the helix (projected to a circle) | |
18 | Double_t fCenterX; // x coordinate of the center of the helix (projected to a circle) | |
19 | Double_t fCenterY; // y coordinate of the center of the helix (projected to a circle) | |
108615fc | 20 | Bool_t fFromMainVertex; // true if tracks origin is the main vertex, otherwise false |
21 | ||
22 | Int_t fRowRange[2]; //Subsector where this track was build | |
23 | Int_t fSector; //Sector # where this track was build | |
24 | ||
25 | //data from momentum fit | |
26 | Int_t fQ; //charge measured fit | |
27 | ||
28 | //track parameters: | |
eeddc64d | 29 | Double_t fTanl; //tan of dipangle |
30 | Double_t fPsi; //azimuthal angle of the momentum | |
31 | Double_t fPt; //transverse momentum | |
32 | Double_t fLength; //length of track (s) | |
108615fc | 33 | |
b1ed0288 | 34 | Double_t fPterr; //errors for fast hough |
35 | Double_t fPsierr; //errors for fast hough | |
36 | Double_t fZ0err; //errors for fast hough | |
37 | Double_t fTanlerr; //errors for fast hough | |
108615fc | 38 | |
eeddc64d | 39 | Double_t fPhi0; //azimuthal angle of the first point |
40 | Double_t fR0; //radius of the first point | |
41 | Double_t fZ0; //z coordinate of the first point (fFirstPoint[2]) | |
42 | ||
b1ed0288 | 43 | Double_t fFirstPoint[3]; //first point |
44 | Double_t fLastPoint[3]; //last point | |
45 | Double_t fPoint[3]; //point | |
108615fc | 46 | Double_t fPointPsi; //azimuthal angle of the momentum at Point |
eeddc64d | 47 | |
b1ed0288 | 48 | Bool_t fIsPoint; //Helix crosses the X-plane |
108615fc | 49 | Bool_t fIsLocal; //Track given in local coordinates. |
50 | ||
b1ed0288 | 51 | Float_t fPID; //pid |
52 | UInt_t fHitNumbers[159]; //Array of hit numbers for this track | |
53 | ||
54 | Bool_t IsPoint(Bool_t ispoint) {fIsPoint = ispoint;return fIsPoint;} | |
108615fc | 55 | |
108615fc | 56 | public: |
57 | ||
58 | AliL3Track(); | |
59 | virtual ~AliL3Track(); | |
60 | ||
61 | virtual void Set(AliL3Track* track); | |
203925a9 | 62 | virtual Int_t Compare(const AliL3Track *track) const; |
6d24ee88 | 63 | virtual void CalculateHelix(); |
108615fc | 64 | |
6d4f1901 | 65 | Bool_t CalculateReferencePoint(Double_t angle,Double_t radius=132);//Calculate Reference Point |
108615fc | 66 | Bool_t CalculateEdgePoint(Double_t angle);//Calculate crossing point with line |
eeddc64d | 67 | Bool_t CalculatePoint(Double_t xplane); //Calculate crossing point with X-plane |
108615fc | 68 | Bool_t IsPoint() {return fIsPoint;} |
3e87ef69 | 69 | Double_t GetCrossingAngle(Int_t padrow,Int_t slice=-1); |
0391971c | 70 | Bool_t GetCrossingPoint(Int_t padrow,Float_t *xyz); |
0615a438 | 71 | Double_t GetDistance(Double_t /*x0*/,Double_t /*x1*/){return 0;} |
3e87ef69 | 72 | void UpdateToFirstPoint(); |
73 | ||
b1ed0288 | 74 | Float_t GetPID() const {return fPID;} |
de3c3890 | 75 | void SetPID(Float_t pid) {fPID=pid;} |
76 | ||
6d4f1901 | 77 | void GetClosestPoint(AliL3Vertex *vertex,Double_t &closest_x,Double_t &closest_y,Double_t &closest_z); |
203925a9 | 78 | void Rotate(Int_t slice,Bool_t tolocal=kFALSE); |
b1ed0288 | 79 | Bool_t IsLocal() const {return fIsLocal;} |
3e87ef69 | 80 | |
108615fc | 81 | // getter |
b1ed0288 | 82 | Double_t GetFirstPointX() const {return fFirstPoint[0];} |
83 | Double_t GetFirstPointY() const {return fFirstPoint[1];} | |
84 | Double_t GetFirstPointZ() const {return fFirstPoint[2];} | |
85 | Double_t GetLastPointX() const {return fLastPoint[0];} | |
86 | Double_t GetLastPointY() const {return fLastPoint[1];} | |
87 | Double_t GetLastPointZ() const {return fLastPoint[2];} | |
88 | ||
89 | Double_t GetPointPsi() const {return fPointPsi;} | |
90 | Double_t GetPointX() const {return fPoint[0];} | |
91 | Double_t GetPointY() const {return fPoint[1];} | |
92 | Double_t GetPointZ() const {return fPoint[2];} | |
108615fc | 93 | |
94 | Double_t GetPt() const {return fPt;} | |
95 | Double_t GetTgl() const {return fTanl;} | |
108615fc | 96 | Double_t GetPsi() const {return fPsi;} |
eeddc64d | 97 | Double_t GetPhi0() const {return fPhi0;} |
108615fc | 98 | Double_t GetR0() const {return fR0;} |
2d04dcbf | 99 | Double_t GetZ0() const {return fFirstPoint[2];} |
b1ed0288 | 100 | //const Double_t GetZ0() const {return fZ0;} |
2d04dcbf | 101 | |
eeddc64d | 102 | Double_t GetKappa() const {return fKappa;} |
103 | Double_t GetRadius() const {return fRadius;} | |
104 | Double_t GetCenterX() const {return fCenterX;} | |
105 | Double_t GetCenterY() const {return fCenterY;} | |
108615fc | 106 | |
203925a9 | 107 | Int_t GetNHits() const {return fNHits;} |
108615fc | 108 | Int_t GetNumberOfPoints() const {return fNHits;} |
eeddc64d | 109 | Bool_t ComesFromMainVertex() const {return fFromMainVertex;} |
108615fc | 110 | |
eeddc64d | 111 | Double_t GetPx() const {return fPt*cos(fPsi);} |
112 | Double_t GetPy() const {return fPt*sin(fPsi);} | |
113 | Double_t GetPz() const {return fPt*fTanl;} | |
108615fc | 114 | |
eeddc64d | 115 | Double_t GetP() const; |
116 | Double_t GetPseudoRapidity() const; | |
203925a9 | 117 | //Double_t GetEta() const; |
108615fc | 118 | Double_t GetRapidity() const; |
119 | ||
eeddc64d | 120 | Int_t GetCharge() const {return fQ;} |
108615fc | 121 | Int_t GetMCid() const {return fMCid;} |
eeddc64d | 122 | Double_t GetLength() const {return fLength;} |
108615fc | 123 | |
eeddc64d | 124 | Int_t GetFirstRow() const {return fRowRange[0];} |
108615fc | 125 | Int_t GetLastRow() const {return fRowRange[1];} |
de3c3890 | 126 | Int_t GetSector() const {return fSector;} |
108615fc | 127 | |
128 | UInt_t *GetHitNumbers() {return fHitNumbers;} | |
129 | ||
130 | // setter | |
203925a9 | 131 | |
108615fc | 132 | void SetMCid(Int_t f) {fMCid = f;} |
133 | void SetFirstPoint(Double_t f,Double_t g,Double_t h) {fFirstPoint[0]=f; fFirstPoint[1]=g; fFirstPoint[2]=h;} | |
134 | void SetLastPoint(Double_t f,Double_t g,Double_t h) {fLastPoint[0]=f; fLastPoint[1]=g; fLastPoint[2]=h;} | |
135 | ||
136 | void SetHits(Int_t nhits,UInt_t *hits) {memcpy(fHitNumbers,hits,nhits*sizeof(UInt_t));} | |
137 | ||
138 | void SetPhi0(Double_t f) {fPhi0 = f;} | |
139 | void SetPsi(Double_t f) {fPsi = f;} | |
140 | void SetR0(Double_t f) {fR0 = f;} | |
141 | void SetTgl(Double_t f) {fTanl =f;} | |
2d04dcbf | 142 | //void SetZ0(Double_t f) {fZ0 = f;} |
143 | void SetZ0(Double_t f) {fFirstPoint[2] = f;} | |
108615fc | 144 | void SetPt(Double_t f) {fPt = f;} |
145 | void SetLength(Double_t f) {fLength = f;} | |
146 | void SetPterr(Double_t f) {fPterr = f;} | |
147 | void SetPsierr(Double_t f) {fPsierr = f;} | |
148 | void SetZ0err(Double_t f) {fZ0err = f;} | |
149 | void SetTglerr(Double_t f) {fTanlerr = f;} | |
150 | void SetKappa(Double_t f) {fKappa = f;} | |
151 | ||
152 | void SetNHits(Int_t f) {fNHits = f;} | |
153 | ||
eeddc64d | 154 | void SetRowRange(Int_t f,Int_t g) {fRowRange[0]=f; fRowRange[1]=g;} |
155 | void SetSector(Int_t f) {fSector = f;} | |
108615fc | 156 | |
eeddc64d | 157 | void SetRadius(Double_t f) {fRadius = f;} |
158 | void SetCenterX(Double_t f) {fCenterX = f;} | |
159 | void SetCenterY(Double_t f) {fCenterY = f;} | |
108615fc | 160 | |
eeddc64d | 161 | void SetCharge(Int_t f) {fQ = f;} |
108615fc | 162 | |
eeddc64d | 163 | void ComesFromMainVertex(Bool_t f) {fFromMainVertex = f;} |
108615fc | 164 | |
2d04dcbf | 165 | ClassDef(AliL3Track,1) //Base track class |
108615fc | 166 | }; |
167 | ||
168 | #endif |