Quark--antiquark states and their radiative transitions in terms of the   spectral integral equation. {\Huge II.} Charmonia
  In the precedent paper of the authors (hep-ph/0510410), the $b\bar b$ states were treated in the framework of the spectral integral equation, together with simultaneous calculations of radiative decays of the considered bottomonia. In the present paper, such a study is carried out for the charmonium $(c\bar c)$ states. We reconstruct the interaction in the $c\bar c$-sector on the basis of data for the charmonium levels with $J^{PC}=0^{-+}$, $1^{--}$, $0^{++}$, $1^{++}$, $2^{++}$, $1^{+-}$ and radiative transitions $\psi(2S)\to\gamma\chi_{c0}(1P)$, $\gamma\chi_{c1}(1P)$, $\gamma\chi_{c2}(1P)$, $\gamma\eta_{c}(1S)$ and $\chi_{c0}(1P)$, $\chi_{c1}(1P)$, $\chi_{c2}(1P)\to\gamma J/\psi$. The $c\bar c$ levels and their wave functions are calculated for the radial excitations with $n\le 6$. Also, we determine the $c\bar c$ component of the photon wave function using the $e^+e^-$ annihilation data: $e^+e^- \to J/\psi(3097)$, $\psi(3686)$, $\psi(3770)$, $\psi(4040)$, $ \psi(4160)$, $\psi(4415)$ and perform the calculations of the partial widths of the two-photon decays for the $n=1$ states: $\eta_{c0}(1S)$, $\chi_{c0}(1P)$, $\chi_{c2}(1P)\to\gamma\gamma$, and $n=2$ states: $\eta_{c0}(2S)\to\gamma\gamma$, $\chi_{c0}(2P)$, $\chi_{c2}(2P)\to \gamma\gamma$. We discuss the status of the recently observed $c\bar c$ states X(3872) and Y(3941): according to our results, the X(3872) can be either $\chi_{c1}(2P)$ or $\eta_{c2}(1D)$, while Y(3941) is $\chi_{c2}(2P)$. 
V.V. Anisovich, L.G. Dakhno, M.A. Matveev, V.A. Nikonov and A.V.   Sarantsev