116ebfadef785c4425b9436606a471fde85fdb2d
[u/mrichter/AliRoot.git] / PWG2 / FORWARD / analysis2 / AddTaskCentraldNdeta.C
1 /**
2  * @file   AddTaskCentraldNdeta.C
3  * @author Christian Holm Christensen <cholm@nbi.dk>
4  * @date   Fri Jan 28 10:22:26 2011
5  * 
6  * @brief Script to add a multiplicity task for the central
7  *        @f$\eta@f$ region
8  * 
9  * 
10  */
11 AliAnalysisTask*
12 AddTaskCentraldNdeta(const char* trig="INEL", 
13                      Double_t vzMin=-10, 
14                      Double_t vzMax=10, 
15                      Float_t centLow=0, 
16                      Float_t centHigh=100)
17 {
18   // analysis manager
19   AliAnalysisManager* mgr = AliAnalysisManager::GetAnalysisManager();
20   
21   // Make our object.  2nd argumenent is absolute max Eta 
22   // 3rd argument is absolute max Vz
23   AliCentraldNdetaTask* task = new AliCentraldNdetaTask("Central");
24   task->SetVertexRange(vzMin, vzMax);
25   task->SetTriggerMask(trig);
26   // task->SetUseShapeCorrection(false);
27   task->AddCentralityBin( 0,  0); // All bin - integrate over centrality
28   task->AddCentralityBin( 0,  5);
29   task->AddCentralityBin( 5, 10);
30   task->AddCentralityBin(10, 20);
31   task->AddCentralityBin(20, 30);
32   task->AddCentralityBin(30, 40);
33   task->AddCentralityBin(40, 50);
34   task->AddCentralityBin(50, 60);
35   task->AddCentralityBin(60,100);
36   mgr->AddTask(task);
37
38   // create containers for input/output
39   AliAnalysisDataContainer *sums = 
40     mgr->CreateContainer("CentralSums", TList::Class(), 
41                          AliAnalysisManager::kOutputContainer, 
42                          AliAnalysisManager::GetCommonFileName());
43   AliAnalysisDataContainer *output = 
44     mgr->CreateContainer("CentralResults", TList::Class(), 
45                          AliAnalysisManager::kParamContainer, 
46                          AliAnalysisManager::GetCommonFileName());
47   
48   // connect input/output
49   mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer());
50   mgr->ConnectOutput(task, 1, sums);
51   mgr->ConnectOutput(task, 2, output);
52
53   return task;
54 }
55
56   
57 //________________________________________________________________________
58 //
59 // EOF
60 //