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 [565 KiB] Thumbnail [330 KiB] *.C file 

Fig1b.pdf [24 KiB] HiDef png [570 KiB] Thumbnail [332 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 [291 KiB] Thumbnail [215 KiB] *.C file 

Fig2b.pdf [23 KiB] HiDef png [306 KiB] Thumbnail [229 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 [403 KiB] Thumbnail [253 KiB] *.C file 

Fig3b.pdf [28 KiB] HiDef png [420 KiB] Thumbnail [277 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 [329 KiB] Thumbnail [224 KiB] *.C file 

Fig4b.pdf [28 KiB] HiDef png [351 KiB] Thumbnail [253 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 [172 KiB] Thumbnail [166 KiB] *.C file 

Fig5b.pdf [14 KiB] HiDef png [102 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 [109 KiB] Thumbnail [118 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 [112 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 [139 KiB] *.C file 

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

Fig8c.pdf [14 KiB] HiDef png [120 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 [875 KiB] Thumbnail [326 KiB] *.C file 

Fig9b.pdf [15 KiB] HiDef png [709 KiB] Thumbnail [287 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 [97 KiB] Thumbnail [58 KiB] *.C file 

Fig10b.pdf [13 KiB] HiDef png [87 KiB] Thumbnail [52 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 [98 KiB] Thumbnail [59 KiB] *.C file 

Fig11b.pdf [13 KiB] HiDef png [88 KiB] Thumbnail [52 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 [166 KiB] Thumbnail [162 KiB] *.C file 

Fig12b.pdf [14 KiB] HiDef png [130 KiB] Thumbnail [137 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 [281 KiB] *.C file 

Fig14b.pdf [16 KiB] HiDef png [585 KiB] Thumbnail [268 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 [66 KiB] HiDef png [74 KiB] Thumbnail [36 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 [72 KiB] HiDef png [85 KiB] Thumbnail [37 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 [72 KiB] HiDef png [85 KiB] Thumbnail [37 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 [60 KiB] HiDef png [154 KiB] Thumbnail [75 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 [59 KiB] HiDef png [66 KiB] Thumbnail [33 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 [51 KiB] HiDef png [78 KiB] Thumbnail [40 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 [51 KiB] HiDef png [80 KiB] Thumbnail [40 KiB] tex code 
Created on 06 March 2021.