Skip to content
Snippets Groups Projects
Commit aed0c56e authored by Wouter Deconinck's avatar Wouter Deconinck
Browse files

Add logQ2 plot for large minQ2 simulations

parent bcae96ed
No related branches found
No related tags found
No related merge requests found
This commit is part of merge request !177. Comments created here will be created in the context of that merge request.
...@@ -82,6 +82,12 @@ int dis_electrons(const std::string& config_name) ...@@ -82,6 +82,12 @@ int dis_electrons(const std::string& config_name)
.Define("Q2_da", "InclusiveKinematicsDA.Q2") .Define("Q2_da", "InclusiveKinematicsDA.Q2")
.Define("Q2_sigma", "InclusiveKinematicsSigma.Q2") .Define("Q2_sigma", "InclusiveKinematicsSigma.Q2")
.Define("Q2_esigma", "InclusiveKinematicseSigma.Q2") .Define("Q2_esigma", "InclusiveKinematicseSigma.Q2")
.Define("Q2_sim", "log(Q2_sim)")
.Define("Q2_el", "log(Q2_el)")
.Define("Q2_jb", "log(Q2_jb)")
.Define("Q2_da", "log(Q2_da)")
.Define("Q2_sigma", "log(Q2_sigma)")
.Define("Q2_esigma", "log(Q2_esigma)")
.Define("Q2_el_res", combinatorial_diff_ratio, {"Q2_sim", "Q2_el"}) .Define("Q2_el_res", combinatorial_diff_ratio, {"Q2_sim", "Q2_el"})
.Define("Q2_jb_res", combinatorial_diff_ratio, {"Q2_sim", "Q2_jb"}) .Define("Q2_jb_res", combinatorial_diff_ratio, {"Q2_sim", "Q2_jb"})
.Define("Q2_da_res", combinatorial_diff_ratio, {"Q2_sim", "Q2_da"}) .Define("Q2_da_res", combinatorial_diff_ratio, {"Q2_sim", "Q2_da"})
...@@ -107,6 +113,12 @@ int dis_electrons(const std::string& config_name) ...@@ -107,6 +113,12 @@ int dis_electrons(const std::string& config_name)
auto h_Q2_da = d0.Histo1D({"h_Q2_da", "; GeV^2; counts", 100, -5, 25}, "Q2_da"); auto h_Q2_da = d0.Histo1D({"h_Q2_da", "; GeV^2; counts", 100, -5, 25}, "Q2_da");
auto h_Q2_sigma = d0.Histo1D({"h_Q2_sigma", "; GeV^2; counts", 100, -5, 25}, "Q2_sigma"); auto h_Q2_sigma = d0.Histo1D({"h_Q2_sigma", "; GeV^2; counts", 100, -5, 25}, "Q2_sigma");
auto h_Q2_esigma = d0.Histo1D({"h_Q2_esigma", "; GeV^2; counts", 100, -5, 25}, "Q2_esigma"); auto h_Q2_esigma = d0.Histo1D({"h_Q2_esigma", "; GeV^2; counts", 100, -5, 25}, "Q2_esigma");
auto h_logQ2_sim = d0.Histo1D({"h_logQ2_sim", "; GeV^2; counts", 100, -1, 4}, "logQ2_sim");
auto h_logQ2_el = d0.Histo1D({"h_logQ2_el", "; GeV^2; counts", 100, -1, 4}, "logQ2_el");
auto h_logQ2_jb = d0.Histo1D({"h_logQ2_jb", "; GeV^2; counts", 100, -1, 4}, "logQ2_jb");
auto h_logQ2_da = d0.Histo1D({"h_logQ2_da", "; GeV^2; counts", 100, -1, 4}, "logQ2_da");
auto h_logQ2_sigma = d0.Histo1D({"h_logQ2_sigma", "; GeV^2; counts", 100, -1, 4}, "logQ2_sigma");
auto h_logQ2_esigma = d0.Histo1D({"h_logQ2_esigma", "; GeV^2; counts", 100, -1, 4}, "logQ2_esigma");
auto h_Q2_el_res = d0.Histo1D({"h_Q2_el_res", "; ; counts", 100, -1, 1}, "Q2_el_res"); auto h_Q2_el_res = d0.Histo1D({"h_Q2_el_res", "; ; counts", 100, -1, 1}, "Q2_el_res");
auto h_Q2_jb_res = d0.Histo1D({"h_Q2_jb_res", "; ; counts", 100, -1, 1}, "Q2_jb_res"); auto h_Q2_jb_res = d0.Histo1D({"h_Q2_jb_res", "; ; counts", 100, -1, 1}, "Q2_jb_res");
auto h_Q2_da_res = d0.Histo1D({"h_Q2_da_res", "; ; counts", 100, -1, 1}, "Q2_da_res"); auto h_Q2_da_res = d0.Histo1D({"h_Q2_da_res", "; ; counts", 100, -1, 1}, "Q2_da_res");
...@@ -237,12 +249,12 @@ int dis_electrons(const std::string& config_name) ...@@ -237,12 +249,12 @@ int dis_electrons(const std::string& config_name)
c.cd(); c.cd();
gPad->SetLogx(false); gPad->SetLogx(false);
gPad->SetLogy(true); gPad->SetLogy(true);
auto& h1 = *h_Q2_sim; auto& h1 = *h_logQ2_sim;
auto& h2 = *h_Q2_el; auto& h2 = *h_logQ2_el;
auto& h3 = *h_Q2_jb; auto& h3 = *h_logQ2_jb;
auto& h4 = *h_Q2_da; auto& h4 = *h_logQ2_da;
auto& h5 = *h_Q2_sigma; auto& h5 = *h_logQ2_sigma;
auto& h6 = *h_Q2_esigma; auto& h6 = *h_logQ2_esigma;
// histogram style // histogram style
h1.SetLineColor(common_bench::plot::kMpBlue); h1.SetLineColor(common_bench::plot::kMpBlue);
h1.SetLineWidth(2); h1.SetLineWidth(2);
...@@ -292,7 +304,7 @@ int dis_electrons(const std::string& config_name) ...@@ -292,7 +304,7 @@ int dis_electrons(const std::string& config_name)
tptr1 = t1.AddText("e#Sigma method"); tptr1 = t1.AddText("e#Sigma method");
tptr1->SetTextColor(common_bench::plot::kMpCyan); tptr1->SetTextColor(common_bench::plot::kMpCyan);
t1.Draw(); t1.Draw();
c.Print(fmt::format("{}_Q2_panels.png", output_prefix).c_str()); c.Print(fmt::format("{}_logQ2_panels.png", output_prefix).c_str());
} }
// Q2 comparison (overlays) // Q2 comparison (overlays)
...@@ -301,12 +313,12 @@ int dis_electrons(const std::string& config_name) ...@@ -301,12 +313,12 @@ int dis_electrons(const std::string& config_name)
c.cd(); c.cd();
gPad->SetLogx(false); gPad->SetLogx(false);
gPad->SetLogy(true); gPad->SetLogy(true);
auto& h1 = *h_Q2_sim; auto& h1 = *h_logQ2_sim;
auto& h2 = *h_Q2_el; auto& h2 = *h_logQ2_el;
auto& h3 = *h_Q2_jb; auto& h3 = *h_logQ2_jb;
auto& h4 = *h_Q2_da; auto& h4 = *h_logQ2_da;
auto& h5 = *h_Q2_sigma; auto& h5 = *h_logQ2_sigma;
auto& h6 = *h_Q2_esigma; auto& h6 = *h_logQ2_esigma;
// histogram style // histogram style
h1.SetLineColor(common_bench::plot::kMpBlue); h1.SetLineColor(common_bench::plot::kMpBlue);
h1.SetLineWidth(2); h1.SetLineWidth(2);
...@@ -350,7 +362,7 @@ int dis_electrons(const std::string& config_name) ...@@ -350,7 +362,7 @@ int dis_electrons(const std::string& config_name)
tptr1 = t1.AddText("e#Sigma method"); tptr1 = t1.AddText("e#Sigma method");
tptr1->SetTextColor(common_bench::plot::kMpCyan); tptr1->SetTextColor(common_bench::plot::kMpCyan);
t1.Draw(); t1.Draw();
c.Print(fmt::format("{}_Q2_overlays.png", output_prefix).c_str()); c.Print(fmt::format("{}_logQ2_overlays.png", output_prefix).c_str());
} }
// Q2 resolution (panels) // Q2 resolution (panels)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment