The decay $Hto ell^+ell^- gamma$, $ell=e,mu$, receives contributions from $Hto Z[to ell^+ell^-] gamma$ and a non-resonant contribution, both of which are loop-induced. We describe how one can separate these sub-processes in a gauge-independent way, define the decay rate $Gamma(Hto Zgamma)$, and extract the latter from differential $Hto ell^+ell^- gamma$ branching ratios. For $ell=mu$ also the tree decay rate, which is driven by the muon Yukawa coupling, is important. We propose kinematic cuts optimized to separate the three contributions, paving the way to the milestones (i) discovery of $Hto Z gamma$, (ii) discovery of $Hto mu^+mu^- left. gammaright|_{rm tree}$, and (iii) quantification of new physics in both the effective $H$-$Z$-$gamma$ and non-resonant $H$-$ell^+$-$ell^-$-$gamma$ couplings.