The production of $B^+$, $B^0$ and $\Lambda_b^0$ hadrons is studied in protonlead collisions at a centreofmass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}=8.16 {\rm TeV}$ recorded with the LHCb detector at the LHC. The measurement uses a dataset corresponding to an integrated luminosity of $12.2\pm0.3 \mathrm{nb}^{1}$ for the case where the proton beam is projected into the LHCb detector (corresponding to measuring hadron production at positive rapidity) and $18.6\pm0.5 \mathrm{nb}^{1}$ for the lead beam projected into the LHCb detector (corresponding to measuring hadron production at negative rapidity). Nuclear effects are probed through doubledifferential crosssections, forwardtobackward crosssection ratios and nuclear modification factors of the beauty hadrons. The doubledifferential crosssections are measured as a function of the beautyhadron transverse momentum and rapidity in the nucleonnucleon centreofmass frame. Forwardtobackward crosssection ratios and nuclear modification factors indicate a significant nuclear suppression at positive rapidity. The ratio of $\Lambda_b^0$ over $B^0$ production crosssections is reported and is consistent with the corresponding measurement in $pp$ collisions.
Invariant mass distribution of $ B ^+ $ candidates reconstructed in the $ B ^+ \rightarrow \overline{ D } {}^0 \pi ^+ $ decay for (left) $ p\mkern 1mu\mathrm{Pb}$ and (right) $\mathrm{Pb}\mkern 1.5mu p$ collisions, with the fit result superimposed. The solid blue line, the solid green line, the crossshaded area, the brown shaded area and the red shaded area represent the total fit, the signal component, the partially reconstructed background, the combinatorial background and $ B ^+ \rightarrow \overline{ D } {}^0 K ^+ $ decays, respectively. 
Fig1a.pdf [24 KiB] HiDef png [571 KiB] Thumbnail [333 KiB] *.C file 

Fig1b.pdf [24 KiB] HiDef png [581 KiB] Thumbnail [335 KiB] *.C file 

Invariant mass distribution of $ B ^+ $ candidates reconstructed in the $ B ^+ \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} K ^+ $ decay for (left) $ p\mkern 1mu\mathrm{Pb}$ and (right) $\mathrm{Pb}\mkern 1.5mu p$ collisions, with the fit result superimposed. The solid blue line, the solid green line, the crossshaded area, the brown shaded area and the red shaded area represent the total fit, the signal component, the partially reconstructed background, the combinatorial background and $ B ^+ \rightarrow { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu} \pi ^+ $ decays, respectively. 
Fig2a.pdf [25 KiB] HiDef png [299 KiB] Thumbnail [215 KiB] *.C file 

Fig2b.pdf [23 KiB] HiDef png [311 KiB] Thumbnail [230 KiB] *.C file 

Invariant mass distribution of $ B ^0 $ candidates reconstructed in the $ B ^0 \rightarrow D ^ \pi ^+ $ decay for (left) $ p\mkern 1mu\mathrm{Pb}$ and (right) $\mathrm{Pb}\mkern 1.5mu p$ collisions, with the fit result superimposed. The solid blue line, the solid green line, the crossshaded area, the brown shaded area and the red shaded area represent the total fit, the signal component, the partially reconstructed background, the combinatorial background and $ B ^0 \rightarrow D ^ K ^+ $ decays, respectively. 
Fig3a.pdf [26 KiB] HiDef png [411 KiB] Thumbnail [255 KiB] *.C file 

Fig3b.pdf [28 KiB] HiDef png [425 KiB] Thumbnail [279 KiB] *.C file 

Invariant mass distribution of $\Lambda ^0_ b $ candidates reconstructed in the $\Lambda ^0_ b \rightarrow \Lambda ^+_ c \pi ^ $ decay for (left) $ p\mkern 1mu\mathrm{Pb}$ and (right) $\mathrm{Pb}\mkern 1.5mu p$ collisions, with the fit result superimposed. The solid blue line, the solid green line, the crossshaded area, the brown shaded area and the red shaded area represent the total fit, the signal component, the partially reconstructed background, the combinatorial background and $\Lambda ^0_ b \rightarrow \Lambda ^+_ c K ^ $ decays, respectively. 
Fig4a.pdf [27 KiB] HiDef png [337 KiB] Thumbnail [226 KiB] *.C file 

Fig4b.pdf [28 KiB] HiDef png [355 KiB] Thumbnail [255 KiB] *.C file 

Production crosssection of $ B ^+ $ mesons as a function of (left) $ p_{\mathrm{T}}$ in $y$ bins and (right) $y$ integrated over $ p_{\mathrm{T}}$ . The vertical bars (boxes) show statistical (total) uncertainties. 
Fig5a.pdf [15 KiB] HiDef png [173 KiB] Thumbnail [166 KiB] *.C file 

Fig5b.pdf [14 KiB] HiDef png [101 KiB] Thumbnail [108 KiB] *.C file 

Production crosssection of $ B ^0 $ mesons as function of (left) $ p_{\mathrm{T}}$ in $y$ bins and (right) $y$ integrated over $ p_{\mathrm{T}}$ . The vertical bars (boxes) show statistical (total) uncertainties. 
Fig6a.pdf [15 KiB] HiDef png [178 KiB] Thumbnail [171 KiB] *.C file 

Fig6b.pdf [14 KiB] HiDef png [108 KiB] Thumbnail [117 KiB] *.C file 

Production crosssection of $\Lambda ^0_ b $ baryons as a function of (left) $ p_{\mathrm{T}}$ in $y$ bins and (right) $y$ integrated over $ p_{\mathrm{T}}$ . The vertical bars (boxes) show statistical (total) uncertainties. 
Fig7a.pdf [15 KiB] HiDef png [179 KiB] Thumbnail [171 KiB] *.C file 

Fig7b.pdf [14 KiB] HiDef png [111 KiB] Thumbnail [118 KiB] *.C file 

Production crosssection ratios of $\Lambda ^0_ b $ baryons over $ B ^0 $ mesons and of $ B ^0 $ mesons over $ B ^+ $ mesons (top) as a function of $y$ integrated over $ p_{\mathrm{T}}$ and as a function of $ p_{\mathrm{T}}$ for (bottom left) $2.5<y<3.5$ and (bottom right) $3.5<y<2.5$. The vertical bars (boxes) show statistical (total) uncertainties. 
Fig8a.pdf [14 KiB] HiDef png [132 KiB] Thumbnail [138 KiB] *.C file 

Fig8b.pdf [14 KiB] HiDef png [119 KiB] Thumbnail [121 KiB] *.C file 

Fig8c.pdf [14 KiB] HiDef png [119 KiB] Thumbnail [122 KiB] *.C file 

Forwardbackward ratio, $R_{\rm FB}$, for $ B ^+ $ mesons as a function of (left) $ p_{\mathrm{T}}$ and (right) $y$ in protonlead collisions compared with HELAConia calculations using different nPDF sets. For the data points, the vertical bars (boxes) represent the statistical (total) uncertainties. The rapidity range $2.5< y <3.5$ represents the common range between the $ p\mkern 1mu\mathrm{Pb}$ and $\mathrm{Pb}\mkern 1.5mu p$ data samples studied in this analysis. 
Fig9a.pdf [15 KiB] HiDef png [876 KiB] Thumbnail [324 KiB] *.C file 

Fig9b.pdf [15 KiB] HiDef png [711 KiB] Thumbnail [286 KiB] *.C file 

Forwardbackward ratio, $R_{\rm FB}$, of $ B ^0 $ mesons as a function of (left) $ p_{\mathrm{T}}$ and as a function of (right) $y$ in protonlead collisions. The vertical bars (boxes) represent the statistical (total) uncertainties. 
Fig10a.pdf [14 KiB] HiDef png [124 KiB] Thumbnail [134 KiB] *.C file 

Fig10b.pdf [13 KiB] HiDef png [111 KiB] Thumbnail [120 KiB] *.C file 

Forwardbackward ratio, $R_{\rm FB}$, of $\Lambda ^0_ b $ baryons as a function of (left) $ p_{\mathrm{T}}$ and (right) $y$ in protonlead collisions. The vertical bars (boxes) represent the statistical (total) uncertainties. 
Fig11a.pdf [14 KiB] HiDef png [125 KiB] Thumbnail [135 KiB] *.C file 

Fig11b.pdf [13 KiB] HiDef png [112 KiB] Thumbnail [120 KiB] *.C file 

Forwardbackward ratio, $R_{\rm FB}$, of (red) $ B ^+ $, (blue) $ B ^0 $ mesons and (green) $\Lambda ^0_ b $ baryons as a function of (left) $ p_{\mathrm{T}}$ and (right) $y$ in protonlead collisions. The vertical bars (boxes) represent the statistical (total) uncertainties. Data points are shifted horizontally for better visibility. 
Fig12a.pdf [15 KiB] HiDef png [165 KiB] Thumbnail [162 KiB] *.C file 

Fig12b.pdf [14 KiB] HiDef png [129 KiB] Thumbnail [136 KiB] *.C file 

Ratio of nuclear modification factors, $R^{\Lambda ^0_ b / B ^0 }_{ p\mkern 1mu\mathrm{Pb} }$, as a function of (left) $ p_{\mathrm{T}}$ and (right) $y$ in $ p\mkern 1mu\mathrm{Pb} $ and $\mathrm{Pb}\mkern 1.5mu p $ collisions. The vertical bars (boxes) represent the statistical (total) uncertainties. 
Fig13a.pdf [14 KiB] HiDef png [146 KiB] Thumbnail [149 KiB] *.C file 

Fig13b.pdf [14 KiB] HiDef png [113 KiB] Thumbnail [123 KiB] *.C file 

Nuclear modification factor, $R_{ p\mkern 1mu\mathrm{Pb} }$, for $ B ^+ $ mesons as function of (top) $y$ and as a function of $ p_{\mathrm{T}}$ in (bottom left) $ p\mkern 1mu\mathrm{Pb} $ and (bottom right) $\mathrm{Pb}\mkern 1.5mu p $ compared with HELAConia calculations using different nPDF sets as well as with the measurement of $R_{ p\mkern 1mu\mathrm{Pb} }$ for nonprompt $ { J \mskip 3mu/\mskip 2mu\psi \mskip 2mu}$ production. For the data points, the vertical bars (boxes) represent the statistical (total) uncertainties. 
Fig14a.pdf [17 KiB] HiDef png [649 KiB] Thumbnail [280 KiB] *.C file 

Fig14b.pdf [16 KiB] HiDef png [586 KiB] Thumbnail [266 KiB] *.C file 

Fig14c.pdf [16 KiB] HiDef png [364 KiB] Thumbnail [227 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2018048.gif Thumbnail 
Signal yields in the range $2< p_{\mathrm{T}} <20\text{ Ge V /}c $ and $1.5< y <3.5$ ($4.5< y <2.5$) for $ p\mkern 1mu\mathrm{Pb}$ ( $\mathrm{Pb}\mkern 1.5mu p$ ) collisions. Uncertainties are statistical only. 
Table_1.pdf [50 KiB] HiDef png [88 KiB] Thumbnail [42 KiB] tex code 

Summary of systematic uncertainties (in %) for the measured crosssections for different decay modes in $ p\mkern 1mu\mathrm{Pb}$ . The ranges correspond to the minimum and maximum values over the $ p_{\mathrm{T}}$ and $ y$ bins of the measurement. 
Table_2.pdf [56 KiB] HiDef png [102 KiB] Thumbnail [48 KiB] tex code 

Summary of systematic uncertainties (in %) for the measured crosssections for different decay modes in $\mathrm{Pb}\mkern 1.5mu p$ . The ranges correspond to the minimum and maximum values over the $ p_{\mathrm{T}}$ and $ y$ bins of the measurement. 
Table_3.pdf [56 KiB] HiDef png [102 KiB] Thumbnail [48 KiB] tex code 

Differential crosssections of $ B ^+ $ , $ B ^0$ and $\Lambda ^0_ b $ production in bins of $ p_{\mathrm{T}}$ and $y$, $\frac{{\rm d}^2\sigma}{{\rm d} p_{\mathrm{T}} {\rm d}y}( \upmu\text{b} /[ \text{ Ge V /}c ])$, and in bins of $y$ integrated over $2< p_{\mathrm{T}} <20\text{ Ge V /}c $, $\frac{{\rm d}\sigma}{{\rm d}y} ( \upmu\text{b} )$. The first uncertainty is statistical and the second systematic. 
Table_4.pdf [44 KiB] HiDef png [176 KiB] Thumbnail [92 KiB] tex code 

Forwardbackward ratios, $R_{\rm FB}$, of $ B ^+ $ , $ B ^0$ and $\Lambda ^0_ b $ production in bins of $ p_{\mathrm{T}}$ and integrated over $2.5<y<3.5$. The first uncertainty is statistical and the second systematic. 
Table_5.pdf [42 KiB] HiDef png [76 KiB] Thumbnail [40 KiB] tex code 

Ratios of nuclear modification factors, $R^{\Lambda ^0_ b / B ^0 }_{ p\mkern 1mu\mathrm{Pb} }$, in bins of $ p_{\mathrm{T}}$ and integrated over $2.5<y<3.5$, for $ p\mkern 1mu\mathrm{Pb} $ and $\mathrm{Pb}\mkern 1.5mu p $ samples. The first uncertainty is statistical and the second systematic. 
Table_6.pdf [35 KiB] HiDef png [91 KiB] Thumbnail [48 KiB] tex code 

Nuclear modification factor, $R_{ p\mkern 1mu\mathrm{Pb} }$, of $ B ^+ $ production in $ p\mkern 1mu\mathrm{Pb} $ and $\mathrm{Pb}\mkern 1.5mu p $ collisions, in bins of $ p_{\mathrm{T}}$ for the range $2.5<y<3.5$. The first uncertainty is statistical and the second systematic. 
Table_7.pdf [35 KiB] HiDef png [93 KiB] Thumbnail [48 KiB] tex code 
Created on 19 October 2019.