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