Study of $e^+e^- o gamma omega J/psi$ and Observation of $X(3872) o omega J/psi$

We study the $e^+e^-togammaomega J/psi$ process using $11.6 ~rm fb^{-1}$ $e^+ e^-$ annihilation data taken at center-of-mass energies from $sqrt{s}=4.008~rm GeV$ to $4.600~rm GeV$ with the BESIII detector at the BEPCII storage ring. The $X(3872)$ resonance is observed for the first time in the $omega J/psi$ system with a significance of more than $5sigma$. The relative decay ratio of $X(3872)toomega J/psi$ and $pi^+pi^- J/psi$ is measured to be $mathcal{R}=1.6^{+0.4}_{-0.3}pm0.2$, where the first error is statistical and the second systematic (the same hereafter). The $sqrt{s}$-dependent cross section of $e^+e^-togamma X(3872)$ is also measured and investigated, and it can be described by a single Breit-Wigner resonance, referred to as the $Y(4200)$, with a mass of $4200.6^{+7.9}_{-13.3}pm3.0~{rm MeV}/c^2$ and a width of $115^{+38}_{-26}pm12~{rm MeV}$. In addition, to describe the $omega J/psi$ mass distribution above $3.9~rm GeV/c^2$, we need at least one additional Breit-Wigner resonance, labeled as $X(3915)$, in the fit. The mass and width of the $X(3915)$ are measured to be $3926.4pm2.2pm1.2~rm MeV/c^2$ and $3.8pm7.5pm2.6~rm MeV$, or $3932.6pm8.7pm4.7~rm MeV/c^2$ and $59.7pm15.5pm3.7~rm MeV$, respectively, depending on the fit models. The resonant parameters of the $X(3915)$ agree with those of the $Y(3940)$ in $Bto Komega J/psi$ and of the $X(3915)$ in $gammagammatoomega J/psi$ by the Belle and BABAR experiments within errors.