]>
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 | ||
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 | 15 | class 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 |