STEER/AliKalmanTrack.h

   1 #ifndef ALIKALMANTRACK_H
   2 #define ALIKALMANTRACK_H
   3
   4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   5  * See cxx source for full Copyright notice                               */
   6
   7 /* $Id$ */
   8
   9 //-------------------------------------------------------------------------
  10 //                          Class AliKalmanTrack
  11 //      fixed the interface for the derived reconstructed track classes
  12 //            Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
  13 //-------------------------------------------------------------------------
  14
  15 #include <TObject.h>
  16 #include "AliLog.h"
  17 #include "AliPID.h"
  18
  19 class AliCluster;
  20
  21 class AliKalmanTrack : public TObject {
  22 public:
  23   AliKalmanTrack();
  24   AliKalmanTrack(const AliKalmanTrack &t);
  25
  26   virtual ~AliKalmanTrack(){};
  27   void SetLabel(Int_t lab) {fLab=lab;}
  28   void SetFakeRatio(Float_t ratio) {fFakeRatio=ratio;}
  29
  30   Bool_t   IsSortable() const {return kTRUE;}
  31   Int_t    GetLabel()   const {return fLab;}
  32   Float_t    GetFakeRatio()   const {return fFakeRatio;}
  33   Double_t GetChi2()    const {return fChi2;}
  34   Double_t GetMass()    const {return fMass;}
  35   Int_t    GetNumberOfClusters() const {return fN;}
  36   virtual Int_t GetClusterIndex(Int_t) const { //reserved for AliTracker
  37     AliWarning("Method must be overloaded !\n");
  38     return 0;
  39   }
  40   virtual Double_t GetPIDsignal() const {
  41     AliWarning("Method must be overloaded !\n");
  42     return 0.;
  43   }
  44
  45   virtual Double_t GetDCA(const AliKalmanTrack *p,Double_t &xthis,Double_t &xp) const;
  46   virtual
  47   Double_t PropagateToDCA(AliKalmanTrack *p, Double_t d=0., Double_t x0=0.);
  48   virtual Double_t GetAlpha() const {
  49     AliWarning("Method must be overloaded !\n");
  50     return 0.;
  51   }
  52   virtual Double_t GetSigmaY2() const {
  53     AliWarning("Method must be overloaded !\n");
  54     return 0.;
  55   }
  56   virtual Double_t GetSigmaZ2() const {
  57     AliWarning("Method must be overloaded !\n");
  58     return 0.;
  59   }
  60
  61   virtual Int_t Compare(const TObject *) const {return 0;}
  62
  63   virtual void GetExternalParameters(Double_t&/*xr*/,Double_t/*x*/[5]) const=0;
  64   virtual void GetExternalCovariance(Double_t /*cov*/[15]) const = 0;
  65
  66   virtual Double_t GetPredictedChi2(const AliCluster *) const = 0;
  67   virtual Int_t PropagateTo(Double_t/*xr*/,Double_t/*x0*/,Double_t/*rho*/) = 0;
  68   //virtual Int_t PropagateToVertex(Double_t /*d*/=0., Double_t /*x0*/=0.) = 0;
  69   virtual Int_t Update(const AliCluster*, Double_t /*chi2*/, UInt_t) = 0;
  70
  71   //static Double_t GetConvConst();
  72   static Double_t MeanMaterialBudget(Double_t *start, Double_t *end, Double_t *mparam);
  73
  74   // Time integration (S.Radomski@gsi.de)
  75   void   StartTimeIntegral();
  76   void SetIntegratedLength(Double_t l) {fIntegratedLength=l;}
  77   void SetIntegratedTimes(const Double_t *times);
  78
  79   Bool_t IsStartedTimeIntegral() const {return fStartTimeIntegral;}
  80   void     AddTimeStep(Double_t length);
  81   void GetIntegratedTimes(Double_t *times) const;
  82   Double_t GetIntegratedTime(Int_t pdg) const;
  83   Double_t GetIntegratedLength() const {return fIntegratedLength;}
  84   virtual void GetXYZ(Float_t r[3]) const = 0;
  85
  86   static Double_t GetConvConst();
  87   void SetNumberOfClusters(Int_t n) {fN=n;}
  88 protected:
  89   void     SaveLocalConvConst();
  90   Double_t GetLocalConvConst() const {return fLocalConvConst;}
  91
  92   void External2Helix(Double_t helix[6]) const;
  93
  94   void SetChi2(Double_t chi2) {fChi2=chi2;}
  95   void SetMass(Double_t mass) {fMass=mass;}
  96
  97   Int_t fLab;             // track label
  98   Float_t fFakeRatio;     // fake ratio
  99   Double_t fChi2;         // total chi2 value for this track
 100   Double_t fMass;         // mass hypothesis
 101   Int_t fN;               // number of associated clusters
 102
 103 private:
 104   Double_t fLocalConvConst;   //local conversion "curvature(1/cm) -> pt(GeV/c)"
 105
 106   // variables for time integration (S.Radomski@gsi.de)
 107   Bool_t  fStartTimeIntegral;       // indicator wether integrate time
 108   Double_t fIntegratedTime[AliPID::kSPECIES];       // integrated time
 109   Double_t fIntegratedLength;        // integrated length
 110
 111   ClassDef(AliKalmanTrack,5)    // Reconstructed track
 112 };
 113
 114 #endif
 115
 116