Total $J/\psi$ efficiency as a function of cos$\theta$ for $2 < p_{T}^{J/\psi} < 3$ GeV/c

Total $J/\psi$ efficiency as a function of cos$\theta$ for $3 < p_{T}^{J/\psi} < 4$ GeV/c

Total $J/\psi$ efficiency as a function of cos$\theta$ for $4 < p_{T}^{J/\psi} < 6$ GeV/c

Different contributions to total $J/\psi$ efficiency as a function of cos$\theta$ for $2 < p_{T}^{J/\psi} < 3$ GeV/c

Different contributions to total $J/\psi$ efficiency as a function of cos$\theta$ for $3 < p_{T}^{J/\psi} < 4$ GeV/c

Different contributions to total $J/\psi$ efficiency as a function of cos$\theta$ for $4 < p_{T}^{J/\psi} < 6$ GeV/c

Corrected cos$\theta$ distribution for $2 < p_{T}^{J/\psi} < 3$ GeV/c

Corrected cos$\theta$ distribution for $3 < p_{T}^{J/\psi} < 4$ GeV/c

Corrected cos$\theta$ distribution for $4 < p_{T}^{J/\psi} < 6$ GeV/c

$J/\psi$ polarization parameter $\lambda_{\theta}$ as a function of $J/\psi$ $p_{T}$ for rapidity |y|<1

The invariant yield versus transverse momentum for |y| < 1 in 0-60% centrality in Au+Au collisions (solid circles). The results are compared to high-$p_T$ (3 < $p_T$ < 10 GeV/c) results from STAR [9] (solid squares) and PHENIX data [8] (open squares). Also shown is the yield in 0-60% centrality Cu+Cu collisions for low $p_T$ (solid triangles) and high $p_T$ [48] (open triangles). The models are described in the text [49, 50]. The $J/\psi$ cross section in p+p collisions is also shown at STAR (stars) and PHENIX [51] (diamonds), and the scale is indicated on the right axis.

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 0-60% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 0-20% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 20-40% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 40-60% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ invariant yield and nuclear modification factor as a function of centrality for |y| < 1 in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $J/\psi$ invariant yield $Bd^2N/(2\pi p_Tdydp_T) [(GeV/c)^{-2}]$ and nuclear modification factor as a function of transverse momentum for |y| < 1 in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $J/\psi$ invariant yield and nuclear modification factor for |y| < 1 in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Excitation function at cos(Theta_K+)cm = -0.5

Excitation function at cos(Theta_K+)cm = -0.4

Excitation function at cos(Theta_K+)cm = -0.3

Excitation function at cos(Theta_K+)cm = -0.2

Excitation function at cos(Theta_K+)cm = -0.1

Excitation function at cos(Theta_K+)cm = 0.0

Excitation function at cos(Theta_K+)cm = 0.1

Excitation function at cos(Theta_K+)cm = 0.2

Excitation function at cos(Theta_K+)cm = 0.3

Excitation function at cos(Theta_K+)cm = 0.4

Excitation function at cos(Theta_K+)cm = -0.8

Excitation function at cos(Theta_K+)cm = -0.7

Excitation function at cos(Theta_K+)cm = -0.6

Excitation function at cos(Theta_K+)cm = -0.5

Excitation function at cos(Theta_K+)cm = -0.4

Excitation function at cos(Theta_K+)cm = -0.3

Excitation function at cos(Theta_K+)cm = -0.2

Excitation function at cos(Theta_K+)cm = -0.1

Excitation function at cos(Theta_K+)cm = 0.0

Excitation function at cos(Theta_K+)cm = 0.1

Excitation function at cos(Theta_K+)cm = 0.2

Excitation function at cos(Theta_K+)cm = 0.3

Excitation function at cos(Theta_K+)cm = 0.4

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

The differential cross section of $Z$ boson production scaled by 1/$p_\mathrm{T}^{Z}$, (1/$p_\mathrm{T}^{Z}$) $d\sigma /dp_\mathrm{T}^{Z}$, for $-3<y^\mathrm{*}_{Z}<2$.

The differential cross section of $Z$ boson production scaled by 1/$p_\mathrm{T}^{Z}$, (1/$p_\mathrm{T}^{Z}$) $d\sigma /dp_\mathrm{T}^{Z}$, for $-2<y^\mathrm{*}_{Z}<0$ and $0<y^\mathrm{*}_{Z}<2.

Centrality bias corrected $Z$ boson yields per event for $-3<y^\mathrm{*}_{Z}<2$ scaled by by $\langle N_{coll}\rangle$. (To remove the centrality bias correction each value may be multiplied by the approriate correction value found in arXiv:1412.0976.).

The rapidity differential Z boson yields per event scaled by $\langle N_{coll}\rangle$ for three centrality ranges.

THERMINATOR 2 simulation correlation function for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system in the 0--10\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{-}}\oplus\overline{\Lambda}\mathrm{K^{+}}$ system in the 0--10\% centrality interval.

THERMINATOR 2 simulation correlation function for the $\Lambda\mathrm{K^{-}}\oplus\overline{\Lambda}\mathrm{K^{+}}$ system in the 0--10\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{0}_{S}}\oplus\overline{\Lambda}\mathrm{K^{0}_{S}}$ system in the 0--10\% centrality interval.

THERMINATOR 2 simulation correlation function for the $\Lambda\mathrm{K^{0}_{S}}\oplus\overline{\Lambda}\mathrm{K^{0}_{S}}$ system in the 0--10\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system in the 10--30\% centrality interval.

THERMINATOR 2 simulation correlation function for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system in the 10--30\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{-}}\oplus\overline{\Lambda}\mathrm{K^{+}}$ system in the 10--30\% centrality interval.

THERMINATOR 2 simulation correlation function for the $\Lambda\mathrm{K^{-}}\oplus\overline{\Lambda}\mathrm{K^{+}}$ system in the 10--30\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{0}_{S}}\oplus\overline{\Lambda}\mathrm{K^{0}_{S}}$ system in the 10--30\% centrality interval.

THERMINATOR 2 simulation correlation function for the $\Lambda\mathrm{K^{0}_{S}}\oplus\overline{\Lambda}\mathrm{K^{0}_{S}}$ system in the 10--30\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system in the 30--50\% centrality interval.

THERMINATOR 2 simulation correlation function for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system in the 30--50\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{-}}\oplus\overline{\Lambda}\mathrm{K^{+}}$ system in the 30--50\% centrality interval.

THERMINATOR 2 simulation correlation function for the $\Lambda\mathrm{K^{-}}\oplus\overline{\Lambda}\mathrm{K^{+}}$ system in the 30--50\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{0}_{S}}\oplus\overline{\Lambda}\mathrm{K^{0}_{S}}$ system in the 30--50\% centrality interval.

THERMINATOR 2 simulation correlation function for the $\Lambda\mathrm{K^{0}_{S}}\oplus\overline{\Lambda}\mathrm{K^{0}_{S}}$ system in the 30--50\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system in the 0--10\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{-}}\oplus\overline{\Lambda}\mathrm{K^{+}}$ system in the 0--10\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{0}_{S}}\oplus\overline{\Lambda}\mathrm{K^{0}_{S}}$ system in the 0--10\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system in the 10--30\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{-}}\oplus\overline{\Lambda}\mathrm{K^{+}}$ system in the 10--30\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{0}_{S}}\oplus\overline{\Lambda}\mathrm{K^{0}_{S}}$ system in the 10--30\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system in the 30--50\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{-}}\oplus\overline{\Lambda}\mathrm{K^{+}}$ system in the 30--50\% centrality interval.

Measured correlation function for the $\Lambda\mathrm{K^{0}_{S}}\oplus\overline{\Lambda}\mathrm{K^{0}_{S}}$ system in the 30--50\% centrality interval.

Extracted imaginary and real components of the complex scattering length, $\Im f_{0}$ and $\Re f_{0}$, for the $\Lambda\mathrm{K^{0}_{S}}$ system.

Extracted imaginary and real components of the complex scattering length, $\Im f_{0}$ and $\Re f_{0}$, for the $\Lambda\mathrm{K^{-}}$ system.

Extracted imaginary and real components of the complex scattering length, $\Im f_{0}$ and $\Re f_{0}$, for the $\Lambda\mathrm{K^{+}}$ system.

Extracted effective range of the interaction, $d_{0}$, for the $\Lambda\mathrm{K^{0}_{S}}$ system.

Extracted effective range of the interaction, $d_{0}$, for the $\Lambda\mathrm{K^{-}}$ system.

Extracted effective range of the interaction, $d_{0}$, for the $\Lambda\mathrm{K^{+}}$ system.

The extracted $\lambda_{\mathrm{Fit}}$ and radius parameter for the $\Lambda$K system in the 0--10\% centrality interval.

The extracted $\lambda_{\mathrm{Fit}}$ and radius parameter for the $\Lambda$K system in the 10--30\% centrality interval.

The extracted $\lambda_{\mathrm{Fit}}$ and radius parameter for the $\Lambda$K system in the 30--50\% centrality interval.

Extracted fit $R_{\mathrm{inv}}$ parameters as a function of pair transverse mass ($m_{\mathrm{T}}$) for the $\Lambda$K system in the 0--10\% centrality interval.

Extracted fit $R_{\mathrm{inv}}$ parameters as a function of pair transverse mass ($m_{\mathrm{T}}$) for the $\Lambda$K system in the 10--30\% centrality interval.

Extracted fit $R_{\mathrm{inv}}$ parameters as a function of pair transverse mass ($m_{\mathrm{T}}$) for the $\Lambda$K system in the 30--50\% centrality interval.

Spherical harmonics component $C_{00}$ of the $\Lambda\mathrm{K^{+}}$ correlation function for the 0--10\% centrality interval.

Spherical harmonics component $\Re C_{11}$ of the $\Lambda\mathrm{K^{+}}$ correlation function for the 0--10\% centrality interval.

Correlation functions for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system built using the normal mixed-event reference distribution for the 0--10\% centrality interval.

Correlation functions for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system built using the Stavinskiy method for the 0--10\% centrality interval.

Correlation functions for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system built using the normal mixed-event reference distribution for the 10--30\% centrality interval.

Correlation functions for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system built using the Stavinskiy method for the 10--30\% centrality interval.

Correlation functions for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system built using the normal mixed-event reference distribution for the 30--50\% centrality interval.

Correlation functions for the $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ system built using the Stavinskiy method for the 30--50\% centrality interval.

One-dimensional $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ correlation function for the experimental data in the 0--10\% centrality interval.

One-dimensional $\Lambda\mathrm{K^{+}}\oplus\overline{\Lambda}\mathrm{K^{-}}$ correlation function for the THERMINATOR 2 simulation with impact parameter b = 2 fm.

$\Re C_{11}$ component of a spherical harmonic decomposition for the experimental data in the 0--10\% centrality interval.

$\Re C_{11}$ component of a spherical harmonic decomposition for the THERMINATOR 2 simulation with impact parameter b = 2 fm.

The pair source distribution from the THERMINATOR 2 simulation in the out direction.

The pair source distribution from the THERMINATOR 2 simulation in the side direction.

The pair source distribution from the THERMINATOR 2 simulation in the long direction.

The temporal characteristics of the pair source distribution from the THERMINATOR 2 simulation.

Probing the effect of varying the source shift in the outward direction, $\mu_{\mathrm{out}}$, within the THERMINATOR 2 framework. To achieve this, particle pairs are formed from the simulation, but with altered spatial characteristics achieved by drawing the out, side, and long components from predetermined Gaussian distributions. The sources in all three directions are Gaussians of width 5 fm. The distributions used for the side and long direction are centered at the origin, while the shift in the outward direction, $\mu_{\mathrm{out}}$, is varied. THERMINATOR 2 simulated correlation functions with $\mu_{\mathrm{out}} = 0$ fm.

Probing the effect of varying the source shift in the outward direction, $\mu_{\mathrm{out}}$, within the THERMINATOR 2 framework. To achieve this, particle pairs are formed from the simulation, but with altered spatial characteristics achieved by drawing the out, side, and long components from predetermined Gaussian distributions. The sources in all three directions are Gaussians of width 5 fm. The distributions used for the side and long direction are centered at the origin, while the shift in the outward direction, $\mu_{\mathrm{out}}$, is varied. THERMINATOR 2 simulated correlation functions with $\mu_{\mathrm{out}} = 1$ fm.

Probing the effect of varying the source shift in the outward direction, $\mu_{\mathrm{out}}$, within the THERMINATOR 2 framework. To achieve this, particle pairs are formed from the simulation, but with altered spatial characteristics achieved by drawing the out, side, and long components from predetermined Gaussian distributions. The sources in all three directions are Gaussians of width 5 fm. The distributions used for the side and long direction are centered at the origin, while the shift in the outward direction, $\mu_{\mathrm{out}}$, is varied. THERMINATOR 2 simulated correlation functions with $\mu_{\mathrm{out}} = 3$ fm.

Probing the effect of varying the source shift in the outward direction, $\mu_{\mathrm{out}}$, within the THERMINATOR 2 framework. To achieve this, particle pairs are formed from the simulation, but with altered spatial characteristics achieved by drawing the out, side, and long components from predetermined Gaussian distributions. The sources in all three directions are Gaussians of width 5 fm. The distributions used for the side and long direction are centered at the origin, while the shift in the outward direction, $\mu_{\mathrm{out}}$, is varied. THERMINATOR 2 simulated correlation functions with $\mu_{\mathrm{out}} = 6$ fm.

Spherical harmonics component $C_{00}$ of the $\Lambda\mathrm{K^{+}}$ correlation function for the 0--10\% centrality interval.

Spherical harmonics component $C_{00}$ of the $\Lambda\mathrm{K^{+}}$ correlation function from the THERMINATOR 2 simulation implemented with a shift $\mu_{\mathrm{out}} = 0$ fm in the outward direction.

Spherical harmonics component $C_{00}$ of the $\Lambda\mathrm{K^{+}}$ correlation function from the THERMINATOR 2 simulation implemented with a shift $\mu_{\mathrm{out}} = 1$ fm in the outward direction.

Spherical harmonics component $C_{00}$ of the $\Lambda\mathrm{K^{+}}$ correlation function from the THERMINATOR 2 simulation implemented with a shift $\mu_{\mathrm{out}} = 3$ fm in the outward direction.

Spherical harmonics component $C_{00}$ of the $\Lambda\mathrm{K^{+}}$ correlation function from the THERMINATOR 2 simulation implemented with a shift $\mu_{\mathrm{out}} = 6$ fm in the outward direction.

Spherical harmonics component $\Re C_{11}$ of the $\Lambda\mathrm{K^{+}}$ correlation function for the 0--10\% centrality interval.

Spherical harmonics component $\Re C_{11}$ of the $\Lambda\mathrm{K^{+}}$ correlation function from the THERMINATOR 2 simulation implemented with a shift $\mu_{\mathrm{out}} = 0$ fm in the outward direction.

Spherical harmonics component $\Re C_{11}$ of the $\Lambda\mathrm{K^{+}}$ correlation function from the THERMINATOR 2 simulation implemented with a shift $\mu_{\mathrm{out}} = 1$ fm in the outward direction.

Spherical harmonics component $\Re C_{11}$ of the $\Lambda\mathrm{K^{+}}$ correlation function from the THERMINATOR 2 simulation implemented with a shift $\mu_{\mathrm{out}} = 3$ fm in the outward direction.

Spherical harmonics component $\Re C_{11}$ of the $\Lambda\mathrm{K^{+}}$ correlation function from the THERMINATOR 2 simulation implemented with a shift $\mu_{\mathrm{out}} = 6$ fm in the outward direction.