New option to discard mono-cathod clusters (by changing their resolution) and refit...
[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 // $Id$
7
8 /// \ingroup rec
9 /// \class AliMUONRecoParam
10 /// \brief Class with MUON reconstruction parameters
11 ///
12 //  Author: Philippe Pillot
13
14 #include "AliDetectorRecoParam.h"
15 #include "TString.h"
16 #include <TMath.h>
17
18 class AliMUONRecoParam : public AliDetectorRecoParam
19 {
20  public: 
21   AliMUONRecoParam();
22   virtual ~AliMUONRecoParam();
23   
24   static AliMUONRecoParam *GetLowFluxParam();
25   static AliMUONRecoParam *GetHighFluxParam();
26   static AliMUONRecoParam *GetCosmicParam();
27   static AliMUONRecoParam *GetCalibrationParam();
28   
29   /// set the calibration mode (see GetCalibrationMode() for possible modes)
30   void SetCalibrationMode(Option_t* mode) { fCalibrationMode = mode; fCalibrationMode.ToUpper();}
31   
32   Option_t* GetCalibrationMode() const;
33   
34   /// set the clustering (pre-clustering) mode
35   void      SetClusteringMode(Option_t* mode) {fClusteringMode = mode; fClusteringMode.ToUpper();}
36   /// get the clustering (pre-clustering) mode
37   Option_t* GetClusteringMode() const {return fClusteringMode.Data();}
38   
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      
55   /// set the tracking mode
56   void      SetTrackingMode(Option_t* mode) {fTrackingMode = mode; fTrackingMode.ToUpper();}
57   /// get the tracking mode
58   Option_t* GetTrackingMode() const {return fTrackingMode.Data();}
59   
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   
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   
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;}
81   
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;}
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;}
90   
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   
96   /// set the vertex dispersion (cm) in non bending plane
97   void     SetNonBendingVertexDispersion(Double_t val) {fNonBendingVertexDispersion = val;} 
98   /// return the vertex dispersion (cm) in non bending plane
99   Double_t GetNonBendingVertexDispersion() const {return fNonBendingVertexDispersion;}
100   /// set the vertex dispersion (cm) in bending plane
101   void     SetBendingVertexDispersion(Double_t val) {fBendingVertexDispersion = val;} 
102   /// return the vertex dispersion (cm) in bending plane
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;}
127   
128   /// set the cut in sigma to apply on track during trigger hit pattern search
129   void     SetSigmaCutForTrigger(Double_t val) {fSigmaCutForTrigger = val; fMaxNormChi2MatchTrigger = val*val;} 
130   /// return the cut in sigma to apply on track during trigger hit pattern search
131   Double_t GetSigmaCutForTrigger() const {return fSigmaCutForTrigger;}
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;}
140   
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;} 
146   /// return kTRUE/kFALSE if the tracking of all the possible candidates is switched on/off
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;} 
151   /// return kTRUE/kFALSE if the recovering of tracks being lost during reconstruction is switched on/off
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;} 
156   /// return kTRUE/kFALSE if the fast building of track candidates is switched on/off
157   Bool_t   MakeTrackCandidatesFast() const {return fMakeTrackCandidatesFast;}
158   
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   
164   /// switch on/off the completion of reconstructed track
165   void     ComplementTracks(Bool_t flag) {fComplementTracks = flag;} 
166   /// return kTRUE/kFALSE if completion of the reconstructed track is switched on/off
167   Bool_t   ComplementTracks() const {return fComplementTracks;}
168   
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   
174   /// switch on/off the use of the smoother
175   void     UseSmoother(Bool_t flag) {fUseSmoother = flag;} 
176   /// return kTRUE/kFALSE if the use of the smoother is switched on/off
177   Bool_t   UseSmoother() const {return fUseSmoother;}
178   
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;}
188   /// return an integer where first 5 bits are set to 1 if the corresponding station is requested
189   UInt_t   RequestedStationMask() const;
190   
191   /// set the bypassSt45 value
192   void   BypassSt45(Bool_t st4, Bool_t st5);
193   
194   /// return kTRUE if we should replace clusters in St 4 and 5 by generated clusters from trigger tracks
195   Bool_t BypassSt45() const { return fBypassSt45==45; }
196   
197   /// return kTRUE if we should replace clusters in St 4 by generated clusters from trigger tracks
198   Bool_t BypassSt4() const { return BypassSt45() || fBypassSt45==4 ; }
199   
200   /// return kTRUE if we should replace clusters in St 5 by generated clusters from trigger tracks
201   Bool_t BypassSt5() const { return BypassSt45() || fBypassSt45==5 ; }
202   
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;
207   
208   /// Set Low and High threshold for pedestal mean
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   
215   /// Set Low and High threshold for pedestal sigma 
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; }
247   
248   /// Number of sigma cut we must apply when cutting on adc-ped
249   Double_t ChargeSigmaCut() const { return fChargeSigmaCut; }
250   
251   /// Number of sigma cut we must apply when cutting on adc-ped
252   void ChargeSigmaCut(Double_t value) { fChargeSigmaCut=value; }
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   
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   
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]; }
293   
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
299   virtual void Print(Option_t *option = "") const;
300   
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
313   /// Set the percentage of token lost error we allow
314   virtual void SetTokenLostLimit(Double_t limit) { fTokenLostLimit = limit; }
315
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
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   
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   
340 private:
341   
342   void SetDefaultLimits();
343   
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
355   Double32_t fMaxNonBendingSlope; ///< maximum value of the non bending slope
356   Double32_t fMaxBendingSlope;    ///< maximum value of the bending slope (used only if B = 0)
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
365
366   Double32_t fSigmaCutForImprovement; ///< cut in sigma to apply on cluster (local chi2) during track improvement
367   
368   Double32_t fSigmaCutForTrigger; ///< cut in sigma to apply on track during trigger hit pattern search
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
373   
374   Double32_t fMaxNormChi2MatchTrigger; ///< maximum normalized chi2 of tracking/trigger track matching
375   
376   Double32_t fPercentOfFullClusterInESD; ///< percentage of events for which all cluster info are stored in ESD
377   
378   Bool_t     fCombinedClusterTrackReco; ///< switch on/off the combined cluster/track reconstruction
379   
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   
386   Bool_t     fMakeMoreTrackCandidates; ///< kTRUE to make candidate tracks starting from 1 cluster in each of the stations 4 and 5
387   
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   
394   Bool_t     fSaveFullClusterInESD; ///< kTRUE to save all cluster info (including pads) in ESD
395   
396   /// calibration mode:  GAIN, NOGAIN, GAINCONSTANTCAPA, INJECTIONGAIN
397   TString    fCalibrationMode; ///<\brief calibration mode
398   
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
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   
405   Double32_t fGainA1Limits[2]; ///< Low and High threshold for gain a0 parameter
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
408   Double32_t fHVSt12Limits[2]; ///< DEPRECATED. See fHVLimits
409   Double32_t fHVSt345Limits[2]; ///< DEPRECATED. See fHVLimits
410   Double32_t fPedMeanLimits[2]; ///< Low and High threshold for pedestal mean
411   Double32_t fPedSigmaLimits[2]; ///< Low and High threshold for pedestal sigma
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
419   
420   Bool_t     fRemoveConnectedTracksInSt12; ///< kTRUE to remove tracks sharing cluster in station 1 and 2
421   
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   
425   Bool_t     fSelectTrackOnSlope; ///< select track candidates according to their slope (instead of their impact parameter)
426   
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
430
431   Double32_t fMissingPadFractionLimit; ///< DEPRECATED
432   Double32_t fFractionOfBuspatchOutsideOccupancyLimit; ///< above this limit, we consider we have too many buspatches out of the allowed occupancy range
433
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   
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   
442   Bool_t     fTryRecover; ///< try to recover corrupted raw data
443
444   Double32_t fHVLimit[10]; // HV limit (below which we consider that chamber efficiency is to be considered zero)
445   
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   
450   // functions
451   void SetLowFluxParam();
452   void SetHighFluxParam();
453   void SetCosmicParam();
454   void SetCalibrationParam();
455   
456   ClassDef(AliMUONRecoParam,170) // MUON reco parameters
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...
459 };
460
461 #endif
462