Add ROy

scripts/online_reco.py

@@ -56,11 +56,14 @@ class ZenithDistribution(kp.Module):
         n_ok = n - np.count_nonzero(np.isnan(self.zeniths))
         fontsize = 16
 
+        roy_zeniths = -np.loadtxt("reco.txt", "r")[-self.max_events:]
+        roy_zeniths[np.abs(roy_zeniths) > .99] = np.nan
+        n_roy = len(roy_zeniths)
+
         fig, ax = plt.subplots(figsize=(16, 8))
-        ax.hist(self.zeniths, bins=180)
-        ax.set_title(
-            r"Zenith distribution of JGandalf track reconstructions"
-            "\nbased on %d reconstructed tracks out of %d events" % (n_ok, n))
+        ax.hist(self.zeniths, bins=180, label="JGandalf (last %d events)" % n)
+        ax.hist(roy_zeniths, bins=180, label="ROy (last %d events)" % n_roy)
+        ax.set_title("Zenith distribution of online track reconstructions")
         ax.set_xlabel(r"cos(zenith)", fontsize=fontsize)
         ax.set_ylabel("count", fontsize=fontsize)
         ax.tick_params(labelsize=fontsize)