Comments on “An update of the HLS estimate of the muon g-2” by M. Benayoun et al., EPJC (2013) 73:2453

In a recent paper, Benayoun et al. (Eur. J. Phys. C 73:2453, 2013) use a specific model to compare results on the existing data for the cross section of the process e+e−→π+π− and state conclusions about the inconsistency of the BABAR results with those from the other experiments. We show that a direct model-independent comparison of the data at hand contradicts this claim. Clear discrepancies with the results of Benayoun et al. (Eur. J. Phys. C 73:2453, 2013) are pointed out. As a consequence we do not believe that the lower value and the smaller uncertainty obtained for the prediction of the muon magnetic anomaly are reliable results.


Introduction
The authors of Ref. [1] use a specific model of the cross section for the process e + e − → π + π − in order to compute the corresponding contribution to the muon magnetic anomaly. Specifically they fix parameters of a theoretical model based on the concept of Hidden Local Symmetry (references to be found in Ref. [1]) using the spectral function from τ decays obtained by the Belle collaboration [2] and other less precise inputs from additional processes assuming vector dominance.
In these comments we consider four points relevant to the analysis in Ref. [1] where the model is used to discriminate between the different sets of recent e + e − data, allowing the authors to make statements about their consistency. First, their model for the e + e − cross section relies on the τ spectral function from Belle within the framework of the HLS approach. Such a procedure relies on assumptions regarding the isospin-breaking corrections needed in order to a e-mail: davier@lal.in2p3.fr b e-mail: malaescu@in2p3.fr link the τ and e + e − spectral functions. We provide some comments on this point in Sect. 2. Second, the model is used by the authors to compare the results of the model with existing measurements of the π + π − annihilation process. In the peak ρ region they claim good agreement between the data from KLOE [3,4], CMD-2 [5,6], and SND [7], while pointing out a discrepancy with BABAR [8,9]. As we show in Sect. 3 this conclusion is not consistent with the direct comparison of these data.
Third, still using the model, a specific analysis of the ρ − ω interference pattern in the above experimental data is performed, again showing an even bigger discrepancy between BABAR and the other experiments. However the ρ − ω interference has been already extracted by the experiments (but KLOE) and the published results clearly do not support the claim made by the authors, as demonstrated in Sect. 4.
Fourth, we comment on the values obtained for the hadronic contribution to the muon g − 2 and on the stated increased significance of the discrepancy with the direct measurement [10] which we believe is not warranted.
In the following we address these four issues. For the second and third points, we show that in a model-independent way the direct comparison of the measured cross sections yields results which are at variance with those of Ref. [1].

Isospin breaking corrections
In addition to the standard isospin-breaking (IB) corrections [11][12][13][14][15][16][17][18], M. Benayoun et al. consider the effect of γ − ρ mixing, already proposed previously [19]. The effect of mixing is not shown quantitatively in Ref. [1], but from the results in Ref. [19] one can observe that a major consequence in the ρ region is a mass shift between neutral and charged ρ resonances, although some additional effects are seen in the tails. A mass shift was already considered in the analysis of Ref. [18], together with electromagnetic contributions in the ρ decays, with the conclusion that the isospin-corrected τ spectral function agreed well with the BABAR data and rather well with the Novosibirsk results, while disagreeing with KLOE.
The authors of Ref. [1] reach the opposite conclusion, i.e. good agreement with KLOE, fair agreement with CMD-2 and SND, and discrepancy with BABAR. We point out that even if the corrections used follow the same rules as in Ref. [19], their final results differ significantly. Therefore it seems that the IB corrections as used by the authors are still subject to some debate, especially since they do not provide specific studies of their systematic uncertainties. In these conditions it is probably not safe to use the IB-corrected τ spectral function as an exact replica of the e + e − cross section when comparing to the e + e − data.

Consistency of different data in the peak ρ region
The direct comparison of the cross sections for BABAR and the other experiments has been shown in detail in Ref. [9]. There is indeed a discrepancy with KLOE, but agreement with CMD-2 and SND within the quoted respective systematic uncertainties. One can quantify the observations by considering the ratio of the cross sections measured by each experiment to the BABAR result. Since the energy corresponding to the experimental data points do not match the BABAR bins, we perform the quantitative comparison to the form factor fit of the BABAR data, described in Ref. [9] which nicely interpolates the data.
From inspection of their Fig. 5 the authors of Ref.
[1] point out major differences in the ρ peak region 0.70-0.85 GeV. Thus we restrict the ratio fits to this mass range. We exclude the ρ − ω region 0.776-0.788 GeV, discussed in the next section, as the ratios are very sensitive to the relative mass calibration of the experiments (this effect is not taken into account in Ref. [1]). The ratios to BABAR as a function of ππ mass are displayed in Fig. 1. The fit to a constant is a good description of the respective data, except for KLOE where deviations are seen at lower masses, however not affecting significantly the fit over most of the mass range.
The fitted values for the ratios are given in Table 1. Of course the BABAR value in the first row is very close to one by construction, but we show it to demonstrate that the form factor fit, performed in the full mass range from threshold to 3 GeV, does represent well the BABAR data between 0.7 and 0.85 GeV. Only statistical uncertainties are used in the ratio fits, while systematic uncertainties are quoted separately. All the values are displayed in Fig. 2.
From Table 1 and Fig. 2 we see that, unlike the qualitative claim in Ref. [1], the BABAR data in the ρ peak region is fully compatible with CMD-2 and SND, and disagrees with KLOE. The CMD-2 data also disagrees with KLOE, while SND with a larger uncertainty is consistent with both BABAR and KLOE. Since in Ref. [1] CMD-2 and SND are treated together, it is interesting to quote the ratio for the combined Novosibirsk data (taking into account correlated uncertainties). Its value of 0.999 ± 0.007 shows the same disagreement with KLOE as BABAR. Therefore we conclude that the claim by Benayoun et al. is incorrect. There is no justified argument for discarding the BABAR data in the evaluation of the g − 2 dispersion integral.

Analysis of ρ − ω interference
Except for KLOE, all the experiments have provided fits of the pion form factor with ρ and ω amplitudes, even including isovector resonances at higher mass. In these analyses the parameters describing the ρ resonance and the interference amplitude are fitted, so that the product of the branching ratios B ω→ee · B ω→ππ can be directly obtained. Using the fit results published by the experiments, the values deduced for the product are given in Table 2. The values presented in Ref. [1] are listed for comparison. The agreement is reasonable for CMD-2 and SND, however there is a large Table 1 Ratios of the measured e + e − → π + π − cross sections from the quoted experiments relative to the BABAR form factor fit in the mass range 0.7-0.85 GeV, excluding the ρ − ω interference region.
The quoted uncertainties are statistical ('stat'), systematic ('syst'), and total ('tot'). The χ 2 per degree of freedom (DF) of the ratio fits to a constant are given   ates considerably from the direct value and the quoted uncertainty is one order of magnitude too small, given the statistical accuracy of the BABAR data. Thus the conclusion Fig. 2 Ratios of the measured e + e − → π + π − cross sections from the quoted experiments relative to the BABAR cross section fit in the mass range 0.7-0.85 GeV, excluding the ρ − ω interference region. The quoted error bars include both statistical and systematic uncertainties reached in Ref. [1] that the BABAR data strongly disagree with the Novosibirsk experiments is ill-founded.

The muon magnetic anomaly
The authors of Ref.
[1] present a range of values indicating a discrepancy with the direct measurement [10], with their 'best' estimate given at 4.9σ . This result follows from two facts: the prediction has a lower central value, as the result of discarding the BABAR data, and its uncertainty is reduced because their model brings additional constraints. In comparison the most recent analyses [20,21] are model-independent and use all the available e + e − data, taking into account the existing discrepancies between experiments in order to set realistic uncertainties. In that case the discrepancy in the muon magnetic anomaly is only at 3.3-3.7σ .

Conclusions
We have pointed out serious discrepancies in the modeldependent analysis of Ref. [1] when compared to the original e + e − → π + π − cross section data. Let us summarize our points: • The approach relies crucially on model-dependent isospin-breaking corrections involved in relating the τ spectral function to the e + e − cross section. A detailed study of the corresponding systematic uncertainties, mandatory before using the model to discriminate between various data sets, is lacking. • The direct and quantitative comparison of the different e + e − data sets contradicts the claim made that BABAR data disagree with the other experiments. • The direct BABAR analysis of ρ − ω interference gives a result in agreement with Novosibirsk data, in contradiction with the value given in Ref. [1] which makes no sense and is again used to discriminate against the BABAR data. • As a consequence we do not believe that the lower value and the smaller uncertainty obtained for the prediction of the muon magnetic anomaly are reliable results.