TRD/AliTRDtrack.h

   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
   9 class AliTRDcluster;
  10
  11 class 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
  38 public:
  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 //-----------------------------------------------------------------
  89 class AliTRDseed : public AliTRDtrack {
  90    Float_t fdEdx[200];
  91 public:
  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