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.

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.

$\Upsilon$(nS) polarization parameter $\tilde\lambda$ in the HX frame for $p_T$ of 10 to 15 GeV. The global uncertainties, independent of state and $N_{ch}$ bin, are also indicated.

$\Upsilon$(nS) polarization parameter $\lambda_\vartheta$ in the HX frame for $p_T$ of 15 to 35 GeV. The global uncertainties, independent of state and $N_{ch}$ bin, are also indicated.

$\Upsilon$(nS) polarization parameter $\lambda_\varphi$ in the HX frame for $p_T$ of 15 to 35 GeV. The global uncertainties, independent of state and $N_{ch}$ bin, are also indicated.

$\Upsilon$(nS) polarization parameter $\lambda_{\vartheta\varphi}$ in the HX frame for $p_T$ of 15 to 35 GeV. The global uncertainties, independent of state and $N_{ch}$ bin, are also indicated.

$\Upsilon$(nS) polarization parameter $\tilde\lambda$ in the HX frame for $p_T$ of 15 to 35 GeV. The global uncertainties, independent of state and $N_{ch}$ bin, are also indicated.

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.

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

Differential production cross sections of psi(2S) as a function of pT.

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

integrated production cross section of psi(2S).

psi(2S) over J/psi ratio as a function of pt.

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

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

Differential production cross sections of Upsilon(1S) as a function of pT.

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

integrated production cross section of Upsilon (1S).

integrated production cross section of Upsilon (2S).

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

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

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range < 2.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range < 0.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range < 0.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range < 0.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.4 to 0.8. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.4 to 0.8. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.4 to 0.8. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.8 to 1.2. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.8 to 1.2. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 0.8 to 1.2. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.2 to 1.6. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.2 to 1.6. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.2 to 1.6. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.6 to 2.0. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.6 to 2.0. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 1.6 to 2.0. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 2.0 to 2.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 2.0 to 2.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of PT for the |rapidity| range 2.0 to 2.4. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(1S) production cross sections (times di-muon branching ratio) as a function of rapidity for PT < 50 GeV/c. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(1S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(2S) production cross sections (times di-muon branching ratio) as a function of rapidity for PT < 50 GeV/c. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(2S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The fiducial and acceptance corrected UPSI(3S) production cross sections (times di-muon branching ratio) as a function of rapidity for PT < 50 GeV/c. Note these are integrated cross sections and the acceptance-corrected cross sections assume the UPSI(3S) are unpolarized with the variations due to the 4 extreme polarization scenarios shown in the last 4 columns. The fiducial cross sections do not need to make any assumptions on the polarizations scenarios. The luminosity uncertainty of 4% is not included in the systematic errors.

The ratios of the UPSI(3S) to UPSI(1S) fiducial and acceptance cross sections as a function of PT for the |rapidiy| range < 2.4.

The ratios of the UPSI(2S) to UPSI(1S) fiducial and acceptance cross sections as a function of PT for the |rapidiy| range < 2.4.

The ratios of the UPSI(3S) to UPSI(2S) fiducial and acceptance cross sections as a function of PT for the |rapidiy| range < 2.4.

The acceptance-corrected UPSI(1S) production cross section, integrated and differential in PT, times the respective dimuon branching fraction, integrated over the rapidity range |y| < 1.2. The first (sys) error is the quadratic sum of all systematic uncertainties excluding that associated with the polarization of the UPSI(1S) which is shown separately as the second (sys) value. The latter has been obtained by calculating the acceptance with the +68 and -68% CL as given by the measured UPSI(1S) polarization [see. PRL 110 (2013) 081802].

The acceptance-corrected UPSI(2S) production cross section, integrated and differential in PT, times the respective dimuon branching fraction, integrated over the rapidity range |y| < 1.2. The first (sys) error is the quadratic sum of all systematic uncertainties excluding that associated with the polarization of the UPSI(2S) which is shown separately as the second (sys) value. The latter has been obtained by calculating the acceptance with the +68 and -68% CL as given by the measured UPSI(2S) polarization [see. PRL 110 (2013) 081802].

The acceptance-corrected UPSI(3S) production cross section, integrated and differential in PT, times the respective dimuon branching fraction, integrated over the rapidity range |y| < 1.2. The first (sys) error is the quadratic sum of all systematic uncertainties excluding that associated with the polarization of the UPSI(3S) which is shown separately as the second (sys) value. The latter has been obtained by calculating the acceptance with the +68 and -68% CL as given by the measured UPSI(3S) polarization [see. PRL 110 (2013) 081802].

Distribution of Lambda-Phi in the CS frame for Y(1S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta-Phi in the CS frame for Y(1S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta-Phi in the CS frame for Y(1S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Tilde in the CS frame for Y(1S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Tilde in the CS frame for Y(1S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta in the HX frame for Y(1S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta in the HX frame for Y(1S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Phi in the HX frame for Y(1S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Phi in the HX frame for Y(1S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta-Phi in the HX frame for Y(1S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta-Phi in the HX frame for Y(1S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Tilde in the HX frame for Y(1S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Tilde in the HX frame for Y(1S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta in the PX frame for Y(1S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta in the PX frame for Y(1S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Phi in the PX frame for Y(1S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Phi in the PX frame for Y(1S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta-Phi in the PX frame for Y(1S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta-Phi in the PX frame for Y(1S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Tilde in the PX frame for Y(1S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Tilde in the PX frame for Y(1S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta in the CS frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta in the CS frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Phi in the CS frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Phi in the CS frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta-Phi in the CS frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta-Phi in the CS frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Tilde in the CS frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Tilde in the CS frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta in the HX frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta in the HX frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Phi in the HX frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Phi in the HX frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta-Phi in the HX frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta-Phi in the HX frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Tilde in the HX frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Tilde in the HX frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta in the PX frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta in the PX frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Phi in the PX frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Phi in the PX frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta-Phi in the PX frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta-Phi in the PX frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Tilde in the PX frame for Y(2S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Tilde in the PX frame for Y(2S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta in the CS frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta in the CS frame for Y(3S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Phi in the CS frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Phi in the CS frame for Y(3S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta-Phi in the CS frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta-Phi in the CS frame for Y(3S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Tilde in the CS frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Tilde in the CS frame for Y(3S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta in the HX frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta in the HX frame for Y(3S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Phi in the HX frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Phi in the HX frame for Y(3S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta-Phi in the HX frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta-Phi in the HX frame for Y(3S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Tilde in the HX frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Tilde in the HX frame for Y(3S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta in the PX frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta in the PX frame for Y(3S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Phi in the PX frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Phi in the PX frame for Y(3S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Theta-Phi in the PX frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Theta-Phi in the PX frame for Y(3S) production in the |y| range 0.6-1.2.

Distribution of Lambda-Tilde in the PX frame for Y(3S) production in the |y| range 0.0-0.6.

Distribution of Lambda-Tilde in the PX frame for Y(3S) production in the |y| range 0.6-1.2.

Estimated number of signal events for Y(1S) production.

Estimated number of signal events for Y(2S) production.

Estimated number of signal events for Y(3S) production.

Double differential cross section for UPSI(1S) production times the dimuon branching fraction as a function of PT for the rapidity region 3.0-3.5. The second systematic (sys) error is due to the unknown polarisation of the UPSI(1S).

Double differential cross section for UPSI(1S) production times the dimuon branching fraction as a function of PT for the rapidity region 3.5-4.0. The second systematic (sys) error is due to the unknown polarisation of the UPSI(1S).

Double differential cross section for UPSI(1S) production times the dimuon branching fraction as a function of PT for the rapidity region 4.0-4.5. The second systematic (sys) error is due to the unknown polarisation of the UPSI(1S).

Double differential cross section for UPSI(2S) production times the dimuon branching fraction as a function of PT for the rapidity region 2.0-2.5. The second systematic (sys) error is due to the unknown polarisation of the UPSI(2S).

Double differential cross section for UPSI(2S) production times the dimuon branching fraction as a function of PT for the rapidity region 2.5-3.0. The second systematic (sys) error is due to the unknown polarisation of the UPSI(2S).

Double differential cross section for UPSI(2S) production times the dimuon branching fraction as a function of PT for the rapidity region 3.0-3.5. The second systematic (sys) error is due to the unknown polarisation of the UPSI(2S).

Double differential cross section for UPSI(2S) production times the dimuon branching fraction as a function of PT for the rapidity region 3.5-4.0. The second systematic (sys) error is due to the unknown polarisation of the UPSI(2S).

Double differential cross section for UPSI(2S) production times the dimuon branching fraction as a function of PT for the rapidity region 4.0-4.5. The second systematic (sys) error is due to the unknown polarisation of the UPSI(2S).

Double differential cross section for UPSI(3S) production times the dimuon branching fraction as a function of PT for the rapidity region 2.0-2.5. The second systematic (sys) error is due to the unknown polarisation of the UPSI(3S).

Double differential cross section for UPSI(3S) production times the dimuon branching fraction as a function of PT for the rapidity region 2.5-3.0. The second systematic (sys) error is due to the unknown polarisation of the UPSI(3S).

Double differential cross section for UPSI(3S) production times the dimuon branching fraction as a function of PT for the rapidity region 3.0-3.5. The second systematic (sys) error is due to the unknown polarisation of the UPSI(3S).

Double differential cross section for UPSI(3S) production times the dimuon branching fraction as a function of PT for the rapidity region 3.5-4.0. The second systematic (sys) error is due to the unknown polarisation of the UPSI(3S).

Double differential cross section for UPSI(3S) production times the dimuon branching fraction as a function of PT for the rapidity region 4.0-4.5. The second systematic (sys) error is due to the unknown polarisation of the UPSI(3S).

Ratios of cross section with respect to the UPSI(1S) cross section as a function of PT in the rapidity range 2.0-4.4. The second systematic (sys) error is due to the unknown polarisation of the three states.