An empirical model to determine the hadronic resonance contributions to B ¯ 0 → K ¯ ∗ 0 μ + μ - transitions.

A method for analysing the hadronic resonance contributions in B ¯ 0 → K ¯ ∗ 0 μ + μ - decays is presented. This method uses an empirical model that relies on measurements of the branching fractions and polarisation amplitudes of final states involving J PC = 1 - - resonances, relative to the short-distance component, across the full dimuon mass spectrum of B ¯ 0 → K ¯ ∗ 0 μ + μ - transitions. The model is in good agreement with existing calculations of hadronic non-local effects. The effect of this contribution to the angular observables is presented and it is demonstrated how the narrow resonances in the q 2 spectrum provide a dramatic enhancement to C P -violating effects in the short-distance amplitude. Finally, a study of the hadronic resonance effects on lepton universality ratios, R K ( ∗ ) , in the presence of new physics is presented.

Introduction

Decays with a transition receive contributions predominantly from loop-level, flavour changing neutral current transitions. These transitions are mediated by heavy (short-distance) particles and are suppressed in the Standard Model (SM). Over the last few years, discrepancies have emerged when comparing measurements of the properties of decays to SM predictions [1-10]. Global analyses of these decays imply that there might be a new vector current which is destructively interfering with the SM contribution to the decay, producing inconsistency with the SM at the 4–5 [11-17].

In this paper, the possibility that hadronic resonances are interfering with the short-distance amplitude and mimicking physics beyond the SM is considered. This is because in addition to the short-distance contribution to decays, the same final state can be obtained through non-local transitions, where denotes a quark-anti-quark pair. An example of such a decay is the decay , where the meson decays into two leptons.1 As the decay rate of this process is two orders of magnitude larger than its short-distance counterpart, sizeable interference effects are possible far from the mass.

The approach presented in this paper models the hadronic contributions originating from charm and light quark resonances as Breit–Wigner amplitudes. This approach is inspired by Refs. [18, 19] and is used to describe the hadronic resonances across the full dimuon mass spectrum of decays. The LHCb collaboration performed a measurement of the interference between the non-local and short-distance components of decays by modelling the hadronic resonance contributions as Breit–Wigner amplitudes [20]. The level of interference was found to be small and the measurement of the short-distance component was found to be compatible with that of previous interpretations.

These non-local contributions are difficult to calculate and to date there is no consensus as to whether the deviations seen in global analyses can be explained by the these intermediate hadronic contributions, or by physics beyond the SM. Differentiating between these two hypotheses is of prime importance for confirming the existence and subsequently characterising phenomena not predicted by the SM. More detailed discussions on this point can be found in Refs. [18, 19, 21–28].

Due to the more complex amplitude structure of the decay, for each resonant final state there are three relative phases and magnitudes that need to be determined instead of one in the case of the decay. Existing measurements of the branching fractions of and decays, together with measurements of their polarisation amplitudes [29-32] can be used to assess the impact of these decays to the observables of the process, up to a single overall phase per resonance that needs to be determined through a simultaneous fit to both the short-distance and non-local components in the final state. In the absence of such a measurement, scanning over all possible values for the global phase for each resonant final state, results in a prediction of the range of hadronic effects that can be compared to more formal calculations. The angular distribution of the decay is sensitive to the strong-phases of non-local contributions, particularly through the observables and . This sensitivity allows for a data-driven extraction of the non-local parameters of the proposed model.

The level of violation in decays such as depends on weak- and strong-phase differences with interfering processes, such as . Therefore, a model for the strong phases of the non-local contributions to transitions, offers new insight on both the kinematic regions where violation might be enhanced, as well as what the level of enhancement could be.

An increasingly large part of the discrepancy in transitions is being driven by tests of lepton universality in decays [3, 33, 34]. These deviations cannot be explained by hadronic effects (the meson, for example, decays equally often to electrons and muons). Although a significant deviation from lepton-universality would be a clear indication of physics beyond the SM, the precise characterisation of the new physics model still depends on the treatment of hadronic contributions. The angular distribution of decays is critical in order to both determine the size of the new physics contribution, as well as to distinguish between models with left- or right-handed currents giving rise to new vector and axial-vector couplings.

This paper is organised as follows: Section 2 describes the model of the non-local contributions as well as the experimental inputs; Section 3 presents the comparison of the model to existing calculations; Section 4 shows how current model uncertainties impact both -averaged and -violating observables of decays, as well as the expected precision of the observables using the data that is expected from the LHCb experiment by the end of Run 2 of the LHC; finally in Section 5 there is a discussion of the impact of the non-local contributions in and transitions in the presence of lepton-universality violating physics.

The model

The differential decay rate of transitions, where the is a P-wave state and ignoring scalar or timelike contributions to the dimuon system, depends on eight independent observables [35]. Each of these observables is made up of bilinear combinations of six complex amplitudes representing the three polarisation states of the for both the left- and right-handed chirality of the dilepton system. The expression for the differential decay rate in terms of the angular observables and their subsequent definition in terms of amplitudes, can be found in Ref. [36]. The decay amplitudes are written in terms of the complex valued Wilson Coefficients , and , encoding short distance effects, and the dependent form-factors , given in Ref. [15], that express the matrix elements of the operators involved in these decays. The coefficient corresponds to the coupling strength of the vector current operator, to the axial-vector current operator and to the electromagnetic dipole operator. A detailed review of these decays, including the operator definitions and the numerical values of the Wilson Coefficients in the SM, can be found in Ref. [36]. The decay amplitudes in the transversity basis and assuming a narrow can be written as where , and are the masses of the B-meson, -meson, and lepton respectively, denotes the mass of the dimuon system squared, , andIn the above expressions, and for the remainder of this analysis, contributions from right handed Wilson Coefficients have been omitted. These are numerically small or zero in the SM and are not currently favoured by global analyses of processes. Following Ref. [15], the form factors are writtenwhere the z function is given bywith and . The parameters are taken from Ref. [15] and the coefficients including their correlations are taken from a combined fit to light-cone sum rule calculations and Lattice QCD results given in Refs. [15, 37].

The functions describe the non-local hadronic contributions to the amplitudes and are given by a simplistic empirical parametrisation inspired by the procedure of Refs. [18-20]. In particular where the sum in the above expressions represents the coherent sum of vector meson resonant amplitudes with and the magnitude and phase of each resonant amplitude relative to . The exact normalisation of the parameters is shown in Appendix A and is chosen such the integral of the sum of the squared magnitudes of the amplitude of a given amplitude produce the correct experimental branching fraction. Similarly, the parameters and need to be determined from experimental measurements. In this analysis, the central values of these parameters are set to zero, unless otherwise specified.

The dependence of each resonant amplitude is given by . As indicated from the analysis of the dimuon mass spectrum in decays from Ref. [20], the resonances considered in this analysis are the , , , , , and . Contributions from light-quark resonances are expected to be either CKM- or loop-suppressed compared to final states occurring through charmonium resonances. As experiments accumulate more data, additional broad light-quark states will start becoming statistically significant and can be easily incorporated in the model. For simplicity, is modelled by a relativistic Breit–Wigner function given bywhere and are the pole mass and natural width of the resonance and their values are taken from Ref. [38]. The running width is given bywhere p is the momentum of the muons in the rest frame of the dimuon system evaluated at q, and is the momentum evaluated at the mass of the resonance.

This isobar approach, although not rigorous, it provides a model for the strong phase variation of the amplitude across the full spectrum. This variation can result in sizeable effects even far from the pole of the resonances as discussed in Sect. 4.

It is customary that for each helicity amplitude, the expressions of the non-local components are recast as shifts to the Wilson coefficient , referred to as . This convention is particularly useful for comparisons with formal predictions of the non-local contributions.

Measurements of decays, where V denotes any state, are only sensitive to relative phases of the three transversity amplitudes. Therefore, the convention used in previous measurements of these modes is such that phases and are defined relative to . Using this convention, the remaining phase difference of each resonant polarisation amplitude relative to the corresponding short-distance one, is given by .

Experimental input

In order to assess the impact of the resonances appearing in the dimuon spectrum of decays, knowledge of the resonance parameters and appearing in Eqs. 7–9 is required. The amplitude analyses of and transitions performed by the LHCb, BaBar and Belle collaborations [30, 31, 39] constrain the relative phases and magnitudes of the transversity amplitudes of the resonant decay modes. Combined with the measured branching fractions of these decays by the Belle experiment [30, 32], the parameters and are determined up to an overall phase, , relative to the short-distance amplitude for the decay. Similarly, the amplitude components of transitions have been determined up to an overall phase, through the amplitude analyses and branching fraction measurements given in Refs. [40-42].

For the decay , the magnitude of the total decay amplitude is set using the world average branching fraction of this transition [38, 43, 44]. As no amplitude analysis of this mode has been performed, the relative phases and magnitudes of the transversity amplitudes are taken to be the same as those of the decay. As the overall contribution of the is expected to be small, this assumption will not impact the main conclusions of this study.

No measurements exist for final states involving the , and resonances, denoted as . To estimate the contributions of these final states, the relative phases and magnitudes of the transversity amplitudes are taken from the amplitude analysis of decays. An approximate value of the branching fraction of each of the modes is obtained by scaling the measured branching fraction of the decay , with , by the known ratio of branching fractions between and decays, with , given in Ref. [20]. The values used for the relative amplitudes and phases for each resonant contribution are summarised in Table 1.

Table 1

Summary of the input values used to model the non-local amplitude components . The input values rely on measurements given in Refs. [20, 29–32, 40–44]. The phases are measured relative to . As the measurements are given for the decay of the meson, in order to convert to the decay of the , the phase given in the table above must be shifted by

Mode	(ηj‖,θj‖[rad])	(ηj⊥,θj⊥[rad])	ηj0
B^0→ρ^0K^∗0	(1.5, 2.6)	(1.9, 2.6)	5.1×10^-1
B^0→ϕK^∗0	(2.5×10^+1,2.6)	(3.2×10^+1,2.6)	1.0×10^+1
B^0→J/ψK^∗0	(4.9×10^+3,-2.9)	(6.5×10^+3,2.9)	7.1×10^+3
B^0→ψ(2S)K^∗0	(5.3×10^+2,-2.8)	(8.1×10^+2,2.8)	9.6×10^+2
B^0→ψ(3770)K^∗0	(9.3×10^-1,-2.9)	(1.5, 2.9)	1.7
B^0→ψ(4040)K^∗0	(2.9×10^-1,-2.9)	(5.6×10^-1,2.9)	6.0×10^-1
B^0→ψ(4160)K^∗0	(8.3×10^-1,-2.9)	(2.0, 2.9)	1.8

Model comparisons

The study presented in Ref. [28] provides a prediction of the non-local charm loop contribution to decays. It relies on QCD light-cone sum rule calculations of matrix elements for and extrapolated to larger through a hadronic dispersion relation. The extrapolation uses input from experimental measurements of the rate and amplitude structure of and decays. As this calculation does not account for the factorisable next-to-leading order corrections to the charm loop, all phases of the non-local relative to the short-distance amplitudes are set to zero.

Figure 1 shows the parametrisation of the non-local contributions in the invariant amplitude basis of decays given in Ref. [28]. The relation of this amplitude basis to the helicity basis is also given in Ref. [28]. The predictions using the model described in Sect. 2, where only the contributions from the and resonances are considered, are shown for comparison. The free phases and appearing in Eqs. 7–9 are both set to 0 or . As a consistency check, the model presented in this paper is also shown, with the phases of all transversity amplitudes set to zero. The parameters and also appearing in Eqs. 7–9 are chosen such that they are broadly consistent with the values of Ref. [15] and the predictions of Ref. [28], with and . Ignoring all phases of the transversity amplitudes of and decays, the model of described in this analysis is consistent to that of Ref. [28]. However, accounting for the measured relative phases in the resonant decay amplitudes results in large differences between the two models. The level of disagreement depends on the value of the free phases and . The effect of the non-local charm contributions in Ref. [28] are known to move the central value of predictions of angular observables such as further away from experimental measurements [16]. However, this effect is only true due to the fact that the analysis of Ref. [28] did not account for the phases of the resonant amplitudes. An assessment of the impact of the phases on the angular observables is discussed in Sect. 4.

Fig. 1

The non-local contributions to the invariant amplitudes as a function of . The prediction using the model discussed in Sect. 2 is shown, where only the contributions from the and resonances are considered. The free phases and are both set to 0 (cyan solid line) or (cyan dashed-dotted line). The prediction where all phases of the and appearing in Eqs. 7–9 are set to zero is also depicted (black solid line), alongside the prediction from Ref. [28] (magenta band)

Building on the ideas of Ref. [28], a recent analysis presented in Ref. [21] provides a prediction of the non-local charm contribution that is valid up to a . This prediction also makes use of experimental measurements of and decays. In contrast to Ref. [28], the calculations of the non-local contributions are performed at to next-to-leading order in . The parametrisation is given by a z-expansion truncated after the second order as in Eq. (5). Figure 2 shows both the real and imaginary parts of the non-local contributions to decays presented in Ref. [21]. As the correlations between the z-expansion parameters are not provided, only the central values of the predictions are shown. The phase convention used in Ref. [21] is such that the transversity amplitudes of the and decays are related to those presented in this study through . The model described in Sect. 2, where only the contributions from the and resonances are considered, is in qualitative agreement with that of Ref. [21] for the following parameter choice: , , and . The small level of disagreement observed in the imaginary part of the amplitudes at low is due to the choice of setting , with smaller values giving a better agreement.

Fig. 2

The non-local contributions to the transversity amplitudes as a function of . The real (left) and imaginary (right) components are shown separately. The prediction from Ref. [21] is shown (magenta points). Predictions using the model discussed in Sect. 2, where only the contributions from the and resonances are considered, are overlaid for different choices of the phases and (cyan lines). See text for further details

To conclude, the simplistic model of the non-local contributions to decays presented in this paper is in good agreement with existing models, provided appropriate choices of , , and . For the latter, a larger value is required to match the predictions of Ref. [21], compared to Ref. [28]. The expressions of have sufficient freedom to capture the dependence of formal theory predictions in the range . In addition, in contrast to current predictions, the model of can naturally accommodate hadronic contributions from states composed of light quarks such as the and , as well as resonances appearing in the region , where denotes the mass of the D-meson. This is due to the use of Breit–Wigner functions to approximate the resonant contributions, that experiments can easily adopt.

Effect on angular observables

Using the model of described in Sect. 2, the effect of the hadronic resonance contributions on the angular observables of decays can be estimated. Figure 3 shows the distribution of the angular observables , , and [45, 46] in the SM. The observable exhibits a particularly large dependence on the strong phases, demonstrating that measurements of the angular distribution of decays can be used to determine the phases of the hadronic resonances. Therefore, this observable can be used to separate short-distance from the non-local contributions, as only the non-local part has a strong-phase difference. The remaining -averaged observables can be found in Appendix B. Definitions of these observables can be found for instance in Ref. [47]. As the phase of all the resonant final states appearing in Table 1 are unknown, all possible variations of phases are considered. The uncertainties arising from the combined light-cone sum rules and lattice QCD calculations of form factors are accounted for using the covariance matrix provided in Ref. [15]. The predictions of these observables using flavio [48] are also shown for comparison. The lack of knowledge of the phase results in a large uncertainty for the prediction of , diluting the sensitivity of this observable to the effects of physics beyond the SM. However, for the choice of that results in a non-local charm contribution that is compatible with the latest prediction presented in Ref. [21] and is shown in Fig. 2), the tension of the prediction with the measured value of cannot be explained solely through hadronic effects.

Fig. 3

Distributions of the angular observables , , and as a function of for regions below (left) and above (right) the open charm threshold (cyan). Specific choices are highlighted for (hatched band) and (dark band). The measured values of the observables from Ref. [49] are also shown (black points). The theoretical predictions (magenta band) using flavio [48] are shown for comparison

Sensitivity to violation

The model of the hadronic resonance contributions to decays described in this paper provides a prediction for the strong phase differences involved in these transitions. Direct violation will arise when there are interfering amplitudes that have different weak phases as well as different strong phases, as discussed within the context of and decays in Refs. [24, 50]. Therefore, it is interesting to study the effect that potential weak phases beyond the SM have on angular observables such as the direct asymmetry , defined aswhere and correspond to the partial widths of the decays and respectively, as well as the so-called -odd angular observables , defined for instance in Ref. [36]. The method is similar to what is discusssed in Ref. [51] for semileptonic charm decays.

Figure 4 shows the observables , and for of all resonances set to , 0, and . To illustrate the effect that the model of the strong phase differences have in the presence of new weak phases, two new physics models are considered which are compatible with existing experimental constraints. One with , and one with both and [52, 53]. The notation denotes the new physics contribution to the corresponding Wilson Coefficient. In both these models, all other Wilson Coefficients are set to their SM values. It is clear that the non-local contribution enhances -violating effects in these decays, with the level of this enhancement depending on the value of the unknown phase . As it can be seen, there is a huge effect in the vicinity of the resonances, thus giving sensitivity to an imaginary component of in a way which have not been considered before. The only other viable way to gain sensitivity would be through a time dependent analysis of the or the with the decaying to the eigenstate . In contrast, -violating effects arising through a weak phases appearing in the Wilson coefficient , are best constrained from measurements of decays [53].

Fig. 4

Distribution of observables , and as a function of , for of all resonances set to , 0, and . Two new physics models are considered, one with (left), and one with , (right)

Expected experimental precision

The experimental sensitivity to the phases between the short-range and hadronic resonance contributions to is determined using simulated decays that include contributions from both short-distance and non-local components. The size of this sample corresponds to the approximate number of decays expected2 to be collected by the LHCb experiment by the end of Run2 of the LHC [49]. The decays are generated with the parameters , and set to zero. The S-wave contribution to decays is accounted for using the angular terms and amplitude expressions as a function of the invariant mass of the system given in Refs. [54, 55]. In addition to the S-wave component for the short-range amplitude, S-wave components are introduced with an amplitude and phase (, ), for the and the resonances, based on the measurements given in Refs. [29, 30]. The overall effect of the S-wave contribution to the remaining resonances is considered to be negligible and is therefore ignored. In this study, all Wilson Coefficients are assumed to be real.

In order to ascertain the statistical precision on the non-local contribution, the detector resolution in needs to be accounted for by smearing the spectrum of the simulated events. For simplicity, a Gaussian resolution function is used with a width based on the RMS value of the dimuon mass resolution provided in Ref. [20], and converted into a resolution in . As the resolution in the helicity angles are far better than the variations in the angular distributions, any resolution effect in angles can be ignored; the sharp shape of the , and resonances mean that a similar argument is not valid for the distribution.

A four dimensional maximum likelihood fit is performed to the , , and distributions of the decays in this sample. Both the non-local parameters, including and , as well as the Wilson Coefficients and are left to vary in the fit. The form factor parameters however are fixed to their central values given in Ref. [15]. The resulting covariance matrix is used to ascertain the statistical precision on . Based on the assessment of the systematic uncertainties in Ref. [20], the dominant source of experimental uncertainty is expected to be statistical in nature. However, the presence of tetra-quark states appearing in and decays [30, 56] will impact the determination of the non-local parameters. Although the effect is expected to be small, an accurate assessment of the effect is beyond the scope of this study.

The statistical precision on the angular observables is estimated by generating values for the non-local parameters of , according to a multivariate Gaussian distribution centred at the values used to simulate the decays, with a covariance matrix obtained from the resulting fit to the simulated data. These values are then propagated to the angular observables in order to obtain their 68% confidence interval as a function of . Figure 5 shows the statistical precision to , , and in the SM, where the non-local parameters are given by Table 1 with . The equivalent plots for the remaining -averaged observables can be found in Appendix C.

Fig. 5

Predictions of the observables , , and in the SM using the expected post-fit precision of the non-local parameters at the end of Run2 of the LHC. A sample of simulated decays that include contributions from both short-distance and non-local components, is used to determine the parameters of . The decays are simulated in the SM, with the parameters , and set to zero. The 68% confidence intervals are shown for the statistical uncertainty (cyan band) and the combination of the statistical uncertainty with the form-factor uncertainties (magenta band) given in Ref. [15]

By the end of Run2 of the LHC, the dominant theoretical uncertainty of the angular observables in the region  , will be due to the knowledge of the form-factors, rather than the non-local components. Future runs of the LHC will result in an even larger number of decays. Therefore, it will, in a fit that combines the experimental data and the form factor uncertainties [57], be possible to use experimental data to further constrain Wilson Coefficients, as well as improve the precision of form factors and non-local contributions from charm and light quark resonances.

Hadronic resonance effects in tests of lepton universality

Recent tests of lepton universality in decays have revealed hints of non-universal new physics entering in the dimuon Wilson Coefficient [11, 12, 58, 59]. The level of this potential new physics effect is compatible with the observed anomalies in the amplitude analyses and branching fraction measurements of transitions. Lepton universality tests rely on measurements such as the ratios of branching fractions between decays with muons and electrons in the final state. The observables and are defined asHadronic effects in decays are lepton universal and observables such as and can be predicted precisely in the SM, due to the cancellation of hadronic uncertainties. Therefore, any significant deviation between measurements and predictions of these quantities is a clear sign of physics beyond the SM. However, in the presence of new physics effects that enter through the Wilson Coefficient , the cancellation of hadronic uncertainties is no longer exact. Consequently, in order to determine the exact nature of any potential new physics model, an accurate determination of the non-local contributions in decays is essential.

The model of discussed in Sect. 2 is used to provide a prediction for that accounts for the residual dependence on the unknown phases . Figure 6 summarises this prediction in models with values of between -0.5 and -2.0, as suggested by global analyses of transitions. The confidence interval for is determined by considering the full variation of the unknown phases . The residual form factor uncertainty is found to be subdominant compared to the variation of the phase. A prediction for is also provided, which uses the long distance contributions measured in Ref. [20] with the 68% confidence interval determined by treating the measured non-local parameters as uncorrelated. It can be seen that when the experimental data is used for measuring the phase of the non-local contribution, the residual uncertainty becomes very small. It is worth noting that for , there is no dependence on the unknown phase . Tabulated values of these predictions can be found in Appendix D. In the presence of new physics entering the Wilson coefficient , a modest variation of with the unknown phase is observed. However, this variation is around 6 times smaller than the estimated uncertainty of in the presence of lepton non-universal effects suggested by Ref. [11].

Fig. 6

Predictions of at large recoil (hatched magenta) and low recoil (hatched cyan), and at large recoil (solid burgundy) for different values of . The values at low recoil are identical to those at large recoil and thus not shown. The interval for is determined using the model described in Sect. 2, considering the full variation of the unknown phases . In contrast the 68% confidence interval of the prediction is obtained using the measured non-local contributions in decays [20]

Conclusions

An empirical model to describe the hadronic resonance contributions in transitions that relies on measurements of the branching fractions and polarisation amplitudes of decays, is presented. For a particular choice of the relative phases between the short-distance component and the hadronic amplitudes, this model was found to be in good agreement with more formal predictions such as those of Refs. [21, 28]. The approach of this paper can naturally accommodate broad hadronic contributions from states such as the , the and charm-resonances above the open charm threshold, which can be inserted into experimental analyses of decays.

The lack of knowledge of the longitudinal phase differences between and decays results in a larger uncertainty on the predictions of the angular observables of decays compared to current approaches. A measurement of these phases is critical as it will reduce the uncertainty in the determination of the Wilson Coefficients.

In addition, the resonant contributions to the decay provide large strong-phase differences that enhance sensitivity to CP violating effects. In this way, there is no need to rely on a time dependent analysis to a eigenstate. For the method to be exploited, it is required to have a model of the strong phase differences between short- and non-local contributions to transitions as proposed here.

In the SM, observables such as and are independent of hadronic uncertainties. However, in the presence of non-universal effects in transitions, these observables receive uncertainties from both the form-factor calculations and the interference between short- and non-local amplitudes. Using the models described in Ref. [20] and in this paper, predictions for and are provided for various choices of the Wilson coefficient . In order to maximise the potential of observables such as as a way of characterising the exact physics model behind potential lepton-universality violating effects, a measurement of the non-local contributions in decays is crucial. The data sample that will be collected by the LHCb experiment by the end of Run2 of the LHC will allow for a simultaneous amplitude analysis of both short-distance and non-local contributions to decays across the full spectrum of the decay. The model described in this paper, allows for a precise determination of both of these components.

Table 2

Predictions of and at large and low recoil for different values of . The interval for is determined using the model described in Sect. 2, considering the full variation of the unknown phases . The uncertainty due to the residual form factor dependence is found to be subdominant. In contrast, the 68% confidence interval of the prediction is obtained using the measured long distance contributions in decays [20]

Observable	C9μNP=-0.5	C9μNP=-1.0	C9μNP=-1.5	C9μNP=-2.0
	1.1<q^2<6.0 GeV^2/c^4
R_K^∗	[0.902, 0.912]	[0.827, 0.850]	[0.769, 0.808]	[0.727, 0.784]
R_K	[0.888, 0.889]	[0.792, 0.794]	[0.712, 0.718]	[0.651, 0.658]
	15<q^2<19 GeV^2/c^4
R_K^∗	[0.889, 0.894]	[0.796, 0.806]	[0.719, 0.735]	[0.658, 0.680]
R_K	[0.888, 0.889]	[0.792, 0.794]	[0.712, 0.718]	[0.651, 0.658]