]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - MUON/AliMUONReconstructor.cxx
Added protections
[u/mrichter/AliRoot.git] / MUON / AliMUONReconstructor.cxx
index b3225cb3a93f06c69f0667ad1283a5183d87422b..ad85be65841851d567ba68e11b2a115aa8afbfdf 100644 (file)
 ///
 /// Implementation of AliReconstructor for MUON subsystem.
 ///
-/// 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);
+/// The clustering mode and the associated parameters can be changed through the
+/// AliMUONRecoParam object set in the reconstruction macro or read from the CDB
+/// (see methods in AliMUONRecoParam.h file for details)
 ///
 /// Valid modes are :
 ///
 ///
 /// TRIGGERDISABLE : disable the treatment of MUON trigger
 ///
-/// USEFASTTRKDECODER : makes the digit maker class use the high performance decoder
-///                  AliMUONTrackerDDLDecoder instead of AliMUONPayloadTracker.
-///
-/// USEFASTTRGDECODER : makes the digit maker class use the high performance decoder
-///                  AliMUONTriggerDDLDecoder instead of AliMUONPayloadTrigger.
-///
-/// USEFASTDECODERS : makes the digit maker class use the high performance decoders
-///                  AliMUONTrackerDDLDecoder and AliMUONTriggerDDLDecoder.
+/// ENABLEERRORLOGGING : enable error logging (this activates also warnings in RawStreamTracker)
 ///
 /// \author Laurent Aphecetche, Subatech
 //-----------------------------------------------------------------------------
 #include "AliMUONPreClusterFinder.h"
 #include "AliMUONPreClusterFinderV2.h"
 #include "AliMUONPreClusterFinderV3.h"
-#include "AliMUONRecoParam.h"
 #include "AliMUONSimpleClusterServer.h"
 #include "AliMUONTracker.h"
 #include "AliMUONTriggerCircuit.h"
-#include "AliMUONTriggerCrateStore.h"
 #include "AliMUONTriggerStoreV1.h"
 #include "AliMUONVClusterFinder.h"
 #include "AliMUONVClusterServer.h"
 #include "AliMpArea.h"
 #include "AliMpCDB.h"
 #include "AliMpConstants.h"
+#include "AliMpDDLStore.h"
+#include "AliMpSegmentation.h"
 
-#include "AliRecoParam.h"
 #include "AliRawReader.h"
 #include "AliCDBManager.h"
 #include "AliCodeTimer.h"
 ClassImp(AliMUONReconstructor)
 /// \endcond 
 
-AliMUONRecoParam* AliMUONReconstructor::fgRecoParam = 0x0; // reconstruction parameters
-
 //_____________________________________________________________________________
 AliMUONReconstructor::AliMUONReconstructor() : 
 AliReconstructor(),
-fCrateManager(0x0),
 fDigitMaker(0x0),
 fTransformer(new AliMUONGeometryTransformer()),
 fDigitStore(0x0),
@@ -144,10 +132,19 @@ fCalibrationData(0x0),
 fDigitCalibrator(0x0),
 fClusterServer(0x0),
 fTriggerStore(0x0),
-fTrackStore(0x0)
+fTrackStore(0x0),
+fClusterStore(0x0)
 {
   /// normal ctor
 
+  AliDebug(1,"");
+
+  // Unload and delete old mapping
+  AliCDBManager::Instance()->UnloadFromCache("MUON/Calib/Mapping");
+  AliCDBManager::Instance()->UnloadFromCache("MUON/Calib/DDLStore");
+  delete AliMpDDLStore::Instance();
+  delete AliMpSegmentation::Instance();
+
   // Load mapping
   if ( ! AliMpCDB::LoadDDLStore() ) {
     AliFatal("Could not access mapping from OCDB !");
@@ -162,44 +159,22 @@ fTrackStore(0x0)
 AliMUONReconstructor::~AliMUONReconstructor()
 {
   /// dtor
+
+  AliDebug(1,"");
+
   delete fDigitMaker;
   delete fDigitStore;
   delete fTransformer;
-  delete fCrateManager;
   delete fTriggerCircuit;
   delete fCalibrationData;
   delete fDigitCalibrator;
   delete fClusterServer;
   delete fTriggerStore;
   delete fTrackStore;
-}
+  delete fClusterStore;
 
-//_____________________________________________________________________________
-const AliMUONRecoParam* AliMUONReconstructor::GetRecoParam()
-{
-  /// get reconstruction parameters
-  
-  if (!fgRecoParam) {
-    
-    // get reconstruction parameters from AliRecoParam if any
-    TObjArray *recoParams = AliRecoParam::Instance()->GetRecoParam("MUON");
-    
-    if (recoParams) {
-      
-      fgRecoParam = (AliMUONRecoParam*) recoParams->Last();
-      
-    } else {
-      
-      // initialize reconstruction parameters if not already done
-      cout<<"W-AliMUONReconstructor::GetRecoParam: Reconstruction parameters not initialized - Use default one"<<endl;
-      fgRecoParam = AliMUONRecoParam::GetLowFluxParam();
-      AliRecoParam::Instance()->RegisterRecoParam(fgRecoParam);
-      
-    }
-    
-  }
-  
-  return fgRecoParam;
+  delete AliMpSegmentation::Instance(false);
+  delete AliMpDDLStore::Instance(false);
 }
 
 //_____________________________________________________________________________
@@ -261,17 +236,6 @@ AliMUONReconstructor::ConvertDigits(AliRawReader* rawReader, TTree* digitsTree)
   }
 }
 
-//_____________________________________________________________________________
-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
@@ -282,25 +246,15 @@ AliMUONReconstructor::CreateDigitMaker() const
   AliCodeTimerAuto("")
 
   TString option = GetOption();
-  Bool_t enableErrorLogging = kTRUE;
-  Bool_t useFastTrackerDecoder = kFALSE;
-  Bool_t useFastTriggerDecoder = kFALSE;
-  if (option.Contains("USEFASTTRKDECODER"))
-  {
-    useFastTrackerDecoder = kTRUE;
-  }
-  if (option.Contains("USEFASTTRGDECODER"))
-  {
-    useFastTriggerDecoder = kTRUE;
-  }
-  if (option.Contains("USEFASTDECODERS"))
+  
+  Bool_t enableErrorLogging = kFALSE;
+
+  if (option.Contains("ENABLEERRORLOGGING"))
   {
-    useFastTrackerDecoder = kTRUE;
-    useFastTriggerDecoder = kTRUE;
+    enableErrorLogging = kTRUE;
   }
-  fDigitMaker = new AliMUONDigitMaker(
-      enableErrorLogging, useFastTrackerDecoder, useFastTriggerDecoder
-    );
+
+  fDigitMaker = new AliMUONDigitMaker(enableErrorLogging);
   option.ToUpper();
   if ( option.Contains("SAVEDIGITS" ))
     {
@@ -333,9 +287,10 @@ AliMUONReconstructor::CreateTracker() const
 
   AliMUONTracker* tracker(0x0);
   
-  if ( ! AliMUONReconstructor::GetRecoParam()->CombineClusterTrackReco() )
+  if ( ! GetRecoParam()->CombineClusterTrackReco() )
   {
-    tracker = new AliMUONTracker(0x0,
+    tracker = new AliMUONTracker(GetRecoParam(),
+                                0x0,
                                  *DigitStore(),
                                  fDigitMaker,
                                  fTransformer,
@@ -343,7 +298,8 @@ AliMUONReconstructor::CreateTracker() const
   }
   else
   {
-    tracker = new AliMUONTracker(fClusterServer,
+    tracker = new AliMUONTracker(GetRecoParam(),
+                                fClusterServer,
                                  *DigitStore(),
                                  fDigitMaker,
                                  fTransformer,
@@ -482,7 +438,7 @@ AliMUONReconstructor::CreateCalibrator() const
 
   TString calibMode = GetRecoParam()->GetCalibrationMode();
 
-  fDigitCalibrator = new AliMUONDigitCalibrator(*fCalibrationData,calibMode.Data());
+  fDigitCalibrator = new AliMUONDigitCalibrator(*fCalibrationData,GetRecoParam(),calibMode.Data());
 }
 
 //_____________________________________________________________________________
@@ -531,7 +487,7 @@ AliMUONReconstructor::FillTreeR(AliMUONVTriggerStore* triggerStore,
   Bool_t ok(kFALSE);
   Bool_t alone(kTRUE); // is trigger the only info in TreeR ?
   
-  if ( ! AliMUONReconstructor::GetRecoParam()->CombineClusterTrackReco() )
+  if ( ! GetRecoParam()->CombineClusterTrackReco() )
   {
     alone = kFALSE; // we'll get both tracker and trigger information in TreeR
   }
@@ -545,16 +501,17 @@ AliMUONReconstructor::FillTreeR(AliMUONVTriggerStore* triggerStore,
     }
   }
 
-  AliMUONVClusterStore* clusterStore(0x0);
-  
   if ( !alone )
   {
-    clusterStore = new AliMUONClusterStoreV2;
+    if (!fClusterStore)
+    {
+      fClusterStore = new AliMUONClusterStoreV2;
+    }
     
     CreateClusterServer();
     
     TIter next(DigitStore()->CreateIterator());
-    fClusterServer->UseDigits(next);
+    fClusterServer->UseDigits(next,DigitStore());
 
     AliMpArea area;
     
@@ -562,21 +519,22 @@ AliMUONReconstructor::FillTreeR(AliMUONVTriggerStore* triggerStore,
     
     for ( Int_t i = 0; i < AliMpConstants::NofTrackingChambers(); ++i ) 
     {
-      if (AliMUONReconstructor::GetRecoParam()->UseChamber(i))
+      if (GetRecoParam()->UseChamber(i))
       {
-        if ( i >= 6 && AliMUONReconstructor::GetRecoParam()->BypassSt45() ) continue;
+        if ( ( i == 6 || i == 7 )  && GetRecoParam()->BypassSt4() ) continue;
+        if ( ( i == 8 || i == 9 )  && GetRecoParam()->BypassSt5() ) continue;
         
-        fClusterServer->Clusterize(i,*clusterStore,area);
+        fClusterServer->Clusterize(i,*fClusterStore,area,GetRecoParam());
       }
     }
     
-    Bool_t cok = clusterStore->Connect(clustersTree,alone);
+    Bool_t cok = fClusterStore->Connect(clustersTree,alone);
     
     if (!cok) AliError("Could not connect clusterStore to clusterTree");
     
-    AliDebug(1,Form("Number of clusters found = %d",clusterStore->GetSize()));
+    AliDebug(1,Form("Number of clusters found = %d",fClusterStore->GetSize()));
     
-    StdoutToAliDebug(1,clusterStore->Print());
+    StdoutToAliDebug(1,fClusterStore->Print());
   }
          
   if (ok) // at least one type of branches created successfully
@@ -584,7 +542,7 @@ AliMUONReconstructor::FillTreeR(AliMUONVTriggerStore* triggerStore,
     clustersTree.Fill();
   }
   
-  delete clusterStore;
+  if (fClusterStore) fClusterStore->Clear();
 }
 
 //_____________________________________________________________________________