Measurements of two- and multi-particle angular correlations in pp collisions at sqrt(s) = 5, 7, and 13 TeV are presented as a function of charged-particle multiplicity. The data, corresponding to integrated luminosities of 1.0 inverse picobarn (5 TeV), 6.2 inverse picobarns (7 TeV), and 0.7 inverse picobarns (13 TeV), were collected using the CMS detector at the LHC. The second-order (v[2]) and third-order (v[3]) azimuthal anisotropy harmonics of unidentified charged particles, as well as v[2] of K0 short and Lambda/anti-Lambda particles, are extracted from long-range two-particle correlations as functions of particle multiplicity and transverse momentum. For high-multiplicity pp events, a mass ordering is observed for the v[2] values of charged hadrons (mostly pions), K0 short, and Lambda/anti-Lambda, with lighter particle species exhibiting a stronger azimuthal anisotropy signal below pt of about 2 GeV/c. For 13 TeV data, the v[2] signals are also extracted from four- and six-particle correlations for the first time in pp collisions, with comparable magnitude to those from two-particle correlations. These observations are similar to those seen in pPb and PbPb collisions, and support the interpretation of a collective origin for the observed long-range correlations in high-multiplicity pp collisions.

The study of (anti-)deuteron production in pp collisions has proven to be a powerful tool to investigate the formation mechanism of loosely bound states in high energy hadronic collisions. In this paper the production of (anti-)deuterons is studied as a function of the charged particle multiplicity in inelastic pp collisions at $\sqrt{s}=13$ TeV using the ALICE experiment. Thanks to the large number of accumulated minimum bias events, it has been possible to measure (anti-)deuteron production in pp collisions up to the same charged particle multiplicity ($\rm{d} N_{ch}/\rm{d}\eta\sim26$) as measured in p-Pb collisions at similar centre-of-mass energies. Within the uncertainties, the deuteron yield in pp collisions resembles the one in p-Pb interactions, suggesting a common formation mechanism behind the production of light nuclei in hadronic interactions. In this context the measurements are compared with the expectations of coalescence and Statistical Hadronisation Models (SHM).

Particle production sensitive to non-factorizable and non-perturbative processes that contribute to the underlying event associated with a high transverse momentum ($p_{T}$) jet in proton+proton collisions at $\sqrt{s}$=200 GeV is studied with the STAR detector. Each event is divided into three regions based on the azimuthal angle with respect to the highest-$p_{T}$ jet direction: in the leading jet direction ("Toward"), opposite to the leading jet ("Away"), and perpendicular to the leading jet ("Transverse"). In the Transverse region, the average charged particle density is found to be between 0.4 and 0.6 and the mean transverse momentum, $\langle p_{T}\rangle$, between 0.5-0.7 GeV/$c$ for particles with $p_{T}$$>$0.2 GeV/$c$ at mid-pseudorapidity ($|\eta|$$<$1) and jet $p_{T}$$>$15 GeV/$c$. Both average particle density and $\langle p_{T}\rangle$ depend weakly on the leading jet $p_{T}$. Closer inspection of the Transverse region hints that contributions to the underlying event from initial- and final-state radiation are significantly smaller in these collisions than at the higher energies, up to 13 TeV, recorded at the LHC. Underlying event measurements associated with a high-$p_{T}$ jet will contribute to our understanding of QCD processes at hard and soft scales at RHIC energies, as well as provide constraints to modeling of underlying event dynamics.

The production cross-sections of $\Upsilon(1S)$, $\Upsilon(2S)$ and $\Upsilon(3S)$ mesons in proton-proton collisions at $\sqrt{s}$= 13 TeV are measured with a data sample corresponding to an integrated luminosity of $277 \pm 11$ $\rm pb^{-1}$ recorded by the LHCb experiment in 2015. The $\Upsilon$ mesons are reconstructed in the decay mode $\Upsilon\to\mu^{+}\mu^{-}$. The differential production cross-sections times the dimuon branching fractions are measured as a function of the $\Upsilon$ transverse momentum, $p_{\rm T}$, and rapidity, $y$, over the range $0 < p_{\rm T}< 30$ GeV/c and $2.0 < y < 4.5$. The ratios of the cross-sections with respect to the LHCb measurement at $\sqrt{s}$= 8 TeV are also determined. The measurements are compared with theoretical predictions based on NRQCD.

The azimuthal anisotropy Fourier coefficients ($v_n$) in 8.16 TeV pPb data are extracted via long-range two-particle correlations as a function of event multiplicity and compared to corresponding results in pp and PbPb collisions. Using a four-particle cumulant technique, $v_n$ correlations are measured for the first time in pp and pPb collisions. The $v_2$ and $v_4$ coefficients are found to be positively correlated in all collision systems. For high multiplicity pPb collisions an anticorrelation of $v_2$ and $v_3$ is observed, with a similar correlation strength as in PbPb data at the same multiplicity. The new correlation results strengthen the case for a common origin of the collectivity seen in pPb and PbPb collisions in the measured multiplicity range.

Production cross-sections of prompt charm mesons are measured using data from $pp$ collisions at the LHC at a centre-of-mass energy of $5\,$TeV. The data sample corresponds to an integrated luminosity of $8.60\pm0.33\,$pb$^{-1}$ collected by the LHCb experiment. The production cross-sections of $D^0$, $D^+$, $D_s^+$, and $D^{*+}$ mesons are measured in bins of charm meson transverse momentum, $p_{\text{T}}$, and rapidity, $y$. They cover the rapidity range $2.0<y<4.5$ and transverse momentum ranges $0 < p_{\text{T}} < 10\, \text{GeV}/c$ for $D^0$ and $D^+$ and $1 < p_{\text{T}} < 10\, \text{GeV}/c$ for $D_s^+$ and $D^{*+}$ mesons. The inclusive cross-sections for the four mesons, including charge-conjugate states, within the range of $1 < p_{\text{T}} < 8\, \text{GeV}/c$ are determined to be \sigma(pp\rightarrow D^0 X) = 1004 \pm 3 \pm 54\,\mu\text{b} \sigma(pp\rightarrow D^+ X) = 402 \pm 2 \pm 30\,\mu\text{b} \sigma(pp\rightarrow D_s^+ X) = 170 \pm 4 \pm 16\,\mu\text{b} \sigma(pp\rightarrow D^{*+} X)= 421 \pm 5 \pm 36\,\mu\text{b} where the uncertainties are statistical and systematic, respectively.

A search for supersymmetry is presented based on proton-proton collision events containing identified hadronically decaying top quarks, no leptons, and an imbalance $p_\mathrm{T}^\text{miss}$ in transverse momentum. The data were collected with the CMS detector at the CERN LHC at a center-of-mass energy of 13 TeV, and correspond to an integrated luminosity of 35.9 fb$^{-1}$. Search regions are defined in terms of the multiplicity of bottom quark jet and top quark candidates, the $p_\mathrm{T}^\text{miss}$, the scalar sum of jet transverse momenta, and the $m_{\mathrm{T2}}$ mass variable. No statistically significant excess of events is observed relative to the expectation from the standard model. Lower limits on the masses of supersymmetric particles are determined at 95% confidence level in the context of simplified models with top quark production. For a model with direct top squark pair production followed by the decay of each top squark to a top quark and a neutralino, top squark masses up to 1020 GeV and neutralino masses up to 430 GeV are excluded. For a model with pair production of gluinos followed by the decay of each gluino to a top quark-antiquark pair and a neutralino, gluino masses up to 2040 GeV and neutralino masses up to 1150 GeV are excluded. These limits extend previous results.

Differential production cross sections of J/$\psi$ and $\psi$(2S) charmonium and $\Upsilon$(nS) (n = 1, 2, 3) bottomonium states are measured in proton-proton collisions at $\sqrt{s} =$ 13 TeV, with data collected by the CMS detector at the LHC, corresponding to an integrated luminosity of 2.3 fb$^{-1}$ for the J/$\psi$ and 2.7 fb$^{-1}$ for the other mesons. The five quarkonium states are reconstructed in the dimuon decay channel, for dimuon rapidity $|y| <$ 1.2. The double-differential cross sections for each state are measured as a function of $y$ and transverse momentum, and compared to theoretical expectations. In addition, ratios are presented of cross sections for prompt $\psi$(2S) to J/$\psi$, $\Upsilon$(2S) to $\Upsilon$(1S), and $\Upsilon$(3S) to $\Upsilon$(1S) production.

Measurements of the cross-section for producing \bquark quarks in the reaction $pp\to b\bar{b} X$ are reported in 7 and 13 TeV collisions at the LHC as a function of the pseudorapidity $\eta$ in the range $2<\eta<5$ covered by the acceptance of the LHCb experiment. The measurements are done using semileptonic decays of $b$-flavored hadrons decaying into a ground-state charmed hadron in association with a muon. The cross-sections in the covered $\eta$ range are $72.0\pm 0.3\pm6.8~\mu$b and $144\pm 1\pm 21~\mu$b for 7 and 13 TeV. The ratio is $2.00\pm0.02\pm0.26$, where the quoted uncertainties are statistical and systematic, respectively. The agreement with theoretical expectation is good at 7 TeV, but differs somewhat at 13 TeV. The measured ratio of cross-sections is larger at lower $\eta$ than the model prediction.

Using a data sample corresponding to an integrated luminosity of $2.0\,fb^{-1}$, collected by the LHCb experiment, the production of the $\eta_c(1S)$ state in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13 \text{ TeV}$ is studied in the rapidity range ${2.0 < y < 4.5}$ and in the transverse momentum range ${6.5 < p_{T} < 14.0\text{ GeV}}$. The cross-section for prompt production of $\eta_c(1S)$ mesons relative to that of the $J/\psi$ meson is measured using the ${p\bar{p}}$ decay mode and is found to be ${\sigma_{\eta_c(1S)}/\sigma_{J/\psi} = 1.69 \pm 0.15 \pm 0.10 \pm 0.18}$. The quoted uncertainties are, in order, statistical, systematic and due to uncertainties on the branching fractions of the ${J/\psi\to p \bar{p}}$ and ${\eta_c\to p \bar{p}}$ decays. The prompt $\eta_c(1S)$ production cross-section is determined to be ${\sigma_{\eta_c(1S)} = 1.26 \pm 0.11\pm 0.08 \pm 0.14 \,\mu b}$, where the last uncertainty includes that on the ${J/\psi}$ meson cross-section. The ratio of the branching fractions of $b$-hadron decays to the $\eta_c(1S)$ and ${J/\psi}$ states is measured to be ${\mathcal{B}_{b\to\eta_c X}/\mathcal{B}_{b\to J/\psi X} = 0.48 \pm 0.03 \pm 0.03 \pm 0.05}$, where the last uncertainty is due to those on the branching fractions of the ${J/\psi \to p \bar{p}}$ and ${\eta_c\to p \bar{p}}$ decays. The difference between the ${J/\psi}$ and $\eta_c(1S)$ masses is also determined to be ${113.0 \pm 0.7 \pm 0.1\text{ MeV}}$, which is the most precise single measurement of this quantity to date.

The WZ production cross section in proton-proton collisions at sqrt(s) = 13 TeV is measured with the CMS experiment at the LHC using a data sample corresponding to an integrated luminosity of 2.3 inverse femtobarns. The measurement is performed in the leptonic decay modes WZ to l nu l' l', where l, l'= e, mu. The measured cross section for the range 60 < m[l'l'] < 120 GeV is sigma(pp to WZ) = 39.9 +/- 3.2 (stat) +2.9/-3.1 (syst) +/- 0.4 (theo) +/- 1.3 (lumi) pb, consistent with the standard model prediction.

Bose-Einstein correlations of charged hadrons are measured over a broad multiplicity range, from a few particles up to about 250 reconstructed charged hadrons in proton-proton collisions at $\sqrt{s} =$ 13 TeV. The results are based on data collected using the CMS detector at the LHC during runs with a special low-pileup configuration. Three analysis techniques with different degrees of dependence on simulations are used to remove the non-Bose-Einstein background from the correlation functions. All three methods give consistent results. The measured lengths of homogeneity are studied as functions of particle multiplicity as well as average pair transverse momentum and mass. The results are compared with data from both CMS and ATLAS at $\sqrt{s} =$ 7 TeV, as well as with theoretical predictions.

A measurement of the W$^+$W$^-$ boson pair production cross section in proton-proton collisions at $\sqrt{s} =$ 13 TeV is presented. The data used in this study are collected with the CMS detector at the CERN LHC and correspond to an integrated luminosity of 35.9 fb$^{-1}$. The W$^+$W$^-$ candidate events are selected by requiring two oppositely charged leptons (electrons or muons). Two methods for reducing background contributions are employed. In the first one, a sequence of requirements on kinematic quantities is applied allowing a measurement of the total production cross section: 117.6 $\pm$ 6.8 pb, which agrees well with the theoretical prediction. Fiducial cross sections are also reported for events with zero or one jet, and the change in the zero-jet fiducial cross section with the jet transverse momentum threshold is measured. Normalized differential cross sections are reported within the fiducial region. A second method for suppressing background contributions employs two random forest classifiers. The analysis based on this method includes a measurement of the total production cross section and also a measurement of the normalized jet multiplicity distribution in W$^+$W$^-$ events. Finally, a dilepton invariant mass distribution is used to probe for physics beyond the standard model in the context of an effective field theory, and constraints on the presence of dimension-6 operators are derived.

The production of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13$ TeV is studied with the \lhcb detector. Cross-section measurements are performed as a function of the transverse momentum $p_\mathrm{T}$ and the rapidity $y$ of the $J/\psi$ meson in the region $p_\mathrm{T}<14\mathrm{GeV}/c$ and $2.0<y<4.5$, for both prompt $J/\psi$ mesons and $J/\psi$ mesons from $b$-hadron decays. The production cross-sections integrated over the kinematic coverage are $15.03\pm 0.03\pm 0.94$ $\mu$b for prompt $J/\psi$ and $2.25\pm 0.01\pm 0.14$ $\mu$b for $J/\psi$ from $b$-hadron decays, assuming zero polarization of the $J/\psi$ meson. The first uncertainties are statistical and the second systematic. The cross-section reported for $J/\psi$ mesons from $b$-hadron decays is used to extrapolate to a total $b\bar{b}$ cross-section. The ratios of the cross-sections with respect to $\sqrt{s}=8$ TeV are also determined.

Production cross-sections of prompt charm mesons are measured with the first data from $pp$ collisions at the LHC at a centre-of-mass energy of $13\,\mathrm{TeV}$. The data sample corresponds to an integrated luminosity of $4.98 \pm 0.19\,\mathrm{pb}^{-1}$ collected by the LHCb experiment. The production cross-sections of $D^{0}$, $D^{+}$, $D_{s}^{+}$, and $D^{*+}$ mesons are measured in bins of charm meson transverse momentum, $p_{\mathrm{T}}$, and rapidity, $y$, and cover the range $0 < p_{\mathrm{T}} < 15\,\mathrm{GeV}/c$ and $2.0 < y < 4.5$. The inclusive cross-sections for the four mesons, including charge conjugation, within the range of $1 < p_{\mathrm{T}} < 8\,\mathrm{GeV}/c$ are found to be \begin{equation} \sigma(pp \to D^{0} X) = 2072 \pm 2 \pm 124\,\mu\mathrm{b}\\ \sigma(pp \to D^{+} X) = 834 \pm 2 \pm \phantom{1}78\,\mu\mathrm{b}\\ \sigma(pp \to D_{s}^{+} X) = 353 \pm 9 \pm \phantom{1}76\,\mu\mathrm{b}\\ \sigma(pp \to D^{*+} X) = 784 \pm 4 \pm \phantom{1}87\,\mu\mathrm{b} \end{equation} where the uncertainties are due to statistical and systematic uncertainties, respectively.

A combination of four searches for new physics involving signatures with at least one photon and large missing transverse momentum, motivated by generalized models of gauge-mediated supersymmetry (SUSY) breaking, is presented. All searches make use of proton-proton collision data at $\sqrt{s}=$ 13 TeV, which were recorded with the CMS detector at the LHC in 2016, and correspond to an integrated luminosity of 35.9 fb$^{-1}$. Signatures with at least one photon and large missing transverse momentum are categorized into events with two isolated photons, events with a lepton and a photon, events with additional jets, and events with at least one high-energy photon. No excess of events is observed beyond expectations from standard model processes, and limits are set in the context of gauge-mediated SUSY. Compared to the individual searches, the combination extends the sensitivity to gauge-mediated SUSY in both electroweak and strong production scenarios by up to 100 GeV in neutralino and chargino masses, and yields the first CMS result combining various SUSY searches in events with photons at $\sqrt{s}=$ 13 TeV.

This article reports measurements characterizing the Underlying Event (UE) associated with hard scatterings at midrapidity in pp collisions at $\sqrt{s}=13$ TeV. The hard scatterings are identified by the leading particle, the charged particle with the highest transverse momentum ($p_{\rm T}^{\rm leading}$) in the event. Charged-particle number and summed transverse-momentum densities are measured in different azimuthal regions defined with respect to the leading particle direction: Toward, Transverse, and Away. The Toward and Away regions contain the fragmentation products of the hard scatterings in addition to the UE contribution, whereas particles in the Transverse region are expected to originate predominantly from the UE. The study is performed as a function of $p_{\rm T}^{\rm leading}$ with three different $p_{\rm T}$ thresholds for the associated particles, $p_{\rm T}^{\rm min} >$ 0.15, 0.5, and 1.0 GeV/$c$. The charged-particle density in the Transverse region rises steeply for low values of $p_{\rm T}^{\rm leading}$ and reaches a plateau. The results confirm the trend that the charged-particle density in the Transverse region shows a stronger increase with $\sqrt{s}$ than the inclusive charged-particle density at midrapidity. The UE activity is increased by approximately 20% when going from 7 to 13 TeV. The plateau in the Transverse region ($5 < p_{\rm T}^{\rm leading} < ~ 40$ GeV/$c$ ) is further characterized by the probability distribution of its charged-particle multiplicity normalized to its average value (relative transverse activity, $R_{T}$) and the mean transverse momentum as a function of $R_{T}$. Experimental results are compared to model calculations using PYTHIA 8 and EPOS LHC. The overall agreement between models and data is within 30%. These measurements provide new insights on the interplay between hard scatterings and the associated UE in pp collisions.

The production cross-section of $J/\psi$ pairs is measured using a data sample of $pp$ collisions collected by the LHCb experiment at a centre-of-mass energy of $\sqrt{s} = 13 \,{\mathrm{TeV}}$, corresponding to an integrated luminosity of $279 \pm 11 \,{\mathrm{pb^{-1}}}$. The measurement is performed for $J/\psi$ mesons with a transverse momentum of less than $10 \,{\mathrm{GeV}}/c$ in the rapidity range $2.0<y<4.5$. The production cross-section is measured to be $15.2 \pm 1.0 \pm 0.9 \,{\mathrm{nb}}$. The first uncertainty is statistical, and the second is systematic. The differential cross-sections as functions of several kinematic variables of the $J/\psi$ pair are measured and compared to theoretical predictions.

A measurement of the production cross-section of Z bosons in pp collisions at $\sqrt{s} = 13$ TeV is presented using dimuon and dielectron final states in LHCb data. The cross-section is measured for leptons with pseudorapidities in the range $2.0 < \eta < 4.5$, transverse momenta $p_\text{T} > 20$ GeV and dilepton invariant mass in the range $60<m(\ell\ell)<120$ GeV. The integrated cross-section from averaging the two final states is \begin{equation*}\sigma_{\text{Z}}^{\ell\ell} = 194.3 \pm 0.9 \pm 3.3 \pm 7.6\text{ pb,}\end{equation*} where the first uncertainty is statistical, the second is due to systematic effects, and the third is due to the luminosity determination. In addition, differential cross-sections are measured as functions of the Z boson rapidity, transverse momentum and the angular variable $\phi^*_\eta$.

The production cross-sections of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=5$ TeV are measured using a data sample corresponding to an integrated luminosity of $9.13\pm0.18~\text{pb}^{-1}$, collected by the LHCb experiment. The cross-sections are measured differentially as a function of transverse momentum, $p_{\text{T}}$, and rapidity, $y$, and separately for $J/\psi$ mesons produced promptly and from beauty hadron decays (nonprompt). With the assumption of unpolarised $J/\psi$ mesons, the production cross-sections integrated over the kinematic range $0<p_{\text{T}}<20~\text{GeV}/c$ and $2.0<y<4.5$ are $8.154\pm0.010\pm0.283~\mu\text{b}$ for prompt $J/\psi$ mesons and $0.820\pm0.003\pm0.034~\mu\text{b}$ for nonprompt $J/\psi$ mesons, where the first uncertainties are statistical and the second systematic. These cross-sections are compared with those at $\sqrt{s}=8$ TeV and $13$ TeV, and are used to update the measurement of the nuclear modification factor in proton-lead collisions for $J/\psi$ mesons at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\text{NN}}}=5$ TeV. The results are compared with theoretical predictions.

A search for supersymmetry is presented based on multijet events with large missing transverse momentum produced in proton-proton collisions at a center-of-mass energy of sqrt(s) = 13 TeV. The data, corresponding to an integrated luminosity of 35.9 inverse femtobarns, were collected with the CMS detector at the CERN LHC in 2016. The analysis utilizes four-dimensional exclusive search regions defined in terms of the number of jets, the number of tagged bottom quark jets, the scalar sum of jet transverse momenta, and the magnitude of the vector sum of jet transverse momenta. No evidence for a significant excess of events is observed relative to the expectation from the standard model. Limits on the cross sections for the pair production of gluinos and squarks are derived in the context of simplified models. Assuming the lightest supersymmetric particle to be a weakly interacting neutralino, 95% confidence level lower limits on the gluino mass as large as 1800 to 1960 GeV are derived, and on the squark mass as large as 960 to 1390 GeV, depending on the production and decay scenario.

Three of the most significant measured deviations from standard model predictions, the enhanced decay rate for B $\to$ D$^{(*)}\tau\nu$, hints of lepton universality violation in B $\to$ K$^{(*)}\ell\ell$ decays, and the anomalous magnetic moment of the muon, can be explained by the existence of leptoquarks (LQs) with large couplings to third-generation quarks and masses at the TeV scale. The existence of these states can be probed at the LHC in high energy proton-proton collisions. A novel search is presented for pair production of LQs coupled to a top quark and a muon using data at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$, recorded by the CMS experiment. No deviation from the standard model prediction has been observed and scalar LQs decaying exclusively into t$\mu$ are excluded up to masses of 1420 GeV. The results of this search are combined with those from previous searches for LQ decays into t$\tau$ and b$\nu$, which excluded scalar LQs below masses of 900 and 1080 GeV. Vector LQs are excluded up to masses of 1190 GeV for all possible combinations of branching fractions to t$\mu$, t$\tau$ and b$\nu$. With this analysis, all relevant couplings of LQs with an electric charge of -1/3 to third-generation quarks are probed for the first time.

Transverse momentum spectra of charged pions, kaons, and protons are measured in proton-proton collisions at sqrt(s) = 13 TeV with the CMS detector at the LHC. The particles, identified via their energy loss in the silicon tracker, are measured in the transverse momentum range of pt ~ 0.1-1.7 GeV and rapidities abs(y) < 1. The pt spectra and integrated yields are compared to previous results at smaller sqrt(s) and to predictions of Monte Carlo event generators. The average pt increases with particle mass and charged particle multiplicity of the event. Comparisons with previous CMS results at sqrt(s) = 0.9, 2.76, and 7 TeV show that the average pt and the ratios of hadron yields feature very similar dependences on the particle multiplicity in the event, independently of the center-of-mass energy of the pp collision.

A search for a heavy resonance decaying into a top quark and antiquark ($\mathrm{t\bar{t}}$) pair is performed using proton-proton collisions at $\sqrt{s} =$ 13 TeV. The search uses the data set collected with the CMS detector in 2016, which corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The analysis considers three exclusive final states and uses reconstruction techniques that are optimized for top quarks with high Lorentz boosts, which requires the use of nonisolated leptons and jet substructure techniques. No significant excess of events relative to the expected yield from standard model processes is observed. Upper limits on the production cross section of heavy resonances decaying to a $\mathrm{t\bar{t}}$ pair are calculated. Limits are derived for a leptophobic topcolor Z' resonance with widths of 1, 10, and 30%, relative to the mass of the resonance, and exclude masses up to 3.80, 5.25, and 6.65 TeV, respectively. Kaluza-Klein excitations of the gluon in the Randall-Sundrum model are excluded up to 4.55 TeV. To date, these are the most stringent limits on $\mathrm{t\bar{t}}$ resonances.

This Letter presents the observation of the rare Z boson decay Z $\to\psi\ell^+\ell^-$. Here, $\psi$ represents contributions from direct J/$\psi$ and $\psi$(2S) $\to$ J/$\psi X$, $\ell^+\ell^-$ is a pair of electrons or muons, and the J/$\psi$ meson is detected via its decay to $\mu^+\mu^-$. The sample of proton-proton collision data, collected by the CMS experiment at the LHC at a center-of-mass energy of 13 TeV, corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The signal is observed with a significance in excess of 5 standard deviations. After subtraction of the $\psi$(2S) $\to$ J/$\psi X$ contribution, the ratio of the branching fraction of the exclusive decay Z $\to\psi\ell^+\ell^-$ to the decay Z $\to\mu^+\mu^-\mu^+\mu^-$ within a fiducial phase space is measured to be $\mathcal{B}($Z $\to\psi\ell^+\ell^-) / \mathcal{B}($Z $\to\mu^+\mu^-\mu^+\mu^-) =$ 0.67 $\pm$ 0.18 (stat) $\pm$ 0.05 (syst).