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 */
7 /// \class AliMUONRecoParam
8 /// \brief Class with MUON reconstruction parameters
10 // Author: Philippe Pillot
12 #include "AliDetectorRecoParam.h"
15 class AliMUONRecoParam : public AliDetectorRecoParam
19 virtual ~AliMUONRecoParam();
21 static AliMUONRecoParam *GetLowFluxParam();
22 static AliMUONRecoParam *GetHighFluxParam();
24 /// set the calibration mode (see GetCalibrationMode() for possible modes)
25 void SetCalibrationMode(Option_t* mode) { fCalibrationMode = mode; fCalibrationMode.ToUpper();}
27 Option_t* GetCalibrationMode() const;
29 /// set the clustering (pre-clustering) mode
30 void SetClusteringMode(Option_t* mode) {fClusteringMode = mode; fClusteringMode.ToUpper();}
31 /// get the clustering (pre-clustering) mode
32 Option_t* GetClusteringMode() const {return fClusteringMode.Data();}
34 /// set the tracking mode
35 void SetTrackingMode(Option_t* mode) {fTrackingMode = mode; fTrackingMode.ToUpper();}
36 /// get the tracking mode
37 Option_t* GetTrackingMode() const {return fTrackingMode.Data();}
39 /// switch on/off the combined cluster/track reconstruction
40 void CombineClusterTrackReco(Bool_t flag) {fCombinedClusterTrackReco = flag;}
41 /// return kTRUE/kFALSE if the combined cluster/track reconstruction is on/off
42 Bool_t CombineClusterTrackReco() const {return fCombinedClusterTrackReco;}
44 /// save all cluster info (including pads) in ESD, for the given percentage of events
45 void SaveFullClusterInESD(Bool_t flag, Double_t percentOfEvent = 100.) {fSaveFullClusterInESD = flag;
46 fPercentOfFullClusterInESD = (fSaveFullClusterInESD) ? percentOfEvent : 0.;}
47 /// return kTRUE/kFALSE depending on whether we save all cluster info in ESD or not
48 Bool_t SaveFullClusterInESD() const {return fSaveFullClusterInESD;}
49 /// return the percentage of events for which all cluster info are stored in ESD
50 Double_t GetPercentOfFullClusterInESD() const {return fPercentOfFullClusterInESD;}
52 /// set the minimum value (GeV/c) of momentum in bending plane
53 void SetMinBendingMomentum(Double_t val) {fMinBendingMomentum = val;}
54 /// return the minimum value (GeV/c) of momentum in bending plane
55 Double_t GetMinBendingMomentum() const {return fMinBendingMomentum;}
56 /// set the maximum value (GeV/c) of momentum in bending plane
57 void SetMaxBendingMomentum(Double_t val) {fMaxBendingMomentum = val;}
58 /// return the maximum value (GeV/c) of momentum in bending plane
59 Double_t GetMaxBendingMomentum() const {return fMaxBendingMomentum;}
60 /// set the maximum value of the non bending slope
61 void SetMaxNonBendingSlope(Double_t val) {fMaxNonBendingSlope = val;}
62 /// return the maximum value of the non bending slope
63 Double_t GetMaxNonBendingSlope() const {return fMaxNonBendingSlope;}
65 /// set the vertex dispersion (cm) in non bending plane (used for original tracking only)
66 void SetNonBendingVertexDispersion(Double_t val) {fNonBendingVertexDispersion = val;}
67 /// return the vertex dispersion (cm) in bending plane (used for original tracking only)
68 Double_t GetNonBendingVertexDispersion() const {return fNonBendingVertexDispersion;}
69 /// set the vertex dispersion (cm) in non bending plane (used for original tracking only)
70 void SetBendingVertexDispersion(Double_t val) {fBendingVertexDispersion = val;}
71 /// return the vertex dispersion (cm) in bending plane (used for original tracking only)
72 Double_t GetBendingVertexDispersion() const {return fBendingVertexDispersion;}
74 /// set the maximum distance to the track to search for compatible cluster(s) in non bending direction
75 void SetMaxNonBendingDistanceToTrack(Double_t val) {fMaxNonBendingDistanceToTrack = val;}
76 /// return the maximum distance to the track to search for compatible cluster(s) in non bending direction
77 Double_t GetMaxNonBendingDistanceToTrack() const {return fMaxNonBendingDistanceToTrack;}
78 /// set the maximum distance to the track to search for compatible cluster(s) in bending direction
79 void SetMaxBendingDistanceToTrack(Double_t val) {fMaxBendingDistanceToTrack = val;}
80 /// return the maximum distance to the track to search for compatible cluster(s) in bending direction
81 Double_t GetMaxBendingDistanceToTrack() const {return fMaxBendingDistanceToTrack;}
83 /// set the cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
84 void SetSigmaCutForTracking(Double_t val) {fSigmaCutForTracking = val;}
85 /// return the cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
86 Double_t GetSigmaCutForTracking() const {return fSigmaCutForTracking;}
88 /// switch on/off the track improvement and keep the default cut in sigma to apply on cluster (local chi2)
89 void ImproveTracks(Bool_t flag) {fImproveTracks = flag;}
90 /// switch on/off the track improvement and set the cut in sigma to apply on cluster (local chi2)
91 void ImproveTracks(Bool_t flag, Double_t sigmaCut) {fImproveTracks = flag; fSigmaCutForImprovement = sigmaCut;}
92 /// return kTRUE/kFALSE if the track improvement is switch on/off
93 Bool_t ImproveTracks() const {return fImproveTracks;}
94 /// return the cut in sigma to apply on cluster (local chi2) during track improvement
95 Double_t GetSigmaCutForImprovement() const {return fSigmaCutForImprovement;}
97 /// set the cut in sigma to apply on track during trigger hit pattern search
98 void SetSigmaCutForTrigger(Double_t val) {fSigmaCutForTrigger = val;}
99 /// return the cut in sigma to apply on track during trigger hit pattern search
100 Double_t GetSigmaCutForTrigger() const {return fSigmaCutForTrigger;}
102 /// set the maximum normalized chi2 of tracking/trigger track matching
103 void SetMaxNormChi2MatchTrigger(Double_t val) {fMaxNormChi2MatchTrigger = val;}
104 /// return the maximum normalized chi2 of tracking/trigger track matching
105 Double_t GetMaxNormChi2MatchTrigger() const {return fMaxNormChi2MatchTrigger;}
107 /// switch on/off the tracking of all the possible candidates (track only the best one if switched off)
108 void TrackAllTracks(Bool_t flag) {fTrackAllTracks = flag;}
109 /// return kTRUE/kFALSE if the tracking of all the possible candidates is switched on/off
110 Bool_t TrackAllTracks() const {return fTrackAllTracks;}
112 /// switch on/off the recovering of tracks being lost during reconstruction
113 void RecoverTracks(Bool_t flag) {fRecoverTracks = flag;}
114 /// return kTRUE/kFALSE if the recovering of tracks being lost during reconstruction is switched on/off
115 Bool_t RecoverTracks() const {return fRecoverTracks;}
117 /// switch on/off the fast building of track candidates (assuming linear propagation between stations 4 and 5)
118 void MakeTrackCandidatesFast(Bool_t flag) {fMakeTrackCandidatesFast = flag;}
119 /// return kTRUE/kFALSE if the fast building of track candidates is switched on/off
120 Bool_t MakeTrackCandidatesFast() const {return fMakeTrackCandidatesFast;}
122 /// switch on/off the building of track candidates starting from 1 cluster in each of the stations 4 and 5
123 void MakeMoreTrackCandidates(Bool_t flag) {fMakeMoreTrackCandidates = flag;}
124 /// return kTRUE/kFALSE if the building of extra track candidates is switched on/off
125 Bool_t MakeMoreTrackCandidates() const {return fMakeMoreTrackCandidates;}
127 /// switch on/off the completion of reconstructed track
128 void ComplementTracks(Bool_t flag) {fComplementTracks = flag;}
129 /// return kTRUE/kFALSE if completion of the reconstructed track is switched on/off
130 Bool_t ComplementTracks() const {return fComplementTracks;}
132 /// switch on/off the use of the smoother
133 void UseSmoother(Bool_t flag) {fUseSmoother = flag;}
134 /// return kTRUE/kFALSE if the use of the smoother is switched on/off
135 Bool_t UseSmoother() const {return fUseSmoother;}
137 /// switch on/off a chamber in the reconstruction
138 void UseChamber(Int_t iCh, Bool_t flag) {if (iCh >= 0 && iCh < 10) fUseChamber[iCh] = flag;}
139 /// return kTRUE/kFALSE whether the chamber must be used or not
140 Bool_t UseChamber(Int_t iCh) const {return (iCh >= 0 && iCh < 10) ? fUseChamber[iCh] : kFALSE;}
142 /// request or not at least one cluster in the station to validate the track
143 void RequestStation(Int_t iSt, Bool_t flag) {if (iSt >= 0 && iSt < 5) fRequestStation[iSt] = flag;}
144 /// return kTRUE/kFALSE whether at least one cluster is requested in the station to validate the track
145 Bool_t RequestStation(Int_t iSt) const {return (iSt >= 0 && iSt < 5) ? fRequestStation[iSt] : kFALSE;}
147 /// set the bypassSt45 value
148 void BypassSt45(Bool_t value) { fBypassSt45 = value; }
149 /// return kTRUE if we should replace clusters in St 4 and 5 by generated clusters from trigger tracks
150 Bool_t BypassSt45() const { return fBypassSt45; }
153 virtual void Print(Option_t *option = "") const;
157 /// clustering mode: NOCLUSTERING, PRECLUSTER, PRECLUSTERV2, PRECLUSTERV3, COG, <pre>
158 /// SIMPLEFIT, SIMPLEFITV3, MLEM:DRAW, MLEM, MLEMV2, MLEMV3 </pre>
159 TString fClusteringMode; ///< \brief name of the clustering (+ pre-clustering) mode
161 /// tracking mode: ORIGINAL, KALMAN
162 TString fTrackingMode; ///< \brief name of the tracking mode
164 Double32_t fMinBendingMomentum; ///< minimum value (GeV/c) of momentum in bending plane
165 Double32_t fMaxBendingMomentum; ///< maximum value (GeV/c) of momentum in bending plane
166 Double32_t fMaxNonBendingSlope; ///< maximum value of the non bending slope
168 Double32_t fNonBendingVertexDispersion; ///< vertex dispersion (cm) in non bending plane (used for original tracking only)
169 Double32_t fBendingVertexDispersion; ///< vertex dispersion (cm) in bending plane (used for original tracking only)
171 Double32_t fMaxNonBendingDistanceToTrack; ///< maximum distance to the track to search for compatible cluster(s) in non bending direction
172 Double32_t fMaxBendingDistanceToTrack; ///< maximum distance to the track to search for compatible cluster(s) in bending direction
174 Double32_t fSigmaCutForTracking; ///< cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking
176 Double32_t fSigmaCutForImprovement; ///< cut in sigma to apply on cluster (local chi2) during track improvement
178 Double32_t fSigmaCutForTrigger; ///< cut in sigma to apply on track during trigger hit pattern search
180 Double32_t fMaxNormChi2MatchTrigger; ///< maximum normalized chi2 of tracking/trigger track matching
182 Double32_t fPercentOfFullClusterInESD; ///< percentage of events for which all cluster info are stored in ESD
184 Bool_t fCombinedClusterTrackReco; ///< switch on/off the combined cluster/track reconstruction
186 Bool_t fTrackAllTracks; ///< kTRUE to track all the possible candidates; kFALSE to track only the best ones
188 Bool_t fRecoverTracks; ///< kTRUE to try to recover the tracks getting lost during reconstruction
190 Bool_t fMakeTrackCandidatesFast; ///< kTRUE to make candidate tracks assuming linear propagation between stations 4 and 5
192 Bool_t fMakeMoreTrackCandidates; ///< kTRUE to make candidate tracks starting from 1 cluster in each of the stations 4 and 5
194 Bool_t fComplementTracks; ///< kTRUE to try to complete the reconstructed tracks by adding missing clusters
196 Bool_t fImproveTracks; ///< kTRUE to try to improve the reconstructed tracks by removing bad clusters
198 Bool_t fUseSmoother; ///< kTRUE to use the smoother to compute track parameters/covariances and local chi2 at each cluster (used for Kalman tracking only)
200 Bool_t fSaveFullClusterInESD; ///< kTRUE to save all cluster info (including pads) in ESD
202 /// calibration mode: GAIN, NOGAIN, GAINCONSTANTCAPA
203 TString fCalibrationMode; ///<\brief calibration mode
205 Bool_t fBypassSt45; ///< kTRUE to use trigger tracks to generate "fake" clusters in St 4 and 5
207 Bool_t fUseChamber[10]; ///< kTRUE to use the chamber i in the tracking algorithm
209 Bool_t fRequestStation[5]; ///< kTRUE to request at least one cluster in station i to validate the track
212 void SetLowFluxParam();
213 void SetHighFluxParam();
216 ClassDef(AliMUONRecoParam,3) // MUON reco parameters