def plot_pt_j2(self):
dpt_j2_file = "/diphoton/j2pt.dat"
dpt_j2_plotter = RivetPlotter(
rivet_path=self.rivet_dir,
dat_file=dpt_j2_file,
)
xlabel = r"$p_{T,j_2}$ [GeV]"
if self.rescaling == "On":
ylabel = r"d$\sigma/$d$p_{T,j_2}$ [fb GeV$^{-1}$]"
elif self.rescaling == "XS":
ylabel = r"$1/\sigma$ d$\sigma/$d$p_{T,j_2}$ [GeV$^{-1}$]"
if self.legs == 6:
ylim = (1.2e-7, 1.8e-1)
else:
ylim = None
dpt_j2_fig = dpt_j2_plotter.plot_errors(
xlabel=xlabel,
ylabel=ylabel,
rescaling=self.rescaling,
training_reruns=self.training_reruns,
xlim=(0, 500),
ylim=ylim,
)
dpt_j2_fig.savefig(self.save_dir + "dpt_j2.pdf", bbox_inches="tight")
def plot_m_yy(self):
dm_yy_file = self.rivet_dir + "mass.dat"
dm_yy_plotter = RivetPlotter(dm_yy_file)
xlabel = r"$m_{\gamma_1\gamma_2}$"
if self.rescaling == "On":
ylabel = r"d$\sigma/$d$m_{\gamma_1\gamma_2}$ [fb]"
elif self.rescaling == "XS":
ylabel = r"$1/\sigma$ d$\sigma/$d$m_{\gamma_1\gamma_2}$"
if self.legs == 5:
ylim = (11e-7, 3)
else:
ylim = None
dm_yy_fig = dm_yy_plotter.plot(
xlabel=xlabel,
ylabel=ylabel,
ylim=ylim,
rescaling=self.rescaling,
)
dm_yy_fig.savefig(self.save_dir + "dm_yy.png",
dpi=self.dpi,
bbox_inches="tight")
def plot_m_yy(self):
dm_yy_file = "/diphoton/mass.dat"
dm_yy_plotter = RivetPlotter(
rivet_path=self.rivet_dir,
dat_file=dm_yy_file,
)
xlabel = r"$m_{\gamma_1\gamma_2}$"
if self.rescaling == "On":
ylabel = r"d$\sigma/$d$m_{\gamma_1\gamma_2}$ [fb]"
elif self.rescaling == "XS":
ylabel = r"$1/\sigma$ d$\sigma/$d$m_{\gamma_1\gamma_2}$"
if self.legs == 5:
ylim = (11e-7, 3)
else:
ylim = None
dm_yy_fig = dm_yy_plotter.plot_errors(
xlabel=xlabel,
ylabel=ylabel,
ylim=ylim,
training_reruns=self.training_reruns,
rescaling=self.rescaling,
)
dm_yy_fig.savefig(self.save_dir + "dm_yy.pdf", bbox_inches="tight")
def plot_phi_jj(self):
dphi_jj_file = "/diphoton/dphijj.dat"
dphi_jj_plotter = RivetPlotter(rivet_path=self.rivet_dir,
dat_file=dphi_jj_file)
xlabel = r"$\Delta\phi_{j_1j_2}$ [rad]"
if self.rescaling == "On":
ylabel = r"d$\sigma/$d$\Delta\phi_{j_1j_2}$ [fb rad$^{-1}$]"
elif self.rescaling == "XS":
ylabel = r"$1/\sigma$ d$\sigma/$d$\Delta\phi_{j_1j_2}$ [rad$^{-1}$]"
dphi_jj_fig = dphi_jj_plotter.plot_errors(
xlabel=xlabel,
ylabel=ylabel,
training_reruns=self.training_reruns,
rescaling=self.rescaling,
)
dphi_jj_fig.savefig(self.save_dir + "dphi_jj.pdf", bbox_inches="tight")
def plot_phi_jj(self):
dphi_jj_file = self.rivet_dir + "dphijj.dat"
dphi_jj_plotter = RivetPlotter(dphi_jj_file)
xlabel = r"$\Delta\phi_{j_1j_2}$ [rad]"
if self.rescaling == "On":
ylabel = r"d$\sigma/$d$\Delta\phi_{j_1j_2}$ [fb rad$^{-1}$]"
elif self.rescaling == "XS":
ylabel = r"$1/\sigma$ d$\sigma/$d$\Delta\phi_{j_1j_2}$ [rad$^{-1}$]"
dphi_jj_fig = dphi_jj_plotter.plot(
xlabel=xlabel,
ylabel=ylabel,
rescaling=self.rescaling,
)
dphi_jj_fig.savefig(self.save_dir + "dphi_jj.png",
dpi=self.dpi,
bbox_inches="tight")
def plot_eta_yy(self):
deta_yy_file = self.rivet_dir + "etayy.dat"
deta_yy_plotter = RivetPlotter(deta_yy_file)
xlabel = r"$\Delta\eta_{\gamma_1\gamma_2}$"
if self.rescaling == "On":
ylabel = r"d$\sigma/$d$\Delta\eta_{\gamma_1\gamma_2}$ [fb]"
elif self.rescaling == "XS":
ylabel = r"$1/\sigma$ d$\sigma/$d$\Delta\eta_{\gamma_1\gamma_2}$"
deta_yy_fig = deta_yy_plotter.plot(
xlabel=xlabel,
ylabel=ylabel,
rescaling=self.rescaling,
)
deta_yy_fig.savefig(self.save_dir + "deta_yy.png",
dpi=self.dpi,
bbox_inches="tight")
def plot_r_jy(self):
dr_jy_file = self.rivet_dir + "rsepjy.dat"
dr_jy_plotter = RivetPlotter(dr_jy_file)
xlabel = r"$R_{j_1\gamma_1}$"
if self.rescaling == "On":
ylabel = r"d$\sigma/$d$R_{j_1\gamma_1}$ [fb]"
elif self.rescaling == "XS":
ylabel = r"$1/\sigma$ d$\sigma/$d$R_{j_1\gamma_1}$"
dr_jy_fig = dr_jy_plotter.plot(
xlabel=xlabel,
ylabel=ylabel,
xlim=(0.2, 5),
rescaling=self.rescaling,
)
dr_jy_fig.savefig(self.save_dir + "dr_jy.png",
dpi=self.dpi,
bbox_inches="tight")
def plot_pt_j2(self):
dpt_j2_file = self.rivet_dir + "j2pt.dat"
dpt_j2_plotter = RivetPlotter(dpt_j2_file)
xlabel = r"$p_{T,j_2}$ [GeV]"
if self.rescaling == "On":
ylabel = r"d$\sigma/$d$p_{T,j_2}$ [fb GeV$^{-1}$]"
elif self.rescaling == "XS":
ylabel = r"$1/\sigma$ d$\sigma/$d$p_{T,j_2}$ [GeV$^{-1}$]"
dpt_j2_fig = dpt_j2_plotter.plot(
xlabel=xlabel,
ylabel=ylabel,
rescaling=self.rescaling,
xlim=500,
)
dpt_j2_fig.savefig(self.save_dir + "dpt_j2.png",
dpi=self.dpi,
bbox_inches="tight")
def plot_eta_yy(self):
deta_yy_file = "/diphoton/etayy.dat"
deta_yy_plotter = RivetPlotter(
rivet_path=self.rivet_dir,
dat_file=deta_yy_file,
)
xlabel = r"$\Delta\eta_{\gamma_1\gamma_2}$"
if self.rescaling == "On":
ylabel = r"d$\sigma/$d$\Delta\eta_{\gamma_1\gamma_2}$ [fb]"
elif self.rescaling == "XS":
ylabel = r"$1/\sigma$ d$\sigma/$d$\Delta\eta_{\gamma_1\gamma_2}$"
deta_yy_fig = deta_yy_plotter.plot_errors(
xlabel=xlabel,
ylabel=ylabel,
training_reruns=self.training_reruns,
rescaling=self.rescaling,
)
deta_yy_fig.savefig(self.save_dir + "deta_yy.pdf", bbox_inches="tight")
def plot_r_jy(self):
dr_jy_file = "/diphoton/rsepjy.dat"
dr_jy_plotter = RivetPlotter(
rivet_path=self.rivet_dir,
dat_file=dr_jy_file,
)
xlabel = r"$R_{j_1\gamma_1}$"
if self.rescaling == "On":
ylabel = r"d$\sigma/$d$R_{j_1\gamma_1}$ [fb]"
elif self.rescaling == "XS":
ylabel = r"$1/\sigma$ d$\sigma/$d$R_{j_1\gamma_1}$"
dr_jy_fig = dr_jy_plotter.plot_errors(
xlabel=xlabel,
ylabel=ylabel,
xlim=(0.2, 5),
training_reruns=self.training_reruns,
rescaling=self.rescaling,
)
dr_jy_fig.savefig(self.save_dir + "dr_jy.pdf", bbox_inches="tight")