Abstract: Context. Ground-based gamma-ray astronomy is still a rather young field of research, with strong historical connections to particle physics. This is why most observations are conducted by experiments with proprietary data and analysis software, as is usual in the particle physics field. However, in recent years, this paradigm has been slowly shifting toward the development and use of open-source data formats and tools, driven by upcoming observatories such as the Cherenkov Telescope Array (CTA). In this context, a community-driven, shared data format (the gamma-astro-data-format, or GADF) and analysis tools such as Gammapy and ctools have been developed. So far, these efforts have been led by the Imaging Atmospheric Cherenkov Telescope community, leaving out other types of ground-based gamma-ray instruments. Aims. We aim to show that the data from ground particle arrays, such as the High-Altitude Water Cherenkov (HAWC) observatory, are also compatible with the GADF and can thus be fully analyzed using the related tools, in this case, Gammapy. Methods. We reproduced several published HAWC results using Gammapy and data products compliant with GADF standard. We also illustrate the capabilities of the shared format and tools by producing a joint fit of the Crab spectrum including data from six different gamma-ray experiments. Results. We find excellent agreement with the reference results, a powerful confirmation of both the published results and the tools involved. Conclusions. The data from particle detector arrays such as the HAWC observatory can be adapted to the GADF and thus analyzed with Gammapy. A common data format and shared analysis tools allow multi-instrument joint analysis and effective data sharing. To emphasize this, a sample of Crab nebula event lists is made public with this paper. Because of the complementary nature of pointing and wide-field instruments, this synergy will be distinctly beneficial for the joint scientific exploitation of future observatories such as the Southern Wide-field Gamma-ray Observatory and CTA.

Abstract: The High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory surveys the very high-energy sky in the 300 GeV to >100 TeV energy range. HAWC has detected two blazars above 11 sigma, Markarian 421 (Mrk 421) and Markarian 501 (Mrk 501). The observations are comprised of data taken in the period between 2015 June and 2018 July, resulting in similar to 1038 days of exposure. In this work, we report the time-averaged spectral analyses for both sources, above 0.5 TeV. Taking into account the flux attenuation due to the extragalactic background light, the intrinsic spectrum of Mrk 421 is described by a power law with an exponential energy cutoff with index alpha = 2.26 +/- (0.12)(stat)((+0.17)(-0.2))(sys) and energy cutoff E-c = 5.1 +/- (1.6)(stat)((+1.4)(-2.5))(sys) TeV, while the intrinsic spectrum of Mrk 501 is better described by a simple power law with index alpha = 2.61 +/- (0.11)(stat)((+)(0.01)(-0.07))(sys). The maximum energies at which the Mrk 421 and Mrk 501 signals are detected are 9 and 12 TeV, respectively. This makes these some of the highest energy detections to date for spectra averaged over years-long timescales. Since the observation of gamma radiation from blazars provides information about the physical processes that take place in their relativistic jets, it is important to study the broadband spectral energy distributions (SEDs) of these objects. For this purpose, contemporaneous data in the gamma-ray band to the X-ray range, and literature data in the radio to UV range, were used to build time-averaged SEDs that were modeled within a synchrotron-self Compton leptonic scenario.