University of Birmingham Search for CP violation in b0 pK and b0 p decays

A search for CP violation in Λ 0 b → pK − and Λ 0 b → p π − decays is presented using a sample of pp collisions collected with the LHCb detector and corresponding to an integrated luminosity of 3.0fb − 1 . The CP -violating asymmetries are measured to be A pK − CP = − 0 . 020 ± 0 . 013 ± 0 . 019 and A p π − CP = − 0 . 035 ± 0 . 017 ± 0 . 020, and their difference A pK − CP − A p π − CP = 0 . 014 ± 0 . 022 ± 0 . 010, where the ﬁrst uncertainties are statistical and the second systematic. These are the most precise measurements of such asymmetries to date. by 3 .


Introduction
The non-invariance of weak interactions under the combined application of charge conjugation (C ) and parity ( P ) transformations is accommodated within the Standard Model by the Cabibbo-Kobayashi-Maskawa mechanism [1,2]. The violation of the CP symmetry was discovered in neutral-kaon decays [3], and later observed with B 0 [4][5][6][7][8][9][10][11][12], B + [13] and B 0 s mesons [12,14]. First evidence for CP violation in the b-baryon sector was found more recently [15]. The decays Λ 0 b → pK − and Λ 0 b → pπ − are mediated by the same quark-level transitions contributing to charmless two-body B 0 and B 0 s decays to charged pions and kaons, where nonzero values of the CP asymmetries are well established [14]. The inclusion of charge-conjugate processes is implied throughout.
Predictions for the CP asymmetries in the decays of the Λ 0 b baryon to two-body charmless final states pK − or pπ − range from a few percent in the generalised factorisation approach [16,17] up to approximately 30% within the perturbative quantumchromodynamics formalism [18]. The only measurements of these quantities available to date were performed by the CDF Collaboration [12]. The asymmetries were found to be compatible with zero within an uncertainty of 8 to 9%. This Letter reports on a search for CP violation in Λ 0 b → pK − and Λ 0 b → pπ − decays, using pp-collision data collected with the LHCb detector at centre-of-mass energies of 7 and 8 TeV and corresponding to 3.0 fb −1 of integrated luminosity. The CP asymmetry is defined as where is the partial width of the given decay, with f ≡ pK − (pπ − ) and f ≡ pK + (pπ + ). In addition, the difference of the two CP asymmetries, CP , is also reported. As the main systematic uncertainties cancel in the difference, this quantity will become useful with the increasing size of the data sample.
The Letter is organised as follows. After a brief introduction on the detector, trigger and simulation in Sec. 2, the formalism needed to relate the physical CP asymmetries to the experimental measurements is presented in Sec. 3. Then, the event selection and the invariant-mass fit are described in Secs. 4 and 5, respectively. The determination of instrumental asymmetries and systematic uncertainties is discussed in Sec. 6. Finally, results are given and conclusions are drawn in Sec. 7.

Detector, trigger and simulation
The LHCb detector [19,20] 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 [21], 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 [22] 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 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 [23]. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers [24]. The online event selection is performed by a trigger [25], which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction.
Simulated events are used to study the modelling of various mass line shapes. In the simulation, proton-proton collisions are generated using Pythia [26] with a specific LHCb configuration [27]. Decays of hadronic particles are described by EvtGen [28], in which final-state radiation is generated using Photos [29]. The interaction of the generated particles with the detector and its response are implemented using the Geant4 toolkit [30] as described in Ref. [31].

Formalism
The CP asymmetries of Λ 0 b → pK − and Λ 0 b → pπ − decays are approximated as the sums of various experimental quantities where A f raw is the measured raw asymmetry between the yields its charge conjugate; A h D is the asymmetry between the detection efficiencies for particle h and its charge conjugate, with h = p, K − or π − ; the symbol A f PID stands for the asymmetry between the particle-identification (PID) efficiencies for the final states f and f ; The raw asymmetry is defined as where N denotes the observed signal yield for the given decay, obtained in this analysis by means of extended binned maximumlikelihood fits to the pK − and pπ − invariant-mass spectra. The proton, kaon and pion detection asymmetries are defined as where ε rec is the total efficiency to reconstruct the given particle, excluding PID. Such asymmetries are mostly due to the different interaction cross-sections of particles and antiparticles with the detector material. The kaon and pion detection asymmetries are measured using charm-meson control samples employing the procedures described in Refs. [32,33]. The kaon detection asymmetry is obtained by subtracting the raw asymmetries of the D + → K 0 S π + and D + → K − π + π + decay modes and correcting for the K 0 ( A K 0 D ) [32] and pion detection asymmetries. The latter is measured from the ratio of partially to fully reconstructed D * + → The proton detection asymmetry is obtained from simulated events.
The PID asymmetries are measured from large calibration samples and are defined as where ε f ( f ) PID is the PID efficiency for a final state f ( f ) given a set of PID requirements.
The Λ 0 b production asymmetry is defined as where σ denotes the inclusive production cross-section in the LHCb acceptance. The production asymmetry is taken as an external input, following Ref. [34].
Finally, asymmetries may arise if the hardware and software trigger used to collect data do not have the same efficiencies on oppositely charged particles. These effects are estimated through various data-driven techniques, as described in Sec. 6.

Event selection
The event selection starts with the reconstruction of b hadrons formed by two oppositely charged tracks with p T > 1 GeV/c, inconsistent with originating from any PV and required to form a common vertex. Each b-hadron candidate needs to have a transverse momentum greater than 1.2 GeV/c and an invariant mass, computed assigning the pion mass to both daughter tracks, in the range between 4.8 and 5.8 GeV/c 2 . Finally, each b-hadron candidate is required to be consistent with originating from a PV.
Particle-identification selection criteria are applied to divide the data sample into mutually exclusive subsamples corresponding to the final-state hypotheses pK − , pK + , pπ − , pπ + , K + π − , K − π + , K + K − and π + π − . The latter four combinations are selected to study the background due to two-body B decays, where one or both final-state particles are misidentified.
The event selection is further refined using a boosted decision tree (BDT) classifier [35,36] to reject combinatorial background. This algorithm combines the information from several input quantities to obtain a discriminant variable used to classify the b-hadron candidates as signal or background. The following properties of the final-state particles are used as input variables: the transverse momentum of the b-hadron decay products, the logarithms of their χ 2 IP values, where χ 2 IP is defined as the difference in the vertex-fit χ 2 of a given PV reconstructed with and without the candidate under consideration, the quality of the common vertex fit of the two tracks and the distance of closest approach between the two tracks. The BDT also exploits the following properties of the b-hadron candidate: the transverse momentum, the χ 2 IP quantity, and the logarithm of the flight distance with respect to the associated PV, defined as that with the smallest χ 2 IP with respect to the b-hadron candidate. The BDT is trained using simulated signal decays and combinatorial background events from data in the high-mass sideband.
The selection criteria on the BDT classifier and the PID variables are optimised separately for the Λ 0 b → pK − and Λ 0 b → pπ − decays.
Two different selections, denoted hereafter as S pK − and S pπ − , are aimed at obtaining the best statistical sensitivity on each of the two CP asymmetries. Common PID requirements are used for the final states containing only kaons and pions. Multiple candidates are present in less than 0.05% of the events in the final sample. Only one candidate is accepted for each event on the basis of a reproducible pseudorandom sequence.

Invariant-mass fit
For each final-state hypothesis, namely pK − , pK + , pπ − , pπ + , K + π − , K − π + , K + K − and π + π − , the invariant-mass distribution of selected candidates is modelled by an appropriate probability density function. These models are used to perform a simultaneous fit to the eight invariant-mass spectra and determine at once the yields of all two-body b-hadron decays contributing to the spectra. Three categories are considered for the background: combinatorial, due to random association of tracks; partially reconstructed, due to multibody b-hadron decays with one or more particles not reconstructed; and cross-feed, arising from other two-body b-hadron decays where one or both final-state particles are misidentified. The model used to describe each signal is obtained by convolving the sum of two Gaussian functions with common mean, accounting for mass-resolution effects, with a power-law function that accounts for final-state photon radiation effects. The powerlaw distribution is taken from analytical quantum-electrodynamics calculations [37] and the correctness of the model is checked against simulated events generated with Photos [29]. The parameter governing the tail due to final-state photon radiation effects is different for each decay mode. This model describes well the invariant-mass distributions predicted by the simulation.
The combinatorial background is modelled using exponential functions. The partially reconstructed background is parameterised using ARGUS functions [38] convolved with the sum of two Gaussian functions with zero mean values, whose relative fraction and widths are in common with the signal model. Finally, the crossfeed background is modelled using simulated two-body b-hadron decays and a kernel estimation method [39]. The cross-feed back-ground yields are set to the corresponding two-body b-hadron decay yields, determined by the simultaneous fit, multiplied by appropriate PID efficiency ratios. The efficiencies for a given PID requirement are obtained from large calibration samples of D * + → D 0 (→ K − π + )π + , Λ → pπ − and Λ + c → pK − π + decays, with the aid of simulated events in the case of protons to account for phasespace regions not covered by the calibration samples (about 20% of the protons from signal decays). The efficiencies are determined in bins of particle momentum, pseudorapidity and track multiplicity, as the performances of the RICH detectors depend on such variables. They are then averaged over the momentum and pseudorapidity distributions of the final-state particles and over the distribution of track multiplicity in selected events. After the application of the optimal BDT and PID requirements, an extended binned maximum-likelihood fit with a bin width of 5 MeV/c 2 is performed simultaneously to the eight two-body invariant-mass spectra for each of the two selections, S pK − and S pπ − . The m pK − and m pπ − invariant-mass distributions are shown in Fig. 1, with the results of the fits superimposed. The values of the raw asymmetries and of the signal yields obtained from the fits to the candidates passing the respective S pK − or S pπ − selec-

Instrumental asymmetries and systematic uncertainties
The determination of the instrumental asymmetries introduced in Eqs. (2) and (3) is crucial to obtain the CP asymmetries, as described in Sec. 3.
The kaon detection asymmetry is determined as a function of the kaon momentum, following the approach developed in Ref. [32] and subtracting A K 0 D = (0.054 ± 0.014)% [32] and the pion detection asymmetry. The momentum-dependent values are then weighted with the background-subtracted [40] momentum distribution of kaons from Λ 0 b → pK − decays to obtain A K − D = (−0.76 ± 0.23)%, where the dominant uncertainty is due to the finite size of the samples used. The pion detection asymmetry is obtained in an analogous way, adopting the approach of Ref. [33], and is determined to be A π − D = (0.13 ± 0.11)%. A different approach is followed for the proton detection asymmetry, since no measurement of this quantity is available to date. Simulated events are used to obtain the reconstruction efficiency defined as the number of reconstructed over generated decays, in bins of proton momentum. Then, according to Eq. (5), an asymmetry is defined and weights are computed from the background-subtracted [40] proton-momentum distributions of Λ 0 b → pK − and Λ 0 b → pπ − decays. The proton detection asymmetries for both decays are found to be equal, consistent with the fact that the kinematics of the protons for the two decays do not exhibit significant differences, as shown in Fig. 2. The common value is where the first uncertainty is due to the finite amount of simulated events and the second is associated to the knowledge of the material budget of the LHCb detector. The third uncertainty is due to the assumptions made on the proton and antiproton crosssections used in the computation. The PID asymmetries are calculated using calibration samples with the aid of simulation to account for the limited phase-space coverage of the protons from Λ → pπ − and Λ + c → pK − π + decays. The dominant uncertainty comes from different PID performances in data and simulation in the phase-space region where simulated events are used. This discrepancy has been studied us- Asymmetries related to different trigger efficiencies for the charge-conjugated final states may arise. The efficiency for a charged hadron to be responsible for the affirmative decision of the hardware trigger is determined as a function of transverse momentum, separately for positively and negatively charged particles, using a sample of Λ 0 b → Λ + c (→ pK − π + )π − decays. These efficiencies are used to determine the charge asymmetry introduced by the hardware trigger for the signal candidates that fire it. The charge asymmetry introduced by the hardware trigger for candidates that are retained independently of whether or not they are responsible for an affirmative hardware-trigger decision is determined studying a sample of B 0 → K + π − decays [25]. The asymmetry of the software trigger is also studied Several sources of systematic uncertainties associated with the fit model are investigated. The alternative models used to determine systematic uncertainties associated with the choices of the invariant-mass shapes consist in turn of: adding a Gaussian function to the invariant-mass resolution model used for signals and cross-feed backgrounds to account for long tails due to candidates with a poor determination of the final-state particles momenta; changing the value of the parameter governing the final-state photon radiation effects according to its uncertainty; substituting the exponential function used to model the combinatorial background with a linear function and removing the partially reconstructed background component by rejecting candidates with m pK − (m pπ − ) lower than 5.5 GeV/c 2 . When testing alternative models, 250 pseudoexperiments are generated according to the baseline fit model and using as input the central values of the baseline results. Fits are performed to each of the generated samples using the baseline model and then the alternative models. The mean and the root mean square of the distribution of the difference between the raw asymmetries determined by the two sets of fits are added in quadrature and the resulting value is taken as a systematic uncertainty.
A different approach is adopted to assess systematic uncertainties related to the knowledge of the PID efficiencies. Samples are generated using the baseline fit model and results. The baseline fit model is then fitted 250 times to the generated samples, varying the PID efficiencies according to their uncertainties, which are mainly driven by the choice of the binning scheme used to divide the phase-space. These uncertainties are assessed by changing the baseline binning scheme with alternative schemes and computing again the efficiencies. The largest root mean square of the raw asymmetry distributions is taken as a systematic uncertainty.
The systematic uncertainties due to the fit model choice, PID efficiencies determination and instrumental asymmetries measurement, along with the total uncertainty obtained as the quadratic sum of the individual contributions, are reported in Table 1.

Results and conclusions
Using in Eqs. (2) and (3) the values of the raw asymmetries reported in Sec. 5 and those of the instrumental and production asymmetries reported in Sec. 6, the following CP asymmetries are obtained A pK − CP = −0.020 ± 0.013 ± 0.019, A pπ − CP = −0.035 ± 0.017 ± 0.020, where the first uncertainties are statistical and the second systematic. The correlation between A pK − CP and A pπ − CP is found to be 0.5. No evidence for CP violation is observed.
A quantity that is independent from the proton detection and Λ 0 b production asymmetries is obtained by taking the difference The statistical and systematic correlations between the raw asymmetries, the PID asymmetries and the detection asymmetries are taken into account when propagating the uncertainty to A CP , obtaining A CP = 0.014 ± 0.022 ± 0.010, where the first uncertainty is statistical and the second systematic. These results represent the world's best measurements to date, with much improved precision with respect to previous CDF determinations [12].