-AliAnalysisTaskJetResponse* AddTaskJetResponse(Char_t* type = "clusters", Char_t* jf = "FASTKT", Float_t radius = 0.4, UInt_t filterMask = 256 , Float_t ptTrackMin = 0.15, Int_t iBack = 1, Int_t eventClassMin = 1, Int_t eventClassMax = 4){
+AliAnalysisTaskJetResponse* AddTaskJetResponse(Char_t* type = "clusters", Char_t* jf = "FASTKT", Float_t radius = 0.4, UInt_t filterMask = 256 , Float_t ptTrackMin = 0.15, Int_t iBack = 1, Int_t eventClassMin = 1, Int_t eventClassMax = 4, Int_t reactionPlaneBin = -1){
Printf("adding task jet response\n");
//task->SetJetDeltaEta(0.2);
//task->SetJetDeltaPhi(0.2);
+ task->SetReactionPlaneBin(reactionPlaneBin);
+
+ task->SetJetPtMin(0.); // min jet pt is implicit a cut on delta pT!!
+ //task->SetJetEtaMin(-.4); // reduced tracking acceptance in
+ //task->SetJetEtaMax(.4); // LHC10h pass1
+
+ //task->SetNMatchJets(1); // leading jets only
+ task->SetEventClassMode(1); //0: centrality (default), 1: multiplicity
+
mgr->AddTask(task);
+
+ if(reactionPlaneBin>=0){
+ suffix2 += Form("_rp%d", reactionPlaneBin);
+ }
+
AliAnalysisDataContainer *coutputJetResponse = mgr->CreateContainer(
Form("jetresponse_%s%s", type,suffix2.Data()), TList::Class(), AliAnalysisManager::kOutputContainer,
Form("%s:PWG4_JetResponse_%s%s", AliAnalysisManager::GetCommonFileName(), type, suffix2.Data()));