A search for charmless fourbody decays of $\Lambda_{b}^{0}$ and $\Xi_{b}^{0}$ baryons with a proton and three charged mesons (either kaons or pions) in the final state is performed. The data sample used was recorded in 2011 and 2012 with the LHCb experiment and corresponds to an integrated luminosity of 3 fb$^{1}$. Six decay modes are observed, among which $\Lambda_{b}^{0} \to pK^{}\pi^{+}\pi^{}$, $\Lambda_{b}^{0} \to pK^{}K^{+}K^{}$, $\Xi_{b}^{0} \to pK^{}\pi^{+}\pi^{}$ and $\Xi_{b}^{0} \to pK^{} \pi^{+}K^{}$ are established for the first time. Their branching fractions (including the ratio of hadronisation fractions in the case of the $\Xi_{b}^{0}$ baryon) are determined relative to the $\Lambda_{b}^{0} \to \Lambda_{c}^{+}\pi^{}$ decay.
Results of the fit to the $ p \pi ^\pi ^+\pi ^$ candidate mass spectrum with (top) linear and (bottom) logarithmic scales. The different components employed in the fit are indicated in the legend. The $\Lambda ^0_ b \rightarrow $ 5body legend describes two components, the radiative partially reconstructed background $\Lambda ^0_ b \rightarrow p \pi ^\eta ^{\prime} $ and the partially reconstructed background $\Lambda ^0_ b \rightarrow p \pi ^\pi ^+\pi ^\pi ^0$ where a $\pi^0$ is not reconstructed. The latter has a lowermass endpoint. 
Lb2ppi[..].pdf [49 KiB] HiDef png [538 KiB] Thumbnail [147 KiB] *.C file 

Results of the fit to the $ p K ^\pi ^+\pi ^$ candidate mass spectrum with (top) linear and (bottom) logarithmic scales. The different components employed in the fit are indicated in the legend. The $\Lambda ^0_ b \rightarrow $ 5body legend describes two components, the radiative partially reconstructed background $\Lambda ^0_ b \rightarrow p K ^\eta ^{\prime} $ and the partially reconstructed background $\Lambda ^0_ b \rightarrow p K ^\pi ^+\pi ^\pi ^0$ where a $\pi^0$ is not reconstructed. The latter has a lowermass endpoint. 
Lb2pKp[..].pdf [51 KiB] HiDef png [599 KiB] Thumbnail [159 KiB] *.C file 

Results of the fit to the $ p K ^ K ^+\pi ^$ candidate mass spectrum with (top) linear and (bottom) logarithmic scales. The different components employed in the fit are indicated in the legend. The $\Lambda ^0_ b \rightarrow $ 5body legend describes two components where a $\pi^0$ is not reconstructed, the partially reconstructed background $\Lambda ^0_ b \rightarrow p K ^\pi ^+\pi ^\pi ^0$ where a pion is misidentified as a kaon and the partially reconstructed background $\Lambda ^0_ b \rightarrow p K ^ K ^+\pi ^\pi ^0$ . 
Lb2pKK[..].pdf [59 KiB] HiDef png [707 KiB] Thumbnail [178 KiB] *.C file 

Results of the fit to the $ p K ^\pi ^+ K ^$ candidate mass spectrum with (top) linear and (bottom) logarithmic scales. The different components employed in the fit are indicated in the legend. 
Xib02p[..].pdf [53 KiB] HiDef png [674 KiB] Thumbnail [174 KiB] *.C file 

Results of the fit to the $ p K ^ K ^+ K ^$ candidate mass spectrum with (top) linear and (bottom) logarithmic scales. The different components employed in the fit are indicated in the legend. The $\Lambda ^0_ b \rightarrow $ 5body legend includes two decays, partially reconstructed $\Lambda ^0_ b \rightarrow p K ^ K ^+ K ^\gamma $ and $\Lambda ^0_ b \rightarrow p K ^ K ^+ K ^\pi ^0$ , where the $\gamma $ and $\pi^0$ are not reconstructed. 
Lb2pKK[..].pdf [52 KiB] HiDef png [598 KiB] Thumbnail [156 KiB] *.C file 

Results of the fit to the $\Lambda ^+_ c \pi^$ candidate mass spectrum on linear scale. The different components employed in the fit are indicated in the legend. 
Lb2pKp[..].pdf [20 KiB] HiDef png [191 KiB] Thumbnail [59 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2017034.gif Thumbnail 
Signal yields for each decay mode, determined by summing the fitted yields in each year of data taking. The signal ($S$) to background ($B$, adding all sources) ratios in an invariant mass window, covering $\pm$ 3 times the measured signal widths, are provided. The corresponding invariant mass ranges are reported in the fourth column. 
Table_1.pdf [54 KiB] HiDef png [107 KiB] Thumbnail [35 KiB] tex code 

Ratios of the normalisation decay mode efficiencies, relative to the signal decay mode as used in Eq. 4, for (first row) 2011 and (second row) 2012. The last column shows the efficiency of the veto of charm and charmonium backgrounds (applied to the signal mode only), as discussed in the text. Since the $\Xi ^0_ b \rightarrow p K ^ K ^+ K ^$ decay mode is not observed, the veto efficiency is determined with the simulated data sample. The difference between the simulation value and the average veto efficiency measured on other $\Xi ^0_ b $ modes is reported in the table as the uncertainty. 
Table_2.pdf [54 KiB] HiDef png [183 KiB] Thumbnail [30 KiB] tex code 

Yields and efficiencies of each signal decay with the statistical uncertainty, and systematic uncertainties related to the fit model and the efficiency determination, for the 2011 data samples. 
Table_3.pdf [76 KiB] HiDef png [76 KiB] Thumbnail [12 KiB] tex code 

Yields and efficiencies of each signal decay with the statistical uncertainty, and systematic uncertainties related to the fit model and the efficiency determination, for the 2012 data samples. 
Table_4.pdf [76 KiB] HiDef png [76 KiB] Thumbnail [13 KiB] tex code 

Measurements of the $R$ ratio from the (first row) 2011 and the (second row) 2012 data samples for $\Lambda ^0_ b $ decay modes expressed in percent as well as their combination. The three uncertainties are statistical, systematic related to the fit model and systematic related to the efficiency, respectively. The consistency of the two determinations for each year, denoted $\Delta$, is quantified as the ratio of the signed difference of the central values over the quadratic sum of the related uncertainties. 
Table_5.pdf [53 KiB] HiDef png [94 KiB] Thumbnail [28 KiB] tex code 

Measurements of the $R$ ratio from the (first row) 2011 and the (second row) 2012 data samples for $\Xi ^0_ b $ decay modes expressed in per mil as well as their combination. The three uncertainties are statistical, systematic related to the fit model and systematic related to the efficiency, respectively. The consistency of the two determinations for each year, denoted $\Delta$, is quantified as the ratio of the signed difference of the central values over the quadratic sum of the related uncertainties. 
Table_6.pdf [53 KiB] HiDef png [128 KiB] Thumbnail [22 KiB] tex code 
Created on 23 February 2019.Citation count from INSPIRE on 23 February 2019.