+ fCorrelEventplaneMCDATA = new TH2F("fCorrelEventplaneMCDATA","fCorrelEventplaneMCDATA",40, -2.*TMath::Pi(), 2.*TMath::Pi(), 40, -2.*TMath::Pi(), 2.*TMath::Pi());
+ fCorrelEventplaneMCDATA->GetXaxis()->SetTitle("#phi (event plane)");
+ fCorrelEventplaneMCDATA->GetYaxis()->SetTitle("#phi (MC event plane)");
+ fCorrelEventplaneMCDATA->Sumw2();
+
+ // cross check for event plane resolution
+ fEPDistCent = new TH2F("fEPDistCent","fEPDistCent",20, -1.*TMath::Pi(), TMath::Pi(), fCentralityNbins-1, fBinsCentrality);
+ fEPDistCent->GetXaxis()->SetTitle("#phi (#Psi_{EP})");
+ fEPDistCent->GetYaxis()->SetTitle("Centrality");
+ fEPDistCent->Sumw2();
+
+ fPhiCent = new TH2F("fPhiCent","fPhiCent",200, -2.*TMath::Pi(), 2.*TMath::Pi(), fCentralityNbins-1, fBinsCentrality);
+ fPhiCent->GetXaxis()->SetTitle("#phi");
+ fPhiCent->GetYaxis()->SetTitle("Centrality");
+ fPhiCent->Sumw2();
+
+ fPcosEPCent = new TProfile("fPcosEPCent","fPcosEPCent", 100,0,100);
+ fPcosEPCent->GetXaxis()->SetTitle("Centrality");
+ fPcosEPCent->GetYaxis()->SetTitle("#LT cos 2 #Psi_{EP} #GT");
+ fPcosEPCent->Sumw2();
+
+ fPsinEPCent = new TProfile("fPsinEPCent","fPsinEPCent", 100,0,100);
+ fPsinEPCent->GetXaxis()->SetTitle("Centrality");
+ fPsinEPCent->GetYaxis()->SetTitle("#LT sin 2 #Psi_{EP} #GT");
+ fPsinEPCent->Sumw2();
+
+ fPcosPhiCent = new TProfile("fPcosPhiCent","fPcosPhiCent", 100,0,100);
+ fPcosPhiCent->GetXaxis()->SetTitle("Centrality");
+ fPcosPhiCent->GetYaxis()->SetTitle("#LT cos 2 #phi #GT");
+ fPcosPhiCent->Sumw2();
+
+ fPsinPhiCent = new TProfile("fPsinPhiCent","fPsinPhiCent", 100,0,100);
+ fPsinPhiCent->GetXaxis()->SetTitle("Centrality");
+ fPsinPhiCent->GetYaxis()->SetTitle("#LT sin 2 #phi #GT");
+ fPsinPhiCent->Sumw2();
+