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

TL;DR

This study reports the first observation of the $X(3872)$ resonance in the $ extomega J/ extomega$ system, measures its decay ratio, and investigates related resonances in $e^+e^-$ collisions at BESIII, revealing new insights into exotic charmonium-like states.

Contribution

First observation of $X(3872)$ in $ extomega J/ extomega$ system and detailed measurement of its properties and related resonances at BESIII.

Findings

$X(3872)$ observed with >5σ significance in $ extomega J/ extomega$.

Measured decay ratio $ ext{R}=1.6^{+0.4}_{-0.3} ext{(stat)} ext{±}0.2 ext{(sys)}$.

Identified resonances $Y(4200)$ and $X(3915)$ with measured masses and widths.

Abstract

We study the $e^+e^-\to\gamma\omega 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 $5\sigma$. The relative decay ratio of $X(3872)\to\omega 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^-\to\gamma 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.4\pm2.2\pm1.2~\rm MeV/c^2$ and $3.8\pm7.5\pm2.6~\rm MeV$, or $3932.6\pm8.7\pm4.7~\rm MeV/c^2$ and $59.7\pm15.5\pm3.7~\rm MeV$, respectively, depending on the fit models. The resonant parameters of the $X(3915)$ agree with those of the $Y(3940)$ in $B\to K\omega J/\psi$ and of the $X(3915)$ in $\gamma\gamma\to\omega J/\psi$ by the Belle and BABAR experiments within errors.

Tables

Table 1. Table 1: The masses (in MeV/ c 2 superscript 𝑐 2 c^{2} ) and widths (in MeV) of the X ​ ( 3872 ) 𝑋 3872 X(3872) ,  X ​ ( 3915 ) 𝑋 3915 X(3915) , and X ​ ( 3960 ) 𝑋 3960 X(3960) resonances from the fit. The numbers in brackets represent the fit scenario without the X ​ ( 3960 ) 𝑋 3960 X(3960) . The uncertainties are statistical only.

	Mass	Width
$X\left(3872\right)$	$3873.3\pm 1.1\left(3872.8\pm 1.2\right)$	$1.2\left(1.2\right)$
$X\left(3915\right)$	$3926.4\pm 2.2$ $\left(3932.6\pm 8.7\right)$	$3.8\pm 7.5$ $\left(59.7\pm 15.5\right)$
$X\left(3960\right)$	$3963.7\pm 5.5$	$33.3\pm 34.2$

Table 2. Table 2: Summary of the systematic uncertainties for X ​ ( 3872 ) 𝑋 3872 X(3872) / X ​ ( 3915 ) 𝑋 3915 X(3915) / X ​ ( 3960 ) 𝑋 3960 X(3960) mass and width measurements. The numbers in the brackets correspond to the fit scenario with only the X ​ ( 3872 ) 𝑋 3872 X(3872) and X ​ ( 3915 ) 𝑋 3915 X(3915) as signal PDF.

Source	Mass (MeV/$c^2$)	Width (MeV)
Absolute mass scale	0.8/0.8 (0.8)/0.8	-/-/-
Background shape	0.3/0.4 (4.5)/0.5	-/2.5 (3.6)/8.3
Resolution	0.0/0.8 (0.7)/0.8	-/0.7 (0.3)/0.1
Fit model	0.5/-/-	-/-/-
Total	1.0/1.2 (4.7)/1.3	-/2.6 (3.7)/8.3

Full text

The $X(3872)$ resonance was first observed by the Belle experiment bellex , and confirmed by the CDF CDFx , D0 D0x , BABAR babarx , LHCb LHCbx , and BESIII Collaborations bes3x . Its unusual properties do not accommodate with a charmonium state, and thus, the $X(3872)$ resonance is widely explained as an unconventional meson candidate review . Since the $X(3872)$ mass is near the $\bar{D^{0}}D^{*0}$ mass threshold, it is often interpreted as a hadronic molecule by theoretical models molecule . The hadronic molecule model predicts that the decay of $X(3872)\to\omega J/\psi$ is sensitive to its internal structure, and a precise measurement of the decay rate would help to determine the ratio of various components that contribute to the $X(3872)$ wave function. While the decay $X(3872)\to\pi^{+}\pi^{-}J/\psi$, where $\pi^{+}\pi^{-}$ is found to be dominated by a $\rho^{0}$ CDF-pp , violates the isospin symmetry in the strong interaction, the $X(3872)\to\omega J/\psi$ decay process conserves isospin symmetry, and thus such a decay provides an excellent metric for probing its isospin-violation effect. Previously, the Belle and BABAR Collaborations only reported less than $5\sigma$ evidences for the $X(3872)\to\omega J/\psi$ decay wjpsi . A solid observation is still lacking and necessary for improved interpretation of this first experimentally observed state potentially composed of four quarks.

The BESIII Collaboration recently reported evidence for the radiative transition $Y(4260)\to\gamma X(3872)$ in $X(3872)\to\pi^{+}\pi^{-}J/\psi$ mode bes3x . A charged charmoniumlike state $Z_{c}(3900)$, which is a good candidate for a four-quark state four-quark , was observed near $\sqrt{s}=4.26$ GeV by BESIII bes3-zc and Belle belle-zc , and later confirmed with CLEO-c’s data at $\sqrt{s}=4.17$ GeV cleo-zc . All these observations show potential connections among the $X(3872)$, $Y(4260)$ and $Z_{c}(3900)$ resonances, and strongly hint towards a common underlying nature for them. At the moment, more supportive experimental observation for the transition process $Y(4260)\to\gamma X(3872)$ is needed to establish these connections.

The $Y(3940)$ resonance was observed by the Belle Collaboration belle-y3940 and confirmed by the BABAR Collaboration babar-y3940 in $B\to K\omega J/\psi$. Later on, both Belle and BABAR reported observations of the $X(3915)$ resonance in $\gamma\gamma\to\omega J/\psi$ process x3915 , and it was suggested to be the same resonance as the $Y(3940)$ by the Particle Data Group (PDG) pdg . The underlying nature of the $X(3915)$ is still unclear. It was once considered as a candidate for the $\chi_{c0}(2P)$ charmonium state. However, such kind of assignment was challenged by a recent Belle observation chic0p . Other interpretations, such as a tetraquark tetra-x3915 or a hadronic molecule molecule-x3915 are proposed for the $X(3915)$. Morever, a theoretical calculation predicted a $1^{++}$ tetraquark with mass near 3.95 GeV/$c^{2}$ tetra-liuyr . To make the situation more clear, it is important to provide additional data on the $X(3915)$.

In this Letter, we report the study of the process $e^{+}e^{-}\to\gamma\omega J/\psi$, with $J/\psi\to\ell^{+}\ell^{-}$ ($\ell=e,~{}\mu$) and $\omega\to\pi^{+}\pi^{-}\pi^{0}(\pi^{0}\to\gamma\gamma)$, using data samples collected with the BESIII detector Ablikim:2009aa . We search for the $X(3872)$ and $X(3915)$ resonances in the $\omega J/\psi$ system and study the $\sqrt{s}$-dependent production cross section, $\sigma[e^{+}e^{-}\to\gamma X(3872)]$. The $e^{+}e^{-}$ center-of-mass (CM) energies of the data sets range from $\sqrt{s}=4.008$ to 4.600 GeV (c.f. Supplemental Material xsec-CM ), with a total integrated luminosity of about 11.6 fb*-1*.

The BESIII detector is described in detail elsewhere Ablikim:2009aa ; etof . geant4 geant based Monte Carlo (MC) simulation samples are used to optimize the event selection criteria, determine the detection efficiency, and estimate backgrounds. For the signal process, we generate $e^{+}e^{-}\to\gamma X(3872)/X(3915)\to\gamma\omega J/\psi$ MC events, with $J/\psi\to\ell^{+}\ell^{-}$ ($\ell=e,~{}\mu$) and $\omega\to\pi^{+}\pi^{-}\pi^{0}(\pi^{0}\to\gamma\gamma)$ at each CM energy corresponding to data. The $X(3872)/X(3915)\to\omega J/\psi$ decay is described with the phase-space model from evtgen evtgen . Initial-state-radiation (ISR) is simulated with kkmc kkmc . The maximum ISR photon energy is set to correspond to the 3.90 GeV/$c^{2}$ production threshold of the $\gamma X(3872)$ system. The final-state-radiation (FSR) from charged final-state particles are handled with photos photos .

Events with four charged tracks with net zero charge are selected. For each charged track, the polar angle in the multilayer drift chamber must satisfy $|\cos\theta|<0.93$, and the point of closest approach to the $e^{+}e^{-}$ interaction point must be within $\pm 10$ cm in the beam direction and within $1$ cm in the plane perpendicular to the beam direction. Since the $\pi^{\pm}$ from $\omega$ decay and $\ell^{\pm}$ from $J/\psi$ decay are kinematically well separated, charged tracks with momenta larger than 1.0 GeV/$c$ in the laboratory frame are assumed to be $\ell^{\pm}$, and the ones with momenta less than 1.0 GeV/$c$ are assumed to be $\pi^{\pm}$. The energy deposition of charged tracks in the electromagnetic calorimeter (EMC) is used to separate $e$ and $\mu$. For $\mu^{\pm}$ candidates, the deposited energy in the EMC are required to be less than 0.35 GeV, while for $e^{\pm}$, it is required to be larger than 1.1 GeV.

Showers identified as good photon candidates must satisfy fiducial and shower-quality requirements. The minimum EMC energy is 25 MeV for barrel showers ($|\cos\theta|<0.80$) and 50 MeV for end-cap showers ($0.86<|\cos\theta|<0.92$). To eliminate showers produced by charged particles, a photon must be separated by at least 20 degrees from any charged track in the EMC. The time information from the EMC is also used to suppress electronic noise and energy deposits unrelated to the event. At least three good photon candidates are required in each event.

To improve the momentum and energy resolutions and to reduce backgrounds, a five-constraint (5C) kinematic fit is applied to an event with the hypothesis $e^{+}e^{-}\to\gamma\pi^{+}\pi^{-}\pi^{0}\ell^{+}\ell^{-}$, which constrains the sum of four momentum of the final-state particles to the initial colliding beams, and the mass of two photon combinations to the $\pi^{0}$ world average mass pdg . The $\chi^{2}$ over number of degree of freedom (ndf) of the kinematic fit is required to be less than $100/5$. When there are ambiguities due to multi-combinations or multi-photon candidates in one event, we choose the combination with the smallest $\chi^{2}$.

Background events such as $e^{+}e^{-}\to\pi^{+}\pi^{-}\psi(3686)/\pi^{0}\pi^{0}\psi(3686)\to\pi^{+}\pi^{-}\pi^{0}\pi^{0}J/\psi$ with one photon candidate missing would also pass the previously described event selection. To remove these backgrounds, we require $|M^{\rm recoil}(\pi^{+}\pi^{-})-m[\psi(3686)]|>8$ MeV/$c^{2}$ and $|M(\pi^{+}\pi^{-}J/\psi)-m[\psi(3686)]|>7$ MeV/$c^{2}$, where $M^{\rm recoil}(\pi^{+}\pi^{-})=\sqrt{(P_{e^{+}e^{-}}-P_{\pi^{+}\pi^{-}})^{2}}$, and $m[\psi(3686)]$ is the mass of the $\psi(3686)$ according to Ref. pdg . Other background events, such as $e^{+}e^{-}\to\eta^{\prime}J/\psi\to\gamma\omega J/\psi$, have the same event topology as the signal. Their contribution can be effectively vetoed by rejecting events satisfying both $0.93<M(\gamma\omega)<0.97$ GeV/$c^{2}$ and $M(\omega J/\psi)>3.9$ GeV/$c^{2}$.

After imposing the above requirements, clear peaks from $J/\psi$ and $\omega$ decays are seen in the $\ell^{+}\ell^{-}$ and $\pi^{+}\pi^{-}\pi^{0}$ invariant mass distributions, as shown in Fig. 1. The $\eta$ peak in the right panel of Fig. 1 comes from $e^{+}e^{-}\to\eta J/\psi$ and $\gamma_{\rm ISR}\psi(3686)\to\gamma_{\rm ISR}\eta J/\psi$ processes. To identify signal candidates that involve the $J/\psi$ resonances, we select events within an invariant mass window of $3.07<M(\ell^{+}\ell^{-})<3.14$ GeV/$c^{2}$, referred to as the $J/\psi$-mass window. Non-$J/\psi$ background events are selected within the two sidebands $2.97<M(\ell^{+}\ell^{-})<3.04$ GeV/$c^{2}$ or $3.17<M(\ell^{+}\ell^{-})<3.24$ GeV/$c^{2}$.

The difference between the mass of $X(3872)$ and $J/\psi$ pdg is about 775 MeV/$c^{2}$, which is slightly lower than the world average mass of the $\omega$. A consequence is an asymmetric $M(\pi^{+}\pi^{-}\pi^{0})$ distribution around the $\omega$ resonance, as can be seen in the right panel of Fig. 1. To accommodate for this effect, the $\omega$ mass window is defined as $0.72<M(\pi^{+}\pi^{-}\pi^{0})<0.81$ GeV/$c^{2}$, and its mass sideband as $0.61<M(\pi^{+}\pi^{-}\pi^{0})<0.70$ GeV/$c^{2}$ or $0.83<M(\pi^{+}\pi^{-}\pi^{0})<0.92$ GeV/$c^{2}$. We fitted both the $M(\ell^{+}\ell^{-})$ and $M(\pi^{+}\pi^{-}\pi^{0})$ distributions, and normalized the data of the sidebands according to the fit results.

Figure 2 shows the $M(\omega J/\psi)$ x-mass distribution from the full data set. A signal peak consistent with the $X(3872)$ resonance is observed. In addition, there are evident structures above 3.9 GeV/$c^{2}$. There are irreducible $e^{+}e^{-}\to\omega\chi_{c0}$ background events that produce a broad structure in the $M(\omega J/\psi)$ distribution. Such kind of background is well understood and can be reproduced by the MC simulation at BESIII wcc0 . Other possible backgrounds come from continuum events, such as $e^{+}e^{-}\to\gamma\omega\pi^{+}\pi^{-}$, $\gamma\pi^{+}\pi^{-}\pi^{0}J/\psi$, $\gamma\pi^{+}\pi^{-}\pi^{0}\pi^{+}\pi^{-}$ etc. They are estimated by analyzing the $J/\psi$ and $\omega$ mass sidebands data.

An unbinned maximum-likelihood fit is performed to the $M(\omega J/\psi)$ mass distribution. In the fit, we use as the signal probability-density-function (PDF), the incoherent sum of three Breit-Wigner (BW) resonances (denoted as $X(3872)$, $X(3915)$ and $X(3960)$, respectively), each convolved with a Gaussian resolution function. The $X(3872)$ width is set to 1.2 MeV pdg . The shape and yield of the $e^{+}e^{-}\to\omega\chi_{c0}$ background component are fixed to the results of the MC simulation. Contribution from other backgrounds is parameterized as a linear shape. The upper panel of Fig. 2 shows the fit results (numerical results are listed in Table 1), and the extracted $X(3872)$ mass agrees with its world average value within errors. The obtained $X(3872)$ signal events yield is $N_{\rm sig}=45\pm 9\pm 3$. The statistical significance of the $X(3872)$ resonance is estimated to be $5.7\sigma$, by comparing the likelihood difference with or without the $X(3872)$ in the fit, $\Delta(-2\ln\mathcal{L})=40.8$, and by taking the change of ndf ($\Delta{\rm ndf}=3$) into account. Possible systematic effects on the $X(3872)$ signal significance, including background shape, $\omega\chi_{c0}$ background normalization, $X(3872)$ intrinsic width and mass resolution are investigated, and no sign for a decreased $X(3872)$ significance is observed. The statistical significance of $X(3915)$ and $X(3960)$ are estimated to be $3.1\sigma$ and $3.4\sigma$ only.

As an alternative choice, we fit the $M(\omega J/\psi)$ mass distribution only with the $X(3872)$ and $X(3915)$ resonances as signal PDF. The $e^{+}e^{-}\to\omega\chi_{c0}$ background is handled in the same way as before. The contribution from other backgrounds is parameterized as a linear function and has been fixed to the result from fitting to the data of the $J/\psi$- and $\omega$-mass sidebands. The bottom panel of Fig. 2 shows the fit results (c.f. Table 1), and the number of fitted $X(3872)$ signal events is $N_{\rm sig}=40\pm 8\pm 3$. The statistical significance of $X(3872)$ is estimated to be $5.2\sigma$, and found to be larger than $5.1\sigma$ after considering systematic effects from $\omega\chi_{c0}$ and linear background normalization, $X(3872)$ intrinsic width and mass resolution. The statistical significance of $X(3915)$ is estimated to be $6.9\sigma$. We test the significance between these two fit scenarios, and find they only differ by $2.5\sigma$.

The production cross section of $e^{+}e^{-}\to\gamma X(3872)$ times the branching fraction $\mathcal{B}[X(3872)\to\omega J/\psi]$ at each CM energy is calculated as $\sigma\cdot\mathcal{B}[X(3872)\to\omega J/\psi]=\frac{N^{\rm sig}}{\mathcal{L}\epsilon(1+\delta)\mathcal{B}}$, where $N^{\rm sig}$ is the number of $X(3872)$ signal events, $\mathcal{L}$ is the integrated luminosity, $\epsilon$ is the detection efficiency, $\mathcal{B}$ is the product of branching fractions for $J/\psi\to\ell^{+}\ell^{-}$ and $\omega\to\pi^{+}\pi^{-}\pi^{0}(\pi^{0}\to\gamma\gamma)$, and $1+\delta$ is the ISR radiative correction factor, which is calculated using the kkmc program kkmc . The ISR photon energy distribution is obtained by an iterative procedure using the line shape $\sigma[e^{+}e^{-}\to\gamma X(3872)]$ measured in this study to replace the default one of kkmc. The left panel of Fig. 3 shows the measured $\sigma\cdot\mathcal{B}[X(3872)\to\omega J/\psi]$. Using the same analysis method as described in Ref. bes3x and the radiative correction factor in this study, $\sigma\cdot\mathcal{B}[X(3872)\to\pi^{+}\pi^{-}J/\psi]$ is measured as well. Our result agrees with and supersedes the earlier published BESIII measurement bes3x , as shown in the right panel of Fig. 3. All the numerical results can be found in Supplemental Materials xsec-CM .

A simultaneous maximum-likelihood fit is performed to both the $\sigma\cdot\mathcal{B}[X(3872)\to\omega J/\psi]$ and the $\sigma\cdot\mathcal{B}[X(3872)\to\pi^{+}\pi^{-}J/\psi]$ distributions. We use a single BW resonance, denoted as $Y(4200)$, with free mass and width as PDF. A free parameter $\mathcal{R}=\frac{\mathcal{B}[X(3872)\to\omega J/\psi]}{\mathcal{B}[X(3872)\to\pi^{+}\pi^{-}J/\psi]}$ is used to describe the relative decay rate of $X(3872)\to\omega J/\psi$ and $\pi^{+}\pi^{-}J/\psi$, which is common for every CM energy. The fit gives $M[Y(4200)]=4200.6^{+7.9}_{-13.3}$ MeV/$c^{2}$, $\Gamma[Y(4200)]=115^{+38}_{-26}$ MeV, $\Gamma^{ee}\cdot\mathcal{B}[Y(4200)\to\gamma X(3872)]\cdot\mathcal{B}[X(3872)\to\pi^{+}\pi^{-}J/\psi]=(4.5^{+1.1}_{-0.8})\times 10^{-2}$ eV and $\mathcal{R}=1.6^{+0.4}_{-0.3}$, where $\Gamma^{ee}$ is the electronic partial width of the $Y(4200)$. Here, all the uncertainties are statistical only.

The systematic uncertainty for $X(3872)$, $X(3915)$, and $X(3960)$ mass and width measurements come from the uncertainties in the absolute mass scale, background and resolution effects. The $e^{+}e^{-}\to\gamma_{\rm ISR}\psi(3686)\to\gamma_{\rm ISR}\eta J/\psi$ events with the same event selection (except the $\omega$ mass window is replaced by the $\eta$ mass window) are used as a control sample to calibrate the mass scale. The measured $\psi(3686)$ mass is $3685.4\pm 0.4$ MeV/$c^{2}$, and the difference to the $\psi(3686)$ world average mass is 0.8 MeV/$c^{2}$. Backgrounds are varied from a linear shape to a second-order polynomial or by $\pm 1\sigma$ for the linear component, and varied by $\pm 1\sigma$ for the $\omega\chi_{c0}$ component in the fit. The differences in the mass and width measurements with respect to the nominal results are taken as a systematic uncertainty. The systematic uncertainty of resolution is estimated by varying the Gaussian parameters of the resolution response function by $\pm 1\sigma$ in the signal PDF. In both fit scenarios (with and without the $X(3960)$), the $X(3872)$ mass difference 0.5 MeV/$c^{2}$ is taken as a systematic uncertainty due to the fit model. All these contributions are summarized in Table 2, and the total uncertainty is calculated by adding the independent contributions in quadrature.

The systematic uncertainty for the $e^{+}e^{-}\to\gamma X(3872)$ cross section measurement mainly comes from uncertainties in the luminosity measurements, detection efficiency, signal extraction, radiative correction and branching fractions. The integrated luminosities of each data set are measured with large-angle Bhabha scattering events, with an uncertainty of 1.0% lum . The tracking efficiency is estimated to be 1% per track from a study of the control sample $J/\psi\to p\bar{p}\pi^{+}\pi^{-}$. The uncertainty due to the photon reconstruction is studied using the $J/\psi\to\pi^{+}\pi^{-}\pi^{0}$ events, and is found to be 1% for the radiative photon gamma-error . An additional systematic uncertainty of 1% is assigned to the efficiency of $\pi^{0}$ reconstruction by studying $\psi(3686)\to\pi^{0}\pi^{0}J/\psi$ and $e^{+}e^{-}\to\omega\pi^{0}$ events. In our event selection, a 5C kinematic fit is used, and the systematic uncertainty related to the kinematic fit is estimated to be 0.8% by using a helix correction method as discussed in Ref. KF .

The number of $X(3872)$ signal events is extracted by fitting the $M(\omega J/\psi)$ distribution, and the difference between the two fit scenarios is 9.5%. The $X(3872)$ intrinsic width is fixed to 1.2 MeV in the signal PDF. Varying the width from 50 keV to 1.2 MeV results in a 5% difference for the $X(3872)$ signal yield. The systematic uncertainty of the $\omega\chi_{c0}$ background is estimated by varying the normalization by $\pm 1\sigma$, which will cause a difference of 0.9% in the $X(3872)$ signal yield. The remaining background is parameterized as a linear function. Varying the background shape from linear to a second-order polynomial or the normalization by $\pm 1\sigma$ will cause a 3.1% difference for the $X(3872)$ signal yield.

We iterate the cross section measurement until the value of $(1+\delta)\epsilon$ changes by at most 1% from the previous iteration, and 1% is taken as a systematic uncertainty due to ISR radiative correction. The systematic uncertainty related to the $J/\psi$-mass window cut is 1.6% bes3x . The branching fraction uncertainties of $J/\psi\to\ell^{+}\ell^{-}$, $\omega\to\pi^{+}\pi^{-}\pi^{0}$ and $\pi^{0}\to\gamma\gamma$ are 0.6%, 0.8% and 0.04% pdg , respectively.

The total systematic uncertainty is calculated to be 12.3% by adding all contributions in quadrature.

The systematic uncertainty for the $Y(4200)$ parameters mainly comes from the uncertainties related to the $e^{+}e^{-}$ CM energy measurement, the parameterization of the fit model, and the cross section measurement. The CM energy of each data set is measured with dimuon events, with $\pm 0.8$ MeV uncertainty ecm . Such kind of common uncertainty will shift the $Y(4200)$ line shape globally, and thus, propagate to the $Y(4200)$ mass linearly. In the fit to cross section, the $Y(4200)$ resonance is parameterized as a BW with a constant full width. We also use a BW with a phase-space dependent full width, $\Gamma\frac{\Phi(\sqrt{s})}{\Phi(M)}$, and the difference is 2.8 MeV/$c^{2}$ for the mass, 12 MeV for the width, and 6.5% for $\Gamma^{ee}$. The cross section data measured in $X(3872)\to\omega J/\psi$ and $\pi^{+}\pi^{-}J/\psi$ channels are fitted simultaneously. The common uncertainties of cross section measurements in both channels, including luminosity, tracking, photon detection, radiative correction, kinematic fit, $X(3872)$ intrinsic width, $J/\psi$ mass window, and $J/\psi\to\ell^{+}\ell^{-}$ branching fraction, will propagate to $\Gamma^{ee}$ linearly, i.e. 6.9%. The uncommon ones, including $\pi^{0}$, background, fit model and $\omega\to\pi^{+}\pi^{-}\pi^{0}(\pi^{0}\to\gamma\gamma)$ branching fraction, will affect the $\mathcal{R}$ measurement, and the total contribution is 10.9%, by adding them in quadrature.

In summary, we have studied the $e^{+}e^{-}\to\gamma\omega J/\psi$ process with 11.6 fb*-1* data at the BESIII experiment. For the first time, the $X(3872)\to\omega J/\psi$ decay was firmly observed with more than $5\sigma$ significance, and the $X(3872)$ mass was measured to be $3873.3\pm 1.1\pm 1.0$ MeV/$c^{2}$. The relative decay ratio for $X(3872)\to\omega J/\psi$ and $\pi^{+}\pi^{-}J/\psi$ is measured to be $\mathcal{R}=1.6^{+0.4}_{-0.3}\pm 0.2$, which agrees well with previous measurements within errors wjpsi . These measurements provide important input for the hadronic molecule interpretation for the $X(3872)$ resonance molecule .

To describe the $M(\omega J/\psi)$ distribution above 3.9 GeV/$c^{2}$, we need at least one additional BW resonance $X(3915)$. Its mass and width are measured to be $3926.4\pm 2.2\pm 1.2$ MeV/$c^{2}$ and $3.8\pm 7.5\pm 2.6$ MeV; or $3932.6\pm 8.7\pm 4.7$ MeV/$c^{2}$ and $59.7\pm 15.5\pm 3.7$ MeV, depending on the fit models.

The $e^{+}e^{-}\to\gamma X(3872)$ production cross section is measured at the CM energies between 4.008 and 4.600 GeV xsec-CM . We studied the $\sqrt{s}$-dependent cross section line shape of $e^{+}e^{-}\to\gamma X(3872)$, and find it can be described by a single BW resonance $Y(4200)$. A simultaneous fit to the $X(3872)\to\omega J/\psi$ and $\pi^{+}\pi^{-}J/\psi$ cross section data gives its mass $M[Y(4200)]=4200.6^{+7.9}_{-13.3}\pm 3.0$ MeV/$c^{2}$, and width $\Gamma[Y(4200)]=115^{+38}_{-26}\pm 12$ MeV, which agree with the $\psi(4160)$ pdg or the $Y(4220)$ observed by BESIII in $\pi^{+}\pi^{-}J/\psi$ bes3-ppjpsi and $\pi^{+}\pi^{-}h_{c}$ bes3-pphc within errors. The measured $e^{+}e^{-}\to\gamma X(3872)$ cross section provides useful information for the $D\bar{D^{*}}$ hadronic molecule calculation as described in Ref. zhaoq-model .

The BESIII collaboration thanks the staff of BEPCII and the IHEP computing center for their strong support. This work is supported in part by National Key Basic Research Program of China under Contract No. 2015CB856700; National Natural Science Foundation of China (NSFC) under Contracts Nos. 11335008, 11425524, 11625523, 11635010, 11735014; the Chinese Academy of Sciences (CAS) Large-Scale Scientific Facility Program; the CAS Center for Excellence in Particle Physics (CCEPP); Joint Large-Scale Scientific Facility Funds of the NSFC and CAS under Contracts Nos. U1532257, U1532258, U1732263; CAS Key Research Program of Frontier Sciences under Contracts Nos. QYZDJ-SSW-SLH003, QYZDJ-SSW-SLH040; 100 Talents Program of CAS; INPAC and Shanghai Key Laboratory for Particle Physics and Cosmology; German Research Foundation DFG under Contract No. Collaborative Research Center CRC 1044; Istituto Nazionale di Fisica Nucleare, Italy; Koninklijke Nederlandse Akademie van Wetenschappen (KNAW) under Contract No. 530-4CDP03; Ministry of Development of Turkey under Contract No. DPT2006K-120470; National Science and Technology fund; The Knut and Alice Wallenberg Foundation (Sweden) under Contract No. 2016.0157; The Swedish Research Council; U. S. Department of Energy under Contracts Nos. DE-FG02-05ER41374, DE-SC-0010118, DE-SC-0012069; University of Groningen (RuG) and the Helmholtzzentrum fuer Schwerionenforschung GmbH (GSI), Darmstadt.