Studies of mu-pair and pi-pair production at the electron-positron low   energy colliders
  Predictions for the radiative return with the muon pair and pion pair final state from KKMC and PHOKHARA Monte Carlo programs are compared and discussed. The case of muon pairs is well understood, especially of the initial state radiation (ISR), where three different second order calculations agree very well. The case of the final state radiation (FSR) requires more tests. Matrix element in KKMC of the EEX type with the incomplete second order NLL corrections is not good enough for the radiative return at $Q^2<1$GeV with the precision requirement better than 1%. A method of extending the superior CEEX-type matrix element in KKMC to the pion pair final state is described. 
S. Jadach