Measurements of CP asymmetries in charmless four-body Λb0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varLambda } ^0_{b} $$\end{document} and Ξb0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varXi } ^0_{b} $$\end{document} decays

A search for CP\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C\!P$$\end{document} violation in charmless four-body decays of Λb0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varLambda } ^0_{b} $$\end{document} and Ξb0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varXi } ^0_{b} $$\end{document} baryons with a proton and three charged mesons in the final state is performed. To cancel out production and detection charge-asymmetry effects, the search is carried out by measuring the difference between the CP\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C\!P$$\end{document} asymmetries in a charmless decay and in a decay with an intermediate charmed baryon with the same particles in the final state. The data sample used was recorded in 2011 and 2012 with the LHCb detector and corresponds to an integrated luminosity of 3fb-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3 \text { fb} ^{-1} $$\end{document}. A total of 18 CP\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C\!P$$\end{document} asymmetries are considered, either accounting for the full phase space of the decays or exploring specific regions of the decay kinematics. No significant CP\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C\!P$$\end{document}-violation effect is observed in any of the measurements.


Introduction
All measurements of CP violation performed so far are consistent with the predictions of the Standard Model (SM) [1]. Nonvanishing CP-violating asymmetries have been observed in the decays of both K and B mesons [2]. In contrast, CP violation has not been observed in baryon decays, although some indications for nonvanishing CP asymmetries in bflavoured baryon decays have been reported by the LHCb collaboration [3][4][5][6].
The abundant production of Λ 0 b and Ξ 0 b baryons 1 in proton-proton collisions at the Large Hadron Collider (LHC) gives the LHCb experiment the opportunity to study multibody charmless decays of b-flavoured baryons. In particular, Λ 0 b and Ξ 0 b baryon decays to charmless four-body final states were observed by the LHCb collaboration and their branching fractions measured [7]. Their large yields enable measurements of CP-violating asymmetries to be performed with a precision at the level of a few percent.
This search follows the successful path of the observation of large CP-violating asymmetries in multibody charmless decays of charged and neutral B mesons by LHCb [8][9][10][11]. These decays proceed simultaneously through the chargedcurrent b → u transition and neutral-current b → s, d transitions, and the resulting interference exhibits a weakphase difference. Furthermore, and analogously to the aforementioned charmless multibody B-meson decays, charmless multibody decays of b-flavoured baryons contain rich resonance structures, both in the two-or three-body baryonic invariant-mass spectra (i.e. pK − , pπ − , pπ + , pπ + π − and pK − π + ) and in the two-or three-body nonbaryonic ones (i.e. the π + π − , K ± π ∓ , K + K − , π + π − π ± and K ± π + π − ). Consequently, CP asymmetries might be enhanced due to the strong-phase differences induced by the interference patterns between these transitions in the mass regions around resonances. The charmless b-baryon decays studied in this paper are hence well suited for a potential first observation of CP violation in the baryon sector. However, the presence of these strong phases, that are difficult to predict, would make a potential observation of CP violation difficult to interpret in terms of the weak phase of the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [12,13].
This work focuses on a search for CP violation in X 0 b → phh h charmless decays, where X 0 b stands either for Λ 0 b or Ξ 0 b and h ( , ) stand either for a pion or a kaon. Six decays are studied, The CP asymmetry is defined as where (X 0 b → f ) is the partial width of the given decay. The CP asymmetry measurement relies on counting the number of reconstructed particle and antiparticle decays and includes therefore experimental charge-asymmetric effects such as Table 1 Four-body charmless and charmed decays considered in this analysis. The difference of CP-asymmetries measured for the charmless modes and for the control channels results in ΔA CP measurements. For each observable, the choice of the control channel is aiming at cancelling at first order production and detection asymmetries. Given the data samples at hand, it is not possible to meet both criteria for the signal decay Ξ 0 b → pK − π + K − : the choice of the Cabibbo-favoured decay Ξ 0 b → (Ξ + c → pK − π + )π − as a control channel requires in turn to correct the corresponding ΔA CP for the kaon-detection asymmetry

Charmless mode
Control channel the track detection efficiency or b-baryon production asymmetries. They are cancelled out to first order by comparing the CP asymmetries of the signal modes to those of charmed decays that lead to the same or very similar final states and for which no measurable CP violation is expected in the SM [14]. The decays Λ 0 b → (Λ + c → pK − π + )π − , Λ 0 b → (Λ + c → pπ − π + )π − and Ξ 0 b → (Ξ + c → pK − π + ) π − are thus reconstructed with the same selection as the corresponding charmless signals. The CP-violating observable considered in this work is then referred to as ΔA CP ≡ A CP no-c −A CP c , where A CP no-c (A CP c ) is the asymmetry measured in the charmless (charmed) decays. The decays of interest are reported in Table 1.
In addition to ΔA CP measurements integrated over all of the four-body phase space, specific regions of the space are studied in order to search for local CP asymmetries.
The same final states have been used by the LHCb experiment to search for CP violation using triple product asymmetries [5,15]. The latter technique and the ΔA CP measurements exhibit different sensitivity to CP violation [16], which makes the two approaches complementary.

Detector and data set
The analysis is performed using pp collision data recorded with the LHCb detector, corresponding to an integrated luminosity of 1.0 fb −1 at a centre-of-mass energy of 7 TeV in 2011 and 2.0 fb −1 at a centre-of-mass energy of 8 TeV in 2012. The LHCb detector [17,18] is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 5, designed for the study of particles containing b or c quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction region, a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes placed downstream of the magnet. The tracking system provides a measurement of the momentum, p, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200 GeV/c. The minimum distance of a track to a primary vertex (PV), the impact parameter (IP), is measured with a resolution of (15+29/ p T ) µm, where p T is the component of the momentum transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov detectors. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers.
Simulation is used to investigate backgrounds from other b-hadron decays and also to study the detection and reconstruction efficiencies of the signals. In the simulation, pp collisions are generated using Pythia [19,20] with a specific LHCb configuration [21]. Decays of hadronic particles are described by EvtGen [22] in which final-state radiation is generated using Photos [23]. The interactions of the generated particles with the detector, and its response, are implemented using the Geant4 toolkit [24,25] as described in Ref. [26].

Trigger and selection requirements
The selection follows most of the strategy described in Ref. [7]. The online event selection is performed by a trigger [27] that consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, in which all charged particles with p T > 500 (300) MeV/c are reconstructed for 2011 (2012) data. At the hardware-trigger stage, events are required to include a muon or a dimuon with high transverse momentum or a hadron, photon or electron with high transverse energy. The software trigger reconstructs charged particles with transverse momentum p T > 500 (300) MeV/c for 2011 (2012) data and requires a two-, three-or four-track secondary vertex with significant displacement from all primary pp interaction vertices. At least one charged particle must have transverse momentum p T > 1.7 (1.6) GeV/c for 2011 (2012) data and be inconsistent with originating from any PV. A multivariate algorithm [28] is used for the identification of secondary vertices consistent with the decay of a b hadron. In the offline selection, trigger signals are associated with reconstructed particles. Selection requirements can therefore be made on the trigger selection itself and on whether the decision was due to the signal candidate, other particles produced in the pp collision, or a combination of both.
The events passing the trigger requirements are filtered in two stages. Initial requirements are applied to further reduce the size of the data sample before a multivariate selection is applied. Selection requirements based on topological variables, such as the flight distance of the b baryon, are used as the main discriminants. In order to preserve the phase space of the decays of interest, only loose requirements are placed on the transverse momenta of the decay products, Neutral b-baryon candidates, hereafter denoted as X 0 b , are formed from a proton candidate selected with particle identification (PID) requirements and three additional charged tracks. When more than one PV is reconstructed, the X 0 b candidate is associated to the PV with the smallest value of χ 2 IP , where χ 2 IP is the difference in χ 2 of a given PV reconstructed with and without the considered candidate. Each of the four tracks of the final state is required to have χ 2 IP > 16 and 3 < p < 100 GeV/c. Beyond 100 GeV/c, there is little pion/kaon/proton discrimination. The X 0 b candidates are then required to form a vertex with a fit quality χ 2 vtx < 20 and to be significantly separated from any PV with χ 2 FD > 50, where χ 2 FD is the square of the flight-distance significance. To remove backgrounds from higher-multiplicity decays, the difference in χ 2 vtx when adding any other track must be greater than 4. The X 0 b candidates must have a transverse momentum p T (X 0 b ) greater than 1.5 GeV/c and an invariant mass within the range 5340 < m( phh h ) < 6400 MeV/c 2 . They are further required to be consistent with originating from a PV, quantified by both χ 2 IP < 16 and the cosine of the angle θ DIR between the reconstructed momentum of the b hadron and the vector defined by the associated PV and the decay vertex be greater than 0.999. Finally, PID requirements are applied to provide discrimination between kaons and pions in order to assign the candidates to one of the five different final-state hypotheses pπ − π + π − , pK − π + π − , pK − K + π − , pK − π + K − and pK − K + K − .
There are three main categories of background that contribute significantly in the selected invariant-mass regions: the so-called signal cross-feed background, resulting from a misidentification of one or more final-state particles in a charmless baryon decay, which can therefore be reconstructed as another charmless decay with a different mass hypothesis; the charmless decays of neutral B mesons to final states containing four charged mesons, where a pion or a kaon is misidentified as a proton; and the combinatorial background, which results from a random association of unrelated tracks. The pion and kaon PID requirements, that define mutually exclusive samples, are optimised to reduce the cross-feed background, and hence to maximise the significance of the signal. The charmless B-meson decays are identified by reconstructing the invariant-mass distri-butions of candidates using the pion or kaon mass instead of the proton mass hypothesis, in the high-mass sidebands defined as m sideband < m( phh h ) < 6400 MeV/c 2 , where m sideband = 5680 MeV/c 2 for pπ − π + π − and pK − K + π − final states, and m sideband = 5840 MeV/c 2 for pK − π + π − , pK − π + K − and pK − K + K − final states. This background contribution is reduced by the optimisation of the proton PID requirement.
To reject combinatorial background, multivariate discriminants based on a boosted decision tree (BDT) [29] with the AdaBoost algorithm [30] have been designed. Candidates from simulated Λ 0 b → pπ − π + π − decays and the high-mass sideband are used as the signal and background training samples, respectively. This high-mass sideband region is chosen such that the sample is free of cross-feed background. The samples are divided into two data-taking periods and further subdivided into two equally sized subsamples. Each subsample is then used to train an independent discriminant. The BDT trained on one subsample is used to select candidates from the other subsample, in order to avoid a possible bias in the selection.
The BDTs have the following quantities as inputs: p T , η, χ 2 IP , χ 2 FD , cos θ DIR , and χ 2 vtx of the X 0 b candidate; the smallest change in the b-baryon χ 2 vtx when adding any other track from the event; the sum of the χ 2 IP of the four tracks of the final state; and the p T asymmetry where p cone T is the transverse component of the vector sum of all particle momenta inside a cone around the b-baryon candidate direction, of radius R ≡ δη 2 + δφ 2 = 1.5, where δη and δφ are the difference in pseudorapidity and azimuthal angle (expressed in radians) around the beam direction, between the momentum vector of the track under consideration and that of the b-hadron candidate. The distribution of p asym T for the signal candidates is enhanced towards high values. The BDT output is determined to be uncorrelated with the position in phase space of the decays of interest. The selection requirement placed on the output of the BDTs is optimised for the six decays of interest by minimising the uncertainties on the CP-asymmetry differences.
A number of background contributions consisting of fully reconstructed b-baryon decays into the two-body Λ + c h, Ξ + c h, three-body Dph or (cc) ph combinations, where (cc) represents a charmonium resonance, may produce the same final state as the signal. Hence, they have similar invariant-mass distribution of the b-baryon candidate as the signal along with a similar selection efficiency. The presence of a misidentified hadron in the D, Λ + c and Ξ + c decay also produces peaking background under the signal. Therefore, the following decay channels are explicitly reconstructed under the relevant par-ticle hypotheses and vetoed by means of a requirement on the resulting invariant mass, in all spectra: The decays of other possible broad charmonium resonances to π + π − and K + K − are retained as potential interfering amplitudes with the charmless amplitudes under study.
The same set of trigger, PID and BDT requirements is applied to the control modes to cancel out most of the systematic effects related to the selection criteria. Candidates whose pK − π + or pπ − π + invariant mass is in the range [2213, 2313] MeV/c 2 for Λ + c and [2437, 2497] MeV/c 2 for Ξ + c , are retained as control channels candidates. Events outside these intervals belong to the corresponding signal spectrum, again ensuring statistically independent samples for the simultaneous fit.
The fraction of events containing more than one candidate is below the percent level. The candidate to be retained in each event is chosen randomly and reproducibly.

Simultaneous fit
A simultaneous unbinned extended maximum-likelihood fit is performed to the invariant-mass distributions of the bhadron candidates under each of the mass hypotheses for the signal and control channel final-state tracks. The data samples are split according to the charge of the proton and to the year of data taking. Furthermore, data are split according to the hardware trigger conditions, in order to correct raw measurements for charge-asymmetric trigger efficiencies. The components of the model include, in addition to signal decays, partially reconstructed five-body X 0 b decays, signal and background cross-feeds, four-and five-body decays of B mesons and combinatorial background. The independent data samples obtained for each final state are fitted simultaneously. For each sample, the likelihood is expressed as where N j is the number of events related to the component j and P j,i is the probability distribution function for component j evaluated at the mass of the candidate i.

Fit model
The signal decays are modelled as the sum of two Crystal Ball (CB) functions [31] that share peak positions and widths but have independent power-law tails on opposite sides of the peak. The Λ 0 b mass parameter is free in the fit and shared among the Λ 0 b decays. The difference between the fitted Ξ 0 b and Λ 0 b masses is also a shared parameter and is constrained to the value reported in Ref. [2] by using a Gaussian function.
The width parameter for Λ 0 b → pK − π + π − decays measured in the 2012 data-taking sample is found to be 16.47 ± 0.22 MeV/c 2 and is chosen as reference. The ratio of the experimental widths of the signal decay functions is constrained using Gaussian prior probability distributions multiplying the likelihood function, with parameters obtained from a fit to simulated events. The other parameters of the CB components are obtained from a simultaneous fit to simulated samples, and are fixed to those values in the fits to the data.
The cross-feed backgrounds are modelled by the sum of two CB functions, whose parameters are determined from simulated samples weighted to match the performances of the particle identification algorithm as measured in the data. All cases resulting from the misidentification of either one or two of the final-state particles are considered. The yield of each misidentified decay is constrained to the yield of the corresponding correctly identified decay and the known misidentification probabilities. These constraints are implemented using Gaussian prior probability distributions multiplying the likelihood function. Their mean values are obtained from the ratio of selection efficiencies and their widths include uncertainties originating from the finite size of the simulated events samples as well as the systematic uncertainties related to the determination of the PID efficiencies.
The backgrounds resulting from four-or five-body decays of B mesons are identified in each spectrum by a dedicated fit to the candidates in the high-mass sideband, reconstructed under the hypothesis of the kaon mass for the proton candidates. The relative yield of each decay is then constrained in the simultaneous fit from its observed abundance in the high-mass sidebands. The invariant-mass distributions are modelled by the sum of two CB functions, whose parameters are determined from simulation.
Partially reconstructed backgrounds where a neutral pion is not reconstructed, such as Λ 0 b , Ξ 0 b → phh h π 0 , are modelled by means of generalised ARGUS functions [32] convolved with a Gaussian resolution function. The Gaussian width is taken as the signal Λ 0 b → pK − π + π − width parameter. The parameters of the ARGUS function are shared among all invariant-mass spectra and are determined directly from the fit, except for the threshold, which is given by m(X b ) − m(π 0 ). Partially reconstructed decays with a missing photon such as Λ 0 b → pπ − η and Λ 0 b → pK − η decays, with η → π + π − γ , are modelled separately using the same functional form but where the parameters are fixed from simulation. The Λ 0 b → pK − π + π − π 0 decay modes where a charged pion is misidentified as a kaon can significantly contribute to the pK − K + π − and pK − π + K − spectra. They are modelled with an empirical function determined from the partially reconstructed background candidates in the control channel. Finally, the combinatorial background is modelled by a linear function whose slope is shared among the invariantmass spectra.

The ensemble of measurements
The following three categories of measurements have been considered a priori (before any evaluation of the data) to search for global and local effects of CP violation.
• CP asymmetries are measured, considering the whole selected phase space of the decay candidates. • CP asymmetries are also measured in the phase-space region of low invariant mass on the baryonic pair (i.e. pπ ± or pK − ) and low invariant mass on the pairing of the two other tracks. The ensemble of measurements that are performed with this phase-space selection is hereafter referred to as LBM (Low 2 × 2-Body Mass) measurements. The invariant mass of the baryonic pair is required to be lower than 2 GeV/c 2 while the invariant-mass requirements on the two remaining tracks depends on whether it is a π + π − pair, a K ± π ∓ or a K K pair. These values are chosen to include several known resonances, in particular f 0 (1500) resonance for π + π − , the broad scalar K * 0 (1430) 0 resonance for K + π − and the f 2 (1525) resonance for K + K − . Only the modes with the largest signal yields are considered, The two-body low-mass distributions are displayed in Fig. 1. Several resonant structures are observed, and correspond to baryon resonances like Λ(1520), Δ(1232) ++ and N (1520) or meson resonances like K * (892) 0 , ρ(770) 0 or φ(1020). This phase-space selection focuses therefore on low-invariant-mass resonances (both mesonic and baryonic) as well as low-invariant-mass nonresonant components of the amplitudes. The latter have been shown to generate large CP-violating asymmetries in analogous B-meson decays [8].
• CP asymmetries are measured for regions of the phase space that contain specific quasi-two-body decays, where in the latter only the high pK − mass region is selected. Invariantmass requirements for these measurements are reported in Table 2. Only the narrower baryons or the well-known baryon and meson resonances have been considered, with the noticeable exception of the a 1 (1260) meson. Although the a 1 (1260) meson is a broad resonance, the analogous B-meson decay B 0 → a 1 (1260) ± π ∓ has been studied at the B-factories [33,34] and could serve as a benchmark comparison in the interpretation of the results obtained for the Λ 0 b → pa 1 (1260) − decay.

Corrections for experimental detection asymmetries and related systematic uncertainties
Tracking reconstruction, trigger selection and particle identification requirements can generate charge-dependent selection efficiencies of the decays of interest. Most of these charge-dependent effects are however cancelled out in the ΔA CP observables, up to the kinematical differences between signal and control channels. The remaining impact is addressed by evaluating corrections to the ΔA CP observables. These correction factors are determined from calibration samples as discussed in this Section. Systematic uncertainties are estimated for each correction factor and propagated to the ΔA CP measurements. A summary of the systematic uncertainties is reported in Table 3 for all modes.
• Tracking detection efficiency: differences between the interactions of oppositely charged pions, kaons or protons in the material of the spectrometer induce detection charge asymmetries. The difference in π ± and K ± tracking efficiency has been quantified with calibration samples, as a function of the transverse momentum of the tracks [35,36]. The simulated signal and control channels kinematics is used to weight the simulation track efficiency in order to match the hadron detection efficiencies as measured in those calibration samples. The ΔA CP values (as measured by the difference of signal yields) are then corrected for these efficiencies and the uncertainty on the detection efficiency determination itself is propagated as a systematic uncertainty to the final ΔA CP measurements, taking into account the correlation between signal and control channel induced by the use of the same calibration samples. The systematic uncertainty arises from the size of the simulated samples used in the weighting, the statistical and systematic uncertainties on the charge asymmetry of the data calibration samples and the knowledge of the kinematical distributions generated in the simulated samples. The latter is determined by taking two different kinematic configurations of the final state (saturated by quasi two-body modes on one hand and phase-space decay on the other hand) and using the difference as the systematic estimate. The difference between the p and p particles is not measured decays, respectively. Structures corresponding to the K * (892) 0 , φ(1020) and ρ(770) 0 resonances are visible. The red lines correspond to the invariant-mass requirements applied to the selection of the main quasi-two-or three-body decays analysed to date. Simulation is used to obtain the reconstruction efficiencies as a function of the momentum of the proton or antiproton track. An additional systematic uncertainty related to the knowledge of the material budget in the simulation is added, as reported in Ref. [37]. The proton detection correction follows the same procedure as π ± and K ± detection asymmetry correction.
• The same methodology is used to correct for the difference of triggering efficiency between oppositely charged hadrons of the signal candidate, at the hardware stage of the trigger system. The trigger asymmetry effects are  Table 3 Systematic uncertainties for each decay mode. The uncertainties related to the kaon and proton detection asymmetry, the difference of triggering efficiency, the PID asymmetries and the production asymmetry are respectively reported as σ K , σ p , σ L0 , σ PID and σ A P Decay mode Absolute uncertainties (%) quantified as a function of the transverse momentum of the tracks of interest, by studying the triggering efficiency of K − and π + from the decay D 0 → K − π + [27] and protons from Λ 0 b → (Λ + c → pK − π + )π − decays. • The production asymmetry can depend on the kinematical properties of the reconstructed X 0 b candidates, though the actual dependence has not been observed yet [38]. Differences between signal and control channel X 0 b candidates kinematics would reflect in an incomplete cancel-lation of the production asymmetry in the ΔA CP observable. This effect has been estimated by considering the Λ 0 b production asymmetry measured in Ref. [38] as a function of its p T and pseudorapidity.
• The PID requirements set on the tracks of the final state can induce asymmetries. Efficiencies for the final-state particles are determined from Λ + c decays selected in data, and are parameterised by their momentum and electric charge. The correction factors to apply to the value of ΔA CP are here again determined by performing a weighting of the simulated signal and control channel events to match the efficiencies measured in the data. The uncertainties coming from the finite size of the calibration samples are propagated as a systematic uncertainty for the final ΔA CP measurements.
The first three corrections on the value of ΔA CP are found to be at the few per mille level, commensurate with their uncertainties. The lattermost source is dominating the systematic uncertainty budget, and can reach the percent level. The correction factors are however consistent with zero. The design of the fit model and the simultaneous fit strategy allow the direct measurement of the combinatorial background and the B-meson decay asymmetries. No significant asymmetries are observed and the results are presented in Sect. 7. Systematic uncertainties can be induced by the fit model and the fit complexity and it is evaluated by means of pseudoexperiments reproducing the nominal fit results. No significant biases are obtained. The normalised residuals of the signal yields are computed and the uncertainties on their pull mean value are propagated as a systematic uncertainty to each relevant ΔA CP measurement. The largest uncertainty is determined to be at the level of few 10 −4 , hence negligible in comparison to the aforementioned systematic uncertainty estimate.

Fit results
The results of the simultaneous fits to the five experimental spectra split by year of data taking, magnet polarity and trigger conditions are discussed in this section. The fit results are reported for each final state in the following subsections, and the summary of the measured yields is reported in Table 4.
• pπ − π + π − final state: Figures 2 and 3 show the results of the simultaneous fits to the invariant-mass spectra of the pπ − π + π − spectra for the inclusive, LBM and quasi two-body measurements. The high-mass region of the pπ − π + π − spectrum is only populated by either Bmeson decays or combinatorial background. The good agreement between the data and the fit model, especially in this region, validates the chosen modelling of these components. The same comment is in order for the fit in the different phase-space regions. The combinatorial component becomes negligible in the quasi two-body case. • pK − π + π − final state: Figures 4 and 5 show the results of the simultaneous fits to the pK − π + π − mass spectrum for the inclusive, LBM and quasi two-body measurements. The fit model provides also in this case a sat- isfactory description of the data, despite the very different background contributions depending on the phase-space selection. Raw asymmetries at the level of several percent are observed. • pK − K + K − final state: Figure 6 shows the results of the simultaneous fits to the reconstructed pK − K + K − mass spectrum for the inclusive, LBM and quasi two-body measurements. Negligible raw asymmetries are obtained. • pK − K + π − and pK − π + K − final states: The simultaneous fit results for the two remaining final states are shown in Fig. 7. The result of the fit for the control channel Ξ 0 b → (Ξ + c → pK − π + )π − is also displayed and shows a good description of the spectrum. This control channel is used to account for the production asymmetry of the Ξ 0 b modes.

Measurements of CP asymmetries and concluding remarks
Five charmless final states of Λ 0 b and Ξ 0 b four-body hadronic decays are examined in this paper. Specific regions of their phase space have been selected to search for local CP asym-   b → five-body legend includes two components: the partially reconstructed Λ 0 b → pπ − η and Λ 0 b → pπ − π + π − π 0 decays where a γ or π 0 is not reconstructed. The latter has a lower-mass endpoint metries in addition to the integrated CP-asymmetry. A total of eighteen measurements of CP asymmetries are reported in this paper.

LHCb
A simple counting experiment allows the measurement of a CP asymmetry up to the corrections due to instrumental and b-baryon production asymmetries. These corrections are mitigated by establishing the differences (denoted ΔA CP ) between the raw A CP values of the signals and those of the decay modes with intermediate charmed baryons comprising the same final-state particles. The asymmetries ΔA CP are further corrected for residual experimental charge asymmetries due to kinematic differences between signal and con-trol modes. The integrated ΔA CP asymmetry differences are measured to be   b → five-body legend includes two components: the partially reconstructed Λ 0 b → pπ − η and Λ 0 b → pπ − π + π − π 0 decays where a γ or π 0 is not reconstructed. The latter has a lower-mass endpoint

5-body
The measurements for the two-body low invariant-mass regions are Finally, the measurements for the quasi two-body decays are In all cases the first uncertainties are statistical and the second systematic. No significant CP violation is observed. The ΔA CP measurements for the independent samples of the two   Invariant pK − π + π − mass distributions with the results of the fit superimposed: (first row) full phase space, (second row) LBM and (third row) Λ 0 b → (Λ + c → pK − π + )π − control channel. The two columns correspond to the charge-conjugate final states: (left) baryon, (right) antibaryon. The different components employed in the fit are indicated in the legends. The Λ 0 b → five-body legend includes two components: the partially reconstructed Λ 0 b → pK − η and Λ 0 b → pK − π + π − π 0 decays where a γ or π 0 is not reconstructed. The latter has a lower-mass endpoint magnet polarities, the two categories of trigger requirements and the two distinct data-taking samples are found to be consistent. In addition, the measured asymmetries for the combinatorial background in all spectra are consistent with zero. The background contributions coming from B-meson decays (that could potentially exhibit nonzero CP violation) are also consistent with null asymmetries.
In a previous analysis, the LHCb collaboration reported evidence for CP violation in a specific region of the phase space of the decay Λ 0 b → pπ − π + π − , by measuring tripleproduct asymmetries [5]. By contrast, in the present analysis, no indication of a significant CP-violating asymmetry is obtained with the same data sample, providing complementary insights about the origin of this potential CP-symmetry breaking effect. The quest for the first observation of CP violation in baryon decays continues. LHCb Run 2 data provides about five times larger yields allowing for a more sensitive search of smaller CP-violating effects.

LHCb
Data Fit

LHCb
Data Fit   Invariant (first row) pK − K + π − and (second row) pK − π + K − mass distributions, with the results of the fit superimposed. The two bottom plots are the results of the fit to the Ξ 0 b → (Ξ + c → pK − π + )π − control channel. The two columns correspond to the charge-conjugate final states: (left) baryon, (right) antibaryon. The different components employed in the fit are indicated in the legends. The Λ 0 b → five-body legend includes two components where a π 0 is not reconstructed: the partially reconstructed background Λ 0 b → pK − π + π − π 0 where a pion is misidentified as a kaon and the partially reconstructed background Λ 0 b → pK − K + π − π 0 the plots in this publication as well as in supplementary materials are