Like in the generate_pixel_mask and cluster_hits function of hit_analysis, plotting is done separately from the analysis, with an independent plotting function from plot_utils. I guess this should be the general approach of handling plotting since it is much more conveinient to have separate functions in case one needs to plot something externally. Maybe we can implement this in the future for all functions of the analysis.