correct settings for POI and RP change name macro
authorsnelling <snelling@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 13 Feb 2009 19:01:49 +0000 (19:01 +0000)
committersnelling <snelling@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 13 Feb 2009 19:01:49 +0000 (19:01 +0000)
PWG2/FLOW/AliFlowCommon/AliFlowAnalysisWithMCEventPlane.cxx
PWG2/FLOW/AliFlowCommon/AliFlowAnalysisWithMCEventPlane.h
PWG2/FLOW/AliFlowCommon/AliQCumulantsFunctions.cxx
PWG2/FLOW/AliFlowTasks/AliAnalysisTaskMCEventPlane.cxx
PWG2/FLOW/macros/compareFlowResults.C [moved from PWG2/FLOW/macros/CompareFlowResults.C with 84% similarity]

index f068704..5bc6d3b 100644 (file)
@@ -52,7 +52,7 @@ ClassImp(AliFlowAnalysisWithMCEventPlane)
    fCommonHistsRes(NULL),
    fHistProFlow(NULL),
    fHistRP(NULL),
-   fHistProIntFlowRP(NULL),
+   fHistProIntFlow(NULL),
    fHistProDiffFlowPtRP(NULL),
    fHistProDiffFlowEtaRP(NULL),
    fHistProDiffFlowPtPOI(NULL),
@@ -125,11 +125,11 @@ void AliFlowAnalysisWithMCEventPlane::Init() {
   fHistRP->SetYTitle("Counts");
   fHistList->Add(fHistRP);
   
-  fHistProIntFlowRP = new TProfile("fHistProIntFlowRP","fHistProIntFlowRP",1,0.,1.);
-  fHistProIntFlowRP->SetLabelSize(0.06);
-  (fHistProIntFlowRP->GetXaxis())->SetBinLabel(1,"v_{n}{2}");
-  fHistProIntFlowRP->SetYTitle("");
-  fHistList->Add(fHistProIntFlowRP);
+  fHistProIntFlow = new TProfile("fHistProIntFlow","fHistProIntFlow",1,0.,1.);
+  fHistProIntFlow->SetLabelSize(0.06);
+  (fHistProIntFlow->GetXaxis())->SetBinLabel(1,"v_{n}{2}");
+  fHistProIntFlow->SetYTitle("");
+  fHistList->Add(fHistProIntFlow);
  
   fHistProDiffFlowPtRP = new TProfile("fHistProDiffFlowPtRP","fHistProDiffFlowPtRP",iNbinsPt,dPtMin,dPtMax);
   fHistProDiffFlowPtRP->SetXTitle("P_{t}");
@@ -194,8 +194,8 @@ void AliFlowAnalysisWithMCEventPlane::Make(AliFlowEventSimple* anEvent, Double_t
             dv2 = TMath::Cos(2*(dPhi-aRP));
            dPt = pTrack->Pt();
            dEta = pTrack->Eta();
-            //integrated flow (RP):
-            fHistProIntFlowRP->Fill(0.,dv2);
+            //no-name int. flow (to be improved):
+            fHistProIntFlow->Fill(0.,dv2);
             //differential flow (Pt, RP):
             fHistProDiffFlowPtRP->Fill(dPt,dv2,1.);
             //differential flow (Eta, RP):
@@ -208,8 +208,6 @@ void AliFlowAnalysisWithMCEventPlane::Make(AliFlowEventSimple* anEvent, Double_t
            dv2 = TMath::Cos(2*(dPhi-aRP));
            dPt = pTrack->Pt();
            dEta = pTrack->Eta();
-           //fill histogram
-           fHistProFlow->Fill(dPt,dv2);//to be removed 
            //differential flow (Pt, POI):
             fHistProDiffFlowPtPOI->Fill(dPt,dv2,1.);
             //differential flow (Eta, POI):
@@ -232,83 +230,106 @@ void AliFlowAnalysisWithMCEventPlane::Finish() {
   Int_t iNbinsPt  = AliFlowCommonConstants::GetNbinsPt();  
   Int_t iNbinsEta = AliFlowCommonConstants::GetNbinsEta(); 
          
-  //RP:
-  //integrated flow:
-  Double_t dVRP = fHistProIntFlowRP->GetBinContent(1);  
-  Double_t dErrVRP = fHistProIntFlowRP->GetBinError(1);//to be improved (treatment of errors for non-Gaussian distribution needed!)  
-  fCommonHistsRes->FillIntegratedFlowRP(dVRP,dErrVRP);
+  // no-name int. flow (to be improved):
+  Double_t dV = fHistProIntFlow->GetBinContent(1);  
+  Double_t dErrV = fHistProIntFlow->GetBinError(1); // to be improved (treatment of errors for non-Gaussian distribution needed!)  
+  // fill no-name int. flow (to be improved):
+  fCommonHistsRes->FillIntegratedFlow(dV,dErrV);
+  cout<<"dV{MC} is       "<<dV<<" +- "<<dErrV<<endl;
   
-  //differential flow (Pt): 
+  //RP:
+  TH1F* fHistPtRP = fCommonHists->GetHistPtInt(); // to be improved (change "int" and "diff" to RP and POI in common control histos)
+  Double_t dYieldPtRP = 0.;
+  Double_t dVRP = 0.;
+  Double_t dErrVRP = 0.;
+  Double_t dSumRP = 0.;
+  //differential flow (RP, Pt): 
   Double_t dvPtRP = 0.;           
   Double_t dErrvPtRP = 0.;
-  for(Int_t b=0;b<iNbinsPt;b++)
+  for(Int_t b=1;b<iNbinsPt;b++)
   {
-   dvPtRP    = fHistProDiffFlowPtRP->GetBinContent(b+1);
-   dErrvPtRP = fHistProDiffFlowPtRP->GetBinError(b+1);//to be improved (treatment of errors for non-Gaussian distribution needed!)
-   fCommonHistsRes->FillDifferentialFlowPtRP(b+1, dvPtRP, dErrvPtRP);
+   dvPtRP    = fHistProDiffFlowPtRP->GetBinContent(b);
+   dErrvPtRP = fHistProDiffFlowPtRP->GetBinError(b);//to be improved (treatment of errors for non-Gaussian distribution needed!)
+   fCommonHistsRes->FillDifferentialFlowPtRP(b, dvPtRP, dErrvPtRP);
+   if(fHistPtRP){
+       //integrated flow (RP)
+       dYieldPtRP = fHistPtRP->GetBinContent(b);
+       dVRP += dvPtRP*dYieldPtRP;
+       dSumRP += dYieldPtRP;
+       //error on integrated flow
+       dErrVRP += dYieldPtRP*dYieldPtRP*dErrvPtRP*dErrvPtRP;
+      }
   }
+  if (dSumRP != 0. ) {
+    dVRP /= dSumRP;  //because pt distribution should be normalised
+    dErrVRP /= (dSumRP*dSumRP);
+    dErrVRP = TMath::Sqrt(dErrVRP); 
+  }
+  // fill integrated flow (RP):
+  fCommonHistsRes->FillIntegratedFlowRP(dVRP,dErrVRP);
+  cout<<"dV{MC} (RP) is  "<<dVRP<<" +- "<<dErrVRP<<endl;
   
-  //differential flow (Eta): 
+  //differential flow (RP, Eta): 
   Double_t dvEtaRP = 0.;           
   Double_t dErrvEtaRP = 0.;
-  for(Int_t b=0;b<iNbinsEta;b++)
+  for(Int_t b=1;b<iNbinsEta;b++)
   {
-   dvEtaRP    = fHistProDiffFlowEtaRP->GetBinContent(b+1);
-   dErrvEtaRP = fHistProDiffFlowEtaRP->GetBinError(b+1);//to be improved (treatment of errors for non-Gaussian distribution needed!)
-   fCommonHistsRes->FillDifferentialFlowEtaRP(b+1, dvEtaRP, dErrvEtaRP);
+   dvEtaRP    = fHistProDiffFlowEtaRP->GetBinContent(b);
+   dErrvEtaRP = fHistProDiffFlowEtaRP->GetBinError(b);//to be improved (treatment of errors for non-Gaussian distribution needed!)
+   fCommonHistsRes->FillDifferentialFlowEtaRP(b, dvEtaRP, dErrvEtaRP);
   }
                                                                                                                                    
   //POI:
-  TH1F* fHistPtDiff = fCommonHists->GetHistPtDiff();
-  Double_t dYieldPt = 0.;
-  Double_t dV = 0.;
-  Double_t dErrV = 0.;
-  Double_t dSum = 0.;
-  Double_t dv2proPt = 0.;
-  Double_t dErrdifcombPt = 0.; 
-  Double_t dv2proEta = 0.;
-  Double_t dErrdifcombEta = 0.;   
+  TH1F* fHistPtPOI = fCommonHists->GetHistPtDiff(); // to be improved (change "int" and "diff" to RP and POI in common control histos)
+  Double_t dYieldPtPOI = 0.;
+  Double_t dVPOI = 0.;
+  Double_t dErrVPOI = 0.;
+  Double_t dSumPOI = 0.;
+  Double_t dv2proPtPOI = 0.;
+  Double_t dErrdifcombPtPOI = 0.; 
+  Double_t dv2proEtaPOI = 0.;
+  Double_t dErrdifcombEtaPOI = 0.;   
   //Pt:
   if(fHistProFlow && fHistProDiffFlowPtPOI) {//to be removed (fHistProFlow)
-    for(Int_t b=0;b<iNbinsPt;b++){
+    for(Int_t b=1;b<iNbinsPt;b++){
       //dv2pro = fHistProFlow->GetBinContent(b);//to be removed
       //dErrdifcomb = fHistProFlow->GetBinError(b);//to be removed
-      dv2proPt = fHistProDiffFlowPtPOI->GetBinContent(b+1);
-      dErrdifcombPt = fHistProDiffFlowPtPOI->GetBinError(b+1);//to be improved (treatment of errors for non-Gaussian distribution needed!)
+      dv2proPtPOI = fHistProDiffFlowPtPOI->GetBinContent(b);
+      dErrdifcombPtPOI = fHistProDiffFlowPtPOI->GetBinError(b);//to be improved (treatment of errors for non-Gaussian distribution needed!)
       //fill TH1D
-      fCommonHistsRes->FillDifferentialFlow(b+1, dv2proPt, dErrdifcombPt);//to be removed
-      fCommonHistsRes->FillDifferentialFlowPtPOI(b+1, dv2proPt, dErrdifcombPt); 
-      if (fHistPtDiff){
+      fCommonHistsRes->FillDifferentialFlowPtPOI(b, dv2proPtPOI, dErrdifcombPtPOI); 
+      if (fHistPtPOI){
        //integrated flow (POI)
-       dYieldPt = fHistPtDiff->GetBinContent(b+1);
-       dV += dv2proPt*dYieldPt;
-       dSum += dYieldPt;
+       dYieldPtPOI = fHistPtPOI->GetBinContent(b);
+       dVPOI += dv2proPtPOI*dYieldPtPOI;
+       dSumPOI += dYieldPtPOI;
        //error on integrated flow
-       dErrV += dYieldPt*dYieldPt*dErrdifcombPt*dErrdifcombPt;
+       dErrVPOI += dYieldPtPOI*dYieldPtPOI*dErrdifcombPtPOI*dErrdifcombPtPOI;
       }
     }//end of for(Int_t b=0;b<iNbinsPt;b++)  
   } else { cout<<"fHistProFlow is NULL"<<endl; }
-  if (dSum != 0. ) {
-    dV /= dSum;  //because pt distribution should be normalised
-    dErrV /= (dSum*dSum);
-    dErrV = TMath::Sqrt(dErrV); 
+  if (dSumPOI != 0. ) {
+    dVPOI /= dSumPOI;  //because pt distribution should be normalised
+    dErrVPOI /= (dSumPOI*dSumPOI);
+    dErrVPOI = TMath::Sqrt(dErrVPOI); 
   }
-  cout<<"dV{MC} is "<<dV<<" +- "<<dErrV<<endl;
-  fCommonHistsRes->FillIntegratedFlow(dV,dErrV);//to be removed 
-  fCommonHistsRes->FillIntegratedFlowPOI(dV,dErrV);
+  cout<<"dV{MC} (POI) is "<<dVPOI<<" +- "<<dErrVPOI<<endl;
+
+  fCommonHistsRes->FillIntegratedFlowPOI(dVPOI,dErrVPOI);
   
   //Eta:
   if(fHistProDiffFlowEtaPOI)
   {
-   for(Int_t b=0;b<iNbinsEta;b++)
+   for(Int_t b=1;b<iNbinsEta;b++)
    {
-    dv2proEta = fHistProDiffFlowEtaPOI->GetBinContent(b+1);
-    dErrdifcombEta = fHistProDiffFlowEtaPOI->GetBinError(b+1);//to be improved (treatment of errors for non-Gaussian distribution needed!)
+    dv2proEtaPOI = fHistProDiffFlowEtaPOI->GetBinContent(b);
+    dErrdifcombEtaPOI = fHistProDiffFlowEtaPOI->GetBinError(b);//to be improved (treatment of errors for non-Gaussian distribution needed!)
     //fill common hist results:
-    fCommonHistsRes->FillDifferentialFlowEtaPOI(b+1, dv2proEta, dErrdifcombEta); 
+    fCommonHistsRes->FillDifferentialFlowEtaPOI(b, dv2proEtaPOI, dErrdifcombEtaPOI); 
    }
   }   
-                         
+  
+  cout<<endl;                    
   cout<<".....finished"<<endl;
 }
 
index 91a5860..e7a584d 100644 (file)
@@ -68,8 +68,8 @@ class AliFlowAnalysisWithMCEventPlane {
    TH1F*     GetHistRP()                  {return this->fHistRP; } 
    void      SetHistRP(TH1F* aHistRP)     {this->fHistRP = aHistRP; }
    
-   TProfile* GetHistProIntFlowRP()                            {return this->fHistProIntFlowRP; } 
-   void      SetHistProIntFlowRP(TProfile* aHistProIntFlowRP) {this->fHistProIntFlowRP = aHistProIntFlowRP; }
+   TProfile* GetHistProIntFlow()                          {return this->fHistProIntFlow; } 
+   void      SetHistProIntFlow(TProfile* aHistProIntFlow) {this->fHistProIntFlow = aHistProIntFlow; }
    
    TProfile* GetHistProDiffFlowPtRP()                               {return this->fHistProDiffFlowPtRP; } 
    void      SetHistProDiffFlowPtRP(TProfile* aHistProDiffFlowPtRP) {this->fHistProDiffFlowPtRP = aHistProDiffFlowPtRP; } 
@@ -105,7 +105,7 @@ class AliFlowAnalysisWithMCEventPlane {
    TProfile*    fHistProFlow;       //
    TH1F*        fHistRP;            //
    
-   TProfile*    fHistProIntFlowRP;     //profile used to calculate the integrated flow of RP particles
+   TProfile*    fHistProIntFlow;     //profile used to calculate the integrated flow of RP particles
    TProfile*    fHistProDiffFlowPtRP;  //profile used to calculate the differential flow (Pt) of RP particles 
    TProfile*    fHistProDiffFlowEtaRP; //profile used to calculate the differential flow (Eta) of RP particles 
    
index bd999f6..6cd3a1e 100644 (file)
@@ -159,12 +159,12 @@ void AliQCumulantsFunctions::Calculate()
  Double_t two   = fQCorr->GetBinContent(1);  //<<2>>_{n|n}
  Double_t four  = fQCorr->GetBinContent(11); //<<4>>_{n,n|n,n}
  Double_t six   = fQCorr->GetBinContent(24); //<<6>>_{n,n,n|n,n,n}
- Double_t eight = fQCorr->GetBinContent(31);   //<<8>>_{n,n,n,n|n,n,n,n}
+ Double_t eight = fQCorr->GetBinContent(31); //<<8>>_{n,n,n,n|n,n,n,n}
   
  //errors of 2-, 4- and 6-particle azimuthal correlation:
  Double_t twoErr   = fQCorr->GetBinError(1);  //sigma_{<<2>>_{n|n}} 
  Double_t fourErr  = fQCorr->GetBinError(11); //sigma_{<<4>>_{n,n|n,n}} 
- Double_t sixErr   = fQCorr->GetBinError(21); //sigma_{<<6>>_{n,n,n|n,n,n}}
+ Double_t sixErr   = fQCorr->GetBinError(24); //sigma_{<<6>>_{n,n,n|n,n,n}}
  //Double_t eightErr = fQCorr->GetBinError(); //sigma_{<<8>>_{n,n,n,n|n,n,n,n}}
  
  //covariances of multi-particle correlations
@@ -208,7 +208,7 @@ void AliQCumulantsFunctions::Calculate()
  cout<<endl;
  cout<<"**************************************"<<endl;
  cout<<"**************************************"<<endl;
- cout<<"flow estimates from Q-cumulants (RP):"<<endl;
+ cout<<"flow estimates from Q-cumulants :"<<endl;
  cout<<endl;
  
  Double_t vn2=0.,vn4=0.,vn6=0.,vn8=0.;
@@ -220,9 +220,8 @@ void AliQCumulantsFunctions::Calculate()
   fIntRes->SetBinContent(1,vn2);
   fIntRes->SetBinError(1,sd2);
   //common histograms:
-  fchr2nd->FillIntegratedFlow(vn2,sd2);//to be removed
+  fchr2nd->FillIntegratedFlow(vn2,sd2);
   fchr2nd->FillChi(vn2*pow(AvM,0.5));//to be removed
-  fchr2nd->FillIntegratedFlowRP(vn2,sd2);
  }else{
   cout<<" v_"<<n<<"{2} = Im"<<endl;
  }          
@@ -233,9 +232,8 @@ void AliQCumulantsFunctions::Calculate()
   fIntRes->SetBinContent(2,vn4);
   fIntRes->SetBinError(2,sd4);
   //common histograms:
-  fchr4th->FillIntegratedFlow(vn4,sd4);//to be removed
+  fchr4th->FillIntegratedFlow(vn4,sd4);
   fchr4th->FillChi(vn4*pow(AvM,0.5));//to be removed
-  fchr4th->FillIntegratedFlowRP(vn4,sd4);
  }else{
   cout<<" v_"<<n<<"{4} = Im"<<endl;
  }
@@ -246,9 +244,8 @@ void AliQCumulantsFunctions::Calculate()
   fIntRes->SetBinContent(3,vn6);
   fIntRes->SetBinError(3,sd6);
   //common histograms:
-  fchr6th->FillIntegratedFlow(vn6,sd6);//to be removed
+  fchr6th->FillIntegratedFlow(vn6,sd6);
   fchr6th->FillChi(vn6*pow(AvM,0.5));//to be removed
-  fchr6th->FillIntegratedFlowRP(vn6,sd6);
  }else{
   cout<<" v_"<<n<<"{6} = Im"<<endl;
  }
@@ -258,9 +255,8 @@ void AliQCumulantsFunctions::Calculate()
   fIntRes->SetBinContent(4,vn8);
   fIntRes->SetBinError(4,sd8);
   //common histograms:
-  fchr8th->FillIntegratedFlow(vn8,sd8);//to be removed
+  fchr8th->FillIntegratedFlow(vn8,sd8);
   fchr8th->FillChi(vn8*pow(AvM,0.5));//to be removed
-  fchr8th->FillIntegratedFlowRP(vn8,sd8);
  }else{
   cout<<" v_"<<n<<"{8} = Im"<<endl;
  }
@@ -271,39 +267,95 @@ void AliQCumulantsFunctions::Calculate()
  cout<<endl; 
 //--------------------------------------------------------------------------------------------------------- 
  
-//---------------------------------------------------------------------------------------------------------    
+
+
+//---------------------------------------------------------------------------------------------------------
 //differential flow (RP)
+Int_t nBinsPtRP = fbinPt2p1n1nRP->GetNbinsX();
+Int_t nBinsEtaRP = fbinEta2p1n1nRP->GetNbinsX();
 
 //Pt:
 Double_t secondOrderQCumulantDiffFlowPtRP = 0.;
 Double_t fourthOrderQCumulantDiffFlowPtRP = 0.;
 
-Int_t nBinsPtRP = fbinPt2p1n1nRP->GetNbinsX();
+Double_t dVn2ndRP=0.,dSd2ndRP=0.,dDiffvn2ndRP=0.,dYield2ndRP=0.,dSum2ndRP=0.;
+Double_t dVn4thRP=0.,dSd4thRP=0.,dDiffvn4thRP=0.,dYield4thRP=0.,dSum4thRP=0.;
 
 for(Int_t bb=1;bb<nBinsPtRP+1;bb++)
 {
+ //QC{2}
  if(fbinPt2p1n1nRP->GetBinEntries(bb)>0.&&vn2!=0)
  {
+  //cout<<"bin = "<<bb<<" : "<<(fch2nd->GetHistPtDiff())->GetBinContent(bb)<<endl;
+  //cout<<endl;
   secondOrderQCumulantDiffFlowPtRP = fbinPt2p1n1nRP->GetBinContent(bb);
   fDiffRes2nd->SetBinContent(bb,secondOrderQCumulantDiffFlowPtRP/vn2);
   //common histogram:
-  fchr2nd->FillDifferentialFlowPtRP(bb,secondOrderQCumulantDiffFlowPtRP/vn2, 0.);//to be improved (errors)
+  fchr2nd->FillDifferentialFlowPtRP(bb,secondOrderQCumulantDiffFlowPtRP/vn2, 0.);//to be improved (errors) && bb or bb+1
+  //-------------------------------------------------------------
+  //integrated flow (RP, Pt, 2nd order):
+  dDiffvn2ndRP=(fchr2nd->GetHistDiffFlowPtRP())->GetBinContent(bb);
+  dYield2ndRP=(fch2nd->GetHistPtInt())->GetBinContent(bb);
+  dVn2ndRP+=dDiffvn2ndRP*dYield2ndRP;
+  dSum2ndRP+=dYield2ndRP;
+  //-------------------------------------------------------------
  }
+ //QC{4]
  if(fbinPt4p1n1n1n1nRP->GetBinEntries(bb)>0.&&vn4!=0.)
  {
   fourthOrderQCumulantDiffFlowPtRP = fbinPt4p1n1n1n1nRP->GetBinContent(bb)-2.*fbinPt2p1n1nRP->GetBinContent(bb)*pow(vn2,2.);
   fDiffRes4th->SetBinContent(bb,-1.*fourthOrderQCumulantDiffFlowPtRP/pow(vn4,3.));
   //common histogram:
   fchr4th->FillDifferentialFlowPtRP(bb,-1.*fourthOrderQCumulantDiffFlowPtRP/pow(vn4,3.), 0.);//to be improved (errors)
+  //-------------------------------------------------------------
+  //integrated flow (RP, Pt, 4th order):
+  dDiffvn4thRP=(fchr4th->GetHistDiffFlowPtRP())->GetBinContent(bb);
+  dYield4thRP=(fch4th->GetHistPtInt())->GetBinContent(bb);
+  dVn4thRP+=dDiffvn4thRP*dYield4thRP;
+  dSum4thRP+=dYield4thRP;
+  //-------------------------------------------------------------
  }
-}      
+}  
+
+//integrated flow estimates from Q-cumulants (RP):
+cout<<endl;
+cout<<"**************************************"<<endl;
+cout<<"**************************************"<<endl;
+cout<<"flow estimates from Q-cumulants (RP):"<<endl;
+cout<<endl;
+//storing the final results for integrated flow (RP):
+//QC{2}
+if(dSum2ndRP && fchr2nd)
+{
+ dVn2ndRP/=dSum2ndRP;
+ fchr2nd->FillIntegratedFlowRP(dVn2ndRP,0.);//to be improved (errors)
+ cout<<" v_"<<n<<"{2} = "<<dVn2ndRP<<" +/- "<<dSd2ndRP<<endl;
+}else 
+ {
+  cout<<" v_"<<n<<"{2} = Im"<<endl;
+ }
+
+//QC{4}
+if(dSum4thRP && fchr4th)
+{
+ dVn4thRP/=dSum4thRP;
+ fchr4th->FillIntegratedFlowRP(dVn4thRP,0.);//to be improved (errors)
+ cout<<" v_"<<n<<"{4} = "<<dVn4thRP<<" +/- "<<dSd4thRP<<endl;
+}else
+ {
+  cout<<" v_"<<n<<"{4} = Im"<<endl;
+ }
+
+cout<<endl;
+cout<<"   nEvts = "<<nEvts<<", AvM = "<<AvM<<endl;
+cout<<"**************************************"<<endl;
+cout<<"**************************************"<<endl;
+cout<<endl;  
 
 //Eta:
 Double_t secondOrderQCumulantDiffFlowEtaRP = 0.;
 Double_t fourthOrderQCumulantDiffFlowEtaRP = 0.;
 
-Int_t nBinsEtaRP = fbinEta2p1n1nRP->GetNbinsX();
-
 for(Int_t bb=1;bb<nBinsEtaRP+1;bb++)
 {
  if(fbinEta2p1n1nRP->GetBinEntries(bb)>0.&&vn2!=0)
@@ -320,8 +372,13 @@ for(Int_t bb=1;bb<nBinsEtaRP+1;bb++)
   //common histogram:
   fchr4th->FillDifferentialFlowEtaRP(bb,-1.*fourthOrderQCumulantDiffFlowEtaRP/pow(vn4,3.), 0.);//to be improved (errors)
  }
-}      
-
+}    
+//---------------------------------------------------------------------------------------------------       
+        
+        
+        
+                      
+//---------------------------------------------------------------------------------------------------------
 //differential flow (POI)
 Int_t nBinsPtPOI = fbinPt2p1n1nPOI->GetNbinsX();
 Int_t nBinsEtaPOI = fbinEta2p1n1nPOI->GetNbinsX();
@@ -330,10 +387,8 @@ Int_t nBinsEtaPOI = fbinEta2p1n1nPOI->GetNbinsX();
 Double_t secondOrderQCumulantDiffFlowPtPOI = 0.;
 Double_t fourthOrderQCumulantDiffFlowPtPOI = 0.;
 
-Double_t dVn2ndPOI=0.,dSd2ndPOI=0.,dDiffvn2nd=0.,dYield2nd=0.,dSum2nd=0.;
-Double_t dVn4thPOI=0.,dSd4thPOI=0.,dDiffvn4th=0.,dYield4th=0.,dSum4th=0.;
-
-Double_t dVn2ndPOIEta=0.,dDiffvn2ndEta=0.,dYield2ndEta=0.,dSum2ndEta=0.;//to be removed
+Double_t dVn2ndPOI=0.,dSd2ndPOI=0.,dDiffvn2ndPOI=0.,dYield2ndPOI=0.,dSum2ndPOI=0.;
+Double_t dVn4thPOI=0.,dSd4thPOI=0.,dDiffvn4thPOI=0.,dYield4thPOI=0.,dSum4thPOI=0.;
 
 for(Int_t bb=1;bb<nBinsPtPOI+1;bb++)
 {
@@ -348,10 +403,10 @@ for(Int_t bb=1;bb<nBinsPtPOI+1;bb++)
   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;
+  dDiffvn2ndPOI=(fchr2nd->GetHistDiffFlowPtPOI())->GetBinContent(bb);
+  dYield2ndPOI=(fch2nd->GetHistPtDiff())->GetBinContent(bb);
+  dVn2ndPOI+=dDiffvn2ndPOI*dYield2ndPOI;
+  dSum2ndPOI+=dYield2ndPOI;
   //-------------------------------------------------------------
  }
  //QC{4]
@@ -363,16 +418,14 @@ for(Int_t bb=1;bb<nBinsPtPOI+1;bb++)
   fchr4th->FillDifferentialFlowPtPOI(bb,-1.*fourthOrderQCumulantDiffFlowPtPOI/pow(vn4,3.), 0.);//to be improved (errors)
   //-------------------------------------------------------------
   //integrated flow (POI, Pt, 4th order):
-  dDiffvn4th=(fchr4th->GetHistDiffFlowPtPOI())->GetBinContent(bb);
-  dYield4th=(fch4th->GetHistPtDiff())->GetBinContent(bb);
-  dVn4thPOI+=dDiffvn4th*dYield4th;
-  dSum4th+=dYield4th;
+  dDiffvn4thPOI=(fchr4th->GetHistDiffFlowPtPOI())->GetBinContent(bb);
+  dYield4thPOI=(fch4th->GetHistPtDiff())->GetBinContent(bb);
+  dVn4thPOI+=dDiffvn4thPOI*dYield4thPOI;
+  dSum4thPOI+=dYield4thPOI;
   //-------------------------------------------------------------
  }
 }      
 
-//---------------------------------------------------------------------------------------------------
-//integrated flow estimates from Q-cumulants (POI):
 cout<<endl;
 cout<<"**************************************"<<endl;
 cout<<"**************************************"<<endl;
@@ -380,9 +433,9 @@ cout<<"flow estimates from Q-cumulants (POI):"<<endl;
 cout<<endl;
 //storing the final results for integrated flow (POI):
 //QC{2}
-if(dSum2nd && fchr2nd)
+if(dSum2ndPOI && fchr2nd)
 {
- dVn2ndPOI/=dSum2nd;
+ dVn2ndPOI/=dSum2ndPOI;
  fchr2nd->FillIntegratedFlowPOI(dVn2ndPOI,0.);//to be improved (errors)
  cout<<" v_"<<n<<"{2} = "<<dVn2ndPOI<<" +/- "<<dSd2ndPOI<<endl;
 }else 
@@ -391,9 +444,9 @@ if(dSum2nd && fchr2nd)
  }
 
 //QC{4}
-if(dSum4th && fchr4th)
+if(dSum4thPOI && fchr4th)
 {
- dVn4thPOI/=dSum4th;
+ dVn4thPOI/=dSum4thPOI;
  fchr4th->FillIntegratedFlowPOI(dVn4thPOI,0.);//to be improved (errors)
  cout<<" v_"<<n<<"{4} = "<<dVn4thPOI<<" +/- "<<dSd4thPOI<<endl;
 }else
@@ -406,8 +459,6 @@ cout<<"   nEvts = "<<nEvts<<", AvM = "<<AvM<<endl;
 cout<<"**************************************"<<endl;
 cout<<"**************************************"<<endl;
 cout<<endl;  
-//---------------------------------------------------------------------------------------------------
-
 
 //Eta:
 Double_t secondOrderQCumulantDiffFlowEtaPOI = 0.;
@@ -462,475 +513,6 @@ 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.);
- Double_t second=pow(fourErr,2.);
- Double_t third=-8.*two*cov24;
- cout<<endl;
- cout<<endl;
- cout<<endl;
- cout<<"       1st = "<<(1./16.)*first/(pow(-fourthOrderQCumulant,1.5))<<endl;
- cout<<"       2nd = "<<(1./16.)*second/(pow(-fourthOrderQCumulant,1.5))<<endl;
- cout<<"       3rd = "<<(1./16.)*third/(pow(-fourthOrderQCumulant,1.5))<<endl;
- cout<<endl;
- cout<<"   nEvts = "<<nEvts<<", AvM = "<<AvM<<endl;
- cout<<endl;
- cout<<endl;
- cout<<"c-c = "<<pow((1./16.)*((first+second+third)/(pow(-fourthOrderQCumulant,1.5))),0.5)<<endl; 
- */
-           
-             
- /*
- //<Q-vector components>
- Double_t AvQx  = fQVector->GetBinContent(1); //<Q_x>
- Double_t AvQy  = fQVector->GetBinContent(2); //<Q_y>
- Double_t AvQ2x = fQVector->GetBinContent(3); //<(Q_x)^2>
- Double_t AvQ2y = fQVector->GetBinContent(4); //<(Q_y)^2>
- Double_t AvQ2  = fQVector->GetBinContent(5); //<|Q|^2>
- Double_t AvQ4  = fQVector->GetBinContent(6); //<|Q|^4>
- Double_t AvQ_2n2 = fQVector->GetBinContent(7); //<|Q_2n|^2>
- Double_t ReQ2nQnstarQnstar = fQVector->GetBinContent(8); //<Re[Q_n^2 * Q_2n^*]> 
- Double_t ImQ2nQnstarQnstar = fQVector->GetBinContent(9); //<Im[Q_n^2 * Q_2n^*]>
- Double_t AvQ_3n2 = fQVector->GetBinContent(10); //<|Q_3n|^2>
- Double_t AvQ_4n2 = fQVector->GetBinContent(11); //<|Q_4n|^2>
- Double_t AvQ_5n2 = fQVector->GetBinContent(12); //<|Q_5n|^2>
- Double_t ReQ3nQ2nstarQnstar = fQVector->GetBinContent(13); //Re<Q_{3n} Q_{2n}^* Q_{n}^*>
- Double_t ReQ3nQnstarQnstarQnstar = fQVector->GetBinContent(14); //Re<Q_{3n} Q_{n}^* Q_{n}^* Q_{n}^*>
- Double_t HereQ2nQnQ2nstarQnstar = fQVector->GetBinContent(15); //<|Q_{2n}|^2 |Q_{n}|^2>
- Double_t HereQ2nQnQnstarQnstarQnstar = fQVector->GetBinContent(16); //Re<Q_{2n} Q_{n} Q_{n}^* Q_{n}^* Q_{n}^*>
- Double_t AvQ6 = fQVector->GetBinContent(17); //<|Q|^6>
-*/
- /*
- //direct particle correlations
- Double_t twoDirect_n_n = fDirect->GetBinContent(1);
- Double_t fourDirect = fDirect->GetBinContent(2);
- Double_t threeDirect_2n_n_n = fDirect->GetBinContent(3);
- Double_t threeDirect_3n_2n_n = fDirect->GetBinContent(4);
- Double_t fourDirect_3n_n_n_n = fDirect->GetBinContent(5);
- Double_t fourDirect_2n_n_2n_n = fDirect->GetBinContent(6);
- Double_t fiveDirect_2n_n_n_n_n = fDirect->GetBinContent(7);
- Double_t sixDirect_n_n_n_n_n_n = fDirect->GetBinContent(8);
- Double_t twoDirect_2n_2n = fDirect->GetBinContent(9);
- Double_t twoDirect_3n_3n = fDirect->GetBinContent(10);
- Double_t twoDirect_4n_4n = fDirect->GetBinContent(11);
- Double_t twoDirect_5n_5n = fDirect->GetBinContent(12);
-*/ 
-/* 
- //CORRELATIONS
- //two particle correlations
- Double_t two_n_n=(AvQ2-AvM)/(AvM*(AvM-1));
- Double_t two_2n_2n=(AvQ_2n2-AvM)/(AvM*(AvM-1));
- Double_t two_3n_3n=(AvQ_3n2-AvM)/(AvM*(AvM-1));
- Double_t two_4n_4n=(AvQ_4n2-AvM)/(AvM*(AvM-1));
- Double_t two_5n_5n=(AvQ_5n2-AvM)/(AvM*(AvM-1));
-
-
-
-
-
- //four particle correlations:
- //<4>_{n,n|n,n}
- Double_t four_n_n_n_n = (AvQ4+AvQ_2n2-4.*(AvM-2.)*AvQ2-2.*ReQ2nQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))+2./((AvM-1)*(AvM-2));
- //<4>_{3n|n,n,n}
- //Double_t four_3n_n_n_n=(ReQ3nQnstarQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))-3.*(three_3n_2n_n+three_2n_n_n)/(AvM-3.)-(3.*two+3.*two_2n+two_3n)/((AvM-2.)*(AvM-3.))-1./((AvM-1.)*(AvM-2.)*(AvM-3.)); OK!!!
- //Double_t four_3n_n_n_n_=(ReQ3nQnstarQnstarQnstar-3.*ReQ3nQ2nstarQnstar-3.*ReQ2nQnstarQnstar)/(AvM*(AvM-1.)*(AvM-2.)*(AvM-3.))+ (2.*AvQ_3n2 + 3.*AvQ_2n2 + 6.*AvQ2-6.*AvM ) /(AvM*(AvM-1.)*(AvM-2.)*(AvM-3.)); OK!!! (final version)
-
-
- //<4>_{2n,n|2n,n}
- //Double_t four_2n_n_2n_n = (HereQ2nQnQ2nstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))-2.*(three_2n_n_n+three_3n_2n_n)/(AvM-3.)-((AvM+1.)*two+AvM*two_2n+two_3n)/((AvM-2.)*(AvM-3.))-(AvM)/((AvM-1)*(AvM-2)*(AvM-3)); OK!!!
- //Double_t four_2n_n_2n_n_Alternative = (HereQ2nQnQ2nstarQnstar-2.*ReQ3nQ2nstarQnstar-2.*ReQ2nQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))-((AvM-5)*AvQ2+(AvM-4)*AvQ_2n2-AvQ_3n2)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))+(AvM-6)/((AvM-1)*(AvM-2)*(AvM-3)); OK!!! (final version)
-  
- //<6>_{n,n,n|n,n,n} 
- //Double_t six_n_n_n_n_n_n=AvQ6/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) - 3.*(2.*five_2n_n_n_n_n)/(AvM-5) - (2.*four_3n_n_n_n+9.*(AvM-2.)*four+9.*four_2n_n_2n_n)/((AvM-4)*(AvM-5)) - 3.*(2.*(3.*AvM-5.)*three_2n_n_n+2.*three_3n_2n_n)/((AvM-3)*(AvM-4)*(AvM-5)) - ((18.*AvM*AvM-45.*AvM+33.)*two+3.*(3.*AvM-4.)*two_2n+two_3n)/((AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) - (6.*AvM*AvM-9.*AvM+4.)/((AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)); OK!!!
- //Double_t six_n_n_n_n_n_n = (AvQ6+9.*HereQ2nQnQ2nstarQnstar-6.*HereQ2nQnQnstarQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) + 4.*(ReQ3nQnstarQnstarQnstar-3.*ReQ3nQ2nstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) + 2.*(9.*(AvM-4.)*ReQ2nQnstarQnstar+2.*AvQ_3n2)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) - 9.*(AvQ4+AvQ_2n2)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-5)) + (18.*AvQ2 )/(AvM*(AvM-1)*(AvM-3)*(AvM-4)) - (6.)/((AvM-1)*(AvM-2)*(AvM-3)); OK!!! (final version)
-      
-       
-        
-
-
-*/
-
- /*
- cout<<"<4>_{3n,n,n,n} correlation from Q-vector = "<<four_3n_n_n_n<<endl;AvQ2
- cout<<"<4>_{3n,n,n,n} correlation directly      = "<<fourDirect_3n_n_n_n<<endl;
- cout<<endl;
- cout<<"<4>_{2n,n,2n,n} correlation from Q-vector = "<<four_2n_n_2n_n<<endl;
- cout<<"<4>_{2n,n,2n,n} correlation directly      = "<<fourDirect_2n_n_2n_n<<endl;
- cout<<endl;
- cout<<"<5>_{2n,n,n,n,n} correlation from Q-vector = "<<five_2n_n_n_n_n<<endl;
- cout<<"<5>_{2n,n,n,n,n} correlation directly      = "<<fiveDirect_2n_n_n_n_n<<endl;
- cout<<endl;
- cout<<"<6>_{n,n,n,n,n,n} correlation from Q-vector = "<<six_n_n_n_n_n_n<<endl;
- cout<<"<6>_{n,n,n,n,n,n} correlation directly      = "<<sixDirect_n_n_n_n_n_n<<endl;
- cout<<endl;
- */
-
-  /*
-   
-//Q-CUMULANTS
-Double_t fourthOrderQCumulant = (AvQ4+AvQ_2n2-2.*ReQ2nQnstarQnstar)/(AvM*(AvM-1.)*(AvM-2.)*(AvM-3.)) - (2.*AvQ2*AvQ2)/(AvM*AvM*(AvM-1.)*(AvM-1.)) - (8.*AvQ2)/(AvM*(AvM-1.)*(AvM-1.)*(AvM-3.)) + 2./((AvM-1.)*(AvM-1.)*(AvM-2.));    
-     
-Double_t sixthOrderQCumulant = (AvQ6+9.*HereQ2nQnQ2nstarQnstar-6.*HereQ2nQnQnstarQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) + 4.*(ReQ3nQnstarQnstarQnstar - 3.*ReQ3nQ2nstarQnstar+AvQ_3n2)/(AvM*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)*(AvM-5)) + 18.*(4.*AvM+(AvM-5.)*AvQ2)*ReQ2nQnstarQnstar/(AvM*AvM*(AvM-1)*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-5)) - 9.*(4.*AvM+(AvM-5.)*AvQ2)*(AvQ4+AvQ_2n2)/(AvM*AvM*(AvM-1)*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-5)) + 36.*(5.*AvM-11.)*AvQ2/(AvM*(AvM-1)*(AvM-1)*(AvM-1)*(AvM-2)*(AvM-3)*(AvM-4)) + 12.*AvQ2*AvQ2*(6.*AvM-3.*AvQ2+AvM*AvQ2)/(AvM*AvM*AvM*(AvM-1)*(AvM-1)*(AvM-1)*(AvM-3)) - 24./((AvM-1)*(AvM-1)*(AvM-1)*(AvM-2)*(AvM-3));
-          
-           
-                           
-                
-                 
-                   
- //cout<<endl;
- //cout<<"should be the same? "<<threeDirect_3n_2n_n<<" "<<three_2n_n_nTemp<<endl;
- //cout<<endl;
- //cout<<endl;
- cout<<"***************************"<<endl;
- cout<<endl;
- cout<<"multiplicity = "<<AvM<<endl;
- cout<<endl;
- Double_t two=(AvQ2-AvM)/(AvM*(AvM-1));
- Double_t two_2n=(AvQ_2n2-AvM)/(AvM*(AvM-1));
- Double_t two_3n=(AvQ_3n2-AvM)/(AvM*(AvM-1));
- Double_t two_4n=(AvQ_4n2-AvM)/(AvM*(AvM-1));
- Double_t two_5n=(AvQ_5n2-AvM)/(AvM*(AvM-1));
-
- //cout<<"two's = "<<two<<" "<<two_2n<<" "<<two_3n<<" "<<two_4n<<" "<<two_5n<<endl;
- cout<<endl;
- Double_t four=(2.*AvM*(AvM-3.)+AvQ4-4.*(AvM-2.)*AvQ2-2.*ReQ2nQnstarQnstar+AvQ_2n2)/(AvM*(AvM-1)*(AvM-2)*(AvM-3));
- Double_t three_2n_n_n=(ReQ2nQnstarQnstar-AvM-2*AvM*(AvM-1)*two-AvM*(AvM-1)*two_2n)/(AvM*(AvM-1)*(AvM-2));
- //<3>_{3n,2n,n}  
- Double_t three_3n_2n_n=(ReQ3nQ2nstarQnstar-AvM*(AvM-1)*(two+two_2n+two_3n)-AvM)/(AvM*(AvM-1)*(AvM-2));
- //<4>_{3n,n,n,n}  
- //Double_t four_3n_n_n_n=(ReQ3nQnstarQnstarQnstar-3.*AvM*(AvM-1)*(AvM-2)*(three_2n_n_n+three_3n_2n_n)-AvM*(AvM-1)*(3.*two+3.*two_2n+two_3n)-AvM)/(AvM*(AvM-1)*(AvM-2)*(AvM-3));
- Double_t four_3n_n_n_n=(ReQ3nQnstarQnstarQnstar)/(AvM*(AvM-1)*(AvM-2)*(AvM-3))-3.*(three_3n_2n_n+three_2n_n_n)/(AvM-3.)-(3.*two+3.*two_2n+two_3n)/((AvM-2.)*(AvM-3.))-1./((AvM-1.)*(AvM-2.)*(AvM-3.));
-
- //<4>_{2n,n,2n,n}  
- Double_t four_2n_n_2n_n=(HereQ2nQnQ2nstarQnstar-2.*AvM*(AvM-1)*(AvM-2)*(three_2n_n_n+three_3n_2n_n)-AvM*(AvM-1)*((AvM+1)*two+AvM*two_2n+two_3n)-AvM*AvM)/(AvM*(AvM-1)*(AvM-2)*(AvM-3));
-  
-
-  
-    
-      
-  
- cout<<"<2> correlation from Q-vector           = "<<two<<endl;
- cout<<"<2> correlation directly                = "<<twoDirect_n_n<<endl;
- cout<<endl;
- cout<<"<4> correlation from Q-vector           = "<<four<<endl;
- cout<<"<4> correlation directly                = "<<fourDirect<<endl;
- cout<<endl;
- cout<<"<3>_{2n,n,n} correlation from Q-vector  = "<<three_2n_n_n<<endl;
- cout<<"<3>_{2n,n,n} correlation directly       = "<<threeDirect_2n_n_n<<endl;
- cout<<endl;
- cout<<"<3>_{3n,2n,n} correlation from Q-vector = "<<three_3n_2n_n<<endl;
- cout<<"<3>_{3n,2n,n} correlation directly      = "<<threeDirect_3n_2n_n<<endl;
- cout<<endl;
- cout<<"<4>_{3n,n,n,n} correlation from Q-vector = "<<four_3n_n_n_n<<endl;
- cout<<"<4>_{3n,n,n,n} correlation directly      = "<<fourDirect_3n_n_n_n<<endl;
- cout<<endl;
- cout<<"<4>_{2n,n,2n,n} correlation from Q-vector = "<<four_2n_n_2n_n<<endl;
- cout<<"<4>_{2n,n,2n,n} correlation directly      = "<<fourDirect_2n_n_2n_n<<endl;
- cout<<endl;
- cout<<"<5>_{2n,n,n,n,n} correlation from Q-vector = "<<five_2n_n_n_n_n<<endl;
- cout<<"<5>_{2n,n,n,n,n} correlation directly      = "<<fiveDirect_2n_n_n_n_n<<endl;
- cout<<endl;
- cout<<"<6>_{n,n,n,n,n,n} correlation from Q-vector = "<<six_n_n_n_n_n_n<<endl;
- cout<<"<6>_{n,n,n,n,n,n} correlation directly      = "<<sixDirect_n_n_n_n_n_n<<endl;
- cout<<endl;
-
-
-
-
- //have so far: <2>, <2>_{2n}, <3>_{2n,n,n}, <4>
-
-
- //cout<<"temp re = "<<ReQ2nQnstarQnstar<<endl;
- //cout<<"temp im = "<<ImQ2nQnstarQnstar<<endl;
- //cout<<AvQ2<<" "<<sqrt(AvQ4)<<endl;
- cout<<endl;
-
- if(AvQ2>AvM){
-  cout<<"v_2{2} = "<<100*sqrt(two)<<"%"<<endl;
- }else{
-  cout<<"v_2{2} = Im"<<endl;
- }
-
- if(four>0){
-  cout<<"v_2{4} = "<<100*pow(four,1./4.)<<"%"<<endl;
- }else{
-  cout<<"v_2{4} = Im"<<endl;
- }
-
-
-
-
- cout<<"AvM = "<<AvM<<endl;
- cout<<endl; 
- cout<<"*************************************"<<endl;
- cout<<"*************************************"<<endl;
- cout<<"flow estimates from Q-cumulants:"<<endl;
- Double_t cumulant2Q=two;
- Double_t cumulant4Q=four-2.*two*two;
- Double_t cumulant6Q=12.*two*two*two-9.*two*four+six_n_n_n_n_n_n;
-
- //cout<<endl;
- if(cumulant2Q>0.){
-  cout<<"v_2{2} = "<<100*pow(cumulant2Q,1./2.)<<"%"<<endl;
- }else{
-  cout<<"v_2{2} = Im"<<endl;
- }
- if(cumulant4Q<0.){
-  cout<<"v_2{4} = "<<100*pow(-cumulant4Q,1./4.)<<"%"<<endl;
- }else{
-  cout<<"v_2{4} = Im"<<endl;
- }
- if(cumulant6Q>0.){
-  cout<<"v_2{6} = "<<100*pow((1./4.)*cumulant6Q,1./6.)<<"%"<<endl;
- }else{
-  cout<<"v_2{6} = Im"<<endl;
- }
- cout<<"*************************************"<<endl;
- cout<<"*************************************"<<endl;
- cout<<endl;
- cout<<endl; 
- cout<<"*************************************"<<endl;
- cout<<"*************************************"<<endl;
- cout<<"flow estimates from Q-cumulants"<<endl; 
- cout<<"without multiplicities fluctuations:"<<endl;
- Double_t cumulant2QnoFluct=fQCorr->GetBinContent(1);
- Double_t cumulant4QnoFluct=fQCorr->GetBinContent(11)-2.*fQCorr->GetBinContent(1)*fQCorr->GetBinContent(1);
- Double_t cumulant6QnoFluct=12.*fQCorr->GetBinContent(1)*fQCorr->GetBinContent(1)*fQCorr->GetBinContent(1) - 9.*fQCorr->GetBinContent(1)*fQCorr->GetBinContent(11)+fQCorr->GetBinContent(21);
-
- //cout<<endl;
- if(cumulant2QnoFluct>0.){
-  cout<<"v_2{2} = "<<100*pow(cumulant2QnoFluct,1./2.)<<"%"<<endl;
-  fIntRes->SetBinContent(1,100*pow(cumulant2QnoFluct,1./2.));
- }else{
-  cout<<"v_2{2} = Im"<<endl;
- }
- if(cumulant4QnoFluct<0.){
-  cout<<"v_2{4} = "<<100*pow(-cumulant4QnoFluct,1./4.)<<"%"<<endl;
-  fIntRes->SetBinContent(2,100*pow(-cumulant4QnoFluct,1./4.));
- }else{
-  cout<<"v_2{4} = Im"<<endl;
- }
- if(cumulant6QnoFluct>0.){
-  cout<<"v_2{6} = "<<100*pow((1./4.)*cumulant6QnoFluct,1./6.)<<"%"<<endl;
-  fIntRes->SetBinContent(3,100*pow((1./4.)*cumulant6QnoFluct,1./6.));
- }else{
-  cout<<"v_2{6} = Im"<<endl;
- }
- cout<<"*************************************"<<endl;
- cout<<"*************************************"<<endl;
- cout<<endl;
-
-  
-   
-    
- cout<<"*****************************************************************"<<endl;
- cout<<"*****************************************************************"<<endl;
- cout<<"direct (nested loops) correlations vs correlations from Q-vectors"<<endl;
- cout<<"for ARBITRARY multiplicity:"<<endl;
- cout<<endl;
- cout<<"avarage multiplicity = "<<AvM<<endl;
- cout<<endl;
- cout<<"<2>_{n|n} correlation from Q-vector         = "<<fQCorr->GetBinContent(1)<<endl;
- cout<<"<2>_{n|n} correlation nested loops          = "<<fDirect->GetBinContent(1)<<endl;
- cout<<endl;
- cout<<"<4>_{n,n|n,n} correlation from Q-vector     = "<<fQCorr->GetBinContent(11)<<endl;
- cout<<"<4>_{n,n|n,n} correlation nested loops      = "<<fDirect->GetBinContent(2)<<endl; 
- cout<<endl;
- cout<<"<6>_{n,n,n|n,n,n} correlation from Q-vector = "<<fQCorr->GetBinContent(21)<<endl;
- cout<<"<6>_{n,n,n|n,n,n} correlation nested loops  = "<<fDirect->GetBinContent(8)<<endl; 
- cout<<"*****************************************************************"<<endl;
- cout<<"*****************************************************************"<<endl;
-     
-      
-       
-        
-         
-           
- cout<<"****** BDO ******"<<endl;
- if(BDO4<0.){
-  cout<<"v_2{4} = "<<100*pow(-BDO4,1./4.)<<"%"<<endl;
- }else{
-  cout<<"v_2{4} = Im"<<endl;
- }
- if(BDO6>0.){
-  cout<<"v_2{6} = "<<100*pow((1./4.)*BDO6,1./6.)<<"%"<<endl;
- }else{
-  cout<<"v_2{6} = Im"<<endl;
- }
-
- cout<<endl;
- cout<<"fourth order Q-cumulant should be the same and BDO4? "<<fourthOrderQCumulant<<" "<<cumulant4Q<<" "<<BDO4<<endl;
- cout<<endl;
- cout<<endl;
- cout<<"sixth order Q-cumulant should be the same and BDO6? "<<sixthOrderQCumulant<<" "<<cumulant6Q<<" "<<BDO6<<endl;
- cout<<endl;
- cout<<"*************************************"<<endl;
- cout<<"*************************************"<<endl;
- cout<<endl;
-  cout<<"*********************"<<endl;
- cout<<"multiplicity = "<<AvM<<endl;
- cout<<"correlations: "<<endl;
- cout<<"direct = "<<twoDirect_n_n<<endl;
- cout<<"approx = "<<two_n_n<<endl;
- cout<<"weight = "<<fQCorr->GetBinContent(1)<<endl;
- cout<<"*********************"<<endl;
- cout<<"test = "<<fQCorr->GetBinContent(11)<<" "<<four<<endl;
-
-*/
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
 Bool_t bNestedLoopsResults=kFALSE;
 if(bNestedLoopsResults) 
 { 
index 0ddf454..444ff75 100644 (file)
@@ -334,8 +334,8 @@ void AliAnalysisTaskMCEventPlane::Terminate(Option_t *)
     TProfile *pHistProFlow = dynamic_cast<TProfile*> 
       (fListHistos->FindObject("FlowPro_VPt_MCEP"));
 
-    TProfile *pHistProIntFlowRP = dynamic_cast<TProfile*> 
-      (fListHistos->FindObject("fHistProIntFlowRP")); 
+    TProfile *pHistProIntFlow = dynamic_cast<TProfile*> 
+      (fListHistos->FindObject("fHistProIntFlow")); 
                                
     TProfile *pHistProDiffFlowPtRP = dynamic_cast<TProfile*> 
       (fListHistos->FindObject("fHistProDiffFlowPtRP")); 
@@ -349,11 +349,11 @@ void AliAnalysisTaskMCEventPlane::Terminate(Option_t *)
     TProfile *pHistProDiffFlowEtaPOI = dynamic_cast<TProfile*> 
       (fListHistos->FindObject("fHistProDiffFlowEtaPOI"));                             
 
-    if (pCommonHists && pCommonHistResults && pHistProFlow && pHistProIntFlowRP && pHistProDiffFlowPtRP && pHistProDiffFlowEtaRP && pHistProDiffFlowPtPOI && pHistProDiffFlowEtaPOI) {
+    if (pCommonHists && pCommonHistResults && pHistProFlow && pHistProIntFlow && pHistProDiffFlowPtRP && pHistProDiffFlowEtaRP && pHistProDiffFlowPtPOI && pHistProDiffFlowEtaPOI) {
       fMcTerm->SetCommonHists(pCommonHists);
       fMcTerm->SetCommonHistsRes(pCommonHistResults);
       fMcTerm->SetHistProFlow(pHistProFlow);
-      fMcTerm->SetHistProIntFlowRP(pHistProIntFlowRP);
+      fMcTerm->SetHistProIntFlow(pHistProIntFlow);
       fMcTerm->SetHistProDiffFlowPtRP(pHistProDiffFlowPtRP);      
       fMcTerm->SetHistProDiffFlowEtaRP(pHistProDiffFlowEtaRP);  
       fMcTerm->SetHistProDiffFlowPtPOI(pHistProDiffFlowPtPOI);      
similarity index 84%
rename from PWG2/FLOW/macros/CompareFlowResults.C
rename to PWG2/FLOW/macros/compareFlowResults.C
index d56a7db..ca64b24 100644 (file)
@@ -1,6 +1,6 @@
 //type of analysis can be: ESD, AOD, MC, ESDMC0, ESDMC1
 //const TString type = "ESD"; 
-void CompareFlowResults(TString type="ESD")
+void compareFlowResults(TString type="ESD")
 {
   //load needed libraries:
   gSystem->AddIncludePath("-I$ROOTSYS/include");
@@ -24,9 +24,9 @@ void CompareFlowResults(TString type="ESD")
   //             set here which plots will be shown by default
   //==================================================================================
   //RP = particles used to determine the reaction plane
-  Bool_t plotIntFlowRP = kFALSE;     //integrated flow RP
-  Bool_t plotDiffFlowPtRP = kFALSE;  //differential flow (Pt,RP)
-  Bool_t plotDiffFlowEtaRP = kFALSE; //differential flow (Eta,RP)
+  Bool_t plotIntFlowRP = kTRUE;     //integrated flow RP
+  Bool_t plotDiffFlowPtRP = kTRUE;  //differential flow (Pt,RP)
+  Bool_t plotDiffFlowEtaRP = kTRUE; //differential flow (Eta,RP)
   //POI = particle of interest
   Bool_t plotIntFlowPOI = kTRUE;     //integrated flow POI
   Bool_t plotDiffFlowPtPOI = kTRUE;   //differential flow (Pt,POI)
@@ -39,7 +39,6 @@ void CompareFlowResults(TString type="ESD")
   //==================================================================================
   //                         accessing output files
   //==================================================================================
-  
   //open the output files:
   TString inputFileNameMCEP = "outputMCEPanalysis";
   TFile* fileMCEP = NULL;
@@ -84,7 +83,7 @@ void CompareFlowResults(TString type="ESD")
   Int_t meshStyle = 1001;
   Int_t meshColor = kRed-10;
   
-  //marker style and color (int. flow) 
+  //marker style and color (plots for int. flow) 
   Int_t markerStyle = 21;
   Int_t markerColor = kRed-3;
   //==================================================================================
@@ -120,8 +119,8 @@ void CompareFlowResults(TString type="ESD")
   //booking the graph to store flow values and errors from all methods:  
   Double_t x[nMethods] = {0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5, 11.5};
   Double_t xError[nMethods] = {0.};
-  Double_t flowValue[nMethods] = {0.};//to be removed
-  Double_t flowError[nMethods] = {0.};//to be removed
+  Double_t flowValue[nMethods] = {0.};
+  Double_t flowError[nMethods] = {0.};
   Double_t flowValueRP[nMethods] = {0.};
   Double_t flowErrorRP[nMethods] = {0.};
   Double_t flowValuePOI[nMethods] = {0.};
@@ -138,8 +137,8 @@ void CompareFlowResults(TString type="ESD")
       mcepCommonHist    = dynamic_cast<AliFlowCommonHist*> (pListMCEP->FindObject("AliFlowCommonHistMCEP"));
       mcepCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListMCEP->FindObject("AliFlowCommonHistResultsMCEP"));
       if(mcepCommonHistRes) {
-       flowValue[0] = (mcepCommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
-       flowError[0] = (mcepCommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[0] = (mcepCommonHistRes->GetHistIntFlow())->GetBinContent(1);
+       flowError[0] = (mcepCommonHistRes->GetHistIntFlow())->GetBinError(1);
        flowValueRP[0] = (mcepCommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
        flowErrorRP[0] = (mcepCommonHistRes->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[0] = (mcepCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
@@ -158,8 +157,8 @@ void CompareFlowResults(TString type="ESD")
       lyz1CommonHist = dynamic_cast<AliFlowCommonHist*> (pListLYZ1->FindObject("AliFlowCommonHistLYZ1"));
       lyz1CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListLYZ1->FindObject("AliFlowCommonHistResultsLYZ1"));
       if(lyz1CommonHistRes) {
-       flowValue[10] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
-       flowError[10] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[10] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinContent(1);
+       flowError[10] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinError(1);
        flowValueRP[10] = (lyz1CommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
        flowErrorRP[10] = (lyz1CommonHistRes->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[10] = (lyz1CommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
@@ -190,8 +189,8 @@ void CompareFlowResults(TString type="ESD")
       lyzepCommonHist = dynamic_cast<AliFlowCommonHist*> (pListLYZEP->FindObject("AliFlowCommonHistLYZEP"));
       lyzepCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListLYZEP->FindObject("AliFlowCommonHistResultsLYZEP"));
       if(lyzepCommonHistRes) {
-       flowValue[11] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
-       //flowError[11] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[11] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinContent(1);
+       //flowError[11] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinError(1);
        flowValueRP[11] = (lyzepCommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
        //flowErrorRP[11] = (lyzepCommonHistRes->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[11] = (lyzepCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
@@ -210,8 +209,8 @@ void CompareFlowResults(TString type="ESD")
       fqdCommonHist = dynamic_cast<AliFlowCommonHist*> (pListFQD->FindObject("AliFlowCommonHistFQD"));
       fqdCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListFQD->FindObject("AliFlowCommonHistResultsFQD"));
       if(fqdCommonHistRes) {
-       flowValue[9] = (fqdCommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
-       flowError[9] = (fqdCommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[9] = (fqdCommonHistRes->GetHistIntFlow())->GetBinContent(1);
+       flowError[9] = (fqdCommonHistRes->GetHistIntFlow())->GetBinError(1);
        flowValueRP[9] = (fqdCommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
        flowErrorRP[9] = (fqdCommonHistRes->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[9] = (fqdCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
@@ -233,8 +232,8 @@ void CompareFlowResults(TString type="ESD")
       gfcCommonHist = dynamic_cast<AliFlowCommonHist*> (pListGFC->FindObject("AliFlowCommonHistGFC"));
       gfcCommonHistRes2 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults2ndOrderGFC"));
       if(gfcCommonHistRes2) {
-       flowValue[1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);//to be removed
-       flowError[1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);
+       flowError[1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinError(1);
        flowValueRP[1] = (gfcCommonHistRes2->GetHistIntFlowRP())->GetBinContent(1);
        flowErrorRP[1] = (gfcCommonHistRes2->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[1] = (gfcCommonHistRes2->GetHistIntFlowPOI())->GetBinContent(1);
@@ -242,8 +241,8 @@ void CompareFlowResults(TString type="ESD")
       }
       gfcCommonHistRes4 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults4thOrderGFC"));
       if(gfcCommonHistRes4) {
-       flowValue[3] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinContent(1);//to be removed
-       flowError[3] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[3] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinContent(1);
+       flowError[3] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinError(1);
        flowValueRP[3] = (gfcCommonHistRes4->GetHistIntFlowRP())->GetBinContent(1);
        flowErrorRP[3] = (gfcCommonHistRes4->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[3] = (gfcCommonHistRes4->GetHistIntFlowPOI())->GetBinContent(1);
@@ -251,8 +250,8 @@ void CompareFlowResults(TString type="ESD")
       }
       gfcCommonHistRes6 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults6thOrderGFC"));
       if(gfcCommonHistRes6) {
-       flowValue[5] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinContent(1);//to be removed
-       flowError[5] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[5] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinContent(1);
+       flowError[5] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinError(1);
        flowValueRP[5] = (gfcCommonHistRes6->GetHistIntFlowRP())->GetBinContent(1);
        flowErrorRP[5] = (gfcCommonHistRes6->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[5] = (gfcCommonHistRes6->GetHistIntFlowPOI())->GetBinContent(1);
@@ -260,8 +259,8 @@ void CompareFlowResults(TString type="ESD")
       }
       gfcCommonHistRes8 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults8thOrderGFC"));
       if(gfcCommonHistRes8)  {
-       flowValue[7] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinContent(1);//to be removed
-       flowError[7] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[7] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinContent(1);
+       flowError[7] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinError(1);
        flowValueRP[7] = (gfcCommonHistRes8->GetHistIntFlowRP())->GetBinContent(1);
        flowErrorRP[7] = (gfcCommonHistRes8->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[7] = (gfcCommonHistRes8->GetHistIntFlowPOI())->GetBinContent(1);
@@ -287,8 +286,8 @@ void CompareFlowResults(TString type="ESD")
       qcCommonHist2 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist2ndOrderQC"));
       qcCommonHistRes2 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults2ndOrderQC"));
       if(qcCommonHistRes2)  {
-       flowValue[2] = (qcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);//to be removed
-       //flowError[2] = (qcCommonHistRes2->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[2] = (qcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);
+       //flowError[2] = (qcCommonHistRes2->GetHistIntFlow())->GetBinError(1);
        flowValueRP[2] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinContent(1);
        //flowErrorRP[2] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[2] = (qcCommonHistRes2->GetHistIntFlowPOI())->GetBinContent(1);
@@ -297,8 +296,8 @@ void CompareFlowResults(TString type="ESD")
       qcCommonHist4 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist4thOrderQC"));
       qcCommonHistRes4 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults4thOrderQC"));
       if(qcCommonHistRes4) {
-       flowValue[4] = (qcCommonHistRes4->GetHistIntFlow())->GetBinContent(1);//to be removed
-       //flowError[4] = (qcCommonHistRes4->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[4] = (qcCommonHistRes4->GetHistIntFlow())->GetBinContent(1);
+       //flowError[4] = (qcCommonHistRes4->GetHistIntFlow())->GetBinError(1);
        flowValueRP[4] = (qcCommonHistRes4->GetHistIntFlowRP())->GetBinContent(1);
        //flowErrorRP[4] = (qcCommonHistRes4->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[4] = (qcCommonHistRes4->GetHistIntFlowPOI())->GetBinContent(1);
@@ -307,8 +306,8 @@ void CompareFlowResults(TString type="ESD")
       qcCommonHist6 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist6thOrderQC"));
       qcCommonHistRes6 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults6thOrderQC"));
       if(qcCommonHistRes6) {
-       flowValue[6] = (qcCommonHistRes6->GetHistIntFlow())->GetBinContent(1);//to be removed
-       //flowError[6] = (qcCommonHistRes6->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[6] = (qcCommonHistRes6->GetHistIntFlow())->GetBinContent(1);
+       //flowError[6] = (qcCommonHistRes6->GetHistIntFlow())->GetBinError(1);
        flowValueRP[6] = (qcCommonHistRes6->GetHistIntFlowRP())->GetBinContent(1);
        //flowErrorRP[6] = (qcCommonHistRes6->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[6] = (qcCommonHistRes6->GetHistIntFlowPOI())->GetBinContent(1);
@@ -317,8 +316,8 @@ void CompareFlowResults(TString type="ESD")
       qcCommonHist8 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist8thOrderQC"));
       qcCommonHistRes8 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults8thOrderQC"));
       if(qcCommonHistRes8)  {
-       flowValue[8] = (qcCommonHistRes8->GetHistIntFlow())->GetBinContent(1);//to be removed
-       //flowError[8] = (qcCommonHistRes8->GetHistIntFlow())->GetBinError(1);//to be removed
+       flowValue[8] = (qcCommonHistRes8->GetHistIntFlow())->GetBinContent(1);
+       //flowError[8] = (qcCommonHistRes8->GetHistIntFlow())->GetBinError(1);
        flowValueRP[8] = (qcCommonHistRes8->GetHistIntFlowRP())->GetBinContent(1);
        //flowErrorRP[8] = (qcCommonHistRes8->GetHistIntFlowRP())->GetBinError(1);
        flowValuePOI[8] = (qcCommonHistRes8->GetHistIntFlowPOI())->GetBinContent(1);
@@ -327,8 +326,7 @@ void CompareFlowResults(TString type="ESD")
     }
   }        
   
-  //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
-  //              !!!!  to be removed  !!!!
+  //no-name:
   Double_t dMax=flowValue[0]+flowError[0];
   Double_t dMin=flowValue[0]-flowError[0];
   
@@ -338,7 +336,6 @@ void CompareFlowResults(TString type="ESD")
       if(dMin>flowValue[i]-flowError[i]) dMin=flowValue[i]-flowError[i];
     } 
   }  
-  //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
   
   //RP:
   Double_t dMaxRP=flowValueRP[0]+flowErrorRP[0];
@@ -362,23 +359,18 @@ void CompareFlowResults(TString type="ESD")
     } 
   }  
  
-  //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
-  //                    !!!!  to be removed  !!!!
+  //no-name:
   TGraph* flowResults = new TGraphErrors(nMethods, x, flowValue, xError, flowError);
   flowResults->SetMarkerStyle(markerStyle);
   flowResults->SetMarkerColor(markerColor);
-  //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
   
   //RP:
-  TGraph* flowResultsRP = new TGraphErrors(nMethods, x, flowValueRP, xError, flowErrorRP);
-  
+  TGraph* flowResultsRP = new TGraphErrors(nMethods, x, flowValueRP, xError, flowErrorRP); 
   flowResultsRP->SetMarkerStyle(markerStyle);
   flowResultsRP->SetMarkerColor(markerColor);
   
   //POI:
   TGraph* flowResultsPOI = new TGraphErrors(nMethods, x, flowValuePOI, xError, flowErrorPOI);
-  
   flowResultsPOI->SetMarkerStyle(markerStyle);
   flowResultsPOI->SetMarkerColor(markerColor);
   
@@ -386,46 +378,38 @@ void CompareFlowResults(TString type="ESD")
   
   //----------------------------------------------------------------------------------
   //cosmetics: mesh for MC error bands (integrated flow)
-  TGraph* pMesh = NULL;//to be removed
+  TGraph* pMesh = NULL;
   TGraph* pMeshRP = NULL;
   TGraph* pMeshPOI = NULL;
   
   if(intFlowAll && mcepCommonHistRes) {
-    //Double_t valueMC = intFlowAll->GetBinContent(1);
-    //Double_t errorMC = intFlowAll->GetBinError(1);  
-    //Int_t nPts       = intFlowAll->GetNbinsX();     
-    
     Int_t nPts       = nMethods;
-    Double_t valueMC = flowValue[0];//to be removed
-    Double_t errorMC = flowError[0];//to be removed  
+    Double_t valueMC = flowValue[0];
+    Double_t errorMC = flowError[0]; 
     Double_t valueMCRP = flowValueRP[0];
     Double_t errorMCRP = flowErrorRP[0]; 
     Double_t valueMCPOI = flowValuePOI[0];
     Double_t errorMCPOI = flowErrorPOI[0]; 
     
-    pMesh = new TGraph(nPts);//to be removed
+    pMesh = new TGraph(nPts);
     pMeshRP = new TGraph(nPts);
     pMeshPOI = new TGraph(nPts);
     
-    //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
-    //             !!!! to be removed !!!!
+    //no-name:
     pMesh->SetPoint(1,0,valueMC+errorMC);
     pMesh->SetPoint(2,nPts+1,valueMC+errorMC);
     pMesh->SetPoint(3,nPts+1,valueMC-errorMC);
     pMesh->SetPoint(4,0,valueMC-errorMC);
-    pMesh->SetPoint(5,0,valueMC+errorMC);
-    
+    pMesh->SetPoint(5,0,valueMC+errorMC);    
     pMesh->SetFillStyle(meshStyle);
     pMesh->SetFillColor(meshColor);
-    //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
-    
+     
     //RP:
     pMeshRP->SetPoint(1,0,valueMCRP+errorMCRP);
     pMeshRP->SetPoint(2,nPts+1,valueMCRP+errorMCRP);
     pMeshRP->SetPoint(3,nPts+1,valueMCRP-errorMCRP);
     pMeshRP->SetPoint(4,0,valueMCRP-errorMCRP);
-    pMeshRP->SetPoint(5,0,valueMCRP+errorMCRP);
-    
+    pMeshRP->SetPoint(5,0,valueMCRP+errorMCRP);   
     pMeshRP->SetFillStyle(meshStyle);
     pMeshRP->SetFillColor(meshColor);
     
@@ -434,8 +418,7 @@ void CompareFlowResults(TString type="ESD")
     pMeshPOI->SetPoint(2,nPts+1,valueMCPOI+errorMCPOI);
     pMeshPOI->SetPoint(3,nPts+1,valueMCPOI-errorMCPOI);
     pMeshPOI->SetPoint(4,0,valueMCPOI-errorMCPOI);
-    pMeshPOI->SetPoint(5,0,valueMCPOI+errorMCPOI);
-    
+    pMeshPOI->SetPoint(5,0,valueMCPOI+errorMCPOI);    
     pMeshPOI->SetFillStyle(meshStyle);
     pMeshPOI->SetFillColor(meshColor);     
   }
@@ -448,8 +431,6 @@ void CompareFlowResults(TString type="ESD")
   TPaveText *textDefault = new TPaveText(0.05,0.77,0.95,0.90,"NDC");
   textDefault->SetTextFont(72);
   textDefault->SetTextSize(0.08);
-  //textDefault->SetLineColor(kFALSE);
-  //textDefault->SetShadowColor(kFALSE);
   
   TString *entryDefaultAvM = new TString("Average Multiplicity");
   TString *entryDefaultAnd = new TString("and"); 
@@ -459,12 +440,10 @@ void CompareFlowResults(TString type="ESD")
   textDefault->AddText(entryDefaultAnd->Data());
   textDefault->AddText(entryDefaultNumOfEvts->Data());
   
-  //results:
-  TPaveText *textResults = new TPaveText(0.05,0.12,0.95,0.70,"NDC");//to be removed
-  textResults->SetTextFont(72);//to be removed
-  textResults->SetTextSize(0.06);//to be removed
-  //textResults->SetLineColor(kFALSE);
-  //textResults->SetShadowColor(kFALSE);
+  //results (no-name):
+  TPaveText *textResults = new TPaveText(0.05,0.12,0.95,0.70,"NDC");
+  textResults->SetTextFont(72);
+  textResults->SetTextSize(0.06);
   
   //results (RP):
   TPaveText *textResultsRP = new TPaveText(0.05,0.12,0.95,0.70,"NDC");
@@ -476,8 +455,7 @@ void CompareFlowResults(TString type="ESD")
   textResultsPOI->SetTextFont(72);
   textResultsPOI->SetTextSize(0.06);
   
-  //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
-  //      !!!! to be removed !!!!!!                                      
+  //no-name:              
   TString *entryMC    = new TString("MC ...... ");
   TString *entryGFC   = new TString("GFC ..... "); 
   TString *entryQC2   = new TString("QC{2} ... ");
@@ -487,7 +465,6 @@ void CompareFlowResults(TString type="ESD")
   TString *entryFQD   = new TString("FQD ..... "); 
   TString *entryLYZ1  = new TString("LYZ ..... "); 
   TString *entryLYZEP = new TString("LYZEP ... "); 
-  //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
   
   //RP: 
   TString *entryMCRP    = new TString("MC ...... ");
@@ -511,114 +488,96 @@ void CompareFlowResults(TString type="ESD")
   TString *entryLYZ1POI  = new TString("LYZ ..... "); 
   TString *entryLYZEPPOI = new TString("LYZEP ... "); 
   
-  //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
-  //                  !!!! to be removed !!!!
+  //no-name:
   Double_t avMultMC=0.;
   Long_t nEvtsMC=0;
-  
   Double_t avMultGFC=0.;
   Long_t nEvtsGFC=0;
-  
   Double_t avMultQC2=0., avMultQC4=0., avMultQC6=0., avMultQC8=0.;
   Long_t nEvtsQC2=0, nEvtsQC4=0, nEvtsQC6=0, nEvtsQC8=0;
-  
   Double_t avMultFQD=0.;
   Long_t nEvtsFQD=0;
-  
   Double_t avMultLYZ1=0.;
   Long_t nEvtsLYZ1=0;
-  
   Double_t avMultLYZEP=0.;
   Long_t nEvtsLYZEP=0;
-  //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
   
   //RP:
   Double_t avMultMCRP=0.;
   Long_t nEvtsMCRP=0;
-  
   Double_t avMultGFCRP=0.;
   Long_t nEvtsGFCRP=0;
-  
   Double_t avMultQC2RP=0., avMultQC4RP=0., avMultQC6RP=0., avMultQC8RP=0.;
   Long_t nEvtsQC2RP=0, nEvtsQC4RP=0, nEvtsQC6RP=0, nEvtsQC8RP=0;
-  
   Double_t avMultFQDRP=0.;
   Long_t nEvtsFQDRP=0;
-  
   Double_t avMultLYZ1RP=0.;
   Long_t nEvtsLYZ1RP=0;
-  
   Double_t avMultLYZEPRP=0.;
   Long_t nEvtsLYZEPRP=0;
   
   //POI:
   Double_t avMultMCPOI=0.;
   Long_t nEvtsMCPOI=0;
-  
   Double_t avMultGFCPOI=0.;
   Long_t nEvtsGFCPOI=0;
-  
   Double_t avMultQC2POI=0., avMultQC4POI=0., avMultQC6POI=0., avMultQC8POI=0.;
   Long_t nEvtsQC2POI=0, nEvtsQC4POI=0, nEvtsQC6POI=0, nEvtsQC8POI=0;
-  
   Double_t avMultFQDPOI=0.;
   Long_t nEvtsFQDPOI=0;
-  
   Double_t avMultLYZ1POI=0.;
   Long_t nEvtsLYZ1POI=0;
-  
   Double_t avMultLYZEPPOI=0.;
   Long_t nEvtsLYZEPPOI=0;
   
   //MC:  
   if(mcepCommonHist) {
-    avMultMC = (mcepCommonHist->GetHistMultInt())->GetMean();//to be removed
-    nEvtsMC  = (mcepCommonHist->GetHistMultInt())->GetEntries();//to be removed
+    avMultMC = (mcepCommonHist->GetHistMultInt())->GetMean();
+    nEvtsMC  = (mcepCommonHist->GetHistMultInt())->GetEntries();
     avMultMCRP = (mcepCommonHist->GetHistMultInt())->GetMean();
     nEvtsMCRP  = (mcepCommonHist->GetHistMultInt())->GetEntries();
     avMultMCPOI = (mcepCommonHist->GetHistMultDiff())->GetMean();
     nEvtsMCPOI  = (mcepCommonHist->GetHistMultDiff())->GetEntries();
   }
-  //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
-  //   !!!! to be removed !!!!
+  
   if(entryMC) {   
     entryMC->Append("M = ");
     (*entryMC)+=(Long_t)avMultMC;
     entryMC->Append(", N = ");
     (*entryMC)+=(Long_t)nEvtsMC;
   }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
   if(entryMCRP) {   
     entryMCRP->Append("M = ");
     (*entryMCRP)+=(Long_t)avMultMCRP;
     entryMCRP->Append(", N = ");
     (*entryMCRP)+=(Long_t)nEvtsMCRP;
   }
- if(entryMCPOI) {   
+  
+  if(entryMCPOI) {   
    entryMCPOI->Append("M = ");
    (*entryMCPOI)+=(Long_t)avMultMCPOI;
    entryMCPOI->Append(", N = ");
    (*entryMCPOI)+=(Long_t)nEvtsMCPOI;
- }
+  }
   
  //GFC:
  if(gfcCommonHist) {
-   avMultGFC = (gfcCommonHist->GetHistMultInt())->GetMean();//to be removed
-   nEvtsGFC  = (gfcCommonHist->GetHistMultInt())->GetEntries();//to be removed
+   avMultGFC = (gfcCommonHist->GetHistMultInt())->GetMean();
+   nEvtsGFC  = (gfcCommonHist->GetHistMultInt())->GetEntries();
    avMultGFCRP = (gfcCommonHist->GetHistMultInt())->GetMean();
    nEvtsGFCRP  = (gfcCommonHist->GetHistMultInt())->GetEntries();
    avMultGFCPOI = (gfcCommonHist->GetHistMultDiff())->GetMean();
    nEvtsGFCPOI  = (gfcCommonHist->GetHistMultDiff())->GetEntries();
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- //   !!!! to be removed !!!! 
  if(entryGFC) { 
    entryGFC->Append("M = ");
    (*entryGFC)+=(Long_t)avMultGFC;
    entryGFC->Append(", N = ");
    (*entryGFC)+=(Long_t)nEvtsGFC;
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
  if(entryGFCRP) { 
    entryGFCRP->Append("M = ");
    (*entryGFCRP)+=(Long_t)avMultGFCRP;
@@ -634,15 +593,14 @@ void CompareFlowResults(TString type="ESD")
  
  //QC:
  if(qcCommonHist2) {
-   avMultQC2 = (qcCommonHist2->GetHistMultInt())->GetMean();//to be removed
-   nEvtsQC2  = (qcCommonHist2->GetHistMultInt())->GetEntries();//to be removed
+   avMultQC2 = (qcCommonHist2->GetHistMultInt())->GetMean();
+   nEvtsQC2  = (qcCommonHist2->GetHistMultInt())->GetEntries();
    avMultQC2RP = (qcCommonHist2->GetHistMultInt())->GetMean();
    nEvtsQC2RP  = (qcCommonHist2->GetHistMultInt())->GetEntries();
    avMultQC2POI = (qcCommonHist2->GetHistMultDiff())->GetMean();
    nEvtsQC2POI  = (qcCommonHist2->GetHistMultDiff())->GetEntries();
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- //   !!!! to be removed !!!!   
  if(entryQC2)
  { 
   entryQC2->Append("M = ");
@@ -650,7 +608,7 @@ void CompareFlowResults(TString type="ESD")
   entryQC2->Append(", N = ");
   (*entryQC2)+=(Long_t)nEvtsQC2;
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
  if(entryQC2RP)
  { 
   entryQC2RP->Append("M = ");
@@ -658,6 +616,7 @@ void CompareFlowResults(TString type="ESD")
   entryQC2RP->Append(", N = ");
   (*entryQC2RP)+=(Long_t)nEvtsQC2RP;
  } 
  if(entryQC2POI)
  { 
   entryQC2POI->Append("M = ");
@@ -668,15 +627,14 @@ void CompareFlowResults(TString type="ESD")
 
  if(qcCommonHist4)
  {
-  avMultQC4 = (qcCommonHist4->GetHistMultInt())->GetMean();//to be removed
-  nEvtsQC4  = (qcCommonHist4->GetHistMultInt())->GetEntries();//to be removed
+  avMultQC4 = (qcCommonHist4->GetHistMultInt())->GetMean();
+  nEvtsQC4  = (qcCommonHist4->GetHistMultInt())->GetEntries();
   avMultQC4RP = (qcCommonHist4->GetHistMultInt())->GetMean();
   nEvtsQC4RP  = (qcCommonHist4->GetHistMultInt())->GetEntries();
   avMultQC4POI = (qcCommonHist4->GetHistMultDiff())->GetMean();
   nEvtsQC4POI  = (qcCommonHist4->GetHistMultDiff())->GetEntries();
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- //   !!!! to be removed !!!!  
  if(entryQC4)
  {
   entryQC4->Append("M = ");
@@ -684,7 +642,7 @@ void CompareFlowResults(TString type="ESD")
   entryQC4->Append(", N = ");
   (*entryQC4)+=(Long_t)nEvtsQC4;
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
  if(entryQC4RP)
  {
   entryQC4RP->Append("M = ");
@@ -692,6 +650,7 @@ void CompareFlowResults(TString type="ESD")
   entryQC4RP->Append(", N = ");
   (*entryQC4RP)+=(Long_t)nEvtsQC4RP;
  }
  if(entryQC4POI)
  {
   entryQC4POI->Append("M = ");
@@ -702,15 +661,14 @@ void CompareFlowResults(TString type="ESD")
    
  if(qcCommonHist6)
  {
-  avMultQC6 = (qcCommonHist6->GetHistMultInt())->GetMean();//to be removed
-  nEvtsQC6  = (qcCommonHist6->GetHistMultInt())->GetEntries();//to be removed
+  avMultQC6 = (qcCommonHist6->GetHistMultInt())->GetMean();
+  nEvtsQC6  = (qcCommonHist6->GetHistMultInt())->GetEntries();
   avMultQC6RP = (qcCommonHist6->GetHistMultInt())->GetMean();
   nEvtsQC6RP  = (qcCommonHist6->GetHistMultInt())->GetEntries();
   avMultQC6POI = (qcCommonHist6->GetHistMultDiff())->GetMean();
   nEvtsQC6POI  = (qcCommonHist6->GetHistMultDiff())->GetEntries();
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- //   !!!! to be removed !!!!   
  if(entryQC6)
  {  
   entryQC6->Append("M = ");
@@ -718,7 +676,7 @@ void CompareFlowResults(TString type="ESD")
   entryQC6->Append(", N = ");
   (*entryQC6)+=(Long_t)nEvtsQC6;
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
  if(entryQC6RP)
  {  
   entryQC6RP->Append("M = ");
@@ -726,6 +684,7 @@ void CompareFlowResults(TString type="ESD")
   entryQC6RP->Append(", N = ");
   (*entryQC6RP)+=(Long_t)nEvtsQC6RP;
  }
  if(entryQC6POI)
  {  
   entryQC6POI->Append("M = ");
@@ -736,15 +695,14 @@ void CompareFlowResults(TString type="ESD")
    
  if(qcCommonHist8)
  {
-  avMultQC8 = (qcCommonHist8->GetHistMultInt())->GetMean();//to be removed
-  nEvtsQC8  = (qcCommonHist8->GetHistMultInt())->GetEntries();//to be removed
+  avMultQC8 = (qcCommonHist8->GetHistMultInt())->GetMean();
+  nEvtsQC8  = (qcCommonHist8->GetHistMultInt())->GetEntries();
   avMultQC8RP = (qcCommonHist8->GetHistMultInt())->GetMean();
   nEvtsQC8RP  = (qcCommonHist8->GetHistMultInt())->GetEntries();
   avMultQC8POI = (qcCommonHist8->GetHistMultDiff())->GetMean();
   nEvtsQC8POI  = (qcCommonHist8->GetHistMultDiff())->GetEntries();    
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- //   !!!! to be removed !!!!    
+  
  if(entryQC8)
  {
   entryQC8->Append("M = ");
@@ -752,7 +710,7 @@ void CompareFlowResults(TString type="ESD")
   entryQC8->Append(", N = ");
   (*entryQC8)+=(Long_t)nEvtsQC8;
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
  if(entryQC8RP)
  {
   entryQC8RP->Append("M = ");
@@ -760,6 +718,7 @@ void CompareFlowResults(TString type="ESD")
   entryQC8RP->Append(", N = ");
   (*entryQC8RP)+=(Long_t)nEvtsQC8RP;
  }
  if(entryQC8POI)
  {
   entryQC8POI->Append("M = ");
@@ -771,15 +730,14 @@ void CompareFlowResults(TString type="ESD")
  //FQD:
  if(fqdCommonHist)
  {
-  avMultFQD = (fqdCommonHist->GetHistMultInt())->GetMean();//to be removed
-  nEvtsFQD  = (fqdCommonHist->GetHistMultInt())->GetEntries();//to be removed
+  avMultFQD = (fqdCommonHist->GetHistMultInt())->GetMean();
+  nEvtsFQD  = (fqdCommonHist->GetHistMultInt())->GetEntries();
   avMultFQDRP = (fqdCommonHist->GetHistMultInt())->GetMean();
   nEvtsFQDRP  = (fqdCommonHist->GetHistMultInt())->GetEntries();
   avMultFQDPOI = (fqdCommonHist->GetHistMultDiff())->GetMean();
   nEvtsFQDPOI  = (fqdCommonHist->GetHistMultDiff())->GetEntries();
  } 
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- //   !!!! to be removed !!!!  
  if(entryFQD)
  {
   entryFQD->Append("M = ");
@@ -787,7 +745,7 @@ void CompareFlowResults(TString type="ESD")
   entryFQD->Append(", N = ");
   (*entryFQD)+=(Long_t)nEvtsFQD;
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
  if(entryFQDRP)
  {
   entryFQDRP->Append("M = ");
@@ -795,6 +753,7 @@ void CompareFlowResults(TString type="ESD")
   entryFQDRP->Append(", N = ");
   (*entryFQDRP)+=(Long_t)nEvtsFQDRP;
  }
  if(entryFQDPOI)
  {
   entryFQDPOI->Append("M = ");
@@ -806,15 +765,14 @@ void CompareFlowResults(TString type="ESD")
  //LYZ1:
  if(lyz1CommonHist)
  {
-  avMultLYZ1 = (lyz1CommonHist->GetHistMultInt())->GetMean();//to be removed
-  nEvtsLYZ1  = (lyz1CommonHist->GetHistMultInt())->GetEntries();//to be removed
+  avMultLYZ1 = (lyz1CommonHist->GetHistMultInt())->GetMean();
+  nEvtsLYZ1  = (lyz1CommonHist->GetHistMultInt())->GetEntries();
   avMultLYZ1RP = (lyz1CommonHist->GetHistMultInt())->GetMean();
   nEvtsLYZ1RP  = (lyz1CommonHist->GetHistMultInt())->GetEntries();
   avMultLYZ1POI = (lyz1CommonHist->GetHistMultDiff())->GetMean();
   nEvtsLYZ1POI  = (lyz1CommonHist->GetHistMultDiff())->GetEntries();
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- //   !!!! to be removed !!!!   
  if(entryLYZ1) 
  {
   entryLYZ1->Append("M = ");
@@ -822,7 +780,7 @@ void CompareFlowResults(TString type="ESD")
   entryLYZ1->Append(", N = ");
   (*entryLYZ1)+=(Long_t)nEvtsLYZ1;
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
  if(entryLYZ1RP) 
  {
   entryLYZ1RP->Append("M = ");
@@ -830,6 +788,7 @@ void CompareFlowResults(TString type="ESD")
   entryLYZ1RP->Append(", N = ");
   (*entryLYZ1RP)+=(Long_t)nEvtsLYZ1RP;
  }
  if(entryLYZ1POI) 
  {
   entryLYZ1POI->Append("M = ");
@@ -841,15 +800,14 @@ void CompareFlowResults(TString type="ESD")
  //LYZEP:
  if(lyzepCommonHist)
  {
-  avMultLYZEP = (lyzepCommonHist->GetHistMultInt())->GetMean();//to be removed
-  nEvtsLYZEP  = (lyzepCommonHist->GetHistMultInt())->GetEntries();//to be removed
+  avMultLYZEP = (lyzepCommonHist->GetHistMultInt())->GetMean();
+  nEvtsLYZEP  = (lyzepCommonHist->GetHistMultInt())->GetEntries();
   avMultLYZEPRP = (lyzepCommonHist->GetHistMultInt())->GetMean();
   nEvtsLYZEPRP  = (lyzepCommonHist->GetHistMultInt())->GetEntries();
   avMultLYZEPPOI = (lyzepCommonHist->GetHistMultDiff())->GetMean();
   nEvtsLYZEPPOI  = (lyzepCommonHist->GetHistMultDiff())->GetEntries();    
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- //   !!!! to be removed !!!!   
  if(entryLYZEP)
  {
   entryLYZEP->Append("M = ");
@@ -857,7 +815,7 @@ void CompareFlowResults(TString type="ESD")
   entryLYZEP->Append(", N = ");
   (*entryLYZEP)+=(Long_t)nEvtsLYZEP;
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
  if(entryLYZEPRP)
  {
   entryLYZEPRP->Append("M = ");
@@ -865,6 +823,7 @@ void CompareFlowResults(TString type="ESD")
   entryLYZEPRP->Append(", N = ");
   (*entryLYZEPRP)+=(Long_t)nEvtsLYZEPRP;
  }
  if(entryLYZEPPOI)
  {
   entryLYZEPPOI->Append("M = ");
@@ -873,8 +832,7 @@ void CompareFlowResults(TString type="ESD")
   (*entryLYZEPPOI)+=(Long_t)nEvtsLYZEPPOI;
  }
 
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- //   !!!! to be removed !!!!   
+ //no-name:
  if(textResults)
  {
   textResults->AddText(entryMC->Data());
@@ -887,7 +845,7 @@ void CompareFlowResults(TString type="ESD")
   textResults->AddText(entryLYZ1->Data());
   textResults->AddText(entryLYZEP->Data());
  }
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
  //RP:
  if(textResultsRP)
  {
@@ -916,10 +874,11 @@ void CompareFlowResults(TString type="ESD")
  }
  //----------------------------------------------------------------------------------
  
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- //                           !!!! to be removed !!!!
  //----------------------------------------------------------------------------------
- //final drawing for integrated flow:
+ //final drawing for integrated flow (no-name):
  TCanvas* intFlowAllCanvas = new TCanvas("Integrated Flow","Integrated Flow",1000,600);
  
  intFlowAllCanvas->Divide(2,1);
@@ -957,9 +916,7 @@ void CompareFlowResults(TString type="ESD")
   textDefault->Draw();
   textResults->Draw();
  }
  //----------------------------------------------------------------------------------
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
 
  
  
@@ -1008,6 +965,9 @@ void CompareFlowResults(TString type="ESD")
  }//end of if(plotIntFlowRP} 
  //----------------------------------------------------------------------------------
  
  //----------------------------------------------------------------------------------
  //final drawing for integrated flow of POI (i.e. of particles of interest):
  if(plotIntFlowPOI)
@@ -1118,38 +1078,6 @@ void CompareFlowResults(TString type="ESD")
  //----------------------------------------------------------------------------------
  
  
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- //               !!!! to be removed !!!!
- //----------------------------------------------------------------------------------
- //cosmetics: Monte Carlo error bands for differential flow
- TGraph* pMeshDiffFlow = NULL;
- if(mcepCommonHistRes)
- {
-  Int_t nPtsDiffFlow = (mcepCommonHistRes->GetHistDiffFlow())->GetNbinsX();
-  Double_t binWidth = (mcepCommonHistRes->GetHistDiffFlow())->GetBinWidth(1);//assuming that all bins have the same width
-       
-  pMeshDiffFlow = new TGraph(2*nPtsDiffFlow+1);
-  
-  Double_t valueMC=0., errorMC=0.;
-  for(Int_t i=1;i<nPtsDiffFlow+1;i++)
-  {
-   valueMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinContent(i);
-   errorMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinError(i);       
-   pMeshDiffFlow->SetPoint(i,(i-0.5)*binWidth,valueMC+errorMC);
-  }    
-  for(Int_t i=nPtsDiffFlow+1;i<2*nPtsDiffFlow+1;i++)
-  {
-   valueMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinContent(2*nPtsDiffFlow+1-i);
-   errorMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinError(2*nPtsDiffFlow+1-i);       
-   pMeshDiffFlow->SetPoint(i,(2*nPtsDiffFlow-i+0.5)*binWidth,valueMC-errorMC); 
-  }
-  pMeshDiffFlow->SetPoint(2*nPtsDiffFlow+1,0.5*binWidth,valueMC+errorMC); 
-  pMeshDiffFlow->SetFillStyle(meshStyle);
-  pMeshDiffFlow->SetFillColor(meshColor);
- } 
- //----------------------------------------------------------------------------------
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
 
  //----------------------------------------------------------------------------------
  //RP:
@@ -1381,8 +1309,6 @@ void CompareFlowResults(TString type="ESD")
  //----------------------------------------------------------------------------------
    
  //MCEP = Monte Carlo Event Plane
- Double_t avMultDiffFlowMC=0.;//to be removed
- Double_t nEvtsDiffFlowMC=0;//to be removed
  Double_t avMultDiffFlowMCRP=0.;
  Double_t nEvtsDiffFlowMCRP=0;
  Double_t avMultDiffFlowMCPOI=0.;
@@ -1391,8 +1317,6 @@ void CompareFlowResults(TString type="ESD")
  {
   if(mcepCommonHistRes)
   {
-   (mcepCommonHistRes->GetHistDiffFlow())->SetMarkerColor(markerColorMC);//to be removed
-   (mcepCommonHistRes->GetHistDiffFlow())->SetMarkerStyle(markerStyleMC);//to be removed
    (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorMC);
    (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleMC);
    (mcepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorMC);
@@ -1404,8 +1328,6 @@ void CompareFlowResults(TString type="ESD")
   } 
   if(mcepCommonHist)
   {
-   avMultDiffFlowMC = (mcepCommonHist->GetHistMultDiff())->GetMean();//to be removed
-   nEvtsDiffFlowMC  = (mcepCommonHist->GetHistMultDiff())->GetEntries();//to be removed
    avMultDiffFlowMCRP = (mcepCommonHist->GetHistMultInt())->GetMean();
    nEvtsDiffFlowMCRP  = (mcepCommonHist->GetHistMultInt())->GetEntries();
    avMultDiffFlowMCPOI = (mcepCommonHist->GetHistMultDiff())->GetMean();
@@ -1424,8 +1346,6 @@ void CompareFlowResults(TString type="ESD")
  {
   if(gfcCommonHistRes2)
   {
-   (gfcCommonHistRes2->GetHistDiffFlow())->SetMarkerColor(markerColorGFC2);//to be removed
-   (gfcCommonHistRes2->GetHistDiffFlow())->SetMarkerStyle(markerStyleGFC2);//to be removed
    (gfcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC2);
    (gfcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC2);
    (gfcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC2);
@@ -1437,8 +1357,6 @@ void CompareFlowResults(TString type="ESD")
   }
   if(gfcCommonHistRes4)
   { 
-   (gfcCommonHistRes4->GetHistDiffFlow())->SetMarkerColor(markerColorGFC4);//to be removed
-   (gfcCommonHistRes4->GetHistDiffFlow())->SetMarkerStyle(markerStyleGFC4);//to be removed
    (gfcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC4);
    (gfcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC4);
    (gfcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC4);
@@ -1450,8 +1368,6 @@ void CompareFlowResults(TString type="ESD")
   }
   if(gfcCommonHistRes6)
   { 
-   (gfcCommonHistRes6->GetHistDiffFlow())->SetMarkerColor(markerColorGFC6);//to be removed
-   (gfcCommonHistRes6->GetHistDiffFlow())->SetMarkerStyle(markerStyleGFC6);//to be removed
    (gfcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC6);
    (gfcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC6);
    (gfcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC6);
@@ -1463,8 +1379,6 @@ void CompareFlowResults(TString type="ESD")
   }
   if(gfcCommonHistRes8)
   { 
-   (gfcCommonHistRes8->GetHistDiffFlow())->SetMarkerColor(markerColorGFC8);//to be removed
-   (gfcCommonHistRes8->GetHistDiffFlow())->SetMarkerStyle(markerStyleGFC8);//to be removed
    (gfcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC8);
    (gfcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC8);
    (gfcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC8);
@@ -1476,8 +1390,6 @@ void CompareFlowResults(TString type="ESD")
   }
   if(gfcCommonHist)
   {
-   avMultDiffFlowGFC = (gfcCommonHist->GetHistMultDiff())->GetMean();//to be removed
-   nEvtsDiffFlowGFC  = (gfcCommonHist->GetHistMultDiff())->GetEntries();//to be removed
    avMultDiffFlowGFCRP = (gfcCommonHist->GetHistMultInt())->GetMean();   
    nEvtsDiffFlowGFCRP  = (gfcCommonHist->GetHistMultInt())->GetEntries();
    avMultDiffFlowGFCPOI = (gfcCommonHist->GetHistMultDiff())->GetMean();
@@ -1486,10 +1398,6 @@ void CompareFlowResults(TString type="ESD")
  }
   
  //QC = Q-cumulants
- Double_t avMultDiffFlowQC2=0., avMultDiffFlowQC4=0.;//to be removed
- Double_t nEvtsDiffFlowQC2=0., nEvtsDiffFlowQC4=0.;//to be removed
- Double_t avMultDiffFlowQC6=0., avMultDiffFlowQC8=0.;//to be removed
- Double_t nEvtsDiffFlowQC6=0., nEvtsDiffFlowQC8=0.;//to be removed
  Double_t avMultDiffFlowQC2RP=0.;
  Double_t nEvtsDiffFlowQC2RP=0.;
  Double_t avMultDiffFlowQC2POI=0.;
@@ -1512,8 +1420,6 @@ void CompareFlowResults(TString type="ESD")
   //QC{2}
   if(qcCommonHistRes2)
   {
-   (qcCommonHistRes2->GetHistDiffFlow())->SetMarkerColor(markerColorQC2);//to be removed
-   (qcCommonHistRes2->GetHistDiffFlow())->SetMarkerStyle(markerStyleQC2);//to be removed
    (qcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC2);
    (qcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC2);
    (qcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC2);
@@ -1525,8 +1431,6 @@ void CompareFlowResults(TString type="ESD")
   }
   if(qcCommonHist2)
   {
-   avMultDiffFlowQC2 = (qcCommonHist2->GetHistMultDiff())->GetMean();//to be removed
-   nEvtsDiffFlowQC2  = (qcCommonHist2->GetHistMultDiff())->GetEntries();//to be removed
    avMultDiffFlowQC2RP = (qcCommonHist2->GetHistMultInt())->GetMean();
    nEvtsDiffFlowQC2RP  = (qcCommonHist2->GetHistMultInt())->GetEntries();
    avMultDiffFlowQC2POI = (qcCommonHist2->GetHistMultDiff())->GetMean();
@@ -1535,8 +1439,6 @@ void CompareFlowResults(TString type="ESD")
   //QC{4}
   if(qcCommonHistRes4)
   {
-   (qcCommonHistRes4->GetHistDiffFlow())->SetMarkerColor(markerColorQC4);//to be removed
-   (qcCommonHistRes4->GetHistDiffFlow())->SetMarkerStyle(markerStyleQC4);//to be removed
    (qcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC4);
    (qcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC4);
    (qcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC4);
@@ -1548,8 +1450,6 @@ void CompareFlowResults(TString type="ESD")
   }
   if(qcCommonHist4)
   {
-   avMultDiffFlowQC4 = (qcCommonHist4->GetHistMultDiff())->GetMean();//to be removed
-   nEvtsDiffFlowQC4  = (qcCommonHist4->GetHistMultDiff())->GetEntries();//to be removed
    avMultDiffFlowQC4RP = (qcCommonHist4->GetHistMultInt())->GetMean();
    nEvtsDiffFlowQC4RP  = (qcCommonHist4->GetHistMultInt())->GetEntries();
    avMultDiffFlowQC4POI = (qcCommonHist4->GetHistMultDiff())->GetMean();
@@ -1558,8 +1458,6 @@ void CompareFlowResults(TString type="ESD")
   //QC{6}
   if(qcCommonHistRes6)
   {
-   (qcCommonHistRes6->GetHistDiffFlow())->SetMarkerColor(markerColorQC6);//to be removed
-   (qcCommonHistRes6->GetHistDiffFlow())->SetMarkerStyle(markerStyleQC6);//to be removed
    (qcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC6);
    (qcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC6);
    (qcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC6);
@@ -1571,8 +1469,6 @@ void CompareFlowResults(TString type="ESD")
   }
   if(qcCommonHist6)
   {
-   avMultDiffFlowQC6 = (qcCommonHist6->GetHistMultDiff())->GetMean();//to be removed
-   nEvtsDiffFlowQC6  = (qcCommonHist6->GetHistMultDiff())->GetEntries();//to be removed
    avMultDiffFlowQC6RP = (qcCommonHist6->GetHistMultInt())->GetMean();
    nEvtsDiffFlowQC6RP  = (qcCommonHist6->GetHistMultInt())->GetEntries();
    avMultDiffFlowQC6POI = (qcCommonHist6->GetHistMultDiff())->GetMean();
@@ -1581,8 +1477,6 @@ void CompareFlowResults(TString type="ESD")
   //QC{8}
   if(qcCommonHistRes8)
   {
-   (qcCommonHistRes8->GetHistDiffFlow())->SetMarkerColor(markerColorQC8);//to be removed
-   (qcCommonHistRes8->GetHistDiffFlow())->SetMarkerStyle(markerStyleQC8);//to be removed
    (qcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC8);
    (qcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC8);
    (qcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC8);
@@ -1594,8 +1488,6 @@ void CompareFlowResults(TString type="ESD")
   }
   if(qcCommonHist8)
   {
-   avMultDiffFlowQC8 = (qcCommonHist8->GetHistMultDiff())->GetMean();//to be removed
-   nEvtsDiffFlowQC8  = (qcCommonHist8->GetHistMultDiff())->GetEntries();//to be removed
    avMultDiffFlowQC8RP = (qcCommonHist8->GetHistMultInt())->GetMean();
    nEvtsDiffFlowQC8RP  = (qcCommonHist8->GetHistMultInt())->GetEntries();
    avMultDiffFlowQC8POI = (qcCommonHist8->GetHistMultDiff())->GetMean();
@@ -1604,8 +1496,6 @@ void CompareFlowResults(TString type="ESD")
  } 
 
  //LYZ2 = Lee-Yang Zeros (2nd run)
- Double_t avMultDiffFlowLYZ2=0.;//to be removed
- Double_t nEvtsDiffFlowLYZ2=0;//to be removed
  Double_t avMultDiffFlowLYZ2RP=0.;
  Double_t nEvtsDiffFlowLYZ2RP=0;
  Double_t avMultDiffFlowLYZ2POI=0.;
@@ -1614,8 +1504,6 @@ void CompareFlowResults(TString type="ESD")
  {
   if(lyz2CommonHistRes)
   {
-   (lyz2CommonHistRes->GetHistDiffFlow())->SetMarkerColor(markerColorLYZ2);//to be removed
-   (lyz2CommonHistRes->GetHistDiffFlow())->SetMarkerStyle(markerStyleLYZ2);//to be removed
    (lyz2CommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorLYZ2);
    (lyz2CommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleLYZ2);
    (lyz2CommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorLYZ2);
@@ -1627,8 +1515,6 @@ void CompareFlowResults(TString type="ESD")
   } 
   if(lyz2CommonHist)
   {
-   avMultDiffFlowLYZ2 = (lyz2CommonHist->GetHistMultDiff())->GetMean();//to be removed
-   nEvtsDiffFlowLYZ2  = (lyz2CommonHist->GetHistMultDiff())->GetEntries();//to be removed
    avMultDiffFlowLYZ2RP = (lyz2CommonHist->GetHistMultInt())->GetMean();
    nEvtsDiffFlowLYZ2RP  = (lyz2CommonHist->GetHistMultInt())->GetEntries();
    avMultDiffFlowLYZ2POI = (lyz2CommonHist->GetHistMultDiff())->GetMean();
@@ -1637,8 +1523,6 @@ void CompareFlowResults(TString type="ESD")
  } 
 
  //LYZEP = Lee-Yang Zeros Event Plane
- Double_t avMultDiffFlowLYZEP=0.;//to be removed
- Double_t nEvtsDiffFlowLYZEP=0;//to be removed
  Double_t avMultDiffFlowLYZEPRP=0.;
  Double_t nEvtsDiffFlowLYZEPRP=0;
  Double_t avMultDiffFlowLYZEPPOI=0.;
@@ -1647,9 +1531,6 @@ void CompareFlowResults(TString type="ESD")
  {
   if(lyzepCommonHistRes)
   {
-   (lyzepCommonHistRes->GetHistDiffFlow())->Scale(0.01);//to be improved
-   (lyzepCommonHistRes->GetHistDiffFlow())->SetMarkerColor(markerColorLYZEP);//to be removed
-   (lyzepCommonHistRes->GetHistDiffFlow())->SetMarkerStyle(markerStyleLYZEP);//to be removed
    (lyzepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorLYZEP);
    (lyzepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleLYZEP);
    (lyzepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorLYZEP);
@@ -1661,8 +1542,6 @@ void CompareFlowResults(TString type="ESD")
   } 
   if(lyzepCommonHist)
   {
-   avMultDiffFlowLYZEP = (lyzepCommonHist->GetHistMultDiff())->GetMean();//to be removed
-   nEvtsDiffFlowLYZEP  = (lyzepCommonHist->GetHistMultDiff())->GetEntries();//to be removed
    avMultDiffFlowLYZEPRP = (lyzepCommonHist->GetHistMultInt())->GetMean();
    nEvtsDiffFlowLYZEPRP  = (lyzepCommonHist->GetHistMultInt())->GetEntries();
    avMultDiffFlowLYZEPPOI = (lyzepCommonHist->GetHistMultDiff())->GetMean();
@@ -1670,180 +1549,7 @@ void CompareFlowResults(TString type="ESD")
   } 
  } 
 
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
- // !!!! to be removed !!!!
- //----------------------------------------------------------------------------------
- //final drawing for differential flow:
- TCanvas* diffFlowAllCanvas = new TCanvas("Differential Flow","Differential Flow",1000,600);
- diffFlowAllCanvas->Divide(2,1);
- //1st pad is for plot:
- (diffFlowAllCanvas->cd(1))->SetPad(0.0,0.0,0.75,1.0);
-  
- if(styleHistPt)
- {
-  styleHistPt->Draw();
- }
-
- if(pMeshDiffFlow)
- {
-  pMeshDiffFlow->Draw("LFSAME");
- }
- //MC 
- if(mcepCommonHistRes)
- { 
-  (mcepCommonHistRes->GetHistDiffFlow())->Draw("E1PSAME");
- }
- //GFC
- if(gfcCommonHistRes2)
- { 
-  (gfcCommonHistRes2->GetHistDiffFlow())->Draw("E1PSAME"); 
- } 
- if(gfcCommonHistRes4)
- { 
-  (gfcCommonHistRes4->GetHistDiffFlow())->Draw("E1PSAME"); 
- } 
- if(gfcCommonHistRes6)
- { 
-  (gfcCommonHistRes6->GetHistDiffFlow())->Draw("E1PSAME"); 
- } 
- if(gfcCommonHistRes8)
- { 
-  (gfcCommonHistRes8->GetHistDiffFlow())->Draw("E1PSAME"); 
- }   
- //QC
- if(qcCommonHistRes2)
- { 
-  (qcCommonHistRes2->GetHistDiffFlow())->Draw("E1PSAME");
- }
- if(qcCommonHistRes4)
- { 
-  (qcCommonHistRes4->GetHistDiffFlow())->Draw("E1PSAME");
- }
- //LYZ2
- if(lyz2CommonHistRes)
- { 
-  (lyz2CommonHistRes->GetHistDiffFlow())->Draw("E1PSAME");
- }
- //LYZEP
- if(lyzepCommonHistRes)
- { 
-  (lyzepCommonHistRes->GetHistDiffFlow())->Draw("E1PSAME");
- }
- //2nd pad is for legend:
- (diffFlowAllCanvas->cd(2))->SetPad(0.75,0.0,1.0,1.0);
-
- TLegend* legendDiffFlow = new TLegend(0.02,0.25,0.97,0.75);
- legendDiffFlow->SetTextFont(72);
- legendDiffFlow->SetTextSize(0.06);
-
- //legend's entries:
- TString *entryDiffMC   = new TString("MC ....... ");
- TString *entryDiffGFC2 = new TString("GFC{2} ... ");
- TString *entryDiffGFC4 = new TString("GFC{4} ... ");
- TString *entryDiffGFC6 = new TString("GFC{6} ... ");
- TString *entryDiffGFC8 = new TString("GFC{8} ... "); 
- TString *entryDiffQC2  = new TString("QC{2} .... ");
- TString *entryDiffQC4  = new TString("QC{4} .... ");
- TString *entryDiffLYZ2 = new TString("LYZ ...... ");
- TString *entryDiffLYZEP = new TString("LYZEP ... ");
- //MC
- if(mcepCommonHistRes)
- {
-  (mcepCommonHistRes->GetHistDiffFlow())->SetFillStyle(meshStyle);
-  (mcepCommonHistRes->GetHistDiffFlow())->SetFillColor(meshColor);
-  entryDiffMC->Append("M = ");
-  (*entryDiffMC)+=(Long_t)avMultDiffFlowMC;
-  entryDiffMC->Append(", N = ");
-  (*entryDiffMC)+=(Long_t)nEvtsDiffFlowMC; 
-  legendDiffFlow->AddEntry(mcepCommonHistRes->GetHistDiffFlow(),entryDiffMC->Data(),"f");
- }
-
- //GFC
- if(gfcCommonHistRes2)
- {
-  entryDiffGFC2->Append("M = ");
-  (*entryDiffGFC2)+=(Long_t)avMultDiffFlowGFC;
-  entryDiffGFC2->Append(", N = ");
-  (*entryDiffGFC2)+=(Long_t)nEvtsDiffFlowGFC; 
-  legendDiffFlow->AddEntry(gfcCommonHistRes2->GetHistDiffFlow(),entryDiffGFC2->Data(),"p");
- }
- if(gfcCommonHistRes4)
- {
-  entryDiffGFC4->Append("M = ");
-  (*entryDiffGFC4)+=(Long_t)avMultDiffFlowGFC;
-  entryDiffGFC4->Append(", N = ");
-  (*entryDiffGFC4)+=(Long_t)nEvtsDiffFlowGFC; 
-  legendDiffFlow->AddEntry(gfcCommonHistRes4->GetHistDiffFlow(),entryDiffGFC4->Data(),"p");
- }
- if(gfcCommonHistRes6)
- {
-  entryDiffGFC6->Append("M = ");
-  (*entryDiffGFC6)+=(Long_t)avMultDiffFlowGFC;
-  entryDiffGFC6->Append(", N = ");
-  (*entryDiffGFC6)+=(Long_t)nEvtsDiffFlowGFC; 
-  legendDiffFlow->AddEntry(gfcCommonHistRes6->GetHistDiffFlow(),entryDiffGFC6->Data(),"p");
- } 
- if(gfcCommonHistRes8)
- {
-  entryDiffGFC8->Append("M = ");
-  (*entryDiffGFC8)+=(Long_t)avMultDiffFlowGFC;
-  entryDiffGFC8->Append(", N = ");
-  (*entryDiffGFC8)+=(Long_t)nEvtsDiffFlowGFC; 
-  legendDiffFlow->AddEntry(gfcCommonHistRes8->GetHistDiffFlow(),entryDiffGFC8->Data(),"p");
- }  
- //QC
- if(qcCommonHistRes2)
- {
-  entryDiffQC2->Append("M = ");
-  (*entryDiffQC2)+=(Long_t)avMultDiffFlowQC2;
-  entryDiffQC2->Append(", N = ");
-  (*entryDiffQC2)+=(Long_t)nEvtsDiffFlowQC2; 
-  legendDiffFlow->AddEntry(qcCommonHistRes2->GetHistDiffFlow(),entryDiffQC2->Data(),"p");
- }
- if(qcCommonHistRes4)
- {
-  entryDiffQC4->Append("M = ");
-  (*entryDiffQC4)+=(Long_t)avMultDiffFlowQC4;
-  entryDiffQC4->Append(", N = ");
-  (*entryDiffQC4)+=(Long_t)nEvtsDiffFlowQC4; 
-  legendDiffFlow->AddEntry(qcCommonHistRes4->GetHistDiffFlow(),entryDiffQC4->Data(),"p");
- }
- //LYZ
- if(lyz2CommonHistRes)
- {
-  entryDiffLYZ2->Append("M = ");
-  (*entryDiffLYZ2)+=(Long_t)avMultDiffFlowLYZ2;
-  entryDiffLYZ2->Append(", N = ");
-  (*entryDiffLYZ2)+=(Long_t)nEvtsDiffFlowLYZ2; 
-  legendDiffFlow->AddEntry(lyz2CommonHistRes->GetHistDiffFlow(),entryDiffLYZ2->Data(),"p");
- }
- //LYZEP
- if(lyzepCommonHistRes)
- {
-  entryDiffLYZEP->Append("M = ");
-  (*entryDiffLYZEP)+=(Long_t)avMultDiffFlowLYZEP;
-  entryDiffLYZEP->Append(", N = ");
-  (*entryDiffLYZEP)+=(Long_t)nEvtsDiffFlowLYZEP; 
-  legendDiffFlow->AddEntry(lyzepCommonHistRes->GetHistDiffFlow(),entryDiffLYZEP->Data(),"p");
- }
-
-
- //drawing finally the legend in the 2nd pad:     
- if(legendDiffFlow)
- {
-  legendDiffFlow->SetMargin(0.15);
-  legendDiffFlow->Draw();
- }
- //----------------------------------------------------------------------------------
- //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
 
  //----------------------------------------------------------------------------------
  //final drawing for differential flow (Pt, RP):
  //set here the results of which methods will be plotted by default:
@@ -2304,7 +2010,7 @@ void CompareFlowResults(TString type="ESD")
  
   //MC 
   if(plotMCPtPOI && mcepCommonHistRes)
-  { 
+  { kTRUE;
    (mcepCommonHistRes->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
   }
   //GFC
@@ -2334,7 +2040,7 @@ void CompareFlowResults(TString type="ESD")
    (qcCommonHistRes4->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
   }
   if(plotQC6PtPOI && qcCommonHistRes6)
-  { 
+  { kTRUE;
    //(qcCommonHistRes6->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
   }
   if(plotQC8PtPOI && qcCommonHistRes8)