1 #ifndef AliMFTTrackerMU_H
2 #define AliMFTTrackerMU_H
4 /* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 //====================================================================================================================================================
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
14 // Contact author: antonio.uras@cern.ch
16 //====================================================================================================================================================
20 #include "AliGeomManager.h"
21 #include "AliESDEvent.h"
22 #include "AliESDMuonTrack.h"
23 #include "AliESDMuonGlobalTrack.h"
24 #include "AliTracker.h"
27 #include "AliMUONTrackExtrap.h"
28 #include "AliMUONTrack.h"
29 #include "AliMUONESDInterface.h"
30 #include "AliMuonForwardTrack.h"
32 //====================================================================================================================================================
34 class AliMFTTrackerMU : public AliTracker {
38 enum {kConverged, kDiverged};
41 virtual ~AliMFTTrackerMU();
43 Int_t LoadClusters(TTree *cf);
44 void UnloadClusters();
45 Int_t Clusters2Tracks(AliESDEvent *event);
47 void SetNPlanesMFT(Int_t nPlanesMFT) { fNPlanesMFT = nPlanesMFT; }
48 void SeparateFrontBackClusters();
50 Int_t FindClusterInPlane(Int_t planeId);
52 void SetMinResearchRadiusAtPlane(Int_t plane, Double_t radius) { if (plane>=0 && plane<fNMaxPlanes) fMinResearchRadiusAtPlane[plane] = radius; }
54 Double_t TryOneCluster(const AliMUONTrackParam &trackParam, AliMFTCluster *cluster);
56 Bool_t IsCorrectMatch(AliMFTCluster *cluster, Int_t labelMC);
58 void GetVertexFromMC();
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;}
67 void AddClustersFromUnderlyingEvent();
68 void AddClustersFromPileUpEvents();
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;
76 AliESDEvent *fESD; //! pointer to the ESD event
79 AliMFTSegmentation *fSegmentation; //!
80 Int_t fNPlanesMFT, fNPlanesMFTAnalyzed;
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
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
88 Bool_t fGlobalTrackingDiverged; // to keep memory of a possible divergence in the global tracking finding
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
94 TClonesArray *fCandidateTracks; //! array of candidate global tracks
96 AliMUONTrack *fMUONTrack; //! muon track being analyzed
98 AliMuonForwardTrack *fCurrentTrack; //! muon extrapolated track being tested
99 AliMuonForwardTrack *fFinalBestCandidate; //! best final candidate (if any)
101 // uncertainty on the x position of the primary vertex (seed for the extrapolation of the MUON tracks)
102 Double_t fXExtrapVertex;
103 Double_t fYExtrapVertex;
104 Double_t fZExtrapVertex;
105 Double_t fXExtrapVertexError;
106 Double_t fYExtrapVertexError;
108 Bool_t fBransonCorrection; // if TRUE, Branson Correction is applied when extrapolating the MUON tracks to the vertex region
110 Double_t fMinResearchRadiusAtPlane[fNMaxPlanes];
114 AliMFTTrackerMU(const AliMFTTrackerMU &tracker);
115 AliMFTTrackerMU & operator=(const AliMFTTrackerMU &tracker);
116 ClassDef(AliMFTTrackerMU,2) //MFT tracker for muon tracks
120 //====================================================================================================================================================