diff --git a/scripts/live_royfit.jl b/scripts/live_royfit.jl index a9f124583405a71c5b072b5d7a1dd7ebfdf242c3..24091914c48642171d4da003602e2a09f2747b5d 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