# First Observation of the Decay 𝑩→𝑱/𝝍⁢ϕ⁢𝑲

## Abstract

We present the first observation of the decay $`BJ/\psi \varphi K`$. Using $`9.6\times 10^6`$ $`B\overline{B}`$ meson pairs collected with the CLEO detector, we have observed 10 fully reconstructed $`BJ/\psi \varphi K`$ candidates, whereas the estimated background is $`0.5\pm 0.2`$ events. We obtain a branching fraction of $`(BJ/\psi \varphi K)=(8.8_{3.0}^{+3.5}[\mathrm{stat}]\pm 1.3[\mathrm{syst}])\times 10^5`$. This is the first observed $`B`$ meson decay requiring the creation of an additional $`s\overline{s}`$ quark pair.

An observation of a $`B`$ meson decay requiring the creation of an additional $`s\overline{s}`$ quark pair in the final state would enhance our understanding of strong interactions in the final states of $`B`$ decays. Previous studies of such processes involved searches for the “lower vertex” $`\overline{B}D_s^+X`$ transitions , however no signal was observed. The decay $`BJ/\psi \varphi K`$  can occur only if an additional $`s\overline{s}`$ quark pair is created in the decay chain besides the quarks produced in the weak $`bc\overline{c}s`$ transition. The $`BJ/\psi \varphi K`$ transition most likely proceeds as a three-body decay (Fig. 1). Another possibility is that the $`BJ/\psi \varphi K`$ decay proceeds as a quasi-two-body decay in which the $`J/\psi `$ and $`\varphi `$ mesons are daughters of a hybrid charmonium state .

We searched for $`B^+J/\psi \varphi K^+`$ and $`B^0J/\psi \varphi K_S^0`$ decays, reconstructing $`J/\psi \mathrm{}^+\mathrm{}^{}`$, $`\varphi K^+K^{}`$, and $`K_S^0\pi ^+\pi ^{}`$. Both $`e^+e^{}`$ and $`\mu ^+\mu ^{}`$ modes were used for the $`J/\psi `$ reconstruction. The data were collected at the Cornell Electron Storage Ring (CESR) with two configurations of the CLEO detector, called CLEO II  and CLEO II.V. The components of the CLEO detector most relevant to this analysis are the charged particle tracking system, the CsI electromagnetic calorimeter, the time-of-flight system, and the muon chambers. In CLEO II, the momenta of charged particles are measured in a tracking system consisting of a 6-layer straw tube chamber, 10-layer precision drift chamber, and 51-layer main drift chamber, all operating inside a 1.5 T solenoidal magnet. The main drift chamber also provides a measurement of the specific ionization loss, $`dE/dx`$, used for particle identification. For CLEO II.V, the innermost wire chamber was replaced with a three-layer silicon vertex detector . Muon identification system consists of proportional counters placed at various depths in the steel absorber.

The results of this search are based upon an integrated luminosity of 9.1 $`\mathrm{fb}^1`$ of $`e^+e^{}`$ data taken at the $`\mathrm{{\rm Y}}(4S)`$ energy and 4.4 $`\mathrm{fb}^1`$ recorded 60 MeV below the $`\mathrm{{\rm Y}}(4S)`$ energy. The simulated event samples used in this analysis were generated using GEANT-based  simulation of the CLEO detector response. Simulated events were processed in a similar manner as the data.

When making requirements on such kinematic variables as invariant mass or energy, we took advantage of well-understood track and photon-shower covariance matrices to calculate the expected resolution for each combination. Therefore we extensively used normalized variables, which allowed uniform candidate selection criteria to be used for the data collected with the CLEO II and CLEO II.V detector configurations.

The normalized invariant mass distributions for the $`J/\psi \mathrm{}^+\mathrm{}^{}`$ signal in data are shown in Fig. 2. We required the normalized invariant mass to be from $`10`$ to $`+3`$ (from $`4`$ to $`+3`$) for the $`J/\psi e^+e^{}`$($`J/\psi \mu ^+\mu ^{}`$) candidates. The resolution in the $`\mathrm{}^+\mathrm{}^{}`$ invariant mass is about 10 MeV$`/c^2`$. To improve the energy and momentum resolution of a $`J/\psi `$ candidate, we performed a fit constraining the mass of each $`J/\psi `$ candidate to the world average value .

Electron candidates were identified based on the ratio of the track momentum to the associated shower energy in the CsI calorimeter and specific ionization loss in the drift chamber. The internal bremsstrahlung in the $`J/\psi e^+e^{}`$ decay as well as the bremsstrahlung in the detector material produce a long radiative tail in the $`e^+e^{}`$ invariant mass distribution and impede efficient $`J/\psi e^+e^{}`$ detection. We recovered some of the bremsstrahlung photons by selecting the photon shower with the smallest opening angle with respect to the direction of the $`e^\pm `$ track evaluated at the interaction point, and then requiring this opening angle to be smaller than $`5^{}`$. The addition of the bremsstrahlung photons resulted in a relative increase of approximately 25% in the $`J/\psi e^+e^{}`$ reconstruction efficiency without adding more background.

For the $`J/\psi \mu ^+\mu ^{}`$ reconstruction, one of the muon candidates was required to penetrate the steel absorber to a depth greater than 3 nuclear interaction lengths. We relaxed the absorber penetration requirement for the second muon candidate if it was not expected to reach a muon chamber either because its energy was too low or because it pointed to a region of the detector not covered by the muon chambers. For these muon candidates we required the ionization signature in the CsI calorimeter to be consistent with that of a muon. Muons typically leave a narrow trail of ionization and deposit approximately 200 MeV of energy in the crystal calorimeter. Hadrons, on the other hand, quite often undergo a nuclear interaction in the CsI crystals that have a depth of 80% of a nuclear interaction length. Compared to imposing the absorber penetration requirement on both muon candidates, this procedure increased the $`J/\psi \mu ^+\mu ^{}`$ reconstruction efficiency by 20% with 80% increase of background.

We required that the charged kaon candidates have $`dE/dx`$ and, if available, time-of-flight measurements that lie within 3 standard deviations of the expected values.

If for the $`BJ/\psi \varphi K`$ decays we assume a uniform Dalitz distribution and isotropic decays of $`J/\psi `$ and $`\varphi `$ mesons, then the expected efficiency of the combined $`dE/dx`$ and time-of-flight selection is approximately 90% per kaon candidate. The $`dE/dx`$ measurements alone provide the $`K/\pi `$ separation of more than 4 standard deviations for 92% of the $`\varphi `$ daughter kaons and for 64% of the “bachelor” kaons from $`B`$ decay. We selected $`\varphi K^+K^{}`$ candidates by requiring the $`K^+K^{}`$ invariant mass to be within 10 MeV/$`c^2`$ of the $`\varphi `$ mass . We did not use the normalized $`K^+K^{}`$ invariant mass because the mass resolution (1.2 MeV/$`c^2`$) is smaller than the $`\varphi `$ width (4.4 MeV) .

The $`K_S^0`$ candidates were selected from pairs of tracks forming well-measured displaced vertices. The resolution in $`\pi ^+\pi ^{}`$ invariant mass is about 4 MeV$`/c^2`$. We required the absolute value of the normalized $`\pi ^+\pi ^{}`$ invariant mass to be less than 4, then we performed a fit constraining the mass of each $`K_S^0`$ candidate to the world average value .

The $`BJ/\psi \varphi K`$ candidates were selected by means of two observables. The first observable is the difference between the energy of the $`B`$ candidate and the beam energy $`\mathrm{\Delta }EE(J/\psi )+E(\varphi )+E(K)E_{\mathrm{beam}}`$. The resolution in $`\mathrm{\Delta }E`$ for the $`BJ/\psi \varphi K`$ candidates is approximately 6 MeV. The second observable is the beam-constrained $`B`$ mass $`M(B)\sqrt{E_{\mathrm{beam}}^2p^2(B)}`$, where $`p(B)`$ is the absolute value of the $`B`$ candidate momentum. The resolution in $`M(B)`$ for the $`BJ/\psi \varphi K`$ candidates is about 2.7 MeV/$`c^2`$; it is dominated by the beam energy spread. The distributions of the $`\mathrm{\Delta }E`$ vs $`M(B)`$ for $`B^+J/\psi \varphi K^+`$ and $`B^0J/\psi \varphi K_S^0`$ are shown in Fig. 3. We used the normalized $`\mathrm{\Delta }E`$ and $`M(B)`$ variables to select the $`BJ/\psi \varphi K`$ candidates and defined the signal region as $`|\mathrm{\Delta }E/\sigma (\mathrm{\Delta }E)|<3`$ and $`|(M(B)M_B)/\sigma (M(B))|<3`$. We observed 8(2) events in the signal region for the $`B^+J/\psi \varphi K^+`$ ($`B^0J/\psi \varphi K_S^0`$) mode. Considering that $`K^0`$ can decay as $`K_S^0`$ or as $`K_L^0`$, and also taking into account $`(K_S^0\pi ^+\pi ^{})`$ and the difference in reconstruction efficiencies, we expect to observe on average 4.3 $`B^+J/\psi \varphi K^+`$ candidates for every $`B^0J/\psi \varphi K_S^0`$ candidate.

The Dalitz plot and the cosine of helicity angle distributions for the 10 $`BJ/\psi \varphi K`$ signal candidates are shown in Figs. 4 and 5. The helicity angle for $`J/\psi \mathrm{}^+\mathrm{}^{}`$ decay is defined as the angle between a lepton momentum in the $`J/\psi `$ rest frame and the $`J/\psi `$ momentum in the $`B`$ rest frame. An analogous definition was used for the $`\varphi K^+K^{}`$ decay. No conclusion can be drawn yet either about the $`J/\psi `$ and the $`\varphi `$ polarizations or about the resonant substructure of the $`BJ/\psi \varphi K`$ decay. If the $`J/\psi `$ and $`\varphi `$ mesons are the products of the hybrid charmonium $`\psi _g`$ decay, then the $`J/\psi \varphi `$ invariant mass is expected to be below the $`DD^{}`$ threshold (4.3 GeV/$`c^2`$) because $`\psi _gDD^{}`$ decay is likely to dominate above the threshold . The $`J/\psi \varphi `$ invariant mass is above 4.3 GeV/$`c^2`$ for all 10 $`BJ/\psi \varphi K`$ candidates thus difavoring the hybrid charmonium dominance scenario.

The background can be divided into two categories. First category is the combinatorial background from $`\mathrm{{\rm Y}}(4S)B\overline{B}`$ and continuum non-$`B\overline{B}`$ events. Second category is the background from non-resonant $`BJ/\psi K^+K^{}K`$ decays.

The combinatorial background from $`\mathrm{{\rm Y}}(4S)B\overline{B}`$ events was estimated using a sample of simulated events approximately 32 times the data sample; events containing a $`BJ/\psi K^+K^{}K`$ decay were excluded. We estimated the background from $`\mathrm{{\rm Y}}(4S)B\overline{B}`$ decays to be $`0.25_{0.08}^{+0.10}`$ events. In addition, we specifically considered $`BJ/\psi K^{}\pi ^+`$ with $`K^{}K\pi ^{}`$ and $`BJ/\psi \rho ^0K`$ decays because the beam-constrained $`B`$ mass distribution for these modes is the same as for the $`BJ/\psi \varphi K`$ decays. Using data and simulated events, we verified that those backgrounds are rendered negligible by the kaon identification, $`\varphi `$ mass, and $`\mathrm{\Delta }E`$ requirements. The combinatorial background from the continuum non-$`B\overline{B}`$ events was estimated using simulated events and the data collected below $`B\overline{B}`$ threshold. We found the continuum background to be negligible.

To estimate the background contribution from the non-resonant $`BJ/\psi K^+K^{}K`$ decays, we reconstructed $`B^+J/\psi K^+K^{}K^+`$ and $`B^0J/\psi K^+K^{}K_S^0`$ candidates in data requiring $`|M(K^+K^{})M_\varphi |>20`$ MeV/$`c^2`$ to exclude $`BJ/\psi \varphi K`$ events. We observed 7 $`BJ/\psi K^+K^{}K`$ candidates with the estimated $`B\overline{B}`$ combinatorial background of 2.8 events. We estimated the mean background from $`BJ/\psi K^+K^{}K`$ decays for the $`BJ/\psi \varphi K`$ signal to be $`0.27_{0.17}^{+0.21}`$ events; we assumed that $`BJ/\psi K^+K^{}K`$ decays according to phase space.

In summary, the estimated total background for the combined $`BJ/\psi \varphi K`$ signal is $`0.52_{0.19}^{+0.23}`$ events.

We evaluated the reconstruction efficiency using a sample of simulated $`BJ/\psi \varphi K`$ decays. We assumed a uniform Dalitz distribution and isotropic decays of $`J/\psi `$ and $`\varphi `$ mesons; these assumptions are consistent with data (Figs. 4 and 5). The reconstruction efficiency, which does not include branching fractions of daughter particle decays, is $`(15.5\pm 0.2)\%`$ for the $`B^+J/\psi \varphi K^+`$ mode and $`(10.3\pm 0.2)\%`$ for the $`B^0J/\psi \varphi K_S^0`$ mode. The reconstruction efficiency is close to zero at the edges of phase space where either $`\varphi `$ or $`K`$ meson is produced nearly at rest in the laboratory frame. Thus, the overall detection efficiency would be much smaller than the above values if the $`BJ/\psi \varphi K`$ decay is dominated by either a $`J/\psi K`$ resonance with a mass around 4.3 GeV/$`c^2`$ or a $`J/\psi \varphi `$ resonance with a mass around 4.8 GeV/$`c^2`$. No such resonances are expected. To assign the systematic uncertainty due to the decay model dependence of the reconstruction efficiency, we generated two additional samples of simulated $`BJ/\psi \varphi K`$ events. One sample was generated with a uniform Dalitz distribution for $`BJ/\psi \varphi K`$ and $`100\%`$ transverse polarization for $`J/\psi `$ and $`\varphi `$. The other sample was generated assuming the $`\varphi `$ and $`K`$ mesons to be daughters of a hypothetical spin-0 resonance with mass 1.7 GeV/$`c^2`$ and width 100 MeV. We estimated the relative systematic uncertainty due to the decay model dependence of the reconstruction efficiency extraction to be 7%.

For the branching fraction calculation we assumed equal production of $`B^+B^{}`$ and $`B^0\overline{B}^0`$ pairs at the $`\mathrm{{\rm Y}}(4S)`$ resonance and $`(B^+J/\psi \varphi K^+)=(B^0J/\psi \varphi K^0)=(BJ/\psi \varphi K)`$. We did not assign any systematic uncertainty due to these two assumptions. We used the world average values of $`(J/\psi \mathrm{}^+\mathrm{}^{})`$, $`(\varphi K^+K^{})`$, and $`(K_S^0\pi ^+\pi ^{})`$ . We used the tables in Ref.  to assign the 68.27% C.L. intervals for the Poisson signal mean given the total number of events observed and the known mean background. The resulting branching fraction is $`(BJ/\psi \varphi K)=(8.8_{3.0}^{+3.5}[\mathrm{stat}]\pm 1.3[\mathrm{syst}])\times 10^5`$.

The systematic error includes the uncertainty in the reconstruction efficiency due to decay modeling plus the uncertainties in track finding, track fitting, lepton and charged-kaon identification, $`K_S^0`$ finding, background subtraction, uncertainty in the number of $`B\overline{B}`$ pairs used for this measurement, statistics of the simulated event samples, and the uncertainties on the daughter branching fractions $`(J/\psi \mathrm{}^+\mathrm{}^{})`$ and $`(\varphi K^+K^{})`$ . We estimated the total relative systematic uncertainty of the $`(BJ/\psi \varphi K)`$ measurement to be 15%.

In conclusion, we have fully reconstructed 10 $`BJ/\psi \varphi K`$ candidates with a total estimated background of 0.5 events. Assuming equal production of $`B^+B^{}`$ and $`B^0\overline{B}^0`$ pairs at the $`\mathrm{{\rm Y}}(4S)`$ resonance and $`(B^+J/\psi \varphi K^+)=(B^0J/\psi \varphi K^0)=(BJ/\psi \varphi K)`$, we have measured $`(BJ/\psi \varphi K)=(8.8_{3.0}^{+3.5}[\mathrm{stat}]\pm 1.3(\mathrm{syst}))\times 10^5`$. This is the first observed $`B`$ meson decay requiring the creation of an additional $`s\overline{s}`$ quark pair.

We gratefully acknowledge the effort of the CESR staff in providing us with excellent luminosity and running conditions. This work was supported by the National Science Foundation, the U.S. Department of Energy, the Research Corporation, the UTPA-Faculty Research Council Program, the Natural Sciences and Engineering Research Council of Canada, the A.P. Sloan Foundation, the Swiss National Science Foundation, and the Alexander von Humboldt Stiftung.