1 AliEPSelectionTask *AddTaskEventPlaneTPC(Bool_t useEtaGap=kFALSE,Float_t etaGap=0.,Bool_t posTPC=kFALSE,TString containername = "EPStat")
3 // Macro to connect an event plane selection task to an existing analysis manager.
5 if(useEtaGap && posTPC){
6 ::Error("AddTaskEventplane", "eta-splitting of events and one side of TPC not possible at same time!");
10 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
12 ::Error("AddTaskEventplane", "No analysis manager to connect to.");
15 // Check the analysis type using the event handlers connected to the analysis manager.
16 //==============================================================================
17 if (!mgr->GetInputEventHandler()) {
18 ::Error("AddTaskEventplane", "This task requires an input event handler");
21 TString inputDataType = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
24 AliEPSelectionTask *eventplaneTask = new AliEPSelectionTask("EventplaneSelection");
25 eventplaneTask->SelectCollisionCandidates(AliVEvent::kMB | AliVEvent::kSemiCentral | AliVEvent::kCentral);
26 if (inputDataType == "AOD"){
27 eventplaneTask->SetInput("AOD");
29 eventplaneTask->SetTrackType("TPC");
30 eventplaneTask->SetUsePtWeight();
31 eventplaneTask->SetUsePhiWeight();
32 eventplaneTask->SetSaveTrackContribution();
34 eventplaneTask->SetSubeventsSplitMethod(AliEPSelectionTask::kEta);
35 eventplaneTask->SetEtaGap(etaGap);
38 if(inputDataType == "AOD")
40 eventplaneTask->SetPersonalAODtrackCuts(128,0.,0.8,0.15,20.);
41 eventplaneTask->SetSubeventsSplitMethod(AliEPSelectionTask::kRandom);
44 AliESDtrackCuts *trackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
45 trackCuts->SetPtRange(0.15,20.);
46 trackCuts->SetEtaRange(0.0,0.8);
47 eventplaneTask->SetPersonalESDtrackCuts(trackCuts);
48 eventplaneTask->SetSubeventsSplitMethod(AliEPSelectionTask::kRandom);
53 mgr->AddTask(eventplaneTask);
55 AliAnalysisDataContainer *cinput0 = mgr->GetCommonInputContainer();
56 //AliAnalysisDataContainer *coutput1 = mgr->CreateContainer(containername,
57 // TList::Class(), AliAnalysisManager::kOutputContainer,
58 // "EventStat_temp.root");
60 mgr->ConnectInput(eventplaneTask, 0, mgr->GetCommonInputContainer());
61 //mgr->ConnectOutput(eventplaneTask,1,coutput1);
63 return eventplaneTask;