vertex from tracks
[u/mrichter/AliRoot.git] / ANALYSIS / DphiAnalysis.C
CommitLineData
bd0e8153 1void 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
15e62203 46 AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
bd0e8153 47
15e62203 48// AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
49// AliAnalysisManager::kOutputContainer, "default");
bd0e8153 50 AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
51 AliAnalysisManager::kOutputContainer, "histos.root");
52
53
54 mgr->ConnectInput (dphiana, 0, cinput1 );
55 mgr->ConnectOutput (dphiana, 1, coutput2 );
56
57 //
58 // Run the analysis
59 //
60 mgr->InitAnalysis();
61 mgr->PrintStatus();
62 mgr->StartAnalysis("mix",chain, 1000);
63}