]> git.uio.no Git - u/mrichter/AliRoot.git/blob - STEER/AliESDMuonTrack.h
Removing osolete variable. Adding title to trigger configuration obj also
[u/mrichter/AliRoot.git] / STEER / AliESDMuonTrack.h
1 #ifndef ALIESDMUONTRACK_H
2 #define ALIESDMUONTRACK_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 //  Class to describe the MUON tracks
10 //  in the Event Summary Data class
11 //  Author: G.Martinez
12
13
14 #include <TMath.h>
15 #include <TMatrixD.h>
16 #include <TDatabasePDG.h>
17
18 #include "AliVParticle.h"
19
20 class AliESDMuonCluster;
21 class TClonesArray;
22 class TLorentzVector;
23
24 class AliESDMuonTrack : public AliVParticle {
25 public:
26   AliESDMuonTrack(); //Constructor
27   virtual ~AliESDMuonTrack(); // Destructor
28   AliESDMuonTrack(const AliESDMuonTrack& esdm);
29   AliESDMuonTrack& operator=(const AliESDMuonTrack& esdm);
30   virtual void Copy(TObject &obj) const;
31
32   virtual void Clear(Option_t* opt = "");
33   
34   // Get and Set methods for data at vertex
35   Double_t GetInverseBendingMomentum(void) const {return fInverseBendingMomentum;}
36   void     SetInverseBendingMomentum(Double_t InverseBendingMomentum) 
37                 {fInverseBendingMomentum = InverseBendingMomentum;}
38   Double_t GetThetaX(void) const {return fThetaX;}
39   void     SetThetaX(Double_t ThetaX) {fThetaX = ThetaX;}
40   Double_t GetThetaY(void) const {return fThetaY;}
41   void     SetThetaY(Double_t ThetaY) {fThetaY = ThetaY;}
42   Double_t GetZ(void) const {return fZ;}
43   void     SetZ(Double_t Z) {fZ = Z;}
44   Double_t GetBendingCoor(void) const {return fBendingCoor;}
45   void     SetBendingCoor(Double_t BendingCoor) {fBendingCoor = BendingCoor;}
46   Double_t GetNonBendingCoor(void) const {return fNonBendingCoor;}
47   void     SetNonBendingCoor(Double_t NonBendingCoor) {fNonBendingCoor = NonBendingCoor;}
48   
49   // Get and Set methods for data at Distance of Closest Approach in the vertex plane
50   Double_t GetInverseBendingMomentumAtDCA(void) const {return fInverseBendingMomentumAtDCA;}
51   void     SetInverseBendingMomentumAtDCA(Double_t InverseBendingMomentum) 
52                 {fInverseBendingMomentumAtDCA = InverseBendingMomentum;}
53   Double_t GetThetaXAtDCA(void) const {return fThetaXAtDCA;}
54   void     SetThetaXAtDCA(Double_t ThetaX) {fThetaXAtDCA = ThetaX;}
55   Double_t GetThetaYAtDCA(void) const {return fThetaYAtDCA;}
56   void     SetThetaYAtDCA(Double_t ThetaY) {fThetaYAtDCA = ThetaY;}
57   Double_t GetBendingCoorAtDCA(void) const {return fBendingCoorAtDCA;}
58   void     SetBendingCoorAtDCA(Double_t BendingCoor) {fBendingCoorAtDCA = BendingCoor;}
59   Double_t GetNonBendingCoorAtDCA(void) const {return fNonBendingCoorAtDCA;}
60   void     SetNonBendingCoorAtDCA(Double_t NonBendingCoor) {fNonBendingCoorAtDCA = NonBendingCoor;}
61   Double_t GetDCA(void) const {return TMath::Sqrt(fNonBendingCoorAtDCA*fNonBendingCoorAtDCA +
62                                                   fBendingCoorAtDCA*fBendingCoorAtDCA);}
63   
64   // Get and Set methods for data at first station
65   Double_t GetInverseBendingMomentumUncorrected(void) const {return fInverseBendingMomentumUncorrected;}
66   void     SetInverseBendingMomentumUncorrected(Double_t InverseBendingMomentum) 
67                 {fInverseBendingMomentumUncorrected = InverseBendingMomentum;}
68   Double_t GetThetaXUncorrected(void) const {return fThetaXUncorrected;}
69   void     SetThetaXUncorrected(Double_t ThetaX) {fThetaXUncorrected = ThetaX;}
70   Double_t GetThetaYUncorrected(void) const {return fThetaYUncorrected;}
71   void     SetThetaYUncorrected(Double_t ThetaY) {fThetaYUncorrected = ThetaY;}
72   Double_t GetZUncorrected(void) const {return fZUncorrected;}
73   void     SetZUncorrected(Double_t Z) {fZUncorrected = Z;}
74   Double_t GetBendingCoorUncorrected(void) const {return fBendingCoorUncorrected;}
75   void     SetBendingCoorUncorrected(Double_t BendingCoor) {fBendingCoorUncorrected = BendingCoor;}
76   Double_t GetNonBendingCoorUncorrected(void) const {return fNonBendingCoorUncorrected;}
77   void     SetNonBendingCoorUncorrected(Double_t NonBendingCoor) {fNonBendingCoorUncorrected = NonBendingCoor;}
78   
79   // Get and Set methods for covariance matrix of data at first station
80   void     GetCovariances(TMatrixD& cov) const;
81   void     SetCovariances(const TMatrixD& cov);
82   void     GetCovarianceXYZPxPyPz(Double_t cov[21]) const;
83   
84   // Get and Set methods for global tracking info
85   Double_t GetChi2(void) const {return fChi2;}
86   void     SetChi2(Double_t Chi2) {fChi2 = Chi2;}
87   UChar_t  GetNHit(void) const {return fNHit;}
88   void     SetNHit(UInt_t NHit) {fNHit = NHit;}
89   
90   // Get and Set methods for trigger matching
91   Int_t    GetMatchTrigger() const;
92   Double_t GetChi2MatchTrigger() const {return fChi2MatchTrigger;}
93   void     SetChi2MatchTrigger(Double_t Chi2MatchTrigger) {fChi2MatchTrigger = Chi2MatchTrigger;}
94   UShort_t GetHitsPatternInTrigCh() const {return fHitsPatternInTrigCh;}
95   void     SetHitsPatternInTrigCh(UShort_t hitsPatternInTrigCh) {fHitsPatternInTrigCh = hitsPatternInTrigCh;}
96   void     SetLocalTrigger(Int_t locTrig) { fLocalTrigger = locTrig; }
97   Int_t    LoCircuit(void) const { return fLocalTrigger & 0xFF;       }
98   Int_t    LoStripX(void) const  { return fLocalTrigger >>  8 & 0x1F; }
99   Int_t    LoStripY(void) const  { return fLocalTrigger >> 13 & 0x0F; }
100   Int_t    LoDev(void)    const  { return fLocalTrigger >> 17 & 0x1F; }
101   Int_t    LoLpt(void)    const  { return fLocalTrigger >> 22 & 0x03; }
102   Int_t    LoHpt(void)    const  { return fLocalTrigger >> 24 & 0x03; }
103   
104   // Get and Set methods for the hit strips pattern in the trigger chambers
105   UShort_t GetTriggerX1Pattern() { return fX1Pattern; }
106   UShort_t GetTriggerY1Pattern() { return fY1Pattern; }
107   UShort_t GetTriggerX2Pattern() { return fX2Pattern; }
108   UShort_t GetTriggerY2Pattern() { return fY2Pattern; }
109   UShort_t GetTriggerX3Pattern() { return fX3Pattern; }
110   UShort_t GetTriggerY3Pattern() { return fY3Pattern; }
111   UShort_t GetTriggerX4Pattern() { return fX4Pattern; }
112   UShort_t GetTriggerY4Pattern() { return fY4Pattern; }
113   void     SetTriggerX1Pattern(UShort_t pat) { fX1Pattern = pat; }
114   void     SetTriggerY1Pattern(UShort_t pat) { fY1Pattern = pat; }
115   void     SetTriggerX2Pattern(UShort_t pat) { fX2Pattern = pat; }
116   void     SetTriggerY2Pattern(UShort_t pat) { fY2Pattern = pat; }
117   void     SetTriggerX3Pattern(UShort_t pat) { fX3Pattern = pat; }
118   void     SetTriggerY3Pattern(UShort_t pat) { fY3Pattern = pat; }
119   void     SetTriggerX4Pattern(UShort_t pat) { fX4Pattern = pat; }
120   void     SetTriggerY4Pattern(UShort_t pat) { fY4Pattern = pat; }
121
122   // Get and Set methods for muon cluster map
123   UInt_t   GetMuonClusterMap() const {return fMuonClusterMap;}
124   void     SetMuonClusterMap(UInt_t muonClusterMap) {fMuonClusterMap = muonClusterMap;}
125   void     AddInMuonClusterMap(Int_t chamber);
126   Bool_t   IsInMuonClusterMap(Int_t chamber) const;
127   
128   // Methods to get, fill and check the array of associated clusters
129   Int_t         GetNClusters() const;
130   TClonesArray& GetClusters() const;
131   void          AddCluster(const AliESDMuonCluster &cluster);
132   Bool_t        ClustersStored() const;
133   
134   // Methods to compute track momentum
135   Double_t Px() const;
136   Double_t Py() const;
137   Double_t Pz() const;
138   Double_t P() const;
139   Bool_t   PxPyPz(Double_t p[3]) const { p[0] = Px(); p[1] = Py(); p[2] = Pz(); return kTRUE; }
140   void     LorentzP(TLorentzVector& vP) const;
141   Double_t PxAtDCA() const;
142   Double_t PyAtDCA() const;
143   Double_t PzAtDCA() const;
144   Double_t PAtDCA() const;
145   Bool_t   PxPyPzAtDCA(Double_t p[3]) const { p[0] = Px(); p[1] = Py(); p[2] = Pz(); return kTRUE; }
146   void     LorentzPAtDCA(TLorentzVector& vP) const;
147   Double_t PxUncorrected() const;
148   Double_t PyUncorrected() const;
149   Double_t PzUncorrected() const;
150   Double_t PUncorrected() const;
151   Bool_t   PxPyPzUncorrected(Double_t p[3]) const { p[0] = Px(); p[1] = Py(); p[2] = Pz(); return kTRUE; }
152   void     LorentzPUncorrected(TLorentzVector& vP) const;
153   
154   // additional methods to comply with AliVParticle
155   Double_t Xv() const {return -999.;} // put reasonable values here
156   Double_t Yv() const {return -999.;} //
157   Double_t Zv() const {return -999.;} //
158   Bool_t   XvYvZv(Double_t x[3]) const { x[0] = Xv(); x[1] = Yv(); x[2] = Zv(); return kTRUE; }  
159   Double_t Pt() const { return TMath::Sqrt(Px()*Px() + Py()*Py()); }
160   Double_t OneOverPt() const { return (Pt() != 0.) ? 1./Pt() : FLT_MAX; }
161   Double_t Phi() const { return TMath::ATan2(Py(), Px()); }
162   Double_t Theta() const { return TMath::ATan2(Pt(), Pz()); }
163   Double_t E() const { return TMath::Sqrt(M()*M() + P()*P()); }
164   Double_t M() const { return TDatabasePDG::Instance()->GetParticle("mu-")->Mass(); }
165   Double_t Eta() const { return -TMath::Log(TMath::Tan(0.5 * Theta()));}
166   Double_t Y() const { return (Pz()/E() != 1.) ? TMath::ATanH(Pz()/E()) : FLT_MAX; }
167   Short_t  Charge() const { return (Short_t)TMath::Sign(1., GetInverseBendingMomentum()); }
168   const Double_t *PID() const { return (Double_t*)0x0; }
169   
170   
171 protected:
172   // parameters at vertex
173   Double32_t fInverseBendingMomentum; ///< Inverse bending momentum (GeV/c ** -1) times the charge 
174   Double32_t fThetaX;                 ///< Angle of track at vertex in X direction (rad)
175   Double32_t fThetaY;                 ///< Angle of track at vertex in Y direction (rad)
176   Double32_t fZ;                      ///< Z coordinate (cm)
177   Double32_t fBendingCoor;            ///< bending coordinate (cm)
178   Double32_t fNonBendingCoor;         ///< non bending coordinate (cm)
179   
180   // parameters at Distance of Closest Approach in the vertex plane
181   Double32_t fInverseBendingMomentumAtDCA; ///< Inverse bending momentum (GeV/c ** -1) times the charge 
182   Double32_t fThetaXAtDCA;                 ///< Angle of track at vertex in X direction (rad)
183   Double32_t fThetaYAtDCA;                 ///< Angle of track at vertex in Y direction (rad)
184   Double32_t fBendingCoorAtDCA;            ///< bending coordinate (cm)
185   Double32_t fNonBendingCoorAtDCA;         ///< non bending coordinate (cm)
186   
187   // parameters at first tracking station
188   Double32_t fInverseBendingMomentumUncorrected; ///< Inverse bending momentum (GeV/c ** -1) times the charge 
189   Double32_t fThetaXUncorrected;                 ///< Angle of track at vertex in X direction (rad)
190   Double32_t fThetaYUncorrected;                 ///< Angle of track at vertex in Y direction (rad)
191   Double32_t fZUncorrected;                      ///< Z coordinate (cm)
192   Double32_t fBendingCoorUncorrected;            ///< bending coordinate (cm)
193   Double32_t fNonBendingCoorUncorrected;         ///< non bending coordinate (cm)
194   
195   /// reduced covariance matrix of UNCORRECTED track parameters, ordered as follow:      <pre>
196   /// [0] =  <X,X>
197   /// [1] =<X,ThetaX>  [2] =<ThetaX,ThetaX>
198   /// [3] =  <X,Y>     [4] =  <Y,ThetaX>     [5] =  <Y,Y>
199   /// [6] =<X,ThetaY>  [7] =<ThetaX,ThetaY>  [8] =<Y,ThetaY>  [9] =<ThetaY,ThetaY>
200   /// [10]=<X,InvP_yz> [11]=<ThetaX,InvP_yz> [12]=<Y,InvP_yz> [13]=<ThetaY,InvP_yz> [14]=<InvP_yz,InvP_yz>  </pre>
201   Double32_t fCovariances[15]; ///< \brief reduced covariance matrix of parameters AT FIRST CHAMBER
202   
203   // global tracking info
204   Double32_t fChi2;                ///< chi2 in the MUON track fit
205   Double32_t fChi2MatchTrigger;    ///< chi2 of trigger/track matching
206   Int_t      fLocalTrigger;        ///< packed local trigger information
207
208   // hit strips pattern in the trigger chambers
209   UShort_t fX1Pattern;             ///< x-strips pattern in st6/ch1
210   UShort_t fY1Pattern;             ///< y-strips pattern in st6/ch1
211   UShort_t fX2Pattern;             ///< x-strips pattern in st6/ch2
212   UShort_t fY2Pattern;             ///< y-strips pattern in st6/ch2
213   UShort_t fX3Pattern;             ///< x-strips pattern in st7/ch1
214   UShort_t fY3Pattern;             ///< y-strips pattern in st7/ch1
215   UShort_t fX4Pattern;             ///< x-strips pattern in st7/ch2
216   UShort_t fY4Pattern;             ///< y-strips pattern in st7/ch2
217   
218   UInt_t     fMuonClusterMap;      ///< Map of clusters in tracking chambers
219   UShort_t   fHitsPatternInTrigCh; ///< Word containing info on the hits left in trigger chambers
220   UChar_t    fNHit;                ///< number of hit in the track
221   
222   mutable TClonesArray* fClusters; ///< Array of clusters attached to the track
223   
224   ClassDef(AliESDMuonTrack,10) // MUON ESD track class 
225 };
226
227 #endif