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("libSTEERBase.so");
   8     gSystem->Load("libAOD.so");
   9     gSystem->Load("libESD.so");
  10     gSystem->Load("libJETAN.so");
  11      //
  12     if (gApplication) gApplication->InitializeGraphics();
  13     // Create the chain
  14     //
  15     TChain* chain = new TChain("esdTree");
  16     chain->Add("/home/morsch/AliRoot/data/data_jets102/AliESDs.root");
  17     chain->Add("/home/morsch/AliRoot/data/data_jets103/AliESDs.root");
  18     //
  19     // Create the analysis manager
  20     //
  21     AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
  22     mgr->SetDebugLevel(10);
  23     //
  24     // Common output service
  25     AliAODHandler* aodHandler   = new AliAODHandler();
  26     aodHandler->SetOutputFileName("aod.root");
  27     mgr->SetOutputEventHandler(aodHandler);
  28     //
  29     // Common MC truth services
  30     AliMCEventHandler* mcHandler = new AliMCEventHandler();
  31     mgr->SetMCtruthEventHandler(mcHandler);
  32
  33     //
  34     // Jet analysis
  35     AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
  36     jetana->SetDebugLevel(10);
  37     mgr->AddTask(jetana);
  38
  39     //
  40     // Create containers for input/output
  41     AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(),
  42                                                              AliAnalysisManager::kInputContainer);
  43
  44     AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
  45                                                               AliAnalysisManager::kOutputContainer, "aod.root");
  46
  47     AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TH1F::Class(),
  48                                                               AliAnalysisManager::kOutputContainer, "histos.root");
  49
  50     mgr->ConnectInput  (jetana,  0, cinput1  );
  51     mgr->ConnectOutput (jetana,  0, coutput1 );
  52     mgr->ConnectOutput (jetana,  1, coutput2 );
  53     //
  54     // Run the analysis
  55     //
  56     mgr->InitAnalysis();
  57     mgr->PrintStatus();
  58     mgr->StartAnalysis("local",chain);
  59 }