Merge branch 'master' of https://git.cern.ch/reps/AliRoot
[u/mrichter/AliRoot.git] / MFT / AliMFTTrackerMU.h
1 #ifndef AliMFTTrackerMU_H
2 #define AliMFTTrackerMU_H
3
4 /* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice                               */
6
7 //====================================================================================================================================================
8 //
9 //                          MFT tracker
10 //
11 // Class for the creation of the "global muon tracks" built from the tracks reconstructed in the 
12 // muon spectrometer and the clusters of the Muon Forward Tracker
13 //
14 // Contact author: antonio.uras@cern.ch
15 //
16 //====================================================================================================================================================
17
18 #include "TTree.h"
19 #include "AliLog.h"
20 #include "AliGeomManager.h"
21 #include "AliESDEvent.h"
22 #include "AliESDMuonTrack.h"
23 #include "AliESDMuonGlobalTrack.h"
24 #include "AliTracker.h"
25 #include "AliRun.h"
26 #include "AliMFT.h"
27 #include "AliMUONTrackExtrap.h"
28 #include "AliMUONTrack.h"
29 #include "AliMUONESDInterface.h"
30 #include "AliMuonForwardTrack.h"
31
32 //====================================================================================================================================================
33
34 class AliMFTTrackerMU : public AliTracker {
35
36 public:
37
38   enum {kConverged, kDiverged};
39
40   AliMFTTrackerMU();
41   virtual ~AliMFTTrackerMU();
42
43   Int_t LoadClusters(TTree *cf);
44   void UnloadClusters();
45   Int_t Clusters2Tracks(AliESDEvent *event);
46
47   void SetNPlanesMFT(Int_t nPlanesMFT) { fNPlanesMFT = nPlanesMFT; }
48   void SeparateFrontBackClusters();
49
50   Int_t FindClusterInPlane(Int_t planeId);
51
52   void SetMinResearchRadiusAtPlane(Int_t plane, Double_t radius) { if (plane>=0 && plane<fNMaxPlanes) fMinResearchRadiusAtPlane[plane] = radius; }
53
54   Double_t TryOneCluster(const AliMUONTrackParam &trackParam, AliMFTCluster *cluster);
55
56   Bool_t IsCorrectMatch(AliMFTCluster *cluster, Int_t labelMC);
57
58   void GetVertexFromMC();
59
60   /// Dummy implementation
61   virtual Int_t PropagateBack(AliESDEvent* /*event*/) {return 0;}
62   /// Dummy implementation
63   virtual Int_t RefitInward(AliESDEvent* /*event*/) {return 0;}
64   /// Dummy implementation
65   virtual AliCluster *GetCluster(Int_t /*index*/) const {return 0;}
66
67
68 protected:
69
70   static const Int_t fNMaxPlanes = AliMFTConstants::fNMaxPlanes;        // max number of MFT planes
71   static const Double_t fRadLengthSi;
72   static const Int_t fMaxNCandidates = 1000;
73
74   AliESDEvent *fESD;   //! pointer to the ESD event
75
76   AliMFT *fMFT;                        //!
77   AliMFTSegmentation *fSegmentation;   //!
78   Int_t fNPlanesMFT, fNPlanesMFTAnalyzed;
79
80   Double_t fSigmaClusterCut;         // to select the clusters in the MFT planes which are compatible with the extrapolated muon track
81   Double_t fScaleSigmaClusterCut;    // to tune the cut on the compatible clusters in case of too many candidates
82
83   Int_t fNMaxMissingMFTClusters;              // max. number of MFT clusters which can be missed in the global fit procedure
84   Bool_t fIsPlaneMandatory[fNMaxPlanes];      // specifies which MFT planes cannot be missed in the global fit procedure
85
86   Bool_t fGlobalTrackingDiverged;    // to keep memory of a possible divergence in the global tracking finding
87
88   TClonesArray *fMFTClusterArray[fNMaxPlanes];         //! array of clusters for the planes of the MFT
89   TClonesArray *fMFTClusterArrayFront[fNMaxPlanes];    //! array of front clusters for the planes of the MFT
90   TClonesArray *fMFTClusterArrayBack[fNMaxPlanes];     //! array of back clusters for the planes of the MFT
91
92   TClonesArray *fCandidateTracks;   //! array of candidate global tracks 
93
94   AliMUONTrack *fMUONTrack;  //! muon track being analyzed
95
96   AliMuonForwardTrack *fCurrentTrack;     //! muon extrapolated track being tested
97   AliMuonForwardTrack *fFinalBestCandidate;     //! best final candidate (if any)
98
99   // uncertainty on the x position of the primary vertex (seed for the extrapolation of the MUON tracks)
100   Double_t fXExtrapVertex;   
101   Double_t fYExtrapVertex;
102   Double_t fZExtrapVertex;
103   Double_t fXExtrapVertexError;
104   Double_t fYExtrapVertexError;
105
106   Bool_t fBransonCorrection;    // if TRUE, Branson Correction is applied when extrapolating the MUON tracks to the vertex region
107
108   Double_t fMinResearchRadiusAtPlane[fNMaxPlanes];
109
110 private:
111
112   AliMFTTrackerMU(const AliMFTTrackerMU &tracker);
113   AliMFTTrackerMU & operator=(const AliMFTTrackerMU &tracker);
114   ClassDef(AliMFTTrackerMU,2)   //MFT tracker for muon tracks
115
116 };
117
118 //====================================================================================================================================================
119
120 #endif 
121