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

Single differential cross-section for UPSI(2S) times the dimuon branching fraction as a function of PT for the rapidity region 2.0-4.5 without normalisation to the bin sizes. The first uncertainty is statistical and the second systematic.

Single differential cross-section for UPSI(3S) times the dimuon branching fraction as a function of PT for the rapidity region 2.0-4.5 without normalisation to the bin sizes. The first uncertainty is statistical and the second systematic.

Single differential cross-section for UPSI(1S) times the dimuon branching fraction as a function of rapidity for PT < 15 GeV without normalisation to the bin sizes. The first uncertainty is statistical and the second systematic.

Single differential cross-section for UPSI(2S) times the dimuon branching fraction as a function of rapidity for PT < 15 GeV without normalisation to the bin sizes. The first uncertainty is statistical and the second systematic.

Single differential cross-section for UPSI(3S) times the dimuon branching fraction as a function of rapidity for PT < 15 GeV without normalisation to the bin sizes. The first uncertainty is statistical and the second systematic.

Ratios of UPSI(2S) to UPSI(1S) and UPSI(3S) to UPSI(1S) cross-sections times the dimuon branching fractions as a function of PT in the rapidity range 2.0-4.5. The first uncertainty is statistical and the second systematic.

Ratios of UPSI(2S) to UPSI(1S) and UPSI(3S) to UPSI(1S) cross-sections times the dimuon branching fractions as a function of rapidity for PT < 15 GeV. The first uncertainty is statistical and the second systematic.

The $p_{\rm T}$ bin width and corrected yield ratios R21 and R31 for $|y| < 0.6$. The statistical and systematic uncertainties in the corrected yield ratios are given as the percentage of the ratio.

The $p_{\rm T}$ bin width and corrected yield ratios R21 and R31 for $0.6 < |y| < 1.2$. The statistical and systematic uncertainties in the corrected yield ratios are given as the percentage of the ratio.

The $p_{\rm T}$ bin width and corrected yield ratios R21 and R31 for $|y| < 1.2$. The statistical and systematic uncertainties in the corrected yield ratios are given as the percentage of the ratio.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon$(1S) state and $|y| <0.6$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(1S) state and $0.6 < |y| < 1.2$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(1S) state and $|y| < 1.2$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(2S) state and $|y| < 0.6$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(2S) state and $0.6 < |y| < 1.2$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(2S) state and $|y| < 1.2$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(3S) state and $|y| < 0.6$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(3S) state and $0.6 < |y| < 1.2$.

The dimuon acceptance $\mathcal{A}$ calculated using the CMS measured polarization and its positive and negative uncertainties $\mathcal{A}(\sigma^{\pm})$, and the values of the acceptance assuming no polarization $\mathcal{A}$(unpol), transverse polarization $\mathcal{A}(T)$, and longitudinal polarization $\mathcal{A}(L)$ for the $\Upsilon $(3S) state and $|y| < 1.2$.

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 ;.

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