ANALYSIS/DphiAnalysis.C

   1 void DphiAnalysis()
   2 {
   3       gSystem->Load("libTree.so");
   4       gSystem->Load("libPhysics.so");
   5       gSystem->Load("libGeom.so");
   6       gSystem->Load("libVMC.so");
   7       gSystem->Load("libANALYSIS.so");
   8       gSystem->Load("libSTEERBase.so");
   9       gSystem->Load("libAOD.so");
  10       gSystem->Load("libESD.so");
  11       gSystem->Load("libANALYSISalice.so");
  12       gSystem->Load("libJETAN.so");
  13
  14      //
  15     if (gApplication) gApplication->InitializeGraphics();
  16     // Create the chain
  17     //
  18     TChain* chain = new TChain("aodTree");
  19     chain->Add("./AliAODs.root");
  20
  21     ///////////////////////////////////////////////////////////////////////////////////
  22     // Create the analysis manager
  23     //
  24     // Input
  25     AliMultiAODInputHandler* inpHandler = new AliMultiAODInputHandler(2);
  26     // Pool
  27     AliEventPoolOTF* pool = new AliEventPoolOTF("event pool", "event pool");
  28     pool->SetTagDirectory(".");
  29     pool->SetMultiplicityBin(0, 100, 1);
  30     pool->Init();
  31
  32     AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
  33     mgr->SetInputEventHandler  (inpHandler);
  34     mgr->SetEventPool(pool);
  35     inpHandler->SetEventPool(pool);
  36
  37
  38     mgr->SetDebugLevel(10);
  39     ///////////////////////////////////////////////////////////////////////////////////
  40     AliAnalysisTaskPhiCorr *dphiana = new AliAnalysisTaskPhiCorr("Phi Correlation Analysis");
  41     dphiana->SetDebugLevel(10);
  42     mgr->AddTask(dphiana);
  43
  44     //
  45     // Create containers for input/output
  46     AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(),
  47                                                              AliAnalysisManager::kInputContainer);
  48
  49     AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
  50                                                               AliAnalysisManager::kOutputContainer, "default");
  51     AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
  52                                                               AliAnalysisManager::kOutputContainer, "histos.root");
  53
  54
  55     mgr->ConnectInput  (dphiana,  0, cinput1  );
  56     mgr->ConnectOutput (dphiana,  1, coutput2 );
  57
  58     //
  59     // Run the analysis
  60     //
  61     mgr->InitAnalysis();
  62     mgr->PrintStatus();
  63     mgr->StartAnalysis("mix",chain, 1000);
  64 }