Observation of the Λb 0 →J/ψΛϕ decay in proton-proton collisions at s=13TeV

The observation of the Λ 0b → J/ ψ Λ φ decay is reported using proton-proton collision data collected at √ s = 13 TeV by the CMS experiment at the LHC in 2018, corresponding to an integrated luminosity of 60 fb − 1 . The ratio of the branching fractions B ( Λ 0b → J/ ψ Λ φ ) / B ( Λ 0b → ψ ( 2S ) Λ ) is measured to be ( 8.26 ± 0.90 (stat) ± 0.68 (syst) ± 0.11 ( B )) × 10 − 2 , where the ﬁrst uncertainty is statistical, the second is systematic, and the last uncertainty reﬂects the uncertainties in the world-average branching fractions of φ and ψ ( 2S ) decays to the reconstructed ﬁnal states.


Introduction
Studies of b baryon decays are of great importance for probing the dynamics of heavy-flavor decay processes. Since the observation of the lightest b baryon Λ 0 b by the UA1 Collaboration [1] at the CERN SppS, followed by extensive studies at the Fermilab Tevatron by the CDF [2-11] and D0 [12][13][14][15][16][17] Collaborations, the ATLAS, CMS, and LHCb experiments have accomplished numerous Λ 0 b baryon studies, made possible by the large production cross section of bb pairs at the CERN LHC. Among these studies are precision mass measurements of the ground and excited states [18,19], as well as lifetime and polarization measurements [20][21][22][23]. Most of these studies have been performed in the Λ 0 b → J/ψΛ decay channel. Recently, an observation of the Λ 0 b baryon decay to an excited charmonium state Λ 0 b → ψ(2S)Λ has been reported by the ATLAS Collaboration [24], while the LHCb Collaboration observed other, higher-multiplicity decays involving charmonium states [25,26]. Decays of the Λ 0 b baryon also proved to be a rich source of exotic spectroscopy, as has been demonstrated by the observation by LHCb [27,28] of new pentaquark states P c (4312) + , P c (4380) + , and P c (4450) + in the invariant mass distribution of the J/ψp system produced in the Λ 0 b → J/ψpK − decay. Further studies of the Λ 0 b baryon decay modes involving charmonium states may shed light on the strong interaction processes in hadronic decays of b baryons and on the production of exotic multiquark states.
The Λ 0 b → J/ψΛφ decay is expected to proceed via the b → ccs process, similarly to the Λ 0 b → J/ψΛ decay, but requires an additional ss pair. Consequently, the measurement of its branching fraction could enhance the understanding of the final-state strong interactions in b baryon decays and test heavy-quark effective theory [30]. In addition, the Λ 0 b → J/ψΛφ decay is a baryonic analog of the B + → J/ψφK + decay, where a rich resonant structure in the J/ψφ system has been observed by several experiments [31][32][33][34]. Therefore, detailed studies of the J/ψφ spectrum produced in baryonic decays may provide an important test for the production of these states. Recently, the existence of a hidden-charm pentaquark spectra was predicted for the J/ψΛ final state [35], which can be investigated in the Λ 0 b → J/ψΛφ decay, once a sufficient number of signal events is accumulated.

The CMS detector
The central feature of the CMS apparatus [36] is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter, each composed of a barrel and two endcap sections. Forward calorimeters extend the pseudorapidity (η) coverage provided by the barrel and endcap detectors. Muons are detected in gas-ionization chambers embedded in the steel flux-return yoke outside the solenoid. The main subdetectors used for the present analysis are the silicon tracker and the muon system.
The silicon tracker measures charged particles within the range |η| < 2.5. During the LHC running period when the data used in this Letter were recorded, the silicon tracker consisted of 1856 silicon pixel and 15 148 silicon strip detector modules. For nonisolated particles with transverse momentum 1 < p T < 10 GeV and |η| < 1.4, the track resolution is typically 1.5% in p T .
Muons are measured within |η| < 2.4, with detection planes made using three technologies: drift tubes, cathode strip chambers, and resistive-plate chambers. Tracks in the muon system are matched to those measured in the silicon tracker. The relative p T resolution is measured to be in the range 0.8-3.0% for muons with p T < 10 GeV used in this analysis, depending on the muon |η| [37].
Events of interest are selected using a two-tiered trigger system [38]. The first level (L1), composed of custom hardware processors, uses information from the calorimeters and muon detectors to select events at a rate up to 100 kHz within a fixed time interval of less than 4 µs. The second level, known as the high-level trigger (HLT), consists of a farm of processors running a version of the full event reconstruction software optimized for fast processing, and reduces the event rate to around 1 kHz before data storage.
A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [36].

Data sample and event selection
The analysis described in this Letter is based on a data sample of proton-proton collisions at a center-of-mass energy of 13 TeV, collected with the CMS detector in 2018 and corresponding to an integrated luminosity of 60 fb −1 .
Data were recorded with a dedicated trigger, optimized for the selection of b hadrons decaying to J/ψ(µ + µ − ) and two additional tracks from the charged hadrons emerging from the decay. The L1 trigger requires two oppositely charged muons with p T of at least 4 GeV, or two muons in the barrel region (|η| < 1.479) without any p T threshold. At the HLT, a J/ψ candidate decaying into a µ + µ − pair displaced from the interaction point is required, along with at least two tracks consistent with the displaced vertex. Each muon p T is required to be at least 4 GeV, while the dimuon p T is required to exceed 6.9 GeV. The J/ψ candidates reconstructed from dimuons are required to have an invariant mass between 2.9 and 3.3 GeV. The threedimensional distance of closest approach of the two muons to each other is required to be less than 0.5 cm. The fitted dimuon vertex is required to have a transverse decay length significance L xy (J/ψ)/σ L xy (J/ψ) > 3, where L xy (J/ψ) and σ L xy (J/ψ) are, respectively, the distance from the common vertex to the beam axis in the transverse plane and its uncertainty. Finally, the dimuon vertex fit probability, calculated using the χ 2 and the number of degrees of freedom of the fit, is required to exceed 10%, while the angle α between the dimuon p T vector and the direction connecting the beam axis and the dimuon vertex in the transverse plane is required to satisfy cos α > 0.9. Given the lack of a dedicated kaon identification, the two additional tracks are assigned a kaon mass hypothesis and required to have p T > 0.8 GeV, |η| < 2.5, and an invariant mass in a range of 0.95-1.30 GeV.
In the subsequent offline analysis we follow closely the selection of Ref. [39]. The p T threshold on the two muon candidates of 4 GeV and the requirement of |η| < 2.4 are kept. Two oppositely charged muon candidates are paired and required to originate from a common vertex. The vertex requirements applied at the HLT are confirmed in the offline selection. Also both muon candidates must match those that triggered the event readout. Dimuon candidates with an invariant mass within 100 MeV, which corresponds to approximately four effective widths, around the J/ψ meson mass M PDG J/ψ are selected (hereafter, M PDG X denotes the world-average mass of hadron X [29]), and the p T of the J/ψ meson is required to exceed 7 GeV.
To reconstruct a Λ 0 b candidate, the J/ψ candidate is combined with two oppositely charged, high-purity [40] tracks, assumed to be kaon candidates, and a Λ candidate. The p T of the tracks is required to exceed 0.8 GeV, and their invariant mass must satisfy 0.99 < M(K + K − ) < 1.05 GeV. The Λ candidates are formed from displaced two-prong vertices under the assumption of the Λ → pπ − decay, as described in Ref. [41]. Daughter particles of the Λ candidate are refitted to a common vertex with their invariant mass constrained to M PDG Λ , and the vertex fit probability is required to exceed 1%. The proton mass is assigned to the higher-momentum daughter track. To select the candidates in the Λ signal region, the following additional requirement is applied: The width of this window is chosen to correspond to approximately three times the effective width of the reconstructed Λ candidates. In addition, the Λ candidate is required to have a transverse momentum in excess of 1 GeV.
As the last step of the reconstruction, a fit to the common vertex of the Λ candidate, the two kaon tracks, and the dimuon pair is performed, with the dimuon mass constrained to M PDG J/ψ ; this vertex is referred to as the Λ 0 b vertex. The kinematic vertex fit probability of the Λ 0 b candidate is required to exceed 1%. The selected candidates are required to have p T (Λ 0 b ) > 10 GeV. Multiple proton-proton interactions in the same or nearby beam crossing (pileup) are present in the data, with an average multiplicity of 32, resulting in multiple reconstructed vertices in an event. The vertex with the lowest three-dimensional angle between the line connecting this vertex with the Λ 0 b vertex and the Λ 0 b candidate momentum is chosen as the primary vertex (PV). The following requirement is used to select Λ 0 b candidates consistent with originating from the PV: cos α(Λ 0 b , PV) > 0.99, where α(Λ 0 b , PV) is the two-dimensional angle in the transverse plane between the Λ 0 b candidate momentum and the vector pointing from the PV to the Λ 0 b vertex. The following requirement on the Λ 0 b vertex displacement is also applied: is the distance between the primary and Λ 0 b vertices in the transverse plane, and σ L xy (Λ 0 b ) is its uncertainty. Candidate decays for the normalization channel Λ 0 b → ψ(2S)Λ, with ψ(2S) → J/ψπ + π − , are selected using the same reconstruction chain. Identical requirements are used to select the J/ψ candidate, π + and π − tracks, and Λ candidate. An additional requirement is placed on the J/ψπ + π − invariant mass, |M(J/ψπ + π − ) − M PDG ψ(2S) | < 15 MeV, to select ψ(2S) candidates, where this window corresponds to approximately three effective widths of a reconstructed ψ(2S) candidate.
In case of multiple Λ 0 b candidates per event, the one with the highest vertex fit probability is chosen for both the signal and normalization channels. There are 18.9 and 7.4% of events with two or more reconstructed candidates for signal and normalization channels, respectively. When there are two or more candidates in an event, the MC simulation predicts that the correct candidate is chosen 84 ± 5 and 93 ± 13% of the time for the signal and normalization channels, respectively.
To calculate the reconstruction efficiency, a study based on simulated signal events for both channels is performed. The events are generated with PYTHIA 8.230 [42]. The Λ 0 b baryon decays are modeled with EVTGEN [43] v1.6.0 for both the Λ 0 b → J/ψΛφ and Λ 0 b → ψ(2S)Λ decay channels, following the three-body phase space model. The events are then passed through a detailed CMS detector simulation based on GEANT4 [44].

Signal yield extraction
The invariant mass distribution of the Λ 0 b → J/ψΛK + K − candidates selected using the strategy described in the previous section is shown in Fig. 1 (left). An unbinned, extended maximumlikelihood fit to a signal plus background hypothesis is performed on this observable and further mass distributions.
The signal is described by a double-Gaussian function with a floating common mean and total normalization, while the two widths and the relative fraction of the two Gaussian functions are fixed to the values obtained from simulation. The double-Gaussian function was chosen as a model that provides the best description of the simulated sample. The background is parameterized by a third-order Bernstein polynomial. The fit results in a signal yield of 380 ± 32 events. The signal significance is calculated to be 9.7 standard deviations in the asymptotic approximation [45], using the profile likelihood ratio of the signal plus background over the background-only hypothesis as the test statistic. Including modeling uncertainties in the signal and background shapes (described in Section 6) results in a reduction of the significance value to 9.4 standard deviations.
There is a bin with the yield significantly higher than the average background level in the left panel of Fig. 1, just below the signal Λ 0 b peak. The local significance of the excess is estimated to be less than three standard deviations. Several cross-checks have been performed to investigate this enhancement. The M(Λ 0 b → J/ψΛK + K − ) distribution with the requirement on the φ candidates to have a mass within 10 MeV of the nominal value shows no significant excess in this bin. The statistical significance of the Λ 0 b signal is 10.3 standard deviations in the asymptotic approximation. A statistically independent data set, collected in 2017, has been examined with the same selection, and no significant excess below the Λ 0 b peak was observed. As a result of these cross-checks, we attribute the excess to a statistical fluctuation.
An unbinned likelihood fit to the M(Λ 0 b → J/ψΛK + K − ) observable is employed to separate the signal and background components statistically, which is then used with the sPlot technique [46] to obtain the M(K + K − ) data distribution corresponding to signal Λ 0 b → J/ψΛK + K − decays. To extract the Λ 0 b → J/ψΛφ decay yield, the background-subtracted M(K + K − ) distribution is fitted with the convolution of a double-Gaussian and relativistic Breit-Wigner functions for the φ signal and a first-order Bernstein polynomial for the nonresonant component. The natural width of the φ meson is fixed to the world-average value [29]. It was checked that the natural width of the φ meson obtained from the fit when it was allowed to float was consistent with the world-average value within the uncertainties. Both widths and the relative fraction of the two Gaussians are fixed to the values obtained from fitting the simulated signal sample. The fit results in a signal yield of 286 ± 29 events. The M(K + K − ) invariant mass distribution, along with the result of the fit, are shown in Fig. 1 (right). represent the data and the curve is the result of the fit. The signal is described by a double-Gaussian function with floating common mean and total normalization, while the individual widths and the relative fraction of the two Gaussians are fixed from the fit to a simulated signal sample. The background is described by a third-order Bernstein polynomial function. The fit results in a signal yield of 884 ± 37 events. The non-ψ(2S) contribution in the Λ 0 b → J/ψπ + π − Λ signal was estimated to be negligible in the selected mass window |M(J/ψπ

Efficiency calculation
The Λ 0 b selection efficiencies in the signal and normalization channels are calculated as the ratio of the numbers of selected to generated events in simulated signal samples. The overall efficiency includes the trigger and reconstruction efficiencies and the detector acceptance. The efficiency in each channel is obtained using the simulated samples described in Section 3. The efficiency ratio, which is used in the branching fraction ratio measurement, is found to be (Λ 0 b → ψ(2S)Λ)/ (Λ 0 b → J/ψΛφ) = 0.363 ± 0.011, where the uncertainty is statistical only and accounts for the limited event counts in the corresponding simulated samples. The p T spectrum of pions from the ψ(2S) → J/ψπ + π − decay in the normalizations channel is softer than the p T spectrum of kaons from the φ → K + K − decay in the signal channel, resulting in an efficiency ratio significantly below unity.

Systematic uncertainties
In this section we discuss various sources of systematic uncertainty contributing to the measurement of the ratio B( decay modes have the same topology, the systematic uncertainties related to the muon and track reconstruction, as well as the trigger efficiency, mostly cancel in the ratio. To test this assumption, simulated samples were compared with background-subtracted data in a number of kinematic distributions. As a result of these studies, an additional systematic uncertainty is assigned to account for the observed difference between data and simulation in the Λ 0 b rapidity distribution for the normalization channel, as well as for the difference in the two-body invariant mass distributions M(J/ψΛ), M(J/ψφ), and M(Λφ) in data and simulation for the signal channel. The latter discrepancy could be caused by a deviation from the pure phase space decay model used in the simulation due to contributions from intermediate resonant states; however, the statistical power of the present data set is insufficient to perform a more detailed investigation. To estimate this systematic uncertainty, the simulated samples were reweighted to match the distributions observed in data. The difference in the efficiency ratio before and after the reweighting is taken as the corresponding systematic uncertainty.
The systematic uncertainty related to the choice of the background model is estimated separately for the signal channel, normalization channel, and φ → K + K − decays. The variation of the background model includes Bernstein polynomials of second and fourth orders, independently for the signal and normalization channels, and an exponential function for the background in the φ → K + K − invariant mass distribution. For the signal channel, an additional background function with a threshold behavior is also tested: (x − x 0 ) β multiplied by the Bernstein polynomials of first and second orders, where x 0 = M PDG J/ψ + M PDG Λ + M PDG φ and the exponent β is allowed to vary freely in the fit. In each case, the maximum deviation in the measured signal yield within the variations of the background model is used as the systematic uncertainty.
Another source of systematic uncertainty is the signal shape modeling in the M(J/ψΛK + K − ), M(ψ(2S)Λ), and M(K + K − ) distributions. This uncertainty is estimated by using alternative signal models whose parameters were obtained by fitting the simulated invariant mass distributions. The variation of signal models includes a triple-Gaussian function and a sum of two Crystal Ball [47] functions for the Λ 0 b → J/ψΛK + K − invariant mass distribution; a sum of two Crystal Ball functions for the Λ 0 b → ψ(2S)Λ channel; and a convolution of a double Crystal Ball [48] and relativistic Breit-Wigner functions for the M(K + K − ) distribution. For each of the variations, the largest deviation in the measured signal yield is taken as the systematic uncertainty.
The next source of systematic uncertainty is the difference in the mass resolution of the Λ 0 b and φ peaks between data and simulation. To estimate this uncertainty, several variations were applied to the resolution functions in the M(J/ψΛK + K − ) and M(ψ(2S)Λ) distributions: only the ratio of the two Gaussian widths was fixed to the one measured in simulation instead of fixing both widths, as in the nominal fit. For the M(K + K − ) distribution, a fit with the fixed ratios of the two Gaussian widths and yields, as measured in simulation, is performed. In each case, the maximum variation in the measured Λ 0 b yield is used as the systematic uncertainty. The difference between data and simulation in the measured Λ 0 b mass resolution for the Λ 0 b → J/ψΛK + K − channel results in the largest systematic uncertainty.
The statistical uncertainty in the efficiency ratio obtained from simulation is also considered as a source of systematic uncertainty. Table 1 summarizes the individual sources of the systematic uncertainty, as well as the overall uncertainty obtained as a quadratic sum of the individual components.

Source
Relative uncertainty (%) Data/simulation difference in the Total systematic uncertainty 8.2

Summary
The observation of the Λ 0 b → J/ψΛφ decay and the measurement of the branching fraction ratio B( is presented using a data sample of proton-proton collisions at √ s = 13 TeV collected in 2018 by the CMS experiment and corresponding to an integrated luminosity of 60 fb −1 . The ratio B( is measured to be (8.26 ± 0.90 (stat) ± 0.68 (syst) ± 0.11(B)) × 10 −2 , where the first uncertainty is statistical, the second is systematic, and the last uncertainty reflects the uncertainties in the world-average branching fractions of φ and ψ(2S) decays to the reconstructed final states. The observation of the Λ 0 b → J/ψΛφ decay opens a window on future searches for new resonances in the J/ψΛ and J/ψφ mass spectra, once a sufficient number of signal events is observed.
institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies:  [6] CDF Collaboration, "Measurement of σ [8] CDF Collaboration, "First measurement of the ratio of branching fractions [11] CDF Collaboration, "Measurement of the branching fraction B(Λ 0 b → Λ + c π − π + π − ) at CDF", Phys. Rev. D 85 (2012) 032003, doi:10.1103/PhysRevD.85.032003, arXiv:1112.3334.