Impact of correlations between aμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_{\mu }$$\end{document} and αQED\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _\text {QED}$$\end{document} on the EW fit

We study the potential impact on the electroweak (EW) fits due to the tensions between the current determinations of the hadronic vacuum polarisation (HVP) contributions to the anomalous magnetic moment of the muon (aμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_{\mu }$$\end{document}), based on either phenomenological dispersion integrals using measured hadronic spectra or on Lattice QCD calculations. The impact of the current tension between the experimental measurement of aμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_{\mu }$$\end{document} and the total theoretical prediction based on the phenomenological calculations of the HVP is also studied. The correlations between the uncertainties of the theoretical predictions of aμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_{\mu }$$\end{document} and of the running of αQED\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _\mathrm{QED}$$\end{document} are taken into account in the studies. We conclude that the impact on the EW fit can be large in improbable scenarios involving global shifts of the full HVP contribution, while it is much smaller if the shift is restricted to a lower mass range and/or if the shift in αQED\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _\text {QED}$$\end{document} is obtained from that in aμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_{\mu }$$\end{document} through appropriate use of the correlations. Indeed, the latter scenarios only imply at most a 2.6/16 increase in the χ2/n.d.f.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\chi ^2/\text {n.d.f.} $$\end{document} of the EW fits and relatively small changes for the resulting fit parameter values.


Introduction
A long-standing discrepancy of about 3-4 standard deviations has been observed between the experimental measurement of the anomalous magnetic moment of the muon (a Exp μ ) [1] and its Standard Model prediction (a SM μ ) [2][3][4][5][6][7][8]. In this comparison, the leading order hadronic vacuum polarisation part (a HVP, LO μ ), derived phenomenologically through dispersion integrals using as input experimental data of e + e − → Hadrons (a HVP, LO μ (Pheno) ), yields the dominant uncertainty of the total theoretical prediction of a μ based on such an approach (a SM (Pheno) μ ). Recently, the BMW collaboration has achieved an unprecedented sub-percent level precision a e-mail: Bogdan.Malaescu@cern.ch b e-mail: Matthias.Schott@cern.ch (corresponding author) for a QCD+QED Lattice calculation of this same contribution (a HVP, LO μ (Lattice) ) [9]. While yielding a reduced tension between the experimental measurement and the theoretical prediction, this new calculation is in tension with the phenomenological one based on dispersion integrals. Recent studies indicate that the latter tension seems to originate from the low energy region (see Ref. [10] and the updated version of Ref. [9]). At the same time, it has been shown in Ref. [11] that the last set of NLO radiative corrections for the pion form factor, not considered previously in the event generators used by the experiments, cannot explain the tension between a Exp μ and a SM μ . Comparisons among the a HVP, LO μ (Lattice) results obtained by various collaborations, as well as with the a HVP, LO μ (Pheno) calculations, have also been performed [12], using in particular a window method with smoothed steps at the boundaries [13][14][15]. The high precision achieved for the recent result of the BMW collaboration motivates its use as reference a HVP, LO μ (Lattice) value in the current study.
It has been advocated in the past that a change of the hadronic spectra (and hence of a HVP, LO μ ) to reduce the tension between the experimental measurement and the theoretical prediction of a μ could introduce tensions in the EW fit [16,17]. More recently, while our work was under completion, it was pointed that a change of the hadronic spectra in the low energy region (below 0.7 GeV) could allow to reduce the tension for a μ without having too strong an impact on the EW fit, 1 although this would be improbable given the current precision of the data [18]. At the same time, in a different study, model-independent bounds were set on the impact that the discrepancy between a HVP, LO μ (Pheno) and a HVP, LO μ (Lattice) can have on the running of α QED to the Z mass (Δα had (M 2 Z )) [19].
We study these aspects taking into account, to our knowledge for the first time, the full correlations between the uncertainties of the HVP contributions to a HVP, LO μ (Pheno) and to Δα had (M 2 Z ). Indeed, these correlations are induced by the use in the two dispersion integrals of the same hadronic spectra, perturbative QCD (pQCD) calculations and narrow resonance contributions. They have been evaluated in Ref. [3], 2 taking into account in particular the correlations of the (statistical and systematic) uncertainties between the different points/bins of a measurement in a given hadronic channel, between different measurements in the same channel, a well as between different channels. This evaluation also fully accounts for the tension between the measurements at the BaBar [22,23] and KLOE [24][25][26][27] experiments in the dominant π + π − (γ ) channel, both through a local re-scaling of the uncertainties by a factor χ 2 /n.d.f. and by taking into account the systematic differences between the two measurements (i.e. comparing the combined dispersion integrals obtained when excluding either BaBar or KLOE).

Description of the EW Fit
The idea of the global EW fit is to compare the state-ofthe-art calculations of EW precision observables with the latest experimental data and thus test the consistency of the Standard Model. The starting point is the EW sector of the Standard Model that can be described by the masses of the EW gauge bosons m V , the mass of the Higgs boson M H , the EW mixing angle θ W , as well as the coupling parameters α em = e 2 /(4π) of the electromagnetic interaction, g and g for the weak interaction as well as the Higgs potential parameter λ. The beauty of the EW theory lies in the predicted relations of its parameters, i.e. the fact that not all of its parameters can be chosen independently from each other. The weak mixing angle, for instance, can be expressed at tree-level as while the mass of the W boson (M W ) is related to the Fermi constant and the fine-structure constant via Hence, at tree level only three free parameters are required. A common choice of the observables, which are used for the predictions, are those with the smallest experimental uncertainties, i.e. the fine structure constant α QED , the Z boson mass M Z and the Fermi constant G F . Knowing these, the observables of the EW sector, in particular M W and sin 2 θ W , can be predicted and confronted with experimental results. However, just using these tree-level relations will lead to immediate incompatibilities with the respective measurements, since higher order EW corrections have to be taken into account. These EW corrections can be formally absorbed into form factors, denoted by κ f Z , ρ f Z 3 and Δr , i.e.
Within the Standard Model, these form factors exhibit a logarithmic dependence on M H , a dependence on quark masses, 4 dominated by a quadratic dependence of the heaviest quark mass m t , and an approximately linear dependence on M Z , α QED and α s . Hence, precise measurements of all observables of the EW sector plus the top quark mass, m t , and α s , allows a test of the consistency of the Standard Model or, alternatively, allows a precise prediction of one observable, when all others are known. This idea of the global EW fit has a long history in particle physics and was performed by several groups in the past, e.g. [28][29][30][31][32][33].
The running of the electromagnetic coupling, α QED , depends crucially on the loop leptonic and hadronic contributions. However, the leptonic and top-quark vacuum polarisation contributions are precisely known or negligible and only the hadronic contribution for the five lighter quarks, Δα (5) had (M 2 Z ) adds significant uncertainties. Hence the electromagnetic coupling α QED is typically replaced by Δα (5) had (M 2 Z ) within the EW fit.
In the following, we use the Gfitter framework [34,35] to evaluate the impact of the Δα (5) had (M 2 Z ) observable in the context of the overall fit. In particular, we indirectly determine Δα (5) had (M 2 Z ) using state of the art measurements of the relevant EW precision observables, but also test its impact on the prediction of other observables such as the W boson mass, M W , the Higgs boson mass, M H , and the effective EW mixing angle, sin 2 (θ e f f ). The Gfitter framework includes for the predictions of M W and sin 2 (θ e f f ) the complete two-loop corrections and allows for a rigorous statistical treatment. For example, it is possible to introduce dependencies among parameters, which can be used to parametrise correlations due to common systematic errors, or to rescale parameter values and errors with newly available results. This is relevant for the study of Δα (5) had (M 2 Z ), as it depends on α S (M 2 Z ). The rescaling mechanism of Gfitter allows to automatically account for arbitrary functional interdependencies between an arbitrary number of parameters [29].

Including the correlations between a μ and α QED in the EW fit
In order to study the impact of the recent a μ -related results on the EW fit, we consider three different approaches. They all involve correlated shifts of the a HVP, LO μ and Δα had (M 2 Z ) values, while taking into account the fact that the kernels involved in the dispersion integrals emphasise lower (higher) energy regions of the hadronic spectra for a HVP, LO μ (Δα had (M 2 Z )). However, the methodology and the underlying assumptions are different for each of the three approaches, which is important in the current context where the source of the tension between the various a μ results is unknown.
The values of the HVP contribution integrals used in this study, computed either through a phenomenological approach or through Lattice QCD, are summarised in Table 1, for either the full HVP contribution or more restricted energy ranges. The latter are starting from the energy threshold (Th.) and go up to either 1 or 1.8 GeV, a region where the sum of 32 exclusive hadronic production channels is used for the phenomenological calculation [3]. It is to be noted that for the a HVP, LO μ dispersion integrals in the phenomenological approach the low energy part dominates for both the central value of the integral and its uncertainty, while for Δα had (M 2 Z ) the high energy regions bring larger contributions. These are direct consequences of the different energy dependencies for the corresponding integration kernels. The correlation coefficient ρ between the uncertainties of the two dispersion integrals (due to the use of the same input hadronic spectra, pQCD and narrow resonance contributions, with different integration kernels) are also indicated, for the various energy ranges that are considered here. It amounts to 44% when computing the dispersion integrals for the full energy range [3] and is further enhanced when considering the contributions from lower mass ranges only. We also note that the full a HVP, LO μ (Pheno) contribution obtained in Ref. [2] through the conservative merging of several results (693.1 ± 4.0) · 10 −10 is similar to the corresponding value from Table 1, in terms of both the central value and uncertainty. In addition, in this study we use the difference between a In the Approach 0 we apply the same scaling factor for the contributions (from some energy range of the hadronic spectrum) to the a HVP, LO μ and Δα had (M 2 Z ) phenomenological values determined from dispersion integrals, in order to reach some "target" value for a HVP, LO μ . This scaling can be modelled e.g. as a change of normalisation of the inclusive hadronic spectrum in the corresponding energy range, which is in this sense similar to the studies done in Refs. [16][17][18]. The EW fit is then performed using as input the shifted Δ α had (M 2 Z ) value, with the corresponding uncertainty. 5 For the Approach 1 the goal is to include a HVP, LO μ in the EW fit, using the information on the correlations between the uncertainties of a HVP, LO μ and Δα had (M 2 Z ). The covariance matrix of the two quantities can be described by a set of two uncertainty components, often called "nuisance parameters" (NPs), each of them being fully correlated between the two quantities, but the two being independent between each-other. There is indeed an infinite number of ways of performing such description of the information in the covariance matrix using two NPs. One set of such NPs that is especially interesting in this case has the format indicated in Table 2, the key point being that NP 1 impacts both quantities, while NP 2 only impacts Δα had (M 2 Z ). One can evaluate the number of standard deviations by which NP 1 has to be shifted, in order for the a HVP, LO μ determined from dispersion integrals to reach some "target" value. The same relative shift of the NP 1 can then be applied to Δα had (M 2 Z ). This shifted value Δ α had (M 2 Z ) is used as input for Gfitter, with the uncertainty provided by the NP 2 (which impacts Δα had (M 2 Z ), but not a HVP, LO μ ). While the full a SM μ prediction, which is directly comparable with a Exp μ , also involves contributions from higher order hadronic loops, hadronic light-by-light scattering, QED and EW effects, for a HVP, LO μ a direct comparison between the phenomenological and Lattice QCD approaches is possible. Without loss of generality, in the current application of the two approaches above, the "target" values of the contribution scaling or uncertainty shift are chosen to bring the a HVP, LO μ (Pheno) contribution derived phenomenologically to the Lattice QCD value a HVP, LO μ (Lattice) , or to bring the a SM μ value to the a Exp μ , or yet to reach these values minus one standard deviation of the corresponding uncertainties. 6 Given the different contribu- 5 We do not apply here a relative scaling of the Δα had (M 2 Z ) uncertainties for Approach 0, because in case such corrections would be necessary for the central values it is not obvious that the uncertainties would be expected to scale accordingly. Furthermore, even in cases when the scaling is applied only to (part of) the range covered by exclusive channels and the scale factor is therefore at the few percent level, the relative impact on the total Δα had (M 2 Z ) uncertainty would be small. 6 For the studies where the "target" value is a Exp μ minus the corresponding uncertainty, we did not include in the definition of this "target" other uncertainty components (e.g. from the light-by-light contribution) involved in the a  and Δα had (M 2 Z ) integrals computed in either the full energy range ("Full HVP") or some restricted region, through either a phenomenological approach using experimen-tal hadronic spectra [3] or with Lattice QCD [9]. Where relevant, ρ indicates the correlation coefficient of the uncertainties of the two phenomenological dispersion integrals tions entering in the various energy ranges involved in the a HVP, LO μ (Pheno) calculation [3], it is difficult to identify a possible effect that would cause a constant global scaling for all of them, although this cannot be fully excluded either. In view also of the indications from Ref. [10] and from the updated version of Ref. [9], it is indeed important to perform the studies of the impact on the EW fit for changes of the hadronic spectra in more restricted energy ranges too. Studies are performed considering scenarios where the contribution scaling or uncertainty shift is done either for the full HVP dispersion integral, or for the sum of the exclusive channels from the energy threshold up to 1.8 GeV, or yet for their contribution up to 1 GeV. 7 These various choices are summarised in Table 3. comparison, exploring the consequences for the EW fit would be too speculative at this stage and remains beyond the goal of our study. We note however that contributions like hadronic light-by-light, impacting a SM μ but not Δα had (M 2 Z ), reduce the correlations between the two and hence the impact of a μ on the EW fit. 7 For the study involving the range between the energy threshold and 1 GeV only the Approach 0 is used, because of the existing correlations between the data uncertainties in this range and in the [1 ; 1.8] GeV interval respectively. Indeed, treating (uncertainties from) the low energy range in Approach 1 independently of the [1 ; 1.8] GeV interval would not be justified, while the coherent treatment of the two intervals would effectively require applying the Approach 1 for the full range up to 1.8 GeV. Note also that the relative uncertainties are also rather similar for the a HVP, LO and by the Lattice QCD value a HVP, LO μ (Lattice) . The correlations between the uncertainties of the a HVP, LO μ (Pheno) and Δα had (M 2 Z ) dispersion integrals for the full energy range [3] are taken into account in the fit. The uncertainty due to the finite precision of α S (entering here through the pQCD part of the dispersion integrals in the phenomenological approach) also impacts other quantities in the EW fit and is therefore treated separately. It amounts to 0.14 · 10 −10 (0.41 · 10 −4 ) for a HVP, LO μ (Pheno) (Δα had (M 2 Z )) and is treated as fully correlated between the two quantities, as well as with the other α S -related uncertainties in the fit. The remaining uncertainties are of 3.96 · 10 −10 for a HVP, LO μ (Pheno) and 0.91 · 10 −4 for Δα had (M 2 Z ), with a correlation coefficient of 47% between the two. In the Lattice QCD calculation, pQCD is used in the range Q 2 > 3 GeV 2 for an a HVP, LO μ contribution that amounts to 0.16 · 10 −10 , while its uncertainty due to the finite precision of α S is negligible [9].
Approach 2 brings a slightly improved treatment of the uncertainties and correlations compared to Approach 1, where the total covariance matrix (including the α S -related uncertainty) has to be used when computing the NP 1 that impacts a HVP, LO μ . There, the α S uncertainty impacting Δα had (M 2 Z ) is treated as a sub-component of NP 2 (and further correlated with other quantities in the EW fit), which effectively decorrelates it from the corresponding uncertainty of a HVP, LO μ . This approximate treatment in Approach 1 is however well justified, given the relatively small contribution of the α S uncertainty to a HVP, LO μ . It is worth noting that in the various scenarios displayed in Table 3 the scaling factors applied in Approach 0 go well beyond the (sub-percent level) systematic uncertainties of the modern experimental measurements of hadronic spectra, used in the phenomenological dispersion integrals. For this reason we also do not consider applying Approach 0 in more restricted energy ranges below 1 GeV, as the resulting scaling factors would be even larger and hence unlikely. Similarly, the shifts of NP 1 (expressed as a number of standard deviations) in Approach 1 are relative large, assuming hence that the Gaussian approximation and the correlation coefficients between the dispersion integrals are still valid in Table 3 Scaling factors of the hadronic spectra in Approach 0, shifts applied to NP 1 (in terms of a number of standard deviations) and the uncertainty to be used in the EW fit σ Δα had (M 2 Z ) (which incorporates NP 2 and the uncertainty from high mass contributions if a restricted range is used for the uncertainty shift, hence including the α S -related uncertainty too) in Approach 1, together with the corresponding modified Δ α had (M 2 Z ) values, for various shifts of the a HVP, LO μ . The latter are achieved through a contribution scaling (Approach 0) or uncertainty shift (Approach 1), applied for various energy ranges of the hadronic spectrum (see text). The "−1σ " following various quantities indicates a subtraction of one standard deviation of the corresponding a HVP, LO μ or a Exp μ uncertainty. For Approach 0 the uncertainty indicated for "Full HVP" in Table 1   . For all these reasons, the current study should not be seen as an attempt to precisely incorporate the a μ inputs into the EW fit, but rather to explore their potential impact under various hypotheses. Indeed, the three approaches (with the various choices listed in Table 3) allow to probe different hypotheses concerning the possible source(s) of the difference between the phenomenological prediction on one side and the recent Lattice QCD result or the experimental measurement on the other, while assessing the corresponding impact on the EW fit. While Approach 0 considers a sim-ple normalisation scaling of the hadronic spectra (and is in this sense comparable with some of the approaches used in Refs. [16][17][18][19]), Approaches 1 and 2 use the information on the experimental uncertainties, with their phase-space dependence and correlations, to guide the evaluation of coherent shifts of a HVP, LO μ and Δα had (M 2 Z ).

Results of the EW Fit
The input parameters of the fit are summarized in Table 4. They include in particular the measurements from the LEP and SLC collaborations, i.e. the mass and width of the Z boson, the hadronic pole cross sections as well as the forward- Table 5 Different input values of used   backward asymmetry parameters. The W boson mass and the top-quark mass are based on measurements at the Tevatron and the LHC, while the Higgs Boson mass is only measured at the latter. In summary, the floating parameters in the global EW fit within the Gfitter program are the coupling parameters Δα (5) had (M 2 Z ) and α S (M 2 Z ), the masses M Z , m c , m b , m t and M H , as well as four theoretical error parameters. In Approach 2, the χ 2 definition is modified to include a HVP, LO μ as an extra free parameter, constrained by both a HVP, LO μ (Pheno) and a HVP, LO μ (Lattice) with the corresponding uncertainties, adding hence one degree of freedom to the fit.
In a first step, we determine the minimal χ 2 of the global EW fit, using various values for Δα (5) had (M 2 Z ) according the different approaches described in Sect. 3. 8 The results are summarized in Table 5. As discussed earlier, in Approach 2 the tension between a HVP, LO μ (Pheno) and a HVP, LO μ (Lattice) induces a contribution to the χ 2 of about 9.3 units.
In a second step, we studied in more detail the a HVP, LO μ (Lattice) − a HVP, LO μ (Pheno) case and used the three different approaches to indirectly determine several selected observables. Technically, this indirect parameter determination is performed by scanning the parameter in a chosen range and calculating the corresponding χ 2 values. The value of χ 2 min is not relevant for the uncertainty estimation, but only its difference relative to the global minimum, Δχ 2 ≡ χ 2 − χ 2 min . The Δχ 2 = 1 and Δχ 2 = 4 profiles define the 1σ and 2σ uncertainties, respectively. The where tensions with the measured values are induced in this scenario, which as discussed above is, however, unlikely. In Approaches 1 and 2, as well as when using Approach 0 with shifts of the HVP contribution applied on more restricted mass ranges, there's less of a change for Δα had (M 2 Z ) and one can conclude that under the corresponding (more realistic) scenarios the impact of the tensions for a μ on the EW fit is small.
The dependence of the predicted value for M H , M W , sin 2 θ l eff and m t on Δα had (M 2 Z ) in the global EW fit is illustrated in Fig. 2. The results of the Approach 0 and Approach 1, applied either for the full HVP contribution or for the range [Th.; 1.8 GeV], for the a HVP, LO μ (Lattice) − a HVP, LO μ (Pheno) case, are also 8 The Δα (5) had (M 2 Z ) values are obtained based on the values in Table 3, after subtracting the contribution of the top quark to the pQCD calculation, which amounts to −0.72 · 10 −4 with a negligible uncertainty. Thirdly, we determine indirectly the value of Δα (5) had (M 2 Z ) (without including any explicit constraint on it in the EW fit) using the other EW observables, including and excluding the Higgs boson mass 9 as well as assuming improved precisions on the EW observables at the end of the high luminosity LHC phase (see Table 4). The corresponding Δχ 2 distributions are shown in Fig. 3 for all three cases, yielding values of Δα (5) had (M 2 Z ) = 0.02716 ± 0.00033 (including m H ), Δα (5) had (M 2 Z ) = 0.02817 ± 0.00087 (excluding m H ) and Δα (5) had (M 2 Z ) = 0.02706 ± 0.00025 (with future measurement precisions), respectively. In addition, we indicate the predicted values of Δα (5) had (M 2 Z ), previously discussed in Sect. 3. The uncertainties on these predictions are estimated 9 The comparison of the χ 2 -distributions when not including the Higgs boson mass in the fit was added to illustrate the impact of the Higgs Boson discovery on the prediction of Δα (5) based on the uncertainty of the "target" value, driven either by the experimental measurement or the Lattice QCD calculation.
Including the constraint on the Higgs boson mass significantly improves the accuracy of the indirect Δα (5) had (M 2 Z ) determination. Then, in all the configurations the tension between the fitted and the predicted Δα had (M 2 Z ) is enhanced, compared to the one for the nominal prediction. However, here also the tension becomes significant only when using the Approach 0 applying a scaling for the full energy range of the hadronic spectrum.

Summary and conclusions
We studied the potential impact on the EW fits of the tensions between the current determinations of the HVP contributions to a μ , based on either phenomenological dispersion integrals of hadronic spectra or respectively on Lattice QCD calculations. Similarly, we also considered the impact of the current tension between the experimental measurement of a μ and its total theoretical prediction based on the phenomenological calculations of the HVP. We considered an approach based on coherent shifts of the hadronic spectra in various mass ranges (comparable with some of the approaches used in Refs. [16][17][18][19]) and two novel approaches that take into account the correlations between the uncertainties of the theoretical predictions of a μ and of the running of α QED in the phenomenological approach. It is found that the impact on the EW fit can be large in scenarios involving global shifts of the full HVP contribution. However, such scenarios seem unrealistic, since they require the same relative shift to be applied in mass ranges of the hadronic contributions where very different methodologies and inputs are being used. Indeed, recent studies where the Δα had (M 2 Z ) contribution is broken into energy intervals also indicate that the differences between the two evaluations, based on the phenomenological and Lattice QCD calculations respectively, seem to originate from the low energy region (see Ref. [10] and the updated version of Ref. [9]). The impact on the EW fit is much smaller if the shift is restricted to a lower mass range (as also noted in Refs. [18,19]) and/or if the shift of the a μ prediction is propagated to α QED following the pattern of the current uncertainties with their full set of correlations, as done in our two novel approaches. In the case of the latter scenarios, addressing the current tensions at the level of a μ would not induce significant tensions in the EW fit, implying at most a 2.6/16 increase in the corresponding χ 2 /n.d.f., while the changes for the resulting fit parameter values are small too. 10 10 In the second version of Ref. [9] a somewhat smaller difference between the phenomenological and Lattice QCD calculations of a HVP, LO μ is observed, reducing even further the potential impact on the EW fit. An improved precision for the experimental measurement of a μ , further precise measurements of the hadronic spectra allowing to hopefully also clarify the tension between BABAR and KLOE in the π + π − channel, as well as other precise Lattice QCD calculations are expected to become available in the future [2]. Beyond the main goal of exploring the possibility of a contribution from new physics in the comparison between the measurement of a μ and its theoretical predictions, the improved precision will allow to better scrutinise the impact of these findings on the EW fit.
Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: The relevant numerical results of the study are already provided in the paper.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .