Numerical values of double-differential yields $d^{2}n/dydp_{T}$ presented in Fig. 10, given in units of $10^{−3} (GeV/c)^{−1}$. The first uncertainty is statistical, while the second one is systematic

Numerical values of double-differential yields $1/m_T d^{2}n/dm_{T}dy$ given in units of $10^{−3}(GeV)^{−2}$ and presented in Fig. 11; the values of $m_T−m_0$ specify the bin centers. The first uncertainty is statistical, while the second one is systematic

Numerical values of double-differential yields $1/m_T d^{2}n/dm_{T}dy$ given in units of $10^{−3}(GeV)^{−2}$ and presented in Fig. 11; the values of $m_T−m_0$ specify the bin centers. The first uncertainty is statistical, while the second one is systematic

Numerical values of double-differential yields $1/m_T d^{2}n/dm_{T}dy$ given in units of $10^{−3}(GeV)^{−2}$ and presented in Fig. 11; the values of $m_T−m_0$ specify the bin centers. The first uncertainty is statistical, while the second one is systematic

Numerical values of the $p_T$-integrated and extrapolated $dn/dy$ distribution presented in Fig. 12. The first uncertainty is statistical, while the second one is systematic. Additionally, the width of the Gaussian fit to the $dn/dy$ distribution, as well as the mean multiplicity of $K^{*}(892)^0$ mesons are shown (see the text for details)

The mean multiplicities $\langle K^{*}(892)^0\rangle$ and the widths of the rapidity distributions $\sigma_y$ obtained from $dn/dy$ distributions (see the text for details). The first uncertainty is statistical and the second systematic

The mean multiplicity of $K^{*}(892)^0$ mesons for 158 GeV/$c$ inelastic p+p interactions compared to theoretical multiplicities obtained within Hadron Gas Models

The mean multiplicities of different particle species measured in nucleus-nucleus collisions at 158$A$ GeV/$c$ by NA49 and NA61/SHINE. The total uncertainties of $\langle K^{*}(892)^0\rangle$, $\langle K^{+}\rangle$ and $\langle K^{-}\rangle$ were taken as the square roots of the sums of squares of statistical and systematic uncertainties. For NA49 p+p data, the $\langle K^{+}\rangle$ and $\langle K^{-}\rangle$ results include statistical uncertainties only($\langle K^{+}\rangle$ =0.2267$\pm$0.0006 and $\langle K^{-}\rangle$=0.1303$\pm$0.0004), whereas systematic uncertainties for total yields were not reported. Therefore, NA61/SHINE $\langle K^{+}\rangle$ and $\langle K^{-}\rangle$ values were used in the $\langle K^{*}(892)^0 \rangle / \langle K^{+} \rangle$ and $\langle K^{*}(892)^{0}\rangle / \langle K^{-} \rangle$ ratios. The numbers of $\langle K^{+} \rangle$ and $\langle K^{−}\rangle$ and their uncertainties for the 5% most central Pb+Pb collisions were multiplied by a factor 262/362 in orderto estimate charged kaon multiplicities in the 23.5% most central Pb+Pb reactions

Parameters of the polynomial fit of net-proton C3/C2 and C4/C2 vs energy in 0-5$\%$ central Au+Au collisions. Polynomial of order 5 and 4 are used to fit C3/C2 and C4/C2 ratio vs energies, respectively.

Cumulant ratios C3/C2 and C4/C2 of net-proton distributions in STAR Au+Au collisions for nine energies from $\sqrt{s_{NN}} = 7.7 - 200 GeV for 0-5$\%$ and 70-80$\%$ centrality.

Cumulant ratios C3/C2 and C4/C2 of net-proton distributions in Au+Au collisions from UrQMD for 0-5$\%$ centrality.

Cumulant ratios C3/C2 and C4/C2 of net-proton distributions from HRG GCE, HRG GEC excluded volume, HRG CE model.

Cumulants and their ratios up to fourth order in Au+Au collisions at $\sqrt{s_{NN}} = 200 GeV using conventional analytical efficiency corrections and unfolding approaches with various cases of distribution of detector response.

Rapidity dependence of cumulant ratio C4/C2 of net-proton distributions for 0-5$\%$ central Au+Au collisions at nine energies from $\sqrt{s_{NN}} = 7.7 - 200 GeV.

Significance of deviations in net-proton C4/C2 (or kappa*sigma^2) for central 0-5$\%$ Au+Au collisions and those from UrQMD, HRG and peripheral 70-80$\%$ data.

Sum of invariant mass spectra of XI-PI-, XI-PI+, XI+PI-, XI+PI+.

Sum of invariant mass spectra of XI-PI-, XI-PI+, XI+PI-, XI+PI+ with mix-background subtracted compared with NA49 result.

<K+>/<PI>+.

Inverse slope parameter T.

Inverse slope parameter T.

Inverse slope parameter T.

Inverse slope parameter T.

K-/PI- at y=0.

K-/PI- at y=0.

<K->/<PI->.

<K->/<PI->.

Data from Figure 32 – STEREO exclusion and exclusion sensitivity contours at 90% C.L. for 179 days reactor-on (phase-I+II). A full graphical presentation can be downloaded at "Resources" for reference.

Data from Figures 33 and 34 – STEREO exclusion and exclusion sensitivity contours at 95% C.L. for 179 days reactor-on (phase-I+II) using the two-dimensional method. A graphical presentation can be downloaded at "Resources" for reference.

Data from Figures 33 and 34 – STEREO exclusion and exclusion sensitivity contours at 95% C.L. for 179 days reactor-on (phase-I+II) using the two-dimensional method. A graphical presentation can be downloaded at "Resources" for reference.

$\Delta \chi^2$ map of phase-I+II data calculated by a raster-scan method. To be used in combination with the $\Delta \chi^2_{\text{crit},x}$ values to generate exclusion contours at $x$% C.L. for phase-I+II data. Additional explanations are given there.

Data from Figure 34 – STEREO exclusion and exclusion sensitivity contours at 95% C.L. for 179 days reactor-on (phase-I+II) using the raster-scan method. A graphical presentation can be downloaded at "Resources" for reference.

Data from Figure 34 – STEREO exclusion and exclusion sensitivity contours at 95% C.L. for 179 days reactor-on (phase-I+II) using the raster-scan method. A graphical presentation can be downloaded at "Resources" for reference.

The $\Delta \chi^2_{\text{crit},x}$ map accounts for the fact that the $\Delta \chi^2$ values of the oscillation fit do not follow a $\chi^2$ distribution with 2 degrees of freedom. Therefore, the $\Delta \chi^2_{\text{crit},x}$ value for x% C.L. of each point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ in the parameter space, $\Delta \chi^2_{\text{crit},x}(\sin^2(2\theta_{ee}), \Delta m^2_{41})$, is determined from the $\Delta \chi^2$ obtained in pseudo-experiments at that point, such that $\Delta \chi^2 \leq \Delta \chi^2_{\text{crit},x}$, for $x$% of the pseudo-experiments. Applying these $\Delta \chi^2_{\text{crit},x}$ values to the $\Delta \chi^2$ map obtained with the data, $x$% C.L. exclusion contours are obtained. The point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ is excluded by the data at $x$% C.L. if $\Delta \chi^2(\sin^2(2\theta_{ee}), \Delta m^2_{41}) > \Delta \chi^2_{\text{crit},x}(\sin^2(2\theta_{ee}), \Delta m^2_{41})$. In order to obtain the $\Delta \chi^2_{\text{crit},x}$ map, $10^4$ pseudo-experiments were generated for each point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ in the parameter space, taking into account all statistical and systematic uncertainties. The $\Delta \chi^2$ value of a pseudo-experiment is calculated by subtracting the $\chi^2$ value of the best-fit in the parameter space from the $\chi^2$ value of the fit at the $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ used in the generation of the pseudo-dataset, where all nuisance parameters are free within their pull terms. When combining the exclusion contours with other experimental data, special care should be exercised. The assumption of a standard $\chi^2$ law instead of the provided $\Delta \chi^2_{\text{crit},x}$ values derived from non-standard $\chi^2$ distributions leads to slightly modified contours. In addition, the contours were derived using a raster-scan in several fixed values of $\Delta m_{41}^2$. While this method is particularly suited to derive exclusion contours, it cannot be used to calculate allowed confidence regions for $\Delta m_{41}^2$ and consequently two-dimensional allowed confidence regions. This is because $\Delta \chi^2$ values are not reflecting the likelihood of individual $\Delta m_{41}^2$ values. Thus, a direct comparison of $\Delta \chi^2$ values across different $\Delta m_{41}^2$ values is not possible in a statistically meaningful way. When generating the exclusion contours with the aforementioned procedure, spurious exclusion regions at low values of $\sin^2(2\theta_{ee})$ can be encountered for some values of $\Delta m^2_{41}$. These should be ignored and are owed to the raster-scan procedure used to generate the maps.

Data from Figure 34 – STEREO exclusion and exclusion sensitivity contours at 95% C.L. for 179 days reactor-on (phase-I+II) using the CLs method. A graphical presentation can be downloaded at "Resources" for reference.

Data from Figure 34 – STEREO exclusion and exclusion sensitivity contours at 95% C.L. for 179 days reactor-on (phase-I+II) using the CLs method. A graphical presentation can be downloaded at "Resources" for reference.

$\Delta \chi^2$ map of phase-I+II data calculated by the two-dimensional method. To be used in combination with the $\Delta \chi^2_{\text{crit},x}$ values of the two-dimensional method to generate exclusion contours at $x$% C.L. for phase-I+II data. Additional explanations are given there.

$\Delta \chi^2$ map of phase-I+II data calculated by the two-dimensional method. To be used in combination with the $\Delta \chi^2_{\text{crit},x}$ values of the two-dimensional method to generate exclusion contours at $x$% C.L. for phase-I+II data. Additional explanations are given there.

The $\Delta \chi^2_{\text{crit},x}$ map accounts for the fact that the $\Delta \chi^2$ values of the oscillation fit do not follow a $\chi^2$ distribution with 2 degrees of freedom. Therefore, the $\Delta \chi^2_{\text{crit},x}$ value for x% C.L. of each point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ in the parameter space, $\Delta \chi^2_{\text{crit},x}(\sin^2(2\theta_{ee}), \Delta m^2_{41})$, is determined from the $\Delta \chi^2$ obtained in pseudo-experiments at that point, such that $\Delta \chi^2 \leq \Delta \chi^2_{\text{crit},x}$, for $x$% of the pseudo-experiments. Applying these $\Delta \chi^2_{\text{crit},x}$ values to the $\Delta \chi^2$ map obtained with the data, $x$% C.L. exclusion contours are obtained. The point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ is excluded by the data at $x$% C.L. if $\Delta \chi^2(\sin^2(2\theta_{ee}), \Delta m^2_{41}) > \Delta \chi^2_{\text{crit},x}(\sin^2(2\theta_{ee}), \Delta m^2_{41})$. In order to obtain the $\Delta \chi^2_{\text{crit},x}$ map, $10^4$ pseudo-experiments were generated for each point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ in the parameter space, taking into account all statistical and systematic uncertainties. The $\Delta \chi^2$ value of a pseudo-experiment is calculated by subtracting the $\chi^2$ value of the best-fit in the parameter space from the $\chi^2$ value of the fit at the $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ used in the generation of the pseudo-dataset, where all nuisance parameters are free within their pull terms. When combining the exclusion contours with other experimental data, special care should be exercised. The assumption of a standard $\chi^2$ law instead of the provided $\Delta \chi^2_{\text{crit},x}$ values derived from non-standard $\chi^2$ distributions leads to modified contours. In addition, for the interpretation of derived allowed confidence regions, one should always consider additional available information, in particular rate information. This is especially advisable considering the expected distribution of best-fit points in a prediction-free shape-only analysis, as discussed in the main publication.

The $\Delta \chi^2_{\text{crit},x}$ map accounts for the fact that the $\Delta \chi^2$ values of the oscillation fit do not follow a $\chi^2$ distribution with 2 degrees of freedom. Therefore, the $\Delta \chi^2_{\text{crit},x}$ value for x% C.L. of each point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ in the parameter space, $\Delta \chi^2_{\text{crit},x}(\sin^2(2\theta_{ee}), \Delta m^2_{41})$, is determined from the $\Delta \chi^2$ obtained in pseudo-experiments at that point, such that $\Delta \chi^2 \leq \Delta \chi^2_{\text{crit},x}$, for $x$% of the pseudo-experiments. Applying these $\Delta \chi^2_{\text{crit},x}$ values to the $\Delta \chi^2$ map obtained with the data, $x$% C.L. exclusion contours are obtained. The point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ is excluded by the data at $x$% C.L. if $\Delta \chi^2(\sin^2(2\theta_{ee}), \Delta m^2_{41}) > \Delta \chi^2_{\text{crit},x}(\sin^2(2\theta_{ee}), \Delta m^2_{41})$. In order to obtain the $\Delta \chi^2_{\text{crit},x}$ map, $10^4$ pseudo-experiments were generated for each point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ in the parameter space, taking into account all statistical and systematic uncertainties. The $\Delta \chi^2$ value of a pseudo-experiment is calculated by subtracting the $\chi^2$ value of the best-fit in the parameter space from the $\chi^2$ value of the fit at the $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ used in the generation of the pseudo-dataset, where all nuisance parameters are free within their pull terms. When combining the exclusion contours with other experimental data, special care should be exercised. The assumption of a standard $\chi^2$ law instead of the provided $\Delta \chi^2_{\text{crit},x}$ values derived from non-standard $\chi^2$ distributions leads to modified contours. In addition, for the interpretation of derived allowed confidence regions, one should always consider additional available information, in particular rate information. This is especially advisable considering the expected distribution of best-fit points in a prediction-free shape-only analysis, as discussed in the main publication.

$\Delta \chi^2$ map of phase-I+II data calculated by the raster-scan method. To be used in combination with the $\Delta \chi^2_{\text{crit},x}$ values of the raster-scan method to generate exclusion contours at $x$% C.L. for phase-I+II data. Additional explanations are given there.

$\Delta \chi^2$ map of phase-I+II data calculated by the raster-scan method. To be used in combination with the $\Delta \chi^2_{\text{crit},x}$ values of the raster-scan method to generate exclusion contours at $x$% C.L. for phase-I+II data. Additional explanations are given there.

The $\Delta \chi^2_{\text{crit},x}$ and $\Delta \chi^2$ values in the raster-scan method serve the same purpose as in the two-dimensional method. The explanations given there also apply here, with the following changes: Instead of the two free oscillation parameters $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$, the raster-scan method fits only $\sin^2(2\theta_{ee})$ as free oscillation parameter while repeating the fit for several fixed values of $\Delta m_{41}^2$. While this method is particularly suited to derive exclusion contours, it cannot be used to calculate allowed confidence regions for $\Delta m_{41}^2$ and consequently two-dimensional allowed confidence regions. This is because $\Delta \chi^2$ values are not reflecting the likelihood of individual $\Delta m_{41}^2$ values. Thus, a direct comparison of $\Delta \chi^2$ values across different $\Delta m_{41}^2$ values is not possible in a statistically meaningful way. Moreover, when generating the exclusion contours with the aforementioned procedure, spurious exclusion regions at low values of $\sin^2(2\theta_{ee})$ can be encountered for some values of $\Delta m^2_{41}$. These should be ignored and are owed to the raster-scan procedure used to generate the maps.

The $\Delta \chi^2_{\text{crit},x}$ and $\Delta \chi^2$ values in the raster-scan method serve the same purpose as in the two-dimensional method. The explanations given there also apply here, with the following changes: Instead of the two free oscillation parameters $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$, the raster-scan method fits only $\sin^2(2\theta_{ee})$ as free oscillation parameter while repeating the fit for several fixed values of $\Delta m_{41}^2$. While this method is particularly suited to derive exclusion contours, it cannot be used to calculate allowed confidence regions for $\Delta m_{41}^2$ and consequently two-dimensional allowed confidence regions. This is because $\Delta \chi^2$ values are not reflecting the likelihood of individual $\Delta m_{41}^2$ values. Thus, a direct comparison of $\Delta \chi^2$ values across different $\Delta m_{41}^2$ values is not possible in a statistically meaningful way. Moreover, when generating the exclusion contours with the aforementioned procedure, spurious exclusion regions at low values of $\sin^2(2\theta_{ee})$ can be encountered for some values of $\Delta m^2_{41}$. These should be ignored and are owed to the raster-scan procedure used to generate the maps.

$\Delta T$ map of phase-I+II data calculated by the CLs method. To be used in combination with the $\Delta T_{\text{crit},x}$ values of the CLs method to generate exclusion contours at $x$% C.L. for phase-I+II data. Additional explanations are given there.

$\Delta T$ map of phase-I+II data calculated by the CLs method. To be used in combination with the $\Delta T_{\text{crit},x}$ values of the CLs method to generate exclusion contours at $x$% C.L. for phase-I+II data. Additional explanations are given there.

The $\Delta T_{\text{crit},x}$ and $\Delta T$ values in the CLs method fulfil a similar purpose as the $\Delta \chi^2_{\text{crit},x}$ and $\Delta \chi^2$ values in the two-dimensional method. The explanations given there also apply here, with the following changes: The CLs method differs from the two-dimensional method as it normalises the confidence level of the oscillation-hypothesis to the confidence level of the null-hypothesis (i.e. no-oscillation-hypothesis), instead of the best-fit. The $\Delta T$ value of a dataset at the parameter space-point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ is computed by subtracting the $\chi^2$ value of the fit under the no-oscillation hypothesis from the $\chi^2$ value of the fit with the oscillation parameters fixed to $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$. With this definition, $\Delta T$ values can be positive or negative as opposed to $\Delta \chi^2$ values. Note that the values used to compare with are $-\Delta T$ and not $\Delta T$ directly. The CLs value of a parameter space-point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ is defined as $\text{CL}_\text{s}(\sin^2(2\theta_{ee}), \Delta m^2_{41}):=(1-p_1)/(1-p_0)$ where $(1-p_0)$ and $(1-p_1)$ are the confidence levels of the dataset, determined from the distributions of $-\Delta T$ created in pseudo-experiments assuming the oscillation parameters $[0,0]$ and $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$, respectively. The point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ is excluded by the data at $x$% C.L. if $\Delta T(\sin^2(2\theta_{ee}), \Delta m^2_{41}) < \Delta T_{\text{crit},x}(\sin^2(2\theta_{ee}), \Delta m^2_{41})$. Note the flipped sign with respect to the two-dimensional method and the raster-scan method. More information on the CLs method can be found at "Resources".

The $\Delta T_{\text{crit},x}$ and $\Delta T$ values in the CLs method fulfil a similar purpose as the $\Delta \chi^2_{\text{crit},x}$ and $\Delta \chi^2$ values in the two-dimensional method. The explanations given there also apply here, with the following changes: The CLs method differs from the two-dimensional method as it normalises the confidence level of the oscillation-hypothesis to the confidence level of the null-hypothesis (i.e. no-oscillation-hypothesis), instead of the best-fit. The $\Delta T$ value of a dataset at the parameter space-point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ is computed by subtracting the $\chi^2$ value of the fit under the no-oscillation hypothesis from the $\chi^2$ value of the fit with the oscillation parameters fixed to $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$. With this definition, $\Delta T$ values can be positive or negative as opposed to $\Delta \chi^2$ values. Note that the values used to compare with are $-\Delta T$ and not $\Delta T$ directly. The CLs value of a parameter space-point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ is defined as $\text{CL}_\text{s}(\sin^2(2\theta_{ee}), \Delta m^2_{41}):=(1-p_1)/(1-p_0)$ where $(1-p_0)$ and $(1-p_1)$ are the confidence levels of the dataset, determined from the distributions of $-\Delta T$ created in pseudo-experiments assuming the oscillation parameters $[0,0]$ and $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$, respectively. The point $[\sin^2(2\theta_{ee}), \Delta m^2_{41}]$ is excluded by the data at $x$% C.L. if $\Delta T(\sin^2(2\theta_{ee}), \Delta m^2_{41}) < \Delta T_{\text{crit},x}(\sin^2(2\theta_{ee}), \Delta m^2_{41})$. Note the flipped sign with respect to the two-dimensional method and the raster-scan method. More information on the CLs method can be found at "Resources".

Map of mean $\Delta \chi^2$ values derived from simulations of phase-I+II and calculated by the two-dimensional method. To be used analogously to the $\Delta \chi^2$ map of phase-I+II data in combination with the $\Delta \chi^2_{\text{crit},x}$ values of the two-dimensional method to generate sensitivity contours at $x$% C.L. for phase-I+II. Additional explanations are given there.

The rho distributions for the mixed-flavour category events corresponding to a signal hypothesis with mH = 261 GeV and mA = 150 GeV.

The rho distributions for the same-flavour category events corresponding to a signal hypothesis with mH = 442 GeV and mA = 193 GeV. The signal is normalised to its theoretical cross section.

The rho distributions for the mixed-flavour category events corresponding to a signal hypothesis with mH = 442 GeV and mA = 193 GeV.

Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction for H(A) -> ZA(H) -> bbbar as a function of mA and mH. The limits are computed using the asymptotic CLs method, combining the ee and μμ channels.

Expected and observed 95% CL upper limits on the signal strength for the Type-II 2HDM benchmark (tan(beta)=1.5, cos(beta-alpha)=0.01) as a function of mA and mH . The limits are computed using the asymptotic CLs method, combining the ee and μμ channels.

Expected and observed 95% CL upper limits on the signal strength for the Type-II 2HDM benchmark (mH = 379 GeV and m A = 172 GeV) as a function of tan(beta) and cos(beta-alpha). The limits are computed using the asymptotic CLs method, combining the ee and μμ channels.

P-value for H(A) -> ZA(H) -> bbbar production as a function of mA and mH. The limits are computed using the asymptotic CLs method, combining the ee and μμ channels.

J/psi nuclear modification in p+Au collisions as a function of centrality and rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au collisions as a function of centrality and rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al MinBias collisions as a function of pT. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Au MinBias collisions as a function of pT. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au MInBias collisions as a function of pT. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Au collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Au collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al most central collisions as a function of pT is directly compared with p+Au collisions (Data listed in Table16) at forward and backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in d+Au and p+Au most central collisions as a function of pT is directly compared at forward and backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column. Please see PRC 87, 034904 for d+Au data points.

J/psi nuclear modification in 3He+Au most central collisions as a function of pT is directly compared to p+Au collisions (Data listed in Table 16) at forward and backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Au collisions as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al collisions as a function of nuclear thickness (T_A). The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Au collisions as a function of nuclear thickness (T_A). The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au collisions as a function of nuclear thickness (T_A). The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi mean pT squared in p+Al collisions as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi mean pT squared in p+Au collisions as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi mean pT squared in 3He+Au collisions as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields as a function of rapidity in small systems. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Al collisions as a function of centrality and rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Au collisions as a function of centrality and apidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in 3He+Au collisions as a function of centrality and rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Al collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Al collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Au collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Au collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in 3He+Au collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in 3He+Au collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

Invariant $p_{T}$-differential spectra of $\Omega- + \bar{\Omega+}$ in p+p and p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in -1.8 < $y_{CM}$ < 0

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in -1.8 < $y_{CM}$ < 0

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in -1.8 < $y_{CM}$ < 0

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0 < $y_{CM}$ < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0 < $y_{CM}$ < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0 < $y_{CM}$ < 1.8

Invariant $p_{T}$-differential spectra of ${K_{0}}^{S}$ in p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

Invariant $p_{T}$-differential spectra of $\Lambda + \bar{\Lambda}$ in p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.3 < |$y_{CM}$| < 0.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.3 < |$y_{CM}$| < 0.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.3 < |$y_{CM}$| < 0.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.8 < |$y_{CM}$| < 1.3

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.8 < |$y_{CM}$| < 1.3

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.8 < |$y_{CM}$| < 1.3

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 1.3 < |$y_{CM}$| < 1.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 1.3 < |$y_{CM}$| < 1.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 1.3 < |$y_{CM}$| < 1.8