Skip to content
Snippets Groups Projects
Commit a05e5028 authored by Kipp Cannon's avatar Kipp Cannon
Browse files

gstlal_inspiral_plotsummary: add likelihood ratio vs SNR scatter plots

parent 6cbb09a5
No related branches found
No related tags found
No related merge requests found
......@@ -1120,6 +1120,18 @@ WHERE
def finish(self):
livetime = far.get_live_time(self.farsegs)
fig, axes = create_plot(x_label = "SNR", y_label = r"$\ln \Lambda$")
axes.semilogx(self.background_snr, self.background_ln_likelihood_ratio, "kx", label = "Background")
axes.legend(loc = "upper left")
axes.set_title(r"$\ln \Lambda$ vs.\ SNR Scatter Plot (Closed Box)")
yield fig, "lr_vs_snr", False
fig, axes = create_plot(x_label = "SNR", y_label = r"$\ln \Lambda$")
axes.semilogx(self.background_snr, self.background_ln_likelihood_ratio, "kx", label = "Background")
axes.semilogx(self.zerolag_snr, self.zerolag_ln_likelihood_ratio, "bx", label = "Zero-lag")
axes.legend(loc = "upper left")
axes.set_title(r"$\ln \Lambda$ vs.\ SNR Scatter Plot")
yield fig, "lr_vs_snr", True
for ln_likelihood_ratio, ifars, is_open_box in [(self.zerolag_ln_likelihood_ratio, self.zerolag_ifar, True), (self.background_ln_likelihood_ratio, self.background_ifar, False)]:
fig, axes = create_farplot(ifars, livetime, is_open_box)
if fig is not None:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment