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.

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.

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 first search for narrow resonances decaying to three well-separated hadronic jets is presented. The search uses proton-proton collision data corresponding to an integrated luminosity of 138 fb$^{-1}$ at $\sqrt{s}$ = 13 TeV, collected at the CERN LHC. No significant deviations from the background predictions are observed between 1.75-9.00 TeV. The results provide the first mass limits on a right-handed boson Z$_{\mathrm{R}}$ decaying to three gluons, an excited quark decaying via a vector boson to three quarks, as well as updated limits on a Kaluza-Klein gluon decaying via a radion to three gluons.

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.

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.

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.

Results are presented from a search for CP violation in top quark pair production, using proton-proton collisions at a center-of-mass energy of 13 TeV. The data used for this analysis consist of final states with two charged leptons collected by the CMS experiment, and correspond to an integrated luminosity of 35.9 fb$^{-1}$. The search uses two observables, $\mathcal{O}_1$ and $\mathcal{O}_3$, which are Lorentz scalars. The observable $\mathcal{O}_1$ is constructed from the four-momenta of the charged leptons and the reconstructed top quarks, while $\mathcal{O}_3$ consists of the four-momenta of the charged leptons and the b quarks originating from the top quarks. Asymmetries in these observables are sensitive to CP violation, and their measurement is used to determine the chromoelectric dipole moment of the top quark. The results are consistent with the expectation from the standard model.

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.

The production of Z bosons associated with jets is measured in pp collisions at $\sqrt{s}$ = 13 TeV with data recorded with the CMS experiment at the LHC corresponding to an integrated luminosity of 36.3 fb$^{-1}$. The multiplicity of jets with transverse momentum $p_\mathrm{T}$$\gt$ 30 GeV is measured for different regions of the Z boson's $p_\mathrm{T}$(Z), from lower than 10 GeV to higher than 100 GeV. The azimuthal correlation $\Delta \phi$ between the Z boson and the leading jet, as well as the correlations between the two leading jets are measured in three regions of $p_\mathrm{T}$(Z). The measurements are compared with several predictions at leading and next-to-leading orders, interfaced with parton showers. Predictions based on transverse-momentum dependent parton distributions and corresponding parton showers give a good description of the measurement in the regions where multiple parton interactions and higher jet multiplicities are not important. The effects of multiple parton interactions are shown to be important to correctly describe the measured spectra in the low $p_\mathrm{T}$(Z) regions.

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.

The production of four top quarks ($\mathrm{t\bar{t}t\bar{t}}$) is studied with LHC proton-proton collision data samples collected by the CMS experiment at a center-of-mass energy of 13 TeV, and corresponding to integrated luminosities of up to 138 fb$^{-1}$. Events that have no leptons (all-hadronic), one lepton, or two opposite-sign leptons (where lepton refers only to prompt electrons or prompt muons) are considered. This is the first $\mathrm{t\bar{t}t\bar{t}}$ measurement that includes the all-hadronic final state. The observed significance of the $\mathrm{t\bar{t}t\bar{t}}$ signal in these final states of 3.9 standard deviations (1.5 expected) provides evidence for $\mathrm{t\bar{t}t\bar{t}}$ production, with a measured cross section of 36 $^{+12}_{-11}$ fb. Combined with earlier CMS results in other final states, the signal significance is 4.0 standard deviations (3.2 expected). The combination returns an observed cross section of 17 $\pm$ 4 (stat) $\pm$ 3 (syst) fb, which is consistent with the standard model prediction.

Protons consist of three valence quarks, two up-quarks and one down-quark, held together by gluons and a sea of quark-antiquark pairs. Collectively, quarks and gluons are referred to as partons. In a proton-proton collision, typically only one parton of each proton undergoes a hard scattering - referred to as single-parton scattering - leaving the remainder of each proton only slightly disturbed. Here, we report the study of double- and triple-parton scatterings through the simultaneous production of three J/$\psi$ mesons, which consist of a charm quark-antiquark pair, in proton-proton collisions recorded with the CMS experiment at the Large Hadron Collider. We observed this process - reconstructed through the decays of J/$\psi$ mesons into pairs of oppositely charged muons - with a statistical significance above five standard deviations. We measured the inclusive fiducial cross section to be 272 $^{+141}_{-104}$ (stat) $\pm$ 17 (syst) fb, and compared it to theoretical expectations for triple-J/$\psi$ meson production in single-, double- and triple-parton scattering scenarios. Assuming factorization of multiple hard-scattering probabilities in terms of single-parton scattering cross sections, double- and triple-parton scattering are the dominant contributions for the measured process.

The first collider search for dark matter arising from a strongly coupled hidden sector is presented and uses a data sample corresponding to 138 fb$^{-1}$, collected with the CMS detector at the CERN LHC, at $\sqrt{s} =$ 13 TeV. The hidden sector is hypothesized to couple to the standard model (SM) via a heavy leptophobic Z' mediator produced as a resonance in proton-proton collisions. The mediator decay results in two "semivisible" jets, containing both visible matter and invisible dark matter. The final state therefore includes moderate missing energy aligned with one of the jets, a signature ignored by most dark matter searches. No structure in the dijet transverse mass spectra compatible with the signal is observed. Assuming the Z' has a universal coupling of 0.25 to the SM quarks, an inclusive search, relevant to any model that exhibits this kinematic behavior, excludes mediator masses of 1.5-4.0 TeV at 95% confidence level, depending on the other signal model parameters. To enhance the sensitivity of the search for this particular class of hidden sector models, a boosted decision tree (BDT) is trained using jet substructure variables to distinguish between semivisible jets and SM jets from background processes. When the BDT is employed to identify each jet in the dijet system as semivisible, the mediator mass exclusion increases to 5.1 TeV, for wider ranges of the other signal model parameters. These limits exclude a wide range of strongly coupled hidden sector models for the first time.

A precise measurement of the invisible width of the Z boson produced in proton-proton collisions at a center-of-mass energy of 13 TeV is presented using data recorded by the CMS experiment at the LHC, corresponding to an integrated luminosity of 36.3 fb$^{-1}$. The result is obtained from a simultaneous fit to kinematic distributions for two data samples of Z boson plus jets: one dominated by Z boson decays to invisible particles and the other by Z boson decays to muon and electron pairs. The invisible width is measured to be 523 $\pm$ 3 (stat) $\pm$ 16 (syst) MeV. This result is the first precise measurement of the invisible width of the Z boson at a hadron collider, and is the single most precise direct measurement to date, competitive with the combined result of the direct measurements from the LEP experiments.

Diboson production in association with jets is studied in the fully leptonic final states, pp $\to$ (Z$\gamma^*$)(Z/$\gamma^*$)+jets $\to$ 2$\ell$2$\ell'$+jets, ($\ell,\ell'$ = e or $\mu$) in proton-proton collisions at a center-of-mass energy of 13 TeV. The data sample corresponds to an integrated luminosity of 138 fb$^{-1}$ collected with the CMS detector at the LHC. Differential distributions and normalized differential cross sections are measured as a function of jet multiplicity, transverse momentum $p_\mathrm{T}$, pseudorapidity $\eta$, invariant mass and $\Delta\eta$ of the highest-$p_\mathrm{T}$ and second-highest-$p_\mathrm{T}$ jets, and as a function of invariant mass of the four-lepton system for events with various jet multiplicities. These differential cross sections are compared with theoretical predictions that mostly agree with the experimental data. However, in a few regions we observe discrepancies between the predicted and measured values. Further improvement of the predictions is required to describe the ZZ+jets production in the whole phase space.

Measurements of jet substructure describing the composition of quark- and gluon-initiated jets are presented. Proton-proton (pp) collision data at $\sqrt{s}$ =13 TeV collected with the CMS detector are used, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Generalized angularities are measured that characterize the jet substructure and distinguish quark- and gluon-initiated jets. These observables are sensitive to the distributions of transverse momenta and angular distances within a jet. The analysis is performed using a data sample of dijet events enriched in gluon-initiated jets, and, for the first time, a Z+jet event sample enriched in quark-initiated jets. The observables are measured in bins of jet transverse momentum, and as a function of the jet radius parameter. Each measurement is repeated applying a "soft drop" grooming procedure that removes soft and large angle radiation from the jet. Using these measurements, the ability of various models to describe jet substructure is assessed, showing a clear need for improvements in Monte Carlo generators.

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.

A search is reported for pairs of light Higgs bosons (H$_1$) produced in supersymmetric cascade decays in final states with small missing transverse momentum. A data set of LHC pp collisions collected with the CMS detector at $\sqrt{s}$ = 13 TeV and corresponding to an integrated luminosity of 138 fb$^{-1}$ is used. The search targets events where both H$_1$ bosons decay into $\mathrm{b\bar{b}}$ pairs that are reconstructed as large-radius jets using substructure techniques. No evidence is found for an excess of events beyond the background expectations of the standard model (SM). Results from the search are interpreted in the next-to-minimal supersymmetric extension of the SM, where a "singlino" of small mass leads to squark and gluino cascade decays that can predominantly end in a highly Lorentz-boosted singlet-like H$_1$ and a singlino-like neutralino of small transverse momentum. Upper limits are set on the product of the squark or gluino pair production cross section and the square of the $\mathrm{b\bar{b}}$ branching fraction of the H$_1$ in a benchmark model containing almost mass-degenerate gluinos and light-flavour squarks. Under the assumption of an SM-like H$_1$$\to$$\mathrm{b\bar{b}}$ branching fraction, H$_1$ bosons with masses in the range 40-120 GeV arising from the decays of squarks or gluinos with a mass of 1200 to 2500 GeV are excluded at 95% confidence level.

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.

Measurements of the second Fourier harmonic coefficient ($v_2$) of the azimuthal distributions of prompt and nonprompt D$^0$ mesons produced in pp and pPb collisions are presented. Nonprompt D$^0$ mesons come from beauty hadron decays. The data samples are collected by the CMS experiment at nucleon-nucleon center-of-mass energies of 13 and 8.16 TeV, respectively. In high multiplicity pp collisions, $v_2$ signals for prompt charm hadrons are reported for the first time, and are found to be comparable to those for light-flavor hadron species over a transverse momentum ($p_\mathrm{T}$) range of 2-6 GeV. Compared at similar event multiplicities, the prompt D$^0$ meson $v_2$ values in pp and pPb collisions are similar in magnitude. The $v_2$ values for open beauty hadrons are extracted for the first time via nonprompt D$^0$ mesons in pPb collisions. For $p_\mathrm{T}$ in the range of 2-5 GeV, the results suggest that $v_2$ for nonprompt D$^0$ mesons are smaller than those for prompt D$^0$ mesons. These new measurements indicate a positive charm hadron $v_2$ in pp collisions and suggest a mass dependence in $v_2$ between charm and beauty hadrons in the pPb system. These results provide insights into the origin of heavy-flavor quark collectivity in small systems.

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.

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 resonances decaying into a W boson and a radion, where the radion decays into two W bosons, is presented. The data analyzed correspond to an integrated luminosity of 138 fb$^{-1}$ recorded in proton-proton collisions with the CMS detector at $\sqrt{s} =$ 13 TeV. One isolated charged lepton is required, together with missing transverse momentum and one or two massive large-radius jets, containing the decay products of either two or one W bosons, respectively. No excess over the background estimation is observed. The results are combined with those from a complementary channel with an all-hadronic final state, described in an accompanying paper. Limits are set on parameters of an extended warped extra-dimensional model. These searches are the first of their kind at the LHC.

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.

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.

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).

Measurements are presented of the differential cross sections for Z bosons produced in proton-proton collisions at $\sqrt{s} =$ 13 TeV and decaying to muons and electrons. The data analyzed were collected in 2016 with the CMS detector at the LHC and correspond to an integrated luminosity of 35.9 fb$^{-1}$. The measured fiducial inclusive product of cross section and branching fraction agrees with next-to-next-to-leading order quantum chromodynamics calculations. Differential cross sections of the transverse momentum $p_\mathrm{T}$, the optimized angular variable $\phi^*_\eta$, and the rapidity of lepton pairs are measured. The data are corrected for detector effects and compared to theoretical predictions using fixed order, resummed, and parton shower calculations. The uncertainties of the measured normalized cross sections are smaller than 0.5% for $\phi^*_\eta <$ 0.5 and for $p_\mathrm{T}^\mathrm{Z} <$ 50 GeV.

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 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 measurement of the double-differential inclusive jet cross section as a function of jet transverse momentum pT and absolute jet rapidity |y| is presented. The analysis is based on proton-proton collisions collected by the CMS experiment at the LHC at a centre-of-mass energy of 13 TeV. The data samples correspond to integrated luminosities of 71 and 44 inverse picobarns for |y| < 3 and 3.2 < |y| < 4.7, respectively. Jets are reconstructed with the anti-kt clustering algorithm for two jet sizes, R, of 0.7 and 0.4, in a phase space region covering jet pT up to 2 TeV and jet rapidity up to |y| = 4.7. Predictions of perturbative quantum chromodynamics at next-to-leading order precision, complemented with electroweak and nonperturbative corrections, are used to compute the absolute scale and the shape of the inclusive jet cross section. The cross section difference in R, when going to a smaller jet size of 0.4, is best described by Monte Carlo event generators with next-to-leading order predictions matched to parton showering, hadronisation, and multiparton interactions. In the phase space accessible with the new data, this measurement provides a first indication that jet physics is as well understood at sqrt(s) = 13 TeV as at smaller centre-of-mass energies.

Results on two-particle angular correlations for charged particles produced in pp collisions at a center-of-mass energy of 13 TeV are presented. The data were taken with the CMS detector at the LHC and correspond to an integrated luminosity of about 270 inverse nanobarns. The correlations are studied over a broad range of pseudorapidity (abs(eta) < 2.4) and over the full azimuth (phi) as a function of charged particle multiplicity and transverse momentum (pt). In high-multiplicity events, a long-range (abs(Delta eta) > 2.0), near-side (Delta phi approximately 0) structure emerges in the two-particle Delta eta-Delta phi correlation functions. The magnitude of the correlation exhibits a pronounced maximum in the range 1.0 < pt < 2.0 GeV/c and an approximately linear increase with the charged particle multiplicity, with an overall correlation strength similar to that found in earlier pp data at sqrt(s) = 7 TeV. The present measurement extends the study of near-side long-range correlations up to charged particle multiplicities of N[ch] approximately 180, a region so far unexplored in pp collisions. The observed long-range correlations are compared to those seen in pp, pPb, and PbPb collisions at lower collision energies.

The top quark pair production cross section is measured for the first time in proton-proton collisions at sqrt(s) = 13 TeV by the CMS experiment at the CERN LHC, using data corresponding to an integrated luminosity of 43 inverse picobarns. The measurement is performed by analyzing events with at least one electron and one muon of opposite charge, and at least two jets. The measured cross section is 746 +/- 58 (stat) +/- 53 (syst) +/- 36 (lumi) pb, in agreement with the expectation from the standard model.

A search for top squark pair production in pp collisions at sqrt(s) = 13 TeV is performed using events with a single isolated electron or muon, jets, and a large transverse momentum imbalance. The results are based on data collected in 2016 with the CMS detector at the LHC, corresponding to an integrated luminosity of 35.9 inverse femtobarns. No significant excess of events is observed above the expectation from standard model processes. Exclusion limits are set in the context of supersymmetric models of pair production of top squarks that decay either to a top quark and a neutralino or to a bottom quark and a chargino. Depending on the details of the model, we exclude top squarks with masses as high as 1120 GeV. Detailed information is also provided to facilitate theoretical interpretations in other scenarios of physics beyond the standard model.

The leptonic and inclusive hadronic decay branching fractions of the W boson are measured using proton-proton collision data collected at $\sqrt{s} =$ 13 TeV by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Events characterized by the production of one or two W bosons are selected and categorized based on the multiplicity and flavor of reconstructed leptons, the number of jets, and the number of jets identified as originating from the hadronization of b quarks. A binned maximum likelihood estimate of the W boson branching fractions is performed simultaneously in each event category. The measured branching fractions of the W boson decaying into electron, muon, and tau lepton final states are (10.83 $\pm$ 0.10)%, (10.94 $\pm$ 0.08)%, and (10.77 $\pm$ 0.21)%, respectively, consistent with lepton flavor universality for the weak interaction. The average leptonic and inclusive hadronic decay branching fractions are estimated to be (10.89 $\pm$ 0.08)% and (67.32 $\pm$ 0.23)%, respectively. Based on the hadronic branching fraction, three standard model quantities are subsequently derived: the sum of squared elements in the first two rows of the Cabibbo-Kobayashi-Maskawa (CKM) matrix $\sum_{ij}\vert V_{ij}\vert^{2}$ = 1.984 $\pm$ 0.021, the CKM element $\vert V_\mathrm{cs}\vert$ = 0.967 $\pm$ 0.011, and the strong coupling constant at the W boson mass scale, $\alpha_\mathrm{S}(m^2_\mathrm{W})$ = 0.095 $\pm$ 0.033.

Central exclusive and semiexclusive production of $\pi^+\pi^-$ pairs is measured with the CMS detector in proton-proton collisions at the LHC at center-of-mass energies of 5.02 and 13 TeV. The theoretical description of these nonperturbative processes, which have not yet been measured in detail at the LHC, poses a significant challenge to models. The two pions are measured and identified in the CMS silicon tracker based on specific energy loss, whereas the absence of other particles is ensured by calorimeter information. The total and differential cross sections of exclusive and semiexclusive central $\pi^+\pi^-$ production are measured as functions of invariant mass, transverse momentum, and rapidity of the $\pi^+\pi^-$ system in the fiducial region defined as transverse momentum $p_\mathrm{T}(\pi)$ $>$ 0.2 GeV and pseudorapidity $|\eta(\pi)|$ $<$ 2.4. The production cross sections for the four resonant channels f$_0(500)$, $\rho^0(770)$, f$_0(980)$, and f$_2(1270)$ are extracted using a simple model. These results represent the first measurement of this process at the LHC collision energies of 5.02 and 13 TeV.

A search for the Higgs boson decaying to two oppositely charged muons is presented using data recorded by the CMS experiment at the CERN LHC in 2016 at a center-of-mass energy $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Data are found to be compatible with the predicted background. For a Higgs boson with a mass of 125.09 GeV, the 95% confidence level observed (background-only expected) upper limit on the production cross section times branching fraction to a pair of muons is found to be 3.0 (2.5) times the standard model expectation. In combination with data recorded at center-of-mass energies $\sqrt{s} =$ 7 and 8 TeV, the background-only expected upper limit improves to 2.2 times the standard model value with a standard model expected significance of 1.0 standard deviations. The corresponding observed upper limit is 2.9 with an observed significance of 0.9 standard deviations. This corresponds to an observed upper limit on the standard model Higgs boson branching fraction to muons of 6.4 $\times$ 10$^{-4}$ and to an observed signal strength of 1.0 $\pm$ 1.0 (stat) $\pm$ 0.1 (syst).

Measurements of the differential cross section for the Drell-Yan process, based on proton-proton collision data at a centre-of-mass energy of 13 TeV, collected by the CMS experiment, are presented. The data correspond to an integrated luminosity of 2.8 (2.3) fb$^{-1}$ in the dimuon (dielectron) channel. The total and fiducial cross section measurements are presented as a function of dilepton invariant mass in the range 15 to 3000 GeV, and compared with the perturbative predictions of the standard model. The measured differential cross sections are in good agreement with the theoretical calculations.

A search for long-lived particles decaying into jets is presented. Data were collected with the CMS detector at the LHC from proton-proton collisions at a center-of-mass energy of 13 TeV in 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The search examines the distinctive topology of displaced tracks and secondary vertices. The selected events are found to be consistent with standard model predictions. For a simplified model in which long-lived neutral particles are pair produced and decay to two jets, pair production cross sections larger than 0.2 fb are excluded at 95% confidence level for a long-lived particle mass larger than 1000 GeV and proper decay lengths between 3 and 130 mm. Several supersymmetry models with gauge-mediated supersymmetry breaking or $R$-parity violation, where pair-produced long-lived gluinos or top squarks decay to several final-state topologies containing displaced jets, are also tested. For these models, in the mass ranges above 200 GeV, gluino masses up to 2300-2400 GeV and top squark masses up to 1350-1600 GeV are excluded for proper decay lengths approximately between 10 and 100 mm. These are the most restrictive limits to date on these models.

A search for dark matter in the form of strongly interacting massive particles (SIMPs) using the CMS detector at the LHC is presented. The SIMPs would be produced in pairs that manifest themselves as pairs of jets without tracks. The energy fraction of jets carried by charged particles is used as a key discriminator to suppress efficiently the large multijet background, and the remaining background is estimated directly from data. The search is performed using proton-proton collision data corresponding to an integrated luminosity of 16.1 fb$^{-1}$, collected with the CMS detector in 2016. No significant excess of events is observed above the expected background. For the simplified dark matter model under consideration, SIMPs with masses up to 100 GeV are excluded and further sensitivity is explored towards higher masses.

A search is performed for neutral non-standard-model Higgs bosons decaying to two muons in the context of the minimal supersymmetric standard model (MSSM). Proton-proton collision data recorded by the CMS experiment at the CERN Large Hadron Collider at a center-of-mass energy of 13 TeV were used, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The search is sensitive to neutral Higgs bosons produced via the gluon fusion process or in association with a $\mathrm{b\overline{b}}$ quark pair. No significant deviations from the standard model expectation are observed. Upper limits at 95% confidence level are set in the context of the $m_\mathrm{h}^{\text{mod+}}$ and phenomenological MSSM scenarios on the parameter $\tan\beta$ as a function of the mass of the pseudoscalar A boson, in the range from 130 to 600 GeV. The results are also used to set a model-independent limit on the product of the branching fraction for the decay into a muon pair and the cross section for the production of a scalar neutral boson, either via gluon fusion, or in association with b quarks, in the mass range from 130 to 1000 GeV.

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}$.

Searches for invisible decays of the Higgs boson are presented. The data collected with the CMS detector at the LHC correspond to integrated luminosities of 5.1, 19.7, and 2.3 inverse femtobarns at centre-of-mass energies of 7, 8, and 13 TeV, respectively. The search channels target Higgs boson production via gluon fusion, vector boson fusion, and in association with a vector boson. Upper limits are placed on the branching fraction of the Higgs boson decay to invisible particles, as a function of the assumed production cross sections. The combination of all channels, assuming standard model production, yields an observed (expected) upper limit on the invisible branching fraction of 0.24 (0.23) at the 95% confidence level. The results are also interpreted in the context of Higgs-portal dark matter models.

A search for heavy neutral leptons (HNLs), the right-handed Dirac or Majorana neutrinos, is performed in final states with three charged leptons (electrons or muons) using proton-proton collision data collected by the CMS experiment at $\sqrt{s} =$ 13 TeV at the CERN LHC. The data correspond to an integrated luminosity of 138 fb$^{-1}$. The HNLs could be produced through mixing with standard model neutrinos $\nu$. For small values of the HNL mass ($\lt$ 20 GeV) and the square of the HNL-$\nu$ mixing parameter (10$^{-7}$-10$^{-2}$), the decay length of these particles can be large enough so that the secondary vertex of the HNL decay can be resolved with the CMS silicon tracker. The selected final state consists of one lepton emerging from the primary proton-proton collision vertex, and two leptons forming a displaced, secondary vertex. No significant deviations from the standard model expectations are observed, and constraints are obtained on the HNL mass and coupling strength parameters, excluding previously unexplored regions of parameter space in the mass range 1-20 GeV and squared mixing parameter values as low as 10$^{-7}$.

The measurement of the cross section for the production of a Z boson, decaying to dielectrons or dimuons, in association with at least one bottom quark jet are performed with proton-proton collision data at $\sqrt{s} =$ 13 TeV. The data sample corresponds to an integrated luminosity of 137 fb$^{-1}$, collected by the CMS experiment at the LHC during 2016-2018. The integrated cross sections for Z + $\ge$ 1 b jet and Z + $\ge$ 2 b jets are reported for the electron, muon, and combined channels. The fiducial cross sections in the combined channel are 6.52 $\pm$ 0.04 (stat) $\pm$ 0.40 (syst) $\pm$ 0.14 (theo) pb for Z + $\ge$ 1 b jet and 0.65 $\pm$ 0.03 (stat) $\pm$ 0.07 (syst) $\pm$ 0.02 (theo) pb for Z + $\ge$ 2 b jets. The differential cross section distributions are measured as functions of various kinematic observables that are useful for precision tests of perturbative quantum chromodynamics predictions. The ratios of integrated and differential cross sections for Z + $\ge$ 2 b jets and Z + $\ge$ 1 b jet processes are also determined. The value of the integrated cross section ratio measured in the combined channel is 0.100 $\pm$ 0.005 (stat) $\pm$ 0.007 (syst) $\pm$ 0.003 (theo). All measurements are compared with predictions from various event generators.

A first observation is presented for the electroweak production of a W boson, a photon, and two jets in proton-proton collisions. The W boson decays are selected by requiring one identified electron or muon and an imbalance in transverse momentum. The two jets are required to have a high dijet mass and a large separation in pseudorapidity. The measurement is based on data collected with the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The observed (expected) significance for this process is 4.9 (4.6) standard deviations. After combining with previously reported CMS results at 8 TeV, the observed (expected) significance is 5.3 (4.8) standard deviations. The cross section for the electroweak W$\gamma$jj production in a restricted fiducial region is measured as 20.4 $\pm$ 4.5 fb and the total cross section for W$\gamma$ production in association with 2 jets in the same fiducial region is 108 $\pm$ 16 fb. All results are in good agreement with recent theoretical predictions. Constraints are placed on anomalous quartic gauge couplings in terms of dimension-8 effective field theory operators.

The first observation of the production of W$^\pm$W$^\pm$ bosons from double parton scattering processes using same-sign electron-muon and dimuon events in proton-proton collisions is reported. The data sample corresponds to an integrated luminosity of 138 fb$^{-1}$ recorded at a center-of-mass energy of 13 TeV using the CMS detector at the CERN LHC. Multivariate discriminants are used to distinguish the signal process from the main backgrounds. A binned maximum likelihood fit is performed to extract the signal cross section. The measured cross section for production of same-sign W bosons decaying leptonically is 80.7 $\pm$ 11.2 (stat) $^{+9.5}_{-8.6}$ (syst) $\pm$ 12.1 (model) fb, whereas the measured fiducial cross section is 6.28 $\pm$ 0.81 (stat) $\pm$ 0.69 (syst) $\pm$ 0.37 (model) fb. The observed significance of the signal is 6.2 standard deviations above the background-only hypothesis.

Double-parton scattering is investigated using events with a Z boson and jets. The Z boson is reconstructed using only the dimuon channel. The measurements are performed with proton-proton collision data recorded by the CMS experiment at the LHC at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$ collected in the year 2016. Differential cross sections of Z + $\geq$ 1 jet and Z + $\geq$ 2 jets are measured with transverse momentum of the jets above 20 GeV and pseudorapidity $|\eta|$$\lt$ 2.4. Several distributions with sensitivity to double-parton scattering effects are measured as functions of the angle and the transverse momentum imbalance between the Z boson and the jets. The measured distributions are compared with predictions from several event generators with different hadronization models and different parameter settings for multiparton interactions. The measured distributions show a dependence on the hadronization and multiparton interaction simulation parameters, and are important input for future improvements of the simulations.

A measurement of jet substructure observables is presented using \ttbar events in the lepton+jets channel from proton-proton collisions at $\sqrt{s}=$ 13 TeV recorded by the CMS experiment at the LHC, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Multiple jet substructure observables are measured for jets identified as bottom, light-quark, and gluon jets, as well as for inclusive jets (no flavor information). The results are unfolded to the particle level and compared to next-to-leading-order predictions from POWHEG interfaced with the parton shower generators PYTHIA 8 and HERWIG 7, as well as from SHERPA 2 and DIRE2. A value of the strong coupling at the Z boson mass, $\alpha_S(m_\mathrm{Z}) = $ 0.115$^{+0.015}_{-0.013}$, is extracted from the substructure data at leading-order plus leading-log accuracy.