Speaker
Description
In this work, we demonstrate the application of the Gammapy Python package for simultaneous multi-wavelength spectral fitting of low level data across multiple energy decades. We have integrated the B13 lepto-hadronic emission code into the Gammapy framework, enabling comprehensive multi-instrument data analysis. Additionally, we incorporate IceCube neutrino detector response functions based on publicly available effective areas and point-spread functions and energy dispersion matrices from the "All-sky point-source IceCube data: years 2008-2018" dataset. This allows a self consistent way to constrain model parameters through neutrino flux predictions or limits. Including neutrino data in our approach grants us the ability to assess the viability of different jet composition scenarios based on a comparison between neutrino count expectation and detection.