From 3fd751d378d0996b5ba1cb5eb9146965e289c0ec Mon Sep 17 00:00:00 2001
From: Tamas Gal <tgal@km3net.de>
Date: Thu, 4 Jul 2019 13:40:28 +0200
Subject: [PATCH] Add royprefit

---
 scripts/live_royfit.jl | 10 ++++++++++
 1 file changed, 10 insertions(+)

diff --git a/scripts/live_royfit.jl b/scripts/live_royfit.jl
index a9f1245..2409191 100755
--- a/scripts/live_royfit.jl
+++ b/scripts/live_royfit.jl
@@ -25,7 +25,9 @@ function main()
         hits = calibrate(event.hits, calib)
         triggered_hits = filter(h->h.triggered, hits)
         dus = sort(unique(map(h->h.du, hits)))
+        triggered_dus = sort(unique(map(h->h.du, triggered_hits)))
         n_dus = length(dus)
+        n_triggered_dus = length(triggered_dus)
         n_doms = length(unique(h->h.dom_id, triggered_hits))
 
         if n_doms < 4
@@ -58,6 +60,14 @@ function main()
 
             savefig("plots/ztplot_roy.png")
         end
+
+        if n_triggered_dus > 2
+            println("Starting multiline fit with $(n_dus) DUs")
+            prefit_track = KM3NeT.prefit(triggered_hits)
+            println(prefit_track)
+            plot(triggered_hits, prefit_track)
+            savefig("plots/ztplot_roy_prefit.png")
+        end
     end
 end
 
-- 
GitLab