include Event Plane (only In-Plane / Out-Of-Plane so far) in Event Mixing

diff --git a/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskBFPsi.cxx b/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskBFPsi.cxx
index 0592e1af1caed5486a612dd47e31e410faaeb561..106f1eccdfb671dcd6ef1ee6c06614dad1000673 100755 (executable)
--- a/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskBFPsi.cxx
+++ b/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskBFPsi.cxx
@@ -345,16 +345,22 @@ void AliAnalysisTaskBFPsi::UserCreateOutputObjects() {
   Int_t trackDepth = fMixingTracks; \r
   Int_t poolsize   = 1000;  // Maximum number of events, ignored in the present implemented of AliEventPoolManager\r
    \r
+  // centrality bins\r
   Double_t centralityBins[] = {0.,1.,2.,3.,4.,5.,6.,7.,8.,9.,10.,15.,20.,25.,30.,35.,40.,45.,50.,55.,60.,65.,70.,75.,80.,90.,100.}; // SHOULD BE DEDUCED FROM CREATED ALITHN!!!\r
   Double_t* centbins        = centralityBins;\r
   Int_t nCentralityBins     = sizeof(centralityBins) / sizeof(Double_t) - 1;\r
   \r
-  // bins for second buffer are shifted by 100 cm\r
+  // Zvtx bins\r
   Double_t vertexBins[] = {-10., -7., -5., -3., -1., 1., 3., 5., 7., 10.}; // SHOULD BE DEDUCED FROM CREATED ALITHN!!!\r
   Double_t* vtxbins     = vertexBins;\r
   Int_t nVertexBins     = sizeof(vertexBins) / sizeof(Double_t) - 1;\r
 \r
-  fPoolMgr = new AliEventPoolManager(poolsize, trackDepth, nCentralityBins, centbins, nVertexBins, vtxbins);\r
+  // Event plane angle (Psi) bins\r
+  Double_t psiBins[] = {0.,45.,135.,215.,305.,360.}; // SHOULD BE DEDUCED FROM CREATED ALITHN!!!\r
+  Double_t* psibins     = psiBins;\r
+  Int_t nPsiBins     = sizeof(psiBins) / sizeof(Double_t) - 1;\r
+\r
+  fPoolMgr = new AliEventPoolManager(poolsize, trackDepth, nCentralityBins, centbins, nVertexBins, vtxbins, nPsiBins, psibins);\r
 \r
   if(fESDtrackCuts) fList->Add(fESDtrackCuts);\r
 \r
@@ -493,10 +499,10 @@ void AliAnalysisTaskBFPsi::UserExec(Option_t *) {
       //    FillCorrelations(). Also nMix should be passed in, so a weight\r
       //    of 1./nMix can be applied.\r
       \r
-      AliEventPool* pool = fPoolMgr->GetEventPool(fCentrality, eventMain->GetPrimaryVertex()->GetZ());\r
+      AliEventPool* pool = fPoolMgr->GetEventPool(fCentrality, eventMain->GetPrimaryVertex()->GetZ(),gReactionPlane);\r
       \r
       if (!pool){\r
-       AliFatal(Form("No pool found for centrality = %f, zVtx = %f", fCentrality, eventMain->GetPrimaryVertex()->GetZ()));\r
+       AliFatal(Form("No pool found for centrality = %f, zVtx = %f, psi = %f", fCentrality, eventMain->GetPrimaryVertex()->GetZ(),gReactionPlane));\r
       }\r
       else{\r
        \r