1 /**************************************************************************
2 * Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
18 //-----------------------------------------------------------------------------
19 /// \class AliMUONRecoParam
21 /// Class with MUON reconstruction parameters
23 /// \author Philippe Pillot
24 //-----------------------------------------------------------------------------
28 #include "AliMUONRecoParam.h"
29 #include "AliMUONPadStatusMaker.h"
31 #include "AliRecoParam.h"
34 #include <Riostream.h>
36 ClassImp(AliMUONRecoParam)
39 //_____________________________________________________________________________
40 AliMUONRecoParam::AliMUONRecoParam()
41 : AliDetectorRecoParam(),
42 fClusteringMode("MLEM"),
43 fTrackingMode("KALMAN"),
44 fMinBendingMomentum(0.),
45 fMaxBendingMomentum(0.),
46 fMaxNonBendingSlope(0.),
48 fNonBendingVertexDispersion(0.),
49 fBendingVertexDispersion(0.),
50 fMaxNonBendingDistanceToTrack(0.),
51 fMaxBendingDistanceToTrack(0.),
52 fSigmaCutForTracking(0.),
53 fSigmaCutForImprovement(0.),
54 fSigmaCutForTrigger(0.),
55 fStripCutForTrigger(0.),
56 fMaxStripAreaForTrigger(0.),
57 fMaxNormChi2MatchTrigger(0.),
58 fPercentOfFullClusterInESD(10.),
59 fCombinedClusterTrackReco(kFALSE),
60 fTrackAllTracks(kFALSE),
61 fRecoverTracks(kFALSE),
62 fMakeTrackCandidatesFast(kFALSE),
63 fMakeMoreTrackCandidates(kFALSE),
64 fComplementTracks(kFALSE),
65 fImproveTracks(kFALSE),
67 fSaveFullClusterInESD(kTRUE),
68 fCalibrationMode("NOGAIN"),
75 SetNameTitle("Dummy","Dummy");
79 //_____________________________________________________________________________
80 AliMUONRecoParam::~AliMUONRecoParam()
85 //_____________________________________________________________________________
87 AliMUONRecoParam::BypassSt45(Bool_t st4, Bool_t st5)
89 /// Set the bypass status
91 if ( st4 && st5 ) fBypassSt45 = 45;
92 else if ( st4 ) fBypassSt45 = 4;
93 else if ( st5 ) fBypassSt45 = 5;
97 //_____________________________________________________________________________
99 AliMUONRecoParam::GetCalibrationMode() const
101 /// Return the calibration mode. Can be :
102 /// NOGAIN : only do pedestal subtraction
103 /// GAIN : do pedestal subtraction, and apply gain correction, but with a
104 /// single capacitance value for all channels
105 /// GAINCONSTANTCAPA : as GAIN, but with a channel-dependent capacitance value
107 return fCalibrationMode.Data();
110 //_____________________________________________________________________________
111 AliMUONRecoParam *AliMUONRecoParam::GetLowFluxParam()
113 /// Return default reconstruction parameters for low flux environment
115 AliMUONRecoParam *param = new AliMUONRecoParam();
116 param->SetLowFluxParam();
121 //_____________________________________________________________________________
122 AliMUONRecoParam *AliMUONRecoParam::GetHighFluxParam()
124 /// Return default reconstruction parameters for high flux environment
126 AliMUONRecoParam *param = new AliMUONRecoParam();
127 param->SetHighFluxParam();
132 //_____________________________________________________________________________
133 AliMUONRecoParam *AliMUONRecoParam::GetCosmicParam()
135 /// Return default reconstruction parameters for high flux environment
137 AliMUONRecoParam *param = new AliMUONRecoParam();
138 param->SetCosmicParam();
143 //_____________________________________________________________________________
144 void AliMUONRecoParam::SetLowFluxParam()
146 /// Set reconstruction parameters for low flux environment
148 SetNameTitle("Low Flux","Low Flux");
149 SetEventSpecie(AliRecoParam::kLowMult);
150 fMinBendingMomentum = 1.;
151 fMaxBendingMomentum = 3000.;
152 fMaxNonBendingSlope = 0.3;
153 fMaxBendingSlope = 0.4;
154 fNonBendingVertexDispersion = 10.;
155 fBendingVertexDispersion = 10.;
156 fMaxNonBendingDistanceToTrack = 1.;
157 fMaxBendingDistanceToTrack = 1.;
158 fSigmaCutForTracking = 6.;
159 fSigmaCutForImprovement = 5.;
160 fSigmaCutForTrigger = 8.;
161 fStripCutForTrigger = 1.;
162 fMaxStripAreaForTrigger = 3.;
163 fMaxNormChi2MatchTrigger = 16.;
164 fCombinedClusterTrackReco = kFALSE;
165 fTrackAllTracks = kTRUE;
166 fRecoverTracks = kTRUE;
167 fMakeTrackCandidatesFast = kFALSE;
168 fMakeMoreTrackCandidates = kFALSE;
169 fComplementTracks = kTRUE;
170 fImproveTracks = kTRUE;
171 fUseSmoother = kTRUE;
172 for (Int_t iCh = 0; iCh < 10; iCh++) {
173 fUseChamber[iCh] = kTRUE;
174 fDefaultNonBendingReso[iCh] = 0.144;
175 fDefaultBendingReso[iCh] = 0.01;
177 for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
182 //_____________________________________________________________________________
183 void AliMUONRecoParam::SetHighFluxParam()
185 /// Set reconstruction parameters for high flux environment
187 SetNameTitle("High Flux","High Flux");
188 SetEventSpecie(AliRecoParam::kHighMult);
189 fMinBendingMomentum = 1.;
190 fMaxBendingMomentum = 3000.;
191 fMaxNonBendingSlope = 0.3;
192 fMaxBendingSlope = 0.4;
193 fNonBendingVertexDispersion = 10.;
194 fBendingVertexDispersion = 10.;
195 fMaxNonBendingDistanceToTrack = 1.;
196 fMaxBendingDistanceToTrack = 1.;
197 fSigmaCutForTracking = 6.;
198 fSigmaCutForImprovement = 5.;
199 fSigmaCutForTrigger = 8.;
200 fStripCutForTrigger = 1.;
201 fMaxStripAreaForTrigger = 3.;
202 fMaxNormChi2MatchTrigger = 16.;
203 fCombinedClusterTrackReco = kFALSE;
204 fTrackAllTracks = kTRUE;
205 fRecoverTracks = kTRUE;
206 fMakeTrackCandidatesFast = kFALSE;
207 fMakeMoreTrackCandidates = kFALSE;
208 fComplementTracks = kTRUE;
209 fImproveTracks = kTRUE;
210 fUseSmoother = kTRUE;
211 for (Int_t iCh = 0; iCh < 10; iCh++) {
212 fUseChamber[iCh] = kTRUE;
213 fDefaultNonBendingReso[iCh] = 0.144;
214 fDefaultBendingReso[iCh] = 0.01;
216 for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
221 //_____________________________________________________________________________
223 AliMUONRecoParam::RequestedStationMask() const
225 /// Get the mask of the requested station, i.e. an integer where
226 /// bit n is set to one if the station n was requested
230 for ( Int_t i = 0; i < 5; ++i )
232 if ( RequestStation(i) ) m |= ( 1 << i );
237 //_____________________________________________________________________________
238 void AliMUONRecoParam::SetCosmicParam()
240 /// Set reconstruction parameters for high flux environment
242 SetNameTitle("Cosmic","Cosmic");
243 SetEventSpecie(AliRecoParam::kCosmic);
244 fMinBendingMomentum = 1.;
245 fMaxBendingMomentum = 10000000.;
246 fMaxNonBendingSlope = 0.3;
247 fMaxBendingSlope = 0.4;
248 fNonBendingVertexDispersion = 10.;
249 fBendingVertexDispersion = 10.;
250 fMaxNonBendingDistanceToTrack = 10.;
251 fMaxBendingDistanceToTrack = 10.;
252 fSigmaCutForTracking = 20.;
253 fSigmaCutForImprovement = 20.;
254 fSigmaCutForTrigger = 8.;
255 fMaxNormChi2MatchTrigger = 16.;
256 fPercentOfFullClusterInESD = 100.;
257 fCombinedClusterTrackReco = kFALSE;
258 fTrackAllTracks = kTRUE;
259 fRecoverTracks = kTRUE;
260 fMakeTrackCandidatesFast = kFALSE;
261 fMakeMoreTrackCandidates = kFALSE;
262 fComplementTracks = kTRUE;
263 fImproveTracks = kTRUE;
264 fUseSmoother = kTRUE;
265 fSaveFullClusterInESD = kTRUE;
266 for (Int_t iCh = 0; iCh < 10; iCh++) {
267 fUseChamber[iCh] = kTRUE;
268 fDefaultNonBendingReso[iCh] = 0.144;
269 fDefaultBendingReso[iCh] = 0.01;
271 for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
276 //_____________________________________________________________________________
277 void AliMUONRecoParam::Print(Option_t *option) const
279 /// print reconstruction parameters
280 /// if option = FULL then print also unused parameters
282 cout<<endl<<"\t------MUON Reconstruction parameters ("<<GetName()<<")------"<<endl;
284 if (IsDefault()) cout<<"\t\t*** Parameters used by default ***"<<endl;
286 cout<<Form("Calibration mode = %s",fCalibrationMode.Data())<<endl;
287 cout<<Form("Clustering mode = %s",fClusteringMode.Data())<<endl;
288 cout<<Form("Tracking mode = %s",fTrackingMode.Data())<<endl;
294 bypass = "stations 4 and 5";
296 else if ( BypassSt4() )
298 bypass = "station 4";
300 else if ( BypassSt5() )
302 bypass = "station 5";
305 if (bypass.Length()) cout << "Will bypass " << bypass.Data() << " (replacing real clusters by generated ones from trigger tracks)" << endl;
307 if (fCombinedClusterTrackReco) cout<<"Combined cluster/track reconstruction: ON"<<endl;
308 else cout<<"Combined cluster/track reconstruction: OFF"<<endl;
310 if (fSaveFullClusterInESD) cout<<Form("Save all cluster info in ESD for %5.2f %% of events",fPercentOfFullClusterInESD)<<endl;
311 else cout<<"Save partial cluster info in ESD"<<endl;
313 cout<<Form("Bending momentum range = [%5.2f,%5.2f]",fMinBendingMomentum,fMaxBendingMomentum)<<endl;
315 cout<<Form("Maximum non bending slope = %5.2f",fMaxNonBendingSlope)<<endl;
317 cout<<Form("Maximum bending slope (used only if B=0) = %5.2f",fMaxBendingSlope)<<endl;
319 if (strstr(fTrackingMode,"ORIGINAL"))
320 cout<<Form("Vertex dispertion = (%5.2f,%5.2f)",fNonBendingVertexDispersion,fBendingVertexDispersion)<<endl;
321 else if (strstr(option,"FULL"))
322 cout<<Form("Vertex dispertion (used for original tracking only) = (%5.2f,%5.2f)",fNonBendingVertexDispersion,fBendingVertexDispersion)<<endl;
324 cout<<Form("Maximum distance to track = (%5.2f,%5.2f)",fMaxNonBendingDistanceToTrack,fMaxBendingDistanceToTrack)<<endl;
326 cout<<Form("Sigma cut for tracking = %5.2f",fSigmaCutForTracking)<<endl;
328 cout<<Form("Sigma cut for trigger hit pattern = %5.2f",fSigmaCutForTrigger)<<endl;
330 cout<<Form("Cut in strips for trigger chamber efficiency = %5.2f",fStripCutForTrigger)<<endl;
332 cout<<Form("Max search area in strips for trigger chamber efficiency = %5.2f",fMaxStripAreaForTrigger)<<endl;
334 if (fTrackAllTracks) cout<<"Track all the possible candidates"<<endl;
335 else cout<<"Track only the best candidates"<<endl;
337 if (strstr(option,"FULL")) {
338 cout<<"Make track candidates assuming linear propagation between stations 4 and 5: ";
339 if (fMakeTrackCandidatesFast) cout<<"ON"<<endl;
340 else cout<<"OFF"<<endl;
341 } else if (fMakeTrackCandidatesFast)
342 cout<<"Make track candidates assuming linear propagation between stations 4 and 5"<<endl;
344 if (strstr(option,"FULL")) {
345 cout<<"Make track candidates starting from 1 cluster in each of the stations 4 and 5: ";
346 if (fMakeMoreTrackCandidates) cout<<"ON"<<endl;
347 else cout<<"OFF"<<endl;
348 } else if (fMakeMoreTrackCandidates)
349 cout<<"Make track candidates starting from 1 cluster in each of the stations 4 and 5"<<endl;
351 if (strstr(option,"FULL")) {
352 cout<<"Try to recover tracks getting lost during tracking: ";
353 if (fRecoverTracks) cout<<"ON"<<endl;
354 else cout<<"OFF"<<endl;
355 } else if (fRecoverTracks)
356 cout<<"Try to recover tracks getting lost during tracking"<<endl;
358 if (strstr(option,"FULL")) {
359 cout<<"Try to complete the reconstructed tracks by adding missing clusters: ";
360 if (fComplementTracks) cout<<"ON"<<endl;
361 else cout<<"OFF"<<endl;
362 } else if (fComplementTracks)
363 cout<<"Try to complete the reconstructed tracks by adding missing clusters"<<endl;
365 if (strstr(option,"FULL")) {
366 cout<<"Try to improve the reconstructed tracks by removing bad clusters: ";
367 if (fImproveTracks) cout<<Form("ON (sigma cut = %5.2f)",fSigmaCutForImprovement)<<endl;
368 else cout<<"OFF"<<endl;
369 } else if (fImproveTracks)
370 cout<<Form("Try to improve the reconstructed tracks by removing bad clusters (sigma cut = %5.2f)",fSigmaCutForImprovement)<<endl;
372 if (strstr(option,"FULL")) {
373 cout<<"Use smoother to compute final track parameters, etc, at each cluster (used for Kalman tracking only): ";
374 if (fUseSmoother) cout<<"ON"<<endl;
375 else cout<<"OFF"<<endl;
376 } else if (fUseSmoother)
377 cout<<"Use smoother to compute final track parameters, etc, at each cluster"<<endl;
379 cout<<Form("Maximum normalized chi2 of tracking/trigger track matching = %5.2f",fMaxNormChi2MatchTrigger)<<endl;
381 Bool_t discardedCh = kFALSE;
384 if (!UseChamber(ch)) {
386 cout<<"Discarded chambers(1..): "<<ch+1;
389 else cout<<" "<<ch+1;
392 if (discardedCh) cout<<endl;
394 Bool_t discardedSt = kFALSE;
397 if (!RequestStation(st)) {
399 cout<<"Not requested stations(1..): "<<st+1;
402 else cout<<" "<<st+1;
405 if (discardedSt) cout<<endl;
407 cout << Form("Pad goodness policy mask is 0x%x",PadGoodnessMask()) << endl;
408 cout << "Which means we reject pads having the condition = " <<
409 AliMUONPadStatusMaker::AsCondition(PadGoodnessMask()).Data() << endl;
411 cout << "The pad limits we are using are :" << endl;
413 cout << Form("%5.0f <= HVSt12 <= %5.0f Volts",HVSt12LowLimit(),HVSt12HighLimit()) << endl;
414 cout << Form("%5.0f <= HVSt345 <= %5.0f Volts",HVSt345LowLimit(),HVSt345HighLimit()) << endl;
415 cout << Form("%7.2f <= Pedestal mean <= %7.2f",PedMeanLowLimit(),PedMeanHighLimit()) << endl;
416 cout << Form("%7.2f <= Pedestal sigma <= %7.2f",PedSigmaLowLimit(),PedSigmaHighLimit()) << endl;
417 cout << Form("%e <= Gain linear term <= %e",GainA1LowLimit(),GainA1HighLimit()) << endl;
418 cout << Form("%e <= Gain quadratic term <= %e",GainA2LowLimit(),GainA2HighLimit()) << endl;
419 cout << Form("%5.0f <= Gain threshold term <= %5.0f",GainThresLowLimit(),GainThresHighLimit()) << endl;
421 cout << Form("And we cut on charge >= %7.2f x ( pedestal sigma ) ",ChargeSigmaCut()) << endl;
423 cout << "chamber non bending resolution = |";
424 for (Int_t iCh = 0; iCh < 10; iCh++) cout << Form(" %6.3f |",fDefaultNonBendingReso[iCh]);
426 cout << "chamber bending resolution = |";
427 for (Int_t iCh = 0; iCh < 10; iCh++) cout << Form(" %6.3f |",fDefaultBendingReso[iCh]);
430 cout<<"\t-----------------------------------------------------"<<endl<<endl;
434 //_____________________________________________________________________________
436 AliMUONRecoParam::SetDefaultLimits()
438 /// Set the default limits and pad goodness policy
440 fHVSt12Limits[0]=1500;
441 fHVSt12Limits[1]=2000;
443 fHVSt345Limits[0]=1500;
444 fHVSt345Limits[1]=2000;
446 fPedMeanLimits[0] = 50;
447 fPedMeanLimits[1] = 1024;
449 fPedSigmaLimits[0] = 0.1;
450 fPedSigmaLimits[1] = 100;
452 fGainA1Limits[0] = 0.1;
453 fGainA1Limits[1] = 10;
455 fGainA2Limits[0] = -1E30;
456 fGainA2Limits[1] = 1E30;
458 fGainThresLimits[0] = 0;
459 fGainThresLimits[1] = 4095;
461 fPadGoodnessMask = 0x8080;
463 fChargeSigmaCut = 4.0;