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 SetVertexError(Double_t xErr, Double_t yErr, Double_t zErr) { fVertexErrorX=xErr; fVertexErrorY=yErr; fVertexErrorZ=zErr; }
  53   void SetMinResearchRadiusAtPlane(Int_t plane, Double_t radius) { if (plane>=0 && plane<fNMaxPlanes) fMinResearchRadiusAtPlane[plane] = radius; }
  54
  55   Double_t TryOneCluster(const AliMUONTrackParam &trackParam, AliMFTCluster *cluster);
  56
  57   /// Dummy implementation
  58   virtual Int_t PropagateBack(AliESDEvent* /*event*/) {return 0;}
  59   /// Dummy implementation
  60   virtual Int_t RefitInward(AliESDEvent* /*event*/) {return 0;}
  61   /// Dummy implementation
  62   virtual AliCluster *GetCluster(Int_t /*index*/) const {return 0;}
  63
  64
  65 protected:
  66
  67   static const Int_t fNMaxPlanes = AliMFTConstants::fNMaxPlanes;        // max number of MFT planes
  68   static const Double_t fRadLengthSi;
  69   static const Int_t fMaxNCandidates = 1000;
  70
  71   AliESDEvent *fESD;   //! pointer to the ESD event
  72
  73   AliMFT *fMFT;                        //!
  74   AliMFTSegmentation *fSegmentation;   //!
  75   Int_t fNPlanesMFT, fNPlanesMFTAnalyzed;
  76
  77   Double_t fSigmaClusterCut;         // to select the clusters in the MFT planes which are compatible with the extrapolated muon track
  78   Double_t fScaleSigmaClusterCut;    // to tune the cut on the compatible clusters in case of too many candidates
  79
  80   Int_t fNMaxMissingMFTClusters;              // max. number of MFT clusters which can be missed in the global fit procedure
  81   Bool_t fIsPlaneMandatory[fNMaxPlanes];      // specifies which MFT planes cannot be missed in the global fit procedure
  82
  83   Bool_t fGlobalTrackingDiverged;    // to keep memory of a possible divergence in the global tracking finding
  84
  85   TClonesArray *fMFTClusterArray[fNMaxPlanes];         //! array of clusters for the planes of the MFT
  86   TClonesArray *fMFTClusterArrayFront[fNMaxPlanes];    //! array of front clusters for the planes of the MFT
  87   TClonesArray *fMFTClusterArrayBack[fNMaxPlanes];     //! array of back clusters for the planes of the MFT
  88
  89   TClonesArray *fCandidateTracks;   //! array of candidate global tracks
  90
  91   AliMUONTrack *fMUONTrack;  //! muon track being analyzed
  92
  93   AliMuonForwardTrack *fCurrentTrack;     //! muon extrapolated track being tested
  94   AliMuonForwardTrack *fFinalBestCandidate;     //! best final candidate (if any)
  95
  96   // uncertainty on the x position of the primary vertex (seed for the extrapolation of the MUON tracks)
  97   Double_t fVertexErrorX;
  98   Double_t fVertexErrorY;
  99   Double_t fVertexErrorZ;
 100
 101   Bool_t fBransonCorrection;    // if TRUE, Branson Correction is applied when extrapolating the MUON tracks to the vertex region
 102
 103   Double_t fMinResearchRadiusAtPlane[fNMaxPlanes];
 104
 105 private:
 106
 107   AliMFTTrackerMU(const AliMFTTrackerMU &tracker);
 108   AliMFTTrackerMU & operator=(const AliMFTTrackerMU &tracker);
 109   ClassDef(AliMFTTrackerMU,1)   //MFT tracker for muon tracks
 110
 111 };
 112
 113 //====================================================================================================================================================
 114
 115 #endif
 116