]> git.uio.no Git - u/mrichter/AliRoot.git/blame - MUON/AliMUONRecoParam.h
Added a commented line for disabling QA
[u/mrichter/AliRoot.git] / MUON / AliMUONRecoParam.h
CommitLineData
3304fa09 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
15d30ed4 12#include "AliDetectorRecoParam.h"
3304fa09 13#include "TString.h"
14
15d30ed4 15class AliMUONRecoParam : public AliDetectorRecoParam
3304fa09 16{
17 public:
18 AliMUONRecoParam();
19 virtual ~AliMUONRecoParam();
20
21 static AliMUONRecoParam *GetLowFluxParam();
22 static AliMUONRecoParam *GetHighFluxParam();
0e894e58 23 static AliMUONRecoParam *GetCosmicParam();
3304fa09 24
9bf6860b 25 /// set the calibration mode (see GetCalibrationMode() for possible modes)
7332f213 26 void SetCalibrationMode(Option_t* mode) { fCalibrationMode = mode; fCalibrationMode.ToUpper();}
9bf6860b 27
28 Option_t* GetCalibrationMode() const;
29
3304fa09 30 /// set the clustering (pre-clustering) mode
7332f213 31 void SetClusteringMode(Option_t* mode) {fClusteringMode = mode; fClusteringMode.ToUpper();}
3304fa09 32 /// get the clustering (pre-clustering) mode
33 Option_t* GetClusteringMode() const {return fClusteringMode.Data();}
34
35 /// set the tracking mode
7332f213 36 void SetTrackingMode(Option_t* mode) {fTrackingMode = mode; fTrackingMode.ToUpper();}
3304fa09 37 /// get the tracking mode
38 Option_t* GetTrackingMode() const {return fTrackingMode.Data();}
39
38bcf0ef 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
0a18ba02 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
0e894e58 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
3304fa09 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;}
9bf6860b 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;}
9f093251 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;}
3304fa09 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;}
fda59e58 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;}
b5270f21 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;}
3304fa09 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;}
38bcf0ef 128 /// return kTRUE/kFALSE if the tracking of all the possible candidates is switched on/off
3304fa09 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;}
38bcf0ef 133 /// return kTRUE/kFALSE if the recovering of tracks being lost during reconstruction is switched on/off
3304fa09 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;}
38bcf0ef 138 /// return kTRUE/kFALSE if the fast building of track candidates is switched on/off
3304fa09 139 Bool_t MakeTrackCandidatesFast() const {return fMakeTrackCandidatesFast;}
140
9bf6860b 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
3304fa09 146 /// switch on/off the completion of reconstructed track
147 void ComplementTracks(Bool_t flag) {fComplementTracks = flag;}
38bcf0ef 148 /// return kTRUE/kFALSE if completion of the reconstructed track is switched on/off
3304fa09 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;}
38bcf0ef 153 /// return kTRUE/kFALSE if the use of the smoother is switched on/off
3304fa09 154 Bool_t UseSmoother() const {return fUseSmoother;}
155
9bf6860b 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
6cac085d 166 /// set the bypassSt45 value
167 void BypassSt45(Bool_t value) { fBypassSt45 = value; }
9bf6860b 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; }
3304fa09 170
9bf6860b 171
172 virtual void Print(Option_t *option = "") const;
3304fa09 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
0e894e58 183 Double32_t fMostProbBendingMomentum; ///< most probable value (GeV/c) of muon momentum in bending plane (used when B = 0)
184
3304fa09 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
9bf6860b 187 Double32_t fMaxNonBendingSlope; ///< maximum value of the non bending slope
9f093251 188 Double32_t fMaxBendingSlope; ///< maximum value of the bending slope (used only if B = 0)
3304fa09 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
fda59e58 197
3304fa09 198 Double32_t fSigmaCutForImprovement; ///< cut in sigma to apply on cluster (local chi2) during track improvement
199
fda59e58 200 Double32_t fSigmaCutForTrigger; ///< cut in sigma to apply on track during trigger hit pattern search
b5270f21 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
fda59e58 205
3304fa09 206 Double32_t fMaxNormChi2MatchTrigger; ///< maximum normalized chi2 of tracking/trigger track matching
207
0a18ba02 208 Double32_t fPercentOfFullClusterInESD; ///< percentage of events for which all cluster info are stored in ESD
209
38bcf0ef 210 Bool_t fCombinedClusterTrackReco; ///< switch on/off the combined cluster/track reconstruction
211
3304fa09 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
9bf6860b 218 Bool_t fMakeMoreTrackCandidates; ///< kTRUE to make candidate tracks starting from 1 cluster in each of the stations 4 and 5
219
3304fa09 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
0a18ba02 226 Bool_t fSaveFullClusterInESD; ///< kTRUE to save all cluster info (including pads) in ESD
3304fa09 227
9bf6860b 228 /// calibration mode: GAIN, NOGAIN, GAINCONSTANTCAPA
de98fdc9 229 TString fCalibrationMode; ///<\brief calibration mode
230
9bf6860b 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
3304fa09 237 // functions
238 void SetLowFluxParam();
239 void SetHighFluxParam();
0e894e58 240 void SetCosmicParam();
3304fa09 241
242
b5270f21 243 ClassDef(AliMUONRecoParam,6) // MUON reco parameters
3304fa09 244};
245
246#endif
247