]> git.uio.no Git - u/mrichter/AliRoot.git/blob - STEER/AliKalmanTrack.h
Boost for assymmetric systems added.
[u/mrichter/AliRoot.git] / 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 //
12 //         Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch 
13 //-------------------------------------------------------------------------
14
15 #include <TObject.h>
16 #include <TVector3.h>
17
18 class AliCluster;
19
20 class AliKalmanTrack : public TObject {
21 public:
22   AliKalmanTrack();
23   AliKalmanTrack(const AliKalmanTrack &t);
24
25   virtual ~AliKalmanTrack(){};
26   void SetLabel(Int_t lab) {fLab=lab;}
27
28   Bool_t   IsSortable() const {return kTRUE;}
29   Int_t    GetLabel()   const {return fLab;}
30   Double_t GetChi2()    const {return fChi2;}
31   Double_t GetMass()    const {return fMass;}
32   Int_t    GetNumberOfClusters() const {return fN;}
33   virtual Int_t GetClusterIndex(Int_t) const { //reserved for AliTracker
34     Warning("GetClusterIndex(Int_t)","Method must be overloaded !\n");
35     return 0;
36   } 
37
38   virtual Double_t GetDCA(const AliKalmanTrack *,Double_t &,Double_t &) const; 
39   virtual 
40   Double_t PropagateToDCA(AliKalmanTrack *, Double_t d=0., Double_t x0=0.); 
41   virtual Double_t GetAlpha() const {
42     Warning("GetAlpha()","Method must be overloaded !\n");
43     return 0.;
44   }
45   virtual Double_t GetSigmaY2() const {
46     Warning("GetSigmaY2()","Method must be overloaded !\n");
47     return 0.;
48   }
49   virtual Double_t GetSigmaZ2() const {
50     Warning("GetSigmaZ2()","Method must be overloaded !\n");
51     return 0.;
52   }
53
54   virtual Int_t Compare(const TObject *) const {return 0;} 
55
56   virtual void GetExternalParameters(Double_t &/*xr*/, Double_t /*x*/[5]) const {}
57   virtual void GetExternalCovariance(Double_t /*cov*/[15]) const {}
58
59   virtual Double_t GetX() const;
60   virtual Double_t GetdEdx() const;
61
62   virtual Double_t GetY() const;
63   virtual Double_t GetZ() const;
64   virtual Double_t GetSnp() const;
65   virtual Double_t GetTgl() const;
66   virtual Double_t Get1Pt() const;
67
68   virtual Double_t Phi() const;
69   virtual Double_t SigmaPhi() const;
70   virtual Double_t Theta() const;
71   virtual Double_t SigmaTheta() const;
72   virtual Double_t Px() const;
73   virtual Double_t Py() const;
74   virtual Double_t Pz() const;
75   virtual Double_t Pt() const;
76   virtual Double_t SigmaPt() const;
77   virtual Double_t P() const;
78   virtual TVector3 Momentum() const;
79
80   virtual Double_t GetPredictedChi2(const AliCluster *) const {return 0.;}
81   virtual 
82     Int_t PropagateTo(Double_t /*xr*/, Double_t /*x0*/, Double_t /*rho*/) {return 0;}
83   virtual Int_t Update(const AliCluster*, Double_t /*chi2*/, UInt_t) {return 0;}
84
85   static void SetConvConst(Double_t cc) {fgConvConst=cc;}
86   Double_t GetConvConst() const {return fgConvConst;}
87
88   static void SetMagneticField(Double_t f) {// f - Magnetic field in T
89     fgConvConst=100/0.299792458/f;
90   }
91   Double_t GetMagneticField() const {return 100/0.299792458/fgConvConst;}
92
93   // Time integration (S.Radomski@gsi.de)
94   void   StartTimeIntegral();
95   Bool_t IsStartedTimeIntegral() const {return fStartTimeIntegral;}
96   void     AddTimeStep(Double_t length);
97   Double_t GetIntegratedTime(Int_t pdg) const;
98   Double_t GetIntegratedLength() const {return fIntegratedLength;}
99   void PrintTime() const;
100   
101
102 protected:
103   void SetChi2(Double_t chi2) {fChi2=chi2;} 
104   void SetMass(Double_t mass) {fMass=mass;}
105   void SetNumberOfClusters(Int_t n) {fN=n;} 
106
107 private: 
108   Int_t fLab;             // track label
109   Double_t fChi2;         // total chi2 value for this track
110   Double_t fMass;         // mass hypothesis
111   Int_t fN;               // number of associated clusters
112
113   static Double_t fgConvConst; //conversion constant cm -> GeV/c
114
115   // variables for time integration (S.Radomski@gsi.de)
116   static const Int_t fgkTypes = 5;  // Number of track types (e,mu,pi,k,p)
117   Bool_t  fStartTimeIntegral;       // indicator wether integrate time
118   Float_t fIntegratedTime[5];       // intgrated time
119   Float_t fIntegratedLength;        // integrated length
120   
121   ClassDef(AliKalmanTrack,2)    // Reconstructed track
122 };
123
124 #endif
125
126