We perform a thorough analysis of the $VP\ensuremath{\gamma}({\ensuremath{\gamma}}^{*})$ and $P\ensuremath{\gamma}\ensuremath{\gamma}({\ensuremath{\gamma}}^{*})$ decays in the resonance chiral theory, where $V$ stands for the vector resonances $\ensuremath{\rho}$, ${K}^{*}$, $\ensuremath{\omega}$, and $\ensuremath{\phi}$; $P$ stands for $\ensuremath{\pi}$, $K$, $\ensuremath{\eta}$, and ${\ensuremath{\eta}}^{\ensuremath{'}}$; and ${\ensuremath{\gamma}}^{*}$ subsequently decays into lepton pairs. Upon imposing QCD short-distance constraints on resonance couplings, the $\ensuremath{\omega}\ensuremath{\rightarrow}\ensuremath{\pi}\ensuremath{\gamma}({\ensuremath{\gamma}}^{*})$, $\ensuremath{\rho}\ensuremath{\rightarrow}\ensuremath{\pi}\ensuremath{\gamma}({\ensuremath{\gamma}}^{*})$, ${K}^{*0}\ensuremath{\rightarrow}{K}^{0}\ensuremath{\gamma}$ processes only depend on one free parameter and $\ensuremath{\pi}\ensuremath{\rightarrow}\ensuremath{\gamma}\ensuremath{\gamma}({\ensuremath{\gamma}}^{*})$ can be completely predicted. The four mixing parameters of the $\ensuremath{\eta}--{\ensuremath{\eta}}^{\ensuremath{'}}$ system, i.e. two mixing angles ${\ensuremath{\theta}}_{8}$, ${\ensuremath{\theta}}_{0}$ and two decay constants ${F}_{8}$, ${F}_{0}$, are determined from radiative decays involving $\ensuremath{\eta}$ or ${\ensuremath{\eta}}^{\ensuremath{'}}$. The higher order low energy constants of the pseudo-Goldstone Lagrangian in the chiral anomaly sector are predicted by integrating out heavy resonances. We also predict the decay widths of $\ensuremath{\rho}\ensuremath{\rightarrow}\ensuremath{\pi}{e}^{+}{e}^{\ensuremath{-}}$, ${\ensuremath{\eta}}^{\ensuremath{'}}\ensuremath{\rightarrow}\ensuremath{\gamma}{e}^{+}{e}^{\ensuremath{-}}$ and $\ensuremath{\phi}\ensuremath{\rightarrow}\ensuremath{\eta}{\ensuremath{\mu}}^{+}{\ensuremath{\mu}}^{\ensuremath{-}}$, which can be compared with the future measurement in these channels.