Browse Source

Updated macros Draw_AcceptanceComparison.C

Peter 2 years ago
parent
commit
585d01c368
1 changed files with 146 additions and 99 deletions
  1. 146 99
      macro/Draw_AcceptanceComparison.C

+ 146 - 99
macro/Draw_AcceptanceComparison.C

@@ -2,6 +2,12 @@
 
 std::pair<double, double> GetFHCalCoord(int iModule);
 
+const std::vector<double> pt_binning_fhcal = {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1., 1.2, 1.4, 1.6, 1.8, 2., 2.5};
+const std::vector<double> pt_binning_tpc = {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1., 1.2, 1.4, 1.6, 1.8, 2., 2.5};
+
+const int N_pt_bins_fhcal = pt_binning_fhcal.size() - 1;
+const int N_pt_bins_tpc = pt_binning_tpc.size() - 1;
+
 void Draw_AcceptanceComparison()
 {
   gROOT->SetStyle("Pub");
@@ -10,24 +16,36 @@ void Draw_AcceptanceComparison()
 
   TFile *fiQv[2], *fiRes[2], *fiFlow[2];
 
-  fiQv[0] = new TFile("/home/peter/Documents/Dataset/good_acceptance/qv_11.5gev.root", "read");
-  fiQv[1] = new TFile("/home/peter/Documents/Dataset/bad_acceptance/qv_11.5gev.root", "read");
+  // fiQv[0] = new TFile("/home/peter/Documents/Dataset/good_acceptance/qv_11.5gev.root", "read");
+  // fiQv[1] = new TFile("/home/peter/Documents/Dataset/bad_acceptance/qv_11.5gev.root", "read");
+
+  // fiRes[0] = new TFile("/home/peter/Documents/Dataset/good_acceptance/res_11.5gev_fit.root", "read");
+  // fiRes[1] = new TFile("/home/peter/Documents/Dataset/bad_acceptance/res_11.5gev_fit.root", "read");
+
+  // fiFlow[0] = new TFile("/home/peter/Documents/Dataset/good_acceptance/flow_11.5gev.root", "read");
+  // fiFlow[1] = new TFile("/home/peter/Documents/Dataset/bad_acceptance/flow_11.5gev.root", "read");
 
-  fiRes[0] = new TFile("/home/peter/Documents/Dataset/good_acceptance/res_11.5gev_fit.root", "read");
-  fiRes[1] = new TFile("/home/peter/Documents/Dataset/bad_acceptance/res_11.5gev_fit.root", "read");
+  // TFile *fiRes_uncorrected = new TFile("/home/peter/Documents/Dataset/bad_acceptance/res_11.5gev_fit_uncorrected.root", "read");
+  // TFile *fiFlow_incorrected = new TFile("/home/peter/Documents/Dataset/bad_acceptance/flow_11.5gev_uncorrected.root", "read");
 
-  fiFlow[0] = new TFile("/home/peter/Documents/Dataset/good_acceptance/flow_11.5gev.root", "read");
-  fiFlow[1] = new TFile("/home/peter/Documents/Dataset/bad_acceptance/flow_11.5gev.root", "read");
+  fiQv[0] = new TFile("/home/peter/Documents/Dataset/mpd_data/7.7gev/good_acceptance/qv.root", "read");
+  fiQv[1] = new TFile("/home/peter/Documents/Dataset/mpd_data/7.7gev/bad_acceptance/qv.root", "read");
 
-  TFile *fiRes_uncorrected = new TFile("/home/peter/Documents/Dataset/bad_acceptance/res_11.5gev_fit_uncorrected.root", "read");
-  TFile *fiFlow_incorrected = new TFile("/home/peter/Documents/Dataset/bad_acceptance/flow_11.5gev_uncorrected.root", "read");
+  fiRes[0] = new TFile("/home/peter/Documents/Dataset/mpd_data/7.7gev/good_acceptance/res_fit.root", "read");
+  fiRes[1] = new TFile("/home/peter/Documents/Dataset/mpd_data/7.7gev/bad_acceptance/res_fit.root", "read");
 
-  const int npid = 3;  // pion, kaon, proton
+  fiFlow[0] = new TFile("/home/peter/Documents/Dataset/mpd_data/7.7gev/good_acceptance/flow.root", "read");
+  fiFlow[1] = new TFile("/home/peter/Documents/Dataset/mpd_data/7.7gev/bad_acceptance/flow.root", "read");
+
+  TFile *fiRes_uncorrected = new TFile("/home/peter/Documents/Dataset/mpd_data/7.7gev/bad_acceptance_no_corr/res_fit.root", "read");
+  TFile *fiFlow_incorrected = new TFile("/home/peter/Documents/Dataset/mpd_data/7.7gev/bad_acceptance_no_corr/flow.root", "read");
+
+  const int npid = 6;  // h+, pion, kaon, proton, h-, h+-
   const int nharm = 2; // v1, v2
   const int nsub = 2;  //TPC eta-gap0 (0), FHCal total (8)
 
-  const int pidId[npid] = {1, 2, 3};
-  const int subId[nsub] = {0, 8};
+  const int pidId[npid] = {0, 1, 2, 3, 4, 999};
+  const int subId[nsub] = {0, 14};
   const TString qvDirNames[3] = {"Raw", "Recentering", "Flattening"};
 
   TH2F *h2TpcAcceptance[2], *h2FHCalAcceptance[2], *hQxCent[2][3][nsub], *hQyCent[2][3][nsub], *hPsiEpCent[2][3][nsub];
@@ -82,7 +100,7 @@ void Draw_AcceptanceComparison()
   {
     for (int j = 0; j < nharm; j++)
     {
-      for (int k = 0; k < npid; k++)
+      for (int k = 0; k < npid - 1; k++)
       {
         p2FlowPt[i][j][k] = (TProfile2D *)fiFlow[i]->Get(Form("FHCalEp/pv%iPtFHCal_pid%i", j + 1, pidId[k]));
         p2FlowPt[i][j][k]->SetName(Form("fi%i_%s", i, p2FlowPt[i][j][k]->GetName()));
@@ -105,6 +123,30 @@ void Draw_AcceptanceComparison()
     }
   }
 
+  for (int i = 0; i < 2; i++)
+  {
+    for (int j = 0; j < nharm; j++)
+    {
+      p2FlowPt[i][j][npid - 1] = (TProfile2D *)p2FlowPt[i][j][0]->Clone(Form("fi%i_pv%iPtFHCal_pid%i", i, j + 1, pidId[npid - 1]));
+      p2FlowPt[i][j][npid - 1]->Add(p2FlowPt[i][j][4]);
+      if (i == 0)
+      {
+        p2FlowPt_uncorrected[j][npid - 1] = (TProfile2D *)p2FlowPt_uncorrected[j][0]->Clone(Form("fi%i_uncorrected_pv%iPtFHCal_pid%i", i, j + 1, pidId[npid - 1]));
+        p2FlowPt_uncorrected[j][npid - 1]->Add(p2FlowPt_uncorrected[j][4]);
+      }
+      if (j == 1)
+      {
+        p2FlowTpcPt[i][npid - 1] = (TProfile2D *)p2FlowTpcPt[i][0]->Clone(Form("fi%i_pv%iPtTpc_gap%i_pid%i", i, j + 1, subId[0], pidId[npid - 1]));
+        p2FlowTpcPt[i][npid - 1]->Add(p2FlowTpcPt[i][4]);
+        if (i == 0)
+        {
+          p2FlowTpcPt_uncorrected[npid - 1] = (TProfile2D *)p2FlowTpcPt_uncorrected[0]->Clone(Form("fi%i_uncorrected_pv%iPtFHCal_gap%i_pid%i", i, j + 1, subId[0], pidId[npid - 1]));
+          p2FlowTpcPt_uncorrected[npid - 1]->Add(p2FlowTpcPt_uncorrected[4]);
+        }
+      }
+    }
+  }
+
   const int Ncent = 8;
   const int NcentFlow = 3;
   const std::pair<int, int> centFlowBins[NcentFlow] = {{1, 2}, {2, 4}, {4, 8}};
@@ -165,22 +207,22 @@ void Draw_AcceptanceComparison()
         for (int k = 0; k < npid; k++)
         {
           pFlowPt[cent][i][j][k] = p2FlowPt[i][j][k]->ProfileX(Form("cent%i_%s", cent, p2FlowPt[i][j][k]->GetName()), centFlowBins[cent].first, centFlowBins[cent].second);
-          pFlowPt[cent][i][j][k]->Rebin(10);
+          pFlowPt[cent][i][j][k] = (TProfile *)pFlowPt[cent][i][j][k]->Rebin(N_pt_bins_fhcal, Form("rebin_%s", pFlowPt[cent][i][j][k]->GetName()), &pt_binning_fhcal[0]);
 
           if (i == 0)
           {
             pFlowPt_uncorrected[cent][j][k] = p2FlowPt_uncorrected[j][k]->ProfileX(Form("cent%i_%s", cent, p2FlowPt_uncorrected[j][k]->GetName()), centFlowBins[cent].first, centFlowBins[cent].second);
-            pFlowPt_uncorrected[cent][j][k]->Rebin(10);
+            pFlowPt_uncorrected[cent][j][k] = (TProfile *)pFlowPt_uncorrected[cent][j][k]->Rebin(N_pt_bins_fhcal, Form("rebin_%s", pFlowPt_uncorrected[cent][j][k]->GetName()), &pt_binning_fhcal[0]);
           }
 
           if (j == 1)
           {
             pFlowTpcPt[cent][i][k] = p2FlowTpcPt[i][k]->ProfileX(Form("cent%i_%s", cent, p2FlowTpcPt[i][k]->GetName()), centFlowBins[cent].first, centFlowBins[cent].second);
-            pFlowTpcPt[cent][i][k]->Rebin(15);
+            pFlowTpcPt[cent][i][k] = (TProfile *)pFlowTpcPt[cent][i][k]->Rebin(N_pt_bins_tpc, Form("rebin_%s", pFlowTpcPt[cent][i][k]->GetName()), &pt_binning_tpc[0]);
             if (i == 0)
             {
               pFlowTpcPt_uncorrected[cent][k] = p2FlowTpcPt_uncorrected[k]->ProfileX(Form("cent%i_%s", cent, p2FlowTpcPt_uncorrected[k]->GetName()), centFlowBins[cent].first, centFlowBins[cent].second);
-              pFlowTpcPt_uncorrected[cent][k]->Rebin(15);
+              pFlowTpcPt_uncorrected[cent][k] = (TProfile *)pFlowTpcPt_uncorrected[cent][k]->Rebin(N_pt_bins_tpc, Form("rebin_%s", pFlowTpcPt_uncorrected[cent][k]->GetName()), &pt_binning_tpc[0]);
             }
           }
         }
@@ -195,25 +237,30 @@ void Draw_AcceptanceComparison()
   // Acceptance
 
   TCanvas *canv_TpcAcceptance = new TCanvas("canv_TpcAcceptance", "canv_TpcAcceptance", 900, 700);
-  TLegend *leg_TpcAcceptance = new TLegend(0.2, 0.2, 0.5, 0.35);
+  TLegend *leg_TpcAcceptance = new TLegend(0.2, 0.2, 0.52, 0.45);
+  TLatex tex1;
   canv_TpcAcceptance->cd();
 
   hTpcAcceptance[0]->SetLineColor(kBlack);
   hTpcAcceptance[1]->SetLineColor(kRed);
 
   hTpcAcceptance[0]->GetYaxis()->SetTitle("dN/d#varphi");
-  hTpcAcceptance[0]->GetXaxis()->SetTitle("#varphi, rad");
+  hTpcAcceptance[0]->GetXaxis()->SetTitle("#varphi (rad)");
 
-  hTpcAcceptance[0]->GetYaxis()->SetRangeUser(0., 1e6);
+  hTpcAcceptance[0]->GetYaxis()->SetRangeUser(0., 5e6);
 
   hTpcAcceptance[0]->Draw();
   hTpcAcceptance[1]->Draw("same");
 
-  leg_TpcAcceptance->AddEntry(hTpcAcceptance[0], "Good Acceptance", "l");
-  leg_TpcAcceptance->AddEntry(hTpcAcceptance[1], "Bad Acceptance", "l");
+  leg_TpcAcceptance->SetBorderSize(0);
+  leg_TpcAcceptance->AddEntry(hTpcAcceptance[0], "Uniform acc.", "l");
+  leg_TpcAcceptance->AddEntry(hTpcAcceptance[1], "Non-uniform acc.", "l");
   leg_TpcAcceptance->Draw();
 
-  canv_TpcAcceptance->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/FHCalAcceptancePhi.png"));
+  tex1.DrawLatexNDC(0.5,0.5, Form("Au+Au at #sqrt{s_{NN}}=7.7 GeV, UrQMD, GEANT4, 0-80%"));
+
+  canv_TpcAcceptance->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/TpcAcceptancePhi.png"));
+  
 
   TCanvas *canv_FHCalAcceptance = new TCanvas("canv_FHCalAcceptance", "canv_FHCalAcceptance", 900, 700);
   TLegend *leg_FHCalAcceptance = new TLegend(0.2, 0.7, 0.5, 0.85);
@@ -231,15 +278,15 @@ void Draw_AcceptanceComparison()
   hFHCalAcceptance[0]->Draw("hist");
   hFHCalAcceptance[1]->Draw("hist same");
 
-  leg_FHCalAcceptance->AddEntry(hFHCalAcceptance[0], "Good Acceptance", "l");
-  leg_FHCalAcceptance->AddEntry(hFHCalAcceptance[1], "Bad Acceptance", "l");
+  leg_FHCalAcceptance->AddEntry(hFHCalAcceptance[0], "Uniform acc.", "l");
+  leg_FHCalAcceptance->AddEntry(hFHCalAcceptance[1], "Non-uniform acc.", "l");
   leg_FHCalAcceptance->Draw();
 
   canv_FHCalAcceptance->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/FHCalAcceptancePhi.png"));
 
   const TString AcceptanceType[2] = {"Good", "Bad"};
 
-  TCanvas *canv_FHCalVisual[2][2]; // good-bad acceptance, L-R
+  TCanvas *canv_FHCalVisual[2][2]; // good-Non-uniform acc., L-R
   const TString SideName[2] = {"L", "R"};
   for (int i = 0; i < 2; i++)
   {
@@ -408,9 +455,9 @@ void Draw_AcceptanceComparison()
   hResTpc[1]->Draw("p hist same");
   // hResTpc_uncorrected->Draw("p hist same");
 
-  leg_resTpc_v2->AddEntry(hResTpc[0], "Good Acceptance", "p");
-  leg_resTpc_v2->AddEntry(hResTpc[1], "Bad Acceptance", "p");
-  // leg_resTpc_v2->AddEntry(hResTpc_uncorrected, "Bad Acceptance (uncorr.)", "p");
+  leg_resTpc_v2->AddEntry(hResTpc[0], "Uniform acc.", "p");
+  leg_resTpc_v2->AddEntry(hResTpc[1], "Non-uniform acc.", "p");
+  // leg_resTpc_v2->AddEntry(hResTpc_uncorrected, "Non-uniform acc. (uncorr.)", "p");
   leg_resTpc_v2->Draw();
 
   canv_resTpc_v2->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/res_v2_TPC.png"));
@@ -440,9 +487,9 @@ void Draw_AcceptanceComparison()
   hResTpc[1]->Draw("p hist same");
   hResTpc_uncorrected->Draw("p hist same");
 
-  leg_resTpc_v2_uncorr->AddEntry(hResTpc[0], "Good Acceptance", "p");
-  leg_resTpc_v2_uncorr->AddEntry(hResTpc[1], "Bad Acceptance", "p");
-  leg_resTpc_v2_uncorr->AddEntry(hResTpc_uncorrected, "Bad Acceptance (uncorr.)", "p");
+  leg_resTpc_v2_uncorr->AddEntry(hResTpc[0], "Uniform acc.", "p");
+  leg_resTpc_v2_uncorr->AddEntry(hResTpc[1], "Non-uniform acc.", "p");
+  leg_resTpc_v2_uncorr->AddEntry(hResTpc_uncorrected, "Non-uniform acc. (uncorr.)", "p");
   leg_resTpc_v2_uncorr->Draw();
 
   canv_resTpc_v2_uncorr->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/res_v2_TPC_uncorrected.png"));
@@ -472,8 +519,8 @@ void Draw_AcceptanceComparison()
   hResv1FHCal[0]->Draw("p hist");
   hResv1FHCal[1]->Draw("p hist same");
 
-  leg_resFHCal_v1->AddEntry(hResv1FHCal[0], "Good Acceptance", "p");
-  leg_resFHCal_v1->AddEntry(hResv1FHCal[1], "Bad Acceptance", "p");
+  leg_resFHCal_v1->AddEntry(hResv1FHCal[0], "Uniform acc.", "p");
+  leg_resFHCal_v1->AddEntry(hResv1FHCal[1], "Non-uniform acc.", "p");
   leg_resFHCal_v1->Draw();
 
   canv_resFHCal_v1->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/res_v1_FHCal.png"));
@@ -503,9 +550,9 @@ void Draw_AcceptanceComparison()
   hResv1FHCal[1]->Draw("p hist same");
   hResv1FHCal_uncorrected->Draw("p hist same");
 
-  leg_resFHCal_v1_uncorr->AddEntry(hResv1FHCal[0], "Good Acceptance", "p");
-  leg_resFHCal_v1_uncorr->AddEntry(hResv1FHCal[1], "Bad Acceptance", "p");
-  leg_resFHCal_v1_uncorr->AddEntry(hResv1FHCal_uncorrected, "Bad Acceptance (uncorr.)", "p");
+  leg_resFHCal_v1_uncorr->AddEntry(hResv1FHCal[0], "Uniform acc.", "p");
+  leg_resFHCal_v1_uncorr->AddEntry(hResv1FHCal[1], "Non-uniform acc.", "p");
+  leg_resFHCal_v1_uncorr->AddEntry(hResv1FHCal_uncorrected, "Non-uniform acc. (uncorr.)", "p");
   leg_resFHCal_v1_uncorr->Draw();
 
   canv_resFHCal_v1_uncorr->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/res_v1_FHCal_uncorrected.png"));
@@ -535,8 +582,8 @@ void Draw_AcceptanceComparison()
   hResv2FHCal[0]->Draw("p hist");
   hResv2FHCal[1]->Draw("p hist same");
 
-  leg_resFHCal_v2->AddEntry(hResv2FHCal[0], "Good Acceptance", "p");
-  leg_resFHCal_v2->AddEntry(hResv2FHCal[1], "Bad Acceptance", "p");
+  leg_resFHCal_v2->AddEntry(hResv2FHCal[0], "Uniform acc.", "p");
+  leg_resFHCal_v2->AddEntry(hResv2FHCal[1], "Non-uniform acc.", "p");
   leg_resFHCal_v2->Draw();
 
   canv_resFHCal_v2->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/res_v2_FHCal.png"));
@@ -566,9 +613,9 @@ void Draw_AcceptanceComparison()
   hResv2FHCal[1]->Draw("p hist same");
   hResv2FHCal_uncorrected->Draw("p hist same");
 
-  leg_resFHCal_v2_uncorr->AddEntry(hResv2FHCal[0], "Good Acceptance", "p");
-  leg_resFHCal_v2_uncorr->AddEntry(hResv2FHCal[1], "Bad Acceptance", "p");
-  leg_resFHCal_v2_uncorr->AddEntry(hResv2FHCal_uncorrected, "Bad Acceptance (uncorr.)", "p");
+  leg_resFHCal_v2_uncorr->AddEntry(hResv2FHCal[0], "Uniform acc.", "p");
+  leg_resFHCal_v2_uncorr->AddEntry(hResv2FHCal[1], "Non-uniform acc.", "p");
+  leg_resFHCal_v2_uncorr->AddEntry(hResv2FHCal_uncorrected, "Non-uniform acc. (uncorr.)", "p");
   leg_resFHCal_v2_uncorr->Draw();
 
   canv_resFHCal_v2_uncorr->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/res_v2_FHCal_uncorrected.png"));
@@ -576,24 +623,24 @@ void Draw_AcceptanceComparison()
   // FLOW
   // pFlowPt[NcentFlow][2][nharm][npid]
 
-  const double x_pt_min[nharm][NcentFlow][npid] = {{{0.2, 0.2, 0.4}, {0.2, 0.2, 0.4}, {0.2, 0.2, 0.4}},
-                                                   {{0.2, 0.2, 0.4}, {0.2, 0.2, 0.4}, {0.2, 0.2, 0.4}}};
-  const double x_pt_max = 2.;
-  const double y_pt_min[nharm][NcentFlow][npid] = {{{-0.029, -0.029, 0.006}, {-0.045, -0.029, 0.006}, {-0.065, -0.075, -0.021}},
-                                                   {{0., 0., 0.}, {0., 0., 0.}, {0., 0., 0.}}};
-  const double y_pt_max[nharm][NcentFlow][npid] = {{{0.019, 0.039, 0.044}, {0., 0.049, 0.044}, {0.0029, 0.049, 0.029}},
-                                                   {{0.079, 0.079, 0.059}, {0.075, 0.095, 0.085}, {0.139, 0.159, 0.099}}};
+  const double x_pt_min[nharm][NcentFlow][npid] = {{{0.2, 0.2, 0.2, 0.4, 0.2, 0.2}, {0.2, 0.2, 0.2, 0.4, 0.2, 0.2}, {0.2, 0.2, 0.2, 0.4, 0.2, 0.2}},
+                                                   {{0.2, 0.2, 0.2, 0.4, 0.2, 0.2}, {0.2, 0.2, 0.2, 0.4, 0.2, 0.2}, {0.2, 0.2, 0.2, 0.4, 0.2, 0.2}}};
+  const double x_pt_max = 2.5;
+  const double y_pt_min[nharm][NcentFlow][npid] = {{{-0.045, -0.045, -0.029, 0.006, -0.045, -0.045}, {-0.065, -0.065, -0.029, 0.006, -0.065, -0.065}, {-0.079, -0.079, -0.025, -0.0021, -0.079, -0.079}},
+                                                   {{-0.003, -0.003, -0.003, -0.003, -0.003, -0.003}, {-0.003, -0.003, -0.003, -0.003, -0.003, -0.003}, {-0.003, -0.003, -0.003, -0.003, -0.003, -0.003}}};
+  const double y_pt_max[nharm][NcentFlow][npid] = {{{0.019, 0.019, 0.039, 0.044, 0.019, 0.019}, {0., 0., 0.049, 0.044, 0., 0.}, {0.0029, 0.0029, 0.049, 0.029, 0.0029, 0.0029}},
+                                                   {{0.079, 0.079, 0.079, 0.059, 0.079, 0.079}, {0.091, 0.091, 0.095, 0.095, 0.091, 0.098}, {0.139, 0.139, 0.159, 0.099, 0.139, 0.139}}};
 
   const double x_y_min = -1.2;
   const double x_y_max = 1.2;
 
-  const double y_y_min[nharm][NcentFlow][npid] = {{{-0.039, -0.039, -0.054}, {-0.069, -0.069, -0.064}, {-0.079, -0.079, -0.034}},
-                                                  {{-0.01, -0.01, -0.01}, {0.02, 0.02, 0.01}, {0.01, 0.01, 0.002}}};
-  const double y_y_max[nharm][NcentFlow][npid] = {{{0.039, 0.039, 0.054}, {0.069, 0.069, 0.064}, {0.079, 0.079, 0.034}},
-                                                  {{0.039, 0.039, 0.039}, {0.045, 0.045, 0.055}, {0.049, 0.049, 0.069}}};
+  const double y_y_min[nharm][NcentFlow][npid] = {{{-0.039, -0.039, -0.039, -0.054, -0.039, -0.039}, {-0.069, -0.069, -0.069, -0.064, -0.069, -0.069}, {-0.079, -0.079, -0.079, -0.034, -0.079, -0.079}},
+                                                  {{-0.01, -0.01, -0.01, -0.01, -0.01, -0.01}, {0.02, 0.02, 0.01, 0.02, 0.02, 0.02}, {0.01, 0.01, 0.002, 0.01, 0.01, 0.01}}};
+  const double y_y_max[nharm][NcentFlow][npid] = {{{0.039, 0.039, 0.054, 0.039, 0.039, 0.039}, {0.069, 0.069, 0.064, 0.069, 0.069, 0.069}, {0.079, 0.079, 0.034, 0.079, 0.079, 0.079}},
+                                                  {{0.039, 0.039, 0.039, 0.039, 0.039, 0.039}, {0.045, 0.045, 0.055, 0.045, 0.045, 0.045}, {0.049, 0.049, 0.069, 0.049, 0.049, 0.049}}};
 
   const TString centFlowName[] = {"0-10%", "10-40%", "40-80%"};
-  const TString pidName[] = {"#pi^{+}", "K^{+}", "p"};
+  const TString pidName[] = {"h^{+}", "#pi^{+}", "K^{+}", "p", "h^{-}", "h#pm"};
 
   TCanvas *canv_flow_FHCal_Pt[NcentFlow][nharm][npid];
   TLegend *leg_flow_FHCal_Pt[NcentFlow][nharm][npid];
@@ -614,22 +661,22 @@ void Draw_AcceptanceComparison()
         leg_flow_FHCal_Pt[cent][i][j] = new TLegend(0.15, 0.72, 0.42, 0.89);
         canv_flow_FHCal_Pt[cent][i][j]->cd();
 
-        pFlowPt[cent][0][i][j]->GetXaxis()->SetTitle("p_{T}, GeV/c");
-        pFlowPt[cent][0][i][j]->GetYaxis()->SetTitle(Form("v_{%i}(FHCal EP)", i + 1));
+        pFlowPt[cent][0][i][j]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
+        pFlowPt[cent][0][i][j]->GetYaxis()->SetTitle(Form("v_{%i}{#Psi_{1,FHCal}}", i + 1));
 
         pFlowPt[cent][0][i][j]->GetXaxis()->SetRangeUser(x_pt_min[i][cent][j], x_pt_max);
         pFlowPt[cent][0][i][j]->GetYaxis()->SetRangeUser(y_pt_min[i][cent][j], 1.2 * y_pt_max[i][cent][j]);
 
-        pFlowPt[cent][0][i][j]->SetMarkerStyle(20);
-        pFlowPt[cent][1][i][j]->SetMarkerStyle(25);
+        pFlowPt[cent][0][i][j]->SetMarkerStyle(25);
+        pFlowPt[cent][1][i][j]->SetMarkerStyle(20);
 
         pFlowPt[cent][0][i][j]->SetMarkerSize(1.6);
         pFlowPt[cent][1][i][j]->SetMarkerSize(1.6);
 
-        pFlowPt[cent][0][i][j]->SetMarkerColor(kBlack);
-        pFlowPt[cent][1][i][j]->SetMarkerColor(kRed);
-        pFlowPt[cent][0][i][j]->SetLineColor(kBlack);
-        pFlowPt[cent][1][i][j]->SetLineColor(kRed);
+        pFlowPt[cent][0][i][j]->SetMarkerColor(kRed);
+        pFlowPt[cent][1][i][j]->SetMarkerColor(kBlack);
+        pFlowPt[cent][0][i][j]->SetLineColor(kRed);
+        pFlowPt[cent][1][i][j]->SetLineColor(kBlack);
 
         pFlowPt_uncorrected[cent][i][j]->SetMarkerStyle(26);
         pFlowPt_uncorrected[cent][i][j]->SetMarkerSize(1.6);
@@ -641,9 +688,9 @@ void Draw_AcceptanceComparison()
         // pFlowPt_uncorrected[cent][i][j]->Draw("p same");
 
         leg_flow_FHCal_Pt[cent][i][j]->SetHeader(Form("%s, %s", centFlowName[cent].Data(), pidName[j].Data()));
-        leg_flow_FHCal_Pt[cent][i][j]->AddEntry(pFlowPt[cent][0][i][j], "Good Acceptance", "p");
-        leg_flow_FHCal_Pt[cent][i][j]->AddEntry(pFlowPt[cent][1][i][j], "Bad Acceptance", "p");
-        // leg_flow_FHCal_Pt[cent][i][j]->AddEntry(pFlowPt_uncorrected[cent][i][j], "Bad Acceptance (uncorr.)", "p");
+        leg_flow_FHCal_Pt[cent][i][j]->AddEntry(pFlowPt[cent][0][i][j], "Uniform acc.", "p");
+        leg_flow_FHCal_Pt[cent][i][j]->AddEntry(pFlowPt[cent][1][i][j], "Non-uniform acc.", "p");
+        // leg_flow_FHCal_Pt[cent][i][j]->AddEntry(pFlowPt_uncorrected[cent][i][j], "Non-uniform acc. (uncorr.)", "p");
         leg_flow_FHCal_Pt[cent][i][j]->Draw();
 
         canv_flow_FHCal_Pt[cent][i][j]->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/v%i_pt_FHCalEp_cent%i_pid%i.png", i + 1, cent, j));
@@ -652,22 +699,22 @@ void Draw_AcceptanceComparison()
         leg_flow_FHCal_Pt_uncorr[cent][i][j] = new TLegend(0.15, 0.72, 0.42, 0.89);
         canv_flow_FHCal_Pt_uncorr[cent][i][j]->cd();
 
-        pFlowPt[cent][0][i][j]->GetXaxis()->SetTitle("p_{T}, GeV/c");
-        pFlowPt[cent][0][i][j]->GetYaxis()->SetTitle(Form("v_{%i}(FHCal EP)", i + 1));
+        pFlowPt[cent][0][i][j]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
+        pFlowPt[cent][0][i][j]->GetYaxis()->SetTitle(Form("v_{%i}{#Psi_{1,FHCal}}", i + 1));
 
         pFlowPt[cent][0][i][j]->GetXaxis()->SetRangeUser(x_pt_min[i][cent][j], x_pt_max);
         pFlowPt[cent][0][i][j]->GetYaxis()->SetRangeUser(y_pt_min[i][cent][j], 1.2 * y_pt_max[i][cent][j]);
 
-        pFlowPt[cent][0][i][j]->SetMarkerStyle(20);
-        pFlowPt[cent][1][i][j]->SetMarkerStyle(25);
+        pFlowPt[cent][0][i][j]->SetMarkerStyle(25);
+        pFlowPt[cent][1][i][j]->SetMarkerStyle(20);
 
         pFlowPt[cent][0][i][j]->SetMarkerSize(1.6);
         pFlowPt[cent][1][i][j]->SetMarkerSize(1.6);
 
-        pFlowPt[cent][0][i][j]->SetMarkerColor(kBlack);
-        pFlowPt[cent][1][i][j]->SetMarkerColor(kRed);
-        pFlowPt[cent][0][i][j]->SetLineColor(kBlack);
-        pFlowPt[cent][1][i][j]->SetLineColor(kRed);
+        pFlowPt[cent][0][i][j]->SetMarkerColor(kRed);
+        pFlowPt[cent][1][i][j]->SetMarkerColor(kBlack);
+        pFlowPt[cent][0][i][j]->SetLineColor(kRed);
+        pFlowPt[cent][1][i][j]->SetLineColor(kBlack);
 
         pFlowPt_uncorrected[cent][i][j]->SetMarkerStyle(26);
         pFlowPt_uncorrected[cent][i][j]->SetMarkerSize(1.6);
@@ -675,13 +722,13 @@ void Draw_AcceptanceComparison()
         pFlowPt_uncorrected[cent][i][j]->SetMarkerColor(kBlue);
 
         pFlowPt[cent][0][i][j]->Draw("p");
-        pFlowPt[cent][1][i][j]->Draw("p same");
+        // pFlowPt[cent][1][i][j]->Draw("p same");
         pFlowPt_uncorrected[cent][i][j]->Draw("p same");
 
         leg_flow_FHCal_Pt_uncorr[cent][i][j]->SetHeader(Form("%s, %s", centFlowName[cent].Data(), pidName[j].Data()));
-        leg_flow_FHCal_Pt_uncorr[cent][i][j]->AddEntry(pFlowPt[cent][0][i][j], "Good Acceptance", "p");
-        leg_flow_FHCal_Pt_uncorr[cent][i][j]->AddEntry(pFlowPt[cent][1][i][j], "Bad Acceptance", "p");
-        leg_flow_FHCal_Pt_uncorr[cent][i][j]->AddEntry(pFlowPt_uncorrected[cent][i][j], "Bad Acceptance (uncorr.)", "p");
+        leg_flow_FHCal_Pt_uncorr[cent][i][j]->AddEntry(pFlowPt[cent][0][i][j], "Uniform acc.", "p");
+        // leg_flow_FHCal_Pt_uncorr[cent][i][j]->AddEntry(pFlowPt[cent][1][i][j], "Non-uniform acc.", "p");
+        leg_flow_FHCal_Pt_uncorr[cent][i][j]->AddEntry(pFlowPt_uncorrected[cent][i][j], "Non-uniform acc. (uncorr.)", "p");
         leg_flow_FHCal_Pt_uncorr[cent][i][j]->Draw();
 
         canv_flow_FHCal_Pt_uncorr[cent][i][j]->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/v%i_pt_FHCalEp_cent%i_pid%i_with_uncorrected.png", i + 1, cent, j));
@@ -701,8 +748,8 @@ void Draw_AcceptanceComparison()
       leg_flow_Tpc_Pt[cent][i] = new TLegend(0.15, 0.72, 0.42, 0.89);
       canv_flow_Tpc_Pt[cent][i]->cd();
 
-      pFlowTpcPt[cent][0][i]->GetXaxis()->SetTitle("p_{T}, GeV/c");
-      pFlowTpcPt[cent][0][i]->GetYaxis()->SetTitle(Form("v_{%i}(Tpc EP)", 2));
+      pFlowTpcPt[cent][0][i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
+      pFlowTpcPt[cent][0][i]->GetYaxis()->SetTitle(Form("v_{%i}{#Psi_{2,TPC}}", 2));
 
       pFlowTpcPt[cent][0][i]->GetXaxis()->SetRangeUser(x_pt_min[1][cent][i], x_pt_max);
       pFlowTpcPt[cent][0][i]->GetYaxis()->SetRangeUser(y_pt_min[1][cent][i], 1.2 * y_pt_max[1][cent][i]);
@@ -728,9 +775,9 @@ void Draw_AcceptanceComparison()
       // pFlowTpcPt_uncorrected[cent][i]->Draw("p same");
 
       leg_flow_Tpc_Pt[cent][i]->SetHeader(Form("%s, %s", centFlowName[cent].Data(), pidName[i].Data()));
-      leg_flow_Tpc_Pt[cent][i]->AddEntry(pFlowTpcPt[cent][0][i], "Good Acceptance", "p");
-      leg_flow_Tpc_Pt[cent][i]->AddEntry(pFlowTpcPt[cent][1][i], "Bad Acceptance", "p");
-      // leg_flow_Tpc_Pt[cent][i]->AddEntry(pFlowTpcPt_uncorrected[cent][i], "Bad Acceptance (uncorr.)", "p");
+      leg_flow_Tpc_Pt[cent][i]->AddEntry(pFlowTpcPt[cent][0][i], "Uniform acc.", "p");
+      leg_flow_Tpc_Pt[cent][i]->AddEntry(pFlowTpcPt[cent][1][i], "Non-uniform acc.", "p");
+      // leg_flow_Tpc_Pt[cent][i]->AddEntry(pFlowTpcPt_uncorrected[cent][i], "Non-uniform acc. (uncorr.)", "p");
       leg_flow_Tpc_Pt[cent][i]->Draw();
 
       canv_flow_Tpc_Pt[cent][i]->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/v%i_pt_TpcEp_cent%i_pid%i.png", 2, cent, i));
@@ -739,22 +786,22 @@ void Draw_AcceptanceComparison()
       leg_flow_Tpc_Pt_uncorr[cent][i] = new TLegend(0.15, 0.72, 0.42, 0.89);
       canv_flow_Tpc_Pt_uncorr[cent][i]->cd();
 
-      pFlowTpcPt[cent][0][i]->GetXaxis()->SetTitle("p_{T}, GeV/c");
-      pFlowTpcPt[cent][0][i]->GetYaxis()->SetTitle(Form("v_{%i}(Tpc EP)", 2));
+      pFlowTpcPt[cent][0][i]->GetXaxis()->SetTitle("p_{T} (GeV/c)");
+      pFlowTpcPt[cent][0][i]->GetYaxis()->SetTitle(Form("v_{%i}{#Psi_{2,TPC}}", 2));
 
       pFlowTpcPt[cent][0][i]->GetXaxis()->SetRangeUser(x_pt_min[1][cent][i], x_pt_max);
       pFlowTpcPt[cent][0][i]->GetYaxis()->SetRangeUser(0.9 * y_pt_min[1][cent][i], 2.2 * y_pt_max[1][cent][i]);
 
-      pFlowTpcPt[cent][0][i]->SetMarkerStyle(20);
-      pFlowTpcPt[cent][1][i]->SetMarkerStyle(25);
+      pFlowTpcPt[cent][0][i]->SetMarkerStyle(25);
+      pFlowTpcPt[cent][1][i]->SetMarkerStyle(20);
 
       pFlowTpcPt[cent][0][i]->SetMarkerSize(1.6);
       pFlowTpcPt[cent][1][i]->SetMarkerSize(1.6);
 
-      pFlowTpcPt[cent][0][i]->SetMarkerColor(kBlack);
-      pFlowTpcPt[cent][1][i]->SetMarkerColor(kRed);
-      pFlowTpcPt[cent][0][i]->SetLineColor(kBlack);
-      pFlowTpcPt[cent][1][i]->SetLineColor(kRed);
+      pFlowTpcPt[cent][0][i]->SetMarkerColor(kRed);
+      pFlowTpcPt[cent][1][i]->SetMarkerColor(kBlack);
+      pFlowTpcPt[cent][0][i]->SetLineColor(kRed);
+      pFlowTpcPt[cent][1][i]->SetLineColor(kBlack);
 
       pFlowTpcPt_uncorrected[cent][i]->SetMarkerStyle(26);
       pFlowTpcPt_uncorrected[cent][i]->SetMarkerSize(1.6);
@@ -762,13 +809,13 @@ void Draw_AcceptanceComparison()
       pFlowTpcPt_uncorrected[cent][i]->SetLineColor(kBlue);
 
       pFlowTpcPt[cent][0][i]->Draw("p");
-      pFlowTpcPt[cent][1][i]->Draw("p same");
+      // pFlowTpcPt[cent][1][i]->Draw("p same");
       pFlowTpcPt_uncorrected[cent][i]->Draw("p same");
 
       leg_flow_Tpc_Pt_uncorr[cent][i]->SetHeader(Form("%s, %s", centFlowName[cent].Data(), pidName[i].Data()));
-      leg_flow_Tpc_Pt_uncorr[cent][i]->AddEntry(pFlowTpcPt[cent][0][i], "Good Acceptance", "p");
-      leg_flow_Tpc_Pt_uncorr[cent][i]->AddEntry(pFlowTpcPt[cent][1][i], "Bad Acceptance", "p");
-      leg_flow_Tpc_Pt_uncorr[cent][i]->AddEntry(pFlowTpcPt_uncorrected[cent][i], "Bad Acceptance (uncorr.)", "p");
+      leg_flow_Tpc_Pt_uncorr[cent][i]->AddEntry(pFlowTpcPt[cent][0][i], "Uniform acc.", "p");
+      // leg_flow_Tpc_Pt_uncorr[cent][i]->AddEntry(pFlowTpcPt[cent][1][i], "Non-uniform acc.", "p");
+      leg_flow_Tpc_Pt_uncorr[cent][i]->AddEntry(pFlowTpcPt_uncorrected[cent][i], "Non-uniform acc. (uncorr.)", "p");
       leg_flow_Tpc_Pt_uncorr[cent][i]->Draw();
 
       canv_flow_Tpc_Pt_uncorr[cent][i]->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/v%i_pt_TpcEp_cent%i_pid%i_with_uncorrected.png", 2, cent, i));
@@ -782,11 +829,11 @@ void Draw_AcceptanceComparison()
       for (int j = 0; j < npid; j++)
       {
         vgraph.clear();
-        pFlowPt[cent][0][i][j]->SetTitle("Good Acceptance");
-        pFlowPt[cent][1][i][j]->SetTitle("Bad Acceptance");
+        pFlowPt[cent][0][i][j]->SetTitle("Uniform acc.");
+        pFlowPt[cent][1][i][j]->SetTitle("Non-uniform acc.");
         vgraph.push_back(ConvertProfile(pFlowPt[cent][0][i][j]));
         vgraph.push_back(ConvertProfile(pFlowPt[cent][1][i][j]));
-        canvRatio = (TCanvas *)DrawTGraph(vgraph, Form("%s, %s", centFlowName[cent].Data(), pidName[j].Data()), 0.79, 1.21, x_pt_min[i][cent][j], x_pt_max, y_pt_min[i][cent][j], 1.2 * y_pt_max[i][cent][j], 0.18, 0.72, 0.45, 0.89);
+        canvRatio = (TCanvas *)DrawTGraph(vgraph, Form("Au+Au at #sqrt{s_{NN}}=7.7 GeV, UrQMD, GEANT4, %s, %s", centFlowName[cent].Data(), pidName[j].Data()), 0.91, 1.09, 0.*x_pt_min[i][cent][j], x_pt_max, y_pt_min[i][cent][j], 1.2 * y_pt_max[i][cent][j], 0.17, 0.67, 0.77, 0.89);
         canvRatio->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/Ratio_v%i_pt_FHCalEp_cent%i_pid%i.png", i + 1, cent, j));
         delete canvRatio;
       }
@@ -798,11 +845,11 @@ void Draw_AcceptanceComparison()
     for (int j = 0; j < npid; j++)
     {
       vgraph.clear();
-      pFlowTpcPt[cent][0][j]->SetTitle("Good Acceptance");
-      pFlowTpcPt[cent][1][j]->SetTitle("Bad Acceptance");
+      pFlowTpcPt[cent][0][j]->SetTitle("Uniform acc.");
+      pFlowTpcPt[cent][1][j]->SetTitle("Non-uniform acc.");
       vgraph.push_back(ConvertProfile(pFlowTpcPt[cent][0][j]));
       vgraph.push_back(ConvertProfile(pFlowTpcPt[cent][1][j]));
-      canvRatio = (TCanvas *)DrawTGraph(vgraph, Form("%s, %s", centFlowName[cent].Data(), pidName[j].Data()), 0.79, 1.21, x_pt_min[1][cent][j], x_pt_max, 0.9 * y_pt_min[1][cent][j], 1.2 * y_pt_max[1][cent][j], 0.18, 0.72, 0.45, 0.89);
+      canvRatio = (TCanvas *)DrawTGraph(vgraph, Form("Au+Au at #sqrt{s_{NN}}=7.7 GeV, UrQMD, GEANT4, %s, %s", centFlowName[cent].Data(), pidName[j].Data()), 0.91, 1.09, 0.*x_pt_min[1][cent][j], x_pt_max, 0.9 * y_pt_min[1][cent][j], 1.2 * y_pt_max[1][cent][j], 0.17, 0.67, 0.77, 0.89);
       canvRatio->SaveAs(Form("/home/peter/Documents/AcceptanceCompare/Compare/Ratio_v%i_pt_TpcEp_cent%i_pid%i.png", 2, cent, j));
       delete canvRatio;
     }