X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=MUON%2FAliMUONReconstructor.cxx;h=a26a8f89bf528cfaae85a3a08852443f85e260c7;hp=ac0815c7c6dfd0996bbc2e762e8bfd1c75c4df72;hb=7e88424f201cac746e949f2b58a4b2ce9c1cfcb6;hpb=88544f7ecb63ae6d29c7e8d792c1412c997678e4 diff --git a/MUON/AliMUONReconstructor.cxx b/MUON/AliMUONReconstructor.cxx index ac0815c7c6d..a26a8f89bf5 100644 --- a/MUON/AliMUONReconstructor.cxx +++ b/MUON/AliMUONReconstructor.cxx @@ -14,26 +14,24 @@ **************************************************************************/ /* $Id$ */ +//----------------------------------------------------------------------------- /// \class AliMUONReconstructor /// /// Implementation of AliReconstructor for MUON subsystem. /// -/// The behavior of the MUON reconstruction can be changed, besides -/// the usual methods found in AliReconstruction (e.g. to disable tracking) -/// by using AliReconstruction::SetOption("MUON",options) -/// where options should be a space separated string. +/// The clustering mode and the associated parameters can be changed by using +/// AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLow(High)FluxParam(); +/// muonRecoParam->Set...(); // see methods in AliMUONRecoParam.h for details +/// AliRecoParam::Instance()->RegisterRecoParam(muonRecoParam); /// -/// Valid options are : -/// -/// SAVEDIGITS : if you want to save in the TreeD the *calibrated* digits -/// that are used for the clustering +/// Valid modes are : /// /// SIMPLEFIT : use the AliMUONClusterFinderSimpleFit clusterizer /// -/// AZ : use the AliMUONClusterFinderAZ clusterizer (default) +/// SIMPLEFITV3 : SIMPLEFIT with preclustering=PRECLUSTERV3 /// -/// MLEM : another implementation of AZ, where preclustering is external -/// MLEMV3 : MLEM with preclustering=PRECLUSTERV2 +/// MLEM : use AliMUONClusterFinderMLEM and AliMUONPreClusterFinder for preclustering (default) +/// MLEMV2 : MLEM with preclustering=PRECLUSTERV2 /// MLEMV3 : MLEM with preclustering=PRECLUSTERV3 /// /// PRECLUSTER : use only AliMUONPreClusterFinder. Only for debug as @@ -46,34 +44,53 @@ /// option either, as center-of-gravity is generally not a good estimate /// of the cluster position... /// +/// PEAKCOG : COG cluster finder around local maxima +/// PEAKFIT : fit around local maxima with up to 3 peaks, COG otherwise +/// /// NOCLUSTERING : bypass completely the clustering stage /// -/// NOSTATUSMAP : disable the computation and usage of the pad status map. Only -/// for debug ! +/// ------ +/// +/// The behavior of the MUON reconstruction can also be changed, besides +/// the usual methods found in AliReconstruction (e.g. to disable tracking) +/// by using AliReconstruction::SetOption("MUON",options) +/// where options should be a space separated string. +/// +/// Valid options are : +/// +/// SAVEDIGITS : if you want to save in the TreeD the *calibrated* digits +/// that are used for the clustering +/// +/// DIGITSTOREV1 : use the V1 implementation of the digitstore +/// DIGITSTOREV2R : use the V2R implementation of the digitstore /// /// NOLOCALRECONSTRUCTION : for debug, to disable local reconstruction (and hence /// "recover" old behavior) /// /// TRIGGERDISABLE : disable the treatment of MUON trigger /// -/// DIGITSTOREV1 : use the V1 implementation of the digitstore -/// DIGITSTOREV2R : use the V2R implementation of the digitstore +/// NOFASTTRKDECODER : makes the digit maker class use the non-high performance decoder +/// AliMUONPayloadTracker instead of AliMUONTrackerDDLDecoder. +/// +/// NOFASTTRGDECODER : makes the digit maker class use the non-high performance decoder +/// AliMUONPayloadTrigger instead of AliMUONTriggerDDLDecoder. +/// +/// NOFASTDECODERS : makes the digit maker class use the non-high performance decoders +/// AliMUONPayloadTracker and AliMUONPayloadTrigger. /// /// \author Laurent Aphecetche, Subatech +//----------------------------------------------------------------------------- #include "AliMUONReconstructor.h" -#include "AliCDBManager.h" -#include "AliLoader.h" -#include "AliLog.h" -#include "AliRunLoader.h" #include "AliMUONCalibrationData.h" #include "AliMUONClusterFinderCOG.h" #include "AliMUONClusterFinderMLEM.h" #include "AliMUONClusterFinderSimpleFit.h" -#include "AliMUONClusterFinderAZ.h" -#include "AliMUONClusterReconstructor.h" +#include "AliMUONClusterFinderPeakCOG.h" +#include "AliMUONClusterFinderPeakFit.h" #include "AliMUONClusterStoreV1.h" +#include "AliMUONClusterStoreV2.h" #include "AliMUONConstants.h" #include "AliMUONDigitCalibrator.h" #include "AliMUONDigitMaker.h" @@ -83,21 +100,29 @@ #include "AliMUONPreClusterFinder.h" #include "AliMUONPreClusterFinderV2.h" #include "AliMUONPreClusterFinderV3.h" +#include "AliMUONSimpleClusterServer.h" #include "AliMUONTracker.h" -#include "AliMUONVTrackStore.h" -#include "AliMUONTriggerChamberEff.h" #include "AliMUONTriggerCircuit.h" -#include "AliMUONTriggerCrateStore.h" #include "AliMUONTriggerStoreV1.h" #include "AliMUONVClusterFinder.h" +#include "AliMUONVClusterServer.h" +#include "AliMUONVTrackStore.h" + +#include "AliMpArea.h" #include "AliMpCDB.h" +#include "AliMpConstants.h" + #include "AliRawReader.h" +#include "AliCDBManager.h" #include "AliCodeTimer.h" +#include "AliLog.h" + #include +#include #include #include #include -//#include "AliCodeTimer.h" + /// \cond CLASSIMP ClassImp(AliMUONReconstructor) /// \endcond @@ -105,36 +130,26 @@ ClassImp(AliMUONReconstructor) //_____________________________________________________________________________ AliMUONReconstructor::AliMUONReconstructor() : AliReconstructor(), -fCrateManager(0x0), fDigitMaker(0x0), fTransformer(new AliMUONGeometryTransformer()), fDigitStore(0x0), fTriggerCircuit(0x0), fCalibrationData(0x0), fDigitCalibrator(0x0), -fClusterReconstructor(0x0), -fClusterStore(0x0), +fClusterServer(0x0), fTriggerStore(0x0), fTrackStore(0x0), -fTrigChamberEff(0x0) +fClusterStore(0x0) { /// normal ctor - // Load geometry data - fTransformer->LoadGeometryData(); - // Load mapping - if ( ! AliMpCDB::LoadMpSegmentation() ) - { + if ( ! AliMpCDB::LoadDDLStore() ) { AliFatal("Could not access mapping from OCDB !"); } - // Load DDL store - if ( ! AliMpCDB::LoadDDLStore() ) - { - AliFatal("Could not access DDL Store from OCDB !"); - } - + // Load geometry data + fTransformer->LoadGeometryData(); } @@ -145,15 +160,13 @@ AliMUONReconstructor::~AliMUONReconstructor() delete fDigitMaker; delete fDigitStore; delete fTransformer; - delete fCrateManager; delete fTriggerCircuit; delete fCalibrationData; delete fDigitCalibrator; - delete fClusterReconstructor; - delete fClusterStore; + delete fClusterServer; delete fTriggerStore; delete fTrackStore; - delete fTrigChamberEff; + delete fClusterStore; } //_____________________________________________________________________________ @@ -169,39 +182,6 @@ AliMUONReconstructor::Calibrate(AliMUONVDigitStore& digitStore) const fDigitCalibrator->Calibrate(digitStore); } -//_____________________________________________________________________________ -void -AliMUONReconstructor::Clusterize(const AliMUONVDigitStore& digitStore, - AliMUONVClusterStore& clusterStore) const -{ - /// Creates clusters from digits. - - TString sopt(GetOption()); - sopt.ToUpper(); - if ( sopt.Contains("NOCLUSTERING") ) return; - - if (!fClusterReconstructor) - { - CreateClusterReconstructor(); - } - - AliCodeTimerAuto(Form("%s::Digits2Clusters(const AliMUONVDigitStore&,AliMUONVClusterStore&)", - fClusterReconstructor->ClassName())) - fClusterReconstructor->Digits2Clusters(digitStore,clusterStore); -} - -//_____________________________________________________________________________ -AliMUONVClusterStore* -AliMUONReconstructor::ClusterStore() const -{ - /// Return (and create if necessary) the cluster container - if (!fClusterStore) - { - fClusterStore = new AliMUONClusterStoreV1; - } - return fClusterStore; -} - //_____________________________________________________________________________ void AliMUONReconstructor::ConvertDigits(AliRawReader* rawReader, @@ -224,6 +204,7 @@ void AliMUONReconstructor::ConvertDigits(AliRawReader* rawReader, TTree* digitsTree) const { /// convert raw data into a digit tree + AliCodeTimerAuto("") Bool_t alone = ( TriggerStore() == 0 ); @@ -240,22 +221,13 @@ AliMUONReconstructor::ConvertDigits(AliRawReader* rawReader, TTree* digitsTree) else { ConvertDigits(rawReader,DigitStore(),TriggerStore()); + AliCodeTimerStart("Fill digits") digitsTree->Fill(); + AliCodeTimerStop("Fill digits") DigitStore()->Clear(); } } -//_____________________________________________________________________________ -AliMUONTriggerCrateStore* -AliMUONReconstructor::CrateManager() const -{ - /// Return (and create if necessary) the trigger crate store - if (fCrateManager) return fCrateManager; - fCrateManager = new AliMUONTriggerCrateStore; - fCrateManager->ReadFromFile(); - return fCrateManager; -} - //_____________________________________________________________________________ void AliMUONReconstructor::CreateDigitMaker() const @@ -265,7 +237,31 @@ AliMUONReconstructor::CreateDigitMaker() const AliCodeTimerAuto("") - fDigitMaker = new AliMUONDigitMaker; + TString option = GetOption(); + Bool_t enableErrorLogging = kTRUE; + Bool_t useFastTrackerDecoder = kTRUE; + Bool_t useFastTriggerDecoder = kTRUE; + if (option.Contains("NOFASTTRKDECODER")) + { + useFastTrackerDecoder = kFALSE; + } + if (option.Contains("NOFASTTRGDECODER")) + { + useFastTriggerDecoder = kFALSE; + } + if (option.Contains("NOFASTDECODERS")) + { + useFastTrackerDecoder = kFALSE; + useFastTriggerDecoder = kFALSE; + } + fDigitMaker = new AliMUONDigitMaker( + enableErrorLogging, useFastTrackerDecoder, useFastTriggerDecoder + ); + option.ToUpper(); + if ( option.Contains("SAVEDIGITS" )) + { + fDigitMaker->SetMakeTriggerDigits(kTRUE); + } } //_____________________________________________________________________________ @@ -281,55 +277,50 @@ AliMUONReconstructor::CreateTriggerCircuit() const } -//_____________________________________________________________________________ -void -AliMUONReconstructor::CreateTriggerChamberEff() const -{ - /// Create (and create if necessary) the trigger chamber efficiency class - if (fTrigChamberEff) return; - - AliCodeTimerAuto("") - - fTrigChamberEff = new AliMUONTriggerChamberEff(fTransformer,fDigitMaker,kTRUE); - //fTrigChamberEff->SetDebugLevel(1); -} - //_____________________________________________________________________________ AliTracker* -AliMUONReconstructor::CreateTracker(AliRunLoader* runLoader) const +AliMUONReconstructor::CreateTracker() const { /// Create the MUONTracker object - /// The MUONTracker is passed the GetOption(), i.e. our own options CreateTriggerCircuit(); CreateDigitMaker(); - CreateTriggerChamberEff(); + CreateClusterServer(); + + AliMUONTracker* tracker(0x0); - AliLoader* loader = runLoader->GetDetectorLoader("MUON"); - if (!loader) + if ( ! AliMUONReconstructor::GetRecoParam()->CombineClusterTrackReco() ) { - AliError("Cannot get MUONLoader, so cannot create MUONTracker"); - return 0x0; + tracker = new AliMUONTracker(0x0, + *DigitStore(), + fDigitMaker, + fTransformer, + fTriggerCircuit); + } + else + { + tracker = new AliMUONTracker(fClusterServer, + *DigitStore(), + fDigitMaker, + fTransformer, + fTriggerCircuit); } - AliMUONTracker* tracker = new AliMUONTracker(loader,fDigitMaker,fTransformer,fTriggerCircuit,fTrigChamberEff); - tracker->SetOption(GetOption()); + return tracker; } //_____________________________________________________________________________ -void -AliMUONReconstructor::CreateClusterReconstructor() const +AliMUONVClusterFinder* +AliMUONReconstructor::CreateClusterFinder(const char* clusterFinderType) { - /// Create cluster reconstructor, depending on GetOption() + /// Create a given cluster finder instance - AliCodeTimerAuto("") + AliCodeTimerAutoGeneral("") - AliDebug(1,""); - AliMUONVClusterFinder* clusterFinder(0x0); - TString opt(GetOption()); + TString opt(clusterFinderType); opt.ToUpper(); if ( strstr(opt,"PRECLUSTERV2") ) @@ -344,6 +335,14 @@ AliMUONReconstructor::CreateClusterReconstructor() const { clusterFinder = new AliMUONPreClusterFinder; } + else if ( strstr(opt,"PEAKCOG") ) + { + clusterFinder = new AliMUONClusterFinderPeakCOG(kFALSE,new AliMUONPreClusterFinder); + } + else if ( strstr(opt,"PEAKFIT") ) + { + clusterFinder = new AliMUONClusterFinderPeakFit(kFALSE,new AliMUONPreClusterFinder); + } else if ( strstr(opt,"COG") ) { clusterFinder = new AliMUONClusterFinderCOG(new AliMUONPreClusterFinder); @@ -372,22 +371,32 @@ AliMUONReconstructor::CreateClusterReconstructor() const { clusterFinder = new AliMUONClusterFinderMLEM(kFALSE,new AliMUONPreClusterFinder); } - else if ( strstr(opt,"AZ") ) - { - clusterFinder = new AliMUONClusterFinderAZ; - } else { - // default is currently AZ - clusterFinder = new AliMUONClusterFinderAZ; + AliErrorClass(Form("clustering mode \"%s\" does not exist",opt.Data())); + return 0x0; } - if ( clusterFinder ) - { - AliInfo(Form("Will use %s for clusterizing",clusterFinder->ClassName())); - } + return clusterFinder; +} + +//_____________________________________________________________________________ +void +AliMUONReconstructor::CreateClusterServer() const +{ + /// Create cluster server + + if ( fClusterServer ) return; + + AliCodeTimerAuto(""); + + AliMUONVClusterFinder* clusterFinder = CreateClusterFinder(GetRecoParam()->GetClusteringMode()); + + if ( !clusterFinder ) return; - fClusterReconstructor = new AliMUONClusterReconstructor(clusterFinder,fTransformer); + AliInfo(Form("Will use %s for clusterizing",clusterFinder->ClassName())); + + fClusterServer = new AliMUONSimpleClusterServer(clusterFinder,*fTransformer); } //_____________________________________________________________________________ @@ -421,14 +430,15 @@ AliMUONReconstructor::CreateCalibrator() const TString opt(GetOption()); opt.ToUpper(); - Bool_t statusMap(kTRUE); if ( strstr(opt,"NOSTATUSMAP") ) { - AliWarning("Disconnecting status map : SHOULD BE USED FOR DEBUG ONLY. NOT FOR PRODUCTION !!!"); - statusMap = kFALSE; + AliWarning("NOSTATUSMAP is obsolete"); } - fDigitCalibrator = new AliMUONDigitCalibrator(*fCalibrationData,statusMap); + + TString calibMode = GetRecoParam()->GetCalibrationMode(); + + fDigitCalibrator = new AliMUONDigitCalibrator(*fCalibrationData,GetRecoParam(),calibMode.Data()); } //_____________________________________________________________________________ @@ -466,7 +476,6 @@ AliMUONReconstructor::DigitStore() const //_____________________________________________________________________________ void AliMUONReconstructor::FillTreeR(AliMUONVTriggerStore* triggerStore, - AliMUONVClusterStore* clusterStore, TTree& clustersTree) const { /// Write the trigger and cluster information into TreeR @@ -476,30 +485,64 @@ AliMUONReconstructor::FillTreeR(AliMUONVTriggerStore* triggerStore, AliDebug(1,""); Bool_t ok(kFALSE); + Bool_t alone(kTRUE); // is trigger the only info in TreeR ? + + if ( ! AliMUONReconstructor::GetRecoParam()->CombineClusterTrackReco() ) + { + alone = kFALSE; // we'll get both tracker and trigger information in TreeR + } + if ( triggerStore ) { - Bool_t alone = ( clusterStore ? kFALSE : kTRUE ); ok = triggerStore->Connect(clustersTree,alone); if (!ok) { AliError("Could not create triggerStore branches in TreeR"); } } - - if ( clusterStore ) + + if ( !alone ) { - Bool_t alone = ( triggerStore ? kFALSE : kTRUE ); - ok = clusterStore->Connect(clustersTree,alone); - if (!ok) + if (!fClusterStore) { - AliError("Could not create triggerStore branches in TreeR"); - } + fClusterStore = new AliMUONClusterStoreV2; + } + + CreateClusterServer(); + + TIter next(DigitStore()->CreateIterator()); + fClusterServer->UseDigits(next); + + AliMpArea area; + + AliDebug(1,Form("Doing full clusterization in local reconstruction using %s ",fClusterServer->ClassName())); + + for ( Int_t i = 0; i < AliMpConstants::NofTrackingChambers(); ++i ) + { + if (AliMUONReconstructor::GetRecoParam()->UseChamber(i)) + { + if ( ( i == 6 || i == 7 ) && AliMUONReconstructor::GetRecoParam()->BypassSt4() ) continue; + if ( ( i == 8 || i == 9 ) && AliMUONReconstructor::GetRecoParam()->BypassSt5() ) continue; + + fClusterServer->Clusterize(i,*fClusterStore,area); + } + } + + Bool_t cok = fClusterStore->Connect(clustersTree,alone); + + if (!cok) AliError("Could not connect clusterStore to clusterTree"); + + AliDebug(1,Form("Number of clusters found = %d",fClusterStore->GetSize())); + + StdoutToAliDebug(1,fClusterStore->Print()); } - + if (ok) // at least one type of branches created successfully { clustersTree.Fill(); } + + if (fClusterStore) fClusterStore->Clear(); } //_____________________________________________________________________________ @@ -521,23 +564,6 @@ AliMUONReconstructor::HasDigitConversion() const } } -//_____________________________________________________________________________ -Bool_t -AliMUONReconstructor::HasLocalReconstruction() const -{ - /// Whether or not we have local reconstruction - TString opt(GetOption()); - opt.ToUpper(); - if ( opt.Contains("NOLOCALRECONSTRUCTION" ) ) - { - return kFALSE; - } - else - { - return kTRUE; - } -} - //_____________________________________________________________________________ void AliMUONReconstructor::Reconstruct(AliRawReader* rawReader, TTree* clustersTree) const @@ -552,87 +578,8 @@ AliMUONReconstructor::Reconstruct(AliRawReader* rawReader, TTree* clustersTree) } ConvertDigits(rawReader,DigitStore(),TriggerStore()); - Clusterize(*(DigitStore()),*(ClusterStore())); - - FillTreeR(TriggerStore(),ClusterStore(),*clustersTree); -} -//_____________________________________________________________________________ -void -AliMUONReconstructor::Reconstruct(AliRunLoader* runLoader) const -{ - /// Reconstruct simulated data - - AliCodeTimerAuto("Reconstruct(AliRunLoader*)") - - AliLoader* loader = runLoader->GetDetectorLoader("MUON"); - if (!loader) - { - AliError("Could not get MUON loader"); - return; - } - - Int_t nEvents = runLoader->GetNumberOfEvents(); - - for ( Int_t i = 0; i < nEvents; ++i ) - { - runLoader->GetEvent(i); - - loader->LoadRecPoints("update"); - loader->CleanRecPoints(); - loader->MakeRecPointsContainer(); - TTree* clustersTree = loader->TreeR(); - - loader->LoadDigits("read"); - TTree* digitsTree = loader->TreeD(); - - Reconstruct(digitsTree,clustersTree); - - loader->UnloadDigits(); - loader->WriteRecPoints("OVERWRITE"); - loader->UnloadRecPoints(); - } -} - -//_____________________________________________________________________________ -void -AliMUONReconstructor::Reconstruct(AliRunLoader* runLoader, AliRawReader* rawReader) const -{ - /// This method is called by AliReconstruction if HasLocalReconstruction()==kFALSE - - AliCodeTimerAuto("AliMUONReconstructor::Reconstruct(AliRunLoader*, AliRawReader*)") - - AliLoader* loader = runLoader->GetDetectorLoader("MUON"); - if (!loader) - { - AliError("Could not get MUON loader"); - return; - } - - Int_t i(0); - - while (rawReader->NextEvent()) - { - runLoader->GetEvent(i++); - - loader->LoadRecPoints("update"); - loader->CleanRecPoints(); - loader->MakeRecPointsContainer(); - TTree* clustersTree = loader->TreeR(); - - loader->LoadDigits("update"); - loader->CleanDigits(); - loader->MakeDigitsContainer(); - TTree* digitsTree = loader->TreeD(); - ConvertDigits(rawReader, digitsTree); - loader->WriteDigits("OVERWRITE"); - - Reconstruct(digitsTree,clustersTree); - - loader->UnloadDigits(); - loader->WriteRecPoints("OVERWRITE"); - loader->UnloadRecPoints(); - } + FillTreeR(TriggerStore(),*clustersTree); } //_____________________________________________________________________________ @@ -642,7 +589,7 @@ AliMUONReconstructor::Reconstruct(TTree* digitsTree, TTree* clustersTree) const /// This method is called by AliReconstruction if HasLocalReconstruction()==kTRUE /// AND HasDigitConversion()==kTRUE -// AliCodeTimerAuto("(TTree*,TTree*)") + AliCodeTimerAuto("") AliDebug(1,""); @@ -716,14 +663,13 @@ AliMUONReconstructor::Reconstruct(TTree* digitsTree, TTree* clustersTree) const // i.e. w/o going through raw data, this will be the case) TIter next(fDigitStore->CreateIterator()); AliMUONVDigit* digit = static_cast(next()); - if (!digit->IsCalibrated()) + if (digit && !digit->IsCalibrated()) { Calibrate(*fDigitStore); } - Clusterize(*fDigitStore,*(ClusterStore())); } - FillTreeR(fTriggerStore,ClusterStore(),*clustersTree); + FillTreeR(fTriggerStore,*clustersTree); } //_____________________________________________________________________________