Measured and expected fiducial cross section for Z$\gamma$ production. Corresponding to Table 1 in the paper.

Data selection efficiency as a function of nuclear recoil energy

Data points used in analysis in log_10(S2)-S1 space

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson $c\tau$ for $m_{\tilde{\omega}}=5~\textrm{GeV}$ and $\mathcal{B}(\tilde{\omega}\to\mu\mu)=0.162$ in the vector portal model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson $c\tau$ for $m_{\tilde{\omega}}=15~\textrm{GeV}$ and $\mathcal{B}(\tilde{\omega}\to\mu\mu)=0.15$ in the vector portal model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson $c\tau$ for $m_{\tilde{\omega}}=20~\textrm{GeV}$ and $\mathcal{B}(\tilde{\omega}\to\mu\mu)=0.15$ in the vector portal model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\textrm{A}'$ $c\tau$ for $m_{\tilde{\pi}_{3}}=1~\textrm{GeV}$ and $m_{\textrm{A}'}=0.33~\textrm{GeV}$ , and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.458$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\textrm{A}'$ $c\tau$ for $m_{\tilde{\pi}_{3}}=4~\textrm{GeV}$ and $m_{\textrm{A}'}=0.4~\textrm{GeV}$ , and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.436$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\textrm{A}'$ $c\tau$ for $m_{\tilde{\pi}_{3}}=2~\textrm{GeV}$ and $m_{\textrm{A}'}=0.67~\textrm{GeV}$, and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.17$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\textrm{A}'$ $c\tau$ for $m_{\tilde{\pi}_{3}}=4~\textrm{GeV}$ and $m_{\textrm{A}'}=1.33~\textrm{GeV}$, and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.305$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\pi}_{3}$ meson $c\tau$ for $m_{\tilde{\pi}_{3}}=1~\textrm{GeV}$ and $m_{\textrm{A}'}=0.33~\textrm{GeV}$ , and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.458$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\pi}_{3}$ meson $c\tau$ for $m_{\tilde{\pi}_{3}}=2~\textrm{GeV}$ and $m_{\textrm{A}'}=0.4~\textrm{GeV}$ , and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.436$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\pi}_{3}$ meson $c\tau$ for $m_{\tilde{\pi}_{3}}=2~\textrm{GeV}$ and $m_{\textrm{A}'}=0.67~\textrm{GeV}$ , and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.17$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\pi}_{3}$ meson $c\tau$ for $m_{\tilde{\pi}_{3}}=4~\textrm{GeV}$ and $m_{\textrm{A}'}=1.33~\textrm{GeV}$, and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.305$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson mass for $c\tau=10~\textrm{mm}$ in the vector portal model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson mass for $c\tau=100~\textrm{mm}$ in the vector portal model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario A, for $\mathrm{A}'$ $c\tau=10~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=10m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario A, for $\mathrm{A}'$ $c\tau=10~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=3m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario A, for $\mathrm{A}'$ $c\tau=100~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=10m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario A, for $\mathrm{A}'$ $c\tau=100~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=3m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario B1, for $\tilde{\pi}_{3}$ $c\tau=10~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=10m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario B1, for $\tilde{\pi}_{3}$ $c\tau=10~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=3m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario B1, for $\tilde{\pi}_{3}$ $c\tau=100~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=10m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario B1, for $\tilde{\pi}_{3}$ $c\tau=100~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=3m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Observed upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson mass and $c\tau$ in the vector portal model. The phase space region that is left out is kinematically forbidden in the model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Observed upper limits at $95\%$ CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as functions of the $\textrm{A}'$ mass and $c\tau$ for $m_{\tilde{\pi}_{3}}=10m_{A'}$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Observed upper limits at $95\%$ CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as functions of the $\textrm{A}'$ mass and $c\tau$ for $m_{\tilde{\pi}_{3}}=3m_{A'}$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Observed upper limits at $95\%$ CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as functions of the $\textrm{A}'$ mass and $c\tau$ for $m_{\tilde{\pi}_{3}}=10m_{A'}$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Observed upper limits at $95\%$ CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as functions of the $\textrm{A}'$ mass and $c\tau$ for $m_{\tilde{\pi}_{3}}=3m_{A'}$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Cutflow efficiencies for representative signal model points in the vector portal.

Cutflow efficiencies for representative signal model points in Scenario A.

Cutflow efficiencies for representative signal model points in Scenario B1.

Ratio between the prompt $\Xi_c^+$ baryons in a multiplicity class to the multiplicity-integrated (INEL $>$ 0) class.

The prompt production yield ratios between $\Xi_c^0$ and $\rm {D}^0$ measured in the same multiplicity classes.

The prompt production yield ratios between $\Xi_c^+$ and $\rm {D}^0$ measured in the same multiplicity classes.

The prompt production yield ratios between $\Xi_c^0$ and $\Lambda_c^+$ measured in the same multiplicity classes.

The prompt production yield ratios between $\Xi_c^+$ and $\Lambda_c^+$ measured in the same multiplicity classes.

The branching-fraction ratio of BR($\Xi_c^0 \rightarrow \Xi^-e^+\nu_e$)/BR($\Xi_c^0 \rightarrow \Xi^-\pi^+$)

$p_\mathrm{{T}}$-differential production yield of $\overline{\Sigma}^{-}$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $\overline{\Sigma}^{+}$/ $\overline{\Sigma}^{-}$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $\overline{\Sigma}^{+}$/ $\overline{\Sigma}^{-}$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $(p+\overline{p})/(\pi^+ + \pi^-)$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{+}/(\pi^+ + \pi^-)$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{-}/(\pi^+ + \pi^-)$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{+}/(p + \overline{p})$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{-}/(p + \overline{p})$ in INEL pp collisions at $\sqrt{s}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $(p+\overline{p})/(\pi^+ + \pi^-)$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.5<y_\mathrm{CMS}<0$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{+}/(\pi^+ + \pi^-)$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{-}/(\pi^+ + \pi^-)$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $(\Lambda + \overline{\Lambda})/(\pi^+ + \pi^-)$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.5<y_\mathrm{CMS}<0$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{+}/(p + \overline{p})$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{-}/(p + \overline{p})$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $(\Lambda + \overline{\Lambda})/(p + \overline{p})$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.5<y_\mathrm{CMS}<0$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{+}/(\Lambda + \overline{\Lambda})$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Ratio of $p_\mathrm{{T}}$-differential production yields of $2\overline{\Sigma}^{-}/(\Lambda + \overline{\Lambda})$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

$p_\mathrm{{T}}$-differential nuclear modification factor $R_\mathrm{{pPb}}$ of $\overline{\Sigma}^{+}$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

$p_\mathrm{{T}}$-differential nuclear modification factor $R_\mathrm{{pPb}}$ of $\overline{\Sigma}^{-}$ in NSD p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y_\mathrm{CMS}|<0.5$.

Correlation matrix for the correlated systematic uncertainties of the forward rapidity $\eta$ meson cross section.

Correlation matrix for the correlated systematic uncertainties of the central rapidity $\eta$ meson cross section.

The invariant yields of $\phi$ and $\omega+\rho$ mesons as a function of $\langle N_{\rm part}\rangle$ in Au+Au collisions. The systematic uncertainties of type-A (uncorrelated) are combined with statistical uncertainties in quadrature and are labeled as stat. Type-B (correlated) systematic uncertainties are listed as sys.

$R_{AA}$ of $\phi$ meson as a function of $p_T$ for in Au+Au collisions at $\sqrt{s_{_{NN}}}=200$~GeV, compared with Cu+Au collisions at $1.2<|y|<2.2$. The systematic uncertainties of type-A (uncorrelated) are combined with statistical uncertainties in quadrature and are labeled as stat. Type-B (correlated) systematic uncertainties are listed as sys.

$R_{AA}$ of $\phi$ meson as a function of $\langle N_{\rm part}\rangle$ for $1.2<|y|<2.2$ in Au+Au collisions at $\sqrt{s_{_{NN}}}=200$~GeV. The systematic uncertainties of type-A (uncorrelated) are combined with statistical uncertainties in quadrature and are labeled as stat. Type-B (correlated) systematic uncertainties are listed as sys.

$\xi$ distributions for different jet $p_T$ bins.

$\xi$ distributions for different jet $p_T$ bins.

The $j_T/p_{T}^{\rm jet}$ distributions for different jet $p_T$ bins.

The $j_T/p_{T}^{\rm jet}$ distributions for different jet $p_T$ bins.

$r$ distributions for different jet $p_T$ bins.