In this paper, we study the processes of e+e− annihilation into two body final states, either two pseudoscalar mesons or one meson with a photon. The hadronic vacuum polarization form factors are calculated within the framework of resonance chiral theory in the energy region of E ≲ 2 GeV, with final state interactions taken into account. A joint analysis on the processes of e+e−→ π+π−, K+K−, $$ {K}_L^0{K}_S^0 $$
K
L
0
K
S
0
, π0γ, and ηγ has been performed, and the latest experimental data are included. Based on the vacuum polarization form factors of these processes, we estimate their contributions to the lowest order of anomalous magnetic moment of the muon, (g – 2)μ. Combined with other contributions from hadronic vacuum polarization and other interactions from the standard model, the discrepancy between theoretical prediction and experimental measurement is ∆aμ = (24.1 ± 5.4) × 10−10, i.e., 4.5σ.