Observation of the Z-> psi l ( + ) l (-) Decay in pp Collisions at root s = 13 TeV The CMS collaboration 2018-1004

BMWFW and FWF (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, FAPERGS, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and NSFC (China); COLCIENCIAS (Colombia); MSES and CSF (Croatia); RPF (Cyprus); SENESCYT (Ecuador); MoER, ERC IUT, and ERDF (Estonia); Academy of Finland, MEC, and HIP (Finland); CEA and CNRS/IN2P3 (France); BMBF, DFG, and HGF (Germany); GSRT (Greece); NKFIA (Hungary); DAE and DST (India); IPM (Iran); SFI (Ireland); INFN (Italy); MSIP and NRF (Republic of Korea); LAS (Lithuania); MOE and UM (Malaysia); BUAP, CINVESTAV, CONACYT, LNS, SEP, and UASLP-FAI (Mexico); MBIE (New Zealand); PAEC (Pakistan); MSHE and NSC (Poland); FCT (Portugal); JINR (Dubna); MON, RosAtom, RAS, RFBR, and NRC KI (Russia); MESTD (Serbia); SEIDI, CPAN, PCTI, and FEDER (Spain); Swiss Funding Agencies (Switzerland); MST (Taipei); ThEPCenter, IPST, STAR, and NSTDA (Thailand); TUBITAK and TAEK (Turkey); NASU and SFFR (Ukraine); STFC (United Kingdom); DOE and NSF (USA).

Although the Z boson was discovered more than 30 years ago [1], only one exclusive decay channel with leptons, Z → 4l [2-6], has been observed apart from the dilepton final states.For radiative dilepton decays, Z → l þ l − γ, where l ¼ e, μ, experiments have reported yields consistent with the standard model, as well as upper limits on the branching fraction for anomalous production [7].No resonant structure in the four-lepton decay has yet been observed.The high rate of Z boson production at the CERN LHC facilitates the study of rare decay channels such as Z → Vγ, Z → Vl þ l − , and Z → VV, where V is a vector meson with J PC ¼ 1 −− [8,9].In this paper, we present the observation of the decay of the Z boson to a final state with a J=ψ meson and two oppositely charged same-flavor leptons.
The Z → Vl þ l − process has been described and studied in various theoretical papers [10][11][12][13][14][15][16].For the case where V ¼ J=ψ, the branching fraction BðZ → J=ψl þ l − Þ is calculable within the standard model.The dominant diagram is the quantum electrodynamics radiative process illustrated in Fig. 1, with the γ Ã -V transition strength derived from the measured V → l þ l − electromagnetic decays [17].The theoretical estimates of the branching fraction cover the range ð6.7-7.7Þ× 10 −7 [10,11].Although this branching fraction is small, the dileptons and vector meson in the final state offer a clean signature.The measurement of this branching fraction is valuable for the calculation of the fragmentation function for a virtual photon to split into a J=ψ meson.Rare Higgs boson decays, such as those to quarkonia [18,19], will become accessible in the future, making it possible to search for nonstandard model signatures in these decays, including, e.g., anomalous couplings or new exotic light states [20].Accurate knowledge of potential backgrounds from Z decays to quarkonia will be essential for these measurements.
This analysis uses proton-proton (pp) collision data recorded by the CMS experiment at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb −1 .We report the observation of the Z → ψl þ l − decay channel, where ψ represents the contributions from direct J=ψ and J=ψ mesons from ψð2SÞ decays, and the J=ψ is detected via its μ þ μ − decay channel.We measure the ratio of the branching fraction of this decay to that of the Z → μ þ μ − μ þ μ − decay, to take advantage of a partial cancellation of systematic uncertainties.
The central feature of the CMS apparatus 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 (ECAL), and a brass and scintillator hadron calorimeter, each composed of a barrel and two endcap sections.Muons are detected in gasionization chambers embedded in the steel flux-return yoke outside the solenoid, in the pseudorapidity range jηj < 2.4 [21].Electrons are reconstructed using information from the ECAL and the tracker, in the jηj < 2.5 range [22].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. [23].
Events of the Z → μ þ μ − μ þ μ − process are simulated with the next-to-leading-order Monte Carlo (MC) generator POWHEG [24], interfaced with PYTHIA 8.175 [25,26] with parameters set by the CUETP8M1 tune [27] for parton showering, hadronization, and the underlying event.The parton distribution functions are taken from the NNPDF 3.0 [28] set.For the Z → J=ψl þ l − signal we use PYTHIA 8.175 (same tune as for Z → μ þ μ − μ þ μ − ) to simulate the production of Z bosons, with an unpolarized phase-space model for the Z → J=ψl þ l − decay.Matrix-element effects are evaluated by comparison with data and treated as systematic uncertainties.The detector response is simulated with a model of the CMS detector implemented in the GEANT4 package [29].We measure the fiducial branching fraction restricted to a region of phase space covered by the acceptance of the measurement, as described below.
The trigger and offline selection criteria closely follow the previous CMS analysis of Z → 4l decays [2][3][4].Triggers requiring one, two, or three charged leptons, with varying p T requirements, are used.The combined efficiency of the triggers, within the acceptance of this analysis defined below, is greater than 99%.
Among the multiple pp collisions within the time resolution of the data acquisition, the primary vertex is taken to be the reconstructed vertex with the largest sum of p 2 T over the physics objects in the event.These objects include jets, clustered using the anti-k T jet finding algorithm [30,31] with the tracks assigned to the vertex as inputs, and the associated missing transverse momentum, taken as the negative vector sum of the ⃗p T of those jets.The primary vertex is required to lie within 24 cm of the center of the detector along the beam axis and 2 cm perpendicular to that axis.Charged particle tracks associated with vertices other than the primary vertex are ignored.
We require all lepton candidate trajectories to pass within 1 (0.5) cm of the primary vertex in the direction along (perpendicular to) the beam axis.The lepton candidates from Z boson decay are required to be isolated from the hadronic activity in the event.To satisfy this requirement, the scalar sum of transverse energy deposits in the calorimeters and the p T of tracks is computed in a cone of radius ΔR ≡ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðΔηÞ 2 þ ðΔϕÞ 2 p ¼ 0.3 in η-ϕ around the lepton trajectory, where ϕ is the azimuthal angle in radians.The sum is corrected for other leptons from Z boson decay that fall within the isolation cone and for the average hadronic activity in an event.The ratio of this corrected sum to the lepton p T is required to be smaller than 0.35.Leptons are required to be separated by ΔR > 0.02ð0.05Þfor same-(different-)flavor pairs.
We select events with two oppositely charged reconstructed muons consistent with the dimuon decay of a J=ψ meson that, in combination with two additional oppositely charged electrons or muons (which we refer to as prompt leptons, l), are consistent with the decay Z → ψl þ l − .Specifically, the invariant mass of the ψ muon pair must satisfy 2.6 < m μ þ μ − < 3.6 GeV and that of the four leptons must satisfy [17].Each of the muons from J=ψ decay are required to have p T > 3.5 GeV and jηj < 2.4, and the p T of the J=ψ candidate must exceed 8.5 GeV.We require the highest-and second-highest-p T prompt leptons to have p T > 30 and 15 GeV, respectively, satisfy jηj < 2.5ð2.4Þ for l ¼ eðμÞ, and have a dilepton invariant mass m l þ l − < 80 GeV.The lepton p T thresholds ensure high trigger efficiency, and the invariant mass requirement suppresses the background from events in which a dilepton from Z boson decay is combined with a dimuon from an uncorrelated J=ψ decay or a nonresonant muon pair.
The four leptons, and separately the two muons from the J=ψ decay, are fitted to common vertices, with each vertex fit required to have a χ 2 probability greater than 5%.The significance of the three-dimensional impact parameter relative to the primary vertex is required to satisfy jd IP =σ IP j < 4 for each lepton, where d IP is the signed distance of closest approach of the lepton track to the primary vertex and σ IP is the associated uncertainty.
Following the application of the selection criteria described above, 29 (18) masses for the candidate events.The signal appears as a concentration of events in the overlap region of the J=ψ meson and Z boson masses.The events outside the central cluster along the Z boson mass band indicate contributions from the Z → ðcontinuum μ þ μ − Þl þ l − decay, and along the J=ψ meson mass band, nonresonant J=ψl þ l − production.
We measure the branching fraction of the here the required mass ranges of the two oppositely charged muon pairs are 4ð40Þ < mðμ þ μ − Þ < 80 GeV, where the 40 GeV threshold applies to the pair with the larger invariant mass.
The signal yield is obtained from unbinned extended maximum-likelihood fits [32] of the distributions in the two invariant mass variables m μ þ μ − and m μ þ μ − l þ l − , separately for the dimuon and dielectron channels.The probability density function (pdf) is a sum of four terms, each of which is a yield parameter multiplying a component pdf of the form with the mean fixed to the J=ψ meson mass [17] and the width as a free parameter of the fit.
with its central value and width fixed to the mass and width of the Z boson [17], convolved with a Gaussian function whose width is a free parameter.The pdfs for the continuum background in each dimension of the fit, representing backgrounds that are both peaking and nonpeaking in the orthogonal dimension, are exponential functions with free decay parameters.The projections in each variable are shown in Fig. 3, along with the pdf components resulting from the fits.
The yields resulting from the fit are 13.0 AE 3.9 events for the Z → ψμ þ μ − mode and 11.2 AE 3.4 events for Z → ψe þ e − , where the uncertainties are statistical only.The yields of the two decay modes agree within uncertainties, as expected, since the reconstruction efficiencies of the prompt electrons and muons in this p T range are similar.The Z → μ þ μ − μ þ μ − reference signal is extracted with a separate extended unbinned maximum-likelihood fit to the m μ þ μ − μ þ μ − distribution, using the same parametrization as for Z → ψl þ l − .The fit yields 250 AE 20 events.
We evaluate the signal significance for both ψμ þ μ − and ψe þ e − by generating random pseudoexperiments with dimuon and four-lepton invariant mass distributions drawn from the background-only pdf and then fitted with the background-only and signal-plus-background hypotheses.From the pseudoexperiments the likelihood ratio of the two hypotheses is calculated and compared with the likelihood ratio of the data.Taking into account the systematic uncertainties (discussed below), the background-only hypothesis is excluded at 4.0 and 4.3 standard deviations for ψμ þ μ − and ψe þ e − , respectively.The combination of the two significances based on the Fisher formalism [33] results in the observation of the Z → ψl þ l − decay mode with a significance of 5.7 standard deviations.
From the observed signal yield we compute a ratio of branching fractions defined over the fiducial phase space of the measurement defined in Table I.The entries consist   I. Definition of the fiducial phase space for the measurement of the ratio of branching fractions.Here, l refers to a prompt muon or electron from the signal decay, or to either of the two muons from the higher invariant-mass pair in the reference-channel decay, and μ refers to a J=ψ daughter or a member of the lower invariant-mass pair in the reference-channel decay.The symbol l 1 (l 2 ) refers to the prompt lepton having the higher (lower) value of p T .The p J=ψ T threshold is applied to the signal and the mðμ AE μ ∓ Þ requirement to the reference channel.

Fiducial requirement
40 < m l þ l − < 80 GeV jηðelectronsÞj < 2.5, jηðmuonsÞj < 2.4 p T ðl 1 ; l 2 ; μ; μÞ > ð30; 15; 3.5; 3.5Þ GeV Signal: p J=ψ T > 8.5 GeV Reference channel: 4 < mðμ AE μ ∓ Þ < 80 GeV of the kinematical requirements of the event selection given above, plus the additional requirement m l þ l − > 40 GeV for the Z → ψl þ l − candidates, which is added to match the selection of the Z → μ þ μ − μ þ μ − candidates and to avoid regions of the decay phase space in which the acceptance is steeply falling.This requirement removes 2 (0) events from the Z → ψe þ e − (Z → ψμ þ μ − ) sample, and 0.95 events from the fitted Z → ψe þ e − yield.The ratio of the fiducial branching fractions for lepton flavor l is i is the signal yield excluding the ψð2SÞ → J=ψX contribution, and The experimental efficiencies to reconstruct events within the fiducial phase space are determined from simulation; combined with the trigger efficiencies given above they are Calculated contributions from ψð2SÞ → J=ψX decays, the dominant feed-down source of J=ψ mesons, are subtracted from the signal yields, since the natural width of the Z boson does not allow the separation of the process ψð2SÞ → J=ψX from direct J=ψ production.The predicted production ratio of Z → J=ψl þ l − to Z → ψð2SÞl þ l − is 3.5 [11].Taking into account the branching fraction of ψð2SÞ to J=ψX [17], the ratio of NðZ Using this scale factor, we subtract 1.9 (1.7) events from the yield, considering them as J=ψ events from ψð2SÞ meson decays.
Since the signal and reference-channel events are recorded with the same triggers, and the topologies of the selected events are similar, many systematic uncertainties cancel in the ratio.The uncertainties in R J=ψl þ l − are shown for the two signal decay modes in columns 2 and 3 of Table II and are combined in quadrature as uncorrelated, unless stated otherwise, in column 4.
Systematic uncertainties arising from the choice of fit model are calculated by varying the pdfs used for the signal (Z and J=ψ) and combinatorial background.Substitution of a double-Gaussian function for the Z boson signal leads to differences in the signal yields of 0.02, 0.05, and 1.88 The corresponding changes from using a first-order polynomial instead of an exponential function for the Z boson combinatorial background are 0.9, 0.1, and 0.4 events.
A similar approach was followed for the J=ψ meson signal and background pdfs.The maximum difference observed in the signal yields resulting from the substitution of the sum of a double-Gaussian and a Crystal Ball [34] function for the signal pdf is 0.6 events for the ψμ þ μ − and 0.2 events for the ψe þ e − final state.The background pdf was replaced by a first-order polynomial to estimate the background model uncertainty, where a difference of 0.2 events is found in both decay modes.
To measure the uncertainty from the fitting procedure, 1000 random pseudosamples were generated with the number of events of each drawn from a Poisson distribution having a mean equal to the number of events observed in the data.The absolute value of the average deviation of the fit yields from the nominal yield is taken as the systematic uncertainty.
The reconstruction efficiencies of the muons from J=ψ decay and prompt leptons (electrons and muons) are checked with Z → μ þ μ − , Z → e þ e − , and J=ψ → μ þ μ − decay data using the "tag-and-probe" method [21,35], as functions of the lepton η and p T .To calculate the systematic uncertainty in R J=ψl þ l − , these efficiencies are varied within their uncertainty, with the uncertainties from the lepton efficiencies treated as correlated in the ratio.In addition, we assign an uncertainty associated with the finite number of MC signal and reference-channel events used to obtain the reconstruction efficiencies.
We test the three-body Z boson decay model implemented in the MC simulation by comparing distributions from the simulation with those from signal-weighted data, obtained from the fit model by the s Plot method [36].The most sensitive observables were found to be the azimuthal separation between the J=ψ candidate and the highest-and second-highest-p T prompt leptons.We apply the observed shape differences to the simulation and reevaluate the TABLE II.The contributions to the systematic uncertainty in the ratio of branching fractions for the prompt muon, prompt electron, and combined samples, in percent.The last row gives the sum in quadrature of all components.reconstruction efficiency to extract the decay model uncertainty.

Source of uncertainty
The uncertainty in the fraction of J=ψ events that potentially originate from ψð2SÞ is propagated from the uncertainty of the NðZ → The total systematic uncertainty of 7.6% for R J=ψl þ l − is calculated by adding the sources of uncertainty given in the last column of Table II in quadrature.
After subtracting the ψð2SÞ feed-down we extract from Eq. ( 1) the branching fraction ratio R J=ψl þ l − , for the phasespace region defined in Table I: Assuming that the factors applied to extrapolate the signal and reference-channel branching fractions from the phase space defined in Table I to the full phase space approximately cancel in the ratio, we use the measured . This estimate is consistent with standard model predictions of ð6.7 AE 0.7Þ × 10 −7 [10] and 7.7 × 10 −7 [11].
The factors that extrapolate the fiducial measurements to the full phase space depend on the Z boson decay matrix element, which determines the angular distributions of the muons coming from the ψ meson and the prompt leptons.Computing those factors assuming that the ψ is transversely or longitudinally polarized in the helicity frame (λ θ ¼ AE1) [37] leads to a full phase space branching fraction ratio that differs by less than 25% from the unpolarized result.
In summary, a new decay mode of the Z boson into a ψ meson, where ψ represents the contributions from direct J=ψ and ψð2SÞ → J=ψX, and an additional pair of leptons (muons or electrons), is observed with a statistical significance greater than 5 standard deviations.Using data from proton-proton collisions collected with the CMS detector at ffiffi ffi s p ¼ 13 TeV, corresponding to an integrated luminosity of 35.9 fb −1 , 13.0 AE 3.9 events of the Z → ψμ þ μ − and 11.2 AE 3.4 events of the Z → ψe þ e − decay are obtained.This is the first observed Z boson decay to a vector meson and two oppositely charged same-flavor leptons.
The ratio of the branching fraction for this decay to the one for the reference channel Z → μ þ μ − μ þ μ − in the fiducial phase space of the measurement, as defined in Table I, after subtracting the ψð2SÞ feed-down, is We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS 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: BMWFW and FWF (Austria);    B 273, 338 (1991).
FIG. 1. Leading-order Feynman diagram for the production of the Z boson and its decay in the Z → J=ψl þ l − channel.
FIG. 3. Invariant mass distributions for the ψ muon pairs (left) and for ψl þ l − (right), for Z → ψμ þ μ − (upper) and Z → ψe þ e − (lower) candidates.In each histogram the data are represented by the points, with the vertical bars showing the statistical uncertainties, and the solid curve is the overall fit to the data.The shaded region corresponds to the signal yield, while the long-dashed lines are the ψ meson signal from the Z boson background (left) and the Z boson signal from the ψ meson background (right).The shortdashed line represents the combinatorial background.
FIG. 2. Distribution of invariant massesm μ þ μ − vs m μ þ μ − l þ l −for the selected candidates.The values in the legend give the numbers of candidates per bin, which are also indicated by the sizes of the open black boxes.
The four terms account for the PHYSICALREVIEW LETTERS 121, 141801 (2018)