+void PlotRefFlowVsM()
+{
+ // Make plots for reference flow vs multiplicity:
+
+ TCanvas *cRefFlowVsM = new TCanvas("cRefFlowVsM","Reference Flow");
+ cRefFlowVsM->Divide(2,2);
+
+ TLegend *lRefFlowVsM = new TLegend(0.1,0.7,0.33,0.9);
+ lRefFlowVsM->SetFillStyle(0);
+ //lRefFlowVsM->SetHeader(" minClustersTpcRP");
+
+ TString refFlowVsMFlag[4] = {"v_{2}{2,QC}","v_{2}{4,QC}","v_{2}{6,QC}","v_{2}{8,QC}"};
+
+ for(Int_t co=0;co<4;co++) // cumulant order
+ {
+ cRefFlowVsM->cd(co+1);
+ // Style histrogram:.q
+ TH1D *styleHist = (TH1D*)StyleHist(refFlowVsMFlag[co].Data(),co)->Clone();
+ if(styleHist)
+ {
+ styleHist->GetXaxis()->SetRangeUser(refFlowVsMxRange[0],refFlowVsMxRange[1]);
+ styleHist->GetYaxis()->SetRangeUser(refFlowVsMyRange[0],refFlowVsMyRange[1]);
+ styleHist->Draw();
+ }
+ // Plot first the meshes for reference flow out of their own pad:
+ TGraph *rfMeshOutOf = GetErrorMesh(refFlowVsM[0][0][refFlowVsMeshOrder]); // to be improved - hardwired 0, now it works only for the "current" file and for QC
+ if(rfMeshOutOf)
+ {
+ rfMeshOutOf->SetFillColor(refFlowVsMMeshColor);
+ rfMeshOutOf->Draw("lfsame");
+ if(refFlowVsMeshOrder==co){lRefFlowVsM->AddEntry(rfMeshOutOf,Form("v_{2}{%i,QC} stat. error",2*(refFlowVsMeshOrder+1)),"f");}
+ }
+ // Plot on top of the the meshes for reference flow with markers in its own pad:
+ if(refFlowVsM[0][0][co])
+ {
+ refFlowVsM[0][0][co]->SetMarkerStyle(refFlowVsMMarkerStyle[co]);
+ refFlowVsM[0][0][co]->SetMarkerColor(refFlowVsMMarkerColor[co]);
+ refFlowVsM[0][0][co]->Draw("e1same");
+ lRefFlowVsM->AddEntry(refFlowVsM[0][0][co],Form("v_{2}{%i,QC}",2*(co+1)),"p");
+ }
+ } // end of for(Int_t co=0;co<4;co++) // cumulant order
+
+ // Draw a common legend in a 1st pad:
+ cRefFlowVsM->cd(1);
+ lRefFlowVsM->Draw("same");
+
+} // end of void PlotRefFlowVsM()
+
+// =====================================================================================
+
+void PlotCumulantsVsM()