update test macros
authorjthaeder <jthaeder@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 27 Aug 2010 14:17:06 +0000 (14:17 +0000)
committerjthaeder <jthaeder@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 27 Aug 2010 14:17:06 +0000 (14:17 +0000)
HLT/JET/macros/HLTJetReconstruction.C
HLT/JET/macros/readJets.C
HLT/JET/macros/runChain.sh
HLT/JET/macros/runTask.sh
HLT/JET/macros/tasks/ConfigJetAnalysisHLT.C
HLT/JET/macros/tasks/ConfigJetAnalysisHLTMC.C
HLT/JET/macros/tasks/JetAnalysisManagerHLT.C
HLT/JET/macros/tasks/JetAnalysisManagerHLTKine.C [deleted file]
HLT/JET/macros/tasks/JetAnalysisManagerHLTMC.C [new file with mode: 0644]

index 4d6cc10..7b7a85a 100644 (file)
@@ -107,8 +107,9 @@ void HLTJetReconstruction(Int_t nEvents=1, Bool_t generate=kFALSE, PprRun_t runT
     // ------------------------------------------
     TString publisherId("ESDMCEventPublisher");
     // ------------------------------------------
-
-    TString publisherArg( Form("-entrytype MCFAST -dataspec 0x0000001F -datapath /home/jthaeder/jet/data/HEAD_2009-10-26/FastGen/kPythia6Jets125_150_14TeV/JET-ETA=-0.2,0.2_JET-ET=50,1000_R=0.4_200ev") );
+    // MCFAST - MC - ESD
+    TString publisherArg("-entrytype ESD -dataspec 0x0000001F -datapath /lustre/alice/jthaeder/data/HEAD_2010-07-09/7TeV/pp_Perugia0/014000");
+    
     AliHLTConfiguration ESDMCEventPublisher(publisherId.Data(), "ESDMCEventPublisher", NULL, publisherArg.Data() );
 
     if (!analysisInput.IsNull()) analysisInput+=" ";
@@ -142,15 +143,18 @@ void HLTJetReconstruction(Int_t nEvents=1, Bool_t generate=kFALSE, PprRun_t runT
   // -- Processing Components -- //
   // -                         - //
   // ----------------------------//
+
 #if 1
+
   // ------------------------------------------
   // -- ConeJetFinder
   // ------------------------------------------
   TString jetId("JETConeJet");
   // ------------------------------------------
 
-  TString jetArg( Form("-algorithm FSCSquareCell -leading 1 -coneRadius %.1f -trackCutMinPt 0.0 -seedCutMinPt %.1f -jetCutMinEt %.1f",
-                      coneRadius, cutPtSeed, cutEtJet) );
+  TString jetArg(
+                Form("-algorithm FSCSquareCell -leading 1 -coneRadius %.1f -trackCutMinPt 0.0 -seedCutMinPt %.1f -jetCutMinEt %.1f", coneRadius, cutPtSeed, cutEtJet) 
+                );
   
   AliHLTConfiguration jetCone(jetId.Data(), "JETConeJetFinder", jetInput.Data(), jetArg.Data()); 
   
@@ -158,24 +162,22 @@ void HLTJetReconstruction(Int_t nEvents=1, Bool_t generate=kFALSE, PprRun_t runT
   analysisInput += jetId;
 
 #else
+
   // ------------------------------------------
   // -- FastJetFinder
   // ------------------------------------------
 
-  AliHLTConfiguration jetFinder("JETFastJet_Kt", "JETFastJetFinder",
-                               jetInput.Data(),"-finderType kt");
+  AliHLTConfiguration jetFinder("JETFastJet_Kt", "JETFastJetFinder", jetInput.Data(),"-finderType kt");
     
   if (!writerInput.IsNull()) writerInput+=" ";
   writerInput+="JETFastJet_Kt";
   
-  
-  AliHLTConfiguration jetFinder("JETFastJet_AntiKt", "JETFastJetFinder",
-                               jetInput.Data(),"-finderType antikt");
+  AliHLTConfiguration jetFinder("JETFastJet_AntiKt", "JETFastJetFinder", jetInput.Data(),"-finderType antikt");
   
   if (!writerInput.IsNull()) writerInput+=" ";
   writerInput+="JETFastJet_AntiKt";
 #endif
-  
+
   // ------------------------------------------
   // -- Jet Analysis 
   // ------------------------------------------
@@ -198,8 +200,7 @@ void HLTJetReconstruction(Int_t nEvents=1, Bool_t generate=kFALSE, PprRun_t runT
   TString writerArg( Form("-directory analysis -datafile analyze_%d_%s -write-all-events", nEvents, pprRunName[runType] ));
 
   // -- The RootFileWriter 
-  AliHLTConfiguration rootWriter("RootWriter", "ROOTFileWriter", 
-                                writerInput.Data(), writerArg.Data() );
+  AliHLTConfiguration rootWriter("RootWriter", "ROOTFileWriter", writerInput.Data(), writerArg.Data() );
   
   // --------------------------- //
   // -                         - //
index 676455c..0a8f3f0 100644 (file)
@@ -28,7 +28,7 @@ void readJets( Char_t* fileName = "./analysis/EOR_analyze_100000_kPythia6Jets104
   
   AliHLTJETAnalysisMerge* merge = new AliHLTJETAnalysisMerge();
   merge->Initialize();
-
+    
   merge->AddJets( jets ); 
     
   merge->CreateCanvas();
index 3f8dc04..6728cec 100755 (executable)
@@ -1,5 +1,8 @@
 #/bin/bash
 
+NEVENTS=200
+Pt_BIN=kPythia6Jets104_125
+
 if [ ! -d ./analysis ] ; then
     mkdir analysis
 fi
@@ -9,8 +12,9 @@ ROOTFILES=`find . -maxdepth 1 -name "*.root" `
 if [ -n "${ROOTFILES}" ] ; then
     rm ${ROOTFILES}
 fi
 
-aliroot -b -l -q 'HLTJetReconstruction.C(20,kTRUE)' 2>&1 | tee ChainLog.log
+aliroot -b -l -q HLTJetReconstruction.C'('${NEVENTS}', kFALSE, '${Pt_BIN}')' 2>&1 | tee ChainLog.log
 
 ROOTFILES=`find . -maxdepth 1 -name "*.root" `
 
@@ -18,6 +22,6 @@ if [ -n "${ROOTFILES}" ] ; then
     rm ${ROOTFILES}
 fi
 
-aliroot -l 'readJets.C("./analysis/EOR_analyze_20_kPythia6Jets104_125.root")'
+aliroot -l -b -q readJets.C'("./analysis/EOR_analyze_'${NEVENTS}'_'${Pt_BIN}'.root")'
 
-#valgrind --error-limit=no --leak-check=full --show-reachable=yes aliroot -b -l -q 'HLTJetReconstruction.C(10)' 
+# valgrind --error-limit=no --leak-check=full --show-reachable=yes aliroot -b -l -q 'HLTJetReconstruction.C(10)' 
index 525dbd1..e8748b9 100755 (executable)
@@ -1,7 +1,11 @@
 #/bin/bash
 
-if [ ! -d ./analysis ] ; then
-    mkdir analysis
+if [ "$1" == "MC" ] ; then 
+    MODE=MC 
+elif [ "$1" == "ESD" ] ; then 
+    MODE=ESD
+else    
+    MODE=MC
 fi
 
 ROOTFILES=`find . -maxdepth 1 -name "*.root" `
@@ -10,7 +14,13 @@ if [ -n "${ROOTFILES}" ] ; then
     rm ${ROOTFILES}
 fi
 
-aliroot -b -l -q './tasks/JetAnalysisManagerHLT.C' 2>&1 | tee TaskLog.log
+if [ "$MODE" == "ESD" ] ; then
+    echo " -= Load : JetAnalysisManagerHLT.C =- "
+    aliroot -b -l -q './tasks/JetAnalysisManagerHLT.C' 2>&1 | tee TaskLog.log
+else
+    echo " -= Load : JetAnalysisManagerHLTMC.C =- "
+    aliroot -b -l -q './tasks/JetAnalysisManagerHLTMC.C' 2>&1 | tee TaskLog.log
+fi
 
 ROOTFILES=`find . -maxdepth 1 -name "*.root" `
 
index 9964b90..466c286 100644 (file)
@@ -1,6 +1,6 @@
 /** ---------------------------------------------------------------------
  *  @file   ConfigJetAnalysisHLT.C
- *  @author Jochen Thaeder <thaeder@kip.uni-heidelberg.de>
+ *  @author Jochen Thaeder <jochen@thaeder.de>
  *  @brief  Run HLT cone finder in analysis framework, on ESDs
  *
  *  --------------------------------------------------------------------- 
@@ -48,8 +48,8 @@ AliJetFinder*  ConfigJetAnalysis() {
   jetReaderHeader->SetJetAlgorithm(algorithm);
 
   // Set prt to track cuts
-  jetReaderHeader->SetTrackCuts( trackCuts );
-  jetReaderHeader->SetSeedCuts( seedCuts );
+  jetReaderHeader->SetTrackCuts(trackCuts);
+  jetReaderHeader->SetSeedCuts(seedCuts);
 
   // Set Eta min/max and Phi min/max
   jetReaderHeader->SetFiducialEta( -0.9, 0.9) ;
index 38a7fe9..a65a4c3 100644 (file)
@@ -1,6 +1,6 @@
 /** ---------------------------------------------------------------------
  *  @file   ConfigJetAnalysisHLTMC.C
- *  @author Jochen Thaeder <thaeder@kip.uni-heidelberg.de>
+ *  @author Jochen Thaeder <jochen@thaeder.de>
  *  @brief  Run HLT cone finder in analysis framework, on Kinematics
  *
  *  --------------------------------------------------------------------- 
@@ -48,8 +48,8 @@ AliJetFinder*  ConfigJetAnalysis() {
   jetReaderHeader->SetJetAlgorithm(algorithm);
 
   // Set prt to track cuts
-  jetReaderHeader->SetTrackCuts( trackCuts );
-  jetReaderHeader->SetSeedCuts( seedCuts );
+  jetReaderHeader->SetTrackCuts(trackCuts);
+  jetReaderHeader->SetSeedCuts(seedCuts);
 
   // Set Eta min/max and Phi min/max
   jetReaderHeader->SetFiducialEta( -0.9, 0.9) ;
index 9f7d405..58cb89a 100644 (file)
@@ -31,10 +31,13 @@ void JetAnalysisManagerHLT() {
   //
   // --------------------------------------------------------------------------------  
 
-  gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C");
+  gROOT->LoadMacro("${ALICE_ROOT}/PWG0/CreateESDChain.C");
   TChain* chain = new TChain("esdTree");
   
-  chain->Add("~/jet/data/HEAD_2009-06-04/kPythia6Jets104_125_14TeV/JET-ETA=-0.2,0.2_JET-ET=10,1000_R=0.7_100ev/AliESDs.root");
+  chain->Add("~/jet/data/HEAD_2010-01-08/Gen/kPythia6Jets86_104_14TeV/JET-ETA=-0.2,0.2_JET-ET=10,1000_R=0.4_20ev/AliESDs.root");
+
+
+  //  chain->Add("~/jet/data/HEAD_2009-06-04/kPythia6Jets104_125_14TeV/JET-ETA=-0.2,0.2_JET-ET=10,1000_R=0.7_100ev/AliESDs.root");
   //chain->Add("~/jet/data/HEAD_2009-06-04/kPythia6Jets104_125_14TeV/JET-ETA=-0.2,0.2_JET-ET=10,1000_R=0.7_1ev/AliESDs.root");
 
   // --------------------------------------------------------------------------------  
@@ -52,12 +55,12 @@ void JetAnalysisManagerHLT() {
 
   // MC Truth
   AliMCEventHandler* mcHandler = new AliMCEventHandler();
-    
+
   AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
   mgr->SetInputEventHandler  (inpHandler);
   mgr->SetOutputEventHandler (aodHandler);
   mgr->SetMCtruthEventHandler(mcHandler);
-  mgr->SetDebugLevel(0);
+  mgr->SetDebugLevel(10);
 
   // --------------------------------------------------------------------------------
   //
@@ -70,14 +73,7 @@ void JetAnalysisManagerHLT() {
   jetana->SetConfigFile("./tasks/ConfigJetAnalysisHLT.C");
   jetana->SetNonStdBranch("jetsHLT");
   jetana->SetDebugLevel(10);
-  mgr->AddTask(jetana);  
-
-  // -- HLT FFSC - MC
-  AliAnalysisTaskJets *jetanaMC = new AliAnalysisTaskJets("JetAnalysisHLTMC");
-  jetanaMC->SetConfigFile("./tasks/ConfigJetAnalysisHLTMC.C");
-  jetanaMC->SetNonStdBranch("jetsHLTMC");
-  jetanaMC->SetDebugLevel(0);
-  mgr->AddTask(jetanaMC);
+  mgr->AddTask(jetana);
 
   // --------------------------------------------------------------------------------
   //
@@ -88,22 +84,18 @@ void JetAnalysisManagerHLT() {
   AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
   AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
   AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(), AliAnalysisManager::kOutputContainer, "histos.root");
-  AliAnalysisDataContainer *coutputMC2 = mgr->CreateContainer("histosMC", TList::Class(), AliAnalysisManager::kOutputContainer, "histosMC.root");
-  
+
   mgr->ConnectInput  (jetana,     0, cinput1  );
   mgr->ConnectOutput (jetana,     0, coutput1 );
   mgr->ConnectOutput (jetana,     1, coutput2 );
-    
-  mgr->ConnectInput  (jetanaMC,     0, cinput1  );
-  mgr->ConnectOutput (jetanaMC,     0, coutput1 );
-  mgr->ConnectOutput (jetanaMC,     1, coutputMC2 );
-  
+
+
   // --------------------------------------------------------------------------------  
   //
   // Run the analysis
   //    
   // --------------------------------------------------------------------------------  
-  
+
   mgr->InitAnalysis();
   mgr->PrintStatus();
   mgr->StartAnalysis("local",chain);
diff --git a/HLT/JET/macros/tasks/JetAnalysisManagerHLTKine.C b/HLT/JET/macros/tasks/JetAnalysisManagerHLTKine.C
deleted file mode 100644 (file)
index 9f4a76c..0000000
+++ /dev/null
@@ -1,76 +0,0 @@
-void JetAnalysisManagerHLTKine()
-{
-      gSystem->Load("libCGAL.so");
-      gSystem->Load("${FASTJET}/lib/libfastjet.so");
-
-      gSystem->Load("libTree.so");
-      gSystem->Load("libPhysics.so");
-      gSystem->Load("libGeom.so");
-      gSystem->Load("libVMC.so");
-      gSystem->Load("libANALYSIS.so");
-      gSystem->Load("libSTEERBase.so");
-      gSystem->Load("libAOD.so");
-      gSystem->Load("libESD.so");
-      gSystem->Load("libANALYSISalice.so");
-      gSystem->Load("libJETAN.so");
-
-     //
-    if (gApplication) gApplication->InitializeGraphics();
-    // Create the chain
-    //
-
-    gROOT->LoadMacro("${ALICE_ROOT}/PWG0/CreateESDChain.C");
-    TChain* chain = new TChain("TE");
-    chain->Add("/home/jthaeder/jet/data/HEAD_2009-03-17/FastGen/kPythia6Jets104_125_14TeV/JET-ETA=-0.2,0.2_JET-ET=50,1000_R=0.4_10ev/galice.root");
-
-    /////////////////////////////////////////////////////////////////////////////////// 
-    // Create the analysis manager
-    //
-    // Output
-    AliAODHandler* aodHandler = new AliAODHandler();
-    aodHandler->SetOutputFileName("aod.root");
-    // MC Truth
-    AliMCEventHandler* mcHandler = new AliMCEventHandler();
-    mcHandler->SetReadTR(kFALSE);
-
-    AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
-    mgr->SetOutputEventHandler (aodHandler);
-    mgr->SetMCtruthEventHandler(mcHandler);
-    mgr->SetDebugLevel(10);
-
-    AliAnalysisTaskKineFilter *kinefilter = new AliAnalysisTaskKineFilter("Kine Filter");
-    mgr->AddTask(kinefilter);
-    
-    
-    AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
-    //jetana->SetConfigFile("ConfigJetAnalysisFastjetMC.C");
-    jetana->SetConfigFile("${ALICE_ROOT}/HLT/JET/macros/tasks/ConfigJetAnalysisHLTMC.C");
-    jetana->SetDebugLevel(10);
-    mgr->AddTask(jetana);
-
-    //
-    // Create containers for input/output
-    AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(), 
-                                                            AliAnalysisManager::kInputContainer);
-
-    AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
-                                                             AliAnalysisManager::kOutputContainer, "default");
-
-    AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
-                                                             AliAnalysisManager::kOutputContainer, "histos.root");
-
-    mgr->ConnectInput  (kinefilter,     0, cinput1  );
-    mgr->ConnectOutput (kinefilter,     0, coutput1 );
-
-    mgr->ConnectInput  (jetana,     0, cinput1  );
-    mgr->ConnectOutput (jetana,     0, coutput1 );
-    mgr->ConnectOutput (jetana,     1, coutput2 );
-
-    //
-    // Run the analysis
-    //    
-
-    mgr->InitAnalysis();
-    mgr->PrintStatus();
-    mgr->StartAnalysis("local",chain);
-}
diff --git a/HLT/JET/macros/tasks/JetAnalysisManagerHLTMC.C b/HLT/JET/macros/tasks/JetAnalysisManagerHLTMC.C
new file mode 100644 (file)
index 0000000..a418adb
--- /dev/null
@@ -0,0 +1,148 @@
+void JetAnalysisManagerHLTMC() {
+
+  if ( getenv("FASTJET") ) {
+    gSystem->Load("libCGAL.so");
+    gSystem->Load("libfastjet.so");
+    gSystem->Load("libSISConePlugin.so");
+  }
+
+  gSystem->Load("libTree.so");
+  gSystem->Load("libPhysics.so");
+  gSystem->Load("libGeom.so");
+  gSystem->Load("libVMC.so");
+  gSystem->Load("libANALYSIS.so");
+  gSystem->Load("libSTEERBase.so");
+  gSystem->Load("libAOD.so");
+  gSystem->Load("libESD.so");
+  gSystem->Load("libANALYSISalice.so");
+  gSystem->Load("libJETAN.so");
+  gSystem->Load("libHLTbase.so");
+  gSystem->Load("libAliHLTUtil.so");
+  gSystem->Load("libAliHLTJET.so");
+
+  Int_t debugLevel = 10;
+  
+  // --------------------------------------------------------------------------------  
+    
+  if (gApplication) gApplication->InitializeGraphics();
+
+  // --------------------------------------------------------------------------------  
+  //
+  // Create the chain
+  //
+  // --------------------------------------------------------------------------------  
+
+  gROOT->LoadMacro("${ALICE_ROOT}/PWG0/CreateESDChain.C");
+  TChain* chain = new TChain("TE");
+
+  chain->Add("/home/jthaeder/jet/data/HEAD_2009-10-26/FastGen/kPythia6Jets125_150_14TeV/JET-ETA=-0.2,0.2_JET-ET=50,1000_R=0.4_200ev/galice.root");
+
+  // --------------------------------------------------------------------------------  
+  //
+  // Create the analysis manager
+  //
+  // --------------------------------------------------------------------------------  
+
+  // Output
+  AliAODHandler* aodHandler = new AliAODHandler();
+  aodHandler->SetOutputFileName("aod.root");
+
+  // MC Truth
+  AliMCEventHandler* mcHandler = new AliMCEventHandler();
+  mcHandler->SetReadTR(kFALSE); // for fastgen
+
+  AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
+  mgr->SetOutputEventHandler (aodHandler);
+  mgr->SetMCtruthEventHandler(mcHandler);
+  mgr->SetDebugLevel(debugLevel);
+
+  // --------------------------------------------------------------------------------
+  //
+  // Set Configfiles
+  //
+  // --------------------------------------------------------------------------------
+
+  // -- HLT FFSC - MC
+#if FFSC
+  AliAnalysisTaskJets *taskFFSC = new AliAnalysisTaskJets("JetAnalysisHLTMC_FFSC");
+  taskFFSC->SetConfigFile("./tasks/ConfigJetAnalysisHLTMC.C");
+  taskFFSC->SetNonStdBranch("jetsHLTMC_FFSC");
+  taskFFSC->SetDebugLevel(debugLevel);
+  mgr->AddTask(taskFFSC);
+#endif
+
+  // -- HLT FastJet - Kt
+#if FASTJET_KT
+  if ( getenv("FASTJET") ) {
+    AliAnalysisTaskJets *taskKt = new AliAnalysisTaskJets("JetAnalysisHLTMC_Kt");
+    taskKt->SetConfigFile("./tasks/ConfigJetAnalysisHLTMCKt.C");
+    taskKt->SetNonStdBranch("jetsHLTMC_Kt");
+    taskKt->SetDebugLevel(debugLevel);
+    mgr->AddTask(taskKt);
+  }
+#endif
+
+#if FASTJET_ANTIKT
+  // -- HLT FastJet - AntiKt
+  if ( getenv("FASTJET") ) {
+    AliAnalysisTaskJets *taskAntiKt = new AliAnalysisTaskJets("JetAnalysisHLTMC_AntiKt");
+    taskAntiKt->SetConfigFile("./tasks/ConfigJetAnalysisHLTMCAntiKt.C");
+    taskAntiKt->SetNonStdBranch("jetsHLTMC_AntiKt");
+    taskAntiKt->SetDebugLevel(debugLevel);
+    mgr->AddTask(taskAntiKt);
+  }
+#endif
+
+  // --------------------------------------------------------------------------------
+  //
+  // Create containers for input/output
+  //
+  // --------------------------------------------------------------------------------  
+
+  AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
+  AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
+
+#if FFSC
+  AliAnalysisDataContainer *coutputFFSC = mgr->CreateContainer("histos", TList::Class(), 
+                                                              AliAnalysisManager::kOutputContainer, "histos_MC_FFSC.root");
+  mgr->ConnectInput  (taskFFSC,     0, cinput1  );
+  mgr->ConnectOutput (taskFFSC,     0, coutput1 );
+  mgr->ConnectOutput (taskFFSC,     1, coutputFFSC );
+#endif
+  
+#if FASTJET_KT
+  if ( getenv("FASTJET") ) {
+    AliAnalysisDataContainer *coutputKt = mgr->CreateContainer("histos", TList::Class(), 
+                                                              AliAnalysisManager::kOutputContainer, "histos_MC_KT.root");
+    mgr->ConnectInput  (taskKt,     0, cinput1  );
+    mgr->ConnectOutput (taskKt,     0, coutput1 );
+    mgr->ConnectOutput (taskKt,     1, coutputKt );
+  }
+#endif
+
+#if FASTJET_ANTIKT
+  if ( getenv("FASTJET") ) {
+    AliAnalysisDataContainer *coutputAntiKt = mgr->CreateContainer("histos", TList::Class(), 
+                                                                  AliAnalysisManager::kOutputContainer, "histos_MC_ANTIKT.root");
+    mgr->ConnectInput  (taskAntiKt,     0, cinput1  );
+    mgr->ConnectOutput (taskAntiKt,     0, coutput1 );
+    mgr->ConnectOutput (taskAntiKt,     1, coutputAntiKt );
+  }
+#endif
+
+  //  AliAnalysisTaskKineFilter *kinefilter = new AliAnalysisTaskKineFilter("Kine Filter");
+  //  mgr->AddTask(kinefilter);
+  //  mgr->ConnectInput  (kinefilter,     0, cinput1  );
+  //  mgr->ConnectOutput (kinefilter,     0, coutput1 );
+
+  // --------------------------------------------------------------------------------  
+  //
+  // Run the analysis
+  //    
+  // --------------------------------------------------------------------------------  
+
+  mgr->InitAnalysis();
+  mgr->PrintStatus();
+  mgr->StartAnalysis("local",chain);
+}