The large rate of charmless decays of the \(\psi (3770)\)  can also be accommodated
in theoretical models. Among many theoretical efforts trying to solve the
“\(\rho \pi \)  puzzle” in \(J/\psi \)  and \(\psi (3686)\)  decays observed by the
Mark-II experiment [1]}, the 2S-1D charmonium mixing
scenario [2]} relates a partial width in \(\psi (3770)\)  decays
into any final state to the corresponding partial widths in \(J/\psi \) 
and \(\psi (3686)\)  decays [3]}, as both \(\psi (3686)\)  and \(\psi (3770)\) 
are the mixtures of the 2S and 1D charmonium states. As a result,
a large \(\psi (3770)\rightarrow \text{\it light hadrons}\)  is allowed in this model,
and the measurement of the \(\psi (3770)\)  decays can test the
model predictions. The process \(\psi (3770)\rightarrow p\bar{p}\pi ^0\)  is one of the possible
exclusive channels contributing to the charmless decays of the \(\psi (3770)\) .
