MUON/AliMUONRecoParam.h

   1 #ifndef AliMUONRecoParam_H
   2 #define AliMUONRecoParam_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4 * See cxx source for full Copyright notice                               */
   5
   6 /// \ingroup rec
   7 /// \class AliMUONRecoParam
   8 /// \brief Class with MUON reconstruction parameters
   9 ///
  10 //  Author: Philippe Pillot
  11
  12 #include "AliDetectorRecoParam.h"
  13 #include "TString.h"
  14
  15 class AliMUONRecoParam : public AliDetectorRecoParam
  16 {
  17  public:
  18   AliMUONRecoParam();
  19   virtual ~AliMUONRecoParam();
  20
  21   static AliMUONRecoParam *GetLowFluxParam();
  22   static AliMUONRecoParam *GetHighFluxParam();
  23   static AliMUONRecoParam *GetCosmicParam();
  24
  25   /// set the calibration mode (see GetCalibrationMode() for possible modes)
  26   void SetCalibrationMode(Option_t* mode) { fCalibrationMode = mode; fCalibrationMode.ToUpper();}
  27
  28   Option_t* GetCalibrationMode() const;
  29
  30   /// set the clustering (pre-clustering) mode
  31   void      SetClusteringMode(Option_t* mode) {fClusteringMode = mode; fClusteringMode.ToUpper();}
  32   /// get the clustering (pre-clustering) mode
  33   Option_t* GetClusteringMode() const {return fClusteringMode.Data();}
  34
  35   /// set the tracking mode
  36   void      SetTrackingMode(Option_t* mode) {fTrackingMode = mode; fTrackingMode.ToUpper();}
  37   /// get the tracking mode
  38   Option_t* GetTrackingMode() const {return fTrackingMode.Data();}
  39
  40   /// switch on/off the combined cluster/track reconstruction
  41   void      CombineClusterTrackReco(Bool_t flag) {fCombinedClusterTrackReco = flag;}
  42   /// return kTRUE/kFALSE if the combined cluster/track reconstruction is on/off
  43   Bool_t    CombineClusterTrackReco() const {return fCombinedClusterTrackReco;}
  44
  45   /// save all cluster info (including pads) in ESD, for the given percentage of events
  46   void      SaveFullClusterInESD(Bool_t flag, Double_t percentOfEvent = 100.) {fSaveFullClusterInESD = flag;
  47                                  fPercentOfFullClusterInESD = (fSaveFullClusterInESD) ? percentOfEvent : 0.;}
  48   /// return kTRUE/kFALSE depending on whether we save all cluster info in ESD or not
  49   Bool_t    SaveFullClusterInESD() const {return fSaveFullClusterInESD;}
  50   /// return the percentage of events for which all cluster info are stored in ESD
  51   Double_t  GetPercentOfFullClusterInESD() const {return fPercentOfFullClusterInESD;}
  52
  53   /// set the most probable value (GeV/c) of momentum in bending plane
  54   /// needed to get some "reasonable" corrections for MCS and E loss even if B = 0
  55   void     SetMostProbBendingMomentum(Double_t val) {fMostProbBendingMomentum = val;}
  56   /// return the most probable value (GeV/c) of momentum in bending plane
  57   Double_t GetMostProbBendingMomentum() const {return fMostProbBendingMomentum;}
  58
  59   /// set the minimum value (GeV/c) of momentum in bending plane
  60   void     SetMinBendingMomentum(Double_t val) {fMinBendingMomentum = val;}
  61   /// return the minimum value (GeV/c) of momentum in bending plane
  62   Double_t GetMinBendingMomentum() const {return fMinBendingMomentum;}
  63   /// set the maximum value (GeV/c) of momentum in bending plane
  64   void     SetMaxBendingMomentum(Double_t val) {fMaxBendingMomentum = val;}
  65   /// return the maximum value (GeV/c) of momentum in bending plane
  66   Double_t GetMaxBendingMomentum() const {return fMaxBendingMomentum;}
  67   /// set the maximum value of the non bending slope
  68   void     SetMaxNonBendingSlope(Double_t val) {fMaxNonBendingSlope = val;}
  69   /// return the maximum value of the non bending slope
  70   Double_t GetMaxNonBendingSlope() const {return fMaxNonBendingSlope;}
  71   /// set the maximum value of the bending slope
  72   void     SetMaxBendingSlope(Double_t val) {fMaxBendingSlope = val;}
  73   /// return the maximum value of the bending slope
  74   Double_t GetMaxBendingSlope() const {return fMaxBendingSlope;}
  75
  76   /// set the vertex dispersion (cm) in non bending plane (used for original tracking only)
  77   void     SetNonBendingVertexDispersion(Double_t val) {fNonBendingVertexDispersion = val;}
  78   /// return the vertex dispersion (cm) in bending plane (used for original tracking only)
  79   Double_t GetNonBendingVertexDispersion() const {return fNonBendingVertexDispersion;}
  80   /// set the vertex dispersion (cm) in non bending plane (used for original tracking only)
  81   void     SetBendingVertexDispersion(Double_t val) {fBendingVertexDispersion = val;}
  82   /// return the vertex dispersion (cm) in bending plane (used for original tracking only)
  83   Double_t GetBendingVertexDispersion() const {return fBendingVertexDispersion;}
  84
  85   /// set the maximum distance to the track to search for compatible cluster(s) in non bending direction
  86   void     SetMaxNonBendingDistanceToTrack(Double_t val) {fMaxNonBendingDistanceToTrack = val;}
  87   /// return the maximum distance to the track to search for compatible cluster(s) in non bending direction
  88   Double_t GetMaxNonBendingDistanceToTrack() const {return fMaxNonBendingDistanceToTrack;}
  89   /// set the maximum distance to the track to search for compatible cluster(s) in bending direction
  90   void     SetMaxBendingDistanceToTrack(Double_t val) {fMaxBendingDistanceToTrack = val;}
  91   /// return the maximum distance to the track to search for compatible cluster(s) in bending direction
  92   Double_t GetMaxBendingDistanceToTrack() const {return fMaxBendingDistanceToTrack;}
  93
  94   /// set the cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
  95   void     SetSigmaCutForTracking(Double_t val) {fSigmaCutForTracking = val;}
  96   /// return the cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
  97   Double_t GetSigmaCutForTracking() const {return fSigmaCutForTracking;}
  98
  99   /// switch on/off the track improvement and keep the default cut in sigma to apply on cluster (local chi2)
 100   void     ImproveTracks(Bool_t flag) {fImproveTracks = flag;}
 101   /// switch on/off the track improvement and set the cut in sigma to apply on cluster (local chi2)
 102   void     ImproveTracks(Bool_t flag, Double_t sigmaCut) {fImproveTracks = flag; fSigmaCutForImprovement = sigmaCut;}
 103   /// return kTRUE/kFALSE if the track improvement is switch on/off
 104   Bool_t   ImproveTracks() const {return fImproveTracks;}
 105   /// return the cut in sigma to apply on cluster (local chi2) during track improvement
 106   Double_t GetSigmaCutForImprovement() const {return fSigmaCutForImprovement;}
 107
 108   /// set the cut in sigma to apply on track during trigger hit pattern search
 109   void     SetSigmaCutForTrigger(Double_t val) {fSigmaCutForTrigger = val;}
 110   /// return the cut in sigma to apply on track during trigger hit pattern search
 111   Double_t GetSigmaCutForTrigger() const {return fSigmaCutForTrigger;}
 112   /// set the cut in strips to apply on trigger track during trigger chamber efficiency
 113   void     SetStripCutForTrigger(Double_t val) {fStripCutForTrigger = val;}
 114   /// return the cut in strips to apply on trigger track during trigger chamber efficiency
 115   Double_t GetStripCutForTrigger() const {return fStripCutForTrigger;}
 116   /// set the maximum search area in strips to apply on trigger track during trigger chamber efficiency
 117   void     SetMaxStripAreaForTrigger(Double_t val) {fMaxStripAreaForTrigger = val;}
 118   /// return the maximum search area in strips to apply on trigger track during trigger chamber efficiency
 119   Double_t GetMaxStripAreaForTrigger() const {return fMaxStripAreaForTrigger;}
 120
 121   /// set the maximum normalized chi2 of tracking/trigger track matching
 122   void     SetMaxNormChi2MatchTrigger(Double_t val) {fMaxNormChi2MatchTrigger = val;}
 123   /// return the maximum normalized chi2 of tracking/trigger track matching
 124   Double_t GetMaxNormChi2MatchTrigger() const {return fMaxNormChi2MatchTrigger;}
 125
 126   /// switch on/off the tracking of all the possible candidates (track only the best one if switched off)
 127   void     TrackAllTracks(Bool_t flag) {fTrackAllTracks = flag;}
 128   /// return kTRUE/kFALSE if the tracking of all the possible candidates is switched on/off
 129   Bool_t   TrackAllTracks() const {return fTrackAllTracks;}
 130
 131   /// switch on/off the recovering of tracks being lost during reconstruction
 132   void     RecoverTracks(Bool_t flag) {fRecoverTracks = flag;}
 133   /// return kTRUE/kFALSE if the recovering of tracks being lost during reconstruction is switched on/off
 134   Bool_t   RecoverTracks() const {return fRecoverTracks;}
 135
 136   /// switch on/off the fast building of track candidates (assuming linear propagation between stations 4 and 5)
 137   void     MakeTrackCandidatesFast(Bool_t flag) {fMakeTrackCandidatesFast = flag;}
 138   /// return kTRUE/kFALSE if the fast building of track candidates is switched on/off
 139   Bool_t   MakeTrackCandidatesFast() const {return fMakeTrackCandidatesFast;}
 140
 141   /// switch on/off the building of track candidates starting from 1 cluster in each of the stations 4 and 5
 142   void     MakeMoreTrackCandidates(Bool_t flag) {fMakeMoreTrackCandidates = flag;}
 143   /// return kTRUE/kFALSE if the building of extra track candidates is switched on/off
 144   Bool_t   MakeMoreTrackCandidates() const {return fMakeMoreTrackCandidates;}
 145
 146   /// switch on/off the completion of reconstructed track
 147   void     ComplementTracks(Bool_t flag) {fComplementTracks = flag;}
 148   /// return kTRUE/kFALSE if completion of the reconstructed track is switched on/off
 149   Bool_t   ComplementTracks() const {return fComplementTracks;}
 150
 151   /// switch on/off the use of the smoother
 152   void     UseSmoother(Bool_t flag) {fUseSmoother = flag;}
 153   /// return kTRUE/kFALSE if the use of the smoother is switched on/off
 154   Bool_t   UseSmoother() const {return fUseSmoother;}
 155
 156   /// switch on/off a chamber in the reconstruction
 157   void     UseChamber(Int_t iCh, Bool_t flag) {if (iCh >= 0 && iCh < 10) fUseChamber[iCh] = flag;}
 158   /// return kTRUE/kFALSE whether the chamber must be used or not
 159   Bool_t   UseChamber(Int_t iCh) const {return (iCh >= 0 && iCh < 10) ? fUseChamber[iCh] : kFALSE;}
 160
 161   /// request or not at least one cluster in the station to validate the track
 162   void     RequestStation(Int_t iSt, Bool_t flag) {if (iSt >= 0 && iSt < 5) fRequestStation[iSt] = flag;}
 163   /// return kTRUE/kFALSE whether at least one cluster is requested in the station to validate the track
 164   Bool_t   RequestStation(Int_t iSt) const {return (iSt >= 0 && iSt < 5) ? fRequestStation[iSt] : kFALSE;}
 165
 166   /// set the bypassSt45 value
 167   void BypassSt45(Bool_t value) { fBypassSt45 = value; }
 168   /// return kTRUE if we should replace clusters in St 4 and 5 by generated clusters from trigger tracks
 169   Bool_t BypassSt45() const { return fBypassSt45; }
 170
 171
 172   virtual void Print(Option_t *option = "") const;
 173
 174  private:
 175
 176   /// clustering mode:  NOCLUSTERING, PRECLUSTER, PRECLUSTERV2, PRECLUSTERV3, COG, <pre>
 177   ///                   SIMPLEFIT, SIMPLEFITV3, MLEM:DRAW, MLEM, MLEMV2, MLEMV3   </pre>
 178   TString fClusteringMode; ///< \brief name of the clustering (+ pre-clustering) mode
 179
 180   /// tracking mode: ORIGINAL, KALMAN
 181   TString fTrackingMode; ///< \brief name of the tracking mode
 182
 183   Double32_t fMostProbBendingMomentum; ///< most probable value (GeV/c) of muon momentum in bending plane (used when B = 0)
 184
 185   Double32_t fMinBendingMomentum; ///< minimum value (GeV/c) of momentum in bending plane
 186   Double32_t fMaxBendingMomentum; ///< maximum value (GeV/c) of momentum in bending plane
 187   Double32_t fMaxNonBendingSlope; ///< maximum value of the non bending slope
 188   Double32_t fMaxBendingSlope;    ///< maximum value of the bending slope (used only if B = 0)
 189
 190   Double32_t fNonBendingVertexDispersion; ///< vertex dispersion (cm) in non bending plane (used for original tracking only)
 191   Double32_t fBendingVertexDispersion;    ///< vertex dispersion (cm) in bending plane (used for original tracking only)
 192
 193   Double32_t fMaxNonBendingDistanceToTrack; ///< maximum distance to the track to search for compatible cluster(s) in non bending direction
 194   Double32_t fMaxBendingDistanceToTrack;    ///< maximum distance to the track to search for compatible cluster(s) in bending direction
 195
 196   Double32_t fSigmaCutForTracking; ///< cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
 197
 198   Double32_t fSigmaCutForImprovement; ///< cut in sigma to apply on cluster (local chi2) during track improvement
 199
 200   Double32_t fSigmaCutForTrigger; ///< cut in sigma to apply on track during trigger hit pattern search
 201
 202   Double32_t fStripCutForTrigger; ///< cut in strips to apply on trigger track during trigger chamber efficiency
 203
 204   Double32_t fMaxStripAreaForTrigger; ///< max. search area in strips to apply on trigger track during trigger chamber efficiency
 205
 206   Double32_t fMaxNormChi2MatchTrigger; ///< maximum normalized chi2 of tracking/trigger track matching
 207
 208   Double32_t fPercentOfFullClusterInESD; ///< percentage of events for which all cluster info are stored in ESD
 209
 210   Bool_t     fCombinedClusterTrackReco; ///< switch on/off the combined cluster/track reconstruction
 211
 212   Bool_t     fTrackAllTracks; ///< kTRUE to track all the possible candidates; kFALSE to track only the best ones
 213
 214   Bool_t     fRecoverTracks; ///< kTRUE to try to recover the tracks getting lost during reconstruction
 215
 216   Bool_t     fMakeTrackCandidatesFast; ///< kTRUE to make candidate tracks assuming linear propagation between stations 4 and 5
 217
 218   Bool_t     fMakeMoreTrackCandidates; ///< kTRUE to make candidate tracks starting from 1 cluster in each of the stations 4 and 5
 219
 220   Bool_t     fComplementTracks; ///< kTRUE to try to complete the reconstructed tracks by adding missing clusters
 221
 222   Bool_t     fImproveTracks; ///< kTRUE to try to improve the reconstructed tracks by removing bad clusters
 223
 224   Bool_t     fUseSmoother; ///< kTRUE to use the smoother to compute track parameters/covariances and local chi2 at each cluster (used for Kalman tracking only)
 225
 226   Bool_t     fSaveFullClusterInESD; ///< kTRUE to save all cluster info (including pads) in ESD
 227
 228   /// calibration mode:  GAIN, NOGAIN, GAINCONSTANTCAPA
 229   TString fCalibrationMode; ///<\brief calibration mode
 230
 231   Bool_t fBypassSt45; ///< kTRUE to use trigger tracks to generate "fake" clusters in St 4 and 5
 232
 233   Bool_t     fUseChamber[10]; ///< kTRUE to use the chamber i in the tracking algorithm
 234
 235   Bool_t     fRequestStation[5]; ///< kTRUE to request at least one cluster in station i to validate the track
 236
 237   // functions
 238   void SetLowFluxParam();
 239   void SetHighFluxParam();
 240   void SetCosmicParam();
 241
 242
 243   ClassDef(AliMUONRecoParam,6) // MUON reco parameters
 244 };
 245
 246 #endif
 247