Cross section measurements of $e^+ e^-\to\omega\chi_{c0}$ from $\sqrt{s}=$4.178 to 4.278 GeV

TL;DR

This paper measures the energy-dependent cross section of the process $e^+e^- ightarrow ext{omega} ext{chi}_{c0}$ between 4.178 and 4.278 GeV, revealing a resonant structure and providing new angular distribution data.

Contribution

It presents the first measurement of the angular distribution and confirms a resonant structure in the cross section with improved precision.

Findings

Resonant structure with mass ~4218.5 MeV/c^2 and width ~28.2 MeV.

First extraction of the angular distribution for this process.

Improved accuracy over previous measurements.

Abstract

The cross section of the process $e^+e^- \rightarrow \omega \chi_{c0}$ is measured at center-of-mass energies from $\sqrt{s} =$ 4.178 to 4.278 GeV using a data sample of 7 fb$^{-1}$ collected with the BESIII detector operating at the BEPCII storage ring. The dependence of the cross section on $\sqrt{s}$ shows a resonant structure with mass of $(4218.5\pm1.6(\text{stat.})\pm4.0(\text{syst.}))$ MeV/$c^2$ and width of $(28.2\pm3.9(\text{stat.})\pm1.6(\text{syst.}))$ MeV, respectively. This observation confirms and improves upon the result of a previous study. The angular distribution of the $e^+e^- \rightarrow \omega \chi_{c0}$ process is extracted for the first time.

Tables

Table 1. Table 1: Born cross sections σ B superscript 𝜎 B \sigma^{\rm B} (or upper limits at 90% C.L. σ upper B subscript superscript 𝜎 B upper \sigma^{\rm B}_{\text{upper}} ) for the e + ​ e − → ω ​ χ c ​ 0 → superscript 𝑒 superscript 𝑒 𝜔 subscript 𝜒 𝑐 0 e^{+}e^{-}\to\omega\chi_{c0} reaction at the different center-of-mass energies s 𝑠 \sqrt{s} , together with integrated luminosities ℒ ℒ \mathcal{L} , number of signal events N sig superscript 𝑁 sig N^{\rm sig} , radiative correction factor 1 + δ 𝛿 \delta (s), vacuum polarization factor 1 | 1 − Π | 2 1 superscript 1 Π 2 \frac{1}{|1-\Pi|^{2}} , and efficiency ϵ italic-ϵ \epsilon . The first uncertainties are statistical, and the second systematic.

$\sqrt{s}$ (GeV)	$ℒ$ (pb-1)	$N^{\mathrm{sig}}$	1 + $\delta$(s)	$\frac{1}{{\left|1-\Pi \right|}^2}$	$ϵ$ $(\%)$	${\sigma }^B\left({\sigma }_{\text{upper}}^B\right)$ (pb)
4.178	3194.5	${0.0}_{-0.0}^{+11.6}$	0.63	1.055	24.71
4.189	524.6	$5.4\pm 4.7$	0.64	1.056	24.59
4.199	526.0	$21.5\pm 6.4$	0.66	1.057	25.68	$22.6\pm 6.7\pm 2.6$
4.209	518.0	$27.8\pm 8.4$	0.68	1.057	25.70	$28.8\pm 8.7\pm 4.3$
4.219	514.6	$92.5\pm 11.2$	0.71	1.057	25.52	$93.0\pm 11.3\pm 8.5$
4.236	530.3	$61.3\pm 9.9$	0.80	1.056	25.92	$52.2\pm 8.4\pm 4.7$
4.244	538.1	$21.9\pm 8.0$	0.86	1.055	25.51	$17.4\pm 6.4\pm 2.5$
4.267	531.1	$12.7\pm 9.1$	1.44	1.053	22.10
4.278	175.7	${0.0}_{-0.0}^{+3.0}$	2.68	1.053	17.11

Table 2. Table 2: Relative systematic uncertainties (in % percent \% ) from the different sources. Sources marked with an asterisk have common relative systematic uncertainties for the different center-of-mass energies. Dashes mean that the results are not applicable.

Source / $\sqrt{s}$ (GeV)	4.178	4.189	4.199	4.209	4.219	4.236	4.244	4.267	4.278
Luminosity∗	1.0	1.0	1.0	1.0	1.0	1.0	1.0	1.0	1.0
Tracking efficiency∗	4.0	4.0	4.0	4.0	4.0	4.0	4.0	4.0	4.0
Photon detection∗	2.0	2.0	2.0	2.0	2.0	2.0	2.0	2.0	2.0
Kinematic fit	0.6	0.8	0.6	0.6	0.5	0.6	0.8	2.4	3.8
Radiative correction	14.0	11.2	4.5	3.6	0.2	0.9	1.8	19.8	45.9
Fit range	$-$	9.3	2.3	2.9	2.5	2.0	2.3	8.7	$-$
Signal shape	$-$	1.9	5.1	9.4	0.2	1.5	11.4	1.6	$-$
Background shape	$-$	1.9	3.7	6.5	2.8	1.5	2.3	6.3	$-$
Angular distribution	0.1	0.1	0.3	0.5	0.5	0.8	0.6	0.7	0.9
Branching fraction∗	6.9	6.9	6.9	6.9	6.9	6.9	6.9	6.9	6.9
Sum	16.3	17.0	11.6	14.9	9.2	8.9	14.6	24.2	46.8

Table 3. Table 3: Summary of systematic uncertainties on the resonant parameters. The units for Γ e ​ e ​ ℬ ​ ( ω ​ χ c ​ 0 ) subscript Γ 𝑒 𝑒 ℬ 𝜔 subscript 𝜒 𝑐 0 \Gamma_{ee}\mathcal{B}(\omega\chi_{c0}) , M 𝑀 M , and Γ Γ \Gamma are eV, MeV/ c 2 superscript 𝑐 2 c^{2} , and MeV, respectively.

	${\Gamma }_{ee}ℬ\left(\omega {\chi }_{c0}\right)$	$M$	$\Gamma$
Absolute beam energy	0.1	0.9	0.2
Resonance parametrization	0.1	3.9	1.1
Cross section measurement	0.2	0.3	1.1
Sum	0.3	4.0	1.6

Equations

$$m\cdot (1-(\frac{m}{m_0}{)}^2{)}^p\cdot \text{exp}(k(1-(\frac{m}{m_0}{)}^2))\cdot \theta (m-m_0),$$

$${\sigma }^{\mathrm{B}}(e^{+}{e}^{-}\to \omega {\chi }_{c0})=\frac{N^{\mathrm{sig}}}{\mathcal{L}(1+\delta (s))\frac{1}{|1-\Pi {|}^2}\mathcal{B}ϵ},$$

$$\sigma (\sqrt{s})=\frac{12\pi {\Gamma }_{ee}\mathcal{B}(\omega {\chi }_{c0})\Gamma }{(s-M^2{)}^2+M^2{\Gamma }^2}\times \frac{\Phi (\sqrt{s})}{\Phi (M)},$$

Full text

I I. INTRODUCTION

$Y(4260)$ is the first charmonium-like $Y$ state, which was observed in the process $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}J/\psi$ by the BABAR experiment using an initial-state-radiation (ISR) technique Y4260-babar . This observation was immediately confirmed by the CLEO Y4260-cleo and Belle experiments Y4260-belle in the same process. $Y(4360)$ and $Y(4660)$ were also observed in $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}\psi(3686)$ Y4360-babar ; Y4360-belle . The observation of these $Y$ states has stimulated substantial theoretical discussions on their nature Y-theory1 ; Y-theory2 . These $Y$ states do not fit in the conventional charmonium spectroscopy, so they are good candidates for exotic states, such as hybrid states, tetraquark states and molecule states Y-theory3 . BESIII recently investigated the process $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ using data collected at $\sqrt{s}$ = 4.23 and 4.26 GeV combined with smaller data samples at nearby energies omegachic . An enhancement was found in the cross section around $\sqrt{s}$ = 4.22 GeV, referred to as the $Y(4220)$. Resonance signals were not observed in a subsequent study above 4.4 GeV omegachic2 . Various models the1 ; the2 ; the3 ; the4 ; the5 ; the6 ; the7 are proposed to explain the observed line-shape. Possible scenarios include a missing $\psi(4S)$ state the3 , a contribution from the $\psi(4160)$ state the4 , a tetraquark state the5 , or a molecule state the6 ; the7 . Intriguingly, similar and possibly related structures are also observed in the same energy region for other processes, such as $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}J/\psi$ pipijpsi-bes , $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}\psi(3686)$ pipipsip-bes , $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}h_{c}$ pipihc-bes , and $e^{+}e^{-}\rightarrow\pi^{+}D^{0}D^{*-}+c.c.$ piDDstar .

In this paper, we report a study of the $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ reaction based on the most recent $e^{+}e^{-}$ annihilation data collected with the BESIII detector besiii at nine energy points in the range 4.178 $\leqslant\sqrt{s}\leqslant$ 4.278 GeV, with a total integrated luminosity of about 7 fb*-1*. The $\chi_{c0}$ state is detected via $\chi_{c0}\rightarrow\pi^{+}\pi^{-}/K^{+}K^{-}$, and the $\omega$ is reconstructed via the $\omega\rightarrow\pi^{+}\pi^{-}\pi^{0}$ decay.

II II. BESIII DETECTOR AND MONTE CARLO SIMULATION

The BESIII detector is a magnetic spectrometer besiii located at the Beijing Electron Positron Collider (BEPCII) bepcii . The cylindrical core of the BESIII detector consists of a helium-based multilayer drift chamber (MDC), a plastic scintillator time-of-flight system (TOF), and a CsI(Tl) electromagnetic calorimeter (EMC), which are all enclosed in a superconducting solenoidal magnet providing a 1.0 T magnetic field. The solenoid is supported by an octagonal flux-return yoke with resistive plate counter muon identifier modules interleaved with steel. The acceptance of charged particles and photons is 93% over the $4\pi$ solid angle. The charged-particle momentum resolution at $1~{}{\rm GeV}/c$ is $0.5\%$, and the $\textrm{d}E/\textrm{d}x$ resolution is $6\%$ for the electrons from Bhabha scattering. The EMC measures photon energies with a resolution of $2.5\%$ ($5\%$) at $1$ GeV in the barrel (end cap) region. The time resolution of the TOF barrel part is 68 ps, while that of the end cap part is 110 ps. The end cap TOF system was upgraded in 2015 with multi-gap resistive plate chamber technology, providing a time resolution of 60 ps etof .

Simulated data samples produced with the geant4-based geant4 Monte Carlo (MC) package, which includes the geometric description of the BESIII detector and the detector response, are used to determine the detection efficiency and to estimate the background contributions. The simulation models the beam energy spread and initial state radiation (ISR) in the $e^{+}e^{-}$ annihilations using the generator kkmc KKMC . For the signal we use a MC sample of the $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ process generated according to the measured angular distribution, which is introduced in section VI. The inclusive MC samples consist of the production of open charm processes, the ISR production of vector charmonium(-like) states, and the continuum processes incorporated in kkmc KKMC . The known decay modes are modelled with evtgen ref:evtgen using branching fractions taken from the Particle Data Group (PDG) pdg , and the remaining unknown decays from the charmonium states are generated with lundcharm ref:lundcharm . Final state radiation (FSR) effects from charged final state particles are incorporated $via$ the photos package photos .

III III. EVENT SELECTION

For each charged track, the distance of closest approach to the interaction point (IP) is required to be within $\pm 10$ cm in the beam direction and within 1 cm in the plane perpendicular to the beam direction. The polar angle ($\theta$) of the tracks must be within the fiducial volume of the MDC $(|\cos\theta|<0.93)$. Photons are reconstructed from isolated showers in the EMC, which are at least $10^{\circ}$ away from the nearest charged track. The photon energy is required to be at least 25 MeV in the barrel region $(|\cos\theta|<0.8)$ or 50 MeV in the end cap region $(0.86<|\cos\theta|<0.92)$. To suppress electronic noise and energy depositions unrelated to the event, the time after the collision at which the photon is recorded in the EMC is required to satisfy $0\leq t\leq 700$ ns.

Since the final states of the $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ signal are $\pi^{0}\pi^{+}\pi^{-}\pi^{+}\pi^{-}$ or $\pi^{0}\pi^{+}\pi^{-}K^{+}K^{-}$, candidate events must have four tracks with zero net charge and at least two photons. The tracks with a momentum larger than 1 GeV/$c$ are identified as $\pi/K$ from the decay of the $\chi_{c0}$, whereas lower momentum tracks are considered as pions from $\omega$ decays. Since the tracks from $\omega$ and $\chi_{c0}$ can be separated clearly according to the momentum, the mis-identification rate is negligible. To reduce the background contributions and to improve the mass resolution, a 5C-kinematic fit is performed to both constrain the total four momentum of the final state particles to the total initial four momentum of the colliding beams and to constrain the invariant mass of the two photons from the decay of the $\pi^{0}$ to its nominal mass pdg . If there is more than one candidate in an event, the average multiplicity for signal is 1.09, the one with the smallest $\chi^{2}_{\text{5C}}$ of the kinematic fit is selected. The two track candidates of the decay of the $\chi_{c0}$ are considered to be either a $\pi^{+}\pi^{-}$ or a $K^{+}K^{-}$ pair depending on the $\chi^{2}$ of the 5C-kinematic fit. If $\chi^{2}_{\text{5C}}(\pi^{+}\pi^{-})<\chi^{2}_{\text{5C}}(K^{+}K^{-})$, the two tracks are identified as a $\pi^{+}\pi^{-}$ pair, otherwise, as a $K^{+}K^{-}$ pair. The $\chi^{2}_{\text{5C}}$ of the candidate events is required to be less than 100.

IV IV. BORN CROSS SECTION MEASUREMENT

The correlation between the $\pi^{+}\pi^{-}\pi^{0}$ invariant mass, $M(\pi^{+}\pi^{-}\pi^{0})$, and the $\pi^{+}\pi^{-}/K^{+}K^{-}$ mass, $M(\pi^{+}\pi^{-}/K^{+}K^{-})$, is shown in the top panel in Fig. 1 for data taken at $\sqrt{s}=4.219$ GeV. A high density area can be observed that originates from the $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ channel. The mass range [0.75, 0.81] GeV/$c^{2}$ in $M(\pi^{+}\pi^{-}\pi^{0})$ is defined as the $\omega$ signal region and is indicated by horizontal dashed lines. A sideband in the range [0.60, 0.72] GeV$/c^{2}$ is used to study the non-resonant background. The $\chi_{c0}$ signal region is indicated by the vertical dashed lines and is defined as [3.38, 3.45] GeV/$c^{2}$. The bottom panel of Fig. 1 shows the distribution of $M(\pi^{+}\pi^{-}/K^{+}K^{-})$ for data in the $\omega$ signal region. The shaded (green) histogram corresponds to normalized events in the $\omega$ sideband region.

An unbinned maximum likelihood fit is performed to obtain the signal yields. In the fit, we use the MC-determined shape to describe the $\chi_{c0}$ signal. The background is described with a generalized ARGUS function argus

[TABLE]

where $m_{0}$ is fixed to $(\sqrt{s}-0.75$ GeV), with 0.75 GeV being the lower limit of $M(\pi^{+}\pi^{-}\pi^{0})$, and $p$, $k$ are free parameters. The red solid curve in the bottom panel of Fig. 1 shows the fit result. The data taken at the other center-of-mass energies are analyzed using the same method and the fit results are summarized in Table 1.

The Born cross section is calculated with

[TABLE]

where $N^{\text{sig}}$ is the number of signal events, $\mathcal{L}$ is the integrated luminosity obtained using the same method in Ref. luminosity , 1 + $\delta(s)$ is the radiative correction factor obtained from a Quantum Electrodynamics (QED) calculation QED ; KKMC using the obtained preliminary cross section as input and iterating it until the results converge, $\frac{1}{|1-\Pi|^{2}}$ is the correction factor for vacuum polarization vacuum , $\mathcal{B}$ is the product of branching fractions $\mathcal{B}(\chi_{c0}\rightarrow\pi^{+}\pi^{-}/K^{+}K^{-})\times\mathcal{B}(\omega\rightarrow\pi^{+}\pi^{-}\pi^{0})\times\mathcal{B}(\pi^{0}\rightarrow\gamma\gamma)$, and $\epsilon$ is the event selection efficiency. The Born cross sections (or upper limits at 90% C.L.) at each energy point for $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ are listed in Table 1.

The systematic uncertainty of the Born cross section measurement originates mainly from the luminosity determination, the tracking efficiency, photon detection efficiency, kinematic fit, radiative correction, fit range, signal and background shapes, angular distribution, and the branching fractions for $\mathcal{B}(\chi_{c0}\rightarrow\pi^{+}\pi^{-}/K^{+}K^{-})\times\mathcal{B}(\omega\rightarrow\pi^{+}\pi^{-}\pi^{0})\times\mathcal{B}(\pi^{0}\rightarrow\gamma\gamma)$.

The luminosity is measured with a precision of about $1.0\%$ using the well-known Bhabha scattering process luminosity . The uncertainty in the tracking efficiency is obtained as $1.0\%$ per track using the process $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}K^{+}K^{-}$ omegachic2 . The uncertainty in photon reconstruction is $1.0\%$ per photon, obtained by studying the $J/\psi\rightarrow\rho^{0}\pi^{0}$ decay photon .

The systematic uncertainty due to the kinematic fit is estimated by correcting the helix parameters of charged tracks according to the method described in Ref. helix . The difference between detection efficiencies obtained from MC samples with and without correction is taken as the uncertainty.

The line-shape of the $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ cross section will affect the radiative correction factor and the efficiency. In the nominal results, we use a bifurcated Gaussian function as the line-shape to describe the cross section. The shape is used as input and is iterated until the results converge. To estimate the uncertainty from the radiative correction, we change the line-shape to the Breit-Wigner (BW) function of the $Y(4220)$ omegachic2 . The difference between the results is taken as a systematic uncertainty.

The uncertainty from the fit range is obtained by varying the limits of the fit range by $\pm$0.01 GeV/$c^{2}$. We take the largest difference of the corresponding cross section measurement with respect to the nominal one as the systematic uncertainty. We use the MC-determined shape convolved with a Gaussian function to fit the data as input to get the uncertainty of the signal shape. The difference in the results with respect to the nominal one is taken as the systematic uncertainty. To estimate the systematic uncertainty caused by the background shape, we vary $m_{0}$ by $\pm 0.01$ GeV/$c^{2}$ in the ARGUS function, and take the largest difference in the results as the uncertainty.

The measured angular distribution is used as a model to generate signal events in the MC simulations. The detection efficiency of the $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ reaction will depend upon its angular distribution. We obtained an angular distribution parameter, defined in section VI, of $\alpha=-0.30\pm 0.18(\text{stat.})\pm 0.05(\text{sys.})$. The systematic uncertainty of the efficiency due to uncertainties in the angular distribution is estimated by varying the $\alpha$ value by one standard deviation, the total uncertainty on $\alpha$.

The uncertainty in the product of the branching fractions $\mathcal{B}(\chi_{c0}\rightarrow\pi^{+}\pi^{-}/K^{+}K^{-})\times\mathcal{B}(\omega\rightarrow\pi^{+}\pi^{-}\pi^{0})\times\mathcal{B}(\pi^{0}\rightarrow\gamma\gamma)$ is taken from the uncertainties quoted by the PDG pdg .

Table 2 summarizes all the systematic uncertainties related to the cross section measurements of the $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ process for each center-of-mass energy. The overall systematic uncertainties are obtained by adding all the sources of systematic uncertainties in quadrature assuming they are uncorrelated.

V V. resonant parameter measurement

Figure 2 shows the dressed cross sections ($\sigma=\frac{\sigma^{\rm B}}{|1-\Pi|^{2}}$) for the $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ reaction as a function of center-of-mass energy. The black square points are taken from Refs. omegachic ; omegachic2 , and the blue circular points are from this work. We observe an enhancement in the cross section around 4.22 GeV. By assuming that the $\omega\chi_{c0}$ signals all come from a single resonance, which we label as the $Y(4220)$, with mass $M$ and width $\Gamma$, we fit the cross section data with the following formula convolved with a Gaussian function for the energy spread:

[TABLE]

where $\Phi(\sqrt{s})$ is the two-body phase space factor and $\Gamma_{ee}$ is the electronic width. The fit to all the data in Fig. 2 gives $\Gamma_{ee}\mathcal{B}(\omega\chi_{c0})=(2.5\pm 0.2)$ eV, $M=(4218.5\pm 1.6)$ MeV/$c^{2}$, $\Gamma=(28.2\pm 3.9)$ MeV, where the uncertainties are statistical only. In the fit, the cross sections’ statistical uncertainties are used only. The goodness of fit is $\chi^{2}/\text{ndf}=29/19$, where ndf is the number of degrees of freedom.

The systematic uncertainties on the resonant parameters mainly arise from uncertainties in the absolute beam energy, the parametrization of the BW function, and the cross section measurement. Because the energy spread effect has been considered in the fit, we ignore the systematic uncertainty from energy spread.

Since the uncertainty of the beam energy is about 0.8 MeV, which is obtained using the same method in Ref. ecms , the uncertainty of the resonant parameters caused by the beam energy is estimated by varying $\sqrt{s}$ within 0.8 MeV.

The cross section has been fitted with a BW function having the energy-dependent width $\Gamma=\Gamma^{0}\frac{\Phi(\sqrt{s})}{\Phi(M)}$ in the denominator, where $\Gamma^{0}$ is the nominal width of the resonance, to estimate the uncertainty from parametrization of the BW function. The difference between this fit result and the nominal result is taken as the uncertainty from the parametrization of the BW function.

The systematic uncertainty of the cross section measurement will affect the resonant parameters in the fit and can be divided into two parts. One part comes from the uncorrelated uncertainty among the different center-of-mass energies, and the other part is a common uncertainty. The first part has been considered by including the systematic uncertainty of the cross section in the fit. The difference between the parameters obtained in this fit to those from the nominal fit is taken as the uncertainty. We vary the cross section within the systematic uncertainty coherently for the second part and take the difference between this fit result and the nominal result as the uncertainty. We add the two parts in quadrature assuming they are uncorrelated.

Table 3 summarizes all the systematic uncertainties of the resonant parameters. The total systematic uncertainty is obtained by summing all the sources of systematic uncertainties in quadrature by assuming they are uncorrelated.

VI VI. angular distribution measurement

Both $S$ and $D$-wave contributions are possible in the process $Y(4220)\rightarrow\omega\chi_{c0}$. A measurement of their strengths can be helpful to extract information about the underlying dynamics of the decay process. We therefore performed an angular analysis angular of the relatively high-statistics data samples taken at $\sqrt{s}$ = 4.219, 4.226, and 4.236 GeV (selection of $\sqrt{s}$ = 4.226 GeV data was reported in Ref. omegachic ). The helicity angle, $\theta_{\omega}$, defined by the scattering angle of the $\omega$ with respect to the electron beam in the $e^{+}e^{-}$ center-of-mass frame was reconstructed for each event. Figure 3 shows the bin-by-bin efficiency-corrected events as a function of $\text{cos}\theta_{\omega}$ for the three center-of-mass energies. The signal yield in each of the 10 bins is determined with the same method as that in the cross section measurement, and the detection efficiency in each bin is determined with the signal MC sample. We performed a simultaneous fit using the function $1+\alpha$cos${}^{2}\theta_{\omega}$ with a least-square method, assuming $\alpha$ is common to the three energy points. The red line in Fig. 3 shows the best fit result with $\alpha=-0.30\pm 0.18\pm 0.05$, where the first uncertainty is statistical and the second systematic. The goodness of the fit is $\chi^{2}/$ndf$=31/26$. The fit result indicates evidence for a combination of $S$ and $D-$wave contributions in the $Y(4220)\rightarrow\omega\chi_{c0}$ process, although the statistical significance of this conclusion is only $2\sigma$ compared with a pure $S-$wave contribution of $\alpha=0$.

The systematic uncertainty of $\alpha$ has been estimated by varying the fit range (0.02), the signal (0.01) and background shapes (0.03), and the radiative correction factor (0.03). The uncertainties are indicated in brackets and are determined with the same method described earlier for the cross section measurements. In addition, we estimate an additional source of systematic uncertainty by varying the number of bins. For this, we change the number of bins from 10 to 8, and repeat the process. The difference in $\alpha$ is found to be 0.01. The overall systematic uncertainty (0.05) is obtained by summing all the items of systematic uncertainties in quadrature by assuming they are uncorrelated.

VII VII. Summary

The process $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ has been studied using 9 data samples collected at center-of-mass energies from $\sqrt{s}=$ 4.178 to 4.278 GeV. The $\sqrt{s}$-dependence of the cross section has been measured and the results are listed in Table 1 and are shown in Fig. 2. A clear enhancement is seen around $\sqrt{s}=$ 4.22 GeV which confirms, and statistically improves upon, an earlier observation omegachic . By fitting the $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ cross section with a single resonance, the mass and width for the structure are determined to be $M=(4218.5\pm 1.6(\text{stat.})\pm 4.0(\text{sys.}))$ MeV/$c^{2}$ and $\Gamma=(28.2\pm 3.9(\text{stat.})\pm 1.6(\text{sys.}))$ MeV. The obtained resonance parameters are not compatible with the vector charmonium state $\psi(4160)$, ruling out its possible contribution to the structure the4 . Moreover, we studied the angular distribution of the process $Y(4220)\rightarrow\omega\chi_{c0}$. We measured $\alpha=-0.30\pm 0.18\pm 0.05$, which indicates a combination of $S$ and $D-$wave contributions in the decay.

Figure 4 shows the measured mass and width of the $Y(4220)$ from the different processes. The masses are consistent with each other, while the widths are not. The widths from the processes $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}h_{c}$, $\pi^{+}\pi^{-}\psi(3686)$, and $\pi^{+}D^{0}D^{*-}+c.c.$ are larger than those from the processes $e^{+}e^{-}\rightarrow\omega\chi_{c0}$ and $\pi^{+}\pi^{-}J/\psi$. From these inconsistencies in the width, we cannot draw a conclusion on whether the structure observed in these processes is the same state or whether the inconsistencies are caused by the BW parameterization. Further experimental studies with higher statistics are needed to draw a more reliable conclusion on the nature of this structure.

VIII ACKNOWLEDGMENTS

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. 11625523, 11635010, 11735014; National Natural Science Foundation of China (NSFC) under Contract No. 11835012; National Natural Science Foundation of China (NSFC) under Contract No. 11847028; National Natural Science Foundation of China (NSFC) under Contract No. 11575198; the Chinese Academy of Sciences (CAS) Large-Scale Scientific Facility Program; Joint Large-Scale Scientific Facility Funds of the NSFC and CAS under Contracts Nos. U1532257, U1532258, U1732263, U1832207; 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; Foundation of Henan Educational Committee (No. 19A140015); Nanhu Scholars Program for Young Scholars of Xinyang Normal University; 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 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.