where
\(BW(m)=\frac{\sqrt{12\pi \Gamma _{e^+e^-}\Gamma _{\rm tot}\times {\cal B}(\psi \rightarrow p\bar{p}\pi ^0)}}{m^2-M^2+i M\Gamma _{\rm tot}}\)  is a Breit-Wigner function to describe the
resonance amplitude with a mass \(M\) ,
total width \(\Gamma _{\rm tot}\) , electronic partial width \(\Gamma _{e^+e^-}\) , and
the branching fraction \({\cal B}(\psi \rightarrow p\bar{p}\pi ^0)\) .
The continuum term \(\sigma _{\rm con}\)  is defined the same as before.
In the fits, both the mass and the total width
of the resonances are fixed according to the PDG [1]},
and the product \(\Gamma _{e^+e^-}\times {\cal B}(\psi \rightarrow p\bar{p}\pi ^0)\)  is a free parameter.
The continuum term \(\sigma _{\rm con}\)  and the relative phase between
the continuum and the resonance amplitudes, \(\phi \) , are also float parameters.
