$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 10-20%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 20-40%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 40-60%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 60-80%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 80-100%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (NSD).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 0-5%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 5-10%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 10-20%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 20-40%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 40-60%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 60-80%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 80-100%).

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of K$^{*0}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of K$^{*0}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\Xi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\Omega$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Integrated fiducial cross section for the DY process measured in the muon channel. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of dimuon invariant mass. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Integrated cross section for the DY process measured in the muon channel. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Forward-backward ratios for $15<m_{\mu\mu}<60$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Forward-backward ratios for $60<m_{\mu\mu}<120$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of dimuon invariant mass.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in rapidity in the centre-of-mass frame. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in $p_{\textrm{T}}$. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in $\phi^*$. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

$\Upsilon(1S)$ production cross-section in $p$Pb and Pb$p$, as a function of $y^*$. The uncertainty is the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ production cross-section in $p$Pb and Pb$p$, as a function of $p_{T}$. The uncertainty is the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ production cross-section in $p$Pb and Pb$p$, as a function of $y^*$. The uncertainty is the sum in quadrature of the statistical and systematic components.

Scaled $pp$ differential cross-section in $p_{T}$ at $\sqrt{s_{NN}}=8.16$ TeV. The first uncertainty is statistical, the second is systematic, which includes the systematic uncertainty from the $pp$ measurement and that estimated by changing the interpolation function.

Scaled $pp$ differential cross-section in $y$ at $\sqrt{s_{NN}}=8.16$ TeV. The first uncertainty is statistical, the second is systematic, which includes the systematic uncertainty from the $pp$ measurement and that estimated by changing the interpolation function.

$\Upsilon(1S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(1S)}$, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(1S)}$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV$/c$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(2S)}$, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(2S)}$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV$/c$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ forward-to-backward ratio, $R_{{\rm FB}}^{\Upsilon(1S)}$, as a function of $p_{T}$ integrated over $\vert y^* \vert$ in the range $2.5 < \vert y^* \vert < 4.0$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ forward-to-backward ratio, $R_{{\rm FB}}^{\Upsilon(1S)}$, as a function of $\vert y^* \vert$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ to $\Upsilon(1S)$ ratio, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ to $\Upsilon(1S)$ ratio, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ to non-prompt $J/\psi$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

Ratios of transverse-momentum distributions of charged particles in three intervals of multiplicities to the respective one for inclusive (INEL>0) collisions. The spectra were normalised by the integral prior to division.

Differential production cross sections of $\psi$(2S) as a function of $p_{\rm T}$.

Differential production cross sections of $\psi$(2S) as a function of rapidity.

integrated production cross section of $\psi$(2S).

Differential production cross sections of $\Upsilon$(1S) as a function of $p_{\rm T}$.

Differential production cross sections of $\Upsilon$(1S) as a function of rapidity.

integrated production cross section of $\Upsilon$(1S).

Differential production cross sections of $\Upsilon$(2S) as a function of rapidity.

integrated production cross section of Upsilon (2S).

integrated production cross section of Upsilon (3S).

$\psi$(2S) over J/$\psi$ ratio as a function of $p_{\rm T}$.

$\psi$(2S) over J/$\psi$ ratio as a function of rapidity.

integrated $\psi$(2S) over J/$\psi$ production cross section ratio.

integrated $\Upsilon$(2S) over $\Upsilon$(1S) production cross section ratio.

integrated $\Upsilon$(3S) over $\Upsilon$(1S) production cross section ratio.