/************************************************************************** * Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. * * * * Author: The ALICE Off-line Project. * * Contributors are mentioned in the code where appropriate. * * * * Permission to use, copy, modify and distribute this software and its * * documentation strictly for non-commercial purposes is hereby granted * * without fee, provided that the above copyright notice appears in all * * copies and that both the copyright notice and this permission notice * * appear in the supporting documentation. The authors make no claims * * about the suitability of this software for any purpose. It is * * provided "as is" without express or implied warranty. * **************************************************************************/ // $Id$ //----------------------------------------------------------------------------- /// \class AliMUONRecoParam /// /// Class with MUON reconstruction parameters /// /// \author Philippe Pillot //----------------------------------------------------------------------------- #include "AliMUONRecoParam.h" #include "AliMUONPadStatusMaker.h" #include "AliRecoParam.h" #include "AliLog.h" #include ClassImp(AliMUONRecoParam) //_____________________________________________________________________________ AliMUONRecoParam::AliMUONRecoParam() : AliDetectorRecoParam(), fClusteringMode("MLEM"), fTrackingMode("KALMAN"), fMinBendingMomentum(0.), fMaxBendingMomentum(0.), fMaxNonBendingSlope(0.), fMaxBendingSlope(0.), fNonBendingVertexDispersion(0.), fBendingVertexDispersion(0.), fMaxNonBendingDistanceToTrack(0.), fMaxBendingDistanceToTrack(0.), fSigmaCutForTracking(0.), fSigmaCutForImprovement(0.), fSigmaCutForTrigger(0.), fStripCutForTrigger(0.), fMaxStripAreaForTrigger(0.), fMaxNormChi2MatchTrigger(0.), fPercentOfFullClusterInESD(10.), fCombinedClusterTrackReco(kFALSE), fTrackAllTracks(kFALSE), fRecoverTracks(kFALSE), fMakeTrackCandidatesFast(kFALSE), fMakeMoreTrackCandidates(kFALSE), fComplementTracks(kFALSE), fImproveTracks(kFALSE), fUseSmoother(kFALSE), fSaveFullClusterInESD(kTRUE), fCalibrationMode("NOGAIN"), fBypassSt45(0), fPadGoodnessMask(0), fChargeSigmaCut(4.0), fRemoveConnectedTracksInSt12(kFALSE), fMaxTriggerTracks(0), fMaxTrackCandidates(0), fSelectTrackOnSlope(kFALSE), fMissingPadFractionLimit(0), fFractionOfBuspatchOutsideOccupancyLimit(0) { /// Constructor SetNameTitle("Dummy","Dummy"); for (Int_t iCh = 0; iCh < 10; iCh++) { fUseChamber[iCh] = kTRUE; fDefaultNonBendingReso[iCh] = 0.; fDefaultBendingReso[iCh] = 0.; } for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE; SetDefaultLimits(); } //_____________________________________________________________________________ AliMUONRecoParam::~AliMUONRecoParam() { /// Destructor } //_____________________________________________________________________________ void AliMUONRecoParam::BypassSt45(Bool_t st4, Bool_t st5) { /// Set the bypass status if ( st4 && st5 ) fBypassSt45 = 45; else if ( st4 ) fBypassSt45 = 4; else if ( st5 ) fBypassSt45 = 5; else fBypassSt45 = 0; } //_____________________________________________________________________________ Option_t* AliMUONRecoParam::GetCalibrationMode() const { /// Return the calibration mode. Can be : /// NOGAIN : only do pedestal subtraction /// GAIN : do pedestal subtraction, and apply gain correction, but with a /// single capacitance value for all channels /// INJECTIONGAIN : as GAIN, but with gain values taken as EMELEC factory values /// GAINCONSTANTCAPA : as GAIN, but with a channel-dependent capacitance value return fCalibrationMode.Data(); } //_____________________________________________________________________________ AliMUONRecoParam *AliMUONRecoParam::GetLowFluxParam() { /// Return default reconstruction parameters for low flux environment AliMUONRecoParam *param = new AliMUONRecoParam(); param->SetLowFluxParam(); return param; } //_____________________________________________________________________________ AliMUONRecoParam *AliMUONRecoParam::GetHighFluxParam() { /// Return default reconstruction parameters for high flux environment AliMUONRecoParam *param = new AliMUONRecoParam(); param->SetHighFluxParam(); return param; } //_____________________________________________________________________________ AliMUONRecoParam *AliMUONRecoParam::GetCosmicParam() { /// Return default reconstruction parameters for high flux environment AliMUONRecoParam *param = new AliMUONRecoParam(); param->SetCosmicParam(); return param; } //_____________________________________________________________________________ AliMUONRecoParam *AliMUONRecoParam::GetCalibrationParam() { /// Return default (dummy) reconstruction parameters for calibration environment AliMUONRecoParam *param = new AliMUONRecoParam(); param->SetCalibrationParam(); return param; } //_____________________________________________________________________________ void AliMUONRecoParam::SetLowFluxParam() { /// Set reconstruction parameters for low flux environment SetNameTitle("Low Flux","Low Flux"); SetEventSpecie(AliRecoParam::kLowMult); fMinBendingMomentum = 0.8; fMaxBendingMomentum = 1.e10; fMaxNonBendingSlope = 0.3; fMaxBendingSlope = 0.4; fSelectTrackOnSlope = kFALSE; fNonBendingVertexDispersion = 70.; fBendingVertexDispersion = 70.; fMaxNonBendingDistanceToTrack = 1.; fMaxBendingDistanceToTrack = 1.; fSigmaCutForTracking = 6.; fSigmaCutForImprovement = 5.; fSigmaCutForTrigger = 8.; fStripCutForTrigger = 1.; fMaxStripAreaForTrigger = 3.; fMaxNormChi2MatchTrigger = 16.; fCombinedClusterTrackReco = kFALSE; fTrackAllTracks = kTRUE; fRecoverTracks = kTRUE; fMakeTrackCandidatesFast = kFALSE; fMakeMoreTrackCandidates = kFALSE; fComplementTracks = kTRUE; fImproveTracks = kTRUE; fRemoveConnectedTracksInSt12 = kTRUE; fUseSmoother = kTRUE; for (Int_t iCh = 0; iCh < 10; iCh++) { fUseChamber[iCh] = kTRUE; fDefaultNonBendingReso[iCh] = 0.144; fDefaultBendingReso[iCh] = 0.01; } for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE; fBypassSt45 = 0; fMaxTriggerTracks = 100; fMaxTrackCandidates = 10000; } //_____________________________________________________________________________ void AliMUONRecoParam::SetHighFluxParam() { /// Set reconstruction parameters for high flux environment SetNameTitle("High Flux","High Flux"); SetEventSpecie(AliRecoParam::kHighMult); fMinBendingMomentum = 0.8; fMaxBendingMomentum = 1.e10; fMaxNonBendingSlope = 0.3; fMaxBendingSlope = 0.4; fSelectTrackOnSlope = kFALSE; fNonBendingVertexDispersion = 70.; fBendingVertexDispersion = 70.; fMaxNonBendingDistanceToTrack = 1.; fMaxBendingDistanceToTrack = 1.; fSigmaCutForTracking = 6.; fSigmaCutForImprovement = 5.; fSigmaCutForTrigger = 8.; fStripCutForTrigger = 1.; fMaxStripAreaForTrigger = 3.; fMaxNormChi2MatchTrigger = 16.; fCombinedClusterTrackReco = kFALSE; fTrackAllTracks = kTRUE; fRecoverTracks = kTRUE; fMakeTrackCandidatesFast = kFALSE; fMakeMoreTrackCandidates = kFALSE; fComplementTracks = kTRUE; fImproveTracks = kTRUE; fRemoveConnectedTracksInSt12 = kFALSE; fUseSmoother = kTRUE; for (Int_t iCh = 0; iCh < 10; iCh++) { fUseChamber[iCh] = kTRUE; fDefaultNonBendingReso[iCh] = 0.144; fDefaultBendingReso[iCh] = 0.01; } for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE; fBypassSt45 = 0; fMaxTriggerTracks = 100; fMaxTrackCandidates = 10000; } //_____________________________________________________________________________ void AliMUONRecoParam::SetCosmicParam() { /// Set reconstruction parameters for high flux environment SetNameTitle("Cosmic","Cosmic"); SetEventSpecie(AliRecoParam::kCosmic); fMinBendingMomentum = 0.8; fMaxBendingMomentum = 1.e10; fMaxNonBendingSlope = 0.3; fMaxBendingSlope = 0.4; fSelectTrackOnSlope = kTRUE; fNonBendingVertexDispersion = 170.; fBendingVertexDispersion = 170.; fMaxNonBendingDistanceToTrack = 1.; fMaxBendingDistanceToTrack = 1.; fSigmaCutForTracking = 7.; fSigmaCutForImprovement = 6.; fSigmaCutForTrigger = 8.; fStripCutForTrigger = 1.5; fMaxStripAreaForTrigger = 3.; fMaxNormChi2MatchTrigger = 16.; fPercentOfFullClusterInESD = 100.; fCombinedClusterTrackReco = kFALSE; fTrackAllTracks = kTRUE; fRecoverTracks = kTRUE; fMakeTrackCandidatesFast = kFALSE; fMakeMoreTrackCandidates = kFALSE; fComplementTracks = kTRUE; fImproveTracks = kTRUE; fRemoveConnectedTracksInSt12 = kTRUE; fUseSmoother = kTRUE; fSaveFullClusterInESD = kTRUE; for (Int_t iCh = 0; iCh < 10; iCh++) { fUseChamber[iCh] = kTRUE; fDefaultNonBendingReso[iCh] = 0.4; fDefaultBendingReso[iCh] = 0.4; } fRequestStation[0] = kTRUE; fRequestStation[1] = kTRUE; fRequestStation[2] = kTRUE; fRequestStation[3] = kTRUE; fRequestStation[4] = kTRUE; fBypassSt45 = 0; fPadGoodnessMask = 0x400BE80; // Ped Mean is Zero | Ped Mean Too Low | Ped Mean Too High | Ped Sigma Too Low | Ped Sigma Too High | Ped is missing | HV is missing | manu occupancy too high fMaxTriggerTracks = 100; fMaxTrackCandidates = 10000; SetPedMeanLimits(20, 700); SetManuOccupancyLimits(-1.,0.01); // reject manu above occ=1% SetBuspatchOccupancyLimits(-1,0.01); SetMissingPadFractionLimit(0.1); // 10 % SetFractionOfBuspatchOutsideOccupancyLimit(0.05); // 5 % } //_____________________________________________________________________________ void AliMUONRecoParam::SetCalibrationParam() { /// Set (dummy) reconstruction parameters for calibration environment SetNameTitle("Calibration","Calibration"); SetEventSpecie(AliRecoParam::kCalib); fPedMeanLimits[0] = 5000; fPedMeanLimits[1] = 0; fPadGoodnessMask = 0x8C00; // Pedestal is missing | is too low | too high } //_____________________________________________________________________________ UInt_t AliMUONRecoParam::RequestedStationMask() const { /// Get the mask of the requested station, i.e. an integer where /// bit n is set to one if the station n was requested UInt_t m(0); for ( Int_t i = 0; i < 5; ++i ) { if ( RequestStation(i) ) m |= ( 1 << i ); } return m; } //_____________________________________________________________________________ void AliMUONRecoParam::Print(Option_t *option) const { /// print reconstruction parameters /// if option = FULL then print also unused parameters cout< %5.2f",fMinBendingMomentum)<= %7.2f x ( pedestal sigma ) ",ChargeSigmaCut()) << endl; cout << "Occupancy limits are :" << endl; cout << Form("%e <= Manu occupancy < %7.2f",ManuOccupancyLowLimit(),ManuOccupancyHighLimit()) << endl; cout << Form("%e <= Buspatch occupancy < %7.2f",BuspatchOccupancyLowLimit(),BuspatchOccupancyHighLimit()) << endl; cout << Form("%e <= DE occupancy < %7.2f",DEOccupancyLowLimit(),DEOccupancyHighLimit()) << endl; cout << "'QAChecker' limits" << endl; cout << Form("MissingPadFractionLimit = %5.2f %%",MissingPadFractionLimit()*100.0) << endl; cout << Form("FractionOfBuspatchOutsideOccupancyLimit = %5.2f %%",FractionOfBuspatchOutsideOccupancyLimit()*100.0) << endl; cout << "chamber non bending resolution = |"; for (Int_t iCh = 0; iCh < 10; iCh++) cout << Form(" %6.3f |",fDefaultNonBendingReso[iCh]); cout << endl; cout << "chamber bending resolution = |"; for (Int_t iCh = 0; iCh < 10; iCh++) cout << Form(" %6.3f |",fDefaultBendingReso[iCh]); cout << endl; cout<