]>
Commit | Line | Data |
---|---|---|
4cc511f4 AU |
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" | |
e21f8bf5 | 23 | #include "AliESDMuonGlobalTrack.h" |
4cc511f4 AU |
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 | ||
4cc511f4 AU |
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 | ||
65273893 | 56 | Bool_t IsCorrectMatch(AliMFTCluster *cluster, Int_t labelMC); |
57 | ||
58 | void GetVertexFromMC(); | |
59 | ||
4cc511f4 AU |
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 | ||
7bafb008 | 67 | void AddClustersFromUnderlyingEvent(); |
68 | void AddClustersFromPileUpEvents(); | |
4cc511f4 AU |
69 | |
70 | protected: | |
71 | ||
72 | static const Int_t fNMaxPlanes = AliMFTConstants::fNMaxPlanes; // max number of MFT planes | |
73 | static const Double_t fRadLengthSi; | |
74 | static const Int_t fMaxNCandidates = 1000; | |
75 | ||
76 | AliESDEvent *fESD; //! pointer to the ESD event | |
77 | ||
78 | AliMFT *fMFT; //! | |
79 | AliMFTSegmentation *fSegmentation; //! | |
80 | Int_t fNPlanesMFT, fNPlanesMFTAnalyzed; | |
81 | ||
82 | Double_t fSigmaClusterCut; // to select the clusters in the MFT planes which are compatible with the extrapolated muon track | |
83 | Double_t fScaleSigmaClusterCut; // to tune the cut on the compatible clusters in case of too many candidates | |
84 | ||
85 | Int_t fNMaxMissingMFTClusters; // max. number of MFT clusters which can be missed in the global fit procedure | |
86 | Bool_t fIsPlaneMandatory[fNMaxPlanes]; // specifies which MFT planes cannot be missed in the global fit procedure | |
87 | ||
88 | Bool_t fGlobalTrackingDiverged; // to keep memory of a possible divergence in the global tracking finding | |
89 | ||
90 | TClonesArray *fMFTClusterArray[fNMaxPlanes]; //! array of clusters for the planes of the MFT | |
91 | TClonesArray *fMFTClusterArrayFront[fNMaxPlanes]; //! array of front clusters for the planes of the MFT | |
92 | TClonesArray *fMFTClusterArrayBack[fNMaxPlanes]; //! array of back clusters for the planes of the MFT | |
93 | ||
94 | TClonesArray *fCandidateTracks; //! array of candidate global tracks | |
95 | ||
96 | AliMUONTrack *fMUONTrack; //! muon track being analyzed | |
97 | ||
98 | AliMuonForwardTrack *fCurrentTrack; //! muon extrapolated track being tested | |
99 | AliMuonForwardTrack *fFinalBestCandidate; //! best final candidate (if any) | |
100 | ||
101 | // uncertainty on the x position of the primary vertex (seed for the extrapolation of the MUON tracks) | |
65273893 | 102 | Double_t fXExtrapVertex; |
103 | Double_t fYExtrapVertex; | |
104 | Double_t fZExtrapVertex; | |
105 | Double_t fXExtrapVertexError; | |
106 | Double_t fYExtrapVertexError; | |
4cc511f4 AU |
107 | |
108 | Bool_t fBransonCorrection; // if TRUE, Branson Correction is applied when extrapolating the MUON tracks to the vertex region | |
109 | ||
110 | Double_t fMinResearchRadiusAtPlane[fNMaxPlanes]; | |
111 | ||
112 | private: | |
113 | ||
114 | AliMFTTrackerMU(const AliMFTTrackerMU &tracker); | |
115 | AliMFTTrackerMU & operator=(const AliMFTTrackerMU &tracker); | |
65273893 | 116 | ClassDef(AliMFTTrackerMU,2) //MFT tracker for muon tracks |
4cc511f4 AU |
117 | |
118 | }; | |
119 | ||
120 | //==================================================================================================================================================== | |
121 | ||
122 | #endif | |
123 |