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