Measurements of $C\!P$ violation in the three-body phase space of charmless $B^{\pm}$ decays

The charmless three-body decay modes $B^{\pm} \rightarrow K^{\pm} \pi^{+} \pi^{-}$, $B^{\pm} \rightarrow K^{\pm} K^{+} K^{-}$, $B^{\pm} \rightarrow \pi^{\pm} K^{+} K^{-}$ and $B^{\pm} \rightarrow \pi^{\pm} \pi^{+} \pi^{-}$ are reconstructed using data, corresponding to an integrated luminosity of 3.0\;$\mbox{\,fb}^{-1}$, collected by the LHCb detector. The inclusive $C\!P$ asymmetries of these modes are measured to be $$ A_{C\!P}(B^{\pm} \rightarrow K^{\pm} \pi^{+} \pi^{-})= +0.025 \pm 0.004 \pm 0.004 \pm 0.007 \\ A_{C\!P}(B^{\pm} \rightarrow K^{\pm} K^{+} K^{-}) = -0.036 \pm 0.004 \pm 0.002 \pm 0.007 \\ A_{C\!P}(B^{\pm} \rightarrow \pi^{\pm} \pi^{+} \pi^{-})= +0.058 \pm 0.008 \pm 0.009 \pm 0.007 \\ A_{C\!P}(B^{\pm} \rightarrow \pi^{\pm} K^{+} K^{-})= -0.123 \pm 0.017 \pm 0.012 \pm 0.007 $$ %% where the first uncertainty is statistical, the second systematic, and the third is due to the $C\!P$ asymmetry of the $B^{\pm} \rightarrow J/\psi K^{\pm}$ reference mode. The distributions of these asymmetries are also studied as functions of position in the Dalitz plot and suggest contributions from rescattering and resonance interference processes.


Introduction
The violation of CP symmetry is well established experimentally in the quark sector and, in the Standard Model (SM), is explained by the Cabibbo-Kobayashi-Maskawa [1] matrix through the presence of a single irreducible complex phase.Although the SM is able to describe all CP asymmetries observed experimentally in particle decays, the amount of CP violation within the SM is insufficient to explain the matter-antimatter asymmetry of the universe [2].
The decays of B mesons with three charged charmless mesons in the final state offer interesting opportunities to search for different sources of CP violation, through the study of the signature of these sources in the Dalitz plot.Several theoretical studies modelled the dynamics of the decays in terms of two-body intermediate states, such as ρ(770)K ± or ( ) K * 0 (892)π ± for B ± → K ± π + π − decays, and φ(1020)K ± for B ± → K ± K + K − decays (see e.g.[3]).These intermediate states were identified through amplitude analyses in which a resonant model was assumed.One method of performing such analyses was used by the Belle and the BaBar collaborations and significant CP violation was observed in the intermediate ρ 0 K ± state [4,5] and in the φK ± channel [6].No significant inclusive CP asymmetry (integrated over the Dalitz plot) was found in [4,6].Another method is to measure the CP asymmetry in different regions of the three-body phase space.The LHCb collaboration measured non-zero inclusive CP asymmetries and larger local asymmetries in the decays [8] using a sample corresponding to 1.0 fb −1 of data.These results suggested that final-state interactions may be a contributing factor to CP violation [9,10].
Direct CP violation requires the existence of amplitudes with differences in both their weak and their strong phases.The value of the weak phase can be accessed through interference between tree-level contributions to charmless B decays and other amplitudes (e.g.penguins).The strong phase can originate from three different sources in charmless three-body decays.The first source is related to short-distance processes where the gluon involved in the penguin contribution is timelike, i.e. the momentum transfer satisfies q 2 > 4m 2  i , where m i represents the mass of either the u or the c quark present in the loop diagram [11].This process is similar to that proposed for two-body decays where CP violation is caused by short-distance processes [12].The remaining two sources are related to long-distance effects involving hadron-hadron interactions in the final state.Interference between intermediate states of the decay can introduce large strong-phase differences, and therefore induce local asymmetries in the phase space [9,13,14].Another mechanism is final-state KK ↔ ππ rescattering, which can occur between decay channels having the same flavour quantum numbers [7][8][9][10].Conservation of CP T symmetry constrains hadron rescattering so that the sum of the partial decay widths of all channels with the same final-state quantum numbers related by the scattering matrix must equal that of their charge-conjugated decays [15].The effects of SU(3) flavour symmetry breaking have also been investigated and can explain part of the pattern of CP violation reported by LHCb [9,[15][16][17].
In this paper, the inclusive CP asymmetries of B ± → K ± π + π − , B ± → K ± K + K − , B ± → π ± K + K − and B ± → π ± π + π − decays (henceforth collectively referred to as B ± → h ± h + h − decays) are measured, and local asymmetries in specific regions of the phase space are studied.All asymmetries are measured using the B ± → J/ψ K ± channel, which has similar topology and negligible CP violation, as a reference, thus allowing corrections to be made for production and instrumental asymmetries.We use a sample of proton-proton collisions collected in 2011 (2012) at a centre-of-mass energy of 7 (8) TeV and corresponding to an integrated luminosity of 1.0 (2.0) fb −1 .This analysis supersedes that of [7,8], by using a larger data sample, improved particle identification and a more performant event selection.

LHCb detector and data set
The LHCb detector [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.The tracking system provides a measurement of momentum, p, with a relative uncertainty that varies from 0.4% at low momentum to 0.6% at 100 GeV/c.The minimum distance of a track to a primary vertex, the impact parameter, is measured with a resolution of (15 + 29/p T )µm, where p T is the component of p transverse to the beam, in GeV/c.Charged hadrons are identified using two ring-imaging Cherenkov (RICH) detectors [19].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 [20].
The trigger [21] consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies full event reconstruction.At the hardware trigger stage, events are required to have a muon with high p T or a hadron, photon or electron with high transverse energy in the calorimeters.For hadrons, the transverse energy threshold is 3.5 GeV.In this analysis two partially overlapping categories of events selected by the hardware trigger are considered: events where one of the hadrons from the B ± decay is used in the trigger decision (the "trigger on signal" sample), and events that are triggered by particles other than those hadrons from the B ± decay (the "trigger independent of signal" sample).
At the software trigger stage, events must have at least one good quality track from the signal decay candidate with high p T and a significant displacement from any primary vertex (PV).A secondary vertex, consisting of three good quality tracks that have significant displacements from any PV, is also required.
The magnetic field polarity is reversed regularly during the data taking to reduce any potential bias from charged particle and antiparticle detection asymmetries.The magnetic field bends charged particles in the horizontal plane and the two polarities are referred to as "up" and "down".The fraction of data collected with the magnet down polarity is approximately 60% in 2011, and 52% in 2012.
Possible residual charge-dependent asymmetries, which may originate from left-right differences in detection efficiency, are studied by comparing measurements from data with inverted magnet polarities and found to be negligible.Since the detection and production asymmetries are expected to change between 2011 and 2012 due to different data taking conditions, the analysis is carried out separately for the 2011 and 2012 data and the results combined.
The simulated events are generated using Pythia 8 [22] with a specific LHCb configuration [23].Decays of hadronic particles are produced by EvtGen [24], in which final-state radiation is generated using Photos [25].The interaction of the generated particles with the detector and its response are implemented using the Geant4 toolkit [26] as described in Ref. [27].

Event selection
Since the four B ± signal decay modes considered are topologically and kinematically similar, the same selection criteria are used for each, except for the particle identification requirements, which are specific to each final state.The decay B ± → J/ψ K ± , J/ψ → µ + µ − serves as a control channel for B ± → h ± h + h − decay modes.Since it has negligible CP violation, the raw asymmetry observed in B ± → J/ψ K ± decays is entirely due to production and detection asymmetries.The control channel has a similar topology to the signal and the sample passes the same trigger, kinematic, and kaon particle identification selection as the signal samples.The kaons from B ± → J/ψ K ± decays also have similar kinematic properties in the laboratory frame to those from the B ± → K ± π + π − and In a preselection stage, loose requirements are imposed on the p, p T and the displacement from any PV for the tracks, and on the distance of closest approach between each pair of tracks.The three tracks must form a good quality secondary vertex that has a significant separation from its associated PV.The momentum vector of the reconstructed B ± candidate has to point back to the PV.
A multivariate selection based on a boosted decision tree (BDT) algorithm [29][30][31] is applied to reduce the combinatorial background.The input variables, which are a subset of those used in the preselection, are common to all four decay modes.The BDT is trained using a mixture of simulated signal events as the signal sample, and events reconstructed as B ± → π ± π + π − decays with 5.40 < m(π ± π + π − ) < 5.58 GeV/c 2 as the background sample.
The requirement on the BDT response is chosen to maximise the ratio N S / √ N S + N B , where N S and N B represent the expected number of signal and background candidates, respectively, within an invariant mass window of approximately 40 MeV/c 2 around the signal peak.Since the optimal requirements are similar for the different channels, the same BDT response requirement is chosen for all channels, to simplify the evaluation of the systematic uncertainties.The BDT selection improves the efficiencies for selecting signal events by approximately 50%, compared to the cut-based selection used in Refs.[7] and [8].
Particle identification is used to reduce the cross-feed from other B decays in which hadrons are incorrectly classified.The main source is K → π and π → K misidentification, while p → K and p → π misidentification is negligible.Muons are rejected by a veto applied to each track [32].After the full selection, events with more than one candidate in the range 4.8 < m(B) < 5.8 GeV/c 2 are discarded.This removes approximately 1-2% of candidates.
The B ± → J/ψ K ± control channel is selected using the same criteria as described above, with two exceptions that enhance the selection of J/ψ mesons decaying to two muons: criteria used to identify charged pions are removed and the requirement 3.05 < m(π + π − ) < 3.15 GeV/c 2 is applied.

Determination of signal yields
For each channel the yields and raw asymmetry are extracted from a single simultaneous unbinned extended maximum likelihood fit to the B + and B − invariant mass distribution.The signal components of all four channels are parametrised by a Gaussian function with widths and tails that differ either side of the peak to account for asymmetric effects such as final-state radiation.The means and widths are allowed to vary in the fits, while the tail parameters are fixed to values obtained from simulation.The combinatorial backgrounds are described by exponential functions.The backgrounds due to partially reconstructed four-body B decays are parametrised by an ARGUS function [33] convolved with a Gaussian function.The shapes and yields of peaking backgrounds, i.e. fully reconstructed B decays with at least one misidentified particle in the final state, are obtained from simulation of the relevant decay modes and fixed in the fits.The yields of the peaking and partially reconstructed background components are constrained to be equal for B + and B − decays.
The invariant mass spectra of the four decay modes are shown in Fig. 1.The figure is illustrative only, as the asymmetries are obtained from separate fits of the samples divided by year, trigger selection and magnet polarity, and then combined as described in Sec. 5.
The signal yields obtained for the combined 2011 and 2012 data samples are shown in Table 1.The data samples are larger than those presented in Refs.[7] and [8] due to both an increase in the integrated luminosity and the use of a more efficient selection.
The left panel on each figure shows the B − candidates and the right panel shows the B + candidates.The results of the unbinned maximum likelihood fits are overlaid.The main components of the fits are also shown.
Table 1: Signal yields of charmless three-body B ± decays for the full data set.

Decay mode Yield
5 Inclusive CP asymmetry measurement The CP asymmetry of B ± decays to a final state f ± is defined as where Γ is the partial decay width.To determine the inclusive CP asymmetries, the raw asymmetries measured from the fits are corrected for effects induced by the detector efficiency, interactions of final-state particles with matter, and any asymmetry in the forward production rates between B + and B − mesons.The raw asymmetry, A raw , is written in terms of the B − and B + event yields as where the numbers of signal events N B − and N B + are related to the asymmetries by Here, A P is the B ± meson production asymmetry, N S are the total yields, and ε ± are the average efficiencies for selecting and reconstructing B + and B − decays, respectively.The efficiency is computed on an event-by-event basis and depends on the position in the Dalitz plot.The term A D accounts for residual detection asymmetries, such as differences in interactions of final state particles with the detector material or left-right asymmetries that may not be properly represented in the Monte Carlo.
The three final-state hadrons are treated as the combination of a pair of same-flavour, charge-conjugate hadrons h + h − = π + π − , K + K − , and an unpaired hadron h ± with the same charge as the B ± meson.The detection asymmetry A h D is given in terms of the charge-conjugate detection efficiencies of the unpaired hadron h ± , and the production asymmetry A P is given in terms of the B ± production rates.
The raw asymmetry is expressed in terms of A CP , A P and A h D using Eqs.( 2) and ( 3), For small asymmetries the products are negligible, and the raw asymmetry becomes Throughout this paper, Eq. ( 5) is used in calculating the inclusive asymmetries, as all terms are sufficiently small.For the determination of the asymmetries in regions of the phase space where the raw asymmetries are large, the full formula of Eq. ( 4) is applied.
The four decay channels are divided into two categories according to the flavour of the final-state hadron h ± .For the B ± → K ± π + π − and B ± → K ± K + K − decay channels, the CP asymmetry is expressed in terms of the raw asymmetry and correction terms given by the sum of the B ± production asymmetry and the kaon detection asymmetry, A P and A K D .For the B ± → π ± K + K − and B ± → π ± π + π − decay channels, the pion detection asymmetry A π D is used.The CP asymmetries are calculated as The correction term A ∆ is measured using approximately 265 000 B ± → J/ψ (µ + µ − )K ± decays.The correction is obtained from the raw asymmetry of the B ± → J/ψ K ± mode as using the world average of the CP asymmetry A CP (J/ψ K ± ) = (0.1 ± 0.7)% [28].
The pion detection asymmetry, A π D = (0.00 ± 0.25)%, has been previously measured by LHCb [34] and is consistent with being independent of p and p T .The production asymmetry is obtained from the same sample of B ± → J/ψ K ± decays as A P = A ∆ − A K D , and is consistent with being constant in the interval of momentum measured.Here the kaon interaction asymmetry A K D = (−1.26± 0.18)% is measured in a sample of , where the D * + is produced in the decay of a B meson.The value of A K D is obtained by measuring the ratio of fully to partially reconstructed D * + decays [34].
Since neither the detector efficiencies nor the observed raw asymmetries are uniform across the Dalitz plot, an acceptance correction is applied to the integrated raw asymmetries.This is determined by the ratio of the B − and B + average efficiencies in simulated events, reweighted to reproduce the population of signal data in bins of the Dalitz plot.In addition, to account for the small charge asymmetry introduced by the hadronic hardware trigger, the data are divided into the trigger independent of signal and the trigger on signal samples, as discussed in Sec. 2. The CP asymmetries are calculated using Eqs.( 6) and ( 7), applied to the acceptance-corrected raw asymmetries of the samples collected in each trigger configuration.The inclusive CP asymmetry of each mode is the weighted average of the CP asymmetries for the samples divided by trigger and year of data taking, taking into account the correlation between trigger samples as described in Ref. [35].

Systematic uncertainties and results
Several sources of systematic uncertainty are considered.These include potential mismodellings in the mass fits, the phase-space acceptance corrections and the trigger composition of the samples.
The systematic uncertainties due to the mass fit models are evaluated as the full difference in CP asymmetry resulting from variations of the model.The alternative fits have good quality and describe the data accurately.To estimate the uncertainty due to the choice of the signal mass function, the initial model is replaced by an alternative empirical distribution [36].A systematic uncertainty to account for the use of equal means and widths for B − and B + signal peaks in the default fit is assigned by repeating the fits with these parameters allowed to vary independently.The resulting means and widths are found to agree and the difference in the value of A CP is assigned as a systematic uncertainty.
The systematic uncertainty associated with the peaking background fractions reflects the uncertainties in the expected yields determined from simulation, and the influence of combining 2011 and 2012 simulated samples when determining the fractions in the nominal fit, by repeating the fits with the background fractions obtained for the samples separately.The uncertainty due to background shape is obtained by increasing the width of the Gaussian function according to the observed differences between simulation and data for peaking backgrounds, and allowing the four-body shape to vary in the fit.Similarly, the possibility of non-zero background asymmetries is tested by letting the peaking and four-body-background normalisations vary separately for B − and B + fits.The signal model variations and the background asymmetry are the dominant systematic uncertainties related to the fit procedure.
The systematic uncertainty related to the acceptance correction procedure consists of two parts: the statistical uncertainty on the detection efficiency due to the finite size of the simulated samples, and the uncertainty due to the choice of binning, which is evaluated by varying the binning used in the efficiency correction.
A study is performed to investigate the effect of having different trigger admixtures in the signal and the control channels.The acceptance-corrected CP asymmetries are measured separately for each trigger category and found to agree, therefore no additional systematic uncertainty is assigned.Performing this comparison validates the assumption that the detection asymmetry factorises between the h + h − pair and the h ± , within the statistical precision of the test.
The systematic uncertainties, separated by year, are shown in Table 2, where the total systematic uncertainty is the sum in quadrature of the individual contributions.The uncertainties on A π D and A K D are only considered as systematic uncertainties for B ± → π ± π + π − and B ± → π ± K + K − decays, following Eq.( 6).The systematic uncertainty of the 2011 and 2012 combination is taken to be the greater of these two values.LHCb (d) The visible gaps correspond to the exclusion of the J/ψ (in the B ± → K ± π + π − decay) and D 0 (all plots, except for the B ± → π ± K + K − decay) mesons from the samples.
The results for the integrated CP asymmetries are where the first uncertainty is statistical, the second systematic, and the third is due to the limited knowledge of the CP asymmetry of the B ± → J/ψ K ± reference mode [28].
The significances of the inclusive charge asymmetries, calculated by dividing the central values by the sum in quadrature of the uncertainties, are 2.8 standard deviations (σ) for

CP asymmetry in the phase space
The Dalitz plots distributions in the signal region for the four channels are shown in Fig. 2. For the B ± → K ± K + K − and B ± → π ± π + π − decays, folded Dalitz plots are used.
For a given event, the vertical axis of the Dalitz plot corresponds to the invariant mass squared of the decay with the highest value (m 2 (h + h − ) high ), while the horizontal axis is the invariant mass squared with the lowest value between the two (m 2 (h + h − ) low ).The signal region is defined as the three-body invariant mass region within 34 MeV/c 2 of the fitted mass, except for the B ± → π ± K + K − channel, for which the mass window is restricted to ±17 MeV/c 2 of the peak due to the larger background.The expected background contribution is not subtracted from the data presented in these figures.To improve the resolution, the Dalitz variables are calculated after refitting the candidates with their invariant masses constrained to the known B ± value [28].The events are concentrated in low-mass regions, as expected for charmless decays dominated by resonant contributions.
An inspection of the distribution of candidates from the B + mass sidebands confirms that the background is not uniformly distributed, with combinatorial background events tending to be concentrated at the corners of the phase space, as these are dominated by low-momentum particles.
In addition to the inclusive charge asymmetries, the asymmetries are studied in bins of the Dalitz plots.Figure 3 shows these asymmetries, A N raw ≡ N − −N + N − +N + , where N − and N + are the background-subtracted, efficiency-corrected signal yields for B − and B + decays, respectively.Background subtraction is done via a statistical tool to unfold data distributions called sPlot technique [39] using the B ± candidate invariant mass as discriminating variable.The binning is chosen adaptively, to allow approximately equal populations of the total number of entries (N − + N + ) in each bin.
The A N raw distributions in the Dalitz plots reveal rich structures, which are more evident in the two-body invariant-mass projection plots.These are shown in Figs. 4  and 5 for the region of the ρ resonance in B ± → π ± π + π − and B ± → K ± π + π − decays, respectively.The projections are split according to the sign of cos θ, where θ is the angle between the momenta of the unpaired hadron and the resonance decay product with the same-sign charge.Figure 6 shows the projection onto the low K + K − invariant mass for the B ± → K ± K + K − channel, while Fig. 7 shows the projection into m(K + K − ) for the The dynamic origin of the CP asymmetries seen in Fig. 3 can only be fully understood with an amplitude analysis of these channels.Nevertheless, the projections presented in Figs. 4, 5, 6 and 7 indicate two different sources of CP violation.The first one may be associated with the π + π − ↔ K + K − rescattering strong-phase difference in the region around 1.0 to 1.5 GeV/c 2 [7,8].In this region, there are more B − than B + decays into final states including a π + π − pair (positive CP asymmetry) and more B + than B − into final states that include a K + K − pair (negative CP asymmetry).The second source of CP violation, observed in both B ± → K ± π + π − and B ± → π ± π + π − decays around the ρ(770) mass region, can be attributed to the final-state interference between the S-wave and P-wave in the Dalitz plot.

CP asymmetry induced by rescattering
Previous publications [7,8] showed evidence for possible source of CP violation produced by the long-distance strong phase through π + π − ↔ K + K − rescattering.This interaction plays an important role in S-wave π + π − elastic scattering, as was observed by previous experiments [40,41] in the m(π + π − ) mass region between 1.0 and 1.5 GeV/c 2 .The CP T symmetry requires that the sum of partial widths of a family of final states related to each other by strong rescattering, such as the four channels analysed here, are identical for particle and antiparticle.As a consequence, positive CP asymmetry in some channels implies negative CP asymmetry in other channels of the same family.
The large data samples in the present study allow this effect to become evident, as shown in Figs. 4, 5, 6 and 7. Large asymmetries are observed for all the final states in the region between 1.0 and 1.5 GeV/c 2 .Figure 8 shows the invariant mass distributions for events with m(π + π − ) and m(K + K − ) in this interval, excluding the φ meson mass  region for the B ± → K ± K + K − mode.The measured CP asymmetries corresponding to the figure are given in Table 3. Decays involving a K + K − pair in the final state have a larger CP asymmetry than their partner channels with a π + π − pair.The asymmetries are positive for channels with a π + π − pair and negative for those with a K + K − pair.This indicates that the mechanism of π + π − ↔ K + K − rescattering could play an important role for CP violation in charmless three-body B ± decays.

CP asymmetry due to interference between partial waves
In hadronic three-body decays, there is another long-distance strong-interaction phase, which stems from the amplitudes of intermediate resonances.The Breit-Wigner propagator associated with an intermediate resonance can provide a phase that varies with the resonance mass.There is also a phase related to final-state interactions, associated with each intermediate state that contributes to the same final state.In general, the latter phase      is considered constant within the phase space.This phase also includes any short-distance strong phase.These three sources of strong phases are known to give clear signatures in the Dalitz plane [13,14].The short-distance direct CP violation is proportional to the difference of the magnitude between the positive and negative amplitudes of the resonance, and is therefore proportional to the square of the Breit-Wigner propagator associated with the resonance.
The interference term has two components.One is associated with the real part of the Breit-Wigner propagator and is directly proportional to (m 2 R − s), where m R is the central value of the resonance mass and s is the square of the invariant mass of its decay products.The other component is proportional to the product m R Γ, where Γ is the width of the resonance.The relative proportion of real and imaginary terms of these interference  Table 3: Signal yields and charge asymmetries in the rescattering regions of m(π + π − ) or m(K + K − ) between 1.0 and 1.5 GeV/c 2 .For the charge asymmetries, the first uncertainty is statistical, the second systematic, and the third is due to the CP asymmetry of the B ± → J/ψ K ± reference mode.
The charge asymmetry in B ± → π ± π + π − decays changes sign, as shown in Fig. 4, at a value of m(π + π − ) close to the ρ(770) resonance.This is an indication of the dominance of the long-distance interference effect in this region of the Dalitz plot.Moreover, since this change of sign occurs for both cos θ > 0 and cos θ < 0, the dominant term of the Dalitz interference is inferred to be the real part of the Breit-Wigner propagator.
When the interference between P-wave and S-wave involves two resonances, as for   B ± → K ± π + π − decays, with ρ(770) and f 0 (980) resonances, the dominant component of the real Dalitz CP asymmetry is proportional to a product of the type (m 2 ρ − s)(m 2 f 0 − s), thus having two zeros, as can be seen in Fig. 5.In the cos θ < 0 region there is a zero around the ρ(770) mass and another one around the f 0 (980) meson mass.However, in the region of cos θ > 0, a clear change of sign is only seen around the f 0 (980) mass.
The yield of the f 0 (980) resonance in B ± → K ± π + π − decays depends on cos θ (Fig. 5).The yield around the f 0 (980) mass for cos θ > 0 is almost twice that of the region with cos θ < 0. The magnitude of the CP asymmetry indicates the opposite dependence.Also, the yield around the ρ(770) resonance in the B ± → π ± π + π − decay changes significantly in the two cos θ regions (Fig. 4).The largest yield in the ρ(770) mass region with a small CP asymmetry occurs for cos θ < 0, while there is a large CP asymmetry with fewer events in the ρ(770) mass region for cos θ > 0. One possible explanation is that the fractions of tree and penguin contributions may vary across the phase space [15].Understanding this effect may be important to performing amplitude analyses.
The CP asymmetry around the ρ(770) peak in π + π − invariant mass changes sign depending on the invariant mass and angular distribution.To quantify the CP violation in the region where vector particles contribute without losing sensitivity to the interference asymmetry behaviour, the region is divided into four sectors.Sectors I and III are on the low-  Figure 9 shows the invariant mass distributions for B ± → π ± π + π − decays divided into these four sectors.The charge asymmetries are measured using the same method as for the inclusive asymmetries.Similar measurements of the CP asymmetries are performed for events restricted to regions dominated by the ρ(770) and K * (892) resonances in B ± → K ± π + π − decays, and the φ(1020) resonance in B ± → K ± K + K − decays, with the results given in Table 4.Only the decays involving the ρ(770) resonance have a significant CP asymmetry.The K * (892) charge asymmetry is consistent with zero, as expected and for the φ(1020) resonance, the results are in agreement with the previous LHCb analysis [42].
Table 4: Signal yields and charge asymmetries in the regions dominated by the vector resonances.For the charge asymmetries, the first uncertainty is statistical, the second systematic, and the third is due to the CP asymmetry of the B ± → J/ψ K ± reference mode.

Decay mode
Resonance
The CP asymmetry related to the ρ(770) resonance in the π + π − invariant mass below 1 GeV/c 2 is also reported.The behaviour of this asymmetry, which crosses zero around the ρ(770) mass in both B ± → K ± π + π − and B ± → π ± π + π − modes, indicates a CP asymmetry related to the real part of the long-distance interaction between the S-wave and P-wave contributions to π + π − .Further understanding of the resonance contributions and CP asymmetries in these decays will require amplitude analyses.

Figure 4 :
Figure 4: Projections in bins of the m(π + π − ) low variable of (a, b) the number of B − and B + signal events and (c, d) their difference for B ± → π ± π + π − decays.The plots are restricted to events with (a, c) cos θ < 0 and (b, d) cos θ > 0, with cos θ defined in the text.The yields are acceptance-corrected and background-subtracted.A guide line for zero (horizontal red line) was included on plots (c, d).

Figure 5 :
Figure 5: Projections in bins of the m(π + π − ) variable of (a, b) the number of B − and B + signal events and (c, d) their difference for B ± → K ± π + π − decays.The plots are restricted to events with (a, c) cos θ < 0 and (b, d) cos θ > 0. The yields are acceptance-corrected and background-subtracted.A guide line for zero (horizontal red line) was included on plots (c, d).

Figure 6 :
Figure 6: Projections in bins of the m(K + K − ) low variable of (a, b) the number of B − and B + signal events and(c, d) their difference for B ± → K ± K + K − decays.The inset plots show the φ resonance region of m(K + K − ) low between 1.00 and 1.05 GeV/c 2 , which is excluded from the main plots.The plots are restricted to events with (a, c) cos θ < 0 and (b, d) cos θ > 0. The yields are acceptance-corrected and background-subtracted.A guide line for zero (horizontal red line) was included on plots (c, d).

Figure 7 :
Figure 7: Projections in bins of the m(K + K − ) variable of (a) the number of B − and B + signal events and (b) their difference for B ± → π ± K + K − decays.The yields are acceptance-corrected and background-subtracted.A guide line for zero (horizontal red line) was included on plot (b).

Figure 9 :
Figure 9: Invariant mass distribution of B ± → π ± π + π − candidates restricted to (a) sector I, (b) sector II, (c) sector III and (d) sector IV.The left panel in each figure shows the B − candidates and the right panel shows the B + candidates.

Table 2 :
[34]ematic uncertainties on the measured asymmetries, where the total is the sum in quadrature of the individual contributions.The A π D uncertainty is taken from Ref.[34].
Università di Cagliari, Cagliari, Italy f Università di Ferrara, Ferrara, Italy g Università di Firenze, Firenze, Italy h Università di Urbino, Urbino, Italy i Università di Modena e Reggio Emilia, Modena, Italy j Università di Genova, Genova, Italy k Università di Milano Bicocca, Milano, Italy l Università di Roma Tor Vergata, Roma, Italy m Università di Roma La Sapienza, Roma, Italy n Università della Basilicata, Potenza, Italy o AGH -University of Science and Technology, Faculty of Computer Science, Electronics and Telecommunications, Kraków, Poland p LIFAELS, La Salle, Universitat Ramon Llull, Barcelona, Spain q Hanoi University of Science, Hanoi, Viet Nam r Università di Padova, Padova, Italy s Università di Pisa, Pisa, Italy t Scuola Normale Superiore, Pisa, Italy u Università degli Studi di Milano, Milano, Italy v Politecnico di Milano, Milano, Italy e