A global fit of parameters allows us to pin down the Hidden Local Symmetry (HLS) effective Lagrangian, which we apply for the prediction of the leading hadronic vacuum polarization contribution to the muon g − 2. The latter is dominated by the annihilation channel e + e − → π + π − , for which data are available by scan (CMD-2 & SND) and ISR (KLOE-2008, KLOE-2010 experiments. It is well known that the different data sets are not in satisfactory agreement. In fact it is possible to fix the model parameters without using the π + π − data, by using instead the dipion spectra measured in the τ -decays together with experimental spectra for the π 0 γ , ηγ , π + π − π 0 , K + K − , K 0 K 0 final states, supplemented by specific meson decay properties. Among these, the accepted decay width for ρ 0 → e + e − and the partial widths and phase information for the ω/φ → π + π − transitions, are considered. It is then shown that, relying on this global data set, the HLS model, appropriately broken, allows to predict accurately the pion form factor below 1.05 GeV. It is shown that the data samples provided by CMD-2, SND and KLOE-2010 behave consistently with each other and with the other considered data. Consistency problems with the KLOE-2008 and BaBar data samples are substantiated. "All data" global fits are investigated by applying reweighting the conflicting data sets. Constraining to our best fit, the broken HLS model yields a th μ = (11 659 169.55 + +1.26 −0.59 φ + +0.00 −2.00 τ ± 5.21 th ) 10 −10 associated with a very good global fit probability. Correspondingly, we find that a μ = a exp μ − a th μ exhibits a significance ranging between 4.7 and 4.9σ .