$K^-$ (a), invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\phi$ meson (b) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\Xi^-$ (c) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$K^-$ (a), invariant yields as a function of $m_T-m_0$ for various rapidity regions in 40--60\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\phi$ meson (b) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 40--60\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

Rapidity density distributions of $K^-$ (squares), $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

Rapidity density distributions of $\phi$ meson (circles) $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

Rapidity density distributions of $\Xi^-$ (diamonds) $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

$\phi/K^-$ (a) and $\phi/\Xi^-$ (b) ratio as a function of collision energy, $\sqrt{s_{\mathrm{NN}}}$. The solid black circles show the measurements presented here in 0-10\% centrality bin, while empty markers in black or grey are used for data from various other energies and/or collision systems. The grey solid line represents a THERMUS calculation based on the Grand Canonical Ensemble (GCE) while the dotted lines depict calculations based on the Canonical Ensemble (CE) with different parameters of strangeness correlation radius ($r_c$). The green dashed line, green shaded band and the solid red line show transport model calculations from the public versions $\mathrm{UrQMD}^{1}$, modified $\mathrm{UrQMD}^{2}$ and SMASH, respectively.

Transverse mass $m_T$ spectrum at midrapidity of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c. $m_0$ is the PDG mass of the $\phi$ meson.

Transverse mass $m_T$ spectrum at midrapidity of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c. $m_0$ is the PDG mass of the $\phi$ meson.

Dependence of the slope parameter $T$ of the transverse momentum spectrum on rapidity $y$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c.

Dependence of the slope parameter $T$ of the transverse momentum spectrum on rapidity $y$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c.

Dependence of the width $\sigma_y$ of the rapidity distributions on $p_T$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c.

Dependence of the width $\sigma_y$ of the rapidity distributions on $p_T$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c.

Rapidity spectrum of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c.

Rapidity spectrum of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c.

Rapidity spectrum of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 40 GeV/c.

Energy dependence of ratios of total yields of $\phi$ mesons to mean total yields of pions produced in in minimum bias p+p collisions at CERN SPS energies.

Energy dependence of double ratios of total yields of $\phi$ mesons to mean total yields of pions in central Pb+Pb collisions over minimum bias p+p collisions at CERN SPS energies.

Energy dependence of total yields of $\phi$ mesons produced in minimum bias p+p collisions at CERN SPS energies.

Energy dependence of midrapidity yields of $\phi$ mesons produced in minimum bias p+p collisions at CERN SPS energies.

Widths of rapidity distributions of $\phi$ mesons produced in minimum bias p+p collisions at CERN SPS energies, as a function of beam rapidity.