Electron–positron annihilation into hadrons accompanied by an energetic photon provides a powerful tool to measure the hadronic cross-section across a broad energy range at high-luminosity flavour factories such as DAPHNE, CESR, PEP-II, KEK-B, SuperKEKB, and BESIII. The Monte Carlo event generator \( \texttt{Phokhara} \) has been widely used and simulates this process with next-to-leading order (NLO) accuracy, incorporating both virtual and soft photon corrections for single-photon emission, as well as events with two hard real photons.
In this proceedings, we present recent progress in extending \( \texttt{Phokhara} \) beyond NLO precision. The discussion is divided into two parts: (i) we present the plans for adding resummation and refining the hadronic current by incorporating the Generalised Vector Dominance Model as an improvement over the standard scalar QED approximation; and (ii) we outline the preliminary evaluation of two-loop Feynman integrals that form the building blocks of full NNLO predictions.