1 #ifndef ALIKALMANTRACK_H
2 #define ALIKALMANTRACK_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
9 //-------------------------------------------------------------------------
10 // Class AliKalmanTrack
12 // Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
13 //-------------------------------------------------------------------------
19 class AliKalmanTrack : public TObject {
22 AliKalmanTrack(const AliKalmanTrack &t);
24 virtual ~AliKalmanTrack(){};
25 void SetLabel(Int_t lab) {fLab=lab;}
27 Bool_t IsSortable() const {return kTRUE;}
28 Int_t GetLabel() const {return fLab;}
29 Double_t GetChi2() const {return fChi2;}
30 Double_t GetMass() const {return fMass;}
31 Int_t GetNumberOfClusters() const {return fN;}
32 virtual Int_t GetClusterIndex(Int_t) const { //reserved for AliTracker
33 Warning("GetClusterIndex(Int_t)","Method must be overloaded !\n");
37 virtual Double_t GetDCA(const AliKalmanTrack *,Double_t &,Double_t &) const;
39 Double_t PropagateToDCA(AliKalmanTrack *, Double_t d=0., Double_t x0=0.);
40 virtual Double_t GetAlpha() const {
41 Warning("GetAlpha()","Method must be overloaded !\n");
44 virtual Double_t GetSigmaY2() const {
45 Warning("GetSigmaY2()","Method must be overloaded !\n");
48 virtual Double_t GetSigmaZ2() const {
49 Warning("GetSigmaZ2()","Method must be overloaded !\n");
53 virtual Int_t Compare(const TObject *) const {return 0;}
55 virtual void GetExternalParameters(Double_t &/*xr*/, Double_t /*x*/[5]) const {}
56 virtual void GetExternalCovariance(Double_t /*cov*/[15]) const {}
58 virtual Double_t GetPredictedChi2(const AliCluster *) const {return 0.;}
60 Int_t PropagateTo(Double_t /*xr*/, Double_t /*x0*/, Double_t /*rho*/) {return 0;}
61 virtual Int_t Update(const AliCluster*, Double_t /*chi2*/, UInt_t) {return 0;}
63 static void SetConvConst(Double_t cc) {fgConvConst=cc;}
64 Double_t GetConvConst() const {return fgConvConst;}
66 static void SetMagneticField(Double_t f) {// f - Magnetic field in T
67 fgConvConst=100/0.299792458/f;
69 Double_t GetMagneticField() const {return 100/0.299792458/fgConvConst;}
71 // Time integration (S.Radomski@gsi.de)
72 void StartTimeIntegral();
73 Bool_t IsStartedTimeIntegral() const {return fStartTimeIntegral;}
74 void AddTimeStep(Double_t length);
75 Double_t GetIntegratedTime(Int_t pdg) const;
76 Double_t GetIntegratedLength() const {return fIntegratedLength;}
77 void PrintTime() const;
81 void SetChi2(Double_t chi2) {fChi2=chi2;}
82 void SetMass(Double_t mass) {fMass=mass;}
83 void SetNumberOfClusters(Int_t n) {fN=n;}
86 Int_t fLab; // track label
87 Double_t fChi2; // total chi2 value for this track
88 Double_t fMass; // mass hypothesis
89 Int_t fN; // number of associated clusters
91 static Double_t fgConvConst; //conversion constant cm -> GeV/c
93 // variables for time integration (S.Radomski@gsi.de)
94 static const Int_t fgkTypes = 5; // Number of track types (e,mu,pi,k,p)
95 Bool_t fStartTimeIntegral; // indicator wether integrate time
96 Float_t fIntegratedTime[5]; // intgrated time
97 Float_t fIntegratedLength; // integrated length
99 ClassDef(AliKalmanTrack,2) // Reconstructed track