Observation of X(2370) and search for X(2120) in - BESIII Collaboration

Using a sample of events collected with the BESIII detector, we perform a study of . X (2370) is observed in the invariant-mass distribution with a statistical significance of . Its resonance parameters are measured to be and . The product branching fractions for and are determined to be and , respectively. No evident signal for X (2120) is observed in the invariant-mass distribution. The upper limits for the product branching fractions of and are determined to be and at the 90% confidence level, respectively.


Introduction
Quantum chromodynamics (QCD), a non-Abelian gauge field theory, predicts the existence of new types of hadrons with explicit gluonic degrees of freedom (e.g., glueballs, hybrids) [1 -3]. The search for glueballs is an important field of research in hadron physics. It is, however, challenging since possible mixing of pure glueball states with nearby qq nonet mesons makes the identification of glueballs difficult in both experiment and theory. Lattice QCD (LQCD) predicts the lowest-lying glueballs which are scalar (mass 1.5-1.7 GeV/c 2 ), tensor (mass 2.3-2.4 GeV/c 2 ), and pseudoscalar (mass 2.3-2.6 GeV/c 2 ) [4]. Radiative J/ψ decay is a gluon-rich process and it is therefore regarded as one of the most promising hunting grounds for glueballs [5,6]. Recent LQCD calculations predict that the partial width of J/ψ radiatively decaying into the pure gauge pseudoscalar glueball is 0.0215(74) keV which corresponds to a branching ratio 2.31(80) × 10 −4 [7]. Recently, three states, X (1835), X (2120) and X (2370), were observed in the BESIII experiment in the π + π − η invariant-mass distribution through the decay of J/ψ → γ π + π − η with statistical significances larger than 20σ , 7.2σ and 6.4σ , respec-tively [8]. The measured mass of X (2370) is consistent with the pseudoscalar glueball candidate predicted by LQCD calculations [4]. In the case of a pseudoscalar glueball, the branching fractions of X (2370) decaying into K K η and ππη are predicted to be 0.011 and 0.090 [9], respectively, in accordance with calculations that are based upon the chiral effective Lagrangian. Study on the decays to KK η of the glueball candidate X states is helpful to identify their natures.

Detector and Monte Carlo simulations
The BESIII detector is a magnetic spectrometer [11] located at the Beijing Electron Positron Collider II (BEPCII) [12]. 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 (0.9 T in 2012) 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 4π solid angle. The charged-particle momentum resolution at 1 GeV/c is 0.5%, and the d E/dx 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.
Simulated samples produced with the geant4-based [13] 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 backgrounds. The simulation includes the beam energy spread and initial-state radiation (ISR) in the e + e − annihilations modeled with the generator kkmc [14,15]. The inclusive MC sample consists of the production of the J/ψ resonance, and the continuum processes incorporated in kkmc [14,15]. The known decay modes are modeled with evtgen [16,17] using branching fractions taken from the Particle Data Group [18], and the remaining unknown decays from the charmonium states are generated with lundcharm [19,20]. The final-state radiations (FSR) from charged final-state particles are incorporated with the photos package [21]. Background is stud-ied using a sample of 1.2 × 10 9 simulated J/ψ events. Phase-space (PHSP) MC samples of J/ψ → γ K + K − η and J/ψ → γ K 0 S K 0 S η are generated to describe the nonresonant contribution. To estimate the selection efficiency and to optimize the selection criteria, signal MC events are generated for J/ψ → γ X (2120)/ X (2370) → γ K + K − η and J/ψ → γ X (2120)/ X (2370) → γ K 0 S K 0 S η channel, respectively. The polar angle of the photon in the J/ψ centerof-mass system, θ γ , follows a 1 + cos 2 θ γ function. For the process of η → γρ 0 , ρ 0 → π + π − , a generator taking into account both the ρ-ω interference and the box anomaly is used [22]. The analysis is performed in the framework of the BESIII offline software system (BOSS) [23] incorporating the detector calibration, event reconstruction and data storage.

Event selection
Charged-particle tracks in the polar angle range | cos θ | < 0.93 are reconstructed from hits in the MDC. Tracks (excluding those from K 0 S decays) are selected that extrapolated to be within 10 cm from the interaction point in the beam direction and 1 cm in the plane perpendicular to the beam. The combined information from energy-loss (dE/dx) measurements in the MDC and time in the TOF is used to obtain confidence levels for particle identification (PID) for π , K and p hypotheses. For J/ψ → γ K + K − η decay, each track is assigned to the particle type corresponding to the highest confidence level; candidate events are required to have four charged tracks with zero net charge and with two opposite charged tracks identified as kaons and the other two identified as pions. For the J/ψ → γ K 0 S K 0 S η decay, each track is assumed to be a pion and no PID restrictions are applied; candidate events are required to have six charged tracks with zero net charge. K 0 S candidates are reconstructed from a secondary vertex fit to all π + π − pairs, and each K 0 S candidate is required to satisfy |M π + π − − m K 0 candidates are used as input for the subsequent kinematic fit. Photon candidates are required to have an energy deposition above 25 MeV in the barrel region (| cos θ | < 0.80) and 50 MeV in the end cap (0.86 < | cos θ | < 0.92). To exclude showers from charged tracks, the angle between the shower position and the charged tracks extrapolated to the EMC must be greater than 5 • . A timing requirement in the EMC is used to suppress electronic noise and energy deposits unrelated to the event. At least two (three) photons are required for the η → γρ 0 (η → π + π − η) mode.
For the J/ψ → γ K + K − η (η → γρ 0 ) channel, a fourconstraint (4C) kinematic fit is performed to the hypoth- The dots with error bars correspond to data and the histograms are the results of PHSP MC simulations (arbitrary normalization) esis of J/ψ → γ γ K + K − π + π − by requiring the total energy and each momentum component to be conserved. For events with more than two photon candidates, the combination with the minimum χ 2 4C is selected, and χ 2 4C < 25 is required. Events with |M γ γ − m π 0 | < 30 MeV/c 2 or |M γ γ − m η | < 30 MeV/c 2 are rejected to suppress background containing π 0 or η, where m π 0 and m η are the nominal masses of π 0 and η [18]. A clear η signal is observed in the invariant-mass distribution of γ π + π − (M γ π + π − ), as shown in Fig. 1(a). Candidates of ρ and η are reconstructed from the π + π − and γ π + π − combinations with 0.55 GeV/c 2 < M π + π − < 0.85 GeV/c 2 and |M γ π + π − − m η | < 20 MeV/c 2 , where m η is the nominal mass of η [18], respectively. If there is more than one combination satisfying the selection criteria, the combination with M γ π + π − closest to m η is selected. After applying the above requirements, we obtain the invariant-mass distribution of K + K − η (M K + K − η ) as shown in Fig. 1(b). The peak around 2.98 GeV/c 2 is contributed from the decay of J/ψ → γ η c (η c → K + K − η ), while the peak around the right threshold is mainly from the background events of To reduce background and to improve the mass resolution of the J/ψ → γ K + K − η (η → π + π − η) channel, a five-constraint (5C) kinematic fit is performed whereby the total four momenta of the final-state particles are constrained by the total initial four momentum of the col-liding beams and the invariant mass of the two photons from the decay of η is constrained by its nominal mass. If there are more than three photon candidates, the combination with the minimum χ 2 5C is retained, and χ 2 5C < 45 is required. To suppress background from π 0 → γ γ , |M γ γ − m π 0 | > 30 MeV/c 2 is required for all photon pairs. The η candidates are formed from the π + π − η combination satisfying |M π + π − η − m η | < 15MeV/c 2 , where M π + π − η is the invariant mass of π + π − η, as shown in Fig. 1c. After applying the mass restrictions, we obtain the invariantmass distribution of K + K − η (η → π + π − η) as shown in Fig. 1d.
candidates are subjected to a 4C kinematic fit. For events with more than two photons or two K 0 S candidates, the combination with the smallest χ 2 4C is retained, and χ 2 4C < 45 is required. To suppress background events containing a π 0 or η, events with |M γ γ − m π 0 | < 30 MeV/c 2 or |M γ γ − m η | < 30 GeV/c 2 are rejected. The π + π − invariant mass is required to be in the ρ mass region, 0.55 GeV/c 2 < M π + π − < 0.85 GeV/c 2 , and |M γ π + π − − m η | < 20 MeV/c 2 is applied to select the η signal. If more than one combination of γ π + π − is obtained, the combination with M γ π + π − closest to m η is selected as shown in Fig. 2a. After applying the above requirements, we obtain the K 0 Fig. 2b.
Candidate events of the J/ψ → γ K 0 S K 0 S η (η → π + π − η) channel are subjected to a 5C kinematic fit, which is similar to that for the J/ψ → γ K + K − η (η → π + π − η) mode. If there are more than three photons or more than two K 0 S candidates, only the combination with the minimum χ 2 5C is selected and χ 2 5C <50 is required. To reduce the combinatorial background from π 0 → γ γ events, |M γ γ − m π 0 | > 30 MeV/c 2 is required for all photon pairs. For selecting the η signal, the π + π − η combination satisfying |M π + π − η − m η | < 15 MeV/c 2 is required, as shown in Fig. 2c. After applying the above selection criteria, we obtain the invariant-mass distribu-
A structure near 2.34 GeV/c 2 is observed in the invariantmass distribution of K + K − η and K 0 S K 0 S η . We performed a simultaneous unbinned maximum-likelihood fit to the K + K − η and K 0 S K 0 S η invariant-mass distributions between 2.0 and 2.7 GeV/c 2 , as shown in Fig 3. The signal is represented by an efficiency-weighted non-relativistic Breit-Wigner (BW) function convolved with a double Gaussian function to account for the mass resolution. The mass and width of the BW function are left free in the fit, while the parameters of the double Gaussian function are fixed on the results obtained from the fit of signal MC samples generated with zero width. The non-η background events are described with η sideband data and the yields from these sources are fixed; the J/ψ → K * + K − η + c.c. contributions in the J/ψ → γ K + K − η decay channel are studied as discussed above and the shapes and the yields are fixed in the fit; the contribution from the nonresonant γ KK η production is described by the shape from the PHSP MC sample of J/ψ → γ KK η and its absolute yield is set as a free parameter in the fit; the remaining background is described by a second order Chebychev polynomial function and its parameters are left to be free. In the simultaneous fit, the resonance parameters are free parameters and constrained to be the same for all four channels. The signal ratio for the two η decay modes is fixed with a factor calculated by their branching fractions and efficiencies. The signal ratio between J/ψ → γ X (2370) → γ K + K − η and J/ψ → γ X (2370) → γ K 0 S K 0 S η is a free parameter in the fit. The obtained mass, width and the number of signal events for X (2370) are listed in Table 1. A variety of fits with different fit ranges, η sideband regions and background shapes are performed; after considering the systematic uncertainties like quantum number of X(2370) and the presence of X(2120)*, the smallest statistical significance among these fits is found to be 8.3σ . With the detection efficiencies listed in Table 2, the product branching fractions for J/ψ → γ X (2370), X (2370) → K + K − η and J/ψ → γ X (2370), X (2370) → K 0 S K 0 S η are determined to be (1.79 ± 0.23) × 10 −5 and (1.18 ± 0.32) × 10 −5 , respectively, where the uncertainties are statistical only.
No obvious signal of X (2120) is found in the KK η invariant-mass distribution. We performed a simultaneous unbinned maximum-likelihood fit to the KK η invariantmass distribution in the range of [2.0, 2.7] GeV/c 2 . The signal, X (2120), is described with an efficiency-weighted BW function convolved with a double Gaussian function. The mass and width of the BW function are fixed to previously published BESIII results [8]. The backgrounds are modeled with the same components as used in the fit of X (2370) as mentioned above. The contribution from X (2370) is included in the fit and its mass, width and the numbers of events are set free. The distribution of normalized likelihood values for a series of input signal event yields is taken as the probability density function (PDF) for the expected number of events. The number of events at 90% of the integral of the PDF from zero to the given number of events is defined as the upper limit, N U L , at the 90% confidence level (CL). We repeated this procedure with different signal shape parameters of X (2120) (by varying the values of mass and width with 1σ of the uncertainties cited from [8]), fit ranges, η sideband regions and background shapes, and the maximum upper limit among these cases is selected. The statistical significance of X (2120) is determined to be 2.2σ . To calculate   Fig. 3 The fit result for X (2370) in the invariant-mass distribution of KK η for the decays: and B(J/ψ → γ X (2120) → γ K 0 S K 0 S η ) < 6.15 × 10 −6 , respectively.

Systematic uncertainties
Several sources of systematic uncertainties are considered for the determination of the mass and width of X (2370) and the product branching fractions. These include the efficiency differences between data and MC simulation in the MDC tracking, PID, the photon detection, K 0 S reconstruction, the kinematic fitting, and the mass-window requirements of π 0 , η, ρ and η . Furthermore, uncertainties associated with the fit ranges, the background shapes, the sideband regions, the signal shape parameters of X (2120), intermediate resonance decay branching fractions and the total number of J/ψ events are considered.

Efficiency estimation
The MDC tracking efficiencies of charged pions and kaons are investigated using nearly background-free (clean) con- Table 1 Fit results for the structure around 2.34 GeV/c 2 and 2.12 GeV/c 2 . The superscripts a and b represent the decay modes of X → K + K − η and X → K 0 S K 0 S η , respectively. The uncertainties are statistical only trol samples of J/ψ → ppπ + π − and J/ψ → K 0 S K ± π ∓ [24,25], respectively. The difference in tracking efficiencies between data and MC is 1.0% for each charged pion and kaon. The photon detection efficiency is studied with a clean sample of J/ψ → ρ 0 π 0 [26], and the result shows that the difference of photon detection efficiencies between data and MC simulation is 1.0% for each photon. The systematic uncertainty from K 0 S reconstruction is determined from the control sam- Table 2 Summary of the MC detection efficiencies of the signal yields for the two η modes where the KK η invariant mass is constrained to the applied fitting range between 2.0 and 2.7 GeV/c 2 . The superscripts a and b represent the decay modes of X → K + K − η and X → K 0 S K 0 S η , respectively Decay modes ε η →γρ 0 (%) ε η →π + π − η (%) ples of J/ψ → K * ± K ∓ and J/ψ → φ K 0 S K ± π ∓ , which indicates that the efficiency difference between data and MC is less than 1.5% for each K 0 S . Therefore, 3.0% is taken as the systematic uncertainty for the two For the decay channel of J/ψ → γ K + K − η , the PID has been used to identify the kaons and pions. Using a clean sample of J/ψ → ppπ + π − , the PID efficiency of π + /π − has been studied, which indicates that the π + /π − PID efficiency for data agrees with MC simulation within 1%. The PID efficiency for the kaon is measured with a clean sample of J/ψ → K + K − η. The difference of the PID efficiency between data and MC is less than 1% for each kaon. Hence, in this analysis, four charged tracks are required to be identified as two pions and two kaons, and 4% is taken as the systematic uncertainty associated with the PID.
The systematic uncertainties associated with the kinematic fit are studied with the track helix parameter correction method, as described in Ref. [27]. The differences from those without corrections are taken as systematic uncertainties.
Due to the difference in the mass resolution between data and MC, uncertainties related to the ρ 0 and η mass-window requirements are investigated by smearing the MC simulation to improve the consistency between data and MC simulation. The differences in the detection efficiency before and after smearing are assigned as systematic uncertainties for the ρ 0 and η mass-window requirements. The uncertainties from the π 0 and η mass-window requirements are estimated by varying the mass windows of π 0 and η, and differences in the resulting branching fractions are assigned as the systematic uncertainties of this item.
Furthermore, we considered the effects arising from different quantum numbers of X (2120) and X (2370). We generated J/ψ → γ X (2120) and J/ψ → γ X (2370) decays following a sin 2 θ γ angular distribution. The resulting differences in efficiency from the nominal value are taken as systematic uncertainties.

Fit to the signal
To study the uncertainties from the fit range and η sideband region, the fits are repeated with different fit ranges and sideband regions, the largest differences among these signal yields are taken as systematic uncertainties, respectively. To estimate the uncertainties in the description of various background contributions, we performed alternative fits with third-order Chebychev polynomials modeling the background of the K + K − η and K 0 S K 0 S η channels. The maximum differences in signal yield from the nominal fit are taken as systematic uncertainties. The uncertainties from the background of J/ψ → K * + K − η + c.c. are estimated by absorbing this component into a Chebychev polynomial function, and the differences obtained by using the description with or without the background component of J/ψ → K * + K − η + c.c. are taken as systematic uncertainties. The impact of X (2120) is also considered as a systematic uncertainty in the study of X (2370). The difference between a fit with and without a X (2120) contribution is taken as a systematic uncertainty associated to this item. Table 3 Absolute systematic uncertainties of resonance parameters of mass (M, in MeV/c 2 ) and width ( , in MeV) for X (2370). The items with * are common uncertainties of both η decay modes  Table 4 Systematic uncertainties for determination of the branching fraction of J/ψ → γ X (2370) → γ KK η (in %). The items with * are common uncertainties of both η decay modes. I and II represent the decay modes of η → γρ 0 , ρ 0 → π + π − and η → π + π − η, η → γ γ , respectively

Others
Since no evident structures are observed in the invariantmass distributions of M(K η ), M(K η ) and M(KK ) for the events with a KK η invariant mass within the X (2370) mass region (2.2 GeV/c 2 < M KK η < 2.5 GeV/c 2 ), the systematic uncertainties of the reconstruction efficiency due to the possible intermediate states on the K η ,K η and KK mass spectra are ignored. The uncertainties on the intermediate decay branching fractions of η → γρ 0 → γ π + π − , η → π + π − η, η → γ γ and K 0 S → π + π − are taken from the world average values [18], which are 1.7%, 1.6%, 0.5% and 0.1%, respectively. The systematic uncertainty due to the number of J/ψ events is determined as 0.5% according to Ref. [10].
A summary of all the uncertainties is shown in Tables 3, 4 and 5. The total systematic uncertainties are obtained by adding all individual uncertainties in quadrature, assuming all sources to be independent.
X (2120) and X (2370) are studied via J/ψ → γ K + K − η and J/ψ → γ K 0 S K 0 S η with two η decay modes, respec- Table 5 Systematic uncertainties for the determination of the upper limit of the branching fraction of J/ψ → γ X (2120) → γ KK η (in %). The items with * are common uncertainties of both η decay modes. I and II represent the decay modes of η → γρ 0 , ρ 0 → π + π − and η → π + π − η, η → γ γ , respectively  tively. The measurements from the two η decay modes are, therefore, combined by considering the difference in uncertainties of these two measurements. The combined systematic uncertainties are calculated with the weighted least squares method [28] and the results are shown in Table 6.

Results and summary
Using a sample of 1.31 × 10 9 J/ψ events collected with the BESIII detector, the decays of J/ψ → γ K + K − η and J/ψ → γ K 0 S K 0 S η are investigated using the two η decay modes, η → γρ 0 (ρ 0 → π + π − ) and η → π + π − η(η → γ γ ). X (2370) is observed in the KK η invariant-mass distribution with a statistical significance of 8.3σ . The mass and width are determined to be M X (2370) = 2341.6 ± 6.5 (stat.) ± 5.7 (syst.) MeV/c 2 , X (2370) = 117 ± 10 (stat.) ± 8 (syst.) MeV, which are found to be consistent with those of X (2370) observed in the previous BESIII results [8]. The product branching fractions of B(J/ψ → γ X (2370) → γ K + K − η ) and B(J/ψ → γ X (2370) → γ K 0 S K 0 S η ) are measured to be (1.79 ± 0.23 (stat.) ± 0.65 (syst.)) × 10 −5 and (1.18 ± 0.32 (stat.) ± 0.39 (syst.)) × 10 −5 , respectively. No evident signal for X (2120) is observed in the KK η invariant-mass distribution. For a conservative estimate of the upper limits of the product branching fractions of J/ψ → γ X (2120) → K + K − η and J/ψ → γ X (2120) → K 0 S K 0 S η , the multiplicative uncertainties are considered by convolving the normalized likelihood function with a Gaussian function. The upper limits for product branching fractions at 90% C. L. are determined to be B(J/ψ → γ X (2120) → γ K + K − η ) < 1.49 × 10 −5 and B(J/ψ → γ X (2120) → γ K 0 S K 0 S η ) < 6.38 × 10 −6 . To understand the nature of X (2120) and X (2370), it is critical to measure their spin and parity and to search for them in more decay modes. A partial-wave analysis is needed to measure their masses and widths more precisely, and to determine their spin and parity. This might become possible in the future with the foreseen higher statistics of J/ψ data samples. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.