Study of the X ð 5568 Þ state with semileptonic decays of the B 0 s meson

V.M. Abazov, B. Abbott, B. S. Acharya, M. Adams, T. Adams, J. P. Agnew, G. D. Alexeev, G. Alkhazov, A. Alton, A. Askew, S. Atkins, K. Augsten, V. Aushev, Y. Aushev, C. Avila, F. Badaud, L. Bagby, B. Baldin, D. V. Bandurin, S. Banerjee, E. Barberis, P. Baringer, J. F. Bartlett, U. Bassler, V. Bazterra, A. Bean, M. Begalli, L. Bellantoni, S. B. Beri, G. Bernardi, R. Bernhard, I. Bertram, M. Besançon, R. Beuselinck, P. C. Bhat, S. Bhatia, V. Bhatnagar, G. Blazey, S. Blessing, K. Bloom, A. Boehnlein, D. Boline, E. E. Boos, G. Borissov, M. Borysova, A. Brandt, O. Brandt, M. Brochmann, R. Brock, A. Bross, D. Brown, X. B. Bu, M. Buehler, V. Buescher, V. Bunichev, S. Burdin, C. P. Buszello, E. Camacho-Pérez, B. C. K. Casey, H. Castilla-Valdez, S. Caughron, S. Chakrabarti, K. M. Chan, A. Chandra, E. Chapon, G. Chen, S. W. Cho, S. Choi, B. Choudhary, S. Cihangir, D. Claes, J. Clutter, M. Cooke, W. E. Cooper, M. Corcoran, F. Couderc, M.-C. Cousinou, J. Cuth, D. Cutts, A. Das, G. Davies, S. J. de Jong, E. De La Cruz-Burelo, F. Déliot, R. Demina, D. Denisov, S. P. Denisov, S. Desai, C. Deterre, K. DeVaughan, H. T. Diehl, M. Diesburg, P. F. Ding, A. Dominguez, A. Drutskoy, A. Dubey, L. V. Dudko, A. Duperrin, S. Dutt, M. Eads, D. Edmunds, J. Ellison, V. D. Elvira, Y. Enari, H. Evans, A. Evdokimov, V. N. Evdokimov, A. Fauré, L. Feng, T. Ferbel, F. Fiedler, F. Filthaut, W. Fisher, H. E. Fisk, M. Fortner, H. Fox, J. Franc, S. Fuess, P. H. Garbincius, A. Garcia-Bellido, J. A. García-González, V. Gavrilov, W. Geng, C. E. Gerber, Y. Gershtein, G. Ginther, O. Gogota, G. Golovanov, P. D. Grannis, S. Greder, H. Greenlee, G. Grenier, Ph. Gris, J.-F. Grivaz, A. Grohsjean, S. Grünendahl, M.W. Grünewald, T. Guillemin, G. Gutierrez, P. Gutierrez, J. Haley, L. Han, K. Harder, A. Harel, J. M. Hauptman, J. Hays, T. Head, T. Hebbeker, D. Hedin, H. Hegab, A. P. Heinson, U. Heintz, C. Hensel, I. Heredia-De La Cruz, K. Herner, G. Hesketh, M. D. Hildreth, R. Hirosky, T. Hoang, J. D. Hobbs, B. Hoeneisen, J. Hogan, M. Hohlfeld, J. L. Holzbauer, I. Howley, Z. Hubacek, V. Hynek, I. Iashvili, Y. Ilchenko, R. Illingworth, A. S. Ito, S. Jabeen, M. Jaffré, A. Jayasinghe, M. S. Jeong, R. Jesik, P. Jiang, K. Johns, E. Johnson, M. Johnson, A. Jonckheere, P. Jonsson, J. Joshi, A.W. Jung, A. Juste, E. Kajfasz, D. Karmanov, I. Katsanos, M. Kaur, R. Kehoe, S. Kermiche, N. Khalatyan, A. Khanov, A. Kharchilava, Y. N. Kharzheev, I. Kiselevich, J. M. Kohli, A. V. Kozelov, J. Kraus, A. Kumar, A. Kupco, T. Kurča, V. A. Kuzmin, S. Lammers, P. Lebrun, H. S. Lee, S. W. Lee, W.M. Lee, X. Lei, J. Lellouch, D. Li, H. Li, L. Li, Q. Z. Li, J. K. Lim, D. Lincoln, J. Linnemann, V. V. Lipaev, R. Lipton, H. Liu, Y. Liu, A. Lobodenko, M. Lokajicek, R. Lopes de Sa, R. Luna-Garcia, A. L. Lyon, A. K. A. Maciel, R. Madar, R. Magaña-Villalba, S. Malik, V. L. Malyshev, J. Mansour, J. Martínez-Ortega, R. McCarthy, C. L. McGivern, M. M. Meijer, A. Melnitchouk, D. Menezes, P. G. Mercadante, M. Merkin, A. Meyer, J. Meyer, F. Miconi, N. K. Mondal, M. Mulhearn, E. Nagy, M. Narain, R. Nayyar, H. A. Neal, J. P. Negret, P. Neustroev, H. T. Nguyen, T. Nunnemann, J. Orduna, N. Osman, A. Pal, N. Parashar, V. Parihar, S. K. Park, R. Partridge, N. Parua, A. Patwa, B. Penning, M. Perfilov, Y. Peters, K. Petridis, G. Petrillo, P. Pétroff, M.-A. Pleier, V. M. Podstavkov, A. V. Popov, M. Prewitt, D. Price, N. Prokopenko, J. Qian, A. Quadt, B. Quinn, P. N. Ratoff, I. Razumov, I. Ripp-Baudot, F. Rizatdinova, M. Rominsky, A. Ross, C. Royon, P. Rubinov, R. Ruchti, G. Sajot, A. Sánchez-Hernández, M. P. Sanders, A. S. Santos, G. Savage, M. Savitskyi, L. Sawyer, T. Scanlon, R. D. Schamberger, Y. Scheglov, H. Schellman, M. Schott, C. Schwanenberger, R. Schwienhorst, J. Sekaric, H. Severini, E. Shabalina, V. Shary, S. Shaw, A. A. Shchukin, O. Shkola, V. Simak, P. Skubic, P. Slattery, G. R. Snow, J. Snow, S. Snyder, S. Söldner-Rembold, L. Sonnenschein, K. Soustruznik, J. Stark, N. Stefaniuk, D. A. Stoyanova, M. Strauss, L. Suter, P. Svoisky, M. Titov, V. V. Tokmenin, Y.-T. Tsai, D. Tsybychev, B. Tuchming, C. Tully, L. Uvarov, S. Uvarov, S. Uzunyan, R. Van Kooten, W.M. van Leeuwen, N. Varelas, E.W. Varnes, I. A. Vasilyev, A. Y. Verkheev, L. S. Vertogradov, M. Verzocchi, M. Vesterinen, D. Vilanova, P. Vokac, H. D. Wahl, M. H. L. S. Wang, J. Warchol, G. Watts, M. Wayne, J. Weichert, L. Welty-Rieger, M. R. J. Williams, G.W. Wilson, M. Wobisch, D. R. Wood, T. R. Wyatt, Y. Xie, R. Yamada, S. Yang, T. Yasuda, Y. A. Yatsunenko, W. Ye, Z. Ye, H. Yin, K. Yip, S. W. Youn, J. M. Yu, J. Zennamo, T. G. Zhao, B. Zhou, J. Zhu, M. Zielinski, D. Zieminska, and L. Zivkovic


I. INTRODUCTION
Since the creation of the quark model [1,2] it was understood that exotic mesons containing more than one quark-antiquark pair are possible.However, for exotic mesons containing only the up, down and strange quarks it has been difficult to make a definitive experimental case for such exotic states, although some persuasive arguments have been made (for recent comprehensive discussions of * Deceased.Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license.Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.Funded by SCOAP 3 .
exotic hadrons containing both light and heavy quarks, see Refs.[3][4][5][6]).Multiquark states that contain heavy quarks can be more recognizable owing to the distinctive decay structure of heavy quark hadrons.The 2003 discovery by the Belle experiment [7] of the Xð3872Þ in the channel B AE → K AE Xð→ π þ π − J=ψÞ was the first candidate exotic meson in which heavy flavor quarks participate.This state was subsequently confirmed in several production and decay modes by ATLAS [8], BABAR [9], BES III [10], CDF [11], CMS [12], D0 [13] and LHCb [14] Collaborations.Several additional four-quark candidate exotic mesons have since been found, though in many cases not all experiments have been able to confirm their existence.
Four-quark mesons can be generically categorized as either "molecular states" or tetraquark states of a diquark and an anti-diquark.In the example of the Xð3872Þ, a molecular state interpretation would be a colorless D 0 (c ū) and a colorless DÃ0 (uc) in a loosely bound state.Such a state would be expected to lie close in mass to the D 0 DÃ0 threshold.The tetraquark mode of a colored diquark (cu) and colored anti-diquark (c ū) is more strongly bound by the exchange of gluons and would be expected to have a mass somewhat below the D 0 DÃ0 threshold.In many cases, interpretations of four-quark mesons as pure molecular or tetraquark states are difficult, and more complex mechanisms may be required [4][5][6].The firm identification of multiquark mesons and baryons and the study of their properties are of importance for further understanding of nonperturbative QCD.
Recently the D0 Collaboration presented evidence for a new four-quark candidate that decays to B 0 s π AE where B 0 s decays to J=ψϕ [15].This system would be composed of two quarks and two antiquarks of four different flavors: b, s, u, d, with either a molecular constitution as a loosely bound B 0 d and K AE system or a tightly bound tetraquark such as ðbdÞ − ðs ūÞ, ðbuÞ − ðs dÞ, ðsuÞ − ð b dÞ, or ðsdÞ − ð b ūÞ (because the B 0 s meson is fully mixed, the exact quark antiquark composition cannot be determined).The mass of X AE ð5568Þ is about 200 MeV=c 2 below the B 0 d K AE threshold, thus disfavoring a B 0 d − K AE molecular interpretation.The X AE ð5568Þ was previously reported [15] with a significance of 5.1σ (including systematic uncertainties and the look-elsewhere effect [16]) in the decay X AE ð5568Þ → B 0 s ðJ=ψϕÞπ AE in proton-antiproton collisions at a center of mass energy of 1.96 TeV.The ratio of the number of B 0 s that are from the decay of the X AE ð5568Þ to all B 0 s produced was measured to be ½8.6 AE 1.9ðstatÞ AE 1.4ðsystÞ%.In order to reduce the background, a selection was imposed on the angle between the B 0 s and π AE (the ).Without the cone cut the significance was found to be 3.9σ.In addition to increasing the signal-to-background ratio this cone cut limits backgrounds, such as possible excited states of the B c meson, that are not included in the available simulations.
Multiple checks were carried out to ensure that the cone cut did not create an anomalous signal [15].Varying the cone cut from ΔR max ¼ 0.2 to 0.5 gave stable fitted masses and resulted in no unexpected changes in the result.The invariant mass spectra of the B 0 s candidates and charged tracks with kaon or proton mass hypotheses were checked, and no resonant enhancements in these distributions were found.The invariant mass distribution of B 0 d π AE was also examined with no unexpected resonances or reflections found.Subsequent analyses by the LHCb Collaboration [18] and by the CMS Collaboration [19] have not found evidence for the X AE ð5568Þ in proton-proton interactions at ffiffi ffi s p ¼ 7 and 8 TeV.The CDF Collaboration has recently reported no evidence for X AE ð5568Þ in proton-antiproton collisions at ffiffi ffi s p ¼ 1.96 TeV [20] with different kinematic coverage than that of Ref. [15].
In this article, we present a study of the X AE ð5568Þ in the decay to B 0 s π AE using semileptonic using the full run II integrated luminosity of 10.4 fb −1 in proton-antiproton collisions at a center of mass energy of 1.96 TeV collected with the D0 detector at the Fermilab Tevatron Collider.Charge conjugate states are assumed.Here X includes the unseen neutrino and possibly a photon or π 0 from a D Ã s decay or other hadrons from the B 0 s decay.The decay process is illustrated in Fig. 1.The semileptonic decay channel has a higher branching fraction than the hadronic channel (B 0 s → J=ψϕ).However the presence of the unmeasured neutrino in the final state deteriorates the mass resolution of the signal.Still, a good mass resolution for the X AE ð5568Þ can be obtained in the semileptonic channel for events with a large invariant mass of the μ þ D − s system, yielding a comparable number of selected B 0 s candidates in the two channels.The backgrounds in the semileptonic channel are independent of, but somewhat larger than, those in the hadronic channel.The character of possible reflections of other resonant structures is quite different in the semileptonic and hadronic channels.

FIG. 1. An illustration of the decay
s X in the plane perpendicular to the beam.
Thus observation of the X AE ð5568Þ in the semileptonic decay channel enables an independent confirmation of its existence.We report here the results of the search for the X AE ð5568Þ in the semileptonic channel, as well as a combination of the results in the hadronic and semileptonic channels.

II. D0 DETECTOR
The detector components most relevant to this analysis are the central tracking and the muon systems.The D0 detector has a central tracking system consisting of a silicon microstrip tracker (SMT) and a central fiber tracker (CFT), both located within a 2 T superconducting solenoidal magnet [21,22].The SMT has a design optimized for tracking and vertexing for pseudorapidity of jηj < 3.For charged particles, the resolution on the distance of closest approach as provided by the tracking system is approximately 50 μm for tracks with p T ≈ 1 GeV=c, where p T is the component of the momentum perpendicular to the beam axis.It improves asymptotically to 15 μm for tracks with p T > 10 GeV=c.Preshower detectors and electromagnetic and hadronic calorimeters surround the tracker.A muon system, positioned outside the calorimeter, covering jηj < 2 consists of a layer of tracking detectors and scintillation trigger counters in front of 1.8 T iron toroidal magnets, followed by two similar layers after the toroids [23].

III. EVENT RECONSTRUCTION AND SELECTION
The B 0 s → μ þ D − s X selection requirements have been chosen to optimize the mass resolution of the B 0 s π þ system and to minimize background from random combinations of tracks from muons and charged hadrons.The selection criteria are based on those used in Ref. [24] with the cut on the B 0 s isolation removed and have been selected by maximizing the significance of the signal.
The data were collected with a suite of single and dimuon triggers (approximately 95% of the sample is recorded using single muon triggers).The selection and reconstruction of μ þ D − s decays requires tracks with at least two hits in both the CFT and SMT.
The muon is required to have hits in at least two layers of the muon system, with segments reconstructed both inside and outside the toroid.The muon track segment is required to be matched to a track found in the central tracking system that has transverse momentum 3 < p T < 25 GeV=c.The D − s → ϕπ − ; ϕ → K þ K − decay is selected as follows.The two particles from the ϕ decay are assumed to be kaons and are required to have p T > 1.0 GeV=c, opposite charge and an invariant mass 1.012 < mðK þ K − Þ < 1.03 GeV=c 2 .The charge of the third particle, assumed to be a pion, has to be opposite to that of the muon.This particle is required to have transverse momentum 0.5 < p T < 25 GeV=c.The mass of the three particles must satisfy 1.91 The three tracks are combined to form a common D − s decay vertex using the algorithm described in Ref. [25].The D − s vertex is required to be displaced from the p p primary interaction vertex (PV) in the transverse plane with a significance of at least three standard deviations.The cosine of the angle between the D − s momentum and the vector from the PV to the D − s decay vertex is required to be greater than 0.9.
The trajectories of the muon and D − s candidate are required to be consistent with originating from a common vertex (assumed to be the B 0 s semileptonic decay vertex).The cosine of the angle between the combined μ þ D − s transverse momentum, an approximation of the B 0 s direction, and the direction from the PV to the B 0 s decay vertex has to be greater than 0.95.The B 0 s decay vertex has to be displaced from the PV in the transverse plane with a significance of at least four standard deviations.The transverse momentum of the μ þ D − s system is required to satisfy the condition p T > 10 GeV=c to suppress backgrounds.To minimize the effect of the neutrino in the final state the effective mass is limited to 4.5 The impact parameters (IP) [26] with respect to the PVof the four tracks from the B 0 s decay are required to satisfy the following criteria: the two-dimensional (2D) IPs of the tracks of the muon and the pion from the D − s decay are required to be at least 50 μm to reject tracks emerging promptly from the PV (this requirement is not applied to the tracks associated with the charged kaons since the mass requirements provide satisfactory background suppression).The three-dimensional (3D) IPs of all four tracks are required to be less than 2 cm to suppress combinations with tracks emerging from different p p vertices reconstructed in the same beam crossing. The is shown in Fig. 2, where the invariant mass distribution in data is compared to a fit using a function which includes three terms: a second order polynomial used to describe combinatorial background, a Gaussian used to model the D − peak, and a double Gaussian with similar, but different masses and widths used to model the D − s peak.The selection criteria for the pion in the B 0 s π AE combination have been chosen to match those used in the hadronic analysis.The track representing the pion is required to have transverse momentum 0.5 < p T < 25 GeV=c (the upper limit is applied to reduce background).The pion and the B 0 s candidate are combined to form a vertex that is consistent with the PV.The pion is required to be associated with the PV and have a 2D IP of at most 200 μm and a 3D IP that is less than 0.12 cm.Events with more than 20 B 0 s π AE candidates are rejected.The most likely number of candidates per event is 5.1, and only about 0.1% of the events have more than 20 candidates per event.To improve the resolution of the invariant mass of the B 0 s π AE system we define the invariant mass as mðB 0 [27].We study the mass distribution in the range 5.506 < mðB 0 s π AE Þ < 5.906 GeV=c 2 .When using the hadronic data from Ref. [15] in this paper we use the same mass range as the semileptonic data instead of the slightly shifted mass range used in the original analysis [5.5 < mðB 0 s π AE Þ < 5.9 GeV=c 2 ].The semileptonic data are studied with and without a cone cut which is used to suppress background, in which the angle between the μ þ D − s system and π AE is required to satisfy ΔR ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi Δη 2 þ Δϕ 2 p < 0.3.The resulting invariant mass distributions for the semileptonic channel are shown in Fig. 3.
The selection cuts and resulting kinematics for the hadronic and semileptonic channels are quite similar.The requirement that muons be seen outside the toroids means that the minimum p T for the J=ψ in the hadronic channel is about 4 GeV and about 3 GeV for the single muon in the semileptonic channel.The minimum p T for the additional pion is 0.5 GeV for both the hadronic and semileptonic channels.For both channels, we require the minimum p T ðB 0 s πÞ to be greater than 10 GeV and the average p T ðB 0 s πÞ for events with mðB s πÞ ≈ 5.5 GeV is ≈17 GeV.For both channels the B 0 s π candidates are in the range of −2 < η < 2, and more than half of the events have a muon with jηj > 1.

IV. MONTE CARLO SIMULATION, BACKGROUND MODELING AND PARAMETRIZATION
Monte Carlo (MC) samples are generated using the PYTHIA [28] event generator, modified to use EVTGEN [29] for the decay of hadrons containing b or c quarks.The generated events are processed by the full detector simulation chain.Data events recorded in random beam crossings are overlaid on the MC events to simulate the effect of additional collisions in the same or nearby bunch crossings.The resulting events are then processed with the same reconstruction and selection algorithms as used for data events.
The MC sample for X AE ð5568Þ signal is generated by modifying the mass of the B AE meson and forcing it to decay to B 0 s π AE using an isotropic S-wave decay model.The X AE ð5568Þ is simulated with zero width and zero lifetime.The resulting K þ K − π − and B 0 s π AE invariant mass distributions are shown in Fig. 4 with all selection requirements.
The signal component of the K þ K − π AE invariant mass distribution (Fig. 4a) is modeled by two Gaussian functions and the background by a second-order polynomial.The signal of the mðB 0 s π AE Þ distribution (Fig. 4b) is well modeled with a single Gaussian and the background with a third-order polynomial times an exponential.Using the results of these fits the reconstruction efficiency of the charged pion in the decay GeV=c where the systematic uncertainty represents the expected differences between the reconstruction efficiencies for low-momentum tracks in the data and MC simulation.
It is not possible to create a model of the background that is based only on data.Since the X AE ð5568Þ decays to B 0 s mesons, any data sample that includes B 0 s decays will also include the signal and is unsuitable for modeling the background.Hence, we use MC-generated B 0 s events that result from known particles that have decays that include a B 0 s in the decay chain, combined with data events where the muon has the same sign as the D − s candidate (SS events).MC event generators do not include all possible states as in There are two distinct sources of background in this analysis.The first occurs when an X AE ð5568Þ candidate is reconstructed from a real μ þ and D − s together with a random charged track.This background is modeled using MC samples.
The background MC sample is generated using the PYTHIA inclusive heavy flavor production model, and events are selected that contain at least one muon and a D ∓ s → ϕπ ∓ decay where ϕ → K þ K − .To correct for the difference in lifetimes in the MC simulation and data, a weighting is applied to all nonprompt events in the simulation, based on the generated lifetime of the B candidate, to give the world-average B hadron lifetimes [27].To correct for the effects of the trigger selection and the reconstruction in data, we also weight each MC event so that the transverse momenta of the reconstructed muon and the μ þ D − s system agree with those in the data.The p T distribution of the B 0 s π AE system is altered significantly by the weighting as shown in Fig. 5(a).However, the effect is relatively small for the B 0 s π AE mass distribution as seen in Fig. 5(b).The second source of background is the combinatorial background that occurs when a X AE ð5568Þ candidate is reconstructed from a spurious D − s candidate formed from three random charged tracks that form a vertex.This background is modeled using data events where the muon has the same sign as the D − s candidate (SS events).In Fig. 6(b) we compare the reweighted MC background simulation, smoothed using one iteration of the 353QH algorithm [30], with the SS data for the no cone cut case.These two backgrounds are in good agreement since the χ 2 between them is 50 for 50 bins.We therefore choose to use the MC background shape only, for the data without the cone cut.In Fig. 6(a) we make the same comparison for the data with the cone cut.In this case, χ 2 ¼ 77 for the 50 bins, and we therefore need to model the background shape with a combination of the MC and SS backgrounds.
To construct the background sample for the data with the cone cut the fraction of MC and SS backgrounds need to be determined.This is found by fitting the data with a combination of the MC and SS with the fraction of MC events as a free parameter in the sideband mass range 5.506 < mðB 0 s π AE Þ < 5.55 and 5.650 < mðB 0 s π AE Þ < 5.906 GeV=c 2 .The best agreement is found when the MC fraction is ð62 AE 2Þ%.
We choose the background parametrization for the invariant mass distribution, both with and without the cone cut, to be 5063 GeV=c 2 is the mass threshold.Our baseline choice of Eq. ( 1) gives an equivalently good description of the background as that used in Ref. [15] [Eq.(2)].It has the advantages of having one fewer parameter and being zero at the mass threshold.
Three alternative parametrizations are used to model the background.The first is that used in Ref. [15], where m Δ ¼ m − Δ and Δ ¼ 5.500 GeV=c 2 .The second is the ARGUS function [31] which is specifically constructed to describe background near a threshold, The third alternative model used to fit the background is the MC histogram (or combined MC and SS data) smoothed using one iteration of the 353QH algorithm [30].The ARGUS function is not used as an alternate parametrization in the semileptonic data with the cone cut, because the fit to background is strongly disfavored (the χ 2 of the fit to the MC background is 145 compared with approximately 50 for the alternate functions).The χ 2 per number of degrees of freedom (ndf) for the four representations of the background are shown in Table I.We choose the background description of Eq. ( 1) as the baseline.The alternative functions and the smoothed MC are used to estimate the systematic uncertainty on the background shape.The mðB 0 s π AE Þ background model distribution along with the fit using Eq. ( 1) is presented in Fig. 7.

V. SIGNAL MASS RESOLUTION
We calculate the mass of the B 0 s π AE system using the quantity, Before carrying out the search for the X AE ð5568Þ in the semileptonic channel we ensure that it is an unbiased and precise estimator of the mass of the B 0 s π AE system.This is studied by simulating the two body decay Xð5568Þ AE → B 0 s π AE where B 0 s → μ AE D ∓ s X, starting with a range of input masses mðB 0 s π AE Þ.Following the decay chain B 0 s → μ AE D ∓ s X and forming the invariant masses mðμ AE D ∓ s π AE Þ and mðμ AE D ∓ s Þ are found.Then mðB 0 s π AE Þ is calculated and compared to the input mass mðB 0 s π AE Þ.To evaluate how well the mass approximation works to compensate for the missing neutrino, we model the decay with a toy MC that simulates the virtual with an isotropic distribution of μ and ν in the W boson rest frame.The resulting resolution of a zero width resonance due to the presence of the neutrino is modeled by a Gaussian.The width varies according to mðB 0 s π AE Þ as illustrated by the solid line in Fig. 8.
The mass resolution for the D0 detector of a state decaying into five reconstructed charged particles with a similar kinematic range as in this study is measured using the MC simulation and is given by a Gaussian function of width 3.85 MeV=c 2 .The mðB 0 s π AE Þ resolution function is obtained by convoluting the Gaussian tracking resolution and the smearing resolution resulting from the missing neutrino.The resulting combined resolution, the dashed line in Fig. 8, can be approximated by where m 0 has the same definition as in Eq. ( 1).These studies show that the difference between mðB 0 s π AE Þ and mðB 0 s π AE Þ is less than 1 MeV=c 2 in the search region.This is confirmed with the signal MC sample.

VI. SIGNAL FIT FUNCTION
The X AE ð5568Þ resonance is modeled by a relativistic Breit-Wigner function convolved with a Gaussian detector resolution function given in Eq. ( 5), F sig ðm; m X ; Γ X Þ,  where m X and Γ X are the mass and the width of the resonance.
The fit function has the form where f sig and f bgr are normalization factors.The shape parameters in the background term F bgr are fixed to the values obtained from fitting the MC background distribution (see Fig. 7).
We use the Breit-Wigner parametrization appropriate for an S-wave two-body decay near threshold, The mass-dependent width ΓðmÞ ¼ Γ X • ðq 1 =q 0 Þ, where q 1 and q 0 are the magnitudes of momenta of the B 0 s meson in the rest frame of the B 0 s π AE system at the invariant mass equal to m and m X , respectively.

VII. X AE ð5568Þ SEMILEPTONIC FIT RESULTS
In the fit to the semileptonic data with the cone cut shown in Fig. 9(a), the normalization parameters f sig and f bgr and the Breit-Wigner parameters m X and Γ X are allowed to vary.The fit yields the mass and width of m X ¼ 5566.4 þ3. 4  −2.8MeV=c 2 , Γ X ¼ 2.0 þ9.5 −2.0 MeV=c 2 , the number of signal events, N ¼ 121 þ51 −34 , and a χ 2 ¼ 34.9 for 46 degrees of freedom.The local statistical significance of the signal is defined as ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi −2 lnðL 0 =L max Þ p , where L max and L 0 are likelihood values at the best-fit signal yield and the signal yield fixed to zero obtained from a binned maximum-likelihood fit.The p-value of the background only fit is 2.1 × 10 −5 , and the local statistical significance is 4.3σ.
In the fit to the semileptonic data without the cone cut shown in Fig. 9(b), the mass and width of m X ¼ 5566.7 þ3.6 −3.4 MeV=c 2 , Γ X ¼ 6.0 þ9.5 −6.0 MeV=c 2 , the number of signal events, N ¼ 139 þ51 −63 , and a χ 2 ¼ 30.4 for 46 degrees of freedom.The p-value of the background only fit is 7.7 × 10 −6 , and the local statistical significance is 4.5σ.The fit results, both for the cone cut and no cone cut cases, are given in Table II and for various background parametrizations in Table III.The X AE ð5568Þ parameters for the cone cut and no cone cut cases are consistent.IV) are obtained for the measured values of the mass, width and event yield of the X AE ð5568Þ signal.The dominant uncertainty is due to TABLE II.Results for the fit to the semileptonic data sets (see Fig. 9).

Cone cut No cone cut
Fitted mass, MeV=c (i) the description of the background shape.We evaluate this systematic uncertainty by using the alternative parametrizations of the background, Eqs. ( 2) and ( 3) and the smoothed MC histogram and finding the maximal deviations from the nominal fit.The effect of (ii) the MC weighting is estimated by creating 1000 background samples where the weights have been randomly varied based on the uncertainties in the weighting procedure.
Other sources of systematic uncertainty are evaluated by (iii) varying the energy scale in the MC sample relative to the data by AE1 MeV=c 2 , (iv) varying the mass resolution of the X AE ð5568Þ either by AE1 MeV=c 2 around the mean value, or by using a constant resolution of 11.1 MeV=c 2 obtained from the MC simulation of the X AE ð5568Þ signal, (v) using a P-wave relativistic Breit-Wigner function, and (vi) estimating the shift of the fitted mass peak due to the missing neutrino.
Systematic uncertainties are summarized in Table IV.The uncertainties are added in quadrature separately for positive and negative values to obtain the total systematic uncertainties for each measured parameter.The results including systematic uncertainties are given in Table II.

B. Significance
Since we are seeking to confirm the result presented in Ref. [15] we do not apply a correction for a look elsewhere effect.The systematic uncertainties are treated as nuisance parameters to construct a prior predictive model [27,32] of our test statistic.When the systematic uncertainties are included, the significance of the observed semileptonic signal with the cone cut is 3.2σ (p-value ¼ 1.4 × 10 −3 ).The significance of the semileptonic signal without the cone is 3.4σ (p-value ¼ 6.4 × 10 −4 ).Fitted width, MeV=c 2 6.0 þ9.5 −6.0 6.0 þ12.0 −6.0 6.5 þ8.9 −6.5

C. Closure tests
We have tested the accuracy of the fitting procedure using toy MC event samples constructed with input mass and width of 5568.3 and 21.9 MeV=c 2 , respectively, with the number of input signal events varied in steps of 25 between 75 and 350.At each number of input signal events, 10,000 pseudoexperiments were generated.The signals are modeled with a relativistic Breit-Wigner function convolved with a Gaussian function representing the appropriate detector resolution.The background distribution is based on Eq. ( 1).For each trial the fitting procedure is performed to obtain the mass and width and the number of semileptonic signal events.The results of each set of trials is fitted with a Gaussian to determine the mean and the uncertainty in the number of signal events, the mass and the width (see Table V).The number of fitted signal events vs the number of injected signal events for the semileptonic samples are plotted in Fig. 10.
For the ensembles with a number of input events similar to that observed in data, there is a slight overestimate of the yield and fitted mass, and the width is underestimated.This width reduction is in agreement with the results of the fits to data (Sec.VII) and indicate that the semileptonic data are not sensitive to the width.These effects are accounted for in the calculation of the significance.

D. Comparison with hadronic channel
The measured values of the mass, width, the number of signal events, and significance of the signal for the semileptonic channel and the hadronic channel [15] are given in Table VI.The mass and width of the X AE ð5568Þ for the semileptonic and hadronic channels are consistent taking into account the uncertainties.The observed yields are consistent with coming from a common particle given the number of B 0 s events in the sample and the B 0 s branching ratios.

E. Cross-checks
As a cross-check the B 0 s π AE mass-bin size is set to 5 MeV=c 2 and to 10 MeV=c 2 instead of 8 MeV=c 2 , and the lower edge of the fitted mass range is shifted by 2, 3, 5, and 7 MeV=c 2 .This leads to maximal variations in the mass of þ0.1 −0.6 MeV=c 2 , in the width of þ1.7 −0.9 MeV=c 2 and in the number of signal events þ0 −9 which are small compared to the statistical and systematic uncertainties.
To test the stability of the results, alternative choices are made regarding the fit parameters.In the first, the background fit parameters are allowed to float.The resulting fit is consistent with the nominal fit and the p-value of the background-only fit is 1.7 × 10 −4 corresponding to a statistical significance of 3.8σ (Table VII).The second cross-check fixes the mass and width of the X AE ð5568Þ to the values found in Ref. [15].Again, the resulting fit is consistent with the nominal fit with an increase in the number of signal events due to the increased width of the peak.The p-value of the background-only fit is 1.1 × 10 −4 corresponding to a statistical significance of  3.9σ (Table VII).These cross-checks are also repeated without the cone cut (Table VII).

VIII. PRODUCTION RATIO OF X AE ð5568Þ TO B 0 s
To calculate the production ratio of the X AE ð5568Þ to B 0 s , the number of the B 0 s -mesons needs to be estimated.The fitting of the K þ K − π ∓ mass distribution is described in Sec.IV.The number of D ∓ s mesons extracted from the fit and adjusted for the mass range 1.91 2).The number of μ AE D ∓ s events in the signal sample that are the result of a random combination of a promptly produced D ∓ s and a muon in the event is estimated using events where the muon and the D ∓ smeson have the same sign.The same sign data sample is analyzed using the same model as the opposite sign data with the means and widths of the Gaussians fixed to the values obtained from the opposite sign data.The number of events in the same-sign sample is NðD AE s Þ ¼ 426 AE 61.The mass distributions of the K þ K − π ∓ for opposite and same-sign data are shown in Fig. 2.
The number of B 0 s -meson decays in the semileptonic data is estimated by subtracting the contribution of the promptly produced μ AE D ∓ s events from the overall μ AE D ∓ s sample.A study of the MC background simulations shows that the purity of the resulting sample is 99.5 þ0.5 −1.0 %.We find 6222 AE 141 B 0 s events.Combining these results and using the efficiency for the charged pion in the X(5568) decay (Sec.IV), we obtain a production ratio for the semileptonic data of for our fiducial selection [which includes the requirements where N sl ðX AE ð5568Þ → B 0 s ðslÞπ AE Þ is the number of X AE ð5568Þ decays to B 0 s π AE and N B 0 s ðslÞ is the inclusive number of B 0 s decays, both for semileptonic decays of the B 0 s .This result is similar to the ratio measured in TABLE VI.Fit results obtained in the semileptonic channel and in the hadronic channel (Ref.[15]).In the hadronic channel with no cone cut the mass and width of the X AE ð5568Þ were set to the values found with the cone cut.LEE-look elsewhere effect.

IX. COMBINED SIGNAL EXTRACTION
We now proceed to fit the hadronic and semileptonic data sets simultaneously.The hadronic data set is the same as used in Ref. [15] except that the data are fitted in the mass range 5.506 < mðB 0 s π AE Þ < 5.906 GeV=c 2 instead of 5.500 < mðB 0 s π AE Þ < 5.900 GeV=c 2 .The data selection and background modeling for the hadronic data set are described in detail in Ref. [15].
The fit function has the form where f hðslÞ;sig and f hðslÞ;bgr are normalization factors.The shape parameters in the background terms F hðslÞ;bgr are fixed to the values obtained from fitting the respective background models for the hadronic (h) and semileptonic (sl) samples to Eq. ( 1).The signal shape F hðslÞ;sig is modeled by relativistic Breit-Wigner function convolved with a Gaussian detector resolution function that depends on the data sample.For the semileptonic sample the detector resolution is given by Eq. ( 5), and for the hadronic channel it is 3.85 MeV=c 2 .For the data without the cone cut the combined data are fitted in the range 5.506 < mðB 0 s π AE Þ < 5.706 GeV=c 2 as the hadronic background is not well modeled for mðB 0 s π AE Þ > 5.706 GeV=c 2 [15].The same Breit-Wigner parameters m X and Γ X are used for the hadronic and semileptonic samples.In the fits shown in Fig. 11, the normalization parameters f hðslÞ;sig and f hðslÞ;bgr and the Breit-Wigner parameters m X and Γ X are allowed to vary.Since the fraction of B 0 s events produced by the decay of the X AE ð5568Þ should be essentially the same in the hadronic and semileptonic channels the X AE ð5568Þ event yields (N h and N sl ) are constrained using the parameter which is required to be consistent with the B 0 s -meson production rate in the hadronic and semileptonic channels where N B 0 s ðslÞ ¼ 6222 AE 144, N B 0 s ðhÞ ¼ 5582 AE 100 are the number of semileptonic and hadronic B 0 s decays in the sample.A likelihood penalty of 0.5½ðA sl;h − A sl;h ðB 0 s ÞÞ= ΔA sl;h ðB 0 s Þ 2 is applied where ΔA sl;h ðB 0 s Þ ¼ 0.020 is the uncertainty.This uncertainty includes the statistical uncertainty in the number of B 0 s events and the uncertainties in the relative reconstruction efficiencies and acceptances between the hadronic and semileptonic data.A ratio has been chosen for the constraint as it is well behaved if either of the event yields (N h and N sl ) approaches zero.
The fit results are summarized in Table VIII, and the correlations between the fit parameters are given in Table IX.The correlation of nearly one between N X ðslÞ and N X ðhadÞ is a result of the constraint on the event yields [Eq.(11)].The local statistical significance of the signal is defined as ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi −2 lnðL 0 =L max Þ p , where L max and L 0 are likelihood values at the best-fit signal yield and the signal yield fixed to zero obtained from a binned maximumlikelihood fit.For the cone cut the p-value of the fit to the data with the cone cut is 2.2 × 10 −14 and the local statistical  VIII).The background parametrization function is taken from Eq. (1).significance is 7.6σ.The p-value without the cone cut is 8.2 × 10 −9 , and the local statistical significance is 5.8σ.

A. Systematic uncertainties
The systematic uncertainties of the combined fit are given in Table X.The uncertainty on (i) the background shape descriptions is evaluated by using the alternative parametrizations of the background, Eqs. ( 2) and ( 3) and the smoothed MC histogram independently for the semileptonic and the hadronic channels (16 different fits) and finding the maximal deviations from the nominal fit.
The effect of (ii) the MC weighting for the semileptonic background is estimated by creating 1000 background samples where the weights have been randomly varied based on the uncertainties in the weighting procedure and measuring the standard deviation and bias of the measured values.
The (iii) MC component of the background for the hadronic sample is made up of a mixture of two independent MC samples with different production properties (see Ref. [15]), and the systematic uncertainties due to this are found by varying the composition of this mixture and measuring the standard deviation and bias of the measured values.The (iv) size of the hadronic sidebands is varied using the maximal deviations from the nominal fit to estimate the systematic uncertainty.
The systematic uncertainty due to the (v) fraction of MC and SS data in the semileptonic sample, (vi) the MC and sideband data in the case of the hadronic, is varied independently between the two samples measuring the standard deviation and bias of the measured values.Since the background model for the semileptonic sample without the cone cut only uses the MC background simulation this uncertainty (v) does not apply.
All of the uncertainties due to the modeling of the background are assumed to be independent for the hadronic and semileptonic data samples.
The remaining uncertainties are measured by finding the maximal deviations from the nominal fit for (vii) varying the energy scale in the semileptonic and hadronic MC data samples by AE1 MeV=c 2 in both samples simultaneously; (viii) varying the nominal mass resolution of 3.85 MeV=c 2 for the D0 detector by AE1 MeV=c 2 and þ2 MeV=c 2 in both the hadronic and semileptonic data samples simultaneously; (ix) varying the resolution of the X AE ð5568Þ peak in the semileptonic channel either by AE1 MeV=c 2 around the mean value given by Eq. ( 5) or by using a constant resolution of 11.1 MeV=c 2 for the semileptonic data while the mass resolution in the hadronic channel remains at 3.85 MeV=c 2 ; (x) using a P-wave relativistic Breit-Wigner function for both data sets; (xi) setting the shift of the fitted mass peak in the semileptonic data with respect to the hadronic data due to the missing neutrino to AE1 MeV=c 2 ; and (xii) varying the constraint on the relative number of signal events in hadronic and semileptonic channels [Eq.(11)] between 0.034 and 0.074.The correlation of each of the sources of systematic uncertainty between the hadronic and semileptonic data sets is indicated in Table X.The uncertainties are added in quadrature separately for positive and negative values to obtain the total systematic TABLE VIII.Results for the combined fit to the hadronic and semileptonic data sets (see Fig. 11).

Cone cut
No cone cut Fitted mass, MeV=c 2 5566.9IX.Correlations between the parameters of the combined fit to the hadronic and semileptonic data sets (see Fig. 11).
The yield in the semileptonic channel is N X ðslÞ, the hadronic channel N X ðhÞ, while the fraction of background events is f sl;bgr and f h;bgr , respectively.
small reference value) is measured.The mean number of upcrossings for a reference level of 0.5 is determined, and the global significance is calculated.The resulting significance is consistent with the method described above.The systematic uncertainties are treated as nuisance parameters to construct a prior predictive model [27,32] of our test statistic.When the systematic uncertainties are included, the significance of the observed signal with the cone cut applied for the combined fit is reduced to 6.7σ (pvalue = 1.5 × 10 −11 ), and the significance of the signal without the cone cut is 4.7σ (p-value ¼ 2.0 × 10 −6 ).

C. Closure tests
To test the sensitivity and accuracy of the fitting procedure for the combined signal extraction we repeat the closure tests carried out in Sec.VII C with the following modifications.The size of the associated hadronic signal is set using Eqs.(11) and (12).The appropriate detector resolution is used, Eq. ( 5) for the semileptonic sample and 3.85 MeV=c 2 for the hadronic sample.For each trial the fitting procedure is performed to obtain the mass and width and the number of semileptonic and hadronic signal events.The results of each set of trials is fitted with a Gaussian to determine the mean and the uncertainty in the number of signal events, the mass and the width (see Table XI).The number of fitted signal events vs the number of injected signal events for the semileptonic and hadronic samples is plotted in Fig. 12.These results show excellent agreement between the input and fit parameters.

D. Cross-checks
To test the stability of the results, alternative choices are made regarding the fit parameters (see Table XII).
When no constraint is placed on the ratio of the event yields in the hadronic and semileptonic channels, Eq. ( 11), the results are entirely consistent with the fit with the constraint.
We have also carried out a fit in which two of the systematic effects are treated as nuisance parameters in the fit.We allow a mass shift, Δm, between the hadronic and  semileptonic data with a likelihood penalty of 0.5ðΔm=1 MeV=c 2 Þ 2 .We also allow the overall resolution of the semileptonic signal to vary by Δσ SL with a likelihood penalty of 0.5ðΔσ SL =1 MeV=c 2 Þ 2 .The resultant fit produces a mass, width and event yields that are consistent with the default fit and shifts of Δm ¼ ð0.0 AE 1.4Þ MeV=c 2 and Δσ SL ¼ ð−0.1 AE 1.4Þ MeV=c 2 .The significance of a nonzero width is determined by fitting the data with the width set to zero and comparing it with the fit with no constraint on the width (Table XII).Using the data with the cone cut the p-value of the width being consistent with zero is 5.4 × 10 −6 , and the statistical significance is 4.5σ.The significance without the cone cut is 3.3σ (p-value = 1.1 × 10 −3 ).

X. CONCLUSIONS
We have presented the results of a search for the X AE ð5568Þ → B 0 s π AE with semileptonic decays of the B 0 s meson.The X AE ð5568Þ → B 0 s π AE state reported in the case that B 0 s → J=ψϕ [15] is confirmed for the case that B 0 s → μ ∓ D AE s X, D AE s → ϕπ AE .The analyses of the hadronic and semileptonic data give similar measurements of the mass, width and production ratio of X AE ð5568Þ to a B 0 s meson.The mass and width of this state are measured using a combined fit of both data sets with the cone cut, yielding m ¼ 5566.9 þ3.2 −3.1 ðstatÞ þ0.6 −1.2 ðsystÞ MeV=c 2 , Γ ¼ 18:6 þ7.9 −6.1 ðstatÞ þ3.5 −3.8 ðsystÞ MeV=c 2 .The p-value for the null signal hypothesis to represent the data is 1.5 × 10 −11 (6.7σ).

a
Visitor from Augustana College, Sioux Falls, SD 57197, USA.b Visitor from University of Liverpool, Liverpool L69 3BX, United Kingdom.c Visitor from Deutshes Elektronen-Synchrotron (DESY), Notkestrasse 85, Germany.d Visitor from CONACyT, M-03940 Mexico City, Mexico.e Visitor from SLAC, Menlo Park, CA 94025, USA.f Visitor from University College London, London WC1E 6BT, United Kingdom.g Visitor from Centro de Investigacion en Computacion-IPN, CP 07738 Mexico City, Mexico.h Visitor from Universidade Estadual Paulista, São Paulo, SP 01140, Brazil.i Visitor from Karlsruher Institut für Technologie (KIT)-Steinbuch Centre for Computing (SCC), D-76128 Karlsruhe, Germany.j Visitor from Office of Science, U.S. Department of Energy, Washington, D.C. 20585, USA.k Visitor from American Association for the Advancement of Science, Washington, D.C. 20005, USA.l Visitor from Kiev Institute for Nuclear Research (KINR), Kyiv 03680, Ukraine.m Visitor from University of Maryland, College Park, MD 20742, USA.n Visitor from European Organization for Nuclear Research (CERN), CH-1211 Geneva, Switzerland.o Visitor from Purdue University, West Lafayette, IN 47907, USA.p Visitor from Institute of Physics, Belgrade, Belgrade, Serbia.q Visitor from P.N.Lebedev Physical Institute of the Russian Academy of Sciences, 119991, Moscow, Russia.

FIG. 2 .FIG. 3 .
FIG.2.The K þ K − π AE invariant mass distribution for the μ AE ϕπ ∓ sample (right sign) with the solid curve representing the fit.The lower mass peak is due to the decay D AE → ϕπ AE and the second peak is due to the D AE s decay.The blue histogram below the data points is the invariant mass distribution for the same-sign sample, μ AE ϕπ AE .

FIG. 5 .
FIG. 4. MC simulation ofX AE ð5568Þ → B 0 s π AE where B 0 s → μ þ D − s X and the width of the X AE ð5568Þ is zero.The invariant mass distributions (a) mðK þ K − π − Þ and (b) mðB 0 s π þ Þ are shown.The background in the mðB 0 s π þ Þdistribution is produced by the combination of a random charged track with the B 0 s meson.

FIG. 6 .
FIG.6.The comparison of the mðB 0 s π AE Þ background only distributions (a) with the cone cut and (b) without the cone cut, obtained using the weighted MC (histogram) and from the same sign data samples (points with error bars).The fluctuations in the number of MC events with the cone cut are due to the weighting procedure and the size of the sample.

FIG. 7 .
FIG. 7. The background model produced according to the procedure described in the text is shown along with background function (1) (dotted line) (a) with and (b) without the cone cut.The gray band shows the systematic uncertainties on the background model (see Sec. VII D).

8 .
The resolution for a zero width resonance as a function of mðB 0 s π AE Þ.The solid circles and the solid line show the effect of the missing neutrino, and the open squares and dashed line show the convolution of the resolution due to the missing neutrino convolved with the 3.85 MeV=c 2 detector mass resolution.

FIG. 9 .
FIG. 9.The mðB 0 s π AE Þ distribution (a) with and (b) without the cone cut.The fitting function is superimposed (see text for details).

FIG. 10 .
FIG.10.Results of the toy MC tests of the fitting procedure (black circles) used in the analysis of the semileptonic data with the cone cut.The number of fitted signal events are plotted vs fitted number of injected signal events.The dotted line shows N in ¼ N out .

FIG. 11 .
FIG.11.The mðB 0 s π AE Þ distribution for the hadronic (red squares) and semileptonic (black circles) data with the combined fitting function superimposed (a) with and (b) without the cone cut.(see text for details, the resulting fit parameters are given in TableVIII).The background parametrization function is taken from Eq. (1).
FIG.12.Results of the toy MC tests of the combined sample fitting procedure (black circles) used in the analysis with the cone cut.The number of fitted signal events are plotted vs fitted number of injected signal events for the (a) semileptonic and (b) hadronic samples.The dotted line shows N in ¼ N out and the red line shows the fit to a line.

TABLE I .
Fit results for different parametrizations to the background model.

TABLE III .
Semileptonic data fits for the different background parametrizations.

TABLE IV .
Systematic uncertainties for the X AE ð5568Þ state mass, width and the event yield obtained from the semileptonic data.

TABLE V .
Mean values and uncertainties for fitted number of events, mass and width from Gaussian fits to corresponding distributions from 10,000 trials with the cone cut.Also given is the expected statistical uncertainties on the fitted number of events, ΔðN fit ðslÞÞ, and the expected uncertainties on the measurement of the width, ΔðΓ X Þ MeV=c 2 .A range of signals with 75, 100, 125, 150, 175 and 200 signal events, mass m x ¼ 5568.3MeV=c 2 and width Γ X ¼ 21.9 MeV=c 2 have been simulated.Background parametrization Eq. (1) is used.

TABLE VII .
[15]results for the semileptonic channel using parametrization(1)with the nominal fit, with all parameters free and the mass and width fixed to those of the hadronic channel.Statistical uncertainties only.½8.6 AE 1.9ðstatÞ AE 1.4ðsystÞ% for p T ðJ=ψϕπ AE Þ > 10 GeV=c 2[15].

TABLE XI .
Mean values and uncertainties for fitted number of events, mass and width from Gaussian fits to corresponding distributions from 10,000 trials with the cone cut.Also given is the expected statistical uncertainties on the fitted number of events, ΔðN fit Þ, and the expected uncertainties on the measurement of the width, ΔðΓ X Þ MeV=c 2 .A range of signals with 75, 100, 125, 150, 175 and 200 signal events, mass m x ¼ 5568.3MeV=c 2 and width Γ X ¼ 21.9 MeV=c 2 have been simulated.Background parametrization Eq. (1) is used.