An important role in applications of quantum chromodynamics to hadronic physicsrnis played by charmonium and bottomomium which are built up of a heavy quarkrnand heavy antiquark. In this work, we calculate the mass spectrum of charmoniumrnfor 1P, ..., 3P states of 0++ , 1++ and 1+À€€, as well as for 1S, ..., 4S states of 0ð€€€+, andrn1S, ..., 5D states of 1ð€€€ð€€€ along with the two-photon decay widths of ground and firstrnexcited state of 0++ quarkonia for the process, 0++ ! in the framework of arnQCD motivated Bethe-Salpeter Equation. In this 4_4 BSE framework, the coupledrnSalpeter equations are first shown to decouple for the confining part of interaction,rnunder heavy-quark approximation, and analyically solved, and later the one-gluonexchangerninteraction is perturbatively incorporated leading to mass spectral equationsrnfor various quarkonia. The analytic forms of wave functions obtained are usedrnfor calculation of two-photon decay widths of _c0. Our results are in reasonablernagreement with data (where ever available) and other models. We derive the massrnspectrum of the meson anallytically by decoupling the coupled integral equationrnshowing that the incorporation of coulomb potential purterbatively brings our numericalrnvalue much more closer to the experimental value and lift the degeneracyrnof the spectrum.