]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PWGJE/macros/AddTaskJetShape.C
lhc12g
[u/mrichter/AliRoot.git] / PWGJE / macros / AddTaskJetShape.C
1 AliAnalysisTaskJetShape* AddTaskJetShape(const char* bRec1, const char* bRec2, const char* bBkg1, const char *bBkg2,
2                                          UInt_t filterMask = 272 ,Bool_t kIsPbPb = kFALSE, UInt_t kTriggerMask=0, Int_t eventClassMin = 0, Int_t eventClassMax = 4)
3 {
4
5    Printf("adding task jet shape\n");
6
7    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
8    if(!mgr){
9       ::Error("AddTaskJetShape", "No analysis manager to connect to.");
10       return NULL;
11    }
12    if(!mgr->GetInputEventHandler()){
13       ::Error("AddTaskJetShape", "This task requires an input event handler.");
14       return NULL;
15    }
16
17      
18   
19
20   TString sRec(bRec1);
21   TString sGen(bRec2);
22   TString sRecBkg(bBkg1);
23   TString sGenBkg(bBkg2);
24       
25   AliAnalysisTaskJetShape *task = new AliAnalysisTaskJetShape(Form("JetShape_%s_%s_%d",bRec1,bRec2,kTriggerMask));
26    
27    task->SetBranchNames(sRec,sGen);
28    task->SetBackgroundBranch(sRecBkg, sGenBkg);
29    task->SetOfflineTrgMask(kTriggerMask);
30    task->SetEvtClassMin(eventClassMin);
31    task->SetEvtClassMax(eventClassMax);
32    task->SetCentMin(0.);
33    task->SetCentMax(100.);
34    task->SetFilterMask(filterMask); 
35    task->SetJetPtCorrMin(20.,20);   
36    task->SetPbPb(kIsPbPb);
37    mgr->AddTask(task);
38
39
40    AliAnalysisDataContainer *coutputJetShape = mgr->CreateContainer(Form("pwgjeJetShape_%s_%s_%s_%s_%d",bRec1,bRec2,bBkg1,bBkg2,kTriggerMask), TList::Class(),AliAnalysisManager::kOutputContainer,Form("%s:PWGJE_JetShape_%s_%s_%d",AliAnalysisManager::GetCommonFileName(),bRec1,bRec2,kTriggerMask));
41
42
43    mgr->ConnectInput (task, 0, mgr->GetCommonInputContainer());
44    mgr->ConnectOutput(task, 0, mgr->GetCommonOutputContainer());
45    mgr->ConnectOutput(task, 1, coutputJetShape);
46
47    return task;
48 }