Load pythia libraries.

[u/mrichter/AliRoot.git] / MUON / runReconstruction.C

diff --git a/MUON/runReconstruction.C b/MUON/runReconstruction.C
index 9b3b27f04dd0ab3f03bb415710eda01a6707ad20..b4637b6e7d4720723f8f46653bbf5ad2c65b42a8 100644 (file)
--- a/MUON/runReconstruction.C
+++ b/MUON/runReconstruction.C
@@ -15,23 +15,36 @@
 
 /* $Id$ */
 
-// Macro extracted from MUON test script
-// By Laurent Aphecetche
+/// \ingroup macros
+/// \file runReconstruction.C
+/// \brief Macro for running reconstruction
+///
+/// Macro extracted from the MUON test script
+///
+/// \author Laurent Aphecetche
 
 #if !defined(__CINT__) || defined(__MAKECINT__)
+#include "AliMUONReconstructor.h"
+#include "AliMUONRecoParam.h"
+#include "AliRecoParam.h"
 #include "AliCDBManager.h"
 #include "AliMagFMaps.h"
 #include "AliTracker.h"
 #include "AliReconstruction.h"
 #include <TRandom.h>
+//#include <TObjectTable.h>
 #endif
 
-void runReconstruction(int run, int seed, const char* input, const char* recoptions)
+void runReconstruction(int seed, const char* input, const char* recoptions)
 { 
-  AliCDBManager::Instance()->SetRun(run);
+  AliCDBManager* man = AliCDBManager::Instance();
+  man->SetDefaultStorage("local://$ALICE_ROOT");
+  
   gRandom->SetSeed(seed);
+  
   AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 1, 1., 10., AliMagFMaps::k5kG);
   AliTracker::SetFieldMap(field, kFALSE);
+  
   AliReconstruction* MuonRec = new AliReconstruction("galice.root");
   MuonRec->SetInput(input);
   MuonRec->SetRunVertexFinder(kFALSE);
@@ -42,8 +55,19 @@ void runReconstruction(int run, int seed, const char* input, const char* recopti
   MuonRec->SetNumberOfEventsPerFile(1000);
   MuonRec->SetOption("MUON",recoptions);
   //  MuonRec->SetEventRange(319,319);
+  MuonRec->SetWriteAOD();
+  
+  AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam();
+  muonRecoParam->CombineClusterTrackReco(kTRUE);
+  //muonRecoParam->SetClusteringMode("PEAKFIT");
+  //muonRecoParam->SetClusteringMode("PEAKCOG");
+  muonRecoParam->Print("FULL");
+  AliRecoParam::Instance()->RegisterRecoParam(muonRecoParam);
+  
   MuonRec->Run();
+  
   delete MuonRec;
+  
   //gObjectTable->Print();
 }