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