remoe duplicate QA initialisation and do ESD QA for same detectors as RecPoint QA
[u/mrichter/AliRoot.git] / 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   
72   /// set the vertex dispersion (cm) in non bending plane (used for original tracking only)
73   void     SetNonBendingVertexDispersion(Double_t val) {fNonBendingVertexDispersion = val;} 
74   /// return the vertex dispersion (cm) in bending plane (used for original tracking only)
75   Double_t GetNonBendingVertexDispersion() const {return fNonBendingVertexDispersion;}
76   /// set the vertex dispersion (cm) in non bending plane (used for original tracking only)
77   void     SetBendingVertexDispersion(Double_t val) {fBendingVertexDispersion = val;} 
78   /// return the vertex dispersion (cm) in bending plane (used for original tracking only)
79   Double_t GetBendingVertexDispersion() const {return fBendingVertexDispersion;}
80   
81   /// set the maximum distance to the track to search for compatible cluster(s) in non bending direction
82   void     SetMaxNonBendingDistanceToTrack(Double_t val) {fMaxNonBendingDistanceToTrack = val;} 
83   /// return the maximum distance to the track to search for compatible cluster(s) in non bending direction
84   Double_t GetMaxNonBendingDistanceToTrack() const {return fMaxNonBendingDistanceToTrack;}
85   /// set the maximum distance to the track to search for compatible cluster(s) in bending direction
86   void     SetMaxBendingDistanceToTrack(Double_t val) {fMaxBendingDistanceToTrack = val;} 
87   /// return the maximum distance to the track to search for compatible cluster(s) in bending direction
88   Double_t GetMaxBendingDistanceToTrack() const {return fMaxBendingDistanceToTrack;}
89   
90   /// set the cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
91   void     SetSigmaCutForTracking(Double_t val) {fSigmaCutForTracking = val;} 
92   /// return the cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
93   Double_t GetSigmaCutForTracking() const {return fSigmaCutForTracking;}
94   
95   /// switch on/off the track improvement and keep the default cut in sigma to apply on cluster (local chi2)
96   void     ImproveTracks(Bool_t flag) {fImproveTracks = flag;} 
97   /// switch on/off the track improvement and set the cut in sigma to apply on cluster (local chi2)
98   void     ImproveTracks(Bool_t flag, Double_t sigmaCut) {fImproveTracks = flag; fSigmaCutForImprovement = sigmaCut;} 
99   /// return kTRUE/kFALSE if the track improvement is switch on/off
100   Bool_t   ImproveTracks() const {return fImproveTracks;}
101   /// return the cut in sigma to apply on cluster (local chi2) during track improvement
102   Double_t GetSigmaCutForImprovement() const {return fSigmaCutForImprovement;}
103
104   /// set the cut in sigma to apply on track during trigger hit pattern search
105   void     SetSigmaCutForTrigger(Double_t val) {fSigmaCutForTrigger = val;} 
106   /// return the cut in sigma to apply on track during trigger hit pattern search
107   Double_t GetSigmaCutForTrigger() const {return fSigmaCutForTrigger;}
108   
109   /// set the maximum normalized chi2 of tracking/trigger track matching
110   void     SetMaxNormChi2MatchTrigger(Double_t val) {fMaxNormChi2MatchTrigger = val;} 
111   /// return the maximum normalized chi2 of tracking/trigger track matching
112   Double_t GetMaxNormChi2MatchTrigger() const {return fMaxNormChi2MatchTrigger;}
113   
114   /// switch on/off the tracking of all the possible candidates (track only the best one if switched off)
115   void     TrackAllTracks(Bool_t flag) {fTrackAllTracks = flag;} 
116   /// return kTRUE/kFALSE if the tracking of all the possible candidates is switched on/off
117   Bool_t   TrackAllTracks() const {return fTrackAllTracks;}
118   
119   /// switch on/off the recovering of tracks being lost during reconstruction
120   void     RecoverTracks(Bool_t flag) {fRecoverTracks = flag;} 
121   /// return kTRUE/kFALSE if the recovering of tracks being lost during reconstruction is switched on/off
122   Bool_t   RecoverTracks() const {return fRecoverTracks;}
123   
124   /// switch on/off the fast building of track candidates (assuming linear propagation between stations 4 and 5)
125   void     MakeTrackCandidatesFast(Bool_t flag) {fMakeTrackCandidatesFast = flag;} 
126   /// return kTRUE/kFALSE if the fast building of track candidates is switched on/off
127   Bool_t   MakeTrackCandidatesFast() const {return fMakeTrackCandidatesFast;}
128   
129   /// switch on/off the building of track candidates starting from 1 cluster in each of the stations 4 and 5
130   void     MakeMoreTrackCandidates(Bool_t flag) {fMakeMoreTrackCandidates = flag;} 
131   /// return kTRUE/kFALSE if the building of extra track candidates is switched on/off
132   Bool_t   MakeMoreTrackCandidates() const {return fMakeMoreTrackCandidates;}
133   
134   /// switch on/off the completion of reconstructed track
135   void     ComplementTracks(Bool_t flag) {fComplementTracks = flag;} 
136   /// return kTRUE/kFALSE if completion of the reconstructed track is switched on/off
137   Bool_t   ComplementTracks() const {return fComplementTracks;}
138   
139   /// switch on/off the use of the smoother
140   void     UseSmoother(Bool_t flag) {fUseSmoother = flag;} 
141   /// return kTRUE/kFALSE if the use of the smoother is switched on/off
142   Bool_t   UseSmoother() const {return fUseSmoother;}
143   
144   /// switch on/off a chamber in the reconstruction
145   void     UseChamber(Int_t iCh, Bool_t flag) {if (iCh >= 0 && iCh < 10) fUseChamber[iCh] = flag;}
146   /// return kTRUE/kFALSE whether the chamber must be used or not
147   Bool_t   UseChamber(Int_t iCh) const {return (iCh >= 0 && iCh < 10) ? fUseChamber[iCh] : kFALSE;}
148   
149   /// request or not at least one cluster in the station to validate the track
150   void     RequestStation(Int_t iSt, Bool_t flag) {if (iSt >= 0 && iSt < 5) fRequestStation[iSt] = flag;}
151   /// return kTRUE/kFALSE whether at least one cluster is requested in the station to validate the track
152   Bool_t   RequestStation(Int_t iSt) const {return (iSt >= 0 && iSt < 5) ? fRequestStation[iSt] : kFALSE;}
153   
154   /// set the bypassSt45 value
155   void BypassSt45(Bool_t value) { fBypassSt45 = value; }
156   /// return kTRUE if we should replace clusters in St 4 and 5 by generated clusters from trigger tracks
157   Bool_t BypassSt45() const { return fBypassSt45; }
158   
159   
160   virtual void Print(Option_t *option = "") const;
161   
162  private:
163   
164   /// clustering mode:  NOCLUSTERING, PRECLUSTER, PRECLUSTERV2, PRECLUSTERV3, COG, <pre>
165   ///                   SIMPLEFIT, SIMPLEFITV3, MLEM:DRAW, MLEM, MLEMV2, MLEMV3   </pre>
166   TString fClusteringMode; ///< \brief name of the clustering (+ pre-clustering) mode
167   
168   /// tracking mode: ORIGINAL, KALMAN
169   TString fTrackingMode; ///< \brief name of the tracking mode
170   
171   Double32_t fMostProbBendingMomentum; ///< most probable value (GeV/c) of muon momentum in bending plane (used when B = 0)
172   
173   Double32_t fMinBendingMomentum; ///< minimum value (GeV/c) of momentum in bending plane
174   Double32_t fMaxBendingMomentum; ///< maximum value (GeV/c) of momentum in bending plane
175   Double32_t fMaxNonBendingSlope; ///< maximum value of the non bending slope
176   
177   Double32_t fNonBendingVertexDispersion; ///< vertex dispersion (cm) in non bending plane (used for original tracking only)
178   Double32_t fBendingVertexDispersion;    ///< vertex dispersion (cm) in bending plane (used for original tracking only)
179   
180   Double32_t fMaxNonBendingDistanceToTrack; ///< maximum distance to the track to search for compatible cluster(s) in non bending direction
181   Double32_t fMaxBendingDistanceToTrack;    ///< maximum distance to the track to search for compatible cluster(s) in bending direction
182   
183   Double32_t fSigmaCutForTracking; ///< cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
184
185   Double32_t fSigmaCutForImprovement; ///< cut in sigma to apply on cluster (local chi2) during track improvement
186   
187   Double32_t fSigmaCutForTrigger; ///< cut in sigma to apply on track during trigger hit pattern search
188   
189   Double32_t fMaxNormChi2MatchTrigger; ///< maximum normalized chi2 of tracking/trigger track matching
190   
191   Double32_t fPercentOfFullClusterInESD; ///< percentage of events for which all cluster info are stored in ESD
192   
193   Bool_t     fCombinedClusterTrackReco; ///< switch on/off the combined cluster/track reconstruction
194   
195   Bool_t     fTrackAllTracks; ///< kTRUE to track all the possible candidates; kFALSE to track only the best ones
196   
197   Bool_t     fRecoverTracks; ///< kTRUE to try to recover the tracks getting lost during reconstruction
198   
199   Bool_t     fMakeTrackCandidatesFast; ///< kTRUE to make candidate tracks assuming linear propagation between stations 4 and 5
200   
201   Bool_t     fMakeMoreTrackCandidates; ///< kTRUE to make candidate tracks starting from 1 cluster in each of the stations 4 and 5
202   
203   Bool_t     fComplementTracks; ///< kTRUE to try to complete the reconstructed tracks by adding missing clusters
204   
205   Bool_t     fImproveTracks; ///< kTRUE to try to improve the reconstructed tracks by removing bad clusters
206   
207   Bool_t     fUseSmoother; ///< kTRUE to use the smoother to compute track parameters/covariances and local chi2 at each cluster (used for Kalman tracking only)
208   
209   Bool_t     fSaveFullClusterInESD; ///< kTRUE to save all cluster info (including pads) in ESD
210   
211   /// calibration mode:  GAIN, NOGAIN, GAINCONSTANTCAPA
212   TString fCalibrationMode; ///<\brief calibration mode
213   
214   Bool_t fBypassSt45; ///< kTRUE to use trigger tracks to generate "fake" clusters in St 4 and 5
215   
216   Bool_t     fUseChamber[10]; ///< kTRUE to use the chamber i in the tracking algorithm
217   
218   Bool_t     fRequestStation[5]; ///< kTRUE to request at least one cluster in station i to validate the track
219   
220   // functions
221   void SetLowFluxParam();
222   void SetHighFluxParam();
223   void SetCosmicParam();
224   
225   
226   ClassDef(AliMUONRecoParam,4) // MUON reco parameters
227 };
228
229 #endif
230