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 <AliKalmanTrack.h>
   8 #include <TMath.h>
   9
  10 class AliTRDcluster;
  11 class AliTPCtrack;
  12
  13 const unsigned kMAX_CLUSTERS_PER_TRACK=210;
  14
  15 class AliTRDtrack : public AliKalmanTrack {
  16
  17 // Represents reconstructed TRD track
  18
  19 public:
  20
  21    AliTRDtrack():AliKalmanTrack(){}
  22    AliTRDtrack(const AliTRDcluster *c, UInt_t index, const Double_t xx[5],
  23                const Double_t cc[15], Double_t xr, Double_t alpha);
  24    AliTRDtrack(const AliTRDtrack& t);
  25    AliTRDtrack(const AliKalmanTrack& t, Double_t alpha);
  26
  27    Int_t    Compare(const TObject *o) const;
  28    void     CookdEdx(Double_t low=0.05, Double_t up=0.70);
  29
  30    Double_t GetAlpha() const {return fAlpha;}
  31    Double_t GetC()     const {return fC;}
  32    Int_t    GetClusterIndex(Int_t i) const {return fIndex[i];}
  33    Float_t  GetClusterdQdl(Int_t i) const {return fdQdl[i];}
  34
  35    void     GetCovariance(Double_t cov[15]) const;
  36    Float_t  GetdEdx()  const {return fdEdx;}
  37    Double_t GetEta()   const {return fE;}
  38
  39    void     GetExternalCovariance(Double_t cov[15]) const ;
  40    void     GetExternalParameters(Double_t& xr, Double_t x[5]) const ;
  41
  42    Double_t GetLikelihoodElectron() const { return fLhElectron; };
  43
  44    Double_t Get1Pt()   const {return TMath::Abs(fC*GetConvConst());}
  45    Double_t GetP()     const {
  46      return TMath::Abs(GetPt())*sqrt(1.+GetTgl()*GetTgl());
  47    }
  48    Double_t GetPredictedChi2(const AliTRDcluster*) const ;
  49    Double_t GetPt()    const {return 1./Get1Pt();}
  50    void     GetPxPyPz(Double_t &px, Double_t &py, Double_t &pz) const ;
  51    void     GetGlobalXYZ(Double_t &x, Double_t &y, Double_t &z) const ;
  52    Int_t    GetSeedLabel() const { return fSeedLab; }
  53    Double_t GetSigmaC2()   const {return fCcc;}
  54    Double_t GetSigmaTgl2() const {return fCtt;}
  55    Double_t GetSigmaY2()   const {return fCyy;}
  56    Double_t GetSigmaZ2()   const {return fCzz;}
  57    Double_t GetSnp()  const {return fX*fC - fE;}
  58    Double_t GetTgl()  const {return fT;}
  59    Double_t GetX()    const {return fX;}
  60    Double_t GetY()    const {return fY;}
  61    Double_t GetZ()    const {return fZ;}
  62
  63    Int_t    PropagateTo(Double_t xr,
  64                    Double_t x0=8.72,Double_t rho=5.86e-3,Double_t pm=0.139);
  65    void     ResetCovariance();
  66    Int_t    Rotate(Double_t angle);
  67
  68    void     SetdEdx(Float_t dedx) {fdEdx=dedx;}
  69    void     SetLikelihoodElectron(Float_t l) { fLhElectron = l; };
  70
  71    void     SetSampledEdx(Float_t q, Int_t i) {
  72               Double_t s=GetSnp(), t=GetTgl();
  73               q*= TMath::Sqrt((1-s*s)/(1+t*t));
  74               fdQdl[i]=q;
  75             }
  76
  77    void     SetSeedLabel(Int_t lab) { fSeedLab=lab; }
  78
  79    Int_t    Update(const AliTRDcluster* c, Double_t chi2, UInt_t i);
  80
  81
  82
  83 protected:
  84
  85    Int_t    fSeedLab;     // track label taken from seeding
  86    Float_t  fdEdx;        // dE/dx
  87
  88    Double_t fAlpha;       // rotation angle
  89    Double_t fX;           // running local X-coordinate of the track (time bin)
  90
  91    Double_t fY;           // Y-coordinate of the track
  92    Double_t fZ;           // Z-coordinate of the track
  93    Double_t fC;           // track curvature
  94    Double_t fE;           // C*x0
  95    Double_t fT;           // tangent of the track dip angle
  96
  97    Double_t fCyy;                         // covariance
  98    Double_t fCzy, fCzz;                   // matrix
  99    Double_t fCcy, fCcz, fCcc;             // of the
 100    Double_t fCey, fCez, fCec, fCee;       // track
 101    Double_t fCty, fCtz, fCtc, fCte, fCtt; // parameters
 102
 103    UInt_t  fIndex[kMAX_CLUSTERS_PER_TRACK];  // global indexes of clusters
 104    Float_t fdQdl[kMAX_CLUSTERS_PER_TRACK];   // cluster amplitudes corrected
 105                                              // for track angles
 106
 107    Float_t fLhElectron;    // Likelihood to be an electron
 108
 109    ClassDef(AliTRDtrack,2) // TRD reconstructed tracks
 110
 111 };
 112
 113
 114 #endif