We study the mass spectrum and electromagnetic processes of charmonium system with the spin-dependent potentials fully taking into account in the solution of the Schroedinger equation and the results for the pure scalar and scalar-vector mixing linear confining potentials are compared. It is revealed that the scalar-vector mixing confinement is important for reproducing the mass spectrum and decay widths and the vector component is found to be around 22%, the long-standing discrepancy in M1 radiative transitions of J/ψ and ψ ′ is alleviated by means of the state wave functions obtained via the Hamiltonian with the full spin-dependent potential. This work also intends to identify few of the copious higher charmonium-like states as the cc ones. Particularly, the newly observed X(4160) and X(4350) are assigned as M (2 1 D 2 ) = 4164.9 MeV and M (3 3 P 2 ) = 4352.4 MeV, which strongly favor the J P C = 2 −+ , 2 ++ assignments respectively. The corresponding radiative transitions, leptonic and two-photon decay widths have been also calculated for the further experimental study.Keywords charmonium · potential model · Gaussian expansion method · XYZ mesons 1 IntroductionDue to the impressive increase of experimental results, charmonium (cc) spectroscopy has renewed great interest recently, coming along with the striking disagreement with theoretical expectations [1,2,3]. The unexpected and still-fascinating X(3872) has been joined by more than a dozen other charmonium-like states, while the series of vacancy have been left on the cc list. It is urgent to identify the possible new members of charmonium family from the abundant observations. The QCD inspired potential models have been playing an important role in investigating heavy quarkonium, owning to the presence of large nonperturbative effects in this energy region. Most quark potential models [4,5,6,7,8,9,10,11,12,13,14,15,16] have common ingredients under the nonrelativistic limit, despite some differences in the detailed corrections for relativistic and coupled channel effects, which typically are the Coulomb-like term induced by one-gluon exchange plus the long-range confining potential expected from nonperturbative QCD. Anyway, the nature of confining mechanism has been veiled so far. In the original Cornell model [17,18], it was assumed as Lorentz scalar, which gives a vanishing long-range magnetic contribution and agrees with the flux tube picture of quark