Measured $R_\text{AA}$ for the $\Upsilon(3S)$ state in the 0–90% centrality interval, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty "PP MB" represents the pp luminosity and PbPb $N_\text{MB}$ combined uncertainties, whereas the global uncertainty "PP 3S" corresponds to the uncertainty on the$\Upsilon(3S)$ pp yields.

Measured $R_\text{AA}$ for the $\Upsilon(2S)$ state as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval. The global uncertainty combines the uncertainties of $T_\text{AA}$ , pp luminosity, and PbPb $N_\text{MB}$.

Measured $R_\text{AA}$ for the $\Upsilon(3S)$ state as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval. The global uncertainty combines the uncertainties of $T_\text{AA}$ , pp luminosity, and PbPb $N_\text{MB}$.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ as functions of PbPb collision centrality, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty represents the combined systematic and statistical uncertainties from pp data.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ in the 0–90% centrality interval, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty represents the combined systematic and statistical uncertainties from pp data.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval.

Yields for $\Upsilon(2S)$ mesons in PbPb collisions in centrality $0--90\%$ and $|y| < 2.4$, corrected for acceptance and efficiency, and normalized by the nuclear thickness function $\langle T_{AA} \rangle$ and the number of minimum bias events $N_{MB}$. The numbers are in units of pb$/GeV/c$.

Yields for $\Upsilon(2S)$ mesons in PbPb collisions in $p_\text{T} < 30GeV/c$ and $|y| < 2.4$, corrected for acceptance and efficiency, and normalized by the nuclear thickness function $\langle T_{AA} \rangle$ and the number of minimum bias events $N_{MB}$. The numbers are in units of pb$/GeV/c$.

Yields for $\Upsilon(3S)$ mesons in PbPb collisions in centrality $0--90\%$ and $|y| < 2.4$, corrected for acceptance and efficiency, and normalized by the nuclear thickness function $\langle T_{AA} \rangle$ and the number of minimum bias events $N_{MB}$. The numbers are in units of pb$/GeV/c$.

Yields for $\Upsilon(3S)$ mesons in PbPb collisions in $p_\text{T} < 30GeV/c$ and $|y| < 2.4$, corrected for acceptance and efficiency, and normalized by the nuclear thickness function $\langle T_{AA} \rangle$ and the number of minimum bias events $N_{MB}$. The numbers are in units of pb$/GeV/c$.

Yields for $\Upsilon(1S)$ mesons in PbPb collisions in $p_\text{T} < 30GeV/c$ and $|y| < 2.4$, corrected for acceptance and efficiency, and normalized by the nuclear thickness function $\langle T_{AA} \rangle$ and the number of minimum bias events $N_{MB}$. The numbers are in units of pb$/GeV/c$. The values and global integrated luminosity uncertainty are quoted from data taken from CMS in 2018 LHC PbPb run.

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.

$R_{dAu}$ vs. $y$.

Comparison of our d+Au measurements to the pA measurements from E772. Ratio of $\Upsilon$ production in pA to pp scaled by mass number as a function of mass number. Shown are the 1S and 2S+3S $\Upsilon$ measurements from E772 and our 1S measurement.

Comparison of our d+Au measurements to the pA measurements from E772. Exponent $\alpha$ as a function of $x_{F}$.

Invariant mass distributions of electron pairs in the region $|y_{ee}| < 1.0$ for the centrality selections $30-60\%$.

Invariant mass distributions of electron pairs in the region $|y_{ee}| < 1.0$ for the centrality selections $10-30\%$.

Invariant mass distributions of electron pairs in the region $|y_{ee}| < 1.0$ for the centrality selections $0-10\%$.

Nuclear modification factor for $\Upsilon$(1S+2S+3S), in $|y| < 1.0$ in d+Au and Au+Au collisions as a function of $N_{part}$.

Nuclear modification factor for $\Upsilon$(1S+2S+3S), in $|y| < 0.5$, in d+Au and Au+Au collisions as a function of $N_{part}$.

Nuclear modification factor for $\Upsilon$(1S) in $|y| < 1.0$, in d+Au and Au+Au collisions as a function of $N_{part}$.

Nuclear modification factor for $\Upsilon$(1S) in $|y| < 0.5$, in d+Au and Au+Au collisions as a function of $N_{part}$.

Nuclear modification factor of quarkonium states as a function of binding energy as measured by STAR.

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 section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factors for Υ(1S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.3%) are displayed at unity as empty and filled red black boxes, respectively.

Nuclear modification factors for Υ(1S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.3%) are displayed at unity as empty and filled red black boxes, respectively.

Nuclear modification factors for Υ(2S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.9%) are displayed at unity as empty and filled black black boxes, respectively.

Nuclear modification factors for Υ(2S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.9%) are displayed at unity as empty and filled black black boxes, respectively.

Nuclear modification factors for Υ(1S), Υ(2S) and Υ(3S) meson production in PbPb collisions, integrated over centrality. For the Υ(3S), the upper limit derived on the nuclear modification factor is given.

Nuclear modification factors for Υ(1S), Υ(2S) and Υ(3S) meson production in PbPb collisions, integrated over centrality. For the Υ(3S), the upper limit derived on the nuclear modification factor is given.

The cross section of the process E+ E- --> DEUTBAR X.

The ratio of the cross sections of the processes E+ E- --> DEUTBAR X and E+ E- --> HADRONS.

Fiducial Upsilon(3S) production cross-section, where pT>4 GeV and |eta|<2.3 for both muons, as a function of Upsilon(3S) pT in the Upsilon(3S) rapidity (|y|) bins 0-1.2 and 1.2-2.25. The first uncertainty is statistical, the second is systematic.

Fiducial Upsilon(nS) production cross-section, where pT>4 GeV and |eta|<2.3 for both muons, as a function of Upsilon(nS) rapidity (|y|) in Upsilon(nS) PT ranging from 0 - 70 GeV. The first uncertainty is statistical, the second is systematic.

Inclusive Upsilon(1S) production cross-section as a function of Upsilon(1S) pT in the Upsilon(1S) rapidity (|y|) bins 0-1.2 and 1.2-2.25. The first uncertainty is statistical, the second is systematic and the third encapsulates any possible variation due to spin-alignment from the unpolarised central value.

Inclusive Upsilon(2S) production cross-section as a function of Upsilon(2S) pT in the Upsilon(2S) rapidity (|y|) bins 0-1.2 and 1.2-2.25. The first uncertainty is statistical, the second is systematic and the third encapsulates any possible variation due to spin-alignment from the unpolarised central value.

Inclusive Upsilon(3S) production cross-section as a function of Upsilon(3S) pT in the Upsilon(3S) rapidity (|y|) bins 0-1.2 and 1.2-2.25. The first uncertainty is statistical, the second is systematic and the third encapsulates any possible variation due to spin-alignment from the unpolarised central value.

Inclusive Upsilon(nS) production cross-section as a function of Upsilon(nS) rapidity (|y|) in Upsilon(nS) PT ranging from 0 - 70 GeV. The first uncertainty is statistical, the second is systematic and the third encapsulates any possible variation due to spin-alignment from the unpolarised central value.

Ratio of Upsilon(2S) to Upsilon(1S) inclusive cross-section as a function of Upsilon pT in the Upsilon rapidity (|y|) bins 0-1.2 and 1.2-2.25. The first uncertainty is statistical, the second is systematic.

Ratio of Upsilon(3S) to Upsilon(1S) inclusive cross-section as a function of Upsilon pT in the Upsilon rapidity (|y|) bins 0-1.2 and 1.2-2.25. The first uncertainty is statistical, the second is systematic.

Ratio of Upsilon(2S) and Upsilon(3S) to Upsilon(1S) inclusive cross-section as a function of Upsilon rapidity (|y|) in Upsilon pT 0 - 70 GeV. The first uncertainty is statistical, the second is systematic.

UPSI(3S) differential cross sections for |y| < 2.

UPSI(1S) differential cross sections for |y| < 1.

UPSI(2S) differential cross sections for |y| < 1.

UPSI(3S) differential cross sections for |y| < 1.

UPSI(1S) differential cross sections for 1 < |y| < 2.

UPSI(2S) differential cross sections for 1 < |y| < 2.

UPSI(3S) differential cross sections for 1 < |y| < 2.

UPSI(1S) differential cross sections as a function of rapidity.

Ratio of upsilon cross sectis.on.

UPSI differential cross sections for |y| < 2.

UPSI differential cross sections for |y| < 1.

UPSI differential cross sections for 1< |y| < 2.

Differential cross section per nucleon as a function of transverse momentum for UPSILON production on the P target.

Measurements of the longitudinally polarized fraction in UPSILON(1S) events. Also given is the measured value of ALPHA in the angular distribution 1+ALPHA*COS(THETA*)**2.

SIG*Br(UPSI --> MU+ MU-).

The first (asymmetric) systematic error represents the uncertainty in the decay angular distribution of the form (1 + a*cos(theta*)**2) where a = +1 and -1 respectively. The second is the quadratic sum of the other systematic errors.

The first (asymmetric) systematic error represents the uncertainty in the decay angular distribution of the form (1 + a*cos(theta*)**2) where a = +1 and -1 respectively. The second is the quadratic sum of the other systematic errors.

The first (asymmetric) systematic error represents the uncertainty in the decay angular distribution of the form (1 + a*cos(theta*)**2) where a = +1 and -1 respectively. The second is the quadratic sum of the other systematic errors.