mixing example
[u/mrichter/AliRoot.git] / 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 }