ALIROOT-5420 Changes for CDH v3
[u/mrichter/AliRoot.git] / MUON / AliMUONRecoParam.h
CommitLineData
7d5d0cc5 1#ifndef ALIMUONRECOPARAM_H
2#define ALIMUONRECOPARAM_H
3304fa09 3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4* See cxx source for full Copyright notice */
5
bf4d93eb 6// $Id$
7
3304fa09 8/// \ingroup rec
9/// \class AliMUONRecoParam
10/// \brief Class with MUON reconstruction parameters
11///
12// Author: Philippe Pillot
13
15d30ed4 14#include "AliDetectorRecoParam.h"
3304fa09 15#include "TString.h"
35be7ed7 16#include <TMath.h>
3304fa09 17
15d30ed4 18class AliMUONRecoParam : public AliDetectorRecoParam
3304fa09 19{
20 public:
21 AliMUONRecoParam();
22 virtual ~AliMUONRecoParam();
23
24 static AliMUONRecoParam *GetLowFluxParam();
25 static AliMUONRecoParam *GetHighFluxParam();
0e894e58 26 static AliMUONRecoParam *GetCosmicParam();
c6e702f8 27 static AliMUONRecoParam *GetCalibrationParam();
3304fa09 28
9bf6860b 29 /// set the calibration mode (see GetCalibrationMode() for possible modes)
7332f213 30 void SetCalibrationMode(Option_t* mode) { fCalibrationMode = mode; fCalibrationMode.ToUpper();}
35be7ed7 31
9bf6860b 32 Option_t* GetCalibrationMode() const;
35be7ed7 33
3304fa09 34 /// set the clustering (pre-clustering) mode
7332f213 35 void SetClusteringMode(Option_t* mode) {fClusteringMode = mode; fClusteringMode.ToUpper();}
3304fa09 36 /// get the clustering (pre-clustering) mode
37 Option_t* GetClusteringMode() const {return fClusteringMode.Data();}
38
110edb51 39 /// Get the (truncated) average of sigmas of pedestal measurements, i.e. noise, of pads
40 Double_t AverageNoisePadCharge() const { return fAverageNoisePadCharge; }
41 /// Set the average of sigmas of pedestal measurements, i.e. noise, of pads
42 void AverageNoisePadCharge(Double_t noise) { fAverageNoisePadCharge = noise; }
43
44 /// Get the lowest charge we allow for pads
45 Double_t LowestPadCharge() const { return fChargeSigmaCut*fAverageNoisePadCharge; }
46
47 /// Get the cut applied to cut on cluster charge (the charge is cut if below fClusterChargeCut*LowestPadCharge())
48 Double_t ClusterChargeCut() const { return fClusterChargeCut; }
49 /// Set the cut applied to cut on cluster charge (the charge is cut if below fClusterChargeCut*LowestPadCharge())
50 void ClusterChargeCut(Double_t n) { fClusterChargeCut=n; }
51
52 /// Get the lowest possible cluster charge
53 Double_t LowestClusterCharge() const { return ClusterChargeCut()*LowestPadCharge(); }
54
3304fa09 55 /// set the tracking mode
7332f213 56 void SetTrackingMode(Option_t* mode) {fTrackingMode = mode; fTrackingMode.ToUpper();}
3304fa09 57 /// get the tracking mode
58 Option_t* GetTrackingMode() const {return fTrackingMode.Data();}
59
38bcf0ef 60 /// switch on/off the combined cluster/track reconstruction
61 void CombineClusterTrackReco(Bool_t flag) {fCombinedClusterTrackReco = flag;}
62 /// return kTRUE/kFALSE if the combined cluster/track reconstruction is on/off
63 Bool_t CombineClusterTrackReco() const {return fCombinedClusterTrackReco;}
64
0a18ba02 65 /// save all cluster info (including pads) in ESD, for the given percentage of events
66 void SaveFullClusterInESD(Bool_t flag, Double_t percentOfEvent = 100.) {fSaveFullClusterInESD = flag;
67 fPercentOfFullClusterInESD = (fSaveFullClusterInESD) ? percentOfEvent : 0.;}
68 /// return kTRUE/kFALSE depending on whether we save all cluster info in ESD or not
69 Bool_t SaveFullClusterInESD() const {return fSaveFullClusterInESD;}
70 /// return the percentage of events for which all cluster info are stored in ESD
71 Double_t GetPercentOfFullClusterInESD() const {return fPercentOfFullClusterInESD;}
72
3304fa09 73 /// set the minimum value (GeV/c) of momentum in bending plane
74 void SetMinBendingMomentum(Double_t val) {fMinBendingMomentum = val;}
75 /// return the minimum value (GeV/c) of momentum in bending plane
76 Double_t GetMinBendingMomentum() const {return fMinBendingMomentum;}
77 /// set the maximum value (GeV/c) of momentum in bending plane
78 void SetMaxBendingMomentum(Double_t val) {fMaxBendingMomentum = val;}
79 /// return the maximum value (GeV/c) of momentum in bending plane
80 Double_t GetMaxBendingMomentum() const {return fMaxBendingMomentum;}
6b191dea 81
9bf6860b 82 /// set the maximum value of the non bending slope
83 void SetMaxNonBendingSlope(Double_t val) {fMaxNonBendingSlope = val;}
84 /// return the maximum value of the non bending slope
85 Double_t GetMaxNonBendingSlope() const {return fMaxNonBendingSlope;}
9f093251 86 /// set the maximum value of the bending slope
87 void SetMaxBendingSlope(Double_t val) {fMaxBendingSlope = val;}
88 /// return the maximum value of the bending slope
89 Double_t GetMaxBendingSlope() const {return fMaxBendingSlope;}
3304fa09 90
6b191dea 91 /// switch on/off the track selection according to their slope (instead of their impact parameter)
92 void SelectOnTrackSlope(Bool_t flag) {fSelectTrackOnSlope = flag;}
93 /// return kTRUE/kFALSE if tracks are selected according to their slope/impact parameter
94 Bool_t SelectOnTrackSlope() const {return fSelectTrackOnSlope;}
95
ad6250dc 96 /// set the vertex dispersion (cm) in non bending plane
3304fa09 97 void SetNonBendingVertexDispersion(Double_t val) {fNonBendingVertexDispersion = val;}
ad6250dc 98 /// return the vertex dispersion (cm) in non bending plane
3304fa09 99 Double_t GetNonBendingVertexDispersion() const {return fNonBendingVertexDispersion;}
ad6250dc 100 /// set the vertex dispersion (cm) in bending plane
3304fa09 101 void SetBendingVertexDispersion(Double_t val) {fBendingVertexDispersion = val;}
ad6250dc 102 /// return the vertex dispersion (cm) in bending plane
3304fa09 103 Double_t GetBendingVertexDispersion() const {return fBendingVertexDispersion;}
104
105 /// set the maximum distance to the track to search for compatible cluster(s) in non bending direction
106 void SetMaxNonBendingDistanceToTrack(Double_t val) {fMaxNonBendingDistanceToTrack = val;}
107 /// return the maximum distance to the track to search for compatible cluster(s) in non bending direction
108 Double_t GetMaxNonBendingDistanceToTrack() const {return fMaxNonBendingDistanceToTrack;}
109 /// set the maximum distance to the track to search for compatible cluster(s) in bending direction
110 void SetMaxBendingDistanceToTrack(Double_t val) {fMaxBendingDistanceToTrack = val;}
111 /// return the maximum distance to the track to search for compatible cluster(s) in bending direction
112 Double_t GetMaxBendingDistanceToTrack() const {return fMaxBendingDistanceToTrack;}
113
114 /// set the cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
115 void SetSigmaCutForTracking(Double_t val) {fSigmaCutForTracking = val;}
116 /// return the cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
117 Double_t GetSigmaCutForTracking() const {return fSigmaCutForTracking;}
118
119 /// switch on/off the track improvement and keep the default cut in sigma to apply on cluster (local chi2)
120 void ImproveTracks(Bool_t flag) {fImproveTracks = flag;}
121 /// switch on/off the track improvement and set the cut in sigma to apply on cluster (local chi2)
122 void ImproveTracks(Bool_t flag, Double_t sigmaCut) {fImproveTracks = flag; fSigmaCutForImprovement = sigmaCut;}
123 /// return kTRUE/kFALSE if the track improvement is switch on/off
124 Bool_t ImproveTracks() const {return fImproveTracks;}
125 /// return the cut in sigma to apply on cluster (local chi2) during track improvement
126 Double_t GetSigmaCutForImprovement() const {return fSigmaCutForImprovement;}
35be7ed7 127
fda59e58 128 /// set the cut in sigma to apply on track during trigger hit pattern search
4dd97b60 129 void SetSigmaCutForTrigger(Double_t val) {fSigmaCutForTrigger = val; fMaxNormChi2MatchTrigger = val*val;}
fda59e58 130 /// return the cut in sigma to apply on track during trigger hit pattern search
131 Double_t GetSigmaCutForTrigger() const {return fSigmaCutForTrigger;}
b5270f21 132 /// set the cut in strips to apply on trigger track during trigger chamber efficiency
133 void SetStripCutForTrigger(Double_t val) {fStripCutForTrigger = val;}
134 /// return the cut in strips to apply on trigger track during trigger chamber efficiency
135 Double_t GetStripCutForTrigger() const {return fStripCutForTrigger;}
136 /// set the maximum search area in strips to apply on trigger track during trigger chamber efficiency
137 void SetMaxStripAreaForTrigger(Double_t val) {fMaxStripAreaForTrigger = val;}
138 /// return the maximum search area in strips to apply on trigger track during trigger chamber efficiency
139 Double_t GetMaxStripAreaForTrigger() const {return fMaxStripAreaForTrigger;}
3304fa09 140
3304fa09 141 /// return the maximum normalized chi2 of tracking/trigger track matching
142 Double_t GetMaxNormChi2MatchTrigger() const {return fMaxNormChi2MatchTrigger;}
143
144 /// switch on/off the tracking of all the possible candidates (track only the best one if switched off)
145 void TrackAllTracks(Bool_t flag) {fTrackAllTracks = flag;}
38bcf0ef 146 /// return kTRUE/kFALSE if the tracking of all the possible candidates is switched on/off
3304fa09 147 Bool_t TrackAllTracks() const {return fTrackAllTracks;}
148
149 /// switch on/off the recovering of tracks being lost during reconstruction
150 void RecoverTracks(Bool_t flag) {fRecoverTracks = flag;}
38bcf0ef 151 /// return kTRUE/kFALSE if the recovering of tracks being lost during reconstruction is switched on/off
3304fa09 152 Bool_t RecoverTracks() const {return fRecoverTracks;}
153
154 /// switch on/off the fast building of track candidates (assuming linear propagation between stations 4 and 5)
155 void MakeTrackCandidatesFast(Bool_t flag) {fMakeTrackCandidatesFast = flag;}
38bcf0ef 156 /// return kTRUE/kFALSE if the fast building of track candidates is switched on/off
3304fa09 157 Bool_t MakeTrackCandidatesFast() const {return fMakeTrackCandidatesFast;}
158
9bf6860b 159 /// switch on/off the building of track candidates starting from 1 cluster in each of the stations 4 and 5
160 void MakeMoreTrackCandidates(Bool_t flag) {fMakeMoreTrackCandidates = flag;}
161 /// return kTRUE/kFALSE if the building of extra track candidates is switched on/off
162 Bool_t MakeMoreTrackCandidates() const {return fMakeMoreTrackCandidates;}
163
3304fa09 164 /// switch on/off the completion of reconstructed track
165 void ComplementTracks(Bool_t flag) {fComplementTracks = flag;}
38bcf0ef 166 /// return kTRUE/kFALSE if completion of the reconstructed track is switched on/off
3304fa09 167 Bool_t ComplementTracks() const {return fComplementTracks;}
168
c59f70b9 169 /// remove tracks sharing cluster in stations 1 or 2
170 void RemoveConnectedTracksInSt12(Bool_t flag) {fRemoveConnectedTracksInSt12 = flag;}
171 /// return kTRUE/kFALSE whether tracks sharing cluster in station 1 and 2 must be removed or not
172 Bool_t RemoveConnectedTracksInSt12() const {return fRemoveConnectedTracksInSt12;}
173
3304fa09 174 /// switch on/off the use of the smoother
175 void UseSmoother(Bool_t flag) {fUseSmoother = flag;}
38bcf0ef 176 /// return kTRUE/kFALSE if the use of the smoother is switched on/off
3304fa09 177 Bool_t UseSmoother() const {return fUseSmoother;}
178
9bf6860b 179 /// switch on/off a chamber in the reconstruction
180 void UseChamber(Int_t iCh, Bool_t flag) {if (iCh >= 0 && iCh < 10) fUseChamber[iCh] = flag;}
181 /// return kTRUE/kFALSE whether the chamber must be used or not
182 Bool_t UseChamber(Int_t iCh) const {return (iCh >= 0 && iCh < 10) ? fUseChamber[iCh] : kFALSE;}
183
184 /// request or not at least one cluster in the station to validate the track
185 void RequestStation(Int_t iSt, Bool_t flag) {if (iSt >= 0 && iSt < 5) fRequestStation[iSt] = flag;}
186 /// return kTRUE/kFALSE whether at least one cluster is requested in the station to validate the track
187 Bool_t RequestStation(Int_t iSt) const {return (iSt >= 0 && iSt < 5) ? fRequestStation[iSt] : kFALSE;}
9bdbee64 188 /// return an integer where first 5 bits are set to 1 if the corresponding station is requested
189 UInt_t RequestedStationMask() const;
9bf6860b 190
6cac085d 191 /// set the bypassSt45 value
35be7ed7 192 void BypassSt45(Bool_t st4, Bool_t st5);
193
9bf6860b 194 /// return kTRUE if we should replace clusters in St 4 and 5 by generated clusters from trigger tracks
004a9ccd 195 Bool_t BypassSt45() const { return fBypassSt45==45; }
9bf6860b 196
35be7ed7 197 /// return kTRUE if we should replace clusters in St 4 by generated clusters from trigger tracks
004a9ccd 198 Bool_t BypassSt4() const { return BypassSt45() || fBypassSt45==4 ; }
35be7ed7 199
200 /// return kTRUE if we should replace clusters in St 5 by generated clusters from trigger tracks
004a9ccd 201 Bool_t BypassSt5() const { return BypassSt45() || fBypassSt45==5 ; }
202
32f1b761 203 /// Set HV threshold for chambers (chamberId=0..9, use -1 to set all chambers equal)
204 void SetHVLimit(Int_t chamberId, Double_t ichamber);
205 /// Retrieve HV limit for chamber (chamberId=0..9)
206 Double_t HVLimit(Int_t chamberId) const;
35be7ed7 207
004a9ccd 208 /// Set Low and High threshold for pedestal mean
35be7ed7 209 void SetPedMeanLimits(float low, float high) { fPedMeanLimits[0]=low; fPedMeanLimits[1]=high; }
210 /// Retrieve low limit of ped mean
211 Float_t PedMeanLowLimit() const { return fPedMeanLimits[0]; }
212 /// Retrieve high limit of ped mean
213 Float_t PedMeanHighLimit() const { return fPedMeanLimits[1]; }
214
004a9ccd 215 /// Set Low and High threshold for pedestal sigma
35be7ed7 216 void SetPedSigmaLimits(float low, float high) { fPedSigmaLimits[0]=low; fPedSigmaLimits[1]=high; }
217 /// Retrieve low limit of ped sigma
218 Float_t PedSigmaLowLimit() const { return fPedSigmaLimits[0]; }
219 /// Retrieve high limit of ped sigma
220 Float_t PedSigmaHighLimit() const { return fPedSigmaLimits[1]; }
221
222 /// Set Low and High threshold for gain a0 term
223 void SetGainA1Limits(float low, float high) { fGainA1Limits[0]=low; fGainA1Limits[1]=high; }
224 /// Retrieve low limit of a1 (linear term) gain parameter
225 Float_t GainA1LowLimit() const { return fGainA1Limits[0]; }
226 /// Retrieve high limit of a1 (linear term) gain parameter
227 Float_t GainA1HighLimit() const { return fGainA1Limits[1]; }
228
229 /// Set Low and High threshold for gain a1 term
230 void SetGainA2Limits(float low, float high) { fGainA2Limits[0]=low; fGainA2Limits[1]=high; }
231 /// Retrieve low limit of a2 (quadratic term) gain parameter
232 Float_t GainA2LowLimit() const { return fGainA2Limits[0]; }
233 /// Retrieve high limit of a2 (quadratic term) gain parameter
234 Float_t GainA2HighLimit() const { return fGainA2Limits[1]; }
235
236 /// Set Low and High threshold for gain threshold term
237 void SetGainThresLimits(float low, float high) { fGainThresLimits[0]=low; fGainThresLimits[1]=high; }
238 /// Retrieve low limit on threshold gain parameter
239 Float_t GainThresLowLimit() const { return fGainThresLimits[0]; }
240 /// Retrieve high limit on threshold gain parameter
241 Float_t GainThresHighLimit() const { return fGainThresLimits[1]; }
242
243 /// Set the goodness mask (see AliMUONPadStatusMapMaker)
244 void SetPadGoodnessMask(UInt_t mask) { fPadGoodnessMask=mask; }
245 /// Get the goodness mask
246 UInt_t PadGoodnessMask() const { return fPadGoodnessMask; }
3304fa09 247
170f4046 248 /// Number of sigma cut we must apply when cutting on adc-ped
249 Double_t ChargeSigmaCut() const { return fChargeSigmaCut; }
35be7ed7 250
170f4046 251 /// Number of sigma cut we must apply when cutting on adc-ped
252 void ChargeSigmaCut(Double_t value) { fChargeSigmaCut=value; }
35be7ed7 253
254 /// Set the default non bending resolution of chamber iCh
255 void SetDefaultNonBendingReso(Int_t iCh, Double_t val) {if (iCh >= 0 && iCh < 10) fDefaultNonBendingReso[iCh] = val;}
256 /// Get the default non bending resolution of chamber iCh
257 Double_t GetDefaultNonBendingReso(Int_t iCh) const {return (iCh >= 0 && iCh < 10) ? fDefaultNonBendingReso[iCh] : FLT_MAX;}
258 /// Set the default bending resolution of chamber iCh
259 void SetDefaultBendingReso(Int_t iCh, Double_t val) {if (iCh >= 0 && iCh < 10) fDefaultBendingReso[iCh] = val;}
260 /// Get the default bending resolution of chamber iCh
261 Double_t GetDefaultBendingReso(Int_t iCh) const {return (iCh >= 0 && iCh < 10) ? fDefaultBendingReso[iCh] : FLT_MAX;}
262
89c8d66d 263 /// Set the maximum number of trigger tracks above which the tracking is cancelled
264 void SetMaxTriggerTracks(Int_t maxTriggerTracks) {fMaxTriggerTracks = maxTriggerTracks;}
265 /// Get the maximum number of trigger tracks above which the tracking is cancelled
266 Int_t GetMaxTriggerTracks() const {return fMaxTriggerTracks;}
267
268 /// Set the maximum number of track candidates above which the tracking abort
269 void SetMaxTrackCandidates(Int_t maxTrackCandidates) {fMaxTrackCandidates = maxTrackCandidates;}
270 /// Get the maximum number of track candidates above which the tracking abort
271 Int_t GetMaxTrackCandidates() const {return fMaxTrackCandidates;}
272
7eafe398 273 /// Set the limits for the acceptable manu occupancy
274 void SetManuOccupancyLimits(float low, float high) { fManuOccupancyLimits[0]=low; fManuOccupancyLimits[1]=high; }
275 /// Retrieve low value of manu occupancy limit
276 Float_t ManuOccupancyLowLimit() const { return fManuOccupancyLimits[0]; }
277 /// Retrieve high value of manu occupancy limit
278 Float_t ManuOccupancyHighLimit() const { return fManuOccupancyLimits[1]; }
279
280 /// Set the limits for the acceptable bp occupancy
281 void SetBuspatchOccupancyLimits(float low, float high) { fBuspatchOccupancyLimits[0]=low; fBuspatchOccupancyLimits[1]=high; }
282 /// Retrieve low value of bp occupancy limit
283 Float_t BuspatchOccupancyLowLimit() const { return fBuspatchOccupancyLimits[0]; }
284 /// Retrieve high value of bp occupancy limit
285 Float_t BuspatchOccupancyHighLimit() const { return fBuspatchOccupancyLimits[1]; }
286
287 /// Set the limits for the acceptable DE occupancy
288 void SetDEOccupancyLimits(float low, float high) { fDEOccupancyLimits[0]=low; fDEOccupancyLimits[1]=high; }
289 /// Retrieve low value of DE occupancy limit
290 Float_t DEOccupancyLowLimit() const { return fDEOccupancyLimits[0]; }
291 /// Retrieve high value of DE occupancy limit
292 Float_t DEOccupancyHighLimit() const { return fDEOccupancyLimits[1]; }
35be7ed7 293
ece56eb9 294 /// Set the fraction of buspatches outside the occupancy limits
295 void SetFractionOfBuspatchOutsideOccupancyLimit(float v) { fFractionOfBuspatchOutsideOccupancyLimit = v; }
296 /// Get the fraction of buspatches outside the occupancy limits
297 Float_t FractionOfBuspatchOutsideOccupancyLimit() const { return fFractionOfBuspatchOutsideOccupancyLimit; }
298
7eafe398 299 virtual void Print(Option_t *option = "") const;
35be7ed7 300
b3d57767 301 /// Get the max event size (soft limit)
302 virtual Double_t EventSizeSoftLimit() const { return fEventSizeSoftLimit; }
303
304 /// Get the max event size (hard limit)
305 virtual Double_t EventSizeHardLimit() const { return fEventSizeHardLimit; }
306
307 /// Set the max event size limits
308 virtual void SetEventSizeLimits(Double_t soft, Double_t hard) { fEventSizeSoftLimit=soft; fEventSizeHardLimit=hard; }
309
310 /// Get the percentage of token lost error we allow
311 virtual Double_t TokenLostLimit() const { return fTokenLostLimit; }
312
9074a9a9 313 /// Set the percentage of token lost error we allow
b3d57767 314 virtual void SetTokenLostLimit(Double_t limit) { fTokenLostLimit = limit; }
315
9074a9a9 316 /// Whether or not we try to recover corrupted raw data
317 virtual Bool_t TryRecover() const { return fTryRecover; }
318
319 /// Set the try recover corrupted raw data (use kTRUE only if you know what you are doing. Should be left to kFALSE by default)
320 virtual void TryRecover(Bool_t flag) { fTryRecover = flag; }
321
852fefb4 322 /// Discard or not the mono-cathod clusters by assigning to them different resolutions (use default values)
323 void DiscardMonoCathodClusters(Bool_t flag) { fDiscardMonoCathodClusters = flag; }
324 /// Discard or not the mono-cathod clusters by assigning to them different resolutions (use given values)
325 void DiscardMonoCathodClusters(Bool_t flag, Double_t resNB, Double_t resB) { fDiscardMonoCathodClusters = flag;
326 fMonoCathodClNonBendingRes = resNB; fMonoCathodClBendingRes = resB; }
327 /// Check whether to discard or not the mono-cathod clusters
328 Bool_t DiscardMonoCathodClusters() const { return fDiscardMonoCathodClusters; }
329 /// Get the non-bending resolution of mono-cathod clusters when the non-bending plane is missing
330 Double_t GetMonoCathodClNonBendingRes() const { return fMonoCathodClNonBendingRes; }
331 /// Get the bending resolution of mono-cathod clusters when the bending plane is missing
332 Double_t GetMonoCathodClBendingRes() const { return fMonoCathodClBendingRes; }
333
9074a9a9 334 /// Create object ready to be put in OCDB
335 static TObjArray* Create(const char* settings);
336
337 /// Show what is the OCDB for that run
338 static void Show(Int_t runNumber, const char* ocdbPath="raw://");
339
b3d57767 340private:
35be7ed7 341
342 void SetDefaultLimits();
343
3304fa09 344 private:
345
346 /// clustering mode: NOCLUSTERING, PRECLUSTER, PRECLUSTERV2, PRECLUSTERV3, COG, <pre>
347 /// SIMPLEFIT, SIMPLEFITV3, MLEM:DRAW, MLEM, MLEMV2, MLEMV3 </pre>
348 TString fClusteringMode; ///< \brief name of the clustering (+ pre-clustering) mode
349
350 /// tracking mode: ORIGINAL, KALMAN
351 TString fTrackingMode; ///< \brief name of the tracking mode
352
353 Double32_t fMinBendingMomentum; ///< minimum value (GeV/c) of momentum in bending plane
354 Double32_t fMaxBendingMomentum; ///< maximum value (GeV/c) of momentum in bending plane
9bf6860b 355 Double32_t fMaxNonBendingSlope; ///< maximum value of the non bending slope
9f093251 356 Double32_t fMaxBendingSlope; ///< maximum value of the bending slope (used only if B = 0)
3304fa09 357
358 Double32_t fNonBendingVertexDispersion; ///< vertex dispersion (cm) in non bending plane (used for original tracking only)
359 Double32_t fBendingVertexDispersion; ///< vertex dispersion (cm) in bending plane (used for original tracking only)
360
361 Double32_t fMaxNonBendingDistanceToTrack; ///< maximum distance to the track to search for compatible cluster(s) in non bending direction
362 Double32_t fMaxBendingDistanceToTrack; ///< maximum distance to the track to search for compatible cluster(s) in bending direction
363
364 Double32_t fSigmaCutForTracking; ///< cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
fda59e58 365
3304fa09 366 Double32_t fSigmaCutForImprovement; ///< cut in sigma to apply on cluster (local chi2) during track improvement
367
fda59e58 368 Double32_t fSigmaCutForTrigger; ///< cut in sigma to apply on track during trigger hit pattern search
b5270f21 369
370 Double32_t fStripCutForTrigger; ///< cut in strips to apply on trigger track during trigger chamber efficiency
371
372 Double32_t fMaxStripAreaForTrigger; ///< max. search area in strips to apply on trigger track during trigger chamber efficiency
fda59e58 373
3304fa09 374 Double32_t fMaxNormChi2MatchTrigger; ///< maximum normalized chi2 of tracking/trigger track matching
375
0a18ba02 376 Double32_t fPercentOfFullClusterInESD; ///< percentage of events for which all cluster info are stored in ESD
377
38bcf0ef 378 Bool_t fCombinedClusterTrackReco; ///< switch on/off the combined cluster/track reconstruction
379
3304fa09 380 Bool_t fTrackAllTracks; ///< kTRUE to track all the possible candidates; kFALSE to track only the best ones
381
382 Bool_t fRecoverTracks; ///< kTRUE to try to recover the tracks getting lost during reconstruction
383
384 Bool_t fMakeTrackCandidatesFast; ///< kTRUE to make candidate tracks assuming linear propagation between stations 4 and 5
385
9bf6860b 386 Bool_t fMakeMoreTrackCandidates; ///< kTRUE to make candidate tracks starting from 1 cluster in each of the stations 4 and 5
387
3304fa09 388 Bool_t fComplementTracks; ///< kTRUE to try to complete the reconstructed tracks by adding missing clusters
389
390 Bool_t fImproveTracks; ///< kTRUE to try to improve the reconstructed tracks by removing bad clusters
391
392 Bool_t fUseSmoother; ///< kTRUE to use the smoother to compute track parameters/covariances and local chi2 at each cluster (used for Kalman tracking only)
393
0a18ba02 394 Bool_t fSaveFullClusterInESD; ///< kTRUE to save all cluster info (including pads) in ESD
3304fa09 395
66cdf5b3 396 /// calibration mode: GAIN, NOGAIN, GAINCONSTANTCAPA, INJECTIONGAIN
35be7ed7 397 TString fCalibrationMode; ///<\brief calibration mode
de98fdc9 398
35be7ed7 399 Int_t fBypassSt45; ///< non-zero to use trigger tracks to generate "fake" clusters in St 4 and 5. Can be 0, 4, 5 or 45 only
9bf6860b 400
401 Bool_t fUseChamber[10]; ///< kTRUE to use the chamber i in the tracking algorithm
402
403 Bool_t fRequestStation[5]; ///< kTRUE to request at least one cluster in station i to validate the track
404
35be7ed7 405 Double32_t fGainA1Limits[2]; ///< Low and High threshold for gain a0 parameter
113ad708 406 Double32_t fGainA2Limits[2]; ///< Low and High threshold for gain a1 parameter
407 Double32_t fGainThresLimits[2]; ///< Low and High threshold for gain threshold parameter
32f1b761 408 Double32_t fHVSt12Limits[2]; ///< DEPRECATED. See fHVLimits
409 Double32_t fHVSt345Limits[2]; ///< DEPRECATED. See fHVLimits
113ad708 410 Double32_t fPedMeanLimits[2]; ///< Low and High threshold for pedestal mean
411 Double32_t fPedSigmaLimits[2]; ///< Low and High threshold for pedestal sigma
35be7ed7 412
413 UInt_t fPadGoodnessMask; ///< goodness mask (see AliMUONPadStatusMaker)
414
415 Double32_t fChargeSigmaCut; ///< number of sigma to cut on adc-ped
416
417 Double32_t fDefaultNonBendingReso[10]; ///< default chamber resolution in the non-bending direction
418 Double32_t fDefaultBendingReso[10]; ///< default chamber resolution in the bending direction
170f4046 419
c59f70b9 420 Bool_t fRemoveConnectedTracksInSt12; ///< kTRUE to remove tracks sharing cluster in station 1 and 2
421
89c8d66d 422 Int_t fMaxTriggerTracks; ///< maximum number of trigger tracks above which the tracking is cancelled
423 Int_t fMaxTrackCandidates; ///< maximum number of track candidates above which the tracking abort
424
6b191dea 425 Bool_t fSelectTrackOnSlope; ///< select track candidates according to their slope (instead of their impact parameter)
426
7eafe398 427 Double32_t fManuOccupancyLimits[2]; ///< low and high thresholds for manu occupancy cut
428 Double32_t fBuspatchOccupancyLimits[2]; ///< low and high thresholds for bus patch occupancy cut
429 Double32_t fDEOccupancyLimits[2]; ///< low and high thresholds for DE occupancy cut
ece56eb9 430
b3d57767 431 Double32_t fMissingPadFractionLimit; ///< DEPRECATED
ece56eb9 432 Double32_t fFractionOfBuspatchOutsideOccupancyLimit; ///< above this limit, we consider we have too many buspatches out of the allowed occupancy range
433
110edb51 434 Double32_t fAverageNoisePadCharge; ///< the (truncated, typically at 10%) mean of the sigma of the pedestals, in femto-coulomb
435 Double32_t fClusterChargeCut; ///< the cluster is cut if its charge is below fClusterChargeCut*LowestPadCharge()
436
b3d57767 437 Double32_t fEventSizeSoftLimit; ///< (soft) limit on mean event size per event (KB)
438 Double32_t fEventSizeHardLimit; ///< (hard) limit on mean event size per event (KB)
439
440 Double32_t fTokenLostLimit; ///< limit on the fraction of token lost error per event we allow
441
9074a9a9 442 Bool_t fTryRecover; ///< try to recover corrupted raw data
32f1b761 443
444 Double32_t fHVLimit[10]; // HV limit (below which we consider that chamber efficiency is to be considered zero)
9074a9a9 445
852fefb4 446 Double32_t fDiscardMonoCathodClusters; // assign a different resolution to mono-cathod clusters in the direction of the missing plane
447 Double32_t fMonoCathodClNonBendingRes; // resolution of mono-cathod clusters in the non-bending direction when the non-bending plane is missing
448 Double32_t fMonoCathodClBendingRes; // resolution of mono-cathod clusters in the bending direction when the bending plane is missing
449
3304fa09 450 // functions
451 void SetLowFluxParam();
452 void SetHighFluxParam();
0e894e58 453 void SetCosmicParam();
c6e702f8 454 void SetCalibrationParam();
35be7ed7 455
852fefb4 456 ClassDef(AliMUONRecoParam,170) // MUON reco parameters
b3d57767 457 // we're at 167 not because we had that many versions, but because at some point (version 15->16)
458 // 166 was committed by error, and we did not to go reverse afterwards...
3304fa09 459};
460
461#endif
462