Suppl. Fig. 7. The reconstructed shower angle for all single photon events. The individual signal and background event type categories added together form the unconstrained prediction.

Suppl. Fig. 7. The constrained covariance matrix for the reconstructed shower angle for all single-photon events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Suppl. Fig. 7. The unconstrained covariance matrix for reconstructed shower angle for all single-photon events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 3. The reconstructed number of protons, with a proton KE energy threshold of 35 MeV, for all single-photon events. The individual signal and background event type categories added together form the unconstrained prediction.

Fig. 4, Suppl. Fig. 8. The reconstructed shower energy for zero proton events binned from 0 to 600 MeV. The individual signal and background event type categories added together form the unconstrained prediction.

Fig. 4. The constrained covariance matrix for the reconstructed shower energy for zero proton events binned from 0 to 600 MeV. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 4. The unconstrained covariance matrix for reconstructed shower energy for zero proton events binned from 0 to 600 MeV. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 6. The reconstructed shower energy for zero proton events binned from 150 to 1250 MeV, with LEE prediction. The individual signal and background event type categories added together form the unconstrained prediction. The "LEE" category is the unfolded MiniBooNE LEE model prediction scaled under the neutrino-target nucleon interaction assumption.

Fig. 6. The constrained covariance matrix for the reconstructed shower energy for zero proton events binned from 150 to 1250 MeV. Both with and without the LEE model added. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Fig. 6. The unconstrained covariance matrix for reconstructed shower energy for zero proton events binned from 150 to 1250 MeV. Both with and without the LEE model added. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Suppl. Fig. 10. The reconstructed shower backwards projected distance for zero proton events. The individual signal and background event type categories added together form the unconstrained prediction.

Suppl. Fig. 9. The reconstructed shower angle for zero proton events. Includes LEE model prediction. The individual signal and background event type categories added together form the unconstrained prediction. The "LEE" category is the unfolded MiniBooNE LEE model prediction scaled under the neutrino-target nucleon interaction assumption.

Suppl. Fig. 9. The constrained covariance matrix for the reconstructed shower angle for zero proton events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

Suppl. Fig. 9. The unconstrained covariance matrix for reconstructed shower angle for zero proton events. The matrix shows uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Data statistical uncertainties are not included. An example of how to add Pearson data statistical uncertainties can be found in the example code repository.

1$\gamma$Xp shower energy and angle selection efficiencies for true 1$\gamma$Xp events. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. Energy and Angle values represent bin lower edges.

1$\gamma$0p shower energy and angle selection efficiencies for true single-photon events with no other true particles of any energy. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. Energy and Angle values represent bin lower edges.

Covariance matrix showing uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Pearson data statistical uncertainties have been included, and include small correlations due to events which can be selected by both WC and Pandora. The first four bins are the WC $1\gamma Np$, WC $1\gamma 0p$, Pandora $1\gamma 1p$, and Pandora $1\gamma 0p$ channels, and the next four sets of 16 bins are the NC $\pi^0 Np$, NC $\pi^0 0p$, $\nu_\mu$CC $Np$, and $\nu_\mu$CC $0p$ channels. This corresponds to Fig. 1 of the paper and Fig. 6 of the Supplemental Material.

Constrained WC $1\gamma Np$ energy distribution. There are 6 bins from 0 to 600 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained WC $1\gamma 0p$ energy distribution. There are 14 bins from 0 to 700 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained Pandora $1\gamma 1p$ energy distribution. There are 6 bins from 0 to 600 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained Pandora $1\gamma 0p$ energy distribution. There are 14 bins from 0 to 700 MeV of reconstructed shower energy, with an additional overflow bin. Constrained systematic uncertainties on the predictions are illustrated, and a more detailed covariance matrix is included in the Constrained Energy Distributions Covariance Matrix tab. This corresponds to Fig. 2 of the paper.

Constrained covariance matrix showing constrained uncertainties and correlations between bins due to flux uncertainties, cross-section uncertainties, hadron reinteraction uncertainties, detector systematic uncertainties, Monte-Carlo statistical uncertainties, and dirt (outside cryostat) uncertainties. Pearson data statistical uncertainties have been included, and include small correlations due to events which can be selected by both WC and Pandora. The first seven bins are the WC $1\gamma Np$ reconstructed shower energy from 0-600 MeV with an overflow bin, the next 15 bins are the WC $1\gamma 0p$ reconstructed shower energy from 0-700 MeV with an overflow bin, the next seven bins are the Pandora $1\gamma 1p$ reconstructed shower energy from 0-600 MeV with an overflow bin, and the last 15 bins are the Pandora $1\gamma 0p$ reconstructed shower energy from 0-700 MeV with an overflow bin. This corresponds to Fig. 2 of the paper.

2D efficiency matrix showing the efficiency of the WC $1\gamma Np$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 15 bins from 0 to 1500 MeV in true photon energy. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the WC $1\gamma 0p$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the Pandora $1\gamma 1p$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 15 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the Pandora $1\gamma 0p$ selection with respect to all true NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied. This corresponds to Fig. 9 of the Supplemental Material.

2D efficiency matrix showing the efficiency of the WC $1\gamma 0p$ selection with respect to all true zero-primary-final-state-proton NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied.

2D efficiency matrix showing the efficiency of the Pandora $1\gamma 0p$ selection with respect to all true zero-primary-final-state-proton NC Delta Radiative Decay events which interact in the fiducial volume. The efficiency is shown as a function of true photon energy and angle. There are 10 bins from -1 to 1 in true photon $\cos \theta$ and 30 bins from 0 to 1500 MeV in true photon energy, with lower bin edges indicated in the table. NaN values indicate a "divide by zero" case where zero events were present in the bin before any selection was applied.

$a_2(1320)$ differential cross section in the m=-2 projection split into different reflectivity components, $\frac{d\sigma^+_{-2}}{dt}$ and $\frac{d\sigma^-_{-2}}{dt}$, in bins of four-momentum transfer. The first uncertainty is statistical, the second systematic.

$a_2(1320)$ differential cross section in the m=0 projection split into different reflectivity components, $\frac{d\sigma^+_{0}}{dt}$ and $\frac{d\sigma^-_{0}}{dt}$, in bins of four-momentum transfer. The first uncertainty is statistical, the second systematic.

$a_2(1320)$ differential cross section in the m=+1 projection split into different reflectivity components, $\frac{d\sigma^+_{+1}}{dt}$ and $\frac{d\sigma^-_{+1}}{dt}$, in bins of four-momentum transfer. The first uncertainty is statistical, the second systematic.

$a_2(1320)$ differential cross section in the m=+2 projection split into different reflectivity components, $\frac{d\sigma^+_{+2}}{dt}$ and $\frac{d\sigma^-_{+2}}{dt}$, in bins of four-momentum transfer. The first uncertainty is statistical, the second systematic.

Numerical data for $(K^+ + K^-)/2$ from Figure 1.

Numerical data for $K^0_S$ from Figure 2.

Numerical data for $K^0_S$ from Figure 2.

Numerical data for $(K^+ + K^-)/2$ from Figure 2.

Numerical data for $(K^+ + K^-)/2$ from Figure 2.

Numerical data for $K^0_S$ lifetime from Figure 6.

Numerical data for $K^0_S$ lifetime from Figure 6.

Numerical data for $K^0_S$ $y$-$p_T$ spectrum from Figure 7.

Full experimental (statistical and systematic) correlations (in \%) of the normalized partial decay rates over the $\overline{B}^0\to D^{*+} e^- \bar\nu_e$ and $\overline{B}^0\to D^{*+} \mu^- \bar\nu_\mu$ decays.

Measured partial decay rates $\Delta\Gamma$ (in units of $10^{-15}$ GeV) using the matrix inversion unfolding method

Average of normalized decay rates that are determined using the matrix inversion unfolding method

Full experimental (statistical and systematic) correlations (in \%) of the partial decay rates determined with the matrix inversion unfolding method for the $\overline{B}^0\to D^{*+} e^- \bar\nu_e$ and $\overline{B}^0\to D^{*+} \mu^- \bar\nu_\mu$ decays.

Full experimental (statistical and systematic) correlations (in \%) of the normalized partial decay rates determined using the matrix inversion unfolding method

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^+$ at 13$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^-$ at 13$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $p$ at 19$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\pi^+$ at 19$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\pi^-$ at 19$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^+$ at 19$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^-$ at 19$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $p$ at 30$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\bar{p}$ at 30$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\pi^+$ at 30$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\pi^-$ at 30$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^+$ at 30$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^-$ at 30$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $p$ at 40$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\bar{p}$ at 40$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\pi^+$ at 40$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\pi^-$ at 40$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^+$ at 40$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^-$ at 40$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $p$ at 75$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\bar{p}$ at 75$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\pi^+$ at 75$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\pi^-$ at 75$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^+$ at 75$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^-$ at 75$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $p$ at 150$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\bar{p}$ at 150$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\pi^+$ at 150$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $\pi^-$ at 150$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^+$ at 150$A$ GeV/c

Two-dimensional distributions ($y$ vs. $p_T$ ) of double differential yields of $K^-$ at 150$A$ GeV/c

Rapidity ($y$) spectra of $\pi^+$ at 13$A$ GeV/c

Rapidity ($y$) spectra of $\pi^-$ at 13$A$ GeV/c

Rapidity ($y$) spectra of $\pi^+$ at 19$A$ GeV/c

Rapidity ($y$) spectra of $\pi^-$ at 19$A$ GeV/c

Rapidity ($y$) spectra of $\pi^+$ at 30$A$ GeV/c

Rapidity ($y$) spectra of $\pi^-$ at 30$A$ GeV/c

Rapidity ($y$) spectra of $\pi^+$ at 40$A$ GeV/c

Rapidity ($y$) spectra of $\pi^-$ at 40$A$ GeV/c

Rapidity ($y$) spectra of $\pi^+$ at 75$A$ GeV/c

Rapidity ($y$) spectra of $\pi^-$ at 75$A$ GeV/c

Rapidity ($y$) spectra of $\pi^+$ at 150$A$ GeV/c

Rapidity ($y$) spectra of $\pi^-$ at 150$A$ GeV/c

The parameters of double-Gaussian fit to $\pi^+$ rapidity spectra

The parameters of double-Gaussian fit to $\pi^-$ rapidity spectra

Rapidity spectrum of negatively charged pions obtained with $h^-$ method in central Ar+Sc collisions at 30$A$ GeV/$c$ The results of $h^-$ analysis of Ar+Sc collisions used for comparison are published in reference [17]

Rapidity spectrum of negatively charged pions obtained with $h^-$ method in central Ar+Sc collisions at 150$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^+$ at 13$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^+$ at 19$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^+$ at 30$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^+$ at 40$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^+$ at 75$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^+$ at 150$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^-$ at 13$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^-$ at 19$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^-$ at 30$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^-$ at 40$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^-$ at 75$A$ GeV/$c$

Rapidity dependence of the inverse slope parameter $T$ fitted to $p_T$ spectrum of $K^-$ at 150$A$ GeV/$c$

Rapidity ($y$) spectra of $K^+$ at 13$A$ GeV/c

Rapidity ($y$) spectra of $K^-$ at 13$A$ GeV/c

Rapidity ($y$) spectra of $K^+$ at 19$A$ GeV/c

Rapidity ($y$) spectra of $K^-$ at 19$A$ GeV/c

Rapidity ($y$) spectra of $K^+$ at 30$A$ GeV/c

Rapidity ($y$) spectra of $K^-$ at 30$A$ GeV/c

Rapidity ($y$) spectra of $K^+$ at 40$A$ GeV/c

Rapidity ($y$) spectra of $K^-$ at 40$A$ GeV/c

Rapidity ($y$) spectra of $K^+$ at 75$A$ GeV/c

Rapidity ($y$) spectra of $K^-$ at 75$A$ GeV/c

Rapidity ($y$) spectra of $K^+$ at 150$A$ GeV/c

Rapidity ($y$) spectra of $K^-$ at 150$A$ GeV/c

The parameters of double-Gaussian fit to $K^+$ rapidity spectra

The parameters of double-Gaussian fit to $K^-$ rapidity spectra

Rapidity ($y$) spectra of $p$ at 13$A$ GeV/c

Rapidity ($y$) spectra of $p$ at 19$A$ GeV/c

Rapidity ($y$) spectra of $p$ at 30$A$ GeV/c

Rapidity ($y$) spectra of $\bar{p}$ at 30$A$ GeV/c

Rapidity ($y$) spectra of $p$ at 40$A$ GeV/c

Rapidity ($y$) spectra of $\bar{p}$ at 40$A$ GeV/c

Rapidity ($y$) spectra of $p$ at 75$A$ GeV/c

Rapidity ($y$) spectra of $\bar{p}$ at 75$A$ GeV/c

Rapidity ($y$) spectra of $p$ at 150$A$ GeV/c

Rapidity ($y$) spectra of $\bar{p}$ at 150$A$ GeV/c

The parameters of double-Gaussian fit to $\bar{p}$ rapidity spectra

The energy dependence of the inverse slope parameter $T(K^+)$ of $p_T$ spectra at mid-rapidity of positively charged $K$ mesons for ArSc. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the inverse slope parameter $T(K^+)$ of $p_T$ spectra at mid-rapidity of positively charged $K$ mesons for BeBe. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the inverse slope parameter $T(K^+)$ of $p_T$ spectra at mid-rapidity of positively charged $K$ mesons for ppWorld. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the inverse slope parameter $T(K^+)$ of $p_T$ spectra at mid-rapidity of positively charged $K$ mesons for pp. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the inverse slope parameter $T(K^+)$ of $p_T$ spectra at mid-rapidity of positively charged $K$ mesons for AuAu. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the inverse slope parameter $T(K^+)$ of $p_T$ spectra at mid-rapidity of positively charged $K$ mesons for PbPb. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the inverse slope parameter $T(K^-)$ of $p_T$ spectra at mid-rapidity of negatively charged $K$ mesons for ArSc. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the inverse slope parameter $T(K^-)$ of $p_T$ spectra at mid-rapidity of negatively charged $K$ mesons for BeBe. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the inverse slope parameter $T(K^-)$ of $p_T$ spectra at mid-rapidity of negatively charged $K$ mesons for ppWorld. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the inverse slope parameter $T(K^-)$ of $p_T$ spectra at mid-rapidity of negatively charged $K$ mesons for pp. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the inverse slope parameter $T(K^-)$ of $p_T$ spectra at mid-rapidity of negatively charged $K$ mesons for AuAu. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the inverse slope parameter $T(K^-)$ of $p_T$ spectra at mid-rapidity of negatively charged $K$ mesons for PbPb. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^+/\pi^+$ ratio at mid-rapidity for ArSc. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^+/\pi^+$ ratio at mid-rapidity for BeBe. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^+/\pi^+$ ratio at mid-rapidity for ppWorld. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^+/\pi^+$ ratio at mid-rapidity for pp. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^+/\pi^+$ ratio at mid-rapidity for AuAu. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^+/\pi^+$ ratio at mid-rapidity for PbPb. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^-/\pi^-$ ratio at mid-rapidity for ArSc. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^-/\pi^-$ ratio at mid-rapidity for BeBe. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^-/\pi^-$ ratio at mid-rapidity for ppWorld. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^-/\pi^-$ ratio at mid-rapidity for pp. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^-/\pi^-$ ratio at mid-rapidity for AuAu. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $K^-/\pi^-$ ratio at mid-rapidity for PbPb. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^+\rangle/\langle\pi^+\rangle$ mean multiplicity ratio for ArSc. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^+\rangle/\langle\pi^+\rangle$ mean multiplicity ratio for BeBe. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^+\rangle/\langle\pi^+\rangle$ mean multiplicity ratio for ppWorld. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^+\rangle/\langle\pi^+\rangle$ mean multiplicity ratio for pp. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^+\rangle/\langle\pi^+\rangle$ mean multiplicity ratio for AuAu. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^+\rangle/\langle\pi^+\rangle$ mean multiplicity ratio for PbPb. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^-\rangle/\langle\pi^-\rangle$ mean multiplicity ratio for ArSc. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^-\rangle/\langle\pi^-\rangle$ mean multiplicity ratio for BeBe. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^-\rangle/\langle\pi^-\rangle$ mean multiplicity ratio for ppWorld. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^-\rangle/\langle\pi^-\rangle$ mean multiplicity ratio for pp. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^-\rangle/\langle\pi^-\rangle$ mean multiplicity ratio for AuAu. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

The energy dependence of the $\langle K^-\rangle/\langle\pi^-\rangle$ mean multiplicity ratio for PbPb. The data used for comparison, other than Ar+Sc, are published in references [2-4,13,17,53-64]

Proton rapidity spectra in central Be+Be collisions at $19A$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in central Be+Be collisions at $30A$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in central Be+Be collisions at $40A$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in central Be+Be collisions at $75A$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in central Be+Be collisions at $150A$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in inelastic p+p collisions at $20$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in inelastic p+p collisions at $30$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in inelastic p+p collisions at $40$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in inelastic p+p collisions at $75$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in inelastic p+p collisions at $150$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in central Pb+Pb collisions at $20A$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in central Pb+Pb collisions at $30A$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in central Pb+Pb collisions at $40A$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in central Pb+Pb collisions at $80A$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

Proton rapidity spectra in central Pb+Pb collisions at $150A$ GeV/$c$ beam momentum. The data used for comparison, other than Ar+Sc, are published in references [2,13,17]

System size dependence of the $K^+/\pi^+$ ratio at mid-rapidity in $central$ nucleus-nucleus and inelastic $p$+$p$ interactions obtained at beam momenta of $\approx 150A$ GeV/$c$. The data used for comparison, other than Ar+Sc, are published in references [2,13,17,53].