Published for SISSA by Springer Received: February 19, 2019 Revised: July 6, 2019 Accepted: July 8, 2019 Published: July 17, 2019 What RK and Q5 can tell us about New Physics in b → s`` transitions? JHEP07(2019)096 Marcel Algueró,a,b Bernat Capdevila,a,b Sébastien Descotes-Genon,c Pere Masjuana,b and Joaquim Matiasb,c a Grup de Fı́sica Teòrica (Departament de Fı́sica), Universitat Autònoma de Barcelona, E-08193 Bellaterra, Barcelona, Spain b Institut de Fı́sica d’Altes Energies (IFAE), The Barcelona Institute of Science and Technology, Campus UAB, E-08193 Bellaterra, Barcelona, Spain c Laboratoire de Physique Théorique, UMR 8627, CNRS, Univ. Paris-Sud, Université Paris-Saclay, 91405 Orsay Cedex, France E-mail: malguero@ifae.es, bcapdevila@ifae.es, sebastien.descotes-genon@th.u-psud.fr, masjuan@ifae.es, matias@ifae.es Abstract: The deviations with respect to the Standard Model that are currently observed in b → s`` transitions, or B anomalies, can be interpreted in terms of different New Physics (NP) scenarios within a model-independent effective approach. We identify a set of internal tensions of the fit that require further attention and whose theoretical or experimental nature could be determined with more data. In this landscape of NP, we discuss possible ways to discriminate among favoured NP hypotheses in the short term thanks to current and forthcoming observables. While an update of RK should help to disentangle the type of NP we may be observing (Lepton-Flavour Universality Violating and/or Lepton Flavour Universal), additional observables, in particular Q5 , turn out to be central to determine which NP hypothesis should be preferred. We also analyse the preferences shown by the current global fit concerning various NP hypotheses, using two different tools: the behaviour of the pulls of individual observables under NP scenarios and the directions favoured by approximate quadratic parametrisations of the observables in terms of Wilson coefficients. Keywords: Beyond Standard Model, Effective Field Theories, Heavy Quark Physics ArXiv ePrint: 1902.04900 Open Access, c The Authors. https://doi.org/10.1007/JHEP07(2019)096 Article funded by SCOAP3 . Contents 1 Introduction and motivation 1 2 Inner tensions of the global fit 4 2.1 RK ∗ in the first bin 4 2.2 Bs → φµ+ µ− versus B → K ∗ µ+ µ− 5 2.3 Tensions between large and low recoil in angular observables 7 JHEP07(2019)096 3 Potential of RK (and Q5 ) to disentangle NP hypotheses 8 3.1 Global fits 10 3.2 LFUV fits 14 4 Pulls of individual observables 15 4.1 Observables with a large pullobs i within the Standard Model 19 obs 4.2 Observables with no pulli with respect to the Standard Model 21 5 Sensitivity of observables to different Wilson coefficients 21 5.1 Observables and parametrisations 22 5.2 Observables as conics 23 5.3 Directions favoured by the anomalies 25 6 Conclusions 27 A Dependence of best-fit points on hRK i[1.1,6] and hQ5 i[1.1,6] 28 1 Introduction and motivation The LHCb [1] and Belle [2] collaborations have recently updated the measurements of the Lepton-Flavour Universality Violating (LFUV) ratios RK and RK∗ : B(B → K(∗ )µ+ µ− ) RK(∗ ) = (1.1) B(B → K(∗ )e+ e− ) for various bins in the dilepton invariant mass. The LHCb collaboration observed hints of deviations from the SM between 2 and 3σ in these observables, whereas the Belle mea- surements for RK ∗ , affected by large uncertainties are compatible with the LHCb mea- surements as well as with SM. These deviations as well as others measured in b → sµµ observables have triggered the combined analysis of LFUV and lepton-flavour dependent (LFD) observables performed in ref. [3]. The global analysis presented in the aforementiond reference shows that, using a model-independent approach, the Standard Model (SM) hy- pothesis is disfavoured compared to various hypotheses of New Physics (NP) contributions –1– in b → s`` decays, with pulls w.r.t. the SM ranging from 5.3σ to 5.9σ. Similar results were obtained by other groups using different treatments of hadronic uncertainties and sets of observables [4–11]. These model-independent analyses constrain NP scenarios expressed as contributions to the short-distance Wilson coefficients Ci` in the effective Hamiltonian approach for b → s`` transitions. The first point to address is obviously whether NP has been discovered, but once this is established, it will prove important to determine the specific pattern of NP discovered. Indeed, even if the amount of data obtained up to now for b → sµµ makes sophisticated global fits to several Wilson coefficients possible [3–11], the outcome is still not conclusive enough to draw definite conclusions about the actual pattern of NP. Disentangling the JHEP07(2019)096 realized pattern is an essential guide to build NP models in agreement with these obser- vations. It is therefore usual to limit NP contributions to a few Wilson coefficients, that from now on we will refer as hypotheses, and to build NP models in agreement with these assumptions of the global fits. Most scenarios discussed in the literature assumed that there is NP in muons only, i.e. the LFUV-NP contributions come from allowing the presence of NP in the muon channel and not in the electron one (or it is considered small). Three particularly interesting one- dimensional scenarios have emerged, namely NP in C9µ , in C9µ = −C90 µ and in C9µ = −C10µ , with a larger significance for the first two scenarios and a smaller one for the latter in ref. [3]. On the contrary, we also found that a fit restricted to a subset of mainly LFUV observables exhibits a marginal preference for the C9µ NP = −C NP scenario compared to 10µ the other scenarios. Moreover, we found several interesting two-dimensional scenarios, all involving C9µNP , noticing that the inclusion of small but non-vanishing right-handed currents helps in solving some of the observed deviations with respect to SM expectations. Recently, several works allowed for NP also in the electron channel, but no particular structure was envisaged from these fits by simply taking some of the electronic Wilson coefficients different from zero [7, 8, 11–13]. However, in a recent article [14], we allowed the possibility of a specific structure, namely, that the b → s`` transitions get a common Lepton Flavour Universal (LFU) NP contribution for all charged leptons (electrons, muons and tau leptons). This permitted us to identify new favoured NP hypotheses. This idea was implemented by allowing two NP contributions inside the semileptonic Wilson coefficients: NP V Ci` = Ci` + CiU (1.2) with ` = e, µ, τ and where Ci` V stands for LFUV-NP and C U for LFU-NP contributions. i We distinguished the two contributions by imposing that Cie V = 0. It is important at this point to emphasize the difference between simply allowing the presence of NP also in electrons or allowing the existence of two different kinds of NP contributions (LFU and LFUV). The case of simply allowing NP in the electron channel has been discussed quite extensively in refs. [11, 13] (see also refs. [8, 12] for a smaller subset of scenarios with and without including low-recoil observables). However, our approach of distinguishing LFU- and LFUV-NP structures provides new ideas to model building and extends the possible interpretations of the current fits. Performing the fits with this new setting [3], we obtained our previous results in ref. [4] but also new scenarios different from refs. [8, 11–13]. This –2– can be seen by translating LFU and LFUV contributions into NP contributions to muons and electrons (leaving τ aside at this stage) NP V C9µ = C9µ + C9U , NP C10µ V = C10µ U + C10 , NP C9e = C9U , NP C10e U = C10 . (1.3) This seemingly innocuous redefinition yields interesting consequences, as discussed in ref. [14]. It opens interesting perspectives to explain with different mechanisms the anoma- lies coming purely from the muon sector (like hP50 i[4,6] ) and the ones describing the violation of lepton flavour universality (like hRK i[1.1,6] ). These results were later updated in ref. [3]. Let us notice in particular that this approach is different from all the analyses including JHEP07(2019)096 NP in electrons [7, 8, 11–13] where the muonic NP contribution is not correlated in any way with the electronic one. When translated from one language to the other, the most interesting one- or two- dimensional scenarios in refs. [3, 14] become, for the purely LFUV cases: V NP [Hyp. I] {C9µ } →{C9µ } V V NP NP [Hyp. II] {C9µ = −C10µ } →{C9µ = −C10µ } V [Hyp. III] {C9µ = −C9V0 µ } →{C9µ NP = −C9NP 0µ } V V NP NP [Hyp. IV] {C9µ , C10µ } →{C9µ , C10µ } V [Hyp. V] {C9µ , C9V0 µ } →{C9µ NP NP , C 90 µ } V V NP NP [Hyp. VI] {C9µ , C10 0 µ } →{C9µ , C100 µ } V [Hyp. VII] {C9µ = −C9V0 µ , C10µ V V = C10 NP NP 0 µ } →{C9µ = −C90 µ , C10µ = C100 µ } V [Hyp. VIII] {C9µ , C9V0 µ = −C10 V NP NP 0 µ } →{C9µ , C90 µ = −C100 µ } (1.4) and the scenarios allowing both LFUV and LFU contributions V V [Hyp.IX] {C9µ = −C10µ , C9U = C10 U NP } →{C9µ NP = −C10µ NP NP + 2C9e NP , C9e = C10e } V [Hyp.X] {C9µ , C9U } →{C9µ NP NP , C9e } V V [Hyp.XI] {C9µ = −C10µ , C9U } →{C9µ NP NP = −C10µ NP NP + C9e NP , C10µ , C9e } V V U NP NP NP NP NP [Hyp.XII] {C9µ = −C10µ , C10 } →{C9µ , C10µ = −C9µ + C10e , C10e } V U NP NP NP [Hyp.XIII] {C9µ , C10 } →{C9µ , C10µ = C10e } V U NP NP NP [Hyp.XIV] {C9µ , C10 0 } →{C9µ , C100 µ = C100 e } (1.5) The hypotheses VII and VIII correspond to the hypotheses 1 and 5 in refs. [3, 4]. The hypotheses IX to XIV correspond to the scenarios 6 to 11 in refs. [3, 14] (scenarios 5 and 13 are also interesting in terms of their ability to explain the deviations observed, but they require three or four free parameters and will not be considered in the following). In this situation, it becomes clear that new data will be instrumental to disentangle the different hypotheses. The goal of the present article is to scrutinize the results of the fit from a different perspective to prepare the next step, i.e. to discriminate the most relevant NP scenario among the ones already favoured, complementing our previous works [3, 4, 14, 15]. Currently, the most significant patterns identified exhibit a pull w.r.t. the SM very close –3– to each other (within a range of half a σ). We explore strategies to disentangle different scenarios and to identify the impact of a more precise measurement of RK . We then combine information on RK and Q5 in order to illustrate that RK by itself will not be sufficient to disentangle clearly one or a small subset of scenarios, but that a combination of RK and Q5 can be useful, depending on the (future) measured value.1 In section 2 we discuss the inner tensions of the fit in order to point those observables where further experimental or theoretical work would be required. In section 3 we explore how a forthcoming precise measurement of RK can disentangle or disfavour scenarios as- suming that the statistical error is reduced and the central value stays within 2σ of its JHEP07(2019)096 present value. We analyse it considering two different fits, either with all observables or only the LFUV subset. We also discuss the impact of a measurement of Q5 in relation with its possible measurement by Belle II and LHCb. We then discuss the structure of the current fits, looking more closely at the deviations of some observables in section 4, focus- ing on the change in their pulls depending on the NP scenario considered. In section 5 we discuss the structure of the observables in terms of their Wilson coefficients to determine their sensitivities and the directions preferred by each of the anomalies, before drawing our conclusions. 2 Inner tensions of the global fit In refs. [3, 4], we saw that different NP scenarios involving C9µ NP led to a much better description of the data than the SM, with fits reaching p-values around 60–80% (the SM being around 9%) and providing pulls with respect to the SM above 5σ. The overall agreement is thus already very good within these NP scenarios and from a purely statistical point of view, it should be expected that these fits exhibit slight tensions. It is however interesting to look at these remaining tensions in more detail in order to determine where statistical fluctuations may be reduced with more data or where improved measurements might help to lift the degeneracy among NP scenarios. We focus on three main tensions that we consider particularly relevant in the current global fit. 2.1 RK ∗ in the first bin A first tension related to RK ∗ occurs in the global fit and it proves interesting to consider both RK ∗ and B(B → K ∗ µ+ µ− ) as measured by LHCb in order to understand its nature (see figure 1). Let us first consider the second bin (from 1 to 6 GeV2 ) for RK ∗ . Even though the deficit could be consistent with an excess in the electron channel with respect to the muon one, the study of the corresponding bins of B(B → K ∗ µ+ µ− ) points towards a deficit of muons. The mechanism that explains the deviation with respect to the SM in the long second bin of RK ∗ is consistent with all the deviations that have been observed in other channels and different invariant di-lepton mass square regions. 1 Up to now only the Belle experiment has been able to perform a measurement of Q5 , leading to hQBelle 5 i[1,6]= +0.656 ± 0.485 ± 0.103 [16]. –4– JHEP07(2019)096 Figure 1. RK ∗ (left panel) and B(B → K ∗ µ+ µ− ) (right panel) measured by LHCb. Figures extracted from refs. [23] and [24] respectively. The situation is different for the first bin of RK ∗ , where B(B → K ∗ µ+ µ− ) is clearly compatible with the SM (see figure 1). An excess in the electron channel would then be needed in order to explain the observed deficit in hRK ∗ i[1.1,6] . This difference of mechanism between the first and the second bins of RK ∗ can be understood in two ways: i) a specific NP effect [17, 18] localised at very low q 2 and able to compete with the dominant Wilson coefficient C7 (well determined to be in agreement with the SM expectations from B(B → Xs γ)) [3, 4, 19–22]; ii) some experimental issue in measuring di-electron pairs at very small invariant mass, close to the photon pole. It would be very interesting that LHCb keep on their efforts to understand the systematics in this bin. Interestingly the recent Belle measurement [2] indicates also a low central value in the same bin, even though the large uncertainty affecting the measurement prevents us from drawing any definite conclusion and makes it compatible also with the SM. Another approach to slightly reduce the tension between data and SM in the first bin of RK ∗ through a NP explanation consists in including NP contributions to the b → see channel, in particular, considering right-handed currents affecting electrons, as discussed in ref. [7]. In the scenarios S8-S11 (using the notation of ref. [7]) the prediction of hRK ∗ i[0.045,1.1] is found to be within ∼ 1σ range of the current measurement. This could open a new window to explore the existence of right-handed currents and to explain some of the tensions found, even though more data is required in order to be conclusive. 2.2 Bs → φµ+ µ− versus B → K ∗ µ+ µ− Another tension in the fit concerns the branching ratio for Bs → φµ+ µ− , in particular when compared with the related decay B → K ∗ µ+ µ− . The prediction for the branching ratio B(B → K ∗ µ+ µ− ) involves hadronic form fac- tors to be determined using different theoretical approaches, depending on the di-lepton invariant mass region analysed: at large recoil, one can use light-cone sum-rules based on light-meson distribution amplitudes [25], while lattice form factors are available at low recoil. Due to the difficulty to assess precisely the uncertainties attached to light-cone sum rules, we perform our computation using more conservative results from light-cone sum –5– rules based on B-meson distribution amplitudes [26] with conservative error estimates, ex- ploiting QCD factorisation to restore correlations that were not available in ref. [26]. We checked that our results are compatible with those obtained in ref. [25] and that the two approaches yield very similar results for the fits [4, 6, 27, 28]. A recent update of these form factors is available in ref. [29] using the same approach as ref. [26], adding corrections from higher twists and providing correlations. We will update our results accordingly in a coming publication, but we do not expect very significant changes for the present article, based on our previous studies [15]. For instance, we checked that even if a large reduction of 50% is achieved on the error of the form factors, the resulting uncertainty of key optimized observables like P50 is minor (it would imply a reduction from JHEP07(2019)096 10% to 8% for the anomalous bins of P50 ). On the contrary, a large impact is observed in unprotected observables like branching ratios or Si observables. As our fit is driven by the optimized observables, we expect only minor changes in the outcome of the fits. Contrary to the case of B(B → K ∗ µ+ µ− ), there are no computations available using the B-meson light-cone sum rules of refs. [26, 29] for Bs → φµ+ µ− , and one must rely on the estimates given in ref. [25]. One can see in figure 2 that at low recoil, where lattice form factors are used, the prediction for B(B → K ∗ µ+ µ− ) is expected to be slightly larger than B(Bs → φµ+ µ− ) and indeed data (with large error bars) follows the same trend. On the contrary, in the large-recoil region where the light-cone sum rules results of ref. [25] are used, the SM predictions lead to a larger value for B(Bs → φµ+ µ− ) than for B(B → K ∗ µ+ µ− ). Surprisingly, data shows the opposite trend, which may come from a statistical fluctuation of the data leading to an inversion of the experimental measurements of both modes at large recoil. Alternatively, this issue may signal a problem in the theoretical prediction of the form factors of ref. [25]. Firstly, these predictions are obtained by combining results in different kinematic regions (light-cone sum rules and lattice QCD) which do not fully agree with each other when they are extrapolated: the fit to a common parametrisation over the whole kinematic space leads to a fit with uncertainties that may be artificially small due to these incompatibilities of the inputs. Moreover, the choice of the z-parametrisation [25, 26] used to describe the form factors over the whole kinematic range has interesting properties of convergence, but it may in some cases lead to potential unitarity violations [30]. Finally, another issue that specifically affects B(Bs → φµ+ µ− ) is the Bs -B̄s mixing. As it is well known, Bs -B̄s mixing implies that the time evolution of the Bs meson before its decay involves two mass states with different widths that are linear combinations of the flavour states Bs and B̄s . The current measurements performed at LHCb are integrated over time, and the neat effect of the evolution between the two mass states is a correction of O(∆Γs /Γs ) in the relation between the theoretical computation of the branching ratio and its measurement [31–34]. This effect is taken into account in the global fit [6] as an additional source of uncertainty for the theoretical estimate of the branching ratios. The experimental efficiencies should also be corrected for this effect, which depend on the CP-asymmetry A∆Γ that can also be affected by NP contributions. It should thus be kept free within a large range in the absence of measurements. Neglecting this effect and assuming a SM value for this asymmetry may lead to an underestimation of some systematics on the efficiencies. For instance, ref. [35] showed that this issue can lead to –6– 0.8 SM from B→K * μ+ μ- 0.7 SM from Bs →ϕ μ+ μ- 0.6 0.5 JHEP07(2019)096 0.4 0.3 0.2 q (GeV2) 0 5 10 15 20 2 Figure 2. Theoretical predictions for B(Bs → φµ+ µ− ) and B(B → K ∗ µ+ µ− ) within the SM along with their corresponding experimental measurement. The results at large recoil are presented here only for illustrative purposes and are based on the form factors presented in ref. [25] (these results are not used in our global analyses). The results at low recoil are indeed used in ref. [6] and are based on available lattice QCD inputs for the form factors. an additional systematic effect of 10% in the Bs → µ+ µ− systematics. The impact on efficiencies from NP effects was indeed considered in ref. [36] for Bs → φµ+ µ− by varying C9µ in the underlying physics model used to compute signal efficiencies, leading to a much smaller effect in this case (of a few percent, in line with back-of-the-envelope estimates). 2.3 Tensions between large and low recoil in angular observables We discuss for the first time here a rather different type of tension, concerning the B → K ∗ µ+ µ− angular observables at large and low recoil. On the one hand, we observe that branching ratios exhibit the same discrepancy pattern between theory and experiment at low and large recoil.2 On the other hand, the current deviations at LHCb in P5 require NP contributions with opposite sign in the two kinematic regions. Indeed, the pull between the SM value and the LHCb experimental measurement in P5 [15,19] has the opposite sign (albeit the significance is only 1.2σ) w.r.t. its large-recoil bins, in particular P5 [4,6] and P5 [6,8] . This very slight tension is not there in the case of the Belle data where same-sign deviations are observed, even though the error bars are rather large in this case. For the purposes of illustration, let us consider the NP scenario where there is no LFU contribution and NP occurs only in C9µ V and C V . This is illustrated in figure 3 where 10µ 2 This is true for all b → s modes, apart from the decay Λ0b → Λµ+ µ− , where the experimental errors at low recoil are very large and the normalisation chosen prevents further interpretation [37, 38]. –7– JHEP07(2019)096 Figure 3. 68.3% (left) and 95% (right) CL solutions regions for the observables discussed in the V V main text in the (C9µ , C10µ ) plane. The yellow region corresponds to the overlap region. hP2 i[15,19] is only shown in the left panel. the constraints for these observables (as well as other relevant observables that will be listed below) are shown at 68.3% (left) and 95% (right) CL. One can notice their milder V . hP 0 i sensitivity to C10µ V 0 5 [4,6] (blue region) would prefer a negative C9µ while hP5 i[15,19] (green region) would favour a positive C9µ V at 68.3% CL. Black dots indicate the particular solutions (−1.02, 0) and (−0.45, 0.45) corresponding to the best-fit points of the 1D favoured scenarios in ref. [3] (hypotheses I and II of the present article). We also indicate the constraints from hP2 i[4,6] , hP2 i[15,19] , and hRK i[1.1,6] , hB(B 0 → K ∗0 µ+ µ− )i[15,19] since we believe that they are representative of the set of observables driving our global fit.3 The former pair of observables (hP2 i[4,6] , hP2 i[15,19] ) has a large overlap region compatible with the SM while the latter one (hRK i[1.1,6] , hB(B 0 → K ∗0 µ+ µ− )i[15,19] ) overlaps far from the SM point. While P50 and RK strongly constrain NP solutions, the P2 bins are weakly constraining. Finally, the yellow region in the right panel in figure 3 is the overlap of the regions from the five observables obtained after considering the data regions at 95% CL. In summary, an interesting tension between low- and large-recoil regions for P50 is observed at the 2-sigma level, favouring C9µ contributions of different signs in the two kinematic regions. Although not statistically significant, this inner tension seems to require either different sources of NP or a shift in the data once more statistics is added. 3 Potential of RK (and Q5 ) to disentangle NP hypotheses In this section we discuss the potential impact of the prospective measurements of hRK i[1.1,6] and hQ5 i[1.1,6] on the global fits in order to distinguish NP hypotheses. We perform the following illustrative exercise: we vary the experimental values of hRK i[1.1,6] P1 and P40 observables are known to behave in a more SM-like way than the ones selected here, thus 3 providing weaker constraints. –8– and hQ5 i[1.1,6] within suitable ranges, and we perform fits according to these values taken as actual measurements. First only hRK i[1.1,6] is allowed to vary before we consider the combined impact of hRK i[1.1,6] and hQ5 i[1.1,6] . The ‘pseudo-data’ for hRK i[1.1,6] takes into account the expected increase in statistics soon available for this observable. For this exer- cise we take as the prospective experimental error for hRK i[1.1,6] +0.044: this corresponds √ to a reduction by 2 of the statistical uncertainty in ref. [1]. Indeed, according to the latter reference, this would amount to the inclusion of the data sets of 2017 and 2018 which are said to have the same statistical power as the combined data set of Run 1, 2015 and 2016. For each fit (corresponding to a given hypothesis and set of data), both the pull of the hypothesis w.r.t. the SM (PullSM ) and the best-fit-point (b.f.p) are computed, which JHEP07(2019)096 we plot as functions of either hRK i[1.1,6] or hQ5 i[1.1,6] . Before discussing the results of our analysis, we first state our assumptions: • We follow the same approach as in refs. [3–6, 14]. We consider a set of b → s`` observables measured by different experiments, we determine the experimental and theoretical correlation matrices between these observables assuming Gaussian distri- butions. Under a given NP hypothesis (generally described by one or two parameters added to specific short-distance Wilson coefficients), we build a χ2 from these observ- ables and their correlation matrices that is used to extract the best-fit points and the confidence intervals of the NP parameters, as well as the pull of the NP hypothesis with respect to the SM, within a frequentist framework. • We consider two different kinds of fits with different subsets of observables [3]: on one side, the global fit (or Fit “All”, to all 178 available observables) and on the other one, the LFUV fit, where only the observables measuring LFUV are included (plus constraints coming from radiative decays, leading to 20 observables). When several experiments have measured the same observable, we do not average the results but we include all these measurements in the χ2 taking into account their (theoretical) correlations. • Any variation of the experimental value of hRK i[1.1,6] could manifest itself also in a change in the branching ratios B(B + → K + µ+ µ− ) and/or B(B + → K + e+ e− ). However, the update of hRK i[1.1,6] in ref. [1] has not led to significant changes in these branching ratios, and we will assume that this will also occur in the forthcoming updates, so that we modify only the value of hRK i[1.1,6] • hRK i[1.1,6] is freely varied within a 2σ range from its current experimental value. It represents a good compromise between a high coverage of the true value and a span compatible with our computational means. hQ5 i[1.1,6] is varied within the range [−0.5, 1.0] in order to ensure that we scan over values corresponding to the most relevant NP scenarios (see figure 2 of ref. [14]). • With the increased statistics available at Run 2, it will be possible for experiments to provide more precise determinations of key observables. Therefore, besides the reduction in the error of hRK i[1.1,6] , we assume a guesstimated uncertainty of order 0.1 for hQ5 i[1.1,6] . –9– Global Fits Global Fits 8 8 NP I: C9µ NP NP II: C9µ = −C10µ NP 7 III: C9µ = −C9µ 7 PullSM (σ) PullSM (σ) 6 6 5 5  NP NP IV: C9µ , C10µ  NP  V: C9µ , C 9 µ  NP  JHEP07(2019)096 VI: C9µ , C10µ 4  NP NP  4 VII: C9µ = −C9µ, C10µ = C10µ  NP  VIII: C9µ , C9µ = −C10µ 0.75 0.80 0.85 0.90 0.95 0.75 0.80 0.85 0.90 0.95 RK [1.1,6] RK [1.1,6] Global Fits 8 7 PullSM (σ) 6 5  V  V IX: C9µ = −C10µ , C9U = C10 U  V U X: C9µ, C9  V V  XI: C9µ = −C10µ , C9U 4  V  V U XII: C9µ = −C10µ, C10  V U XIII: C9µ , C10  V U XIV: C9µ , C10 3 0.75 0.80 0.85 0.90 0.95 RK [1.1,6] Figure 4. Global fit: impact of the central value of RK [1.1,6] on the PullSM of the NP scenarios under consideration. The purpose of this analysis is not to provide precise determinations of the pull of the SM and the b.f.p.s for different values of RK [1.1,6] and Q5 [1.1,6] but rather to gain qualitative knowledge on how experimental measurements of these two observables will drive the analyses. This is particularly true for the b.f.p. plots that will provide only the central value but not the confidence intervals obtained for the NP contributions. 3.1 Global fits Figure 4 displays the outcome of the global fit (or fit “All”, involving 178 observables) for the pulls with respect of the SM, assuming different experimental central values of RK [1.1,6] and different NP hypotheses varied according to the procedure described above. In appendix A, figure 13 displays a similar result for the best-fit point of each NP hypothesis. The shaded vertical band in the plots of figures 4 and 13 highlights the current experimental 1σ confidence interval for the LHCb average of RK [1.1,6] . – 10 – Figure 4 illustrates the relevance of hRK i[1.1,6] on the global fits. For all the NP scenarios considered, except for Hyp. III, C9µ V = −C V , we observe that their corresponding 90 µ PullSM undergoes a ∼ 3 − 4σ variation from one end of the range of variation of hRK i[1.1,6] to the other. If we restrict the variation of hRK i[1.1,6] to only 1σ, one can see differences of ∼ 2σ between the two extremes, as expected from the linearity of PullSM on hRK i[1.1,6] seen in the plots. The flatness of the PullSM under the hypothesis III, C9µ V = −C V , can be easily under- 90 µ V = −C V , so stood. The theoretical prediction of hRK i[1.1,6] is insensitive to the value of C9µ 90 µ that it remains constant and equal to 1 to a very high accuracy. Therefore the difference between the theoretical and experimental values of hRK i[1.1,6] does not play any role in JHEP07(2019)096 the minimisation of the χ2 function. As a consequence, the b.f.p. is determined using the other observables of the fit, regardless of the experimental value for hRK i[1.1,6] (see figure 13), and the contribution of hRK i[1.1,6] cancels in the ∆χ2 = χ2SM − χ2min statistic. This explains the observed flat curve for the PullSM , up to small variations linked to the numerical minimisation of the χ2 function. The results in figure 4 show that, for most of the values of hRK i[1.1,6] scanned, it is not possible to fully disentangle all the NP scenarios, with the exception of Hyp III: C9µ V = −C9V0 µ . However, large values of hRK i[1.1,6] (around 0.90 or above) provide the potential to disentangle some of the NP LFUV scenarios. Many scenarios get their significances down to the range ∼ 3.8σ −4.8σ, apart from scenarios with right-handed currents like Hyps. V, VII, VIII. Indeed, if a new measurement of hRK i[1.1,6] is found in better agreement with its SM prediction, this favours right-handed currents for C9V0 µ cancelling the contribution for C9µ V, 0 + but there is still an important number of other tensions (i.e. RK ∗ , P5µ and B(Bs → φµ µ )) − that require NP contributions in order to be explained. Large values of hRK i[1.1,6] would help also to distinguish among NP scenarios featuring both LFUV and LFU NP, separating Hyps. X and XI from the others). We then study the combined influence of hRK i[1.1,6] and hQ5 i[1.1,6] . The value of hQ5 i[1.1,6] is varied as explained above and we repeat the analysis for three different values of hRK i[1.1,6] : its current experimental value and the ends of its 1σ range. Figures 5 and 6 show how the PullSM varies with different experimental values of hQ5 i[1.1,6] and hRK i[1.1,6] . In appendix A figures 14, 15 and 16 show how b.f.p.s are affected, with a rather large variety of behaviours. We observe that most of the hypotheses see their PullSM increase for larger values of hQ5 i[1.1,6] . One can separate the discussion according to the value of hRK i[1.1,6] • hRK i[1.1,6] ' 0.8: the pulls are larger than in the current case. The rather low value of RK disfavours in general scenarios with right-handed currents with respect to scenarios involving only SM vector operators. A large value of hQ5 i[1.1,6] (close to 1) favours C9µNP , whereas a low value supports NP in both C NP and C NP either from 9µ 10µ LFUV NP only (Hyps. II, IV) or from a combination of LFU and LFUV NP (Hyp. XI). • hRK i[1.1,6] ' 0.85: the new determination of hRK i[1.1,6] is then nominally close to its current experimental value but with smaller errors. Therefore, the values for the – 11 – Global Fits RK [1.1,6] = 0.789 (−1σ) Global Fits RK [1.1,6] = 0.789 (−1σ) 8 8 7 7 PullSM (σ) PullSM (σ) 6 6 5 NP  NP NP I: C9µ IV: C9µ , C10µ NP NP  NP  II: C9µ = −C10µ V: C9µ , C9µ NP  NP  4 III: C9µ = −C9µ 5 VI: C9µ , C10µ JHEP07(2019)096  NP NP  VII: C9µ = −C9µ, C10µ = C10µ  NP  VIII: C9µ , C9µ = −C10µ 3 4 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] Global Fits RK [1.1,6] = 0.846 (0σ) Global Fits RK [1.1,6] = 0.846 (0σ) 8 8 7 7 PullSM (σ) PullSM (σ) 6 6 5 NP  NP NP I: C9µ 5 IV: C9µ , C10µ NP NP  NP  II: C9µ = −C10µ V: C9µ , C9µ 4 NP III: C9µ = −C9µ  NP  VI: C9µ , C10µ  NP NP  4 VII: C9µ = −C9µ, C10µ = C10µ  NP  VIII: C9µ , C9µ = −C10µ 3 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] Global Fits RK [1.1,6] = 0.908 (+1σ) Global Fits RK [1.1,6] = 0.908 (+1σ) 8 8 7 7 PullSM (σ) PullSM (σ) 6 6 5 5  NP NP NP I: C9µ IV: C9µ , C10µ  NP  4 NP II: C9µ NP = −C10µ 4 V: C9µ , C9µ  NP  NP III: C9µ = −C9µ VI: C9µ , C10µ  NP NP  VII: C9µ = −C9µ, C10µ = C10µ 3  NP  3 VIII: C9µ , C9µ = −C10µ −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] Figure 5. Global fit: impact of Q5 [1.1,6] on the PullSM of the NP scenarios under consideration for different values of RK [1.1,6] . – 12 – Global Fits RK [1.1,6] = 0.789 (−1σ) Global Fits RK [1.1,6] = 0.846 (0σ) 8 8 7 7 PullSM (σ) PullSM (σ) 6 6  V   V  V IX: C9µ = −C10µ V , C9U = C10 U IX: C9µ = −C10µ , C9U = C10 U  V U  V U X: C9µ, C9 5 X: C9µ, C9  V   V V  XI: C9µ = −C10µ V , C9U XI: C9µ = −C10µ , C9U  V   V  JHEP07(2019)096 V U V XII: C9µ = −C10µ, C10 U XII: C9µ = −C10µ , C10 5  V U  V U XIII: C9µ , C10 XIII: C9µ, C10  V U 4  V U XIV: C9µ , C10 XIV: C9µ , C10 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] Global Fits RK [1.1,6] = 0.908 (+1σ) 8 7 6 PullSM (σ) 5  V V  IX: C9µ = −C10µ , C9U = C10 U  V U X: C9µ , C9 4  V V  XI: C9µ = −C10µ , C9U  V V U  XII: C9µ = −C10µ , C10  V U 3 XIII: C9µ , C10  V U XIV: C9µ , C10 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Figure 6. Global fit: impact of Q5 [1.1,6] on the PullSM of the NP scenarios under consideration for different values of RK [1.1,6] . b.f.p.s are numerically similar to the b.f.p.s reported in refs. [3, 4]. On one hand, low and large values of Q5 [1.1,6] provide both a rather clear separation between Hyps. II, IX, X, XI and Hyps. I, III, XIII, XIV. On the other hand, it does not help to separate a set of hypotheses with LFUV NP only (IV to VIII). • RK [1.1,6]  0.9: the pulls are lower with respect to their present values. Large values of Q5 [1.1,6] allow one to disfavour Hyp II, and low and large values of Q5 [1.1,6] provide still both a rather clear separation among hypotheses combining LFU and LFUV NP (Hyps. IX, X, XI on one hand and Hyps. XIII, XIV on the other). However, it does not help to separate a set of hypotheses with LFUV NP only (IV to VIII). – 13 – LFUV Fits 7 LFUV Fits 7 NP I: C9µ NP NP II: C9µ = −C10µ 6 NP III: C9µ = −C9µ 6 5 5 PullSM (σ) PullSM (σ) 4 4  NP NP 3 IV: C9µ , C10µ  NP  V: C9µ , C9µ 3  NP  JHEP07(2019)096 VI: C9µ , C10µ  NP  2 NP VII: C9µ = −C9µ, C10µ = C10µ  NP  VIII: C9µ , C9µ = −C10µ 2 0.75 0.80 0.85 0.90 0.95 0.75 0.80 0.85 0.90 0.95 RK [1.1,6] RK [1.1,6] LFUV Fits 6 5 PullSM (σ) 4  V V  IX: C9µ = −C10µ , C9U = C10 U 3  V U X: C9µ, C9  V V  XI: C9µ = −C10µ , C9U  V  2 XII: C9µ = −C10µ V U , C10  V U XIII: C9µ, C10  V U XIV: C9µ , C10 1 0.75 0.80 0.85 0.90 0.95 RK [1.1,6] Figure 7. LFUV fit: impact of the central value of RK [1.1,6] on the PullSM of the NP scenarios under consideration. In summary, if the update of RK [1.1,6] is larger than its current value, it can help distinguishing among various NP hypotheses, with a preference for hypotheses involving right-handed currents in C9V µ . A value of RK [1.1,6] similar or smaller than the current value would clearly disfavour the hypothesis C9µ V = −C V , but many other NP hypotheses 9 µ (with LFUV only or with a combination of LFU and LFUV) cannot be separated. The observable Q5 [1.1,6] is an excellent candidate to separate among some of these possibilities. Depending on the situation, low and/or large values of this observables provide a good separation in terms of pulls. 3.2 LFUV fits It is also interesting to address the impact of the observables RK [1.1,6] and Q5 [1.1,6] on the LFUV fits, which currently include all the RK and RK ∗ measurements, the measurements – 14 – of Qi (i = 4, 5) by the Belle collaboration, all the b → sγ observables available, as well as B(B → Xs µ+ µ− ) and B(Bs → µ+ µ− ) [39–41]. We follow the same guidelines as for the global fits. We show in figures 7 the impact of an update of hRK i[1.1,6] on the pull with respect to the Standard Model, whereas the role played by hQ5 i[1.1,6] is indicated by figures 8 and 9. One can also plot the variation of the best-fit point as a function of hRK i[1.1,6] and hQ5 i[1.1,6] , leading to plots similar to the case of the global fit, but the information gained this way is not very illuminating and we refrain from showing them. Most of the features observed in the global fit are also observed in the LFUV fit, although with lower pulls. For different values of hRK i[1.1,6] , the separation among hy- JHEP07(2019)096 potheses can be improved if one measures either low or large values of hQ5 i[1.1,6] . One can in particular notice that the separation between the LFUV NP hypotheses (IV to VIII) seems easier for low values of hQ5 i[1.1,6] compared to the fit “All”. This means that the observables in the LFUV fit are more sensitive to details of this scenario, but this gets compensated by other observables in the fit “All” so that the sensitivity is reduced in the more complete fit. LFUV fits lead us to draw similar conclusions to the ones extracted from the global fits. However they exhibit a stronger clustering of the pulls, especially if only hRK i[1.1,6] is used to discriminate them. Moreover, if we take a given hypothesis with LFUV-NP only, and consider hypotheses obtained by adding further LFU-NP contributions, we ob- tain very similar pulls. Therefore, the only way to distinguish among different LFU-NP hypotheses consists in performing the global fits and having access to both muonic and electronic branching ratios. The addition of hQ5 i[1.1,6] allows for strategies that enable the discrimination of various scenarios in a similar way to the case of the global fits. 4 Pulls of individual observables Discerning among different NP hypotheses using the significance of the Pull SM proved difficult using the current set of data. In section 3 we have discussed how new data could possibly improve the situation and allow us to disentangle the different scenarios. We discuss here a complementary approach based on analysing the current picture for each hypothesis using the correlated pull of individual observables. We consider the individual pull of an observable, which can be written for the ith observable as: q pullobs i = χ2min − χ2min w/o obs i , (4.1) where χ2min and χ2min w/o obs i are the minimal values of the χ2 with and without the ith observable X χ2 = (Oth − Oexp )i Vij−1 (Oth − Oexp )j , ij X χ2w/o obs i = (Oth − Oexp )i Vij−1 (Oth − Oexp )j , (4.2) i6=j – 15 – LFUV Fits RK [1.1,6] = 0.789 (−1σ) LFUV Fits RK [1.1,6] = 0.789 (−1σ) 8 8 7 7 6 6 5 PullSM (σ) PullSM (σ) 4 5 3 NP I: C9µ  NP NP 4 IV: C9µ , C10µ NP NP II: C9µ = −C10µ  NP  2 NP V: C9µ , C9 µ III: C9µ = −C9µ  NP  JHEP07(2019)096 3 VI: C9µ , C10µ  NP  1 VII: C9µ NP = −C9µ, C10µ = C10µ  NP  VIII: C9µ , C9µ = −C10µ 0 2 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] LFUV Fits RK [1.1,6] = 0.846 (0σ) LFUV Fits RK [1.1,6] = 0.846 (0σ) 8 8 7 7 6 6 5 PullSM (σ) PullSM (σ) 5 4 4 3 NP I: C9µ  NP NP IV: C9µ , C10µ NP II: C9µ NP = −C10µ 3  NP  2 NP V: C9µ , C9 µ III: C9µ = −C9µ  NP  VI: C9µ , C10µ 1 2  NP NP  VII: C9µ = −C9µ, C10µ = C10µ  NP  VIII: C9µ , C9µ = −C10µ 0 1 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] LFUV Fits RK [1.1,6] = 0.908 (+1σ) LFUV Fits RK [1.1,6] = 0.908 (+1σ)  NP NP 8 8 IV: C9µ , C10µ  NP  V: C9µ , C9µ  NP  7 7 VI: C9µ , C10µ  NP NP  VII: C9µ = −C9µ, C10µ = C10µ  NP  6 6 VIII: C9µ , C9µ = −C10µ 5 PullSM (σ) PullSM (σ) 5 4 4 3 NP I: C9µ 3 NP NP II: C9µ = −C10µ 2 NP III: C9µ = −C9µ 2 1 1 0 0 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] Figure 8. LFUV fit: impact of Q5 [1.1,6] for different values of RK [1.1,6] on the PullSM of the NP hypotheses under consideration. – 16 – LFUV Fits RK [1.1,6] = 0.789 (−1σ) 8 7 6 PullSM (σ) 5  V V  IX: C9µ = −C10µ , C9U = C10 U  V U X: C9µ, C9 4  V V  XI: C9µ = −C10µ , C9U  V  JHEP07(2019)096 V U XII: C9µ = −C10µ, C10  V U 3 XIII: C9µ , C10  V U XIV: C9µ , C10 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] LFUV Fits RK [1.1,6] = 0.846 (0σ)  V  8 IX: C9µ = −C10µ V , C9U = C10 U  V U X: C9µ, C9  V V  XI: C9µ = −C10µ , C9U 7  V  V U XII: C9µ = −C10µ, C10  V U XIII: C9µ, C10 6  V U XIV: C9µ , C10 PullSM (σ) 5 4 3 2 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] LFUV Fits RK [1.1,6] = 0.908 (+1σ)  V  8 IX: C9µ = −C10µ V , C9U = C10 U  V U X: C9µ, C9  V  7 XI: C9µ = −C10µ V , C9U  V V U  XII: C9µ = −C10µ , C10  V U 6 XIII: C9µ, C10  V U XIV: C9µ , C10 PullSM (σ) 5 4 3 2 1 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Figure 9. LFUV fit: impact of Q5 [1.1,6] for different values of RK [1.1,6] on the PullSM for the NP hypotheses under consideration. – 17 – ○ □ ◇●▼ ■ ▲ ■ ◆ ● ■■□▲◆ ● ▼○ ◇ ■●◇ ▼ □◆ ■○ ▲ ■ ■■ ▲◆□●▼◇○ ▼◇ ○□ ●■ ▲ ■◆ ◆ ▼● ◇ ○□◆▲ ■ ■ ○ ◇●▼▲ ◆■■□ ▲ ■ ■●▲ □ ▼ ○ ◆ ◇ ■ ▼ ▲○◆ ◇● □ ■ ▼ □ ◆◇○▼ ■▲ ● ■ □ ◇○ ◆▼ ▲ ■ ● ■ ○ ◇ ◆○ □▼ ▲■ ● ■ □ ◆■ ◇▼▲●■ ○ □ □ ◆◇○■▼ ▲ ● ■ JHEP07(2019)096 ■ ◇▼○□ ◆ ● ■▲ ◇ 0 1 2 3 4 Figure 10. pullobs defined in eq. (4.1) for each selected observable within different NP hypotheses (from I to VIII) compared to the Standard Model (black squares). where Oth and Oexp are the theoretical prediction and the experimental value for the observable respectively, and Vij−1 corresponds to the covariance matrix element {i, j}. This definition, already used in refs. [42–45], differs from the definition often adopted e.g. in the context of electroweak precision observables [46]: this naive pull is defined as the difference between the experimental value and the theoretical value at the b.f.p., normalised by the uncertainty. As discussed in ref. [47], this definition should be refined. The pull of an observable can be considered as the assessment of the impact of the additional external constraint to the fit given by this observable, which must be assessed using a pull involving the central values and the uncertainties with and without the constraint given by the observable. It can be easily shown that the definition of ref. [47] is equivalent to our definition. Moreover, this definition follows the same approach of pulls as for the comparison of hypotheses and it can be expected to follow a normal law with zero mean and unit width. We stress that it includes correlations, so that it might have a different value from the naive pull in the presence of large correlations among observables (from experimental or theoretical nature). In the following, we will compare the pulls, as defined in eq. (4.1), under the SM hypothesis and under various NP hypotheses for all observables. However, we will focus mainly on the observables yielding pullsobs larger than 1.5. These pullsobs with respect to the SM are expected to be reduced under the various NP hypotheses, but some might remain (while the others disappear), providing interesting insights into the role of the various observables within a given NP hypothesis. The observables under discussion are visualized in figures 10 and 11. Black squares represent the pullobs of the observable i within the SM computed following eq. (4.1), while coloured and empty shapes represent the pullobs of the same observable under different NP hypotheses. – 18 – ● ○ ▼ ■ ▲ ■ ◆ ■◆●▼▲■○ ■ ●◆ ▼ ▲ ○■ ● ■ ◆ ▲ ▼●■ ○ ○ ◆ ■▲▼ ● ■ ■ ■ ○● ▼▲ ◆ ■ ○ ■ ▼● ▲◆ ■ ◆ ■ ◆▲■▼● ○ ○▲▼ ■●◆ ■ ▲ ● ▼○ ▲ ◆■ ■ ● ○▼ ▲ ■◆ ■ ▼ ●○ ▼ ▲ ■◆ ■ ▼■■ ◆▲ ●○ ○ ◆ ▲ ○●▼ ■ ■ JHEP07(2019)096 ■ ○ ◆ ■ ▲ ▼● 0 1 2 3 4 Figure 11. pullobs defined in eq. (4.1) for each selected observable within different NP hypotheses (from IX to XIV) compared to the Standard Model (black squares). 4.1 Observables with a large pullobs i within the Standard Model Looking at figures 10 and 11, the first observable considered is hRK i[1.1,6] : it is very sensitive to LFUV NP, it has a large pullobs with respect to the SM. It has also a small pullobs in all the NP hypotheses considered except for Hyps. III and XI. Remarkably, it has almost no correlation with any other observable in the fit4 [48] and, as a consequence, the fit must satisfy both hRK i[1.1,6] and the rest of the observables in parallel, leading to a potential tension unless the solution for both sectors of the fit is similar. Another observable showing a tension with the SM is hP2 i[0.1,0.98] . Remarkably, none of the hypotheses considered is able to reduce the tension of this observable (see figures 10 and 11) so it remains poorly explained either within the SM or under the NP hypotheses considered. A similar problem is seen in hP50 i[0.1,0.98] even though for this observable the pullobs within the SM is smaller than for hP2 i[0.1,0.98] . Including a NP contribution to C7 (which affects mainly the first bin) does not improve these two observables as they require contributions of opposite signs. Therefore, the small tensions observed in the first bin for these two observables will require more data to be understood. The polarization fraction hFL i[2.5,4] follows a similar trend as hP50 i[0.1,0.98] . Its current measurement shows little deviation with the SM prediction [6] in any of its bins due to the large theoretical errors, but its pullobs is around ∼ 2 under the SM hypothesis due to the correlations with the rest of the observables involved in the fit. Moreover, for all NP hypotheses considered its value for the pullobs is larger than in the SM case. In the case of the observables hP50 i[4,6] and hP50 i[6,8] the tension between the SM predic- tion and the experimental value is somewhat reduced if one introduces NP contributions. The less efficient hypothesis is Hyp. II, C9µ NP = −C NP , and one of the best possibilities 10µ 4 There is little theoretical correlation between RK and the rest of the observables as the hadronic form factors, which are the main source of correlated uncertainty among observables, cancel in RK . The experimental correlation between RK and the B + → K + `` branching ratio is not public and therefore assumed to be zero here. We checked the lack of correlation between hRK i[1.1,6] and the rest of the observables at the level of our covariance matrix. – 19 – 〈�� 〉 ��� 2.0 〈� �� �〉[���] ��� (����) 〈�� 〉(� � μ =-� ��� �� μ ) ��� 〈�� 〉(� ��� μ =-���� μ � �� � =��� =-����) � 〈�� 〉(� ��� μ =-���� μ � �� � =��� =-����) � 1.5 〈�� 〉(� ��� μ =-���� μ � �� � =��� =-����) � 〈�� 〉(� � μ) ��� 〈R K 〉 1.0 JHEP07(2019)096 0.5 -0.8 -0.6 -0.4 -0.2 0.0 0.0 〈P5'〉[4,6] 0.2 0.4 Figure 12. RK [1.1,6] versus P5 [4,6] in three different scenarios: C9µ V V (grey), C9µ V = −C10µ (red), V V U U and C9µ = −C10µ , C9 = C10 (three different values of the LFU contributions, blue, light blue and purple). In each case, the band is obtained by varying the LFU contribution within 1σ. The current experimental values from LHCb are also indicated (green vertical and orange horizontal bands). is Hyp. I, i.e. a single contribution C9µ NP . As shown in figure 12, the latter hypothesis (grey band) better accommodates at the same time both the deviations in P5 [4,6] and NP = −C NP , although the updated experimental value RK [1.1,6] w.r.t. the solution with C9µ 10µ of RK [1.1,6] makes the compatibility only partial, whereas the former hypothesis (red band) is unable to obtain the same agreement. If one favours this scenario C9µ NP = −C NP , 10µ a compatibility level similar to that of C9µ NP can be obtained by adding a LFU contri- bution C9U = C10 U [14] (corresponding to Hyp. IX), since R  V K [1.1,6] = 1 + 0.49C9µ while  U V U 2 V 2 U V P5 [4,6] = −0.81 − 0.27C9 − 0.14C9µ + 0.07(C9 ) + 0.05(C9µ ) + 0.11C9 C9µ [14]. This im- plies that a NP scenario with C9µ V = −C V would then need a LFU universal contribution 10µ C9U = C10 U to be viable. This is illustrated in figure 12 where Hyp. IX is displayed varying the LFU contribution within its 1σ confidence interval according to the global fit (blue, light blue and purple bands). One of the tensions of the fits described in section 2.3 belongs to the low-recoil region. This corresponds mainly to two observables here: P5 [15,19] and B(B 0 → K ∗0 µ+ µ− )[15,19] . For P5 [15,19] , including NP contributions leads to an increase of pullobs in all cases, see figures 10 and 11. The reduction of this tension would require NP contributions of opposite sign to the ones favoured by the rest of the observables of the fit. On the other hand, the branching ratio B(B 0 → K ∗0 µ+ µ− )[15,19] has its pullobs significantly reduced under all NP hypotheses. This shows that the correlations between B(B 0 → K ∗0 µ+ µ− )[15,19] and the rest of the observables allow for the reduction of tensions under the different favoured NP hypotheses. – 20 – Deviations with respect to the SM are also found for the branching ratio of the Bs → φµ+ µ− channel. All NP hypotheses yield pullsobs smaller than their SM counterparts, specially for the low-recoil bin and parallel to its B → K ∗ counterpart. The preferred scenarios considering the three bins of the branching ratio in the Bs → φµ+ µ− channel from figures 10 and 11 are those with C9µ V = −C V . 10µ Finally, the tension for hRK ∗ i[1.1,6] is also alleviated by introducing a NP contribution. Whereas the large bin hRK ∗ i[1.1,6] gets a significant reduction in its pullobs under all NP hypotheses, the first bin [0.1, 0.98] is just slightly improved when compared to the SM. This could point towards an inconsistency between the first bin and the other bins of RK ∗ if one takes into account the issues discussed above regarding the first bin of P2 and JHEP07(2019)096 P50 . We should comment here on ref. [7], where several NP scenarios are discussed in order to explain the low value of hRK ∗ i[0.1,0.98] , including contributions to both b → sµµ and b → see channels as well as considering the presence of right-handed currents. It is possible to connect some of the scenarios from [7] with the description of NP in terms of LFU and LFUV contributions proposed in ref. [14]. Remarkably, the Wilson coefficients obtained in scenario S10 of [7] (following their notation) are in excellent agreement with the results obtained in [14] for the 4D fit with {C9µ V , C V , C U , C U }. 10µ 9 10 Hyp. II works better than Hyp. I in reducing the pullobs of the LFUV observables, as can be easily understood looking at the structure of the observables in terms of their Wilson coefficients provided in ref. [14]: hRK i[1,6] = 1 + 0.23C9µ V − 0.26C V while hR ∗ i 10µ K [1.1,6] = V V 1 + 0.16C9µ − 0.29C10µ . Then, the ratio (hRK i[1,6] − 1)/(hRK ∗ i[1.1,6] − 1) turns out to be 1.4 in Hyp. I and 1.1 in Hyp. II. 4.2 Observables with no pullobs i with respect to the Standard Model The kind of analysis undertaken in the previous subsection can also be conducted starting from observables that currently have very small pullsobs with respect to the SM predictions and see how they are affected once a fit under a NP hypothesis is performed. It is to expect their pullsobs will grow under NP hypotheses unless they are observables insensitive to the NP hypotheses considered. For this purpose we consider a group of observables with pulls between 0 and 0.2 within the SM. After a detailed analysis, we found very small pullsobs for all of them, except hB(B + → K + µ+ µ− )i[0.1,0.98] , under all the NP hypotheses considered, i.e. they are basically NP insensitive. In other words, these observables could be removed from the fit and the rest of the observables would remain unaffected. It is remarkable that again an observable taken in the first bin, hB(B + → K + µ+ µ− )i[0.1,0.98] , has a worse pullobs under all NP hypotheses than under the SM. This goes in the same direction as the issues observed in the first bin of observables of the B → K ∗ µ+ µ− channel such as P2 , P50 and RK ∗ , which points out the need of a better understanding of the low-q 2 region. 5 Sensitivity of observables to different Wilson coefficients A detailed scrutiny of the dependence in Wilson coefficients of b → s`` observables is useful in order to understand tensions in the global fit under various NP hypotheses. – 21 – The complexity of the calculation prevents such systematic study, but we will rely on an approximate parametrisation of the observables in order to identify the sensitivity of interesting observables to different Wilson coefficients/NP hypotheses. 5.1 Observables and parametrisations We use our current programs [6] to generate predictions according to a grid in the space of Wilson parameters and fit the proposed parametrisations to this set of pseudo-data. We collect the coefficients in the vector  NP NP NP NP X = C9µ , C10µ , C90 µ , C100 µ , C9e , C10e , C90 e , C100 e (5.1) JHEP07(2019)096 For simplicity, we consider only real NP contributions to the Wilson coefficients C9,90 ,10,100 for muon and for electron operators: we assume that we can neglect NP contributions to C7 , C70 (based on b → sγ observables in good agreement with the SM) as well as to scalar, pseudoscalar and tensor operators (based on the outcome of the general global fits performed on b → s`` observables which do not favour large contributions to these Wilson coefficients. We also neglect imaginary contributions for the Wilson coefficients (as there are no signals of NP sources of CP violation from the corresponding b → s`` observables). We consider only the set of LHCb observables present in the global fit of ref. [6] with the corresponding binning in q 2 , as well as the Belle observables Q4 and Q5 . NP (` = µ, e) in the range [−1, 1] for i = 10, 90 , 100 , We generate points along a grid for all Ci` but [−2, 0] for i = 9, computing both central values and theoretical uncertainties for all the observables of interest. We can divide the set of observables into three different groups: LFD observables, LFUV differences and LFUV ratios, with three different approximate quadratic parametrisations: • LFD observables: this category includes branching ratios and (optimized and aver- aged) angular observables [49, 50] governed by a lepton-specific b → s`` transition. We use a general quadratic parametrisation in the NP contributions to the Wilson coefficients (i, j = 9, 10, 90 , 100 ): X X NP NP NP O =α+ βi Ci` + γij Ci` Cj` (5.2) i i≥j • LFUV differences: Q1,2,4,5 and δS5,6s observables that measure differences in angular observables between muonic and electronic channels [51], (0 ) (0 ) , Qi = Piµ − Pie δSi = Siµ − Sie (5.3) In this case we use the following parametrisation (i, j = 9, 10, 90 , 100 ) X  NP NP Oµ − O e = α µ − α e + βiµ Ciµ − βie Cie (5.4) i X  NP NP NP NP + γijµ Ciµ Cjµ − γije Cie Cje i≥j – 22 – • LFUV ratios: these type of observables are defined as ratios of integrated branching fractions involving muonic (numerator) and electronic (denominator) final states and we use (i, j = 9, 10, 90 , 100 ): P NP + P NP NP Oµ αµ + i βiµ Ciµ i≥j γijµ Ciµ Cjµ = P NP + P NP NP (5.5) Oe αe + i βie Cie i≥j γije Cie Cje Not all observables can be parametrised in a fully satisfying way over this grid through such a quadratic approximation. The difference is at most of 0.3 for all observables (cen- tral value and uncertainty) over the whole range of scan. Pi (B → K ∗ µ+ µ− ) observables JHEP07(2019)096 turn out to be more difficult to fit with a purely quadratic parametrisation, with differ- ences between the approximate parametrisation and the exact values above 0.15 in some parts of the scan range (this could be expected due to their normalisation, leading to rapidly varying functions when the Wilson coefficients are changed). However, even this approximate parametrisation is able to catch a few interesting aspects of the sensitivity to Wilson coefficients. 5.2 Observables as conics Let us start the discussion with the b → sµµ observables, with a quadratic parametrisation that can be seen as a conic in the space of O = a + (X − X0 )T .M.(X − X0 ) (5.6) The symmetric matrix M can be diagonalised as M = RT .∆.R where R is an orthogonal (rotation) matrix, and ∆ is a diagonal matrix with eigenvalues ordered in (absolute) size. We can define the (rotated) Y -basis Y = R.(X − X0 ) O = a + RT .∆.R (5.7) The largest eigenvalues of ∆ correspond to directions in the Y -basis that are able to change the value of O in the most significant way. These directions correspond to linear combinations of Wilson coefficients to which the observable O is the most sensitive. We can consider the directions corresponding to the largest eigenvalues in absolute value (“dominant ones”) and we can determine whether some directions are similar for different observables. This will indicate that some observables share the same sensitivity to NP for their dominant directions. We perform this comparison considering only the dominant eigenvalue as well as the eigenvalues that are at most 80% of this eigenvalue, 5 and we consider “parallel” directions that correspond to pairs of vectors with an angle with a cosine of 0.9 or more (in absolute value). We perform this comparison for all the bins of the different observables. Due to the large number of directions to consider, we give our conclusions in a compact way, where each observable mentioned means “some of the bins of this observable”. We can observe several classes of observables exhibiting dominant parallel directions: 5 Let us add that the comparison of eigenvalues can be performed only for a given observable (possibly in different bins). The rescaling of an observable by an arbitrary factor will affect the absolute values of the eigenvalues but not their relative importance. – 23 – • observables which are expected to exhibit similar directions due to flavour symmetries B(B 0 → K ∗0 ), B(B + → K ∗+ ), B(Bs → φ) B(B 0 → K 0 ), B(B + → K + ) Pi (B 0 → K ∗0 ), Pi (Bs → φ) (i = 1, 4, 6, FL ) (5.8) • observables from the same process B(B 0 → K ∗0 ), Pi (B 0 → K ∗0 ) (i = 1, 3, 4, 5, 6, 8, FL ) JHEP07(2019)096 B(Bs → φ), Pi (Bs → φ) (i = 1, 4, 6, FL ) Pi (Bs → φ), Pj (Bs → φ) ({i, j} = {FL , 1}, {FL , 4}, {1, 4} Pi (B → K ), Pj (B 0 → K ∗0 ) 0 ∗0 ({i, j} = {FL , 1}, {FL , 2}, {FL , 4}, {FL , 8}, {1, 2}, {1, 3}, {1, 4}, {1, 5}, {1, 8}, {2, 4}, {2, 5}, {2, 6}, {2, 8}, {3, 4}, {3, 8}, {4, 8}, {6, 7}) (5.9) • observables with unexpected correlations B(B + → K + ), B(B 0 → K 0 ), P1 (B 0 → K ∗0 ), P1 (Bs → φ), P3 (B 0 → K ∗0 ) (5.10) If we now restrict the analysis to observables that are currently deviating by more than 2σ, we can easily consider a wider range of directions in each case (going down to eigenvalues that are 30% of the maximal ones), we can provide a slightly more precise description of the correlations • between observables as expected from flavour symmetries: B(B 0 → K ∗0 )[15, 19], B(B + → K ∗+ )[15, 19], B(Bs → φ)[15, 18.8] (5.11) • between different close bins of the same observable P50 (B 0 → K ∗0 )[4, 6], P50 (B 0 → K ∗0 )[6, 8] B(Bs → φ)[2, 5], B(Bs → φ)[5, 8] (5.12) • between observables that are not obviously correlated B(Bs → φ)[2, 5], B(Bs → φ)[5, 8], P50 (B 0 → K ∗0 )[4, 6], B(B 0 → K ∗0 )[15, 19], B(B + → K ∗+ )[15, 19] (5.13) – 24 – The same analysis carried out for LFUV quantities involves the same correlations, which is to be expected as they involve ratio or differences of the same observables considered in the muon and the electron cases. Another comment is in order. We have performed our analysis of the directions in the general case where NP is allowed in C9,10,90 ,100 and we have studied the dominant directions for each observable. Scenarios where some of these NP contributions are assumed to be zero correspond to projections of this analysis on specific hyperplanes. Obviously, if directions are similar in the general C9,10,90 ,100 space, they remain similar after projection on a given hyperplane, but the projection can also lead to additional pairs of observables with (projected) dominant directions that are parallel (although they are not parallel when JHEP07(2019)096 considered in the whole space). 5.3 Directions favoured by the anomalies Up to now, we have considered only the leading directions (corresponding to the main axes of the conics) for the quadratic approximation of the theoretical expression of the b → sµµ observables. This provides information on the sensitivity to specific directions and the fact that some of these directions are common to several observables. However, since we have not compared these predictions to the deviations observed currently, this does not provide information on the specific changes in the Wilson coefficients preferred by the global fit. One can perform a global analysis of the constraints [48], but one can also exploit the above parametrisation to gain some insights to understand better the outcome of the global fit. This can be done in the following way. As indicated above, we can write a b → sµµ observable in the Y -space as O = a + Y T .∆.Y (5.14) where the largest axis of the conic correspond to the 1st coordinate in the Y -space, with smaller and smaller axes correspond to the 2nd, 3rd. . . coordinates. Each observable tends to favour a shift in Y in order to reduce O − Oexp compared to OSM − Oexp . One can represent this problem as the fact that the conic O = Oexp does not pass through the SM point YSM = −R.X0 . This distance corresponds to the NP shift to perform in some of the Wilson coefficients. In principle, one could chose any NP shift in the Wilson coefficients to go from the SM point to the conic section O = Oexp , but it turns out to be interesting to think in terms of the main axes of the conic section and first to determine the distance between the SM point and the conic section along a given axis i. This is obtained by s P Oexp − a − k6=i YSM,k2 (i) δYj = ±δij (5.15) ∆i which corresponds to a NP in the Wilson coefficients δX (i) = X0 + RT .δY (i) . Let us emphasize that this equation does not necessarily have a solution, as it is not necessarily possible to reach the conic section from the SM in all directions (this is easily seen in the case of an ellipsoid, depending on the position of the SM point, inside, outside and close or outside and far away). The above equation illustrates two different aspects of the problem: on the one hand, a large shift in O is achieved along the 1st coordinate through – 25 – Observable Direction i δY (i) Xi 0 ∗0 P2 (B → K ) [0.1,0.98] 1 11.82 ± 5.89 (0.01, 0.05, 0.66, 0.75) 0 ∗0 P2 (B → K ) [0.1,0.98] 2 4.71 ± 3.62 (−0.47, −0.88, 0.02, 0.04) P2 (B 0 → K ∗0 ) [0.1,0.98] 4 7.93 ± 6.30 (0.88, −0.47, 0.02, −0.00) P50 (B 0 ∗0 → K ) [4.0,6.0] 1 1.83 ± 0.24 (0.94, 0.22, −0.26, 0.02) P50 (B 0 ∗0 → K ) [4.0,6.0] 2 1.96 ± 0.28 (−0.29, 0.49, −0.57, 0.59) P50 (B 0 ∗0 → K ) [4.0,6.0] 4 8.85 ± 1.82 (−0.05, −0.63, −0.74, −0.22) P50 (B 0 → K ∗0 ) [6.0,8.0] 1 1.95 ± 0.32 (0.87, 0.43, −0.22, −0.04) JHEP07(2019)096 P50 (B 0 ∗0 → K ) [6.0,8.0] 2 1.84 ± 0.37 (−0.31, 0.40, −0.60, 0.62) 0 ∗0 B(B → K ) [15.,19.] 2 0.79 ± 0.15 (0.14, −0.69, −0.14, 0.70) B(B + → K ∗+ ) [15.,19.] 2 1.18 ± 0.17 (0.14, −0.69, −0.14, 0.70) B(Bs → φ) [2.0,5.0] 2 2.34 ± 0.63 (−0.03, −0.70, 0.02, 0.71) B(Bs → φ) [15.,18.8] 2 0.65 ± 0.16 (−0.36, 0.60, 0.36, −0.62) Table 1. Distance δY (i) from the SM to the LFD conic O − Oexp following the main axis i of the conic section. The last column indicates the vector corresponding to the main axis i in the basis X = (C9µ , C10µ , C90 µ , C100 µ ). a smaller shift in Y (and thus in the Wilson coefficients) than along higher coordinates, but on the other hand, the distance between the SM point and the conic section might be easier to bridge along higher coordinates than along the first coordinate. The shifts δY (i) are thus interesting tools to identify the directions preferred by each observable that deviate significantly from the SM. If we perform this exercise for the observables deviating by more than 2σ and keeping only dominant directions (at least 30% of the largest eigenvalue), we obtain the results in table 1. We computed the distances δY (i) taking into account both experimental and theoretical uncertainties, as indicated in the corresponding table. We also indicate the corresponding directions in terms of X-coordinates. We show only directions which can accommodate the central value experimentally measured using the parametrisation that we used. We see that • P2 is affected by large uncertainties that make it difficult to reach definite conclusions • the 4 branching ratios prefer a single direction, which is essentially along ' C10µ −C100 µ • the two bins in P50 have similar behaviours, with two directions favoured (C9µ with a small C10µ component, or along ' C10µ − C90 µ + C100 µ ) The opposite signs for C100 between the branching ratios and the two bins in P50 make it difficult to use this specific parameter to improve the agreement of all observables with experiment. A better agreement between theory and data for these observables can thus be reached by performing shifts in C9µ and C10µ . Naturally, this very qualitative argument – 26 – does not take into account the remaining observables, which are in good agreement with the SM and constrain also the size of the NP shifts in C9µ and C10µ . However it is interesting to see that this rough analysis in terms of favoured directions supports the outcome of the global fit. A similar analysis can be performed in the case of RK and RK ∗ , but almost all direc- tions are equally favoured in both muon and electron directions, without shedding further light on the preferences for the global fit. 6 Conclusions JHEP07(2019)096 Over the last few years, the rare b → s`` decays have proved particularly interesting, with a large set of deviations from SM expectations. If a first explanation was provided by a large NP contribution to the Wilson coefficient C9µ , the current set of data allows for several different NP scenarios, which feature large NP contributions violating lepton-flavour universality (affecting b → sµµ but not b → see) in connection with the measurement of the LFUV ratios RK and RK ∗ . This does not exhaust the possibilities, and for instance, it is also possible to add NP contributions satisfying lepton-flavour universality (affecting b → sµµ and b → see) in the same way. With this wealth of possibilities, it is important to understand how current and forthcoming data may favour one scenario over others, so that one can pin down the best NP explanation for the whole set of anomalies observed. We have first considered some of the current inner tensions of the fits. We highlighted the situation of RK ∗ in the lowest bin, the tension between low and large recoils for Bs → φµ+ µ− , and a similar issue with B → K ∗ µ+ µ− angular observables. We then discussed the impact of forthcoming measurements of hRK i[1.1,6] and hQ5 i[1.1,6] to disentangle NP hypotheses. We considered various central values for these two measurements and we assumed some reduction in the experimental uncertainties in order to study how pulls w.r.t. SM and best-fit points would evolve. hRK i[1.1,6] alone proves to have only a limited ability to separate the various NP hypotheses: C9µ V = −C V is the only hypothesis strongly 90 µ affected. On the other hand, the combination of hRK i[1.1,6] and hQ5 i[1.1,6] proves much more efficient to separate various favoured hypotheses, either with only LFUV NP contributions or with both LFUV and LFU contributions. We then considered two different approaches to understand the current structure of the global fit. Firstly, we discussed the pulls associated with individual observables. Their value can be compared under the SM hypothesis and under various NP hypotheses, indicating whether they are easier to accommodate once NP contributions are added. As could be expected, no single hypothesis is clearly preferred: under each hypothesis, some of the observables get larger pulls and others smaller ones. Secondly, we looked at approximate quadratic parametrisations of the observables with large pulls in the SM. We could then identify directions in the space of Wilson coefficients to which the observables showed a large sensitivity to any shift. Interestingly, the same directions are often favoured by a large set of observables. We also analysed the shifts preferred by the observables showing large tensions with the SM using the same approach of directions. – 27 – The overall conclusion of our study is that the current measurements of b → s`` observables are not enough to disentangle various NP hypotheses: the tensions with the SM do not point towards an unambiguous NP scenario, and several compete at the same level, with only LFUV-NP contributions or combining LFUV and LFU contributions. An update of hRK i[1.1,6] could prove interesting in confirming (or not) the existence of LFUV NP contributions, but our study shows that in most cases, this update will not help in lifting the degeneracy among the various NP hypotheses. On the other hand, we have shown that the determination of hQ5 i[1.1,6] might bring interesting information that would be complementary to the existing observables and that could help significantly in disentangling the various scenarios (depending on the actual measured value). The use of extended data JHEP07(2019)096 sets, the analyses in different experimental environments and the inclusion of additional observables will all prove particularly important in the coming months in order to cross- check and constrain the dynamics at work in b → s`` transitions. This will prove essential to identify models of New Physics that could ultimately resolve the anomalies currently observed in the b-quark sector. Acknowledgments This work received financial support from the grants CICYT-FEDER-FPA 2014-55613-P and FPA 2017-86989-P [JM, SDG, BC, MA], and from Centro de Excelencia Severo Ochoa SEV-2012-0234 [BC]; from the EU Horizon 2020 program from the grants No. 690575, No. 674896 and No. 692194 [SDG]. J.M. acknowledges financial support from the Catalan ICREA Academia Program. The work of P.M. is supported by the Beatriu de Pinós postdoctoral program of the Government of Catalonia’s Secretariat for Universities and Research of the Ministry of Economy and Knowledge of Spain. A Dependence of best-fit points on hRK i[1.1,6] and hQ5 i[1.1,6] In this appendix, we show how the best-fit point of each NP hypothesis considered varies assuming particular values of hRK i[1.1,6] and hQ5 i[1.1,6] . For the two-dimensional hypothesis, we show in solid and dashed lines of the same colour the variations of the two parameters describing the NP contributions of the hypotheses. – 28 – Global Fits Global Fits −0.30 0.75 0.50 −0.45 0.25 −0.60 0.00 NP IV: C9µ BFP BFP NP I: C9µ NP IV: C10µ −0.75 NP II: C9µ NP = −C10µ −0.25 NP V: C9µ NP III: C9µ = −C9µ V: C9µ −0.50 −0.90 NP VI: C9µ −0.75 VI: C10µ JHEP07(2019)096 −1.05 −1.00 −1.20 −1.25 0.75 0.80 0.85 0.90 0.95 0.75 0.80 0.85 0.90 0.95 RK [1.1,6] RK [1.1,6] Global Fits Global Fits 0.50 V V IX: C9µ = −C10µ 0.50 0.25 IX: C9U = C10 U V X: C9µ 0.25 0.00 X: C9U V V XI: C9µ = −C10µ 0.00 −0.25 XI: C9U BFP BFP −0.25 NP VII: C9µ = −C9µ −0.50 NP VII: C10µ = C10µ −0.50 NP VIII: C9µ −0.75 VIII: C9µ = −C10µ −0.75 −1.00 −1.00 −1.25 −1.25 −1.50 0.75 0.80 0.85 0.90 0.95 0.75 0.80 0.85 0.90 0.95 RK [1.1,6] RK [1.1,6] Global Fits V V 0.2 XII: C9µ = −C10µ U XII: C10 V XIII: C9µ 0.0 U XIII: C10 V XIV: C9µ −0.2 U XIV: C10  BFP −0.4 −0.6 −0.8 −1.0 −1.2 0.75 0.80 0.85 0.90 0.95 RK [1.1,6] Figure 13. Global fit: impact of the central value of RK [1.1,6] on the b.f.p.s of the NP scenarios under consideration. – 29 – Global Fits RK [1.1,6] = 0.789 (−1σ) Global Fits RK [1.1,6] = 0.789 (−1σ) −0.45 0.50 0.25 −0.60 0.00 −0.75 −0.25 NP IV: C9µ −0.90 BFP BFP −0.50 NP IV: C10µ NP V: C9µ −1.05 −0.75 V: C9µ NP VI: C9µ −1.20 I: NP C9µ −1.00 VI: C10µ JHEP07(2019)096 NP NP II: C9µ = −C10µ NP −1.25 −1.35 III: C9µ = −C9µ −1.50 −1.50 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] Global Fits RK [1.1,6] = 0.846 (0σ) Global Fits RK [1.1,6] = 0.846 (0σ) −0.30 0.6 −0.45 0.3 −0.60 0.0 −0.75 NP IV: C9µ −0.3 BFP BFP NP −0.90 IV: C10µ NP V: C9µ −0.6 −1.05 V: C9µ NP VI: C9µ I: NP C9µ −0.9 VI: C10µ −1.20 NP NP II: C9µ = −C10µ NP −1.2 III: C9µ = −C9µ −1.35 −1.5 −1.50 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] Global Fits RK [1.1,6] = 0.908 (+1σ) Global Fits RK [1.1,6] = 0.908 (+1σ) NP IV: C9µ −0.30 0.9 NP IV: C10µ NP 0.6 V: C9µ −0.45 V: C9µ NP 0.3 VI: C9µ −0.60 VI: C10µ 0.0 −0.75 BFP BFP −0.3 −0.90 −0.6 −1.05 NP I: C9µ −0.9 −1.20 NP II: C9µ NP = −C10µ NP III: C9µ = −C9µ −1.2 −1.35 −1.5 −1.50 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] Figure 14. Global fit: impact of Q5 [1.1,6] on the b.f.p.s of the NP hypotheses under consideration for different values of RK [1.1,6] . – 30 – Global Fits RK [1.1,6] = 0.789 (−1σ) Global Fits RK [1.1,6] = 0.789 (−1σ) 0.50 0.6 0.25 0.3 0.00 0.0 −0.25 VII: NP C9µ = −C9µ −0.3 NP VII: C10µ = C10µ BFP BFP −0.50 NP VIII: C9µ −0.6 VIII: C9µ = −C10µ −0.75 −0.9 V V IX: C9µ = −C10µ −1.00 −1.2 IX: C9U = C10 U JHEP07(2019)096 V X: C9µ −1.25 −1.5 X: C9U V V XI: C9µ = −C10µ −1.50 −1.8 XI: C9U −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] Global Fits RK [1.1,6] = 0.846 (0σ) Global Fits RK [1.1,6] = 0.846 (0σ) 0.50 0.25 0.25 0.00 0.00 −0.25 NP −0.25 VII: C9µ = −C9µ NP VII: C10µ = C10µ −0.50 −0.50 BFP BFP NP VIII: C9µ VIII: C9µ = −C10µ −0.75 −0.75 V V −1.00 IX: C9µ = −C10µ −1.00 IX: C9U = C10 U −1.25 V X: C9µ −1.25 X: C9U −1.50 V XI: C9µ V = −C10µ −1.50 XI: C9U −1.75 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] Global Fits RK [1.1,6] = 0.908 (+1σ) Global Fits RK [1.1,6] = 0.908 (+1σ) 0.2 0.25 0.0 0.00 −0.2 −0.25 NP VII: C9µ = −C9µ −0.4 NP VII: C10µ = C10µ −0.50 BFP BFP NP VIII: C9µ −0.6 VIII: C9µ = −C10µ −0.75 −0.8 V V IX: C9µ = −C10µ −1.00 −1.0 IX: C9U = C10 U V X: C9µ −1.25 −1.2 X: C9U V V XI: C9µ = −C10µ −1.50 −1.4 XI: C9U −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Q5[1.1,6] Figure 15. Global fit: impact of Q5 [1.1,6] on the b.f.p.s of the NP hypotheses under consideration for different values of RK [1.1,6] . – 31 – Global Fits RK [1.1,6] = 0.789 (−1σ) 0.25 0.00 −0.25 −0.50 BFP −0.75 V V XII: C9µ = −C10µ −1.00 U XII: C10 JHEP07(2019)096 V XIII: C9µ −1.25 U XIII: C10 V XIV: C9µ −1.50 U XIV: C10  −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Global Fits RK [1.1,6] = 0.846 (0σ) 0.2 0.0 −0.2 −0.4 BFP −0.6 −0.8 V V XII: C9µ = −C10µ −1.0 U XII: C10 V XIII: C9µ −1.2 U XIII: C10 V −1.4 XIV: C9µ U XIV: C10  −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Global Fits RK [1.1,6] = 0.908 (+1σ) 0.2 0.0 −0.2 −0.4 BFP −0.6 −0.8 V V XII: C9µ = −C10µ U XII: C10 −1.0 V XIII: C9µ U −1.2 XIII: C10 V XIV: C9µ U −1.4 XIV: C10  −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 Q5[1.1,6] Figure 16. Global fit: impact of Q5 [1.1,6] on the b.f.p.s of the NP hypotheses under consideration for different values of RK [1.1,6] . – 32 – 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. References [1] LHCb collaboration, Search for lepton-universality violation in B + → K + `+ `− decays, Phys. Rev. Lett. 122 (2019) 191801 [arXiv:1903.09252] [INSPIRE]. [2] Belle collaboration, Test of lepton flavor universality in B → K ∗ `+ `− decays at Belle, JHEP07(2019)096 arXiv:1904.02440 [INSPIRE]. [3] M. Algueró et al., Emerging patterns of new physics with and without lepton flavour universal contributions, arXiv:1903.09578 [INSPIRE]. [4] B. Capdevila et al., Patterns of new physics in b → s`+ `− transitions in the light of recent data, JHEP 01 (2018) 093 [arXiv:1704.05340] [INSPIRE]. [5] S. Descotes-Genon, J. Matias and J. Virto, Understanding the B → K ∗ µ+ µ− anomaly, Phys. Rev. D 88 (2013) 074002 [arXiv:1307.5683] [INSPIRE]. [6] S. Descotes-Genon, L. Hofer, J. Matias and J. Virto, Global analysis of b → s`` anomalies, JHEP 06 (2016) 092 [arXiv:1510.04239] [INSPIRE]. [7] J. Kumar and D. London, New physics in b → se+ e− ?, Phys. Rev. D 99 (2019) 073008 [arXiv:1901.04516] [INSPIRE]. [8] M. Ciuchini et al., New Physics in b → s`+ `− confronts new data on lepton universality, arXiv:1903.09632 [INSPIRE]. [9] J. Aebischer et al., B-decay discrepancies after Moriond 2019, arXiv:1903.10434 [INSPIRE]. [10] A.K. Alok, A. Dighe, S. Gangal and D. Kumar, Continuing search for new physics in b → sµµ decays: two operators at a time, JHEP 06 (2019) 089 [arXiv:1903.09617] [INSPIRE]. [11] A. Arbey et al., Update on the b → s anomalies, arXiv:1904.08399 [INSPIRE]. [12] M. Ciuchini et al., On flavourful easter eggs for new physics hunger and lepton flavour universality violation, Eur. Phys. J. C 77 (2017) 688 [arXiv:1704.05447] [INSPIRE]. [13] T. Hurth, F. Mahmoudi, D. Martinez Santos and S. Neshatpour, Lepton nonuniversality in exclusive b→s`` decays, Phys. Rev. D 96 (2017) 095034 [arXiv:1705.06274] [INSPIRE]. [14] M. Algueró et al., Are we overlooking lepton flavour universal new physics in b → s``?, Phys. Rev. D 99 (2019) 075017 [arXiv:1809.08447] [INSPIRE]. [15] B. Capdevila, S. Descotes-Genon, L. Hofer and J. Matias, Hadronic uncertainties in B → K ∗ µ+ µ− : a state-of-the-art analysis, JHEP 04 (2017) 016 [arXiv:1701.08672] [INSPIRE]. [16] Belle collaboration, Lepton-flavor-dependent angular analysis of B → K ∗ `+ `− , Phys. Rev. Lett. 118 (2017) 111801 [arXiv:1612.05014] [INSPIRE]. [17] A. Datta, J. Kumar, J. Liao and D. Marfatia, New light mediators for the RK and RK ∗ puzzles, Phys. Rev. D 97 (2018) 115038 [arXiv:1705.08423] [INSPIRE]. [18] W. Altmannshofer et al., Light resonances and the low-q2 bin of RK ∗ , JHEP 03 (2018) 188 [arXiv:1711.07494] [INSPIRE]. – 33 – [19] M. Misiak et al., Estimate of B(B̄ → Xs γ) at O(αs2 ), Phys. Rev. Lett. 98 (2007) 022002 [hep-ph/0609232] [INSPIRE]. [20] M. Misiak and M. Steinhauser, NNLO QCD corrections to the B̄ → X(s)γ matrix elements using interpolation in mc , Nucl. Phys. B 764 (2007) 62 [hep-ph/0609241] [INSPIRE]. [21] S. Descotes-Genon, D. Ghosh, J. Matias and M. Ramon, Exploring new physics in the C7-C70 plane, JHEP 06 (2011) 099 [arXiv:1104.3342] [INSPIRE]. [22] Heavy Flavor Averaging Group collaboration, Averages of b-hadron, c-hadron and τ -lepton properties, arXiv:1010.1589 [INSPIRE]. [23] LHCb collaboration, Test of lepton universality with B 0 → K ∗0 `+ `− decays, JHEP 08 JHEP07(2019)096 (2017) 055 [arXiv:1705.05802] [INSPIRE]. [24] LHCb collaboration, Measurements of the S-wave fraction in B 0 → K + π − µ+ µ− decays and the B 0 → K ∗ (892)0 µ+ µ− differential branching fraction, JHEP 11 (2016) 047 [Erratum ibid. 04 (2017) 142] [arXiv:1606.04731] [INSPIRE]. [25] A. Bharucha, D.M. Straub and R. Zwicky, B → V `+ `− in the standard model from light-cone sum rules, JHEP 08 (2016) 098 [arXiv:1503.05534] [INSPIRE]. [26] A. Khodjamirian, T. Mannel, A.A. Pivovarov and Y.M. Wang, Charm-loop effect in B → K (∗) `+ `− and B → K ∗ γ, JHEP 09 (2010) 089 [arXiv:1006.4945] [INSPIRE]. [27] W. Altmannshofer, P. Stangl and D.M. Straub, Interpreting hints for lepton flavor universality violation, Phys. Rev. D 96 (2017) 055008 [arXiv:1704.05435] [INSPIRE]. [28] W. Altmannshofer, C. Niehoff, P. Stangl and D.M. Straub, Status of the B → K ∗ µ+ µ− anomaly after Moriond 2017, Eur. Phys. J. C 77 (2017) 377 [arXiv:1703.09189] [INSPIRE]. [29] N. Gubernari, A. Kokulu and D. van Dyk, B → P and B → V form factors from B-meson light-cone sum rules beyond leading twist, JHEP 01 (2019) 150 [arXiv:1811.00983] [INSPIRE]. [30] S. Gonzàlez-Solı́s and P. Masjuan, Study of B → π`ν` and B + → η (0) `+ ν` decays and determination of |Vub |, Phys. Rev. D 98 (2018) 034027 [arXiv:1805.11262] [INSPIRE]. [31] S. Descotes-Genon, J. Matias and J. Virto, An analysis of Bd,s mixing angles in presence of new physics and an update of B→ K0∗ K̄0∗ , Phys. Rev. D 85 (2012) 034010 [arXiv:1111.4882] [INSPIRE]. [32] K. De Bruyn et al., Branching ratio measurements of Bs decays, Phys. Rev. D 86 (2012) 014027 [arXiv:1204.1735] [INSPIRE]. (∗)± [33] K. De Bruyn et al., Exploring Bs → Ds K ∓ decays in the presence of a sizable width difference ∆Γs , Nucl. Phys. B 868 (2013) 351 [arXiv:1208.6463] [INSPIRE]. [34] S. Descotes-Genon and J. Virto, Time dependence in B → V `` decays, JHEP 04 (2015) 045 [Erratum ibid. 07 (2015) 049] [arXiv:1502.05509] [INSPIRE]. [35] F. Dettori and D. Guadagnoli, On the model dependence of measured Bs -meson branching fractions, Phys. Lett. B 784 (2018) 96 [arXiv:1804.03591] [INSPIRE]. [36] LHCb collaboration, Angular analysis and differential branching fraction of the decay Bs0 → φµ+ µ− , JHEP 09 (2015) 179 [arXiv:1506.08777] [INSPIRE]. [37] LHCb collaboration, Differential branching fraction and angular analysis of Λ0b → Λµ+ µ− decays, JHEP 06 (2015) 115 [Erratum ibid. 09 (2018) 145] [arXiv:1503.07138] [INSPIRE]. [38] LHCb collaboration, Angular moments of the decay Λ0b → Λµ+ µ− at low hadronic recoil, JHEP 09 (2018) 146 [arXiv:1808.00264] [INSPIRE]. – 34 – [39] Belle collaboration, Improved measurement of the electroweak penguin process B → Xs l+ l− , Phys. Rev. D 72 (2005) 092005 [hep-ex/0503044] [INSPIRE]. [40] BaBar collaboration, Measurement of the B → Xs l+ l− branching fraction and search for direct CP-violation from a sum of exclusive final states, Phys. Rev. Lett. 112 (2014) 211802 [arXiv:1312.5364] [INSPIRE]. [41] LHCb collaboration, Measurement of the Bs0 → µ+ µ− branching fraction and effective lifetime and search for B 0 → µ+ µ− decays, Phys. Rev. Lett. 118 (2017) 191801 [arXiv:1703.05747] [INSPIRE]. [42] A. Lenz et al., Anatomy of new physics in B − B̄ mixing, Phys. Rev. D 83 (2011) 036004 [arXiv:1008.1593] [INSPIRE]. JHEP07(2019)096 [43] J. Charles et al., Predictions of selected flavour observables within the standard model, Phys. Rev. D 84 (2011) 033005 [arXiv:1106.4041] [INSPIRE]. [44] A. Lenz et al., Constraints on new physics in B − B̄ mixing in the light of recent LHCb data, Phys. Rev. D 86 (2012) 033008 [arXiv:1203.0238] [INSPIRE]. [45] J. Charles, S. Descotes-Genon, V. Niess and L. Vale Silva, Modeling theoretical uncertainties in phenomenological analyses for particle physics, Eur. Phys. J. C 77 (2017) 214 [arXiv:1611.04768] [INSPIRE]. [46] ALEPH, DELPHI, L3, OPAL, SLD, LEP Electroweak Working Group, SLD Electroweak Group, SLD Heavy Flavour Group collaboration, Precision electroweak measurements on the Z resonance, Phys. Rept. 427 (2006) 257 [hep-ex/0509008] [INSPIRE]. [47] L. Demortier and L. Lyons, Everything you always wanted to know about pulls, CDF Note 5776 (2002). [48] B. Capdevila, U. Laa and G. Valencia, Anatomy of a six-parameter fit to the b → s`+ `− anomalies, Eur. Phys. J. C 79 (2019) 462 [arXiv:1811.10793] [INSPIRE]. [49] J. Matias, F. Mescia, M. Ramon and J. Virto, Complete anatomy of B̄d → K̄ ∗0 (→ Kπ)l+ l− and its angular distribution, JHEP 04 (2012) 104 [arXiv:1202.4266] [INSPIRE]. [50] W. Altmannshofer et al., Symmetries and asymmetries of B → K ∗ µ+ µ− decays in the standard model and beyond, JHEP 01 (2009) 019 [arXiv:0811.1214] [INSPIRE]. [51] B. Capdevila, S. Descotes-Genon, J. Matias and J. Virto, Assessing lepton-flavour non-universality from B → K ∗ `` angular analyses, JHEP 10 (2016) 075 [arXiv:1605.03156] [INSPIRE]. [52] G. D’Amico et al., Flavour anomalies after the RK ∗ measurement, JHEP 09 (2017) 010 [arXiv:1704.05438] [INSPIRE]. [53] L.-S. Geng et al., Towards the discovery of new physics with lepton-universality ratios of b → s`` decays, Phys. Rev. D 96 (2017) 093006 [arXiv:1704.05446] [INSPIRE]. [54] A. Arbey, T. Hurth, F. Mahmoudi and S. Neshatpour, Hadronic and new physics contributions to b → s transitions, Phys. Rev. D 98 (2018) 095027 [arXiv:1806.02791] [INSPIRE]. [55] G. Hiller and I. Nisandzic, RK and RK ∗ beyond the standard model, Phys. Rev. D 96 (2017) 035003 [arXiv:1704.05444] [INSPIRE]. [56] A.K. Alok et al., New Physics in b → sµ+ µ− after the Measurement of RK ∗ , Phys. Rev. D 96 (2017) 095009 [arXiv:1704.07397] [INSPIRE]. [57] C. Bobeth, M. Chrzaszcz, D. van Dyk and J. Virto, Long-distance effects in B → K ∗ `` from analyticity, Eur. Phys. J. C 78 (2018) 451 [arXiv:1707.07305] [INSPIRE]. – 35 –