]>
Commit | Line | Data |
---|---|---|
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 | ||
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" |
16 | ||
15d30ed4 | 17 | class AliMUONRecoParam : public AliDetectorRecoParam |
3304fa09 | 18 | { |
19 | public: | |
20 | AliMUONRecoParam(); | |
21 | virtual ~AliMUONRecoParam(); | |
22 | ||
23 | static AliMUONRecoParam *GetLowFluxParam(); | |
24 | static AliMUONRecoParam *GetHighFluxParam(); | |
0e894e58 | 25 | static AliMUONRecoParam *GetCosmicParam(); |
3304fa09 | 26 | |
9bf6860b | 27 | /// set the calibration mode (see GetCalibrationMode() for possible modes) |
7332f213 | 28 | void SetCalibrationMode(Option_t* mode) { fCalibrationMode = mode; fCalibrationMode.ToUpper();} |
9bf6860b | 29 | |
30 | Option_t* GetCalibrationMode() const; | |
31 | ||
3304fa09 | 32 | /// set the clustering (pre-clustering) mode |
7332f213 | 33 | void SetClusteringMode(Option_t* mode) {fClusteringMode = mode; fClusteringMode.ToUpper();} |
3304fa09 | 34 | /// get the clustering (pre-clustering) mode |
35 | Option_t* GetClusteringMode() const {return fClusteringMode.Data();} | |
36 | ||
37 | /// set the tracking mode | |
7332f213 | 38 | void SetTrackingMode(Option_t* mode) {fTrackingMode = mode; fTrackingMode.ToUpper();} |
3304fa09 | 39 | /// get the tracking mode |
40 | Option_t* GetTrackingMode() const {return fTrackingMode.Data();} | |
41 | ||
38bcf0ef | 42 | /// switch on/off the combined cluster/track reconstruction |
43 | void CombineClusterTrackReco(Bool_t flag) {fCombinedClusterTrackReco = flag;} | |
44 | /// return kTRUE/kFALSE if the combined cluster/track reconstruction is on/off | |
45 | Bool_t CombineClusterTrackReco() const {return fCombinedClusterTrackReco;} | |
46 | ||
0a18ba02 | 47 | /// save all cluster info (including pads) in ESD, for the given percentage of events |
48 | void SaveFullClusterInESD(Bool_t flag, Double_t percentOfEvent = 100.) {fSaveFullClusterInESD = flag; | |
49 | fPercentOfFullClusterInESD = (fSaveFullClusterInESD) ? percentOfEvent : 0.;} | |
50 | /// return kTRUE/kFALSE depending on whether we save all cluster info in ESD or not | |
51 | Bool_t SaveFullClusterInESD() const {return fSaveFullClusterInESD;} | |
52 | /// return the percentage of events for which all cluster info are stored in ESD | |
53 | Double_t GetPercentOfFullClusterInESD() const {return fPercentOfFullClusterInESD;} | |
54 | ||
0e894e58 | 55 | /// set the most probable value (GeV/c) of momentum in bending plane |
56 | /// needed to get some "reasonable" corrections for MCS and E loss even if B = 0 | |
57 | void SetMostProbBendingMomentum(Double_t val) {fMostProbBendingMomentum = val;} | |
58 | /// return the most probable value (GeV/c) of momentum in bending plane | |
59 | Double_t GetMostProbBendingMomentum() const {return fMostProbBendingMomentum;} | |
60 | ||
3304fa09 | 61 | /// set the minimum value (GeV/c) of momentum in bending plane |
62 | void SetMinBendingMomentum(Double_t val) {fMinBendingMomentum = val;} | |
63 | /// return the minimum value (GeV/c) of momentum in bending plane | |
64 | Double_t GetMinBendingMomentum() const {return fMinBendingMomentum;} | |
65 | /// set the maximum value (GeV/c) of momentum in bending plane | |
66 | void SetMaxBendingMomentum(Double_t val) {fMaxBendingMomentum = val;} | |
67 | /// return the maximum value (GeV/c) of momentum in bending plane | |
68 | Double_t GetMaxBendingMomentum() const {return fMaxBendingMomentum;} | |
9bf6860b | 69 | /// set the maximum value of the non bending slope |
70 | void SetMaxNonBendingSlope(Double_t val) {fMaxNonBendingSlope = val;} | |
71 | /// return the maximum value of the non bending slope | |
72 | Double_t GetMaxNonBendingSlope() const {return fMaxNonBendingSlope;} | |
9f093251 | 73 | /// set the maximum value of the bending slope |
74 | void SetMaxBendingSlope(Double_t val) {fMaxBendingSlope = val;} | |
75 | /// return the maximum value of the bending slope | |
76 | Double_t GetMaxBendingSlope() const {return fMaxBendingSlope;} | |
3304fa09 | 77 | |
78 | /// set the vertex dispersion (cm) in non bending plane (used for original tracking only) | |
79 | void SetNonBendingVertexDispersion(Double_t val) {fNonBendingVertexDispersion = val;} | |
80 | /// return the vertex dispersion (cm) in bending plane (used for original tracking only) | |
81 | Double_t GetNonBendingVertexDispersion() const {return fNonBendingVertexDispersion;} | |
82 | /// set the vertex dispersion (cm) in non bending plane (used for original tracking only) | |
83 | void SetBendingVertexDispersion(Double_t val) {fBendingVertexDispersion = val;} | |
84 | /// return the vertex dispersion (cm) in bending plane (used for original tracking only) | |
85 | Double_t GetBendingVertexDispersion() const {return fBendingVertexDispersion;} | |
86 | ||
87 | /// set the maximum distance to the track to search for compatible cluster(s) in non bending direction | |
88 | void SetMaxNonBendingDistanceToTrack(Double_t val) {fMaxNonBendingDistanceToTrack = val;} | |
89 | /// return the maximum distance to the track to search for compatible cluster(s) in non bending direction | |
90 | Double_t GetMaxNonBendingDistanceToTrack() const {return fMaxNonBendingDistanceToTrack;} | |
91 | /// set the maximum distance to the track to search for compatible cluster(s) in bending direction | |
92 | void SetMaxBendingDistanceToTrack(Double_t val) {fMaxBendingDistanceToTrack = val;} | |
93 | /// return the maximum distance to the track to search for compatible cluster(s) in bending direction | |
94 | Double_t GetMaxBendingDistanceToTrack() const {return fMaxBendingDistanceToTrack;} | |
95 | ||
96 | /// set the cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking | |
97 | void SetSigmaCutForTracking(Double_t val) {fSigmaCutForTracking = val;} | |
98 | /// return the cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking | |
99 | Double_t GetSigmaCutForTracking() const {return fSigmaCutForTracking;} | |
100 | ||
101 | /// switch on/off the track improvement and keep the default cut in sigma to apply on cluster (local chi2) | |
102 | void ImproveTracks(Bool_t flag) {fImproveTracks = flag;} | |
103 | /// switch on/off the track improvement and set the cut in sigma to apply on cluster (local chi2) | |
104 | void ImproveTracks(Bool_t flag, Double_t sigmaCut) {fImproveTracks = flag; fSigmaCutForImprovement = sigmaCut;} | |
105 | /// return kTRUE/kFALSE if the track improvement is switch on/off | |
106 | Bool_t ImproveTracks() const {return fImproveTracks;} | |
107 | /// return the cut in sigma to apply on cluster (local chi2) during track improvement | |
108 | Double_t GetSigmaCutForImprovement() const {return fSigmaCutForImprovement;} | |
fda59e58 | 109 | |
110 | /// set the cut in sigma to apply on track during trigger hit pattern search | |
111 | void SetSigmaCutForTrigger(Double_t val) {fSigmaCutForTrigger = val;} | |
112 | /// return the cut in sigma to apply on track during trigger hit pattern search | |
113 | Double_t GetSigmaCutForTrigger() const {return fSigmaCutForTrigger;} | |
b5270f21 | 114 | /// set the cut in strips to apply on trigger track during trigger chamber efficiency |
115 | void SetStripCutForTrigger(Double_t val) {fStripCutForTrigger = val;} | |
116 | /// return the cut in strips to apply on trigger track during trigger chamber efficiency | |
117 | Double_t GetStripCutForTrigger() const {return fStripCutForTrigger;} | |
118 | /// set the maximum search area in strips to apply on trigger track during trigger chamber efficiency | |
119 | void SetMaxStripAreaForTrigger(Double_t val) {fMaxStripAreaForTrigger = val;} | |
120 | /// return the maximum search area in strips to apply on trigger track during trigger chamber efficiency | |
121 | Double_t GetMaxStripAreaForTrigger() const {return fMaxStripAreaForTrigger;} | |
3304fa09 | 122 | |
123 | /// set the maximum normalized chi2 of tracking/trigger track matching | |
124 | void SetMaxNormChi2MatchTrigger(Double_t val) {fMaxNormChi2MatchTrigger = val;} | |
125 | /// return the maximum normalized chi2 of tracking/trigger track matching | |
126 | Double_t GetMaxNormChi2MatchTrigger() const {return fMaxNormChi2MatchTrigger;} | |
127 | ||
128 | /// switch on/off the tracking of all the possible candidates (track only the best one if switched off) | |
129 | void TrackAllTracks(Bool_t flag) {fTrackAllTracks = flag;} | |
38bcf0ef | 130 | /// return kTRUE/kFALSE if the tracking of all the possible candidates is switched on/off |
3304fa09 | 131 | Bool_t TrackAllTracks() const {return fTrackAllTracks;} |
132 | ||
133 | /// switch on/off the recovering of tracks being lost during reconstruction | |
134 | void RecoverTracks(Bool_t flag) {fRecoverTracks = flag;} | |
38bcf0ef | 135 | /// return kTRUE/kFALSE if the recovering of tracks being lost during reconstruction is switched on/off |
3304fa09 | 136 | Bool_t RecoverTracks() const {return fRecoverTracks;} |
137 | ||
138 | /// switch on/off the fast building of track candidates (assuming linear propagation between stations 4 and 5) | |
139 | void MakeTrackCandidatesFast(Bool_t flag) {fMakeTrackCandidatesFast = flag;} | |
38bcf0ef | 140 | /// return kTRUE/kFALSE if the fast building of track candidates is switched on/off |
3304fa09 | 141 | Bool_t MakeTrackCandidatesFast() const {return fMakeTrackCandidatesFast;} |
142 | ||
9bf6860b | 143 | /// switch on/off the building of track candidates starting from 1 cluster in each of the stations 4 and 5 |
144 | void MakeMoreTrackCandidates(Bool_t flag) {fMakeMoreTrackCandidates = flag;} | |
145 | /// return kTRUE/kFALSE if the building of extra track candidates is switched on/off | |
146 | Bool_t MakeMoreTrackCandidates() const {return fMakeMoreTrackCandidates;} | |
147 | ||
3304fa09 | 148 | /// switch on/off the completion of reconstructed track |
149 | void ComplementTracks(Bool_t flag) {fComplementTracks = flag;} | |
38bcf0ef | 150 | /// return kTRUE/kFALSE if completion of the reconstructed track is switched on/off |
3304fa09 | 151 | Bool_t ComplementTracks() const {return fComplementTracks;} |
152 | ||
153 | /// switch on/off the use of the smoother | |
154 | void UseSmoother(Bool_t flag) {fUseSmoother = flag;} | |
38bcf0ef | 155 | /// return kTRUE/kFALSE if the use of the smoother is switched on/off |
3304fa09 | 156 | Bool_t UseSmoother() const {return fUseSmoother;} |
157 | ||
9bf6860b | 158 | /// switch on/off a chamber in the reconstruction |
159 | void UseChamber(Int_t iCh, Bool_t flag) {if (iCh >= 0 && iCh < 10) fUseChamber[iCh] = flag;} | |
160 | /// return kTRUE/kFALSE whether the chamber must be used or not | |
161 | Bool_t UseChamber(Int_t iCh) const {return (iCh >= 0 && iCh < 10) ? fUseChamber[iCh] : kFALSE;} | |
162 | ||
163 | /// request or not at least one cluster in the station to validate the track | |
164 | void RequestStation(Int_t iSt, Bool_t flag) {if (iSt >= 0 && iSt < 5) fRequestStation[iSt] = flag;} | |
165 | /// return kTRUE/kFALSE whether at least one cluster is requested in the station to validate the track | |
166 | Bool_t RequestStation(Int_t iSt) const {return (iSt >= 0 && iSt < 5) ? fRequestStation[iSt] : kFALSE;} | |
167 | ||
6cac085d | 168 | /// set the bypassSt45 value |
169 | void BypassSt45(Bool_t value) { fBypassSt45 = value; } | |
9bf6860b | 170 | /// return kTRUE if we should replace clusters in St 4 and 5 by generated clusters from trigger tracks |
171 | Bool_t BypassSt45() const { return fBypassSt45; } | |
3304fa09 | 172 | |
9bf6860b | 173 | |
174 | virtual void Print(Option_t *option = "") const; | |
3304fa09 | 175 | |
176 | private: | |
177 | ||
178 | /// clustering mode: NOCLUSTERING, PRECLUSTER, PRECLUSTERV2, PRECLUSTERV3, COG, <pre> | |
179 | /// SIMPLEFIT, SIMPLEFITV3, MLEM:DRAW, MLEM, MLEMV2, MLEMV3 </pre> | |
180 | TString fClusteringMode; ///< \brief name of the clustering (+ pre-clustering) mode | |
181 | ||
182 | /// tracking mode: ORIGINAL, KALMAN | |
183 | TString fTrackingMode; ///< \brief name of the tracking mode | |
184 | ||
0e894e58 | 185 | Double32_t fMostProbBendingMomentum; ///< most probable value (GeV/c) of muon momentum in bending plane (used when B = 0) |
186 | ||
3304fa09 | 187 | Double32_t fMinBendingMomentum; ///< minimum value (GeV/c) of momentum in bending plane |
188 | Double32_t fMaxBendingMomentum; ///< maximum value (GeV/c) of momentum in bending plane | |
9bf6860b | 189 | Double32_t fMaxNonBendingSlope; ///< maximum value of the non bending slope |
9f093251 | 190 | Double32_t fMaxBendingSlope; ///< maximum value of the bending slope (used only if B = 0) |
3304fa09 | 191 | |
192 | Double32_t fNonBendingVertexDispersion; ///< vertex dispersion (cm) in non bending plane (used for original tracking only) | |
193 | Double32_t fBendingVertexDispersion; ///< vertex dispersion (cm) in bending plane (used for original tracking only) | |
194 | ||
195 | Double32_t fMaxNonBendingDistanceToTrack; ///< maximum distance to the track to search for compatible cluster(s) in non bending direction | |
196 | Double32_t fMaxBendingDistanceToTrack; ///< maximum distance to the track to search for compatible cluster(s) in bending direction | |
197 | ||
198 | Double32_t fSigmaCutForTracking; ///< cut in sigma to apply on cluster (local chi2) and track (global chi2) during tracking | |
fda59e58 | 199 | |
3304fa09 | 200 | Double32_t fSigmaCutForImprovement; ///< cut in sigma to apply on cluster (local chi2) during track improvement |
201 | ||
fda59e58 | 202 | Double32_t fSigmaCutForTrigger; ///< cut in sigma to apply on track during trigger hit pattern search |
b5270f21 | 203 | |
204 | Double32_t fStripCutForTrigger; ///< cut in strips to apply on trigger track during trigger chamber efficiency | |
205 | ||
206 | Double32_t fMaxStripAreaForTrigger; ///< max. search area in strips to apply on trigger track during trigger chamber efficiency | |
fda59e58 | 207 | |
3304fa09 | 208 | Double32_t fMaxNormChi2MatchTrigger; ///< maximum normalized chi2 of tracking/trigger track matching |
209 | ||
0a18ba02 | 210 | Double32_t fPercentOfFullClusterInESD; ///< percentage of events for which all cluster info are stored in ESD |
211 | ||
38bcf0ef | 212 | Bool_t fCombinedClusterTrackReco; ///< switch on/off the combined cluster/track reconstruction |
213 | ||
3304fa09 | 214 | Bool_t fTrackAllTracks; ///< kTRUE to track all the possible candidates; kFALSE to track only the best ones |
215 | ||
216 | Bool_t fRecoverTracks; ///< kTRUE to try to recover the tracks getting lost during reconstruction | |
217 | ||
218 | Bool_t fMakeTrackCandidatesFast; ///< kTRUE to make candidate tracks assuming linear propagation between stations 4 and 5 | |
219 | ||
9bf6860b | 220 | Bool_t fMakeMoreTrackCandidates; ///< kTRUE to make candidate tracks starting from 1 cluster in each of the stations 4 and 5 |
221 | ||
3304fa09 | 222 | Bool_t fComplementTracks; ///< kTRUE to try to complete the reconstructed tracks by adding missing clusters |
223 | ||
224 | Bool_t fImproveTracks; ///< kTRUE to try to improve the reconstructed tracks by removing bad clusters | |
225 | ||
226 | Bool_t fUseSmoother; ///< kTRUE to use the smoother to compute track parameters/covariances and local chi2 at each cluster (used for Kalman tracking only) | |
227 | ||
0a18ba02 | 228 | Bool_t fSaveFullClusterInESD; ///< kTRUE to save all cluster info (including pads) in ESD |
3304fa09 | 229 | |
9bf6860b | 230 | /// calibration mode: GAIN, NOGAIN, GAINCONSTANTCAPA |
de98fdc9 | 231 | TString fCalibrationMode; ///<\brief calibration mode |
232 | ||
9bf6860b | 233 | Bool_t fBypassSt45; ///< kTRUE to use trigger tracks to generate "fake" clusters in St 4 and 5 |
234 | ||
235 | Bool_t fUseChamber[10]; ///< kTRUE to use the chamber i in the tracking algorithm | |
236 | ||
237 | Bool_t fRequestStation[5]; ///< kTRUE to request at least one cluster in station i to validate the track | |
238 | ||
3304fa09 | 239 | // functions |
240 | void SetLowFluxParam(); | |
241 | void SetHighFluxParam(); | |
0e894e58 | 242 | void SetCosmicParam(); |
3304fa09 | 243 | |
244 | ||
b5270f21 | 245 | ClassDef(AliMUONRecoParam,6) // MUON reco parameters |
3304fa09 | 246 | }; |
247 | ||
248 | #endif | |
249 |