Mean $p_T$ values of the $\Upsilon(1$S$)$, $\Upsilon(2$S$)$, and $\Upsilon(3$S$)$ states with $p_T\,>\,0\,GeV$ and $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $0<p_T<5\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $5<p_T<7\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $7<p_T<9\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $9<p_T<11\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $11<p_T<15\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $15<p_T<20\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $20<p_T<50\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $0<p_T<5\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $5<p_T<7\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $7<p_T<9\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $9<p_T<11\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $11<p_T<15\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $15<p_T<20\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $20<p_T<50\,$GeV range

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "forward" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Forward tracks are those with momentum direction in $\Delta\phi\,<\,\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "transverse" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Transverse tracks are those with momentum direction in $\pi/3\,<\,\Delta\phi\,<\,2\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "backward" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Backward tracks are those with momentum direction in $\Delta\phi\,>\,2\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,0$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,1$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,2$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,>\,2$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.00\,<\,S_T\,<\,0.55$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.55\,<\,S_T\,<\,0.70$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.70\,<\,S_T\,<\,0.85$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.85\,<\,S_T\,<\,1.00$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Efficiency-corrected multiplicity bins used in the $\Upsilon(n$S$)$ ratio analysis and the corresponding mean number of charged particle tracks.

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

Double-differential cross section times the dimuon branching fraction of the Y(2S) meson for different ranges of pT in bins of |y| and for the full |y| < 1.2 range, for the unpolarized decay hypothesis. The global uncertainty in the integrated luminosity of 2.3% is not included in the systematic uncertainties.

Double-differential cross section times the dimuon branching fraction of the Y(3S) meson for different ranges of pT in bins of |y| and for the full |y| < 1.2 range, for the unpolarized decay hypothesis. The global uncertainty in the integrated luminosity of 2.3% is not included in the systematic uncertainties.

Multiplicative scaling factors to obtain the J/psi differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, +0.10, -1) from the unpolarized cross section measurements given in Table 1.

Multiplicative scaling factors to obtain the psi(2S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, +0.03, -1) from the unpolarized cross section measurements given in Table 2.

Multiplicative scaling factors to obtain the Y(1S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

Multiplicative scaling factors to obtain the Y(2S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

Multiplicative scaling factors to obtain the Y(3S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

Ratio of differential cross section times branching fractions of the prompt psi(2S) to J/psi mesons assuming unpolarized productions as a function of pT for |y| < 1.2.

Ratios of differential cross section times branching fractions of the prompt Y(2S) to Y(1S) and Y(3S) to Y(1S) mesons assuming unpolarized productions as a function of pT for |y| < 1.2.

Differential cross-sections multiplied by dimuon branching fractions for $\Upsilon$(1S), $\Upsilon$(2S) and $\Upsilon$(3S) states versus $p_T$ integrated over $y$ between 2.0 and 4.5.

Differential cross-sections multiplied by dimuon branching fractions for $\Upsilon$(1S), $\Upsilon$(2S) and $\Upsilon$(3S) states versus $y$ integrated over $p_T$ from 0 to 15 GeV/c.

Ratios of double-differential cross-sections times dimuon branching fractions for $\Upsilon$(2S) to $\Upsilon$(1S).

Ratios of double-differential cross-sections times dimuon branching fractions for $\Upsilon$(3S) to $\Upsilon$(1S).

Ratios of differential cross-sections times dimuon branching fractions versus $p_T$ over $y$ for $\Upsilon$(2S) to $\Upsilon$(1S) and $\Upsilon$(3S) to $\Upsilon$(1S).

Ratios of differential cross-sections times dimuon branching fractions versus $y$ over $p_T$ for $\Upsilon$(2S) to $\Upsilon$(1S) and $\Upsilon$(3S) to $\Upsilon$(1S).

The cross-section ratios between 13 TeV and 8 TeV in different bins of $p_T$ and $y$ for $\Upsilon$(1S).

The cross-section ratios between 13 TeV and 8 TeV in different bins of $p_T$ and $y$ for $\Upsilon$(2S).

The cross-section ratios between 13 TeV and 8 TeV in different bins of $p_T$ and $y$ for $\Upsilon$(3S).

Ratios of differential cross-sections between 13 TeV and 8 TeV measurements versus $p_T$ over $y$ for $\Upsilon$(1S), $\Upsilon$(2S) and $\Upsilon$(3S).

Ratios of differential cross-sections between 13 TeV and 8 TeV measurements versus $y$ over $p_T$ for $\Upsilon$(1S), $\Upsilon$(2S) and $\Upsilon$(3S).

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

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