1. "Fig1_counts_P8_P302_Source_z100_w009_w478_healpix_o7_24bins.fits"
Data used in the analysis - photon counts in pixels.
We use Healpix pixels with nside = 128 (approximate size is 0.5 deg), ring ordering

2. "Fig2_Fig7_Fig9_Fig10_Profiles_data_source.yaml"
Differential energy flux E^2 dN / dE [GeV / cm^2 s sr] of the data 
in latitude, longitude, and energy bins 
(the file contains the arithmetic center of Gal. Lat. and Gal. Lon. bins
and logarithmic centers of the energy bins).
The energies are in GeV.

3. "Fig3_LowE_diffuse_model_0.3-1.0GeV_flux_source.fits"
Intensity of the diffuse model derived from the low energy data,
E^2 dN / dE [GeV / cm^2 s sr],
in 18 energy bins between 1 GeV and 1 TeV.

4. "Fig4_LowE_residual_0.3-1.0GeV_flux_source.fits"
Residual in the "low energy model" in intensity units,
E^2 dN / dE [GeV / cm^2 s sr],
in 18 energy bins between 1 GeV and 1 TeV.
The maps are in Healpix format with nside = 128, ring ordering.
The energy table is in GeV.

5. 
"Fig5_Rectangles_0.3-1.0GeV_flux_source.fits"
"Fig5_Rectangles_diffuse_model_0.3-1.0GeV_flux_source.fits"
"Fig5_Rectangles_plus_residual_0.3-1.0GeV_flux_source.fits"

The rectangles model of the bubbles, the diffuse model, 
and residual + rectangles model of the bubbles in intensity units,
E^2 dN / dE [GeV / cm^2 s sr],
in 18 energy bins between 1 GeV and 1 TeV.
The maps are in Healpix format with nside = 128, ring ordering.
The energy table is in GeV.


6. "Fig6" folder.
diffuse_model_flux_9years_Source_z100_refit_PS_model.fits
gas_model_flux_9years_Source_z100_refit_PS_model.fits
IC_model_flux_9years_Source_z100_refit_PS_model.fits
ps_model_flux_9years_Source_z100_refit_PS_model.fits
residual_counts_9years_Source_z100_refit_PS_resid_signal_bubbles_counts.fits

Diffuse model (interstellar emission model, i.e., no isotropic contribution), 
gas-correlated emission (pi0 + brems), IC model, contribution of point sources, 
and residual.
In units of dN / dE [1/MeV cm^2 s sr].
The energy table is in units of MeV.


7. Residual spectra in GLat and GLon bins:
Fig7_Fig9_Fig10_Fig11_Profiles_rectangles_plus_residual_source.yaml
Fig7_Fig9_Fig10_Profiles_GALPROP_residual_plus_FBs_plus_gNFW_source.yaml
Fig7_Fig9_Fig10_Profiles_lowE_residual_source.yaml

Spectra of residuals including the Fermi bubbles
for "low energy" model, "rectangles" model, and "Galprop" models.
The data has dimensions (N Lat bins, N Lon bins, N energy bins)
(the files contains the arithmetic center of Gal. Lat. and Gal. Lon. bins
and logarithmic centers of the energy bins).
The energies are in GeV.

Fig8_Finer_profile_lowE_residual_source.yaml
Same as above for the "rectangles" model with more GLon bins.

8. 
Fig10_Logpar_index_at_500GeV.yaml
Results of the log parabola fit in figure 10:
(n - 2) indices at 500 GeV with statistical uncertainties

9.
Fig11_IC_model_of_FB_rectangles_plus_residual_source_l=-5.0_b=-4.0.yaml
and other files starting with "Fig11_IC_model_of_FB_rectangles_plus_residual_source"

IC model of the gamma-ray emission from the FBs.
Parameters:
N_0 - normalization of the energy density of electrons at 1 GeV,
dN / dE dV [1 / GeV cm^3],
assuming that the electrons are near the GC.
gamma and E_cut - parameters of the electron spectrum (power-law with exponential cutoff).
N_0_los - integral of N_0 along 1.5 kpc - the assumed emission region near the GC,
in the units of dN / dE dA [1 / GeV cm^2].
Energies are in GeV.
Spectrum values represent the best-fit gamma-ray spectra in units of 
E^2 dN / dE [GeV / cm^2 s sr]

The center of the GLat and GLon bin is given at the end of the filename: "l=-5.0_b=-4.0"
The width in GLon is 10 deg, in GLat is 4 deg.


Fig11_pi0_model_of_FB_rectangles_plus_residual_source_l=-5.0_b=-4.0.yaml
Hadronic model of the FBs.
Parameters in the file have the same definition as in the IC model.


Fig11_Parametric_model_of_FB_rectangles_plus_residual_source_l=-5.0_b=-4.0.yaml
and similar files:
best-fit to the FB spectrum using power-law with exponential cutoff functions.


10.
Fig12_Min_max_rectangles_plaw_source_cutoff_l=-5.0_b=0.0.yaml
Fig12_Min_max_rectangles_plaw_source_cutoff_l=5.0_b=0.0.yaml

The value of the Min and Max models to the left (GLon = 5) 
and to the right (GLon = -5) of the GC
and the parameters of the fit to the min and max points with the 
power-law + exp cutoff functions.

Spectra are in the units of E^2 dN / dE [GeV / cm^2 s sr]
Energies are in GeV.


11.
Fig13_Summary_electron_spectrum_source_plaw_l=-5.0_b=0.0.yaml
Fig13_Summary_proton_spectrum_source_plaw_l=-5.0_b=0.0.yaml

Spectra for electrons and protons expressed in isotropic flul units
I = c / 4 pi * rho
to compare with the flux of electrons and protons on the Earth.


12.
Fig14_Fig15_max_min_models_CTA_KM3NeT_sensitivities.yaml

Points for CTA and KM3NeT sensitivities.

13.
FigC1_ISRF_energy_density_comparison.yaml

Interstellar radiation field (ISRF) densities volume averaged over 
R < 1.5 kpc and |z| < 0.3 kpc around the GC.
For different models of ISRF.

14.
Folder "irfs"
Instrument response functions for the data used in the analysis:
expcube_P8_P302_Source_z100_w009_w478_P8R2_SOURCE_V6_healpix_o7_24bins.fits
Exposure in 24 energy bins for the main data sample.

psf_P8_P302_Source_z100_w009_w478_P8R2_SOURCE_V6.fits
Point spread function in the direction of the GC in 600 theta bins and 100 energy bins.