The production of $\Upsilon(nS)$ mesons ($n=1,2,3$) in $p$Pb and Pb$p$ collisions at a centreofmass energy per nucleon pair $\sqrt{s_{NN}}=8.16$ TeV is measured by the LHCb experiment, using a data sample corresponding to an integrated luminosity of 31.8 nb$^{1}$. The $\Upsilon(nS)$ mesons are reconstructed through their decays into two oppositesign muons. The measurements comprise the differential production crosssections of the $\Upsilon(1S)$ and $\Upsilon(2S)$ states, their forwardtobackward ratios and nuclear modification factors, performed as a function of the transverse momentum \pt and rapidity in the nucleonnucleon centreofmass frame $y^*$ of the $\Upsilon(nS)$ states, in the kinematic range $p_{\rm{T}}<25$ GeV/$c$ and $1.5<y^*<4.0$ ($5.0<y^*<2.5$) for $p$Pb (Pb$p$) collisions. In addition, production crosssections for $\Upsilon(3S)$ are measured integrated over phase space and the production ratios between all three $\Upsilon(nS)$ states are determined. The measurements are compared to theoretical predictions and suppressions for quarkonium in $p$Pb collisions are observed.
Invariantmass distribution of $\mu^+\mu^$ pairs from the (left) $p$Pb and (right) Pb$p$ samples after the trigger and off\line selections. 
Fig1a.pdf [32 KiB] HiDef png [367 KiB] Thumbnail [292 KiB] *.C file 

Fig1b.pdf [33 KiB] HiDef png [371 KiB] Thumbnail [296 KiB] *.C file 

Total efficiency $\varepsilon_{\mathrm{tot}} $ of the \Y1S meson as a function of its $p_{\mathrm{ T}}$ in different $y^*$ bins in (left) $p$Pb and (right) Pb$p$ collisions. The horizontal locations of the markers are roughly the centroids of the bins, with offsets from centre to aid in readability. 
Fig2a.pdf [17 KiB] HiDef png [248 KiB] Thumbnail [231 KiB] *.C file 

Fig2b.pdf [17 KiB] HiDef png [248 KiB] Thumbnail [234 KiB] *.C file 

Doubledifferential crosssection for the \Y1S meson as a function of $p_{\mathrm{ T}}$ for different values of $y^*$ for the (left) forward $p$Pb and (right) backward Pb$p$ samples. The uncertainties are the sums in quadrature of the statistical and systematic components. The horizontal locations of the markers are roughly the centroids of the bins, with offsets from centre to aid in readability. 
Fig3a.pdf [17 KiB] HiDef png [248 KiB] Thumbnail [234 KiB] *.C file 

Fig3b.pdf [17 KiB] HiDef png [248 KiB] Thumbnail [241 KiB] *.C file 

Crosssection of (left) \Y1S and (right) \Y2S production as a function of $y^*$ integrated over $p_{\mathrm{ T}}$ for the backward (negative $y^*$) and forward (positive $y^*$) samples, compared to the crosssection measured in $pp$, interpolated to $\sqrt{s_{\scriptscriptstyle\rm NN}} =8.16\mathrm{ Te V} $. In this and subsequent figures, the uncertainties shown are the sums in quadrature of the statistical and systematic components. 
Fig4a.pdf [14 KiB] HiDef png [156 KiB] Thumbnail [160 KiB] *.C file 

Fig4b.pdf [15 KiB] HiDef png [163 KiB] Thumbnail [172 KiB] *.C file 

Crosssection of (top) \Y1S and (bottom) \Y2S production as a function of $p_{\mathrm{ T}}$ integrated over $y^*$ for the (left) forward and (right) backward samples compared to the crosssection measured in $pp$, interpolated to $\sqrt{s_{\scriptscriptstyle\rm NN}} =8.16\mathrm{ Te V} $. 
Fig5a.pdf [14 KiB] HiDef png [166 KiB] Thumbnail [174 KiB] *.C file 

Fig5b.pdf [14 KiB] HiDef png [166 KiB] Thumbnail [174 KiB] *.C file 

Fig5c.pdf [14 KiB] HiDef png [164 KiB] Thumbnail [170 KiB] *.C file 

Fig5d.pdf [14 KiB] HiDef png [165 KiB] Thumbnail [172 KiB] *.C file 

Nuclear modification factors of the (left) \Y1S and (right) \Y2S mesons as a function of $y^*$ integrated over $p_{\mathrm{ T}}$ for the forward and backward samples. The bands correspond to the theoretical predictions for the nCTEQ15 and EPPS16 nPDFs sets, and the comovers model as reported in the text. 
Fig6a.pdf [20 KiB] HiDef png [1 MiB] Thumbnail [520 KiB] *.C file 

Fig6b.pdf [20 KiB] HiDef png [1 MiB] Thumbnail [562 KiB] *.C file 

Nuclear modification factors of the (top) \Y1S and (bottom) \Y2S mesons as a function of $p_{\mathrm{ T}}$ integrated over $y^*$ for the (left) forward and (right) backward samples. The bands correspond to the theoretical predictions for the nCTEQ15 and EPPS16 nPDFs sets as reported in the text. 
Fig7a.pdf [15 KiB] HiDef png [1002 KiB] Thumbnail [437 KiB] *.C file 

Fig7b.pdf [16 KiB] HiDef png [673 KiB] Thumbnail [336 KiB] *.C file 

Fig7c.pdf [15 KiB] HiDef png [999 KiB] Thumbnail [441 KiB] *.C file 

Fig7d.pdf [16 KiB] HiDef png [675 KiB] Thumbnail [337 KiB] *.C file 

Forwardbackward ratio for the \Y1S as a function of (left) $p_{\mathrm{ T}}$ integrated over $y^*$ and (right) as a function of $y^*$ integrated over $p_{\mathrm{ T}}$ . The bands correspond to the theoretical calculations for the nCTEQ15 and EPPS16 nPDFs sets as reported in the text. 
Fig8a.pdf [15 KiB] HiDef png [961 KiB] Thumbnail [430 KiB] *.C file 

Fig8b.pdf [14 KiB] HiDef png [1 MiB] Thumbnail [437 KiB] *.C file 

Forwardbackward ratio for the \Y2S compared with theoretical calculations for the nCTEQ15 and EPPS16 nPDFs sets as reported in the text. 
Fig9.pdf [14 KiB] HiDef png [1 MiB] Thumbnail [438 KiB] *.C file 

Ratios between \Y2S and \Y1S crosssections as a function of (top) $p_{\mathrm{ T}}$ integrated over $y^*$, and as function of (bottom) $y^*$ integrated over $p_{\mathrm{ T}}$ , for $p$Pb and Pb$p$ collisions. The bands correspond to the theoretical predictions for the nCTEQ15 and EPPS16 nPDFs sets as reported in the text. 
Fig10a.pdf [16 KiB] HiDef png [1 MiB] Thumbnail [498 KiB] *.C file 

Fig10b.pdf [16 KiB] HiDef png [725 KiB] Thumbnail [343 KiB] *.C file 

Fig10c.pdf [15 KiB] HiDef png [1 MiB] Thumbnail [445 KiB] *.C file 

Double ratios for (left) \Y2S and (right) \Y3S. The bands correspond to the theoretical prediction for the comovers model as reported in the text. 
Fig11a.pdf [16 KiB] HiDef png [339 KiB] Thumbnail [231 KiB] *.C file 

Fig11b.pdf [16 KiB] HiDef png [299 KiB] Thumbnail [209 KiB] *.C file 

Ratio of \Y1S to nonprompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu}$ crosssections as a function of $y^*$ integrated over $p_{\mathrm{ T}}$ , for $p$Pb and Pb$p$ collisions. 
Fig12.pdf [15 KiB] HiDef png [195 KiB] Thumbnail [201 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2018035.gif Thumbnail 
Yields of $\Y1S$, $\Y2S$, $\Y3S$ mesons in $p$Pb and Pb$p$ samples as given by the fit. The uncertainties are statistical only. 
Table_1.pdf [43 KiB] HiDef png [38 KiB] Thumbnail [20 KiB] tex code 

Systematic uncertainties (in percent) on the crosssection measurements. The ranges indicate the minimum and maximum values in different bins, among all \Y nS states. 
Table_2.pdf [26 KiB] HiDef png [126 KiB] Thumbnail [61 KiB] tex code 

Ratio $R(\Y nS)$ in $pp$, $p$Pb, and Pb$p$ samples. The uncertainties are combinations of statistical and systematical components. 
Table_3.pdf [36 KiB] HiDef png [73 KiB] Thumbnail [37 KiB] tex code 

\Y1S production crosssection in $p$Pb, as a function of $ p_{\mathrm{ T}} $ and $y^*$. 
Table_4.pdf [44 KiB] HiDef png [428 KiB] Thumbnail [220 KiB] tex code 

\Y1S production crosssection in Pb$p$, as a function of $ p_{\mathrm{ T}} $ and $y^*$. 
Table_5.pdf [44 KiB] HiDef png [353 KiB] Thumbnail [190 KiB] tex code 

\Y1S production crosssection in $p$Pb and Pb$p$, as a function of $ p_{\mathrm{ T}} $. 
Table_6.pdf [36 KiB] HiDef png [73 KiB] Thumbnail [38 KiB] tex code 

\Y1S production crosssection in $p$Pb and Pb$p$, as a function of $y^*$. 
Table_7.pdf [35 KiB] HiDef png [237 KiB] Thumbnail [104 KiB] tex code 

\Y2S production crosssection in $p$Pb and Pb$p$, as a function of $ p_{\mathrm{ T}} $. 
Table_8.pdf [36 KiB] HiDef png [68 KiB] Thumbnail [37 KiB] tex code 

\Y2S production crosssection in $p$Pb, as a function of $y^*$. 
Table_9.pdf [35 KiB] HiDef png [256 KiB] Thumbnail [109 KiB] tex code 

Scaled $pp$ differential crosssection in $p_{\mathrm{ T}}$ at $\sqrt{s_{\scriptscriptstyle\rm 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. 
Table_10.pdf [36 KiB] HiDef png [77 KiB] Thumbnail [42 KiB] tex code 

Scaled $pp$ differential crosssection in $y$ at $\sqrt{s_{\scriptscriptstyle\rm 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. 
Table_11.pdf [28 KiB] HiDef png [98 KiB] Thumbnail [52 KiB] tex code 

$\Upsilon(1S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(1S)}$, in $p$Pb and Pb$p$ as a function of $ p_{\mathrm{ T}} $ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $5.0 < y^* < 2.5$ for Pb$p$. 
Table_12.pdf [44 KiB] HiDef png [124 KiB] Thumbnail [59 KiB] tex code 

$\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_{\mathrm{ T}} $ in the range $0 < p_{\mathrm{ T}} < 25$ $ {\mathrm{ Ge V /}c} $ . 
Table_13.pdf [49 KiB] HiDef png [202 KiB] Thumbnail [86 KiB] tex code 

$\Upsilon(2S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(2S)}$, in $p$Pb and Pb$p$ as a function of $ p_{\mathrm{ T}} $ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $5.0 < y^* < 2.5$ for Pb$p$. 
Table_14.pdf [44 KiB] HiDef png [128 KiB] Thumbnail [60 KiB] tex code 

$\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_{\mathrm{ T}} $ in the range $0 < p_{\mathrm{ T}} < 25$ $ {\mathrm{ Ge V /}c} $ . 
Table_15.pdf [49 KiB] HiDef png [198 KiB] Thumbnail [86 KiB] tex code 

$\Upsilon(1S)$ forwardtobackward ratio, $R_{{\rm FB}}^{\Upsilon(1S)}$, as a function of $ p_{\mathrm{ T}} $ integrated over $\vert y^* \vert$ in the range $2.5 < \vert y^* \vert < 4.0$. 
Table_16.pdf [43 KiB] HiDef png [199 KiB] Thumbnail [83 KiB] tex code 

$\Upsilon(1S)$ forwardtobackward ratio, $R_{{\rm FB}}^{\Upsilon(1S)}$, as a function of $\vert y^* \vert$ integrated over $ p_{\mathrm{ T}} $ in the range $0 < p_{\mathrm{ T}} < 25$ $ {\mathrm{ Ge V /}c} $ . 
Table_17.pdf [49 KiB] HiDef png [99 KiB] Thumbnail [41 KiB] tex code 

$\Upsilon(2S)$ forwardtobackward ratio, $R_{{\rm FB}}^{\Upsilon(2S)}$, as a function of $ p_{\mathrm{ T}} $ integrated over $\vert y^* \vert$ in the range $2.5 < \vert y^* \vert < 4.0$. 
Table_18.pdf [36 KiB] HiDef png [128 KiB] Thumbnail [63 KiB] tex code 

$\Upsilon(2S)$ forwardtobackward ratio, $R_{{\rm FB}}^{\Upsilon(2S)}$, as a function of $\vert y^* \vert$ integrated over $ p_{\mathrm{ T}} $ in the range $0 < p_{\mathrm{ T}} < 25$ $ {\mathrm{ Ge V /}c} $ . 
Table_19.pdf [34 KiB] HiDef png [247 KiB] Thumbnail [106 KiB] tex code 

$\Upsilon(2S)$ to $\Upsilon(1S)$ ratio, $R(\Y2S)$, in $p$Pb and Pb$p$ as a function of $ p_{\mathrm{ T}} $ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $5.0 < y^* < 2.5$ for Pb$p$. 
Table_18.pdf [36 KiB] HiDef png [128 KiB] Thumbnail [63 KiB] tex code 

$\Upsilon(2S)$ to $\Upsilon(1S)$ ratio, $R(\Y2S)$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $ p_{\mathrm{ T}} $ in the range $0 < p_{\mathrm{ T}} < 25$ $ {\mathrm{ Ge V /}c} $ . 
Table_19.pdf [34 KiB] HiDef png [247 KiB] Thumbnail [106 KiB] tex code 

$\Upsilon(1S)$ to nonprompt $J/\psi$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $ p_{\mathrm{ T}} $ in the range $0 < p_{\mathrm{ T}} < 25$ $ {\mathrm{ Ge V /}c} $ . The quoted uncertainties are the sum in quadrature of the statistical and systematic components. 
Table_20.pdf [36 KiB] HiDef png [178 KiB] Thumbnail [64 KiB] tex code 
Created on 19 October 2019.