HLT/src/AliL3Track.h

   1 // @(#) $Id$
   2
   3 #ifndef ALIL3TRACK_H
   4 #define ALIL3TRACK_H
   5
   6 #include <string.h>
   7
   8 #include "AliL3RootTypes.h"
   9
  10 class AliL3Vertex;
  11
  12 class AliL3Track {
  13
  14  private:
  15
  16   Int_t fNHits;
  17
  18   Int_t fMCid;  //Assigned id from MC data.
  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)
  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:
  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)
  37
  38   Double_t fPterr;
  39   Double_t fPsierr;
  40   Double_t fZ0err;
  41   Double_t fTanlerr;
  42
  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
  47   Double_t fFirstPoint[3];
  48   Double_t fLastPoint[3];
  49   Double_t fPoint[3];
  50   Double_t fPointPsi; //azimuthal angle of the momentum at Point
  51
  52
  53   Bool_t fIsPoint;    //Helix crosses the X-plane
  54   Bool_t IsPoint(Bool_t ispoint) {fIsPoint = ispoint;return fIsPoint;}
  55
  56   Bool_t fIsLocal; //Track given in local coordinates.
  57
  58   UInt_t fHitNumbers[159];  //Array of hit numbers for this track
  59
  60  public:
  61
  62   AliL3Track();
  63   virtual ~AliL3Track();
  64
  65   virtual void Set(AliL3Track* track);
  66   virtual Int_t Compare(const AliL3Track *track) const;
  67   virtual void CalculateHelix();
  68
  69   Bool_t CalculateReferencePoint(Double_t angle,Double_t radius=132);//Calculate Reference Point
  70   Bool_t CalculateEdgePoint(Double_t angle);//Calculate crossing point with line
  71   Bool_t CalculatePoint(Double_t xplane);   //Calculate crossing point with X-plane
  72   Bool_t IsPoint() {return fIsPoint;}
  73   Double_t GetCrossingAngle(Int_t padrow,Int_t slice=-1);
  74   Bool_t GetCrossingPoint(Int_t padrow,Float_t *xyz);
  75   Double_t GetDistance(Double_t x0,Double_t x1){return 0;}
  76   void UpdateToFirstPoint();
  77
  78   void GetClosestPoint(AliL3Vertex *vertex,Double_t &closest_x,Double_t &closest_y,Double_t &closest_z);
  79   void Rotate(Int_t slice,Bool_t tolocal=kFALSE);
  80   Bool_t IsLocal() {return fIsLocal;}
  81
  82   // getter
  83   Double_t GetFirstPointX() {return fFirstPoint[0];}
  84   Double_t GetFirstPointY() {return fFirstPoint[1];}
  85   Double_t GetFirstPointZ() {return fFirstPoint[2];}
  86   Double_t GetLastPointX() {return fLastPoint[0];}
  87   Double_t GetLastPointY() {return fLastPoint[1];}
  88   Double_t GetLastPointZ() {return fLastPoint[2];}
  89
  90   Double_t GetPointPsi() {return fPointPsi;}
  91   Double_t GetPointX() {return fPoint[0];}
  92   Double_t GetPointY() {return fPoint[1];}
  93   Double_t GetPointZ() {return fPoint[2];}
  94
  95   Double_t GetPt() const {return fPt;}
  96   Double_t GetTgl() const {return fTanl;}
  97   Double_t GetPsi() const {return fPsi;}
  98   Double_t GetPhi0() const {return fPhi0;}
  99   Double_t GetR0() const {return fR0;}
 100   Double_t GetZ0() const {return fFirstPoint[2];}
 101   //Double_t GetZ0() const {return fZ0;}
 102
 103   Double_t GetKappa() const {return fKappa;}
 104   Double_t GetRadius() const {return fRadius;}
 105   Double_t GetCenterX() const {return fCenterX;}
 106   Double_t GetCenterY() const {return fCenterY;}
 107
 108   Int_t GetNHits() const {return fNHits;}
 109   Int_t   GetNumberOfPoints()   const {return fNHits;}
 110   Bool_t  ComesFromMainVertex() const {return fFromMainVertex;}
 111
 112   Double_t GetPx() const {return fPt*cos(fPsi);}
 113   Double_t GetPy() const {return fPt*sin(fPsi);}
 114   Double_t GetPz() const {return fPt*fTanl;}
 115
 116   Double_t GetP() const;
 117   Double_t GetPseudoRapidity() const;
 118   //Double_t   GetEta() const;
 119   Double_t   GetRapidity() const;
 120
 121   Int_t GetCharge() const {return fQ;}
 122   Int_t GetMCid() const {return fMCid;}
 123   Double_t GetLength() const {return fLength;}
 124
 125   Int_t GetFirstRow() const {return fRowRange[0];}
 126   Int_t GetLastRow()  const {return fRowRange[1];}
 127
 128   UInt_t *GetHitNumbers() {return fHitNumbers;}
 129
 130   // setter
 131
 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;}
 142   //void SetZ0(Double_t f) {fZ0 = f;}
 143   void SetZ0(Double_t f) {fFirstPoint[2] = f;}
 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
 154   void SetRowRange(Int_t f,Int_t g) {fRowRange[0]=f; fRowRange[1]=g;}
 155   void SetSector(Int_t f) {fSector = f;}
 156
 157   void SetRadius(Double_t f) {fRadius = f;}
 158   void SetCenterX(Double_t f) {fCenterX = f;}
 159   void SetCenterY(Double_t f) {fCenterY = f;}
 160
 161   void SetCharge(Int_t f) {fQ = f;}
 162
 163   void ComesFromMainVertex(Bool_t f) {fFromMainVertex = f;}
 164
 165   ClassDef(AliL3Track,1) //Base track class
 166 };
 167
 168 #endif
 169