Flowchart of the fitting procedure. The leptohadronic fit of the data from each source is performed as a three-step process. The initial step, where the optical and gamma ray broadband features are explained with a leptonic model, requires the most simulations because the vast parameter space is initially unconstrained. The number of CPU-hours quoted refers to the entire sample and is based on the code’s performance of 1.6 s for a single leptonic simulation and ~20 s for a leptohadronic simulation.