Added option to set event selection
authorfnoferin <fnoferin@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 26 Feb 2013 11:00:50 +0000 (11:00 +0000)
committerfnoferin <fnoferin@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 26 Feb 2013 11:00:50 +0000 (11:00 +0000)
PWGCF/FLOW/macros/AddTaskFlowCentralityPIDQC.C
PWGCF/FLOW/macros/AddTaskFlowCentralityPIDSP.C
PWGCF/FLOW/macros/AddTaskPidQC.C
PWGCF/FLOW/macros/AddTaskPidVZEROSP.C

index 75a2ef2..13f41fe 100644 (file)
@@ -7,7 +7,8 @@
 //\r
 /////////////////////////////////////////////////////////////////////////////////////////////\r
 \r
-void AddTaskFlowCentralityPIDQC( Float_t centrMin=10.,\r
+void AddTaskFlowCentralityPIDQC(Int_t centralitysel,\r
+                               Float_t centrMin=10.,\r
                                Float_t centrMax=20.,\r
                                TString fileNameBase="output",\r
                                Bool_t isPID = kTRUE,\r
@@ -297,6 +298,7 @@ void AddTaskFlowCentralityPIDQC( Float_t centrMin=10.,
     AliFMDAnalysisTaskSE *taskfmd = NULL;\r
     if (rptypestr == "FMD") {\r
       taskfmd = new AliFMDAnalysisTaskSE("TaskFMD");\r
+      taskfmd->SelectCollisionCandidates(centralitysel);\r
       mgr->AddTask(taskfmd);\r
       \r
       AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();\r
@@ -319,8 +321,10 @@ void AddTaskFlowCentralityPIDQC( Float_t centrMin=10.,
       taskFE->SetFlow(v1,v2,v3,v4); \r
       taskFE->SetNonFlowNumberOfTrackClones(numberOfTrackClones);\r
       taskFE->SetAfterburnerOn();\r
+      taskFE->SelectCollisionCandidates(centralitysel);\r
     }\r
-  else {taskFE = new AliAnalysisTaskFlowEvent(Form("TaskFlowEvent %s",outputSlotName.Data()),"",doQA); }\r
+  else {taskFE = new AliAnalysisTaskFlowEvent(Form("TaskFlowEvent %s",outputSlotName.Data()),"",doQA); \r
+    taskFE->SelectCollisionCandidates(centralitysel);}\r
   if (ExcludeRegion) {\r
     taskFE->DefineDeadZone(excludeEtaMin, excludeEtaMax, excludePhiMin, excludePhiMax); \r
   }\r
@@ -347,40 +351,47 @@ void AddTaskFlowCentralityPIDQC( Float_t centrMin=10.,
   //===========================================================================\r
   if (SP){\r
     AliAnalysisTaskScalarProduct *taskSP = new AliAnalysisTaskScalarProduct(Form("TaskScalarProduct %s",outputSlotName.Data()),WEIGHTS[0]);\r
+    taskSP->SelectCollisionCandidates(centralitysel);\r
     taskSP->SetRelDiffMsub(1.0);\r
     taskSP->SetApplyCorrectionForNUA(kTRUE);\r
     mgr->AddTask(taskSP);\r
   }\r
   if (LYZ1SUM){\r
     AliAnalysisTaskLeeYangZeros *taskLYZ1SUM = new AliAnalysisTaskLeeYangZeros(Form("TaskLeeYangZerosSUM %s",outputSlotName.Data()),kTRUE);\r
+    taskLYZ1SUM->SelectCollisionCandidates(centralitysel);\r
     taskLYZ1SUM->SetFirstRunLYZ(kTRUE);\r
     taskLYZ1SUM->SetUseSumLYZ(kTRUE);\r
     mgr->AddTask(taskLYZ1SUM);\r
   }\r
   if (LYZ1PROD){\r
     AliAnalysisTaskLeeYangZeros *taskLYZ1PROD = new AliAnalysisTaskLeeYangZeros(Form("TaskLeeYangZerosPROD %s",outputSlotName.Data()),kTRUE);\r
+    taskLYZ1PROD->SelectCollisionCandidates(centralitysel);\r
     taskLYZ1PROD->SetFirstRunLYZ(kTRUE);\r
     taskLYZ1PROD->SetUseSumLYZ(kFALSE);\r
     mgr->AddTask(taskLYZ1PROD);\r
   }\r
   if (LYZ2SUM){\r
     AliAnalysisTaskLeeYangZeros *taskLYZ2SUM = new AliAnalysisTaskLeeYangZeros(Form("TaskLeeYangZerosSUM %s",outputSlotName.Data()),kFALSE);\r
+    taskLYZ2SUM->SelectCollisionCandidates(centralitysel);\r
     taskLYZ2SUM->SetFirstRunLYZ(kFALSE);\r
     taskLYZ2SUM->SetUseSumLYZ(kTRUE);\r
     mgr->AddTask(taskLYZ2SUM);\r
   }\r
   if (LYZ2PROD){\r
     AliAnalysisTaskLeeYangZeros *taskLYZ2PROD = new AliAnalysisTaskLeeYangZeros(Form("TaskLeeYangZerosPROD %s",outputSlotName.Data()),kFALSE);\r
+    taskLYZ2PROD->SelectCollisionCandidates(centralitysel);\r
     taskLYZ2PROD->SetFirstRunLYZ(kFALSE);\r
     taskLYZ2PROD->SetUseSumLYZ(kFALSE);\r
     mgr->AddTask(taskLYZ2PROD);\r
   }\r
   if (LYZEP){\r
     AliAnalysisTaskLYZEventPlane *taskLYZEP = new AliAnalysisTaskLYZEventPlane(Form("TaskLYZEventPlane %s",outputSlotName.Data()));\r
+    taskLYZEP->SelectCollisionCandidates(centralitysel);\r
     mgr->AddTask(taskLYZEP);\r
   }\r
   if (GFC){\r
     AliAnalysisTaskCumulants *taskGFC = new AliAnalysisTaskCumulants(Form("TaskCumulants %s",outputSlotName.Data()),useWeights);\r
+    taskGFC->SelectCollisionCandidates(centralitysel);\r
     taskGFC->SetUsePhiWeights(WEIGHTS[0]); \r
     taskGFC->SetUsePtWeights(WEIGHTS[1]);\r
     taskGFC->SetUseEtaWeights(WEIGHTS[2]); \r
@@ -388,6 +399,7 @@ void AddTaskFlowCentralityPIDQC( Float_t centrMin=10.,
   }\r
   if (QC){\r
     AliAnalysisTaskQCumulants *taskQC = new AliAnalysisTaskQCumulants(Form("TaskQCumulants %s",outputSlotName.Data()),useWeights);\r
+    taskQC->SelectCollisionCandidates(centralitysel);\r
     taskQC->SetUsePhiWeights(WEIGHTS[0]); \r
     taskQC->SetUsePtWeights(WEIGHTS[1]);\r
     taskQC->SetUseEtaWeights(WEIGHTS[2]); \r
@@ -402,6 +414,7 @@ void AddTaskFlowCentralityPIDQC( Float_t centrMin=10.,
   }\r
   if (FQD){\r
     AliAnalysisTaskFittingQDistribution *taskFQD = new AliAnalysisTaskFittingQDistribution(Form("TaskFittingQDistribution %s",outputSlotName.Data()),kFALSE);\r
+    taskFQD->SelectCollisionCandidates(centralitysel);\r
     taskFQD->SetUsePhiWeights(WEIGHTS[0]); \r
     taskFQD->SetqMin(0.);\r
     taskFQD->SetqMax(1000.);\r
@@ -410,10 +423,12 @@ void AddTaskFlowCentralityPIDQC( Float_t centrMin=10.,
   }\r
   if (MCEP){\r
     AliAnalysisTaskMCEventPlane *taskMCEP = new AliAnalysisTaskMCEventPlane(Form("TaskMCEventPlane %s",outputSlotName.Data()));\r
+    taskMCEP->SelectCollisionCandidates(centralitysel);\r
     mgr->AddTask(taskMCEP);\r
   }\r
   if (MH){\r
     AliAnalysisTaskMixedHarmonics *taskMH = new AliAnalysisTaskMixedHarmonics(Form("TaskMixedHarmonics %s",outputSlotName.Data()),useWeights);\r
+    taskMH->SelectCollisionCandidates(centralitysel);\r
     taskMH->SetHarmonic(1); // n in cos[n(phi1+phi2-2phi3)] and cos[n(psi1+psi2-2phi3)]\r
     taskMH->SetNoOfMultipicityBins(10000);\r
     taskMH->SetMultipicityBinWidth(1.);\r
@@ -425,6 +440,7 @@ void AddTaskFlowCentralityPIDQC( Float_t centrMin=10.,
   }  \r
   if (NL){\r
     AliAnalysisTaskNestedLoops *taskNL = new AliAnalysisTaskNestedLoops(Form("TaskNestedLoops %s",outputSlotName.Data()),useWeights);\r
+    taskNL->SelectCollisionCandidates(centralitysel);\r
     taskNL->SetHarmonic(1); // n in cos[n(phi1+phi2-2phi3)] and cos[n(psi1+psi2-2phi3)]\r
     taskNL->SetEvaluateNestedLoopsForRAD(kTRUE); // RAD = Relative Angle Distribution\r
     taskNL->SetEvaluateNestedLoopsForMH(kTRUE); // evalaute <<cos[n(phi1+phi2-2phi3)]>> (Remark: three nested loops)   \r
@@ -628,6 +644,7 @@ void AddTaskFlowCentralityPIDQC( Float_t centrMin=10.,
   if (runQAtask)\r
   {\r
     AliAnalysisTaskQAflow* taskQAflow = new AliAnalysisTaskQAflow(Form("TaskQAflow %s",outputSlotName.Data()));\r
+    taskQAflow->SelectCollisionCandidates(centralitysel);\r
     taskQAflow->SetEventCuts(cutsEvent);\r
     taskQAflow->SetTrackCuts(cutsRP);\r
     taskQAflow->SetFillNTuple(FillQAntuple);\r
index 6ea9e63..da5e14b 100644 (file)
@@ -7,7 +7,8 @@
 //\r
 /////////////////////////////////////////////////////////////////////////////////////////////\r
 \r
-void AddTaskFlowCentralityPIDSP( Float_t centrMin=10.,\r
+void AddTaskFlowCentralityPIDSP(Int_t centralitysel,\r
+                               Float_t centrMin=10.,\r
                                Float_t centrMax=20.,\r
                                TString fileNameBase="output",\r
                                Bool_t isPID = kTRUE,\r
@@ -294,6 +295,7 @@ void AddTaskFlowCentralityPIDSP( Float_t centrMin=10.,
     AliFMDAnalysisTaskSE *taskfmd = NULL;\r
     if (rptypestr == "FMD") {\r
       taskfmd = new AliFMDAnalysisTaskSE("TaskFMD");\r
+      taskfmd->SelectCollisionCandidates(centralitysel);\r
       mgr->AddTask(taskfmd);\r
       \r
       AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();\r
@@ -313,11 +315,13 @@ void AddTaskFlowCentralityPIDSP( Float_t centrMin=10.,
   if(useAfterBurner)\r
     { \r
       taskFE = new AliAnalysisTaskFlowEvent(Form("TaskFlowEvent %s",outputSlotName.Data()),"",doQA,1);\r
+      taskFE->SelectCollisionCandidates(centralitysel);\r
       taskFE->SetFlow(v1,v2,v3,v4); \r
       taskFE->SetNonFlowNumberOfTrackClones(numberOfTrackClones);\r
       taskFE->SetAfterburnerOn();\r
     }\r
-  else {taskFE = new AliAnalysisTaskFlowEvent(Form("TaskFlowEvent %s",outputSlotName.Data()),"",doQA); }\r
+  else {taskFE = new AliAnalysisTaskFlowEvent(Form("TaskFlowEvent %s",outputSlotName.Data()),"",doQA); \r
+      taskFE->SelectCollisionCandidates(centralitysel);}\r
   if (ExcludeRegion) {\r
     taskFE->DefineDeadZone(excludeEtaMin, excludeEtaMax, excludePhiMin, excludePhiMax); \r
   }\r
@@ -344,40 +348,47 @@ void AddTaskFlowCentralityPIDSP( Float_t centrMin=10.,
   //===========================================================================\r
   if (SP){\r
     AliAnalysisTaskScalarProduct *taskSP = new AliAnalysisTaskScalarProduct(Form("TaskScalarProduct %s",outputSlotName.Data()),WEIGHTS[0]);\r
+    taskSP->SelectCollisionCandidates(centralitysel);\r
     taskSP->SetRelDiffMsub(1.0);\r
     taskSP->SetApplyCorrectionForNUA(kTRUE);\r
     mgr->AddTask(taskSP);\r
   }\r
   if (LYZ1SUM){\r
     AliAnalysisTaskLeeYangZeros *taskLYZ1SUM = new AliAnalysisTaskLeeYangZeros(Form("TaskLeeYangZerosSUM %s",outputSlotName.Data()),kTRUE);\r
+    taskLYZ1SUM->SelectCollisionCandidates(centralitysel);\r
     taskLYZ1SUM->SetFirstRunLYZ(kTRUE);\r
     taskLYZ1SUM->SetUseSumLYZ(kTRUE);\r
     mgr->AddTask(taskLYZ1SUM);\r
   }\r
   if (LYZ1PROD){\r
     AliAnalysisTaskLeeYangZeros *taskLYZ1PROD = new AliAnalysisTaskLeeYangZeros(Form("TaskLeeYangZerosPROD %s",outputSlotName.Data()),kTRUE);\r
+    taskLYZ1PROD->SelectCollisionCandidates(centralitysel);\r
     taskLYZ1PROD->SetFirstRunLYZ(kTRUE);\r
     taskLYZ1PROD->SetUseSumLYZ(kFALSE);\r
     mgr->AddTask(taskLYZ1PROD);\r
   }\r
   if (LYZ2SUM){\r
     AliAnalysisTaskLeeYangZeros *taskLYZ2SUM = new AliAnalysisTaskLeeYangZeros(Form("TaskLeeYangZerosSUM %s",outputSlotName.Data()),kFALSE);\r
+    taskLYZ2SUM->SelectCollisionCandidates(centralitysel);\r
     taskLYZ2SUM->SetFirstRunLYZ(kFALSE);\r
     taskLYZ2SUM->SetUseSumLYZ(kTRUE);\r
     mgr->AddTask(taskLYZ2SUM);\r
   }\r
   if (LYZ2PROD){\r
     AliAnalysisTaskLeeYangZeros *taskLYZ2PROD = new AliAnalysisTaskLeeYangZeros(Form("TaskLeeYangZerosPROD %s",outputSlotName.Data()),kFALSE);\r
+    taskLYZ2PROD->SelectCollisionCandidates(centralitysel);\r
     taskLYZ2PROD->SetFirstRunLYZ(kFALSE);\r
     taskLYZ2PROD->SetUseSumLYZ(kFALSE);\r
     mgr->AddTask(taskLYZ2PROD);\r
   }\r
   if (LYZEP){\r
     AliAnalysisTaskLYZEventPlane *taskLYZEP = new AliAnalysisTaskLYZEventPlane(Form("TaskLYZEventPlane %s",outputSlotName.Data()));\r
+    taskLYZEP->SelectCollisionCandidates(centralitysel);\r
     mgr->AddTask(taskLYZEP);\r
   }\r
   if (GFC){\r
     AliAnalysisTaskCumulants *taskGFC = new AliAnalysisTaskCumulants(Form("TaskCumulants %s",outputSlotName.Data()),useWeights);\r
+    taskGFC->SelectCollisionCandidates(centralitysel);\r
     taskGFC->SetUsePhiWeights(WEIGHTS[0]); \r
     taskGFC->SetUsePtWeights(WEIGHTS[1]);\r
     taskGFC->SetUseEtaWeights(WEIGHTS[2]); \r
@@ -385,6 +396,7 @@ void AddTaskFlowCentralityPIDSP( Float_t centrMin=10.,
   }\r
   if (QC){\r
     AliAnalysisTaskQCumulants *taskQC = new AliAnalysisTaskQCumulants(Form("TaskQCumulants %s",outputSlotName.Data()),useWeights);\r
+    taskQC->SelectCollisionCandidates(centralitysel);\r
     taskQC->SetUsePhiWeights(WEIGHTS[0]); \r
     taskQC->SetUsePtWeights(WEIGHTS[1]);\r
     taskQC->SetUseEtaWeights(WEIGHTS[2]); \r
@@ -399,6 +411,7 @@ void AddTaskFlowCentralityPIDSP( Float_t centrMin=10.,
   }\r
   if (FQD){\r
     AliAnalysisTaskFittingQDistribution *taskFQD = new AliAnalysisTaskFittingQDistribution(Form("TaskFittingQDistribution %s",outputSlotName.Data()),kFALSE);\r
+    taskFQD->SelectCollisionCandidates(centralitysel);\r
     taskFQD->SetUsePhiWeights(WEIGHTS[0]); \r
     taskFQD->SetqMin(0.);\r
     taskFQD->SetqMax(1000.);\r
@@ -407,10 +420,12 @@ void AddTaskFlowCentralityPIDSP( Float_t centrMin=10.,
   }\r
   if (MCEP){\r
     AliAnalysisTaskMCEventPlane *taskMCEP = new AliAnalysisTaskMCEventPlane(Form("TaskMCEventPlane %s",outputSlotName.Data()));\r
+    taskMCEP->SelectCollisionCandidates(centralitysel);\r
     mgr->AddTask(taskMCEP);\r
   }\r
   if (MH){\r
     AliAnalysisTaskMixedHarmonics *taskMH = new AliAnalysisTaskMixedHarmonics(Form("TaskMixedHarmonics %s",outputSlotName.Data()),useWeights);\r
+    taskMH->SelectCollisionCandidates(centralitysel);\r
     taskMH->SetHarmonic(1); // n in cos[n(phi1+phi2-2phi3)] and cos[n(psi1+psi2-2phi3)]\r
     taskMH->SetNoOfMultipicityBins(10000);\r
     taskMH->SetMultipicityBinWidth(1.);\r
@@ -422,6 +437,7 @@ void AddTaskFlowCentralityPIDSP( Float_t centrMin=10.,
   }  \r
   if (NL){\r
     AliAnalysisTaskNestedLoops *taskNL = new AliAnalysisTaskNestedLoops(Form("TaskNestedLoops %s",outputSlotName.Data()),useWeights);\r
+    taskNL->SelectCollisionCandidates(centralitysel);\r
     taskNL->SetHarmonic(1); // n in cos[n(phi1+phi2-2phi3)] and cos[n(psi1+psi2-2phi3)]\r
     taskNL->SetEvaluateNestedLoopsForRAD(kTRUE); // RAD = Relative Angle Distribution\r
     taskNL->SetEvaluateNestedLoopsForMH(kTRUE); // evalaute <<cos[n(phi1+phi2-2phi3)]>> (Remark: three nested loops)   \r
@@ -625,6 +641,7 @@ void AddTaskFlowCentralityPIDSP( Float_t centrMin=10.,
   if (runQAtask)\r
   {\r
     AliAnalysisTaskQAflow* taskQAflow = new AliAnalysisTaskQAflow(Form("TaskQAflow %s",outputSlotName.Data()));\r
+    taskQAflow->SelectCollisionCandidates(centralitysel);\r
     taskQAflow->SetEventCuts(cutsEvent);\r
     taskQAflow->SetTrackCuts(cutsRP);\r
     taskQAflow->SetFillNTuple(FillQAntuple);\r
index 3208568..bc34ece 100644 (file)
@@ -1,4 +1,4 @@
-AddTaskPidQC(Float_t etamin=-0.8,Float_t etamax=0.8){\r
+AddTaskPidQC(Int_t centralityselection=AliVEvent::kAny,Float_t etamin=-0.8,Float_t etamax=0.8){\r
   gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskFlowCentralityPIDQC.C");\r
 \r
   const Int_t ncentr = 9;\r
@@ -6,15 +6,15 @@ AddTaskPidQC(Float_t etamin=-0.8,Float_t etamax=0.8){
   Int_t cmax[ncentr]={5,10,20,30,40,50,60,70,80};\r
 \r
   for(Int_t i=0;i < ncentr;i++){\r
-    AddTaskFlowCentralityPIDQC(cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,1,etamin,etamax); // no pid + QA\r
-    AddTaskFlowCentralityPIDQC(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDQC(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDQC(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,-1,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDQC(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDQC(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDQC(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDQC(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,-1,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDQC(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,1,etamin,etamax); // no pid + QA\r
+    AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,-1,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,-1,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
   }\r
 }\r
 \r
index 930843a..5699e98 100644 (file)
@@ -1,4 +1,4 @@
-AddTaskPidVZEROSP(Float_t etamin=-0.8,Float_t etamax=0.8){\r
+AddTaskPidVZEROSP(Int_t centralityselection=AliVEvent::kAny,Float_t etamin=-0.8,Float_t etamax=0.8){\r
   gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskFlowCentralityPIDSP.C");\r
 \r
   const Int_t ncentr = 9;\r
@@ -6,15 +6,15 @@ AddTaskPidVZEROSP(Float_t etamin=-0.8,Float_t etamax=0.8){
   Int_t cmax[ncentr]={5,10,20,30,40,50,60,70,80};\r
 \r
   for(Int_t i=0;i < ncentr;i++){\r
-    AddTaskFlowCentralityPIDSP(cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax); // no pid\r
-    AddTaskFlowCentralityPIDSP(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,-1,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,-1,2,0,etamin,etamax);\r
-    AddTaskFlowCentralityPIDSP(cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax); // no pid\r
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,-1,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,-1,2,0,etamin,etamax);\r
+    AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);\r
   }\r
 }\r
 \r