Off-by-one

reco/royfit.jl

@@ -58,7 +58,7 @@ function main()
             end
             fit = NeRCA.single_du_fit(du_hits, sparams)
             push!(Q, fit.Q)
-            plot!(du_hits, fit, markercolor=colours[idx % length(colours)], label="DU $(du)", max_z=calib.max_z)
+            plot!(du_hits, fit, markercolor=colours[idx % length(colours) + 1], label="DU $(du)", max_z=calib.max_z)
             write_time_residuals("/data/reco_timeres.csv", event, fit.selected_hits, fit)
         end
         if sum(Q) < 200 && n_doms > 6 && n_dus > 1