A measurement is presented of the primary Lund jet plane (LJP) density in inclusive jet production in proton-proton collisions. The analysis uses 138 fb$^{-1}$ of data collected by the CMS experiment at $\sqrt{s}$ = 13 TeV. The LJP, a representation of the phase space of emissions inside jets, is constructed using iterative jet declustering. The transverse momentum $k_\mathrm{T}$ and the splitting angle $\Delta R$ of an emission relative to its emitter are measured at each step of the jet declustering process. The average density of emissions as function of $\ln(k_\mathrm{T}$/GeV) and $\ln(R/\Delta R)$ is measured for jets with distance parameters $R$ = 0.4 or 0.8, transverse momentum $p_\mathrm{T} \gt$ 700 GeV, and rapidity $\vert y\vert \lt $ 1.7. The jet substructure is measured using the charged-particle tracks of the jet. The measured distributions, unfolded to the level of stable particles, are compared with theoretical predictions from simulations and with perturbative quantum chromodynamics calculations. Due to the ability of the LJP to factorize physical effects, these measurements can be used to improve different aspects of the physics modeling in event generators.

A search for dark matter in events with a displaced nonresonant muon pair and missing transverse momentum is presented. The analysis is performed using an integrated luminosity of 138 fb$^{-1}$ of proton-proton (pp) collision data at a center-of-mass energy of 13 TeV produced by the LHC in 2016-2018. No significant excess over the predicted backgrounds is observed. Upper limits are set on the product of the inelastic dark matter production cross section $\sigma$(pp $\to$ A' $\to$$\chi_1$$\chi_2$) and the decay branching fraction $\mathcal{B}$($\chi_2$$\to$$\chi_1 \mu^+ \mu^-$), where A' is a dark photon and $\chi_1$ and $\chi_2$ are states in the dark sector with near mass degeneracy. This is the first dedicated collider search for inelastic dark matter.

An inclusive search for long-lived exotic particles (LLPs) decaying to final states with a pair of muons is presented. The search uses data corresponding to an integrated luminosity of 36.6 fb$^{-1}$ collected by the CMS experiment from the proton-proton collisions at $\sqrt{s}$ = 13.6 TeV in 2022, the first year of Run 3 of the CERN LHC. The experimental signature is a pair of oppositely charged muons originating from a common vertex spatially separated from the proton-proton interaction point by distances ranging from several hundred $\mu$m to several meters. The sensitivity of the search benefits from new triggers for displaced dimuons developed for Run 3. The results are interpreted in the framework of the hidden Abelian Higgs model, in which the Higgs boson decays to a pair of long-lived dark photons, and of an $R$-parity violating supersymmetry model, in which long-lived neutralinos decay to a pair of muons and a neutrino. The limits set on these models are the most stringent to date in wide regions of lifetimes for LLPs with masses larger than 10 GeV.

The J/$\psi$$\to$$\mu^+\mu^-\mu^+\mu^-$ decay has been observed with a statistical significance in excess of five standard deviations. The analysis is based on an event sample of proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS experiment in 2018 and corresponding to an integrated luminosity of 33.6 fb${-1}$. Normalizing to the J/$\psi$$\to$$\mu^+\mu^-$ decay mode leads to a branching fraction [10.1$^{+3.3}_{-2.7}$ (stat) $\pm$ 0.4 (syst) ]$\times$ 10$^{-7}$, a value that is consistent with the standard model prediction.

Using proton-proton collision data corresponding to an integrated luminosity of 140 fb$^{-1}$ collected by the CMS experiment at $\sqrt{s}$ = 13 TeV, the $\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$ decay is observed for the first time, with a statistical significance exceeding 5 standard deviations. The relative branching fraction, with respect to the $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$ decay, is measured to be $\mathcal{B}$($\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$)/$\mathcal{B}$( $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$) = [3.38 $\pm$ 1.02 $\pm$ 0.61 $\pm$ 0.03]%, where the first uncertainty is statistical, the second is systematic, and the third is related to the uncertainties in $\mathcal{B}$($\psi$(2S) $\to$ J/$\psi\pi^+\pi^-$) and $\mathcal{B}$($\Xi^-$ $\to$ $\Lambda\pi^-$).

A combination of the results of several searches for the electroweak production of the supersymmetric partners of standard model bosons, and of charged leptons, is presented. All searches use proton-proton collision data at $\sqrt{s}$ = 13 TeV recorded with the CMS detector at the LHC in 2016-2018. The analyzed data correspond to an integrated luminosity of up to 137 fb$^{-1}$. The results are interpreted in terms of simplified models of supersymmetry. Two new interpretations are added with this combination: a model spectrum with the bino as the lightest supersymmetric particle together with mass-degenerate higgsinos decaying to the bino and a standard model boson, and the compressed-spectrum region of a previously studied model of slepton pair production. Improved analysis techniques are employed to optimize sensitivity for the compressed spectra in the wino and slepton pair production models. The results are consistent with expectations from the standard model. The combination provides a more comprehensive coverage of the model parameter space than the individual searches, extending the exclusion by up to 125 GeV, and also targets some of the intermediate gaps in the mass coverage.

A search for long-lived heavy neutral leptons (HNLs) is presented, which considers the hadronic final state and coupling scenarios involving all three lepton generations in the 2-20 GeV HNL mass range for the first time. Events comprising two leptons (electrons or muons) and jets are analyzed in a data sample of proton-proton collisions, recorded with the CMS experiment at the CERN LHC at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. A novel jet tagger, based on a deep neural network, has been developed to identify jets from an HNL decay using various features of the jet and its constituent particles. The network output can be used as a powerful discriminating tool to probe a broad range of HNL lifetimes and masses. Contributions from background processes are determined from data. No excess of events in data over the expected background is observed. Upper limits on the HNL production cross section are derived as functions of the HNL mass and the three coupling strengths $V_{\ell\mathrm{N}}$ to each lepton generation $\ell$ and presented as exclusion limits in the coupling-mass plane, as lower limits on the HNL lifetime, and on the HNL mass. In this search, the most stringent limit on the coupling strength is obtained for pure muon coupling scenarios; values of $\lvert V_{\mu\mathrm{N}}\rvert^{2}$$\gt $ 5 (4) $\times$ 10$^{-7}$ are excluded for Dirac (Majorana) HNLs with a mass of 10 GeV at a confidence level of 95% that correspond to proper decay lengths of 17 (10) mm.

A search for W' bosons decaying to a top and a bottom quark in final states including an electron or a muon is performed with the CMS detector at the LHC. The analyzed data correspond to an integrated luminosity of 138 fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of 13 Tev. Good agreement with the standard model expectation is observed and no evidence for the existence of the W' boson is found over the mass range examined. The largest observed deviation from the standard model expectation is found for a W' boson mass ($m_\mathrm{W'}$) hypothesis of 3.8 TeV with a relative decay width of 1%, with a local (global) significance of 2.6 (2.0) standard deviations. Upper limits on the production cross sections of W' bosons decaying to a top and a bottom quark are set. Left- and right-handed W' bosons with $m_\mathrm{W'}$ below 3.9 and 4.3 TeV, respectively, are excluded at the 95% confidence level, under the assumption that the new particle has a narrow decay width. Limits are also set for relative decay widths up to 30%. These are the most stringent limits to date on this W' boson decay channel.

A search for the production of long-lived particles in proton-proton collisions at a center-of-mass energy of 13 TeV at the CERN LHC is presented. The search is based on data collected by the CMS experiment in 2016-2018, corresponding to a total integrated luminosity of 137 fb$^{-1}$. This search is designed to be sensitive to long-lived particles with mean proper decay lengths between 0.1 and 1000 $\mu$m, whose decay products produce a final state with at least one displaced vertex and missing transverse momentum. A machine learning algorithm, which improves the background rejection power by more than an order of magnitude, is applied to improve the sensitivity. The observation is consistent with the standard model background prediction, and the results are used to constrain split supersymmetry (SUSY) and gauge-mediated SUSY breaking models with different gluino mean proper decay lengths and masses. This search is the first CMS search that shows sensitivity to hadronically decaying long-lived particles from signals with mass differences between the gluino and neutralino below 100 GeV. It sets the most stringent limits to date for split-SUSY models and gauge-mediated SUSY breaking models with gluino proper decay length less than 6 $\mu$m.

The inclusive jet cross section is measured as a function of jet transverse momentum $p_\mathrm{T}$ and rapidity $y$. The measurement is performed using proton-proton collision data at $\sqrt{s}$ = 5.02 TeV, recorded by the CMS experiment at the LHC, corresponding to an integrated luminosity of 27.4 pb$^{-1}$. The jets are reconstructed with the anti-$k_\mathrm{T}$ algorithm using a distance parameter of $R$ = 0.4, within the rapidity interval $\lvert y\rvert$$\lt$ 2, and across the kinematic range 0.06 $\lt$$p_\mathrm{T}$$\lt$ 1 TeV. The jet cross section is unfolded from detector to particle level using the determined jet response and resolution. The results are compared to predictions of perturbative quantum chromodynamics, calculated at both next-to-leading order and next-to-next-to-leading order. The predictions are corrected for nonperturbative effects, and presented for a variety of parton distribution functions and choices of the renormalization/factorization scales and the strong coupling $\alpha_\mathrm{S}$.

A measurement of the dijet production cross section is reported based on proton-proton collision data collected in 2016 at $\sqrt{s}$ = 13 TeV by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of up to 36.3 fb$^{-1}$. Jets are reconstructed with the anti-$k_\mathrm{T}$ algorithm for distance parameters of $R$ = 0.4 and 0.8. Cross sections are measured double-differentially (2D) as a function of the largest absolute rapidity $\lvert y_\text{max}\rvert$ of the two jets with the highest transverse momenta $p_\mathrm{T}$ and their invariant mass $m_{1,2}$, and triple-differentially (3D) as a function of the rapidity separation $y^*$, the total boost $y_\mathrm{b}$, and either $m_{1,2}$ or the average $p_\mathrm{T}$ of the two jets. The cross sections are unfolded to correct for detector effects and are compared with fixed-order calculations derived at next-to-next-to-leading order in perturbative quantum chromodynamics. The impact of the measurements on the parton distribution functions and the strong coupling constant at the mass of the Z boson is investigated, yielding a value of $\alpha_\mathrm{S}$ = 0.1179 $\pm$ 0.0019.

A search for new physics in top quark production with additional final-state leptons is performed using data collected by the CMS experiment in proton-proton collisions at $\sqrt{s}$ = 13 TeV at the LHC during 2016-2018. The data set corresponds to an integrated luminosity of 138 fb$^{-1}$. Using the framework of effective field theory (EFT), potential new physics effects are parametrized in terms of 26 dimension-six EFT operators. The impacts of EFT operators are incorporated through the event-level reweighting of Monte Carlo simulations, which allows for detector-level predictions. The events are divided into several categories based on lepton multiplicity, total lepton charge, jet multiplicity, and b-tagged jet multiplicity. Kinematic variables corresponding to the transverse momentum ($p_\mathrm{T}$) of the leading pair of leptons and/or jets as well as the $p_\mathrm{T}$ of on-shell Z bosons are used to extract the 95% confidence intervals of the 26 Wilson coefficients corresponding to these EFT operators. No significant deviation with respect to the standard model prediction is found.

A search is presented for new Higgs bosons in proton-proton (pp) collision events in which a same-sign top quark pair is produced in association with a jet, via the pp $\to$ tH/A $\to$ t$\mathrm{\bar{t}}$c and pp $\to$ tH/A $\to$ t$\mathrm{\bar{t}}$u processes. Here, H and A represent the extra scalar and pseudoscalar boson, respectively, of the second Higgs doublet in the generalized two-Higgs-doublet model (g2HDM). The search is based on pp collision data collected at a center-of-mass energy of 13 TeV with the CMS detector at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. Final states with a same-sign lepton pair in association with jets and missing transverse momentum are considered. New Higgs bosons in the 200-1000 GeV mass range and new Yukawa couplings between 0.1 and 1.0 are targeted in the search, for scenarios in which either H or A appear alone, or in which they coexist and interfere. No significant excess above the standard model prediction is observed. Exclusion limits are derived in the context of the g2HDM.

Energy correlators that describe energy-weighted distances between two or three particles in a jet are measured using an event sample of $\sqrt{s}$ = 13 TeV proton-proton collisions collected by the CMS experiment and corresponding to an integrated luminosity of 36.3 fb$^{-1}$. The measured distributions reveal two key features of the strong interaction: confinement and asymptotic freedom. By comparing the ratio of the two measured distributions with theoretical calculations that resum collinear emissions at approximate next-to-next-to-leading logarithmic accuracy matched to a next-to-leading order calculation, the strong coupling is determined at the Z boson mass: $\alpha_\mathrm{S}(m_\mathrm{Z})$ = 0.1229$^{+0.0040}_{-0.0050}$, the most precise $\alpha_\mathrm{S}(m_\mathrm{Z})$ value obtained using jet substructure observables.

A search for partonic collective effects inside jets produced in proton-proton collisions is performed via correlation measurements of charged constituents using the CMS detector at the CERN LHC. The analysis uses data collected at a center-of-mass energy of $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. Jets are reconstructed with the anti-$k_\mathrm{T}$ algorithm with a distance parameter of 0.8 and are required to have transverse momentum greater than 550 GeV and pseudorapidity $\lvert\eta\rvert$$\lt$ 1.6. Two-particle correlations among the charged constituents within the jets are studied as functions of the particles' azimuthal angle and pseudorapidity separations ($\Delta\phi^*$ and $\Delta\eta^*$) in a jet coordinate basis, where constituents' $\eta^*$, $\phi^*$ are defined relative to the direction of the jet. The correlation functions are studied in classes of in-jet charged-particle multiplicity up to $N_\text{ch}^\mathrm{j}$$\approx$ 100. Fourier harmonics are extracted from long-range azimuthal correlation functions to characterize azimuthal anisotropy for $\lvert\Delta\eta^*\rvert$$\gt$ 2. For low-multiplicity jets, the long-range elliptic anisotropic harmonic, $v^*_2$, is observed to decrease with $N_\text{ch}^\mathrm{j}$. This trend is well described by Monte Carlo event generators. However, a rising trend for $v^*_2$ emerges at $N_\text{ch}^\mathrm{j}$$\gtrsim$ 80, hinting at a possible onset of collective behavior, which is not reproduced by the models tested. This observation yields new insights into the dynamics of parton fragmentation processes in the vacuum.

A measurement of the production cross section of a Higgs boson with transverse momentum greater than 250 GeV is presented where the Higgs boson decays to a pair of $\tau$ leptons. It is based on proton-proton collision data collected by the CMS experiment at the CERN LHC at a center-of-mass energy of 13 TeV. The data sample corresponds to an integrated luminosity of 138 fb$^{-1}$. Because of the large transverse momentum of the Higgs boson the $\tau$ leptons from its decays are boosted and produced spatially close, with their decay products overlapping. Therefore, a dedicated algorithm was developed to reconstruct and identify them. The observed (expected) significance of the measured signal with respect to the standard model background-only hypothesis is 3.5 (2.2) standard deviations. The product of the production cross section and branching fraction is measured to be 1.64$^{+0.68}_{-0.54}$ times the standard model expectation. The fiducial differential production cross section is also measured as functions of the Higgs boson and leading jet transverse momenta. This measurement extends the probed large-transverse-momentum region beyond 600 GeV.

A search for ZZ and ZH production in the $\mathrm{b\bar{b}b\bar{b}}$ final state is presented, where H is the standard model (SM) Higgs boson. The search uses an event sample of proton-proton collisions corresponding to an integrated luminosity of 133 fb$^{-1}$ collected at a center-of-mass energy of 13 TeV with the CMS detector at the CERN LHC. The analysis introduces several novel techniques for deriving and validating a multi-dimensional background model based on control samples in data. A multiclass multivariate classifier customized for the $\mathrm{b\bar{b}b\bar{b}}$ final state is developed to derive the background model and extract the signal. The data are found to be consistent, within uncertainties, with the SM predictions. The observed (expected) upper limits at 95% confidence level are found to be 3.8 (3.8) and 5.0 (2.9) times the SM prediction for the ZZ and ZH production cross sections, respectively.

A search for long-lived heavy neutrinos (N) in the decays of \PB mesons produced in proton-proton collisions at $\sqrt{s}$ = 13 TeV is presented. The data sample corresponds to an integrated luminosity of 41.6 fb$^{-1}$ collected in 2018 by the CMS experiment at the CERN LHC, using a dedicated data stream that enhances the number of recorded events containing B mesons. The search probes heavy neutrinos with masses in the range 1 $\lt$$m_\mathrm{N}$$\lt$ 3 GeV and decay lengths in the range 10$^{-2}$$\lt$$c\tau$$\lt$ 10$^{4}$ mm, where $\tau_\mathrm{N}$ is the N proper mean lifetime. Signal events are defined by the signature B $\to$$\ell_\mathrm{B}$NX; N $\to$$\ell^{\pm} \pi^{\mp}$, where the leptons $\ell_\mathrm{B}$ and $\ell$ can be either a muon or an electron, provided that at least one of them is a muon. The hadronic recoil system, X, is treated inclusively and is not reconstructed. No significant excess of events over the standard model background is observed in any of the $\ell^{\pm}\pi^{\mp}$ invariant mass distributions. Limits at 95% confidence level on the sum of the squares of the mixing amplitudes between heavy and light neutrinos, $\vert V_\mathrm{N}\vert^2$, and on $c\tau$ are obtained in different mixing scenarios for both Majorana and Dirac-like N particles. The most stringent upper limit $\vert V_\mathrm{N}\vert^2$ $\lt$ 2.0$\times$10$^{-5}$ is obtained at $m_\mathrm{N}$ = 1.95 GeV for the Majorana case where N mixes exclusively with muon neutrinos. The limits on $\vert V_\mathrm{N}\vert^2$ for masses 1 $\lt$ $m_\mathrm{N}$ $\lt$ 1.7 GeV are the most stringent from a collider experiment to date.

A search for the production of a top quark in association with a photon and additional jets via flavor changing neutral current interactions is presented. The analysis uses proton-proton collision data recorded by the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The search is performed by looking for processes where a single top quark is produced in association with a photon, or a pair of top quarks where one of the top quarks decays into a photon and an up or charm quark. Events with an electron or a muon, a photon, one or more jets, and missing transverse momentum are selected. Multivariate analysis techniques are used to discriminate signal and standard model background processes. No significant deviation is observed over the predicted background. Observed (expected) upper limits are set on the branching fractions of top quark decays: $\mathcal{B}$(t $\to$ u$\gamma$) $\lt$ 0.95 $\times$ 10$^{-5}$ (1.20 $\times$ 10$^{-5}$) and $\mathcal{B}$(t $\to$ c$\gamma$) $\lt$ 1.51 $\times$ 10$^{-5}$ (1.54 $\times$ 10$^{-5}$) at 95% confidence level, assuming a single nonzero coupling at a time. The obtained limit for $\mathcal{B}$(t $\to$ u$\gamma$) is similar to the current best limit, while the limit for $\mathcal{B}$(t $\to$ c$\gamma$) is significantly tighter than previous results.

A search is performed for charged-lepton flavor violating processes in top quark (t) production and decay. The data were collected by the CMS experiment from proton-proton collisions at a center-of-mass energy of 13 TeV and correspond to an integrated luminosity of 138 fb$^{-1}$. The selected events are required to contain one opposite-sign electron-muon pair, a third charged lepton (electron or muon), and at least one jet of which no more than one is associated with a bottom quark. Boosted decision trees are used to distinguish signal from background, exploiting differences in the kinematics of the final states particles. The data are consistent with the standard model expectation. Upper limits at 95% confidence level are placed in the context of effective field theory on the Wilson coefficients, which range between 0.024-0.424 TeV$^{-2}$ depending on the flavor of the associated light quark and the Lorentz structure of the interaction. These limits are converted to upper limits on branching fractions involving up (charm) quarks, t$\to$e$\mu$u (t$\to$e$\mu$c), of 0.032 (0.498)$\times$10$^{-6}$, 0.022 (0.369)$\times$10$^{-6}$, and 0.012 (0.216)$\times$10$^{-6}$ for tensor-like, vector-like, and scalar-like interactions, respectively.

A search for an exotic decay of the Higgs boson to a Z boson and a light pseudoscalar particle (a), decaying to a pair of leptons and a pair of photons, respectively, is presented. The search is based on proton-proton collision data at a center-of-mass energy of $\sqrt{s}$ = 13 TeV, collected with the CMS detector and corresponding to an integrated luminosity of 138 fb$^{-1}$. The analysis probes pseudoscalar masses $m_\mathrm{a}$ between 1 and 30 GeV, leading to two pairs of well-isolated leptons and photons. Upper limits at 95% confidence level are set on the Higgs boson production cross section times its branching fraction to two leptons and two photons. The observed (expected) limits are in the range of 1.1-17.8 (1.7-17.9) fb within the probed $m_\mathrm{a}$ interval. An excess of data above the expected standard model background with a local (global) significance of 2.6 (1.3) standard deviations is observed for a mass hypothesis of $m_\mathrm{a}$ = 3 GeV. Limits on models involving axion-like particles, formulated as an effective field theory, are also reported.

A search is described for the production of a pair of bottom-type vector-like quarks (B VLQs) with mass greater than 1000 GeV. Each B VLQ decays into a b quark and a Higgs boson, a b quark and a Z boson, or a t quark and a W boson. This analysis considers both fully hadronic final states and those containing a charged lepton pair from a Z boson decay. The products of the H $to$ bb boson decay and of the hadronic Z or W boson decays can be resolved as two distinct jets or merged into a single jet, so the final states are classified by the number of reconstructed jets. The analysis uses data corresponding to an integrated luminosity of 138 fb$^{-1}$ collected in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the CMS detector at the LHC from 2016 to 2018. No excess over the expected background is observed. Lower limits are set on the B VLQ mass at 95% confidence level. These depend on the B VLQ branching fractions and are 1570 and 1540 GeV for 100% B $\to$ bH and 100% B $\to$ bZ, respectively. In most cases, the mass limits obtained exceed previous limits by at least 100 GeV.

Measurements of inclusive and normalized differential cross sections of the associated production of top quark-antiquark and bottom quark-antiquark pairs, ttbb, are presented. The results are based on data from proton-proton collisions collected by the CMS detector at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The cross sections are measured in the lepton+jets decay channel of the top quark pair, using events containing exactly one isolated electron or muon and at least five jets. Measurements are made in four fiducial phase space regions, targeting different aspects of the ttbb process. Distributions are unfolded to the particle level through maximum likelihood fits, and compared with predictions from several event generators. The inclusive cross section measurements of this process in the fiducial phase space regions are the most precise to date. In most cases, the measured inclusive cross sections exceed the predictions with the chosen generator settings. The only exception is when using a particular choice of dynamic renormalization scale, $\mu_\mathrm{R}=\frac{1}{2} \prod_{i=\mathrm{t, \bar{t}, b, \bar{b}}} m_{\mathrm{T},i}^{1/4}$, where $m_{\mathrm{T},i}^2=m_i^2+p^2_{\mathrm{T},i}$ are the transverse masses of top and bottom quarks. The differential cross sections show varying degrees of compatibility with the theoretical predictions, and none of the tested generators with the chosen settings simultaneously describe all the measured distributions.

A search is presented for the decay of the 125 GeV Higgs boson (H) to a pair of new light pseudoscalar bosons (a), followed by the prompt decay of each a boson to a bottom quark-antiquark pair, H $\to$ aa $\to$$\mathrm{b\bar{b}b\bar{b}}$. The analysis is performed using a data sample of proton-proton collisions collected with the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. To reduce the background from standard model processes, the search requires the Higgs boson to be produced in association with a leptonically decaying W or Z boson. The analysis probes the production of new light bosons in a 15 $\lt$$m_\mathrm{a}$$\lt$ 60 GeV mass range. Assuming the standard model predictions for the Higgs boson production cross sections for pp $\to$ WH and ZH, model independent upper limits at 95% confidence level are derived for the branching fraction $\mathcal{B}$(H $\to$ aa $\to$ $\mathrm{b\bar{b}b\bar{b}}$). The combined WH and ZH observed upper limit on the branching fraction ranges from 1.10 for $m_\mathrm{a} =$ 20 GeV to 0.36 for $m_\mathrm{a} =$ 60 GeV, complementing other measurements in the $\mu\mu\tau\tau$, $\tau\tau\tau\tau$ and bb$\ell\ell$ ($\ell=$ $\mu$,$\tau$) channels.

The observation of WW$\gamma$ production in proton-proton collisions at a center-of-mass energy of 13 TeV with an integrated luminosity of 138 fb$^{-1}$ is presented. The observed (expected) significance is 5.6 (5.1) standard deviations. Events are selected by requiring exactly two leptons (one electron and one muon) of opposite charge, moderate missing transverse momentum, and a photon. The measured fiducial cross section for WW$\gamma$ is 5.9 $\pm$ 0.8 (stat) $\pm$ 0.8 (syst) $\pm$ 0.7 (modeling) fb, in agreement with the next-to-leading order quantum chromodynamics prediction. The analysis is extended with a search for the associated production of the Higgs boson and a photon, which is generated by a coupling of the Higgs boson to light quarks. The result is used to constrain the Higgs boson couplings to light quarks.

A search for the lepton flavor violating $\tau$$\to$ 3$\mu$ decay is performed using proton-proton collision events at a center-of-mass energy of 13 TeV collected by the CMS experiment at the LHC in 2017-2018, corresponding to an integrated luminosity of 97.7 fb$^{-1}$. Tau leptons produced in both heavy-flavor hadron and W boson decays are exploited in the analysis. No evidence for the decay is observed. The results of this search are combined with an earlier null result based on data collected in 2016 to obtain a total integrated luminosity of 131 fb$^{-1}$. The observed (expected) upper limits on the branching fraction $\mathcal{B}$($\tau$$\to$ 3$\mu$) at confidence levels of 90 and 95% are 2.9 $\times$ 10$^{-8}$ (2.4 $\times$ 10$^{-8}$) and 3.6 $\times$ 10$^{-8}$ (3.0 $\times$ 10$^{-8}$), respectively.

The observation of the production of four top quarks in proton-proton collisions is reported, based on a data sample collected by the CMS experiment at a center-of-mass energy of 13 TeV in 2016-2018 at the CERN LHC and corresponding to an integrated luminosity of 138 fb$^{-1}$. Events with two same-sign, three, or four charged leptons (electrons and muons) and additional jets are analyzed. Compared to previous results in these channels, updated identification techniques for charged leptons and jets originating from the hadronization of b quarks, as well as a revised multivariate analysis strategy to distinguish the signal process from the main backgrounds, lead to an improved expected signal significance of 4.9 standard deviations above the background-only hypothesis. Four top quark production is observed with a significance of 5.6 standard deviations, and its cross section is measured to be 17.7 $^{+3.7}_{-3.5}$ (stat) $^{+2.3}_{-1.9}$ (syst) fb, in agreement with the available standard model predictions.

Measurements of the inclusive and differential fiducial cross sections for the Higgs boson production in the H → ZZ → 4ℓ (ℓ = e, μ) decay channel are presented. The results are obtained from the analysis of proton-proton collision data recorded by the CMS experiment at the CERN LHC at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{−1}$. The measured inclusive fiducial cross section is 2.73 ± 0.26 fb, in agreement with the standard model expectation of 2.86 ± 0.1 fb. Differential cross sections are measured as a function of several kinematic observables sensitive to the Higgs boson production and decay to four leptons. A set of double-differential measurements is also performed, yielding a comprehensive characterization of the four leptons final state. Constraints on the Higgs boson trilinear coupling and on the bottom and charm quark coupling modifiers are derived from its transverse momentum distribution. All results are consistent with theoretical predictions from the standard model.

A study of the anomalous couplings of the Higgs boson to vector bosons, including $CP$-violation effects, has been conducted using its production and decay in the WW channel. This analysis is performed on proton-proton collision data collected with the CMS detector at the CERN LHC during 2016-2018 at a center-of-mass energy of 13 TeV, and corresponds to an integrated luminosity of 138 fb$^{-1}$. The different-flavor dilepton (e$\mu$) final state is analyzed, with dedicated categories targeting gluon fusion, electroweak vector boson fusion, and associated production with a W or Z boson. Kinematic information from associated jets is combined using matrix element techniques to increase the sensitivity to anomalous effects at the production vertex. A simultaneous measurement of four Higgs boson couplings to electroweak vector bosons is performed in the framework of a standard model effective field theory. All measurements are consistent with the expectations for the standard model Higgs boson and constraints are set on the fractional contribution of the anomalous couplings to the Higgs boson production cross section.

A search for dark matter particles is performed using events with a pair of W bosons and large missing transverse momentum. Candidate events are selected by requiring one or two leptons ($\ell =$ electrons or muons). The analysis is based on proton-proton collision data collected at a center-of-mass energy of 13 TeV by the CMS experiment at the LHC and corresponding to an integrated luminosity of 138 fb$^{-1}$. No significant excess over the expected standard model background is observed in the $\ell\nu$qq and 2$\ell$2$\nu$ final states of the W$^+$W$^-$ boson pair. Limits are set on dark matter production in the context of a simplified dark Higgs model, with a dark Higgs boson mass above the W$^+$W$^-$ mass threshold. The dark matter phase space is probed in the mass range 100-300 GeV, extending the scope of previous searches. Current exclusion limits are improved in the range of dark Higgs masses from 160 to 250 GeV, for a dark matter mass of 200 GeV.

A search is presented for baryon number violating interactions in top quark production and decay. The analysis uses data from proton-proton collisions at a center-of-mass energy of 13 TeV, collected with the CMS detector at the LHC with an integrated luminosity of 138 fb$^{-1}$. Candidate events are selected by requiring two oppositely-charged leptons (electrons or muons) and exactly one jet identified as originating from a bottom quark. Multivariate discriminants are used to separate the signal from the background. No significant deviation from the standard model prediction is observed. Upper limits are placed on the strength of baryon number violating couplings. For the first time the production of single top quarks via baryon number violating interactions is studied. This allows the search to set the most stringent constraints to date on the branching fraction of the top quark decay to a lepton, an up-type quark (u or c), and a down-type quark (d, s, or b). The results improve the previous bounds by three to six orders of magnitude based on the fermion flavor combination of the baryon number violating interactions.

A search is presented for a third-generation leptoquark (LQ) coupled exclusively to a $\tau$ lepton and a b quark. The search is based on proton-proton collision data at a center-of-mass energy of 13 TeV recorded with the CMS detector, corresponding to an integrated luminosity of 138 fb$^{-1}$. Events with $\tau$ leptons and a varying number of jets originating from b quarks are considered, targeting the single and pair production of LQs, as well as nonresonant $t$-channel LQ exchange. An excess is observed in the data with respect to the background expectation in the combined analysis of all search regions. For a benchmark LQ mass of 2 TeV and an LQ-b-$\tau$ coupling strength of 2.5, the excess reaches a local significance of up to 2.8 standard deviations. Upper limits at the 95% confidence level are placed on the LQ production cross section in the LQ mass range 0.5-2.3 TeV, and up to 3 TeV for $t$-channel LQ exchange. Leptoquarks are excluded below masses of 1.22-1.88 TeV for different LQ models and varying coupling strengths up to 2.5. The study of nonresonant $\tau\tau$ production through $t$-channel LQ exchange allows lower limits on the LQ mass of up to 2.3 TeV to be obtained.

A test of lepton flavor universality in B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$ and B$^{\pm}$$\to$ K$^{\pm}$e$^+$e$^-$ decays, as well as a measurement of differential and integrated branching fractions of a nonresonant B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$ decay are presented. The analysis is made possible by a dedicated data set of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded in 2018, by the CMS experiment at the LHC, using a special high-rate data stream designed for collecting about 10 billion unbiased b hadron decays. The ratio of the branching fractions $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$) to $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}$e$^+$e$^-$) is determined from the measured double ratio $R$(K) of these decays to the respective branching fractions of the B$^\pm$$\to$ J/$\psi$K$^\pm$ with J/$\psi$$\to$$\mu^+\mu^-$ and e$^+$e$^-$ decays, which allow for significant cancellation of systematic uncertainties. The ratio $R$(K) is measured in the range 1.1 $\lt q^2 \lt$ 6.0 GeV$^2$, where $q$ is the invariant mass of the lepton pair, and is found to be $R$(K) = 0.78$^{+0.47}_{-0.23}$, in agreement with the standard model expectation $R$(K) $\approx$ 1. This measurement is limited by the statistical precision of the electron channel. The integrated branching fraction in the same $q^2$ range, $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$) = (12.42 $\pm$ 0.68) $\times$ 10$^{-8}$, is consistent with the present world-average value and has a comparable precision.

A search for the lepton-flavor violating decay of the Higgs boson and potential additional Higgs bosons with a mass in the range 110-160 GeV to an e$^{\pm}\mu^{\mp}$ pair is presented. The search is performed with a proton-proton collision dataset at a center-of-mass energy of 13 TeV collected by the CMS experiment at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. No excess is observed for the Higgs boson. The observed (expected) upper limit on the e$^{\pm}\mu^{\mp}$ branching fraction for it is determined to be 4.4 (4.7) $\times$ 10$^{-5}$ at 95% confidence level, the most stringent limit set thus far from direct searches. The largest excess of events over the expected background in the full mass range of the search is observed at an e$^{\pm}\mu^{\mp}$ invariant mass of approximately 146 GeV with a local (global) significance of 3.8 (2.8) standard deviations.

A measurement of the Higgs boson (H) production via vector boson fusion (VBF) and its decay into a bottom quark-antiquark pair ($\mathrm{b\bar{b}}$) is presented using proton-proton collision data recorded by the CMS experiment at $\sqrt{s}$ = 13 TeV and corresponding to an integrated luminosity of 90.8 fb$^{-1}$. Treating the gluon-gluon fusion process as a background and constraining its rate to the value expected in the standard model (SM) within uncertainties, the signal strength of the VBF process, defined as the ratio of the observed signal rate to that predicted by the SM, is measured to be $\mu^\text{qqH}_\mathrm{Hb\bar{b}}$ = 1.01 $^{+0.55}_{-0.46}$. The VBF signal is observed with a significance of 2.4 standard deviations relative to the background prediction, while the expected significance is 2.7 standard deviations. Considering inclusive Higgs boson production and decay into bottom quarks, the signal strength is measured to be $\mu^\text{incl.}_\mathrm{Hb\bar{b}}$ = 0.99 $^{+0.48}_{-0.41}$, corresponding to an observed (expected) significance of 2.6 (2.9) standard deviations.

A search is presented for fractionally charged particles with charge below 1$e$, using their small energy loss in the tracking detector as a key variable to observe a signal. The analyzed data set corresponds to an integrated luminosity of 138 fb$^{-1}$ of proton-proton collisions collected at $\sqrt{s}$ = 13 TeV in 2016-2018 at the CERN LHC. This is the first search at the LHC for new particles with charges between $e/$3 and $e$. Masses up to 640 GeV and charges as low as $e/$3 are excluded at 95% confidence level. These are the most stringent limits to date for the considered Drell-Yan-like production mode.

A search for high-mass dimuon resonance production in association with one or more b quark jets is presented. The study uses proton-proton collision data collected with the CMS detector at the LHC corresponding to an integrated luminosity of 138 fb$^{-1}$ at a center-of-mass energy of 13 TeV. Model-independent limits are derived on the number of signal events with exactly one or more than one b quark jet. Results are also interpreted in a lepton-flavor-universal model with Z$'$ boson couplings to a bb quark pair ($g_\mathrm{b}$), an sb quark pair ($g_\mathrm{b}\delta_\mathrm{bs}$), and any same-flavor charged lepton ($g_\ell$) or neutrino pair ($g_\nu$), with $\left|g_{\nu}\right| = \left|g_\ell\right|$. For a Z$'$ boson with a mass $m_{\mathrm{Z}'}$ = 350 GeV (2 TeV) and $\left|\delta_\mathrm{bs}\right|$$\lt$ 0.25, the majority of the parameter space with 0.0057 $\lt$$\left|g_\ell\right|$$\lt$ 0.35 (0.25 $\lt$$\left|g_\ell\right|$$\lt$ 0.43) and 0.0079 $\lt$$\left|g_\mathrm{b}\right|$$\lt$ 0.46 (0.34 $\lt$$\left|g_\mathrm{b}\right|$$\lt$ 0.57) is excluded at 95% confidence level. Finally, constraints are set on a Z$'$ model with parameters consistent with low-energy b $\to$ s$\ell\ell$ measurements. In this scenario, most of the allowed parameter space is excluded for a Z$'$ boson with 350 $\lt m_{\mathrm{Z}'}$ $\lt$ 500 GeV, while the constraints are less stringent for higher $m_{\mathrm{Z}'}$ hypotheses. This is the first dedicated search at the LHC for a high-mass dimuon resonance produced in association with multiple b quark jets, and the constraints obtained on models with this signature are the most stringent to date.

The first observation of electroweak production of same-sign W boson pairs in proton-proton collisions is reported. The data sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$ collected at a center-of-mass energy of 13 TeV with the CMS detector at the LHC. Events are selected by requiring exactly two leptons (electrons or muons) of the same charge, moderate missing transverse momentum, and two jets with a large rapidity separation and a large dijet mass. The observed significance of the signal is 5.5 standard deviations, where a significance of 5.7 standard deviations is expected based on the standard model. The ratio of measured event yields to that expected from the standard model at leading-order is 0.90 $\pm$ 0.22. A cross section measurement in a fiducial region is reported. Bounds are given on the structure of quartic vector boson interactions in the framework of dimension-eight effective field theory operators and on the production of doubly charged Higgs bosons.

The angular distribution of the flavor-changing neutral current decay B$^+$$\to$ K$^+\mu^+\mu^-$ is studied in proton-proton collisions at a center-of-mass energy of 8 TeV. The analysis is based on data collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 20.5 fb$^{-1}$. The forward-backward asymmetry $A_{\mathrm{FB}}$ of the dimuon system and the contribution $F_{\mathrm{H}}$ from the pseudoscalar, scalar, and tensor amplitudes to the decay width are measured as a function of the dimuon mass squared. The measurements are consistent with the standard model expectations.

The effect of event background fluctuations on charged particle jet reconstruction in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV has been measured with the ALICE experiment. The main sources of non-statistical fluctuations are characterized based purely on experimental data with an unbiased method, as well as by using single high $p_{\rm T}$ particles and simulated jets embedded into real Pb-Pb events and reconstructed with the anti-$k_{\rm T}$ jet finder. The influence of a low transverse momentum cut-off on particles used in the jet reconstruction is quantified by varying the minimum track $p_{\rm T}$ between 0.15 GeV/$c$ and 2 GeV/$c$. For embedded jets reconstructed from charged particles with $p_{\rm T} > 0.15$ GeV/$c$, the uncertainty in the reconstructed jet transverse momentum due to the heavy-ion background is measured to be 11.3 GeV/$c$ (standard deviation) for the 10% most central Pb-Pb collisions, slightly larger than the value of 11.0 GeV/$c$ measured using the unbiased method. For a higher particle transverse momentum threshold of 2 GeV/$c$, which will generate a stronger bias towards hard fragmentation in the jet finding process, the standard deviation of the fluctuations in the reconstructed jet transverse momentum is reduced to 4.8-5.0 GeV/$c$ for the 10% most central events. A non-Gaussian tail of the momentum uncertainty is observed and its impact on the reconstructed jet spectrum is evaluated for varying particle momentum thresholds, by folding the measured fluctuations with steeply falling spectra.

Measurements of inclusive charged-hadron transverse-momentum and pseudorapidity distributions are presented for proton-proton collisions at sqrt(s) = 0.9 and 2.36 TeV. The data were collected with the CMS detector during the LHC commissioning in December 2009. For non-single-diffractive interactions, the average charged-hadron transverse momentum is measured to be 0.46 +/- 0.01 (stat.) +/- 0.01 (syst.) GeV/c at 0.9 TeV and 0.50 +/- 0.01 (stat.) +/- 0.01 (syst.) GeV/c at 2.36 TeV, for pseudorapidities between -2.4 and +2.4. At these energies, the measured pseudorapidity densities in the central region, dN(charged)/d(eta) for |eta| < 0.5, are 3.48 +/- 0.02 (stat.) +/- 0.13 (syst.) and 4.47 +/- 0.04 (stat.) +/- 0.16 (syst.), respectively. The results at 0.9 TeV are in agreement with previous measurements and confirm the expectation of near equal hadron production in p-pbar and pp collisions. The results at 2.36 TeV represent the highest-energy measurements at a particle collider to date.

A measurement of the transverse momentum spectra of jets in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV is reported. Jets are reconstructed from charged particles using the anti-$k_{\rm T}$ jet algorithm with jet resolution parameters $R$ of $0.2$ and $0.3$ in pseudo-rapidity $|\eta|<0.5$. The transverse momentum $p_{\rm T}$ of charged particles is measured down to $0.15$ GeV/$c$ which gives access to the low $p_{\rm T}$ fragments of the jet. Jets found in heavy-ion collisions are corrected event-by-event for average background density and on an inclusive basis (via unfolding) for residual background fluctuations and detector effects. A strong suppression of jet production in central events with respect to peripheral events is observed. The suppression is found to be similar to the suppression of charged hadrons, which suggests that substantial energy is radiated at angles larger than the jet resolution parameter $R=0.3$ considered in the analysis. The fragmentation bias introduced by selecting jets with a high $p_{\rm T}$ leading particle, which rejects jets with a soft fragmentation pattern, has a similar effect on the jet yield for central and peripheral events. The ratio of jet spectra with $R=0.2$ and $R=0.3$ is found to be similar in Pb-Pb and simulated PYTHIA pp events, indicating no strong broadening of the radial jet structure in the reconstructed jets with $R<0.3$.

Measurements of the sphericity of primary charged particles in minimum bias proton--proton collisions at $\sqrt{s}=0.9$, 2.76 and 7 TeV with the ALICE detector at the LHC are presented. The observable is linearized to be collinear safe and is measured in the plane perpendicular to the beam direction using primary charged tracks with $p_{\rm T}\geq0.5$ GeV/c in $|\eta|\leq0.8$. The mean sphericity as a function of the charged particle multiplicity at mid-rapidity ($N_{\rm ch}$) is reported for events with different $p_{\rm T}$ scales ("soft" and "hard") defined by the transverse momentum of the leading particle. In addition, the mean charged particle transverse momentum versus multiplicity is presented for the different event classes, and the sphericity distributions in bins of multiplicity are presented. The data are compared with calculations of standard Monte Carlo event generators. The transverse sphericity is found to grow with multiplicity at all collision energies, with a steeper rise at low $N_{\rm ch}$, whereas the event generators show the opposite tendency. The combined study of the sphericity and the mean $p_{\rm T}$ with multiplicity indicates that most of the tested event generators produce events with higher multiplicity by generating more back-to-back jets resulting in decreased sphericity (and isotropy). The PYTHIA6 generator with tune PERUGIA-2011 exhibits a noticeable improvement in describing the data, compared to the other tested generators.

A measurement of the cross section of the associated production of a single top quark and a W boson in final states with a muon or electron and jets in proton-proton collisions at $\sqrt{s}$ = 13 TeV is presented. The data correspond to an integrated luminosity of 36 fb$^{-1}$ collected with the CMS detector at the CERN LHC in 2016. A boosted decision tree is used to separate the tW signal from the dominant $\mathrm{t\bar{t}}$ background, whilst the subleading W+jets and multijet backgrounds are constrained using data-based estimates. This result is the first observation of the tW process in final states containing a muon or electron and jets, with a significance exceeding 5 standard deviations. The cross section is determined to be 89 $\pm$ 4 (stat) $\pm$ 12 (syst) pb, consistent with the standard model.

This publication describes the methods used to measure the centrality of inelastic Pb-Pb collisions at a center-of-mass energy of 2.76 TeV per colliding nucleon pair with ALICE. The centrality is a key parameter in the study of the properties of QCD matter at extreme temperature and energy density, because it is directly related to the initial overlap region of the colliding nuclei. Geometrical properties of the collision, such as the number of participating nucleons and number of binary nucleon-nucleon collisions, are deduced from a Glauber model with a sharp impact parameter selection, and shown to be consistent with those extracted from the data. The centrality determination provides a tool to compare ALICE measurements with those of other experiments and with theoretical calculations.

In this paper measurements are presented of $\rm \pi$$^+$, $\rm \pi$$^-$, K$^+$, K$^-$, p and $\overline{\rm p}$ production at mid-rapidity < 0.5, in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV as a function of centrality. The measurement covers the transverse momentum ($p_{\rm T}$) range from 100, 200, 300 MeV/$c$ up to 3, 3, 4.6 GeV/$c$, for $\rm\pi$, K, and p respectively. The measured $p_{\rm T}$ distributions and yields are compared to expectations based on hydrodynamic, thermal and recombination models. The spectral shapes of central collisions show a stronger radial flow than measured at lower energies, which can be described in hydrodynamic models. In peripheral collisions, the $p_{\rm T}$ distributions are not well reproduced by hydrodynamic models. Ratios of integrated particle yields are found to be nearly independent of centrality. The yield of protons normalized to pions is a factor ~1.5 lower than the expectation from thermal models.

The production cross section of electrons from semileptonic decays of beauty hadrons was measured at mid-rapidity (|y| < 0.8) in the transverse momentum range $1 < p_{\rm T} < 8$ Gev/$c$ with the ALICE experiment at the CERN LHC in pp collisions at a center of mass energy $\sqrt{s} = 7$ TeV using an integrated luminosity of 2.2 nb$^{-1}$. Electrons from beauty hadron decays were selected based on the displacement of the decay vertex from the collision vertex. A perturbative QCD calculation agrees with the measurement within uncertainties. The data were extrapolated to the full phase space to determine the total cross section for the production of beauty quark-antiquark pairs.

A search for the standard model Higgs boson produced in association with a W or a Z boson and decaying a pair of $\tau$ leptons is performed. A data sample of proton-proton collisions collected at $\sqrt{s} =$ 13 TeV by the CMS experiment at the CERN LHC is used, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The signal strength is measured relative to the expectation for the standard model Higgs boson, yielding $\mu =$ 2.5 $^{+1.4} _{-1.3}$. These results are combined with earlier CMS measurements targeting Higgs boson decays to a pair of $\tau$ leptons, performed with the same data set in the gluon fusion and vector boson fusion production modes. The combined signal strength is $\mu =$ 1.24 $^{+0.29} _{-0.27}$ (1.00 $^{+0.24} _{-0.23}$ expected), and the observed significance is 5.5 standard deviations (4.8 expected) for a Higgs boson mass of 125 GeV.

A measurement is presented of the charged hadron multiplicity in hadronic PbPb collisions, as a function of pseudorapidity and centrality, at a collision energy of 2.76 TeV per nucleon pair. The data sample is collected using the CMS detector and a minimum-bias trigger, with the CMS solenoid off. The number of charged hadrons is measured both by counting the number of reconstructed particle hits and by forming hit doublets of pairs of layers in the pixel detector. The two methods give consistent results. The charged hadron multiplicity density dN(ch)/d eta, evaluated at eta=0 for head-on collisions, is found to be 1612 +/- 55, where the uncertainty is dominated by systematic effects. Comparisons of these results to previous measurements and to various models are also presented.

A search for WW production from double-parton scattering processes using same-charge electron-muon and dimuon events is reported, based on proton-proton collision data collected at a center-of-mass energy of 13 TeV. The analyzed data set corresponds to an integrated luminosity of 77.4 fb$^{-1}$, collected using the CMS detector at the LHC in 2016 and 2017. Multivariate classifiers are used to discriminate between the signal and the dominant background processes. A maximum likelihood fit is performed to extract the signal cross section. This leads to the first evidence for WW production via double-parton scattering, with a significance of 3.9 standard deviations. The measured inclusive cross section is 1.41 $\pm$ 0.28 (stat) $\pm$ 0.28 (syst) pb.