Changes to the DEtaDPhi corr fun (making 0 in the center of the bin in the histogram...
[u/mrichter/AliRoot.git] / ANALYSIS / AliPriorsTask.C
CommitLineData
9ab172e7 1
2Bool_t AliPriorsTask()
3{
4
5 TBenchmark benchmark;
6 benchmark.Start("AliPriorsTask");
7
8 AliLog::SetGlobalDebugLevel(0);
9
10 Load() ; //load the required libraries
11
12 TChain * analysisChain = new TChain("esdTree");
13 analysisChain->Add("AliESDs.root");
14
15 Int_t nIter = 5;
16
17 // create the task
18 AliPriorsTask *task = new AliPriorsTask("Task for Priors");
19 Double_t priors[5]={0.2,0.2,0.2,0.2,0.2};
20 task->SetPriors(priors);
21 task->SetNiterations(nIter);
22
23 // Make the analysis manager
24 AliAnalysisManager *mgr = new AliAnalysisManager("TestIteration");
25
26 AliMCEventHandler* mcHandler = new AliMCEventHandler();
27 mgr->SetMCtruthEventHandler(mcHandler);
28
29
30 AliESDInputHandler* esdHandler = new AliESDInputHandler();
31 mgr->SetInputEventHandler(esdHandler);
32
33
34 // Create and connect containers for input/output
15e62203 35 AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
9ab172e7 36 //slot 0 : default output tree (by default handled by AliAnalysisTaskSE)
37 AliAnalysisDataContainer *coutput0 = mgr->CreateContainer("ctree0", TTree::Class(),AliAnalysisManager::kOutputContainer,"output.root");
38 // output TH1I for event counting
39 AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("chist0", TH1I::Class(),AliAnalysisManager::kOutputContainer,"output.root");
40 AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("cgraph0", TH2D::Class(),AliAnalysisManager::kOutputContainer,"output.root");
41
42
43 cinput0->SetData(analysisChain);
44
45 mgr->AddTask(task);
46 mgr->ConnectInput(task,0,cinput0);
47 mgr->ConnectOutput(task,0,coutput0);
48 mgr->ConnectOutput(task,1,coutput1);
49 mgr->ConnectOutput(task,2,coutput2);
50
51
52
53 AliEventPoolLoop* loop = new AliEventPoolLoop(nIter);
54 loop->SetChain(analysisChain);
55 mgr->SetEventPool(loop);
56
57 //mgr->SetDebugLevel(1);
58 printf("READY TO RUN\n");
59 //RUN !!!
60 if (mgr->InitAnalysis()) {
61 mgr->PrintStatus();
62 mgr->StartAnalysis("mixing",analysisChain);
63 //mgr->StartAnalysis("local",analysisChain);
64
65 }
66
67 benchmark.Stop("AliPriorsTask");
68 benchmark.Show("AliPriorsTask");
69
70 return kTRUE ;
71}
72
73void Load() {
74
75 //load the required aliroot libraries
76 gSystem->Load("libANALYSIS") ;
77 gSystem->Load("libANALYSISalice") ;
78
79 //compile online the task class
80 gSystem->SetIncludePath("-I. -I$ALICE_ROOT/include -I$ROOTSYS/include");
81 gROOT->LoadMacro("./AliPriorsTask.cxx++g");
82}