A Unified Multi-Wavelength Data Analysis Workflow with gammapy. Constraining the Broadband Emission of FSRQ OP 313

Abstract

The Flat Spectrum Radio Quasar OP~313 entered an enhanced activity phase in November 2023 and has undergone multiple flares since then which have motivated the organisation of several large multi-wavelength campaigns, including two deep observations from the hard X-ray telescope NuSTAR. The broadband emission from OP~313 during these two observations is investigated under a new unified analysis framework, with data spanning from optical to gamma rays. Traditional methods for analyzing blazar emissions often rely on proprietary software tailored to specific instruments, making it challenging to integrate and interpret data from multi-wavelength campaigns comprehensively. This study demonstrates the feasibility of utilizing gammapy, an open-source Python package, and common data formats originally developed for gamma-ray instrumentation, to perform a consistent multi-instrument analysis. This enables a forward folding approach that fully incorporates source observations, detector responses, and various instrumental and astrophysical backgrounds. The methodology is applied as an example to recent data collected from the distant quasar OP~313. We present a comprehensive data reconstruction and analysis for instruments including the Liverpool Telescope's IO:O detector, Swift-UVOT, Swift-XRT, NuSTAR, and Fermi-LAT. The resulting spectral analysis is validated against the native tools for each instrument. Additionally, a multi-wavelength phenomenological model of the source emission, encompassing optical to gamma-ray bands, is developed, incorporating absorption components across different energy regimes. We introduce and validate a new unified framework for multi-wavelength forward folding data analysis based on gammapy and open data formats, demonstrating its application to spectral data from the quasar OP~313.

arXiv comment17 pages, 10 figures github: https://github.com/mireianievas/gammapy_mwl_workflow zenodo: https://zenodo.org/records/13837637

References in corpus (33)

• Array Programming with NumPy
• Maps of Dust IR Emission for Use in Estimation of Reddening and CMBR Foregrounds
• Astropy: A Community Python Package for Astronomy
• The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
• The Large Area Telescope on the Fermi Gamma-ray Space Telescope Mission
• On the Absorption of X-rays in the Interstellar Medium
• Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions
• The Nuclear Spectroscopic Telescope Array (NuSTAR) Mission
• The Swift Ultra-Violet/Optical Telescope
• HI4PI: A full-sky HI survey based on EBHIS and GASS
• The Sloan Digital Sky Survey Quasar Catalog V. Seventh Data Release
• astroquery: An Astronomical Web-Querying Package in Python
• Extragalactic Background Light Inferred from AEGIS Galaxy SED-type Fractions
• The Pan-STARRS1 Database and Data Products
• Incremental Fermi Large Area Telescope Fourth Source Catalog
• The Eighteenth Data Release of the Sloan Digital Sky Surveys: Targeting and First Spectra from SDSS-V
• Swift observations of the very intense flaring activity of Mrk 421 during 2006: I. Phenomenological picture of electron acceleration and predictions for the MeV/GeV emission
• Stochastic Acceleration and the Evolution of Spectral Distributions in SSC Sources: A Self Consistent Modeling of Blazars' Flares
• NuSTAR Observations of the Bullet Cluster: Constraints on Inverse Compton Emission
• Proton synchrotron radiation of large-scale jets in active galactic nuclei
• Swift observations of the X-ray bright GRB 050315
• An observational determination of the evolving extragalactic background light from the multiwavelength HST/CANDELS survey in the Fermi and CTA era
• Very-high-energy gamma-rays from the Universe's middle age: detection of the z=0.940 blazar PKS 1441+25 with MAGIC
• Measurement of the Extragalactic Background Light using MAGIC and Fermi-LAT gamma-ray observations of blazars up to z = 1
• Gammapy: A Python package for gamma-ray astronomy
• A lower bound on intergalactic magnetic fields from time variability of 1ES 0229+200 from MAGIC and Fermi/LAT observations
• Data Unfolding Methods in High Energy Physics
• Time-dependent lepto-hadronic modeling of the emission from blazar jets with SOPRANO: the case of TXS 0506+056, 3HSP J095507.9+355101 and 3C 279
• agnpy: An open-source python package modelling the radiative processes of jetted active galactic nuclei
• Understanding the spectrum of a distant blazar PKS 1424+240 and its implications
• Hunting extreme BL Lacertae blazars with Fermi-LAT
• Detailed Analysis of the TeV γ-Ray Sources 3HWC J1928+178, 3HWC J1930+188, and the New Source HAWC J1932+192
• The Multiwavelength Picture of GRB 221009A's Afterglow