Zooming in on $B\to K^\ast \ell\ell$ decays at low recoil

We analyse $B\to K^\ast \ell\ell$ decays in the region of low hadronic recoil, where an operator product expansion (OPE) in $1/m_b$ applies. Using a local model for charm contributions based on $e^+ e^- \to hadrons$ against the OPE provides a data-driven method to access the limitations to the OPE's accuracy related to binnings in the dilepton mass. Model-independent fits to $B\to K^\ast \mu \mu $ low recoil angular observables exhibit presently only small sensitivity to different charm models. They give similar results as the fits based on the OPE, and are in agreement with the standard model, but leave also room for new physics. Measurements with resolution small enough to probe charm resonances would be desirable.


I. INTRODUCTION
Rare (semi-)leptonic decays induced by b → s flavor-changing neutral current (FCNC) transitions are highly suppressed in the Standard Model (SM) and therefore sensitive to effects from non-standard interactions. The corresponding exclusive B-meson decays have been investigated by the experimental collaborations LHCb [1,2], CMS [3], CDF [4], Belle [5] and BaBar [6]. Recently, LHCb presented updated results on the full angular distribution of the process B → (K * → Kπ)µµ from the data sample that corresponds to the total integrated luminosity of 3 fb −1 [2]. Further significant improvements in the precision of the measurements are expected in the ongoing LHC run and the LHCb upgrade [7], as well as future machines [8].
To fully exploit the forthcoming measurements it requires sufficient understanding of the long-distance backgrounds within the SM and/or the methods to disentangle them from the short-distance effects that might carry information about beyond-the-standard-model (BSM) physics. The non-perturbative QCD dynamics in the matrix elements of the local quark currents between the initial and final meson states is parameterised by hadronic transition form factors. The latter can be computed in the region of low hadronic recoil, which is the focus of this work, in the framework of Lattice QCD. Recent progress for B → K * and B s → φ transitions has been reported in [9].
Another important irreducible class of long-distance phenomena stems from the resonances that are induced by four-quark operators. The model-independent description of these effects, based on first theory principles, is currently not available and one needs to rely on models, which ideally can then be tested, i.e. compared to data. One such tool is the low recoil Operator Product Expansion (OPE), in which the effects of the non-local matrix elements of the four-quark operators can be computed in terms of local matrix elements in powers of 1/Q [10,11]. Here, the hard scale is provided by Q ∼ ( q 2 , m b ), where q 2 denotes the invariant mass of the dileptons which at low recoil is of the order of the b-quark mass, m b [12].
The QCD equations of motion can be used to derive the improved Isgur-Wise relations [13] between form factors, valid at leading order in 1/m b [14]. Together with the OPE, these relations imply universality, that is, independence on the polarization of the final state hadron, of the transversity amplitudes in the high-q 2 region [14]. This feature enables the construction of observables free of short-distance dependence assuming no significant right-handed currents [14,15]. One can then use these observables to extract ratios of form factors independently of the underlying short-distance physics to be used directly in SM tests [16,17]. Uncertainties due to next-to-leading order 1/m b -corrections to the universality relations turn out to be parametrically suppressed, at percent level [14,15].
In the region above thecc-threshold, the charm loop effects turn into the nonperturbative resonant spectrum in B → K ( * ) distributions, that shows up as peaks from narrow resonances J/ψ, ψ(2S) and "wiggles" for higher 1 −− states above theDD-threshold [18,19]. While the narrow resonances are removed by kinematic cuts and are not directly relevant at low recoil anyway, the wiggles, observed in B + → K + µµ decays [20], constitute a background not captured locally by the OPE. Since the local resonance structure is a non-perturbative effect, it is not revealed at any order of the perturbative OPE. The amount of duality violation was investigated in a toy model in Ref. [11]; one expects that the OPE gives a reasonably good description for binned observables. In view of the increasing precision it is therefore important to understand this quantitatively for given bin position and size.
To assess the performance of the OPE we exploit existing data on B → (K * → Kπ)µµ angular distributions [2] in different binnings [15 − 19] GeV 2 , [15 − 17], [17 − 19] GeV 2 , [15 − 16], . . . , [18 − 19] allowing to zoom in with resolution ∆q 2 = 4, 2 and 1 GeV 2 , respectively. The differential branching fraction is available in the two larger binnings only [21]. As we assume new physics at the electroweak scale or higher, a binning-related effect is due to resonances, not BSM physics. The plan of the paper is as follows: In Section II we give the effective Hamiltonian and B → K * (→ Kπ)µµ angular observables. In Section III we briefly review the low recoil OPE and the Krüger-Sehgal approach [18] modelling resonance distributions locally and to be used as a test-case against the OPE. Section IV is devoted to the details of such tests and gives results of a global fit for resonance parameters. In Section V we present the outcome of the global fit for the BSM Wilson coefficients and provide estimates of OPE uncertainties. We conclude in Section VI. Auxiliary information can be seen in three appendices.
We briefly review the effective Hamiltonian in Section II A and the basics of the B → K * (→ Kπ) angular observables in Section II B, respectively.

A. The effective Hamiltonian
We employ in this work the effective weak Hamiltonian description for b → s transitions where O i and C i denote the dimension-six operators and their Wilson coefficients, respectively. µ is an (arbitrary) renormalization scale and V ij are CKM matrix elements. We use the basis of the four-quark operators O 1,...6 introduced in Ref. [22], i.e. the so called CCM-basis Here, T a denote the generators of QCD and the sums are over active quark flavors q = u, d, s, c, b. The photon (gluon) penguin operators O 7 (O 8 ) and the semileptonic operators O 9,10 are given as with chiral projectors P L,R = (1 ∓ γ 5 )/2. The mass of the b-quark is the running mass in the M S scheme at the scale µ. We neglect the mass of the s-quark as well as the ones of the leptons and CKM-subleading contributions proportional to V ub V * us .

(6)
We adopt the definitions of the angles from [15], that is, θ is the angle between the µ − and the B in the rest frame of the muon pair, θ K is the angle between the kaon and the negative direction of flight of the B in the Kπ-rest frame and φ is the angle between the normals to the planes spanned by the Kπ and µ + µ − pairs in the rest frame of the B. The angular coefficients J i = J i (q 2 ) can be expressed in terms of transversity amplitudes A a i , i.e., the transition amplitudes with specified polarization of the final vector meson, i =⊥, 0, and the lepton pair, a = L, R, see Appendix A. Neglecting the mass of the leptons the number of independent angular coefficients is eight [24]. The angular distribution for the CP-conjugate decay, d 4Γ , can be obtained by replacing in J all angular coefficients J 1,2,3,4,7 → +J 1,2,3,4,7 and J 5,6,8,9 → −J 5,6,8,9 , whereJ i equal J i with the weak phases complex-conjugated [25].
We consider the observables F L , the fraction of the longitudinally polarized K * mesons, and the CP-averaged ratios The forward-backward asymmetry in the lepton angles can be identified as A FB = S 6 . Due to the different definitions of angles and normalization of the J j i the following relations to the conventions used by LHCb [2,26] hold Furthermore, endpoint relations apply, which are based on general grounds [27] and hold irrespective of the underlying electroweak model We consider B → K * µµ decays in the high-q 2 region above the peaking charmonium resonances in the OPE (Section III A) and a phenomenological data-driven test case (Section III B).

A. The high-q 2 OPE
At high q 2 one may exploit the presence of this hard scale to employ an OPE [10] to control quark-loop effects. The corresponding contributions can be absorbed into the effective coefficients of O 7,9 following [14] The functions F can be found in [28], while A, B and C are given in [29]. The function h(q 2 , m 2 q ) specifies the one-loop contributions to the vacuum polarization induced by the quarks and reads with w = 4m 2 q /q 2 , where m q denotes the quark's mass. In the limit of the massless quark one finds h(q 2 , 0) = 8 27 One can then employ the heavy quark expansion and the operator identities of the QCD to derive the improved Isgur-Wise relations between the (axial)-vector and tensor form factors [13]. A simple derivation is found in [14].
After applying these relations one finds that at leading order in 1/m b the transversity amplitudes are functions of the universal linear combinations of the Wilson coefficients C L,R [14], namely where κ = 1 + (α s /(3π)) ln(m 2 b /µ 2 ). Note that the above form of the transversity amplitudes follows from the universality of C eff 9 (q 2 ), i.e., its independence on the polarization of the final vector meson. This is a property of the high-q 2 OPE [13]. The 1/m b -corrections to these relations are parametrically suppressed [15]. The transversity form factors f ⊥, ,0 are defined as the following combinations of the standard form factors A 1,2 (q 2 ) and V (q 2 ) where λ = λ(q 2 , m 2 B , m 2 K * ) denotes the Källén function λ(a, b, c) = a 2 + b 2 + c 2 − 2(ab + ac + bc). The form factor f 0 (q 2 ) is proportional to A 12 (q 2 ) that has been directly computed in Lattice QCD [9].
To ease notation in the remainder of this work we frequently drop the explicit q 2 -dependence for transversity amplitudes, form factors, effective coefficients etc.

B. The Krüger-Sehgal approach
In this section we describe the method which aims at a local description of charm resonances in the high-q 2 region of B → K ( * ) decays using the charm contribution to the self-energy of the photon and a factorization ansatz. The full non-perturbative form of the charm vacuum polarization function can be extracted from data on the e + e − → h i scattering, where h i denotes all produced hadrons in the given kinematic region. This idea was first proposed in [18,30]. Such analysis was recently performed in Ref. [31] for the case B → Kµµ.
The experimentally accessible observable is the ratio of the cross section of e + e − scattering into hadrons normalized to the corresponding cross section of the scattering into muon pairs as the function of the center-of-mass energy s ≡ q 2 , namely We fit for the function R(s) in the interval √ s = 3.7 GeV to √ s = 4.8 GeV using the available data on the e + e − → h i processes from the BES experiment [34,36]. The ratio R(s) is the sum of the resonant and the continuum contributions The explicit form of R res, cont (s) with further details of the fitting procedure can be found in Appendix B. The charm contribution to Eq. (16) is extracted using where R uds = 2.16 is the asymptotic value of the light-quark contributions.
The relevant scattering amplitude can be written as Gauge invariance dictates the form of the photon's self-energy, Π µν (s) = (−g µν q 2 + q µ q ν )Π(s). The optical theorem relates the imaginary part of this amplitude to the total hadronic cross section, which implies For easier comparison with Krüger and Sehgal (KS) [18], we introduce Using Im[Π (µ) ] = 1/(12π) we obtain The charm polarization function h c (s) is defined in such a way to match the perturbative evaluations in Eq. (11), as Then Eqs. (23) and (22) imply Together with Eq. (18) we extract the imaginary part of the function h c (q 2 ) from the fit for the function R(s). We obtain the real part of h c (s) from its imaginary part using the subtracted dispersion relation where the arbitrary subtraction point s 0 and the lower limit of integration t 0 are convincingly chosen in the perturbative regime below the J/ψ-resonance peak and P denotes the principal part. The function h c is shown in Fig. 1. We proceed using a factorization ansatz and absorb Π (KS) into the charm contribution of the effective coefficient of O 9 . The corresponding B → K * matrix element, which includes B → K * (cc) → K * charmonium contributions, can be obtained by replacing the propagating resonances with the self-energy Π µν (q 2 ) Therefore, Here, we explicitly included terms that arise from the perturbative b-and light-quark contributions. a 2 is a combination of Wilson coefficients that accounts for the perturbative charm-loop To be specific, in this work we employ the value obtained at next-to-next-to-leading order (NNLO) at the b-mass scale, a 2 = 0.2 in the numerical analyses. (In the operator basis used in earlier works 3a 2 corresponds to C (0) [19].) Furthermore, we introduced in Eq. (27) a fudge function η c ≡ η c (K * j , q 2 ) that corrects for effects beyond factorization. In general η c is complex-valued and depends on the transversity state of the K * , j =⊥, , 0. For previous usage of fudge factors, see, e.g., [19,[31][32][33]35]. Note, that in principle, a dependence on the decay angles is possible as well: θ -dependence can arise from electromagnetic corrections, while θ K -dependence can arise from the K * beyond the narrow width approximation, however, both of these effects are neglected in this work.

IV. WIGGLES AND NON-UNIVERSALITY
Both wiggles in binned q 2 -distributions and non-universality would signal a breakdown of the OPE. We compare the predictions of the OPE (red curves and boxes with form factors from [9]) to data (black) in Figs. 2 and 3, zooming in from 2 GeV 2 bins (plots to the left) to finer resolution with 1 GeV 2 bins (plots to the right). Quite generally one expects an onset of resonance structure, consistent with the measured R-ratio [36], see also Figs. 1 and 8. As the branching ratio has not been measured with resolution smaller than 2 GeV 2 bins we only show this binning in Fig. 4. From these figures one cannot draw firm conclusions on observing a resonance structure in any of the observables due to the limited experimental precision. While currently resonance effects are not noticeable in the 2 GeV 2 bins, the alternating patterns in the 1 GeV 2 bins, however, may be hinting at such as structure. Further data with improved precision is required to clarify this point.
In addition, the data have to meet the endpoint predictions (9) irrespective of BSM contributions. A significant violation of Eq. (9) would, for instance, point to underestimated backgrounds other than from K * → Kπ. In particular with 1 GeV 2 bins data on S 3,4,5 are presently in mild conflict at ∼ 1 − 2σ with the endpoint relations. Note, however, that the endpoint bin is challenged by the dying statistics and needs to be viewed with a grain of salt, see Fig. 2.
In Section IV A we discuss different classes of angular observables according to their sensitivity to short-distance physics and resonance parameters. As wiggles and non-universality are both effects beyond the OPE, yet need to be measured, therefore, and only therefore, we use the phenomenological KS-approach as an efficient parameterization of local spectra. In Section IV B we work out phenomenological constraints on the resonance parameters.
A. Short-distance freedom and short-distance sensitivity The universal feature of the OPE-amplitudes (13) enables the construction of observables in which the dependence on the short-distance coefficients, C L,R , cancels [14]. F L , S 3 , S 4 belong to this class of short-distance free observables, The angular observables FL, S3 and S4 in the OPE for 2 GeV 2 bins (plots to the left) and 1 GeV 2 bins (plots to the right) shown as red boxes versus data (black) from LHCb [2]. Systematic and statistical uncertainties are added in quadrature. The light-shaded red bands illustrate the OPE for infinitesimal binning. Form factors are taken from [9]. The binned observables approach the continuous functions in the limit of infinitesimal bin width.
which are defined in terms of transversity amplitudes as follows, respectively, where the differential decay rate is given as Inserting the transversity amplitudes (13) into Eq. (29), one finds that the dependence on Wilson coefficients in the limit of the infinitesimal bin width cancels The measured observables correspond to the binned values of the angular coefficients, i.e., J i = bin J i dq 2 such that in (31) products of type f i f j , for i = 0, , ⊥ are integrated as bin ρ 1 f i f j dq 2 , where Since the effective coefficients that follow from the OPE are slowly varying functions of q 2 the resulting binning effect is small [16]. By the same argument which makes F L , S 3 , S 4 short-distance insensitive contributions to C eff 9 as in Eq. (27) with universal η c (K * j , q 2 ) drop out in these observables. The good agreement between the data and the OPE shown in Fig. 2 therefore implies that there are no extremely large contributions from non-universal pieces. It also implies constraints on right-handed currents, which could spoil Eq. (31). 2 Such BSM effects, however, would induce shapes essentially flat in q 2 . Another class of observables are short-distance dependent angular observables. These include A FB , S 6 and S 7, 8,9 . The former read in terms of transversity amplitudes Within the OPE (13) these observables reduce to with 2 In the presence of chirality-flipped operators beyond (4) the apparent universality of the short-distance coefficients following from the lowest order OPE, Eq. (13), breaks down to a partial one. Specifically, only the longitudinal and parallel coefficients remain the same. and ρ 1 can be seen in Eq. (32). A FB and S 6 are shown in the SM in Fig. 3, where we employ the SM Wilson coefficients (3), (4), evaluated at NNLO [37,38]. Universality predicts further J 7,8,9 = 0 [14], and consequently which can be explicitly seen from Appendix A.
The branching ratio, shown in Fig. 4 for the smallest available binning, depends on BSM physics and is highly sensitive to wiggles whether universal or not, as no cancellations as in the previously discussed observables can take place. Also in the branching ratio the OPE plus SM is in agreement within 1σ with the data [21]. We learn that in order to maximally probe for local structures and their deviations from binned OPE-results one has to simultaneously fit to BSM coefficients and resonance parameters.

B. Probing resonances
In this section we extract information from data on the η c (K * i , q 2 )-parameters. To begin we note that by means of Lorentz invariance all non-factorizable contributions have to vanish at the endpoint q 2 max [27]. This implies To facilitate a fit already with presently available data we assume constant η c -functions in the entire high-q 2 region. This is, of course, a simplifying working-assumption, however, as we show in Section V, it describes B → K * µµ data well. Similarly, there is support for this from B → Kµµ data, which also gives a good fit for this assumption [31]. With better data one should investigate more general shapes. Let us illustrate how S 7,8,9 are informative for non-universal η c [15], as within our assumptions, approximately, whereC eff 9 equals C eff 9 with the charm contribution removed. Terms quadratic in a 2 have not been spelled out explicitly; they require relative phases in the η c and are mildly suppressed by 3a 2 /C 9 .
The weakest constraints stem from the largest bin size. We stress that the constraints in the left plot of Fig. 5 are obtained without assuming the SM. The corresponding predictions of the fit for the BSM coefficients are discussed in the next section V. For illustrational purposes we show in the plot to the right of Fig. 5 the allowed 1σ contour for η 0 = η and η ⊥ for 1 GeV 2 bins in the SM. The constraints are more tight than in the model-independent fit. We also show the individual 1σ areas of various constraining observables, S 8 , S 9 (magenta area), A FB (blue area) and 2 GeV 2 bins for the branching fraction [21] (green area).
The sensitivity of A FB and the branching fraction to the resonance parameters is illustrated in Fig. 6. (B → K * q 2 -spectra including resonance effects have been given previously in [19,32,33], and recently in [31] using B → Kµµ data.) The sensitivity of S 5 is very similar to the one of A FB and not shown. We recall that in these observables already universal resonance effects do not cancel. Shown are local SM spectra for universal and constant η 0, ,⊥ = ±1. This choice is consistent with the measured B → J/ΨK * and B → Ψ(2S)K * branching ratios, see Appendix C. In addition, we show the impact of non-universality, η ,0 = 1, η ⊥ = −1 (dotted purple curve). The resulting spread for different η j is rather small above ∼ 15 GeV 2 except in the branching ratio, which could be used to detail the charm contribution locally, as in B + → K + µµ decays [20]. Ideally this should be done for each K * transversity state, dΓ j /dq 2 , see Eq. (30).

V. MODEL-INDEPENDENT ANALYSIS
The constraints on the BSM coefficients from the different fits using the three available q 2 -binnings are presented in Fig. 7. We recall that all fits are based on B → K * µµ data at low recoil only. The black dashed contours are from a simultaneous fit to Wilson coefficients and resonance parameters η 0 , η ⊥ as discussed in Section IV B. The red shaded areas are obtained within the OPE. In these plots, only two BSM coefficients are switched on at a time, that is, δC 9 , δC 10 in the upper plots and C 9 , C 10 in the lower plots (Formulae which include right-handed currents can be taken from [15].).
We find that within the OPE, as well as the local charm models, the SM agrees well with the data at the current level of precision. The findings are consistent with the pure low recoil analysis of Ref. [39]. Zooming in from large to small bins, the OPE result undergoes small changes, caused by the binning-dependent experimental uncertainties. With the local KS-model zooming in increases the resolution to charmonium contributions. In the fits for C 9,10 we 2) within the KS-model (simultaneously fitting for C 9,10 and η 0,⊥ ) for (1, 2, 4) GeV 2 bins, respectively. Similar results are obtained for the corresponding fits to C 9,10 . All plots exhibit consistency between local modelling and the OPE. We conclude that within current precision, charm effects appear to be controlled and do not endanger the validity of BSM constrains. To estimate the uncertainties of the OPE predictions for a given binning, we suggest to use the ratios where ρ 1,2 are given in Eqs. (32) and (35) and evaluated with the respective C eff 9 . The k are theory measures of the OPE's binning-related uncertainty. Their relation to the observables is straight-forward in the universality-limit.
The k are worked out in Table I, taking into account the 1σ ranges of η 0,⊥ , C 9,10 of the fit shown in Fig. 5 (41) for different q 2 -bins and 1σ ranges of parameters η0, η ⊥ and C9, 10. The coefficients C 9,10 are set to zero.
As expected, larger bins are better behaved than smaller ones, i.e., have k closer to 1, except for those near the endpoint. The [17 − 19] Table I mostly k < 1. This is driven by Re[h c (q 2 )] < Re[h(q 2 , m 2 c )] above q 2 ∼ 16.4 GeV 2 , see Fig. 1. If the resonance parameters η c would be determined more precisely, the uncertainty on the mismatch between the OPE and the KS-model would shrink. This way, the deviations | k − 1| given here correspond to upper limits.
Note the possible ambiguity that could arise if there is a significant constant or slowly varying contribution to C eff 9 stemming from non-resonant DD-backgrounds. Such effect might be differentiated from a new physics contribution to C 9 only in the case the latter is CP-violating or lepton flavor non-universal. Hadronic backgrounds not captured by the OPE for a given binning cause the KS-fit and the OPE-fit to disagree. We interpret this as an uncertainty of the OPE-fit to the Wilson coefficients. A possible binning-independent uncertainty will be one of the limiting factors to test the SM at low recoil.

VI. CONCLUSIONS
The low recoil region in semileptonic |∆B| = |∆S| = 1 decays is inhabited by wider charm resonances, locally not captured by the OPE. While with current data such resonance patterns may be only at the border of being visible in B → K * µµ decays, see Figs. 2 and 3, in the near future this will be an important background to SM precision tests. In order to understand these effects, which are inaccessible from within the OPE, we use the KS-model [18], cf. Section III B, which does describe resonances locally, as a test-case against the OPE .
Using available B → K * µµ data at low recoil, we performed simultaneous fits to resonance parameters of the KS-model and BSM Wilson coefficients, and compare it to the plain OPE-fit. We find that the resulting constraints are consistent with each other, and consistent with the SM, see Fig. 7. There is room left for sizeable BSM contributions. Let us emphasize the difference between our work and the recent global fits [39][40][41] that used the experimental data for all b → s processes as an input but included only the total low recoil bin for B → K * . We focus solely on the low recoil region of this decay mode and use all available data for this kinematic region including also the smaller bins, since our goal is to scrutinize the local q 2 -shape. Specifically, in the KS-model fit we do not fix the hadronic parameters expressed by the fudge factors but use them as fit parameters together with the short-distance Wilson coefficients.
To estimate the uncertainties of the OPE for a given binning, we use the coefficients k defined in Eq. (41), with current evaluations shown in Table I In the future more precise data with even smaller binning than currently exist would be desirable to determine the resonance parameters more accurately. This will directly influence the estimates in Table I, which measure not only the mismatch between local spectra and the OPE, but include also uncertainties within the KS-model. Refinements of the method, such as less minimal parameterizations for the η c (K * j , q 2 ) and h c (q 2 ) functions, may also be envisaged. As dominant uncertainties within the OPE are due to hadronic form factors, improving their predictions would be desirable, too. 9 16 where we neglected terms proportional to m 2 /q 2 . The full expressions can be found in Ref. [14].
Here we describe the fitting procedure for the ratio R(q 2 ) defined in Eq. (16) using the experimental input on the e + e − → h i cross section from the BES-II experiment [36]. Above theDD threshold the four wide charmonium resonances, ψ(3770), ψ(4040), ψ(4150) and ψ(4415), with quantum numbers J P C = 1 −− , appear in the spectrum. We adopt the fitting procedure from [36] and model the background according to [31].
The transition amplitude of the resonance r into a final state f is modelled by the Breit-Wigner ansatz with a phase δ r and mass m r : Since only three relative phases carry a physical information, we set δ ψ(3770) = 0. The s-dependent decay widths of r → f are given by the formula Here,Γ r denotes a fit parameter specified for every given resonance and Z f ≡ ρP f , whereρ 1GeV. The sum over the orbital angular momenta of the decaying final states and the energy dependent partial wave function B L are given in [36]. The momentum P f of the two body decay of the resonance f into the final mesons with masses m 1 and m 2 is given by the familiar formula The different decay channels are given in Table II. The total hadronic width is the sum over all final states The total width of a charmonium resonance is the sum of the hadronic and leptonic widths  Lepton universality Γ r→e + e − = Γ r→µ + µ − is assumed for the leptonic decay widths of the electron and the muon. The kinematic suppression factors are included for the decay rates of the resonances that involve tau pairs in the final state. The total square of the modulus of the inclusive amplitude of the resonances is the following incoherent sum over the final states: The resonance contribution R res (s) to R(s) (see Eq. (17)) is given by The continuum background is modelled as [31] R with x = 4m 2 D /s. For the light quark ratio R uds and the one including charm, R udsc , we use the predictions from Refs. [42,43]. Specifically, we employ R uds = R(s = (3.73 GeV) 2 ) and R udsc = R(s = (4. The values of the parameters in (B9) used in the fit are given in the Table III. Above the fit interval in the open charm region, we use the Schwinger's O(α s ) result, in the form adopted from [31], where β(s) = 1 − 4m 2 c /s. The final result for h c (q 2 ) is given in Fig. 1  Here we give a brief overview on phenomenological data on B → ψ i K ( * ) decays, where ψ i denotes a generic charmonium 1 −− resonance, within factorization. The matrix element of the B → ψ i K ( * ) decays can be written as where the commonly used color singlet and octet operators read respectively. The Wilson coefficients of these operators read in terms of the ones in the CCM-basis (Eq. (3)) as C (1) = a 2 , given in Eq. (28), and C (8) = 1 3 (−C 1 + 6C 2 ). We employ the value of a 2 at NNLO [37,38] thereby including perturbative corrections to the weak b → ccs vertex.
Assuming factorization, the B → ψ i K ( * ) matrix element reads For κ = 1 this ansatz represents the naive factorization approximation (NFA). Note that the color octet operator does not contribute, e.g., [44]. The dependence on the renormalization scale µ does not cancel between the Wilson coefficients and the effective operators in the ansatz (C3). For further aspects of the factorization ansatz the reader is referred to [45]. The matrix element in factorization can be expressed in terms of charmonium decays constants, which can be extracted from data on Γ(ψ i → ) and form factors, to be evaluated at q 2 = m ψ 2 where η ( ) denotes the K * (ψ i ) polarization vector, k, p the 4-momenta of the K ( * ) ,B mesons, respectively, and q = p − k. Due to · q = 0 the terms proportional to q µ do not contribute. We use f J/ψ = 0.416 ± 0.006 GeV, f ψ(2S) = 0.297 ± 0.003 GeV and f ψ(3770) = 0.100 ± 0.004 GeV [46]. There are several processes for which the measurements reveal deviations from NFA. For instance, the branching fraction of the process B − → K − χ c0 has been observed to be significantly non-vanishing, while the corresponding factorization contribution vanishes due to parity conservation of QCD. As a second example, the branching fraction of B → KJ/ψ also deviates from its NFA value, to be discussed in Section C 1. A possible source of non-factorizable corrections are B-meson decays to D ( * ) s D ( * ) pairs which afterwards re-scatter into Kψ i pairs. The analysis of such effects was undertaken in Ref. [48]. It contains significant theoretical uncertainty related to the vague knowledge of the relevant strongly coupled meson vertices. In the following sections we follow the phenomenological point of view and extract the fudge factors for the processes B → K ( * ) (J/ψ, ψ(2S)), in order to gain some further insights.
Here, we used the form factor f + (m 2 ψi ) evaluated in Lattice QCD in [49], see also [50]. Form factors extrapolated to q 2 = m 2 ψi from light cone sum rules from [35] yield very similar results. The value for ψ(3770) is only given for completeness as this resonance is included in the fit to BES-data. The values in Eq. (C11) reveal an order one deviation from naive factorization, κ = 1. Note, that the corresponding signs (phases) remain undetermined from this extraction. Assuming universal fudge factors for all polarizations of the K * , where κ ψi K * = η c (K * , m 2 ψi ), the B → ψ i K * branching fraction can be written as where and the ratios of the B → K * form factors Here, we employed the B → K * form factors from [9]. Eqs. (C11) and (C17) suggest that non-factorizable corrections for processes involving a K * are smaller than those with a K. Further information on universality can be obtained from data on the polarization fractions of the K * in B → K * (J/ψ, ψ(2S)) decays [27,51,52] |A ⊥ | 2 = 2 c 2 y 2 r , |A | 2 = 2 r , r ≡ (a − bx) 2 + 2(1 + c 2 y 2 ) , which are normalized as |A 0 | 2 + |A | 2 + |A ⊥ | 2 = 1. We compare data to the factorization predictions in Table IV for form factors extrapolated from the lattice [9] and a phenomenological fit (SE2LEL) [17]. Results based on [35] have larger uncertainties and are consistent with both theoretical predictions.
The factorization predictions work for the perpendicular polarization fraction, but exhibit larger spread and uncertainties in the other two, in particular for the J/ψ final state. The spread in theory predictions points to the sensitivity to form factor predictions, which in this intermediate q 2 region need to be extrapolated which brings in additional uncertainties.