$\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.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} p_T(D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(\Upsilon(1S))}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(\Upsilon(1S))}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(D^0)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{\pi}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \left| \Delta \phi \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{\pi}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \left| \Delta \phi \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \left| \Delta y \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \left| \Delta y \right| }$ for $2<y(\Upsilon(1S))<4.5$, for $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} p_T(\Upsilon(1S)D^0) }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} p_T(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(\Upsilon(1S)D^0) }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \mathcal{A}_T }$, where $\mathcal{A}_T = \frac{p_T(\Upsilon(1S))-p_T(D^0)}{p_T(\Upsilon(1S))-p_T(D^0)}$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \mathcal{A}_T }$, where $\mathcal{A}_T = \frac{p_T(\Upsilon(1S))-p_T(D^0)}{p_T(\Upsilon(1S))-p_T(D^0)}$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} m(\Upsilon(1S)D^0) }$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} m(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T$ [pb/(GeV/$c$)] for $2.0 < y < 4.5$.

Differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}y$ [pb] for $p_T < 30$ GeV.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(2S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(1S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(3S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(1S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(3S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(2S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

The ratio of differential cross-sections $( \mathrm{d} \sigma ( pp \to ( \Upsilon(iS) \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T ) / ( \mathrm{d} \sigma ( pp \to ( \Upsilon(jS) \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T )$ for $2.0 < y < 4.5$.

The ratio of differential cross-sections $( \mathrm{d} \sigma ( pp \to ( \Upsilon(iS) \to \mu^+ \mu^- ) X ) / \mathrm{d}y ) / ( \mathrm{d} \sigma ( pp \to ( \Upsilon(jS) \to \mu^+ \mu^- ) X ) / \mathrm{d}y )$ for $p_T < 30$ GeV/$c$.

Double-differential cross-section $\mathrm{d}^2 \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d}y$ [pb/(GeV/$c$)].

Double-differential cross-section $\mathrm{d}^2 \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d}y$ [pb/(GeV/$c$)].

Double-differential cross-section $\mathrm{d}^2 \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d}y$ [pb/(GeV/$c$)].

Differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T$ [pb/(GeV/$c$)] for $2.0 < y < 4.5$.

Differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}y$ [pb] for $p_T < 30$ GeV.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(2S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(1S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(3S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(1S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(3S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(2S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

The ratio of differential cross-sections $( \mathrm{d} \sigma ( pp \to ( \Upsilon(iS) \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T ) / ( \mathrm{d} \sigma ( pp \to ( \Upsilon(jS) \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T )$ for $2.0 < y < 4.5$.

The ratio of differential cross-sections $( \mathrm{d} \sigma ( pp \to ( \Upsilon(iS) \to \mu^+ \mu^- ) X ) / \mathrm{d}y ) / ( \mathrm{d} \sigma ( pp \to ( \Upsilon(jS) \to \mu^+ \mu^- ) X ) / \mathrm{d}y )$ for $p_T < 30$ GeV/$c$.

The ratio of differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T$ [pb/(GeV/$c$)] for $\sqrt{s}=8$ and $\sqrt{s}=7$ TeV in $2.0 < y < 4.5$.

The ratio of differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}y$ for $\sqrt{s}=8$ and $\sqrt{s}=7$ TeV and $p_T < 30$ GeV/$c$.

Fractions $R_{\Upsilon(1S)}^{\chi_{b}(mP)}$ in bins of $p_{\mathrm{T}}^{\Upsilon}$, measured for data collected at $\sqrt{s}= 8$ TeV.

Fractions $R_{\Upsilon(1S)}^{\chi_{b}(mP)}$ in bins of $p_{\mathrm{T}}^{\Upsilon}$, measured for data collected at $\sqrt{s}= 8$ TeV.

Fractions $R_{\Upsilon(1S)}^{\chi_{b}(mP)}$ in bins of $p_{\mathrm{T}}^{\Upsilon}$, measured for data collected at $\sqrt{s}= 8$ TeV.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 2.5-3.0. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 3.0-3.5. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 3.5-4.0. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 4.0-4.5. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The total integrated cross sections for UPSILON 1S 2S and 3S times their dimuon branching ratio in the rapidity range 2.0-4.5 and transverse momentum 0-15 GeV/c.

Differential production cross section for UPSILON 1S 2s and 3S mesons times their dimuon branching fraction as a function of transverse momentum ;.

Differential production cross section for UPSILON 1S 2s and 3S mesons times their dimuon branching fraction as a function of rapidity ;.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 2.0-2.5, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 2.5-3.0, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 3.0-3.5, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 3.5-4.0, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 4.0-4.5, assuming no polarisation.

Ratios of the cross-sections for UPSILON(2S)->MU+MU- and UPSILON(3S)->MU+MU- with respect toe UPSILON(1S)->MU+MU- as a function of PT in the rapidity range 2.0-4.0 ;.

Ratios of the cross-sections for UPSILON(2S)->MU+MU- and UPSILON(3S)->MU+MU- with respect to UPSILON(1S)->MU+MU- as a function of rapidity in the rapidity range 2.0-4.0 ;.