JETAN/JetAnalysisManagerLoc.C

   1 void JetAnalysisManagerLoc()
   2 {
   3     gSystem->Load("libTree.so");
   4     gSystem->Load("libGeom.so");
   5     gSystem->Load("libVMC.so");
   6     gSystem->Load("libANALYSIS.so");
   7     gSystem->Load("libAOD.so");
   8     gSystem->Load("libESD.so");
   9     gSystem->Load("libJETAN.so");
  10      //
  11     if (gApplication) gApplication->InitializeGraphics();
  12     // Create the chain
  13     //
  14     TChain* chain = new TChain("esdTree");
  15     chain->Add("/home/morsch/analysis/AliEn/Interactive/esd/001/AliESDs.root");
  16     chain->Add("/home/morsch/analysis/AliEn/Interactive/esd/002/AliESDs.root");
  17     //
  18     // Make the analysis manager
  19     //
  20     AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
  21     AliAODHandler* aodHandler   = new AliAODHandler();
  22     mgr->SetEventHandler(aodHandler);
  23     mgr-> SetDebugLevel(10);
  24
  25     AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
  26     jetana->SetDebugLevel(10);
  27     mgr->AddTask(jetana);
  28
  29     //
  30     // Create containers for input/output
  31     AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(),
  32                                                              AliAnalysisManager::kInputContainer);
  33
  34     AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
  35                                                               AliAnalysisManager::kOutputContainer, "aod.root");
  36
  37     mgr->ConnectInput  (jetana,  0, cinput1 );
  38     mgr->ConnectOutput (jetana,  0, coutput1 );
  39     //
  40     // Run the analysis
  41     //
  42     mgr->InitAnalysis();
  43     mgr->PrintStatus();
  44     mgr->StartAnalysis("local",chain);
  45 }