Measurement of CP observables in B± → DK± and B± → Dπ± with D → KS0K±π∓\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {K}_{\mathrm{S}}^0{K}^{\pm }{\pi}^{\mp } $$\end{document} decays

Measurements of CP observables in B± → DK± and B± → Dπ± decays are presented, where D represents a superposition of D0 and D¯0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\overline{D}}^0 $$\end{document} states. The D meson is reconstructed in the three-body final states KS0K±π∓\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {K}_{\mathrm{S}}^0{K}^{\pm }{\pi}^{\mp } $$\end{document} and KS0K∓π±\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {K}_{\mathrm{S}}^0{K}^{\mp }{\pi}^{\pm } $$\end{document}. The analysis uses samples of B mesons produced in proton-proton collisions, corresponding to an integrated luminosity of 1.0, 2.0, and 6.0 fb−1 collected with the LHCb detector at centre-of-mass energies of s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sqrt{s} $$\end{document} = 7, 8, and 13 TeV, respectively. These measurements are the most precise to date, and provide important input for the determination of the CKM angle γ.


Introduction
In the Standard Model (SM), CP violation in the hadronic sector is described by the irreducible complex phase of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix [1,2]. This matrix is unitary, which leads to the condition V ud V * ub + V cd V * cb + V td V * tb = 0, where V ij is the CKM matrix element relating quark i to quark j. This condition can be represented by a triangle in the complex plane with internal angles α, β, and γ. The angle γ is defined as γ ≡ arg (−V ud V * ub /V cd V * cb ), which is equal to arg (−V us V * ub /V cs V * cb ) up to O(λ 4 ) ∼ 10 −3 in the Wolfenstein parameterisation [3], where λ is the sine of the Cabibbo angle [1]. Improving knowledge of γ can be achieved in a theoretically clean manner by studying the interference of b → u and b → c transition amplitudes in tree-level b-hadron decays. Such a measurement provides a benchmark against which other flavour observables that are more susceptible to the influence of physics beyond the SM can be compared [4].
A combination of measurements from LHCb currently yields γ = (74.0 +5.0 −5.8 ) • [5,6], which is the most precise determination of γ from a single experiment. The precision is dominated by measurements exploiting the B + → DK + decay, 1 where D indicates a superposition of D 0 and D 0 mesons reconstructed in a final state common to both. To continue improving the precision on γ, independent measurements can be performed 1 The inclusion of charge-conjugate processes is implied throughout, unless otherwise indicated.
-1 -JHEP06(2020)058 using all suitable D meson final states. Several different final states have thus far been analysed at LHCb, including a previous measurement of the singly Cabibbo-suppressed D → K 0 S K − π + and D → K 0 S K + π − modes [7]. These decays are reconstructed in two categories by comparing the charge of the pion produced in the D decay with the charge of the B meson; B + → [K 0 S K + π − ] D h + decays are thus labelled "Same Sign" (SS), and B + → [K 0 S K − π + ] D h + decays are labelled "Opposite Sign" (OS), where h ∈ {π, K}. This paper reports an update to ref. [7], measuring CP observables in B + → DK + and B + → Dπ + decays using the D → K 0 S K + π − and D → K 0 S K − π + final states. Data corresponding to 6.0 fb −1 of integrated luminosity collected between 2015 and 2018 (Run 2) of data taking is used. The Run 1 dataset collected during 2011 and 2012 and corresponding to an integrated luminosity of 3.0 fb −1 is also reprocessed, to benefit from an improved selection as well as a reappraisal of the backgrounds.
In order to interpret interference effects involving multi-body D-decays, it is necessary to account for the amplitude structure of the Dalitz plot. Instead of employing an amplitude model to describe the contributing partial waves, the CLEO collaboration have made measurements of the effective amplitude and phase variation using a sample of quantumcorrelated D decays collected by the CLEO-c experiment [8]. Due to limited sample size, those measurements were performed averaging over large regions of the Dalitz plot, notably defining one of two regions to contain the D → K * (892) + K − mode. In the present work, results are reported for both the K * + and non-K * + regions of the Dalitz plot, respecting the boundary defined by CLEO-c. The use of external CLEO-c results, which were performed across the full Dalitz plot and within the K * + region, avoids the need to introduce a systematic uncertainty resulting from an amplitude model description.
The paper is organised as follows: section 2 presents the observables to be measured and their relationships to the physics parameters of interest; section 3 discusses the aspects of the detector, trigger, and simulation that are relevant for the measurement; sections 4, 5, and 6 describe the candidate selection, the fit to the invariant mass spectra, and the assignment of systematic uncertainties; the observable results are presented in section 7.

Formalism
The SS B + → [K 0 S K + π − ] D K + decay can proceed via the D 0 or D 0 states. As such, the total decay amplitude is given by the sum of two interfering amplitudes, where x represents the Dalitz plot coordinates (m 2 are the D 0 and D 0 decay amplitudes at a specific point in the K 0 S K + π − Dalitz plot [9]. The OS B + → [K 0 S K − π + ] D K + decay also proceeds via both D 0 and D 0 , with a total decay amplitude given by The amplitude ratio r B = |A(B + →D 0 K + )| |A(B + →D 0 K + )| ∼ 0.1 [5,6], and δ B = arg A(B + →D 0 K + ) A(B + →D 0 K + ) is the strong-phase difference between the B decay amplitudes. To calculate the decay rate in a -2 -JHEP06(2020)058 finite region of the Dalitz plot, the integral of the interference term over that region must be known. In ref. [8], measurements of quantum-correlated D decays have been used to determine the amplitude ratio, r D = |A(D 0 →K 0 s K + π − )| |A(D 0 →K 0 s K + π − )| , and the integral of the interference term directly in the form of a coherence factor, κ D , and an average strong phase difference, δ D [10]. The coherence factor is defined as A similar notation also holds for SS and OS B + → Dπ + decays with the replacements r B → r π B and δ B → δ π B , where r π B ∼ 0.015. In each Dalitz region, four decay rates are considered in this analysis [11]: Observables constructed from eq. 2.4 have sensitivity to γ that depends upon the value of the coherence factor, with a higher coherence corresponding to greater sensitivity. The CLEO-c results [8] show high coherence within the K * + region, defined as ±100 MeV/c 2 around the K * + mass; κ D = 0.94 ± 0.12 and δ D = (−16.6 ± 18.4) • are reported. With r D ≈ 0.6 [8], the maximal CP asymmetry that can be expected is 35% in B + → DK + decays, but only 2% in B + → Dπ + decays due to the dissimilarity of r D and r π B . Dedicated measurements in the non-K * + region have not yet been made. Eight yields are measured in this analysis, from which seven ratios are constructed as CP observables; each observable can be related to γ through the decay rates in eq. 2.4. The charge asymmetry is measured in four decay modes, where m ∈ {SS,OS} and h ∈ {π, K}. The ratios of B + → DK + and B + → Dπ + yields, R DK/Dπ m , are determined, and the ratio of SS to OS B + → Dπ + yields, R SS/OS , is also measured. The measurements are reported for the K * + region of the D Dalitz plot as defined above, and outside it; they are not interpreted in terms of γ in this work, as constraints on the B decay hadronic parameters which come from measurements using other D decay modes are necessary at the current level of statistical precision.

Detector and simulation
The LHCb detector [12,13] 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 -3 -

JHEP06(2020)058
vertex detector surrounding the pp interaction region, a 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 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 (RICH) 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.
The online event selection is performed by a trigger, 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. The events considered in the analysis are triggered at the hardware level when either one of the final-state tracks of the signal decay deposits enough energy in the calorimeter system, or when one of the other particles in the event, not reconstructed as part of the signal candidate, fulfils any trigger requirement. At the software stage, it is required that at least one particle should have high p T and high χ 2 IP , where χ 2 IP is defined as the difference in the PV fit χ 2 with and without the inclusion of that particle. A multivariate algorithm [14] is used to identify secondary vertices consistent with being a two-, three-, or four-track b-hadron decay. The PVs are fitted with and without the B candidate tracks, and the PV that gives the smallest χ 2 IP is associated with the B candidate.
Simulated events are used to describe the signal mass shapes and compute efficiencies. In the simulation, pp collisions are generated using Pythia [15] with a specific LHCb configuration [16]. Decays of hadronic particles are described by EvtGen [17], in which final-state radiation is generated using Photos [18]. The interaction of the generated particles with the detector, and its response, are implemented using the Geant4 toolkit [19,20] as described in ref. [21].

Offline selection
Decays of K 0 S mesons to the π + π − final state are reconstructed in two categories, the first containing K 0 S mesons that decay early enough for the pions to be reconstructed in the vertex detector, and the second containing K 0 S mesons that decay later such that track segments of the pions cannot be formed in the vertex detector. These categories are referred to as long and downstream, respectively. The candidates in the long category have better mass, momentum, and vertex resolution than those in the downstream category, but the downstream category contains more candidates and thus both are used. Herein, B + candidates are denoted long or downstream depending on which category of K 0 S candidate is used.

JHEP06(2020)058
The D (K 0 S ) candidates are required to be within ±25 MeV/c 2 (±15 MeV/c 2 ) of the known mass [22], and B + meson candidates with invariant masses in the interval 5080-5700 MeV/c 2 are retained. The kaons and pions originating from both the B + and D decays are required to have p T in the range 0.5-10 GeV/c and p in the range 5-100 GeV/c. These requirements ensure that the tracks are within the kinematic coverage of the RICH detectors, which are used to provide particle identification (PID) information.
A boosted decision tree (BDT) classifier [23] implementing the gradient boost algorithm is employed to achieve further combinatorial background suppression. The BDT is trained using simulated B + → Dh + decays as a proxy for signal and a background sample of candidates in data with invariant masses in the range 5900-7200 MeV/c 2 which are not used in the invariant-mass fit (see section 5). The input to the BDT is a set of features that characterise the signal decay. These features can be divided into two categories: (1) properties of any particle, and (2) properties of composite particles only (the D and B + candidates). Specifically: 1. p, p T , and χ 2 IP ; 2. decay time, flight distance between production and decay vertex, decay vertex quality, radial distance between the decay vertex and the PV, and the angle between the particle's momentum vector and the line connecting the production and decay vertices.
In addition, a feature that estimates the imbalance of p T around the B + candidate momentum vector is also used in the BDT. It is defined as where the sum is taken over tracks inconsistent with originating from the PV that lie within a cone around the B + candidate, excluding tracks used to make the signal candidate. The cone is defined by a circle with a radius of 1.5 units in the plane of pseudorapidity and azimuthal angle expressed in radians. Including the I p T feature in the BDT training gives preference to B + candidates that are isolated from the rest of the event.
Since no PID information is used in the BDT classifier, the efficiency for B + → DK + and B + → Dπ + decays is similar, with insignificant variations arising from small differences in the decay kinematics. The selection requirement applied to the BDT response is optimised by minimising the relative statistical uncertainty on the R SS/OS observable, as measured using the fit described in section 5. PID information from the RICH detectors is used to improve the purity of the B + → DK + samples. A strict PID requirement is applied to the companion kaon in B + → DK + to suppress contamination from B + → Dπ + decays where the companion pion is misidentified as a kaon. The requirement is around 70% efficient, and genuine B + → DK + decays failing the requirement are placed into the B + → Dπ + sample. Less than 0.5% of genuine B + → Dπ + decays pass the kaon PID requirement, and are placed into the B + → DK + sample. This results in a B + → Dπ + background of around 5% relative to the correctly identified B + → DK + signal. Background from the B + → [K 0 S π + π − ] D h + decay, which has a branching fraction around ten -5 -

JHEP06(2020)058
times larger than the signal, is suppressed by placing PID requirements on both the kaon and pion produced in the D decay. For long K 0 S candidates, the square of the flight distance significance with respect to the PV is required to be greater than 100 to suppress background from which peaks at the same invariant mass as the signal, is suppressed by requiring that the flight distance of the D candidate divided by its uncertainty is greater than 2. Where multiple candidates are found in the same event, one candidate is chosen at random, leading to a reduction in the sample size of approximately 2%.
For several quantities used in the selection and analysis of the data, a kinematic fit [24] is imposed on the full B + decay chain. Depending on the quantity being calculated, the D and K 0 S candidates may be constrained to have their known masses [22]. The fit also constrains the B + candidate momentum vector to point towards the associated PV, defined as the PV for which the candidate has the smallest χ 2 IP . These constraints improve the resolution of the calculated quantities, and thus help enhance the separation between signal and background decays. Furthermore, they improve the mass-squared resolution, which is important for identifying the Dalitz region assignment.
The Dalitz plots for selected candidates in the signal region ±25 MeV/c 2 around the B + mass are shown in figure 1; within this region, background from decays involving no charm meson constitute less than 5% of the total sample. The Dalitz coordinates are calculated from the kinematic fit with all mass constraints applied. A band corresponding to the intermediate state, D → K * (892) − K + is visible in each plot.

Invariant-mass fit
In order to measure the CP observables introduced in section 2, an extended binned maximum-likelihood fit to the invariant-mass distributions of the B meson candidates in the range between 5080 MeV/c 2 and 5700 MeV/c 2 is performed. The fit is performed simultaneously to all decay categories, in order to enable sharing of common parameters. A total of 16 categories are included in the fit: (DK, Dπ) × (SS, OS) × (long, downstream) × (B + , B − ). The fit range is between 5080 MeV/c 2 and 5700 MeV/c 2 in the B ± candidate invariant mass. The fit is performed separately for candidates within the K * + region and those outside.
To model the invariant-mass distribution, a total fit probability density function (PDF) is created from several signal and background components. Most of these are modelled using simulated signal and background samples reconstructed as the signal decay and passing all selection requirements. The components are: 1. Signal B + → DK + and B + → Dπ + decays, described by the sum of two Crystal Ball functions [25] with a freely varying common mean and width, and tail parameters fixed from simulation. A single freely varying parameter relates all B + → DK + widths to their B + → Dπ + counterparts. SS and OS decays share all shape parameters, but long and downstream decays have separate freely varying widths due to the  3. Partially reconstructed background from the B + → (D * 0 → D{π 0 /γ})h + , B + → Dh + {π 0 }, and B 0 → (D * − → D{π − })h + decays, where the particle in braces is not reconstructed. These components sit at lower invariant-mass values than the signal, and are described by PDFs constructed from a parabolic function to describe the decay kinematics. This function is convolved with the sum of two Gaussian functions with a common mean in order to describe the detector resolution, as further -7 -JHEP06(2020)058 described in ref. [26]. All shape parameters are fixed from simulation. All partially reconstructed background yields vary freely, but the B 0 → D * − h + component yields are required to be equal in the B − and B + samples; the fast B 0 oscillation renders CP violation effects negligible in this time-integrated measurement.
4. Partially reconstructed background from B 0 s → DK + {π − } decays, where the pion produced in the B 0 s decay is not reconstructed, contributes in the B + → DK + samples. These decays are modelled using a PDF with fixed shape parameters based on the m(DK) distribution observed in ref. [27]. The yield of this component freely varies in the SS and OS samples, but the yields are required to be equal in the B − and B + samples as the fast B 0 s oscillation renders CP violation effects negligible.
5. Charmless background contributions remain in the B + → DK + samples after application of the D-meson flight requirement described in section 4. They are estimated using fits to the B + candidate invariant-mass distributions in data, where candidates falling in the lower sidebands of the D candidate invariant-mass spectra are considered. The charmless contributions are included as fixed-shape Gaussian functions from simulation, with fixed yields as determined by the sideband fits. 6. Backgrounds from particle misidentification, which arise due to the imperfect efficiency of RICH PID requirements applied to companion hadrons in order to separate B + → DK + and B + → Dπ + decays. The efficiencies of the PID requirements are determined using calibration samples of high-purity decays that can be identified without the use of RICH information [28]. Given a PID efficiency K PID ∼ 0.7 for B + → DK + decays, a fixed fraction (1 − K PID ) of the total B + → DK + signal yield is assigned to a PDF in the corresponding B + → Dπ + sample. This component is described by a Crystal Ball function with all shape parameters fixed to those found in simulation; due to the companion hadron misidentification, this component falls below the nominal B + mass. In the same fashion, a small component is included in the B + → DK + sample to model misidentified B + → Dπ + decays, with a yield that freely varies to ∼ 0.4% of the total B + → Dπ + yield. This component is also described by a Crystal Ball function, with all shape parameters fixed to the values found in simulation.
In order to measure CP asymmetries, the detection asymmetries for K ± and π ± mesons must be taken into account. In the fit, a detection asymmetry of (−0.51±0.28)% is assigned for each kaon in the final state, primarily due to the fact that the nuclear interaction length of K − mesons is shorter than that of K + mesons. The value used is computed by comparing the charge asymmetries in D + → K − π + π + and D + → K 0 S π + calibration samples, weighted to match the kinematics of the signal kaons [29]. The equivalent asymmetry for pions is smaller, (−0.06 ± 0.04)% [30]. All measured CP asymmetries are also corrected in the fit for the asymmetry in B ± production, which has a value (+0.14 ± 0.07)% based on measurements of this quantity made in refs. [26] and [30].
To measure the R decays. The efficiency ratio is found to be (DK) (Dπ) = 1.012 ± 0.016, which is employed as a fixed correction term in the fit. A similar correction is applied to the R SS/OS observable, to account for differences in selection efficiency for SS and OS decays caused by efficiency variation across the Dalitz plot. The correction is determined using simulated B + → Dπ + decays and the D → K 0 S K + π − and D → K 0 S K − π + amplitudes measured by LHCb in ref. [31]. The correction is determined in bins across the Dalitz plot, and an average value is calculated to be η = 1.090 ± 0.008 (1.007 ± 0.013) within (outside) the K * + region. s categories are shown together. In table 1, the measured signal yields for each D final state are provided for both the K * + and the non-K * + regions. The fit strategy is validated using pseudoexperiments, and is found to be unbiased for all parameters.

Systematic uncertainties
All of the CP observables measured in this work are constructed as ratios of topologically identical final states. As such, the majority of potential systematic uncertainties cancel with the residual systematic uncertainties detailed here. Small differences in efficiency between B + → DK + and B + → Dπ + decays are corrected using simulation as described in section 5, where the uncertainty on the correction arises due to the finite size of the simulated samples. The correction is varied within its uncertainty to determine the systematic uncertainty. The variation in efficiency across the Dalitz plot causes a difference in the total efficiency of SS and OS decays. An appropriate correction is applied to the R SS/OS observable, with an uncertainty arising from the use of a binned procedure to calculate the average correction.
Several fixed shape parameters are used in the fit, including the signal tail parameters and background PDFs. All fixed shape parameters are determined from fits to simulated samples, and are varied to calculate the propagated systematic uncertainty. Charmless backgrounds are modelled as fixed yield components in the invariant-mass fit. The yields are varied within their respective uncertainties to determine the systematic uncertainty. Each charmless component has a fixed CP asymmetry of zero in the fit; their asymmetries are independently varied according to a Gaussian of width 0.1 to determine the systematic uncertainty. This width chosen to align with the degree of CP asymmetry observed in the charmless background present in measurements of B + → [h + h − ] D h + decays. [26,32].  All measured CP asymmetries are corrected for the B ± production asymmetry as well as for the kaon and pion detection asymmetries where relevant. These corrections are applied as fixed terms in the invariant-mass fit, and are varied within their associated uncertainties to determine the systematic uncertainty. A fixed PID efficiency is used to determine the fraction of B + → DK + signal decays that are misidentified as B + → Dπ + . This efficiency is known within 1% relative uncertainty, and is varied within this range to determine the systematic uncertainty.
The systematic uncertainties for each CP observable, quoted as a percentage of the statistical uncertainty, are listed in tables 2 and 3. The category Eff relates to efficiency corrections, PDF to fixed shape parameters, Cls to charmless background yields and asymmetries, Asym to asymmetry corrections, and PID to the PID efficiency. The total systematic uncertainties are given by the sum in quadrature of each contributing systematic.   Table 2. Systematic uncertainties for the K * + region fit. Uncertainties are quoted as a percentage of the statistical uncertainty for a given observable, and the total uncertainty is given by the sum in quadrature of each contribution.    Table 3. Systematic uncertainties for the non-K * + region fit. Uncertainties are quoted as a percentage of the statistical uncertainty for a given observable, and the total uncertainty is given by the sum in quadrature of each contribution.

Results
The results for the K * + region of the Dalitz plot are The results are in agreement with ref. [7], and all statistical uncertainties are reduced in accordance with the increased signal yields. The systematic uncertainties on each asymmetry are reduced considerably due to improved knowledge of the B ± production asymmetry and the kaon detection asymmetry. The systematic uncertainties on R DK/Dπ SS , R DK/Dπ OS , and R SS/OS are also reduced, due to the use of larger simulated samples. All observables are statistically limited with the current data set. The statistical and systematic correlation matrices for the CP observables are given in appendix A.
A comparison of the K * + region results with the SM expectation is made by calculating the CP observables from the current best-fit values of γ = (74.0 +5.0 −5.8 ) • , δ B = (131.2 +5.1 −5.9 ) • , and r B = (9.89 +0.51 −0.50 )% for B + → DK + decays [5]; no comparison is made using the non-K * + results, since the required charm hadronic parameters have not yet been measured.
For these inputs, the 68% and 95% confidence-level expectation intervals are displayed in figure 6, together with the results presented herein. The dominant uncertainty contribution to the expectation intervals comes from the D-decay parameter inputs. The measurements are found to be compatible with the SM expectation, where the χ 2 per degree of freedom is found to be 1.56 taking into account the uncertainties and correlations of both the measurements and the expected values; the corresponding p-value for rejection of the SM hypothesis is 0.14.

Conclusion
Measurements of CP observables in B + → DK + and B + → Dπ + decays with the D meson decaying to K 0 S K + π − and K 0 S K − π + are performed using LHCb data collected in Run 1 and Run 2. The results are in agreement with the SM, and supersede those of the previous study [7], benefiting from the increased data sample and improved analysis methods. The measurements presented in this paper improve the precision of several of the CP observables used in global fits for γ, which will contribute to improved precision on γ and on the hadronic parameters r B and δ B for these decays. Improved measurements of charm hadronic parameters in both the K * + and non-K * + regions would also benefit the interpretation of these results and the constraints on γ that can be obtained from them.

A Correlation matrices
Statistical and systematic correlation matrices for the seven CP observables are given in tables 4−7, for both the K * (892) ± region and non-K * (892) ± region results.   Table 7. Systematic correlation matrix for the non−K * (892) ± region fit.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.