From 354272be34196964f263af8b3cce66657afc96ab Mon Sep 17 00:00:00 2001
From: Tamas Gal <tgal@km3net.de>
Date: Sat, 9 Mar 2019 13:25:51 +0100
Subject: [PATCH] Update reco script
---
scripts/online_reco.py | 11 ++++++++---
1 file changed, 8 insertions(+), 3 deletions(-)
diff --git a/scripts/online_reco.py b/scripts/online_reco.py
index db31c8a..48352cb 100644
--- a/scripts/online_reco.py
+++ b/scripts/online_reco.py
@@ -57,8 +57,9 @@ class ZenithDistribution(kp.Module):
n_ok = n - np.count_nonzero(np.isnan(self.zeniths))
fontsize = 16
- roy_zeniths = -np.loadtxt("reco.txt")[-self.max_events:]
- roy_zeniths[np.abs(roy_zeniths) > .99] = np.nan
+ reco = np.loadtxt("reco2.txt", delimiter=',')[-self.max_events:]
+ roy_zeniths, roy_Q = -reco[:, 0], reco[:, 1]
+ roy_zeniths = roy_zeniths[roy_Q < 40]
n_roy = len(roy_zeniths)
fig, ax = plt.subplots(figsize=(16, 8))
@@ -66,18 +67,22 @@ class ZenithDistribution(kp.Module):
self.zeniths,
bins=180,
label="JGandalf (last %d events)" % n,
+ # label="JGandalf",
histtype="step",
+ normed=True,
lw=3)
ax.hist(
roy_zeniths,
bins=180,
label="ROy (last %d events)" % n_roy,
+ # label="ROy",
histtype="step",
+ normed=True,
lw=3)
ax.set_title("Zenith distribution of online track reconstructions\n%s"
% datetime.utcnow().strftime("%c"))
ax.set_xlabel(r"cos(zenith)", fontsize=fontsize)
- ax.set_ylabel("count", fontsize=fontsize)
+ ax.set_ylabel("normed count", fontsize=fontsize)
ax.tick_params(labelsize=fontsize)
ax.set_yscale("log")
plt.legend(fontsize=fontsize, loc=2)
--
GitLab