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