Data and analysis products for: Broadband Spectral Energy Distributions of Five LHAASO Ultra-High-Energy Sources Associated with Supernova Remnant and Molecular Cloud Systems

This archive contains the data, analysis code, and figure source behind the accepted ApJ paper "Broadband Spectral Energy Distributions of Five LHAASO Ultra-High-Energy Sources Associated with Supernova Remnant and Molecular Cloud Systems" (AAS81063). Version 1.2 is the post-acceptance release. Contents. The archive includes: (1) the complete Fermi-LAT/LHAASO flux-point table and fitted parameters for Tables 4 and 5; (2) thirty-bin Fermi-LAT profile-likelihood products and the converter that reproduces the SED table byte-for-byte; (3) the controlled PowerLaw/LogParabola/PLSuperExpCutoff4 comparison used for Table 3; (4) the full charged-pion decay-chain neutrino calculation and Table 6 outputs; (5) the band-limited Swift-XRT synchrotron/inverse-Compton calculation; (6) human-readable CSV versions of the principal tables; and (7) the five final vector figure PDFs. Reproduction. "bash code/run_reproduce.sh" refits the single-zone models from enclosed inputs and regenerates the two main-text figures. "python code/run_auxiliary_checks.py" regenerates the neutrino and Swift-XRT products. The Fermi SED converter can be checked with its --verify option and reproducesthe released table byte-for-byte. The model-selection driver requires the public Fermi event data and original ROI state, which are not redistributed; its exact results, configuration, and log are archived. Software. The spectral modelling used Python 3.9 with NumPy 1.26, SciPy 1.13, Matplotlib 3.9, Astropy 6.0 and Naima 0.10.0. Fermi-LAT products used Fermitools 2.5.1 and fermipy 1.4.2 with P8R3_SOURCE_V3, 4FGL-DR4, gll_iem_v07 and iso_P8R3_SOURCE_V3_v1. Auxiliary neutrino/X-ray scripts are compatible with NumPy 1.26 and later. Provenance and caveats. LHAASO values are from Cao et al. (2024, ApJS, 271, 25; VizieR J/ApJS/271/25). Fermi event data and response files are public but are not redistributed. Swift-XRT limits use public LSXPS stacks. IceCube limits are from Abbasi et al. (2023, ApJL, 945, L8). Correlated LHAASO pairs and shared-background Fermi bins are treated as described in the paper. Point-source Swift limits do not constrain extended pulsar-wind nebulae or displaced emission. License. CC BY 4.0.