Studies of the decays B+→pp¯h+ and observation of B+→Λ¯(1520)p

: Dynamics and direct CP violation in three-body charmless decays of charged B mesons to a proton, an antiproton and a light meson (pion or kaon) are studied using data, corresponding to an integrated luminosity of 1.0 fb−1, collected by the LHCb experiment in pp collisions at a center-of-mass energy of 7 TeV. Production spectra are determined as a function of Dalitz-plot and helicity variables. The forward-backward asymmetry of the light meson in the pp¯ rest frame is measured. No significant CP asymmetry in B+→pp¯K+ decay is found in any region of the Dalitz plane. We present the first observation of the decay B+→Λ¯¯(1520)(→K+p¯)p near the K+p¯ threshold and measure B(B+→Λ¯¯(1520)p)=(3.9+1.0−0.9(stat)±0.1(syst)±0.3(BF))×10−7, where BF denotes the uncertainty on secondary branching fractions. Dynamics and direct CP violation in three-body charmless decays of charged B mesons to a proton, an antiproton and a light meson (pion or kaon) are studied using data, corresponding to an integrated luminosity of 1 : 0 fb (cid:1) 1 , collected by the LHCb experiment in pp collisions at a center-of-mass energy of 7 TeV. Production spectra are determined as a function of Dalitz-plot and helicity variables. The forward-backward asymmetry of the light meson in the p (cid:1) p rest frame is measured. No signiﬁcant CP asymmetry in B þ ! p (cid:1) pK þ decay is found in any region of the Dalitz plane. We present the ﬁrst observation of the decay B þ ! (cid:1)(cid:2) ð 1520 Þð! K þ (cid:1) p Þ p near the K þ (cid:1) p threshold and measure B ð B þ ! (cid:1)(cid:2) ð 1520 Þ p Þ ¼ ð 3 : 9 þ 1 : 0 (cid:1) 0 : 9 ð stat Þ (cid:2) 0 : 1 ð syst Þ (cid:2) 0 : 3 ð BF ÞÞ (cid:3) 10 (cid:1) 7 , where BF denotes the uncertainty on secondary branching fractions.

Studies of the decays B þ ! p " ph þ and observation of B þ ! " Ãð1520Þp R. Aaij et al. * (LHCb Collaboration) (Received 23 July 2013;published 19 September 2013) Dynamics and direct CP violation in three-body charmless decays of charged B mesons to a proton, an antiproton and a light meson (pion or kaon) are studied using data, corresponding to an integrated luminosity of 1:0 fb À1 , collected by the LHCb experiment in pp collisions at a center-of-mass energy of 7 TeV. Production spectra are determined as a function of Dalitz-plot and helicity variables. The forwardbackward asymmetry of the light meson in the p " p rest frame is measured. No significant CP asymmetry in B þ ! p " pK þ decay is found in any region of the Dalitz plane. We present the first observation of the decay B þ ! " Ãð1520Þð! K þ " pÞp near the K þ " p threshold and measure BðB þ ! " Ãð1520ÞpÞ ¼ ð3:9 þ1:0 À0:9 ðstatÞ AE 0:1ðsystÞ AE 0:3ðBFÞÞ Â 10 À7 , where BF denotes the uncertainty on secondary branching fractions. DOI: 10.1103/PhysRevD.88.052015 PACS numbers: 13.25.Hw

I. INTRODUCTION
Evidence of inclusive direct CP violation in three-body charmless decays of B þ mesons 1 has recently been found in the modes B þ ! K þ þ À , B þ ! K þ K þ K À , B þ ! þ þ À , and B þ ! K þ K À þ [1,2]. In addition, very large CP asymmetries were observed in the low K þ K À and þ À mass regions, without clear connection to a resonance. The localization of the asymmetries and the correlation of the CP violation between the decays suggest that þ À $ K þ K À rescattering may play an important role in the generation of the strong phase difference needed for such a violation to occur [3,4]. Conservation of CPT symmetry imposes a constraint on the sum of the rates of final states with the same flavor quantum numbers, providing the possibility of entangled long-range effects contributing to the CP violating mechanism [5]. In contrast, h þ h À $ p " p (h ¼ or K throughout the paper) rescattering is expected to be suppressed compared to þ À $ K þ K À , and thus is not expected to play an important role.
The leading quark-level diagrams for the modes B þ ! p " ph þ are shown in Fig. 1. The B þ ! p " pK þ mode is expected to be dominated by the b ! s loop (penguin) transition while the mode B þ ! p " pp þ is likely to be dominated by the b ! u tree decay, which is Cabibbo-Kobayashi-Maskawa matrix suppressed compared to the former. Since the short distance dynamics are similar to that of the B þ ! h þ h þ h À modes, a CP analysis of B þ ! p " ph þ decays could help to clarify the role of long-range scatterings in the CP asymmetries of B þ ! h þ h þ h À decays. First studies were performed at the B factories on the production and dynamics of B þ ! p " ph þ decays [6][7][8]. The results have shown a puzzling opposite behavior of B þ ! p " pK þ and B þ ! p " p þ decays in the asymmetric occupation of the Dalitz plane. Charmonium contributions to the B þ ! p " pK þ decay have been studied by LHCb [9]. This paper reports a detailed study of the dynamics of the B þ ! p " ph þ decays and a systematic search for CP violation, both inclusively and in regions of the Dalitz plane. The charmless region, defined for the invariant mass m p " p < 2:85 GeV=c 2 , is of particular interest. The relevant observables are the differential production spectra of Dalitz-plot variables and the global charge asymmetry A CP , defined as pÞK þ serves as a control channel. The first observation of the decay B þ ! " Ãð1520Þp is presented. Its branching fraction is derived through the ratio of its yield to the measured yield of the B þ ! J=c ð! p " pÞK þ decay.

II. DETECTOR AND SOFTWARE
The LHCb detector [10] 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. The combined tracking system has momentum resolution Áp=p that varies from 0.4% at 5 GeV=c to 0.6% at 100 GeV=c, and impact parameter (IP) resolution of 20 m for tracks with high transverse momentum. Charged hadrons are identified using two ring-imaging Cherenkov detectors (RICH) [11]. Photon, electron and hadron candidates 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.
The trigger [12] 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. Events triggered both on objects independent of the signal, and associated with the signal, are used. In the latter case, the transverse energy of the hadronic cluster is required to be at least 3.5 GeV. The software trigger requires a two-, three-or four-track secondary vertex with a large sum of the transverse momentum, p T , of the tracks and a significant displacement from all primary pp interaction vertices. At least one track must have p T > 1:7 GeV=c, track fit 2 per degree of freedom less than 2, and an impact parameter 2 ( 2 IP ) with respect to any primary interaction greater than 16. The 2 IP is defined as the difference between the 2 of the primary vertex reconstructed with and without the considered track. A multivariate algorithm is used to identify secondary vertices [13].
The simulated pp collisions are generated using PYTHIA 6.4 [14] with a specific LHCb configuration [15]. Decays of hadronic particles are described by EVTGEN [16] in which final state radiation is generated using PHOTOS [17]. The interaction of the generated particles with the detector and its response are implemented using the GEANT4 toolkit [18] as described in Ref. [19]. Nonresonant B þ ! p " ph þ events are simulated, uniformly distributed in phase space, to study the variation of efficiencies across the Dalitz plane, as well as resonant samples such as

III. SIGNAL RECONSTRUCTION AND DETERMINATION
Candidate B þ ! p " ph þ decays are formed by combining three charged tracks, with appropriate mass assignments. The tracks are required to satisfy track fit quality criteria and a set of loose selection requirements on their momenta, transverse momenta, 2 IP , and distance of closest approach between any pair of tracks. The requirement on the momentum of the proton candidates, p > 3 GeV=c, is larger than for the kaon and pion candidates, p > 1:5 GeV=c. The B þ candidates formed by the combinations are required to have p T > 1:7 GeV=c and 2 IP < 10. The distance between the decay vertex and the primary vertex is required to be greater than 3 mm, and the vector formed by the primary and decay vertices must align with the B þ candidate momentum. Particle identification (PID) is applied to the proton, kaon and pion candidates, using combined subdetector information, the main separation power being provided by the RICH system. The PID efficiencies are derived from data calibration samples of kinematically identified pions, kaons and protons originating from the decays D Ãþ ! D 0 ð! K À þ Þ þ and Ã ! p À .
Signal and background are extracted using unbinned extended maximum likelihood fits to the mass of the p " ph þ combinations. The B þ ! p " pK þ signal is modeled by a double Gaussian function. The combinatorial background is represented by a second-order polynomial function. A Gaussian function accounting for a partially reconstructed component from B ! p " pK Ã decays is used. A possible p " p þ cross-feed contribution is included in the fit and is found to be small. An asymmetric Gaussian function with power law tails is used to estimate the uncertainties related to the variation of the signal yield.
In the case of the B þ ! p " p þ decay, the signal yield is smaller and the background is larger. The ranges of the signal and cross-feed parameters are constrained to the values obtained in the simulation within their uncertainties. The signal and the p " pK þ cross-feed contribution are modeled with Gaussian functions. The combinatorial background is represented by a third-order polynomial function.
The B þ ! p " ph þ invariant mass spectra are shown in Fig. 2 To probe the dynamics of the B þ ! p " ph þ decays, differential production spectra are derived as a function of m p " p and cos p , where p is the angle between the charged meson h and the opposite-sign baryon in the rest frame of the p " p system. The p " ph þ invariant mass is fitted in bins of the aforementioned variables and the signal yields are corrected for trigger, reconstruction and selection efficiencies. They are estimated with simulated samples and corrected to account for discrepancies between data and simulation. The signal yields are determined with the fit models described in the previous section, but allowing the combinatorial background parameters to vary. The systematic uncertainties are determined for each bin and include uncertainties related to the PID correction, fit model, trigger efficiency, and the size of the simulated samples. The latter is evaluated from the differences between data and simulation as a function of the Dalitz-plot variables. No trigger-induced distortions are found.
A. Invariant mass of the p " p system The yields and total efficiency for B þ ! p " ph þ in m p " p bins are shown in Tables I and II. The charmonium contributions originate from the decays pK þ mode, and B þ ! J=c ð! p " pÞ þ for the B þ ! p " p þ mode. Before deriving the distributions, the charmonium contributions are unfolded by p þ signal yield, including the J=c mode, B þ ! p " pK þ cross-feed yield, signal efficiency, and relative systematic uncertainty in bins of p " p invariant mass. pK þ signal yield, including the charmonium modes, efficiency and relative systematic uncertainty, in bins of p " p invariant mass. The error on the efficiency includes all the sources of uncertainty.  . Invariant mass distributions of (left) p " pK þ and (right) p " p þ candidates. The points with error bars represent data. The solid black line represents the total fit function. Blue dashed, purple dotted, red long-dashed and green dashed-dotted curves represent the signal, cross-feed, combinatorial background and partially reconstructed background, respectively. D 88, 052015 (2013) performing two dimensional extended unbinned maximum likelihood fits to the p " ph þ and p " p invariant masses. The J=c and c ð2SÞ resonances are modeled by Gaussian functions and the c resonance is modeled by a convolution of Breit-Wigner and Gaussian functions. The nonresonant p " p component and the combinatorial background are modeled by polynomial shapes. Table III shows the yields of contributing charmonium modes. The results are consistent with those reported in Ref. [9].

STUDIES OF THE DECAYS
After unfolding, the efficiency-corrected differential distributions are shown in Fig. 3. An enhancement is observed at low p " p mass both for B þ ! p " pK þ and B þ ! p " p þ , with a more sharply peaked distribution for B þ ! p " p þ . This accumulation of events at low m p " p is a well known feature that has also been observed in different contexts such as Çð1SÞ ! p " p [20], J=c ! p " p [21] and B 0 ! D ðÃÞ0 p " p [22] decays. It appears to be caused by proton-antiproton rescattering and is modulated by the particular kinematics of the decay from which the p " p pair originates [23].

B. Invariant mass squared of the Kp system
The B þ ! p " pK þ signal yield as a function of the Dalitz-plot variable m 2 Kp is considered, where Kp denotes the neutral combinations K À p or K þ " p. Table IV shows the yields and efficiencies, after the charmonium bands have been vetoed in the ranges m p " p 2 ½2:85; 3:15 GeV=c 2 and ½3:60; 3:75 GeV=c 2 . The differential spectrum derived after efficiency correction is shown in Fig. 4. Contrary to the situation for m p " p , the data distribution is in reasonable agreement with the uniform phase space distribution, with some discrepancies in the region m 2 Kp 2 ½4; 12 ðGeV=c 2 Þ 2 .
C. Helicity angle of the p " p system The B þ ! p " ph þ signal yields are considered as a function of cos p . Tables V and VI show the corresponding yields and efficiencies. The differential distributions are shown in Fig. 5.
A clear opposite angular correlation between B þ ! p " pK þ and B þ ! p " p þ decays is observed; the light meson h tends to align with the opposite-sign baryon for B AE ! p " pK AE while it aligns with the same-sign baryon for the B AE ! p " p AE mode. A quark level analysis suggests that the meson should align with the same-sign baryon, since the opposite-sign baryon has larger momentum, being formed by products from the decaying quark [24]. This is in agreement with the angular spectrum of B þ ! p " p þ but not for B þ ! p " pK þ decays.

D. Dalitz plot
From the fits to the B-candidate invariant mass, shown in Fig. 2, signal weights are calculated with the sPlot technique [25] and are used to produce the signal Dalitz-plot distributions shown in Fig. 6. To ease the comparison, the cos p curves corresponding to the boundaries of the eight bins used to make the angular distributions in Fig. 5 are superimposed.
With the exception of the charmonium bands [ c , J=c , c ð2SÞ for B þ ! p " pK þ , and J=c for B þ ! p " p þ ], the structure of the low p " p mass enhancement is very different between B þ ! p " pK þ and B þ ! p " p þ . The B þ ! p " pK þ events are distributed in the middle and lower m 2 Kp half, exhibiting a possible p " p band structure near 4 GeV 2 =c 4 . An enhancement at low m Kp is also observed and is caused to a large extent by a Ãð1520Þ signal, as will be shown in the next section. The B þ ! p " p þ events are mainly clustered in the upper m 2 p half, with also a few events   4 (color online). Efficiency-corrected differential yield as a function of m 2 Kp for B þ ! p " pK þ . The data points are shown with their statistical and total uncertainties. The solid line represents the expectation for a uniform phase space production, normalized to the efficiency-corrected area, for comparison. D 88, 052015 (2013) on the doubly charged top diagonal ðpÞ þþ (near the cos p ¼ À1 boundary). These distributions of events are consistent with the angular distributions and asymmetries reported earlier.

V. MEASUREMENT OF
The raw charge asymmetry is obtained by performing a simultaneous extended unbinned maximum likelihood fit to the B À and B þ samples. The B AE yields are defined as a function of the total yield N and the raw asymmetry, A raw , by The CP asymmetry is then derived after correcting for the B AE production asymmetry A P ðB AE Þ and the kaon detection asymmetry A D ðK AE Þ The correction A Á ¼ A P ðB AE Þ þ A D ðK AE Þ is measured from data with the decay B AE ! J=c ð! p " pÞK AE which is part of the data sample pÞK AE Þ À A CP ðJ=c K AE Þ; (4) where A CP ðJ=c K AE Þ ¼ ð1 AE 7Þ Â 10 À3 [26].
Another correction has been applied to account for the proton-antiproton asymmetry, which exactly cancels for J=c ð! p " pÞK AE but not necessarily in the full phase space of p " pK AE events. This effect has been estimated in simulation studying the difference in the interactions of protons and antiprotons with the detector material between J=c ð! p " pÞK AE and p " pK AE events generated uniformly over phase space. We obtained a m 2 Kp -dependent bias, up to 3% for the highest bin, for A raw .
To measure A raw for charmonium modes, and in particular J=c ð! p " pÞK AE , a two-dimensional ðm B ; m p " p Þ simultaneous fit to the B þ and B À samples is performed. The systematic uncertainties are estimated by varying the fit functions and splitting the data sample according to trigger requirements or magnet polarities, and recombining the results from the subsamples. The procedure is applied to obtain a global value of A CP as well as the variation of the asymmetry as a function of the Dalitz-plot variables. The results are A CP ¼ À0:022 AE 0:031ðstatÞ AE 0:007ðsystÞ for the full p " pK AE spectrum, and A CP ¼ À0:047 AE 0:036ðstatÞ AE 0:007ðsystÞ for the region m p " p < 2:85 GeV=c 2 . Figure 7 shows the variation of A CP as a function of the Dalitz-plot variables.

VI. OBSERVATION OF THE
In the p " pK þ spectrum, near the threshold of the neutral Kp combination, a peak in invariant mass at 1:52 GeV=c 2 is observed, as shown in Fig. 8, corresponding to the " u " d " s resonance " Ãð1520Þ. The possible presence of higher " Ã and " AE resonances may explain the enhancement in the range of ½1:6; 1:7 GeV=c 2 .
To identify the " Ãð1520Þ signal, the B þ signal is analyzed in the region m Kp 2 ½1:44; 1:585 GeV=c 2 . Figure 9 shows the B signal weighted Kp invariant mass, and the expected " Ãð1520Þ shape obtained from a model based on an asymmetric Breit-Wigner function derived from an EVTGEN [16] simulation of the decay B þ ! " Ãð1520Þp, convolved with a Gaussian resolution function, and a second-order polynomial function representing the tail of the non-" Ãð1520Þ B þ ! p " pK þ decays. These shapes are then used in a two-dimensional ðm p " pK þ ; m Kp Þ extended unbinned maximum likelihood fit to obtain the B þ ! " Ãð1520Þp yield. The fit results in NðB þ ! " Ãð1520ÞpÞ ¼ 47 þ12 À11 with a statistical significance of 5.3 standard deviations, obtained by comparing the likelihood at its maximum for the nominal fit and for the background-only hypothesis. Figure 10 shows the projections of the fit for the Kp and p " pK þ invariant masses. To test the robustness of the observation, different representations of the Kp background have been used, combining first-or second-order polynomials and a contribution modeled by a Breit-Wigner function, for which the mean () and width (À) are allowed to vary within the known values of the Ãð1600Þ baryon ( 2 ½1:56; 1:7 GeV=c 2 , À 2 ½0:05; 0:25 GeV=c 2 ). Fits in a wider m Kp range were also considered. In all cases the yield was stable with a statistical significance similar to the nominal fit case.
The branching fraction for the decay B þ ! " Ãð1520Þp is derived from the ratio    D 88, 052015 (2013) where N i is the yield of the decay chain i, and gen denotes the efficiency after geometrical acceptance and simulation requirements. The global selection efficiency sel includes the reconstruction, the trigger, the offline selection, and the particle identification requirements. The ratio of branching fractions obtained is

VII. SUMMARY
Based on a data sample, corresponding to an integrated luminosity of 1:0 fb À1 , collected in 2011 by the LHCb experiment, an analysis of the three-body B þ ! p " ph þ decays (h ¼ K or ) has been performed. The dynamics of the decays has been probed using differential spectra of Dalitz-plot variables and signal-weighted Dalitz plots. The charmless B þ ! p " pK þ decay populates mainly the low m 2 p " p and lower m 2 K þ " p -half regions whereas the B þ ! p " p þ decay has a similar enhancement at low m 2 p " p but with an upper m 2 þ " p -half occupancy. From the occupation pattern of the Dalitz plots, it is likely that the B þ ! p " pK þ decay is primarily driven by p " p rescattering with a secondary contribution from neutral Kp rescattering while the B þ ! p " p þ decay is also dominated by p " p rescattering but with a secondary contribution from doubly charged ðpÞ þþ rescattering, along the lines of the rescattering amplitude analysis performed in Ref. [28]. This difference of behavior is reflected in the values of the forward-backward asymmetry of the light meson in the p " p rest frame pK þ Þ ¼ 0:370 AE 0:018ðstatÞ AE 0:016ðsystÞ; A FB ðp " p þ Þ ¼ À0:392 AE 0:117ðstatÞ AE 0:015ðsystÞ: CP asymmetries for the B þ ! p " pK þ decay have been measured and no significant deviation from zero observed: A CP ¼ À0:047 AE 0:036ðstatÞ AE 0:007ðsystÞ for the charmless region m p " p < 2:85 GeV=c 2 , A CP ð c K AE Þ ¼ 0:046 AE 0:057ðstatÞ AE 0:007ðsystÞ and A CP ðc ð2SÞK AE Þ ¼ À0:002 AE 0:123ðstatÞ AE 0:012ðsystÞ. These measurements are consistent with the current known values, A CP ðB AE ! p " pK AE ; m p " p < 2:85 GeV=c 2 Þ ¼ À0:16 AE 0:07 [26], A CP ð c K AE Þ ¼ À0:16 AE 0:08ðstatÞ AE 0:02ðsystÞ [8],  and A CP ðc ð2SÞK AE Þ ¼ À0:025 AE 0:024 [26]. The absence of any significant charge asymmetry, contrary to the situation for B þ ! h þ h þ h À decays [1,2], may be due to different long-range behavior. Final state interactions in the B þ ! p " ph þ case do not change the nature of the particles, such as p " p ! p " p or ph ! ph, while B þ ! h þ h þ h À modes can be affected by þ À $ K þ K À scattering.

ACKNOWLEDGMENTS
We express our gratitude to our colleagues in the CERN accelerator departments for the excellent performance of the LHC. We thank the technical and administrative staff at the LHCb institutes. We acknowledge support from CERN and from the national agencies: CAPES, CNPq, FAPERJ and FINEP (Brazil); NSFC (China); CNRS/IN2P3 and Region Auvergne (France); BMBF, DFG, HGF and MPG (Germany); SFI (Ireland); INFN (Italy); FOM and NWO (The Netherlands); SCSR (Poland); MEN/IFA (Romania); MinES, Rosatom, RFBR and NRC ''Kurchatov Institute'' (Russia); MinECo, XuntaGal and GENCAT (Spain); SNSF and SER (Switzerland); NAS Ukraine (Ukraine); STFC (United Kingdom); NSF (USA). We also acknowledge the support received from the ERC under FP7. The Tier1 computing centers are supported by IN2P3 (France), KIT and BMBF (Germany), INFN (Italy), NWO and SURF (The Netherlands), PIC (Spain), GridPP (United Kingdom). We are thankful for the computing resources put at our disposal by Yandex LLC (Russia), as well as to the communities behind the multiple open source software packages that we depend on.  D 88, 052015 (2013)