further updates Q cumulants

diff --git a/PWG2/FLOW/AliAnalysisTaskQCumulants.cxx b/PWG2/FLOW/AliAnalysisTaskQCumulants.cxx
index bdc4c45ed4142792db02bd026b7b5a1aef8be347..9db89e70fd1d5cb99b736cdbcf22294d0dc2e75a 100644 (file)
--- a/PWG2/FLOW/AliAnalysisTaskQCumulants.cxx
+++ b/PWG2/FLOW/AliAnalysisTaskQCumulants.cxx
@@ -55,6 +55,7 @@
 #include "AliFlowAnalysisWithQCumulants.h"
 #include "AliFlowCumuConstants.h"
 #include "AliFlowCommonConstants.h"
+#include "AliFlowCommonHist.h"
 #include "AliFlowCommonHistResults.h"
 #include "AliQCumulantsFunctions.h"
 
@@ -323,6 +324,18 @@ void AliAnalysisTaskQCumulants::Terminate(Option_t *)
   //final results for covariances (1st bin <2*4>-<2>*<4>, 2nd bin <2*6>-<2>*<6>, ...)
   TH1D *covariances = dynamic_cast<TH1D*>(fListHistos->FindObject("fCovariances"));
   
+  //common control histograms (taking into account only the events with 2 and more particles)  
+  AliFlowCommonHist *commonHist2nd = dynamic_cast<AliFlowCommonHist*>(fListHistos->FindObject("AliFlowCommonHist2ndOrderQC"));
+  
+  //common control histograms (taking into account only the events with 4 and more particles)  
+  AliFlowCommonHist *commonHist4th = dynamic_cast<AliFlowCommonHist*>(fListHistos->FindObject("AliFlowCommonHist4thOrderQC"));
+  
+  //common control histograms (taking into account only the events with 6 and more particles)  
+  AliFlowCommonHist *commonHist6th = dynamic_cast<AliFlowCommonHist*>(fListHistos->FindObject("AliFlowCommonHist6thOrderQC"));
+  
+  //common control histograms (taking into account only the events with 8 and more particles)  
+  AliFlowCommonHist *commonHist8th = dynamic_cast<AliFlowCommonHist*>(fListHistos->FindObject("AliFlowCommonHist8thOrderQC"));
+  
   //common histograms to store the final results for the 2nd order integrated and differential flow
   AliFlowCommonHistResults *commonHistRes2nd = dynamic_cast<AliFlowCommonHistResults*>(fListHistos->FindObject("AliFlowCommonHistResults2ndOrderQC"));
   
@@ -387,6 +400,11 @@ void AliAnalysisTaskQCumulants::Terminate(Option_t *)
   fQCA->SetDiffFlowResults2nd(diffFlowResults2ndOrder);
   fQCA->SetDiffFlowResults4th(diffFlowResults4thOrder); 
   fQCA->SetCovariances(covariances); 
+  
+  fQCA->SetCommonHists2nd(commonHist2nd); 
+  fQCA->SetCommonHists4th(commonHist4th);
+  fQCA->SetCommonHists6th(commonHist6th);
+  fQCA->SetCommonHists8th(commonHist8th);
 
   fQCA->SetCommonHistsResults2nd(commonHistRes2nd); 
   fQCA->SetCommonHistsResults4th(commonHistRes4th);

diff --git a/PWG2/FLOW/AliFlowAnalysisWithQCumulants.cxx b/PWG2/FLOW/AliFlowAnalysisWithQCumulants.cxx
index 31b6ac97c1f62c2ebfe2434ebbbc6797a2e05b85..b4da5be4f323e68f6dbc2c5998472029627f3526 100644 (file)
--- a/PWG2/FLOW/AliFlowAnalysisWithQCumulants.cxx
+++ b/PWG2/FLOW/AliFlowAnalysisWithQCumulants.cxx
@@ -1500,8 +1500,8 @@ void AliFlowAnalysisWithQCumulants::Make(AliFlowEventSimple* anEvent)
 void AliFlowAnalysisWithQCumulants::Finish()
 {
  //calculate the final results
- AliQCumulantsFunctions finalResults(fIntFlowResultsQC,fDiffFlowResults2ndOrderQC,fDiffFlowResults4thOrderQC,fCovariances,fAvMultIntFlowQC,fQvectorComponents,fQCorrelations, fQProduct,fDirectCorrelations,f2PerPtBin1n1nRP,f2PerPtBin2n2nRP,f3PerPtBin2n1n1nRP,f3PerPtBin1n1n2nRP,f4PerPtBin1n1n1n1nRP, f2PerEtaBin1n1nRP,f2PerEtaBin2n2nRP,f3PerEtaBin2n1n1nRP,f3PerEtaBin1n1n2nRP,f4PerEtaBin1n1n1n1nRP, f2PerPtBin1n1nPOI,f2PerPtBin2n2nPOI,f3PerPtBin2n1n1nPOI,f3PerPtBin1n1n2nPOI,f4PerPtBin1n1n1n1nPOI, f2PerEtaBin1n1nPOI,f2PerEtaBin2n2nPOI,f3PerEtaBin2n1n1nPOI,f3PerEtaBin1n1n2nPOI,f4PerEtaBin1n1n1n1nPOI,fCommonHistsResults2nd, fCommonHistsResults4th,fCommonHistsResults6th,fCommonHistsResults8th);
-         
+ AliQCumulantsFunctions finalResults(fIntFlowResultsQC,fDiffFlowResults2ndOrderQC,fDiffFlowResults4thOrderQC,fCovariances,fAvMultIntFlowQC,fQvectorComponents,fQCorrelations, fQProduct,fDirectCorrelations,f2PerPtBin1n1nRP,f2PerPtBin2n2nRP,f3PerPtBin2n1n1nRP,f3PerPtBin1n1n2nRP,f4PerPtBin1n1n1n1nRP, f2PerEtaBin1n1nRP,f2PerEtaBin2n2nRP,f3PerEtaBin2n1n1nRP,f3PerEtaBin1n1n2nRP,f4PerEtaBin1n1n1n1nRP, f2PerPtBin1n1nPOI,f2PerPtBin2n2nPOI,f3PerPtBin2n1n1nPOI,f3PerPtBin1n1n2nPOI,f4PerPtBin1n1n1n1nPOI, f2PerEtaBin1n1nPOI,f2PerEtaBin2n2nPOI,f3PerEtaBin2n1n1nPOI,f3PerEtaBin1n1n2nPOI,f4PerEtaBin1n1n1n1nPOI,fCommonHists2nd,fCommonHists4th, fCommonHists6th, fCommonHists8th,fCommonHistsResults2nd, fCommonHistsResults4th,fCommonHistsResults6th,fCommonHistsResults8th);
+
  finalResults.Calculate();  
 }
 

diff --git a/PWG2/FLOW/AliFlowAnalysisWithQCumulants.h b/PWG2/FLOW/AliFlowAnalysisWithQCumulants.h
index 6c8816433869c8b9c4a2ce14657f5ca10b2d433a..270e37a65dbf620a877be10944aa7b056d166ce6 100644 (file)
--- a/PWG2/FLOW/AliFlowAnalysisWithQCumulants.h
+++ b/PWG2/FLOW/AliFlowAnalysisWithQCumulants.h
@@ -59,6 +59,18 @@ class AliFlowAnalysisWithQCumulants{
   void SetCovariances(TH1D* const cov) {this->fCovariances = cov;};
   TH1D* GetCovariances() const {return this->fCovariances;};
   
+  void SetCommonHists2nd(AliFlowCommonHist* const ch2nd) {this->fCommonHists2nd = ch2nd;};
+  AliFlowCommonHist* GetCommonHists2nd() const {return this->fCommonHists2nd;};
+  
+  void SetCommonHists4th(AliFlowCommonHist* const ch4th) {this->fCommonHists4th = ch4th;};
+  AliFlowCommonHist* GetCommonHists4th() const {return this->fCommonHists4th;};
+  
+  void SetCommonHists6th(AliFlowCommonHist* const ch6th) {this->fCommonHists6th = ch6th;};
+  AliFlowCommonHist* GetCommonHists6th() const {return this->fCommonHists6th;};
+  
+  void SetCommonHists8th(AliFlowCommonHist* const ch8th) {this->fCommonHists8th = ch8th;};
+  AliFlowCommonHist* GetCommonHists8th() const {return this->fCommonHists8th;};
+  
   void SetCommonHistsResults2nd(AliFlowCommonHistResults* const chr2nd) {this->fCommonHistsResults2nd = chr2nd;};
   AliFlowCommonHistResults* GetCommonHistsResults2nd() const {return this->fCommonHistsResults2nd;};
   
@@ -215,10 +227,10 @@ class AliFlowAnalysisWithQCumulants{
   TProfile*                  f3PerEtaBin1n1n2nPOI;      //<<3'>>_{n,n|2n} per eta-bin
   TProfile*                  f4PerEtaBin1n1n1n1nPOI;    //<<4'>>_{n,n|n,n} per eta-bin  
  
-  AliFlowCommonHist*         fCommonHists2nd;           //common control histograms for 2nd order
-  AliFlowCommonHist*         fCommonHists4th;           //common control histograms for 4th order
-  AliFlowCommonHist*         fCommonHists6th;           //common control histograms for 6th order
-  AliFlowCommonHist*         fCommonHists8th;           //common control histograms for 8th order
+  AliFlowCommonHist*         fCommonHists2nd;           //common control histograms (taking into account only the events with 2 and more particles) 
+  AliFlowCommonHist*         fCommonHists4th;           //common control histograms (taking into account only the events with 4 and more particles) 
+  AliFlowCommonHist*         fCommonHists6th;           //common control histograms (taking into account only the events with 6 and more particles) 
+  AliFlowCommonHist*         fCommonHists8th;           //common control histograms (taking into account only the events with 8 and more particles) 
   
   AliFlowCommonHistResults*  fCommonHistsResults2nd;    //final results for 2nd order int. and diff. flow stored in the common histograms 
   AliFlowCommonHistResults*  fCommonHistsResults4th;    //final results for 4th order int. and diff. flow stored in the common histograms 

diff --git a/PWG2/FLOW/AliQCumulantsFunctions.cxx b/PWG2/FLOW/AliQCumulantsFunctions.cxx
index a43542ca92b8d50e7a014c430aba2b153f030966..eef410901d09efdf6a978388cd9faae7ae863688 100644 (file)
--- a/PWG2/FLOW/AliQCumulantsFunctions.cxx
+++ b/PWG2/FLOW/AliQCumulantsFunctions.cxx
@@ -42,6 +42,7 @@
 #include "AliFlowTrackSimple.h"
 #include "AliFlowAnalysisWithCumulants.h"
 #include "AliFlowCommonConstants.h"
+#include "AliFlowCommonHist.h"
 #include "AliFlowCommonHistResults.h"
 #include "AliQCumulantsFunctions.h"
 
@@ -79,7 +80,11 @@ AliQCumulantsFunctions::AliQCumulantsFunctions():
  fbinEta3p2n1n1nPOI(NULL),
  fbinEta3p1n1n2nPOI(NULL),
  fbinEta4p1n1n1n1nPOI(NULL),  
- fchr2nd(NULL),
+ fch2nd(NULL),//ch = common histogram (control)
+ fch4th(NULL),
+ fch6th(NULL),
+ fch8th(NULL), 
+ fchr2nd(NULL),//chr = common histogram results 
  fchr4th(NULL),
  fchr6th(NULL),
  fchr8th(NULL) 
@@ -92,7 +97,7 @@ AliQCumulantsFunctions::~AliQCumulantsFunctions()
  //destructor
 }
 
-AliQCumulantsFunctions::AliQCumulantsFunctions(TH1D *intRes, TH1D *diffRes2nd, TH1D *diffRes4th, TH1D *covar, TProfile *AvMult, TProfile *QVector, TProfile *QCorr, TProfile *QProd, TProfile *Direct, TProfile *binPt2p1n1nRP, TProfile *binPt2p2n2nRP, TProfile *binPt3p2n1n1nRP, TProfile *binPt3p1n1n2nRP, TProfile *binPt4p1n1n1n1nRP, TProfile *binEta2p1n1nRP, TProfile *binEta2p2n2nRP, TProfile *binEta3p2n1n1nRP, TProfile *binEta3p1n1n2nRP, TProfile *binEta4p1n1n1n1nRP, TProfile *binPt2p1n1nPOI, TProfile *binPt2p2n2nPOI, TProfile *binPt3p2n1n1nPOI, TProfile *binPt3p1n1n2nPOI, TProfile *binPt4p1n1n1n1nPOI, TProfile *binEta2p1n1nPOI, TProfile *binEta2p2n2nPOI, TProfile *binEta3p2n1n1nPOI, TProfile *binEta3p1n1n2nPOI, TProfile *binEta4p1n1n1n1nPOI, AliFlowCommonHistResults *chr2nd, AliFlowCommonHistResults *chr4th, AliFlowCommonHistResults *chr6th, AliFlowCommonHistResults *chr8th):
+AliQCumulantsFunctions::AliQCumulantsFunctions(TH1D *intRes, TH1D *diffRes2nd, TH1D *diffRes4th, TH1D *covar, TProfile *AvMult, TProfile *QVector, TProfile *QCorr, TProfile *QProd, TProfile *Direct, TProfile *binPt2p1n1nRP, TProfile *binPt2p2n2nRP, TProfile *binPt3p2n1n1nRP, TProfile *binPt3p1n1n2nRP, TProfile *binPt4p1n1n1n1nRP, TProfile *binEta2p1n1nRP, TProfile *binEta2p2n2nRP, TProfile *binEta3p2n1n1nRP, TProfile *binEta3p1n1n2nRP, TProfile *binEta4p1n1n1n1nRP, TProfile *binPt2p1n1nPOI, TProfile *binPt2p2n2nPOI, TProfile *binPt3p2n1n1nPOI, TProfile *binPt3p1n1n2nPOI, TProfile *binPt4p1n1n1n1nPOI, TProfile *binEta2p1n1nPOI, TProfile *binEta2p2n2nPOI, TProfile *binEta3p2n1n1nPOI, TProfile *binEta3p1n1n2nPOI, TProfile *binEta4p1n1n1n1nPOI, AliFlowCommonHist *ch2nd, AliFlowCommonHist *ch4th, AliFlowCommonHist *ch6th, AliFlowCommonHist *ch8th, AliFlowCommonHistResults *chr2nd, AliFlowCommonHistResults *chr4th, AliFlowCommonHistResults *chr6th, AliFlowCommonHistResults *chr8th):
  fIntRes(intRes),
  fDiffRes2nd(diffRes2nd),
  fDiffRes4th(diffRes4th),
@@ -121,7 +126,11 @@ AliQCumulantsFunctions::AliQCumulantsFunctions(TH1D *intRes, TH1D *diffRes2nd, T
  fbinEta2p2n2nPOI(binEta2p2n2nPOI),
  fbinEta3p2n1n1nPOI(binEta3p2n1n1nPOI),
  fbinEta3p1n1n2nPOI(binEta3p1n1n2nPOI),
- fbinEta4p1n1n1n1nPOI(binEta4p1n1n1n1nPOI),   
+ fbinEta4p1n1n1n1nPOI(binEta4p1n1n1n1nPOI), 
+ fch2nd(ch2nd),
+ fch4th(ch4th),
+ fch6th(ch6th),
+ fch8th(ch8th),   
  fchr2nd(chr2nd),
  fchr4th(chr4th),
  fchr6th(chr6th),
@@ -313,24 +322,38 @@ for(Int_t bb=1;bb<nBinsEtaRP+1;bb++)
  }
 }      
 
-
 //differential flow (POI)
+Int_t nBinsPtPOI = fbinPt2p1n1nPOI->GetNbinsX();
+Int_t nBinsEtaPOI = fbinEta2p1n1nPOI->GetNbinsX();
 
 //Pt:
 Double_t secondOrderQCumulantDiffFlowPtPOI = 0.;
 Double_t fourthOrderQCumulantDiffFlowPtPOI = 0.;
 
-Int_t nBinsPtPOI = fbinPt2p1n1nPOI->GetNbinsX();
+Double_t dVn2ndPOI=0.,dDiffvn2nd=0.,dYield2nd=0.,dSum2nd=0.;
+
+Double_t dVn2ndPOIEta=0.,dDiffvn2ndEta=0.,dYield2ndEta=0.,dSum2ndEta=0.;//to be removed
 
 for(Int_t bb=1;bb<nBinsPtPOI+1;bb++)
 {
+ //QC{2}
  if(fbinPt2p1n1nPOI->GetBinEntries(bb)>0.&&vn2!=0)
  {
+  //cout<<"bin = "<<bb<<" : "<<(fch2nd->GetHistPtDiff())->GetBinContent(bb)<<endl;
+  //cout<<endl;
   secondOrderQCumulantDiffFlowPtPOI = fbinPt2p1n1nPOI->GetBinContent(bb);
   fDiffRes2nd->SetBinContent(bb,secondOrderQCumulantDiffFlowPtPOI/vn2);
   //common histogram:
-  fchr2nd->FillDifferentialFlowPtPOI(bb,secondOrderQCumulantDiffFlowPtPOI/vn2, 0.);//to be improved (errors)
+  fchr2nd->FillDifferentialFlowPtPOI(bb,secondOrderQCumulantDiffFlowPtPOI/vn2, 0.);//to be improved (errors) && bb or bb+1
+  //-------------------------------------------------------------
+  //integrated flow (POI, Pt, 2nd order):
+  dDiffvn2nd=(fchr2nd->GetHistDiffFlowPtPOI())->GetBinContent(bb);
+  dYield2nd=(fch2nd->GetHistPtDiff())->GetBinContent(bb);
+  dVn2ndPOI+=dDiffvn2nd*dYield2nd;
+  dSum2nd+=dYield2nd;
+  //-------------------------------------------------------------
  }
+ //QC{4]
  if(fbinPt4p1n1n1n1nPOI->GetBinEntries(bb)>0.&&vn4!=0.)
  {
   fourthOrderQCumulantDiffFlowPtPOI = fbinPt4p1n1n1n1nPOI->GetBinContent(bb)-2.*fbinPt2p1n1nPOI->GetBinContent(bb)*pow(vn2,2.);
@@ -340,12 +363,21 @@ for(Int_t bb=1;bb<nBinsPtPOI+1;bb++)
  }
 }      
 
+//storing the final results for integrated flow (POI):
+//QC{2}
+if(dSum2nd && fchr2nd)
+{
+ dVn2ndPOI/=dSum2nd;
+ fchr2nd->FillIntegratedFlowPOI(dVn2ndPOI,0.);
+}
+
+
+
+
 //Eta:
 Double_t secondOrderQCumulantDiffFlowEtaPOI = 0.;
 Double_t fourthOrderQCumulantDiffFlowEtaPOI = 0.;
 
-Int_t nBinsEtaPOI = fbinEta2p1n1nPOI->GetNbinsX();
-
 for(Int_t bb=1;bb<nBinsEtaPOI+1;bb++)
 {
  if(fbinEta2p1n1nPOI->GetBinEntries(bb)>0.&&vn2!=0)
@@ -369,7 +401,35 @@ for(Int_t bb=1;bb<nBinsEtaPOI+1;bb++)
  
 
  
- 
+//xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+//         !!!! to be removed !!!!
+for(Int_t bb=1;bb<nBinsEtaPOI+1;bb++)
+{
+ //integrated flow (POI, Eta, 2nd order):
+  dDiffvn2ndEta=(fchr2nd->GetHistDiffFlowEtaPOI())->GetBinContent(bb);
+  dYield2ndEta=(fch2nd->GetHistEtaDiff())->GetBinContent(bb);
+ // cout<<"bin = "<<bb<<" : "<<dDiffvn2ndEta<<" "<<dYield2ndEta<<endl;
+ // cout<<endl;
+  dVn2ndPOIEta+=dDiffvn2ndEta*dYield2ndEta;
+  dSum2ndEta+=dYield2ndEta;
+
+
+//check why for small number of events there is a mismatch here (should be exactly the same)://to be removed
+//cout<<dSum2nd<<"   to be improved   "<<(fch2nd->GetHistPtDiff())->Integral()<<endl;//to be removed
+}
+if(dSum2ndEta)
+{
+ dVn2ndPOIEta/=dSum2ndEta;
+ //cout<<"POI Pt vs POI Eta:"<<endl;
+ //cout<<dVn2ndPOI<<" vs "<<dVn2ndPOIEta<<endl;
+ //cout<<dSum2nd<<" vs "<<dSum2ndEta<<endl;
+// cout<<"eta = "<<(fchr2nd->GetHistDiffFlowEtaPOI())->GetBinContent(44)<<endl;
+// cout<<"poi = "<<(fchr2nd->GetHistDiffFlowPtPOI())->GetBinContent(14)<<endl;
+// cout<<"+++++++++"<<endl; 
+}
+
+//xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
+
  
  /*
  Double_t first=16.*pow(two*twoErr,2.);

diff --git a/PWG2/FLOW/AliQCumulantsFunctions.h b/PWG2/FLOW/AliQCumulantsFunctions.h
index 5e2c2651fbad76939903f2d456cd9e7fba30c34d..77e7b3ad256248391f13e5f4bda736a75c4e36ed 100644 (file)
--- a/PWG2/FLOW/AliQCumulantsFunctions.h
+++ b/PWG2/FLOW/AliQCumulantsFunctions.h
@@ -35,7 +35,7 @@ class AliQCumulantsFunctions{
  public:
   AliQCumulantsFunctions();
   virtual ~AliQCumulantsFunctions();
-  AliQCumulantsFunctions(TH1D *intRes, TH1D *diffRes2nd, TH1D *diffRes4th, TH1D *covar, TProfile *AvMult, TProfile *QVector, TProfile *QCorr, TProfile *QProd, TProfile *Direct, TProfile *binPt2p1n1nRP, TProfile *binPt2p2n2nRP, TProfile *binPt3p2n1n1nRP, TProfile *binPt3p1n1n2nRP, TProfile *binPt4p1n1n1n1nRP, TProfile *binEta2p1n1nRP, TProfile *binEta2p2n2nRP, TProfile *binEta3p2n1n1nRP, TProfile *binEta3p1n1n2nRP, TProfile *binEta4p1n1n1n1nRP, TProfile *binPt2p1n1nPOI, TProfile *binPt2p2n2nPOI, TProfile *binPt3p2n1n1nPOI, TProfile *binPt3p1n1n2nPOI, TProfile *binPt4p1n1n1n1nPOI, TProfile *binEta2p1n1nPOI, TProfile *binEta2p2n2nPOI, TProfile *binEta3p2n1n1nPOI, TProfile *binEta3p1n1n2nPOI, TProfile *binEta4p1n1n1n1nPOI, AliFlowCommonHistResults *chr2nd, AliFlowCommonHistResults *chr4th, AliFlowCommonHistResults *chr6th, AliFlowCommonHistResults *chr8th);
+  AliQCumulantsFunctions(TH1D *intRes, TH1D *diffRes2nd, TH1D *diffRes4th, TH1D *covar, TProfile *AvMult, TProfile *QVector, TProfile *QCorr, TProfile *QProd, TProfile *Direct, TProfile *binPt2p1n1nRP, TProfile *binPt2p2n2nRP, TProfile *binPt3p2n1n1nRP, TProfile *binPt3p1n1n2nRP, TProfile *binPt4p1n1n1n1nRP, TProfile *binEta2p1n1nRP, TProfile *binEta2p2n2nRP, TProfile *binEta3p2n1n1nRP, TProfile *binEta3p1n1n2nRP, TProfile *binEta4p1n1n1n1nRP, TProfile *binPt2p1n1nPOI, TProfile *binPt2p2n2nPOI, TProfile *binPt3p2n1n1nPOI, TProfile *binPt3p1n1n2nPOI, TProfile *binPt4p1n1n1n1nPOI, TProfile *binEta2p1n1nPOI, TProfile *binEta2p2n2nPOI, TProfile *binEta3p2n1n1nPOI, TProfile *binEta3p1n1n2nPOI, TProfile *binEta4p1n1n1n1nPOI, AliFlowCommonHist *ch2nd, AliFlowCommonHist *ch4th, AliFlowCommonHist *ch6th, AliFlowCommonHist *ch8th, AliFlowCommonHistResults *chr2nd, AliFlowCommonHistResults *chr4th, AliFlowCommonHistResults *chr6th, AliFlowCommonHistResults *chr8th);
  
   void Calculate();
 
@@ -79,6 +79,11 @@ class AliQCumulantsFunctions{
   TProfile *fbinEta3p1n1n2nPOI;   //<<3'>>_{n,n|2n} per eta-bin
   TProfile *fbinEta4p1n1n1n1nPOI; //<<4'>>_{n,n|n,n} per eta-bin  
   
+  AliFlowCommonHist *fch2nd; //common control histograms (taking into account only the events with 2 and more particles) 
+  AliFlowCommonHist *fch4th; //common control histograms (taking into account only the events with 4 and more particles) 
+  AliFlowCommonHist *fch6th; //common control histograms (taking into account only the events with 6 and more particles) 
+  AliFlowCommonHist *fch8th; //common control histograms (taking into account only the events with 8 and more particles) 
+  
   AliFlowCommonHistResults *fchr2nd; //final results for 2nd order int. and diff. flow stored in the common histograms
   AliFlowCommonHistResults *fchr4th; //final results for 4th order int. and diff. flow stored in the common histograms
   AliFlowCommonHistResults *fchr6th; //final results for 6th order int. and diff. flow stored in the common histograms