A search for a heavy Higgs boson in the H to WW and H to ZZ decay channels is reported. The search is based upon proton-proton collision data samples corresponding to an integrated luminosity of up to 5.1 inverse femtobarns at sqrt(s) = 7 TeV and up to 19.7 inverse femtobarns at sqrt(s) = 8 TeV, recorded by the CMS experiment at the CERN LHC. Several final states of the H to WW and H to ZZ decays are analyzed. The combined upper limit at the 95% confidence level on the product of the cross section and branching fraction exclude a Higgs boson with standard model-like couplings and decays in the range 145 < m[H] < 1000 GeV. We also interpret the results in the context of an electroweak singlet extension of the standard model.

A search for neutral Higgs bosons predicted in the minimal supersymmetric standard model (MSSM) for mu+ mu- decay channels is presented. The analysis uses data collected by the CMS experiment at the LHC in proton-proton collisions at centre-of-mass energies of 7 and 8 TeV, corresponding to integrated luminosities of 5.1 and 19.3 inverse femtobarns, respectively. The search is sensitive to Higgs bosons produced through the gluon fusion process or in association with a bb quark pair. No statistically significant excess is observed in the mu+ mu- mass spectrum. Results are interpreted in the framework of several benchmark scenarios, and the data are used to set an upper limit on the MSSM parameter tan(beta) as a function of the mass of the pseudoscalar A boson in the range from 115 to 300 GeV. Model independent upper limits are given for the product of the cross section and branching fraction for gluon fusion and b quark associated production. They are the most stringent limits obtained to date in this channel.

A measurement of the underlying event (UE) activity in proton-proton collisions is performed using events with charged-particle jets produced in the central pseudorapidity region (abs(eta[jet]) < 2) and with transverse momentum 1<= pt[jet] < 100 GeV. The analysis uses a data sample collected at a centre-of-mass energy of 2.76 TeV with the CMS experiment at the LHC. The UE activity is measured as a function of pt[jet] in terms of the average multiplicity and scalar sum of transverse momenta (pt) of charged particles, with abs(eta) < 2 and pt > 0.5 GeV, in the azimuthal region transverse to the highest pt jet direction. By further dividing the transverse region into two regions of smaller and larger activity, various components of the UE activity are separated. The measurements are compared to previous results at 0.9 and 7 TeV, and to predictions of several Monte Carlo event generators, providing constraints on the modelling of the UE dynamics.

A search is described for a Higgs boson decaying into two photons, one of which has an internal conversion to a muon or an electron pair (ll gamma). The analysis is performed using proton-proton collision data recorded with the CMS detector at the LHC at a centre-of-mass energy of 8 TeV, corresponding to an integrated luminosity of 19.7 inverse femtobarns. The events selected have an opposite-sign muon or electron pair and a high transverse momentum photon. No excess above background has been found in the three-body invariant mass range 120 < m[ll gamma] < 150 GeV, and limits have been derived for the Higgs boson production cross section times branching fraction for the decay H to gamma* gamma to ll gamma, where the dilepton invariant mass is less than 20 GeV. For a Higgs boson with m[H] = 125 GeV, a 95% confidence level (CL) exclusion observed (expected) limit is 6.7 (5.9 +2.8/-1.8) times the standard model prediction. Additionally, an upper limit at 95% CL on the branching fraction of H to J/Psi gamma for the 125 GeV Higgs boson is set at 1.5E-3.

A search for pairs of Higgs bosons produced via gluon and vector boson fusion is presented, focusing on the four b quark final state. The data sample consists of proton-proton collisions at a center-of-mass energy of 13 TeV, collected with the CMS detector at the LHC, and corresponds to an integrated luminosity of 138 fb$^{-1}$. No deviation from the background-only hypothesis is observed. A 95% confidence level upper limit on the Higgs boson pair production cross section is observed at 3.9 times the standard model prediction for an expected value of 7.8. Constraints are also set on the modifiers of the Higgs field self-coupling, $\kappa_\lambda$, and of the coupling of two Higgs bosons to two vector bosons, $\kappa_\mathrm{2V}$. The observed (expected) allowed intervals at the 95% confidence level are $-$2.3 $\lt$ $\kappa_\lambda$ $\lt$ 9.4 ($-$5.0 $\lt$ $\kappa_\lambda$ $\lt$ 12.0) and $-$0.1 $\lt$ $\kappa_\mathrm{2V}$ $\lt$ 2.2 ($-$0.4 $\lt$ $\kappa_\mathrm{2V}$ $\lt$ 2.5). These are the most stringent observed constraints to date on the HH production cross section and on the $\kappa_\mathrm{2V}$ coupling.

A search for flavor-changing neutral current interactions of the top quark (t) and the Higgs boson (H) is presented. The search is based on a data sample corresponding to an integrated luminosity of 137 fb$^{-1}$ recorded by the CMS experiment at the LHC in proton-proton collisions at $\sqrt{s}$ = 13 TeV. Events containing exactly one lepton (muon or electron) and at least three jets, among which at least two are identified as originating from the hadronization of a bottom quark, are analyzed. A set of deep neural networks is used for kinematic event reconstruction, while boosted decision trees distinguish the signal from the background events. No significant excess over the background predictions is observed, and upper limits on the signal production cross sections are extracted. These limits are interpreted in terms of top quark decay branching fractions ($\mathcal{B}$) to the Higgs boson and an up (u) or a charm quark (c). Assuming one nonvanishing extra coupling at a time, the observed (expected) upper limits at 95% confidence level are $\mathcal{B}$(t $\to$ Hu) $\lt$ 0.079 (0.11)% and $\mathcal{B}$(t $\to$ Hc) $\lt$ 0.094 (0.086)%.

Inclusive and differential cross sections of single top quark production in association with a Z boson are measured in proton-proton collisions at a center-of-mass energy of 13 TeV with a data sample corresponding to an integrated luminosity of 138 fb$^{-1}$ recorded by the CMS experiment. Events are selected based on the presence of three leptons, electrons or muons, associated with leptonic Z boson and top quark decays. The measurement yields an inclusive cross section of 87.9 $_{-7.3}^{+7.5}$ (stat) $_{-6.0}^{+7.3}$ (syst) fb for a dilepton invariant mass greater than 30 GeV, in agreement with standard model (SM) calculations and the most precise determination to date. The ratio between the cross sections for the top quark and the top antiquark production in association with a Z boson is measured as 2.37 $_{-0.42}^{+0.56}$ (stat) ${}_{-0.13}^{+0.27}$ (syst). Differential measurements at parton and particle levels are performed for the first time. Several kinematic observables are considered to study the modeling of the process. Results are compared to theoretical predictions with different assumptions on the source of the initial-state b quark and found to be in agreement, within the uncertainties. Additionally, the spin asymmetry, which is sensitive to the top quark polarization, is determined from the differential distribution of the polarization angle at parton level to be 0.54 $\pm$ 0.16 (stat) $\pm$ 0.06 (syst), in agreement with SM predictions.

A measurement of the jet mass distribution in hadronic decays of Lorentz-boosted top quarks is presented. The measurement is performed in the lepton+jets channel of top quark pair production ($\mathrm{t\bar{t}}$) events, where the lepton is an electron or muon. The products of the hadronic top quark decay are reconstructed using a single large-radius jet with transverse momentum greater than 400 GeV. The data were collected with the CMS detector at the LHC in proton-proton collisions and correspond to an integrated luminosity of 138 fb$^{-1}$. The differential $\mathrm{t\bar{t}}$ production cross section as a function of the jet mass is unfolded to the particle level and is used to extract the top quark mass. The jet mass scale is calibrated using the hadronic W boson decay within the large-radius jet. The uncertainties in the modelling of the final state radiation are reduced by studying angular correlations in the jet substructure. These developments lead to a significant increase in precision, and a top quark mass of 173.06 $\pm$ 0.84 GeV.

A search is presented for vector-like T and B quark-antiquark pairs produced in proton-proton collisions at a center-of-mass energy of 13 TeV. Data were collected by the CMS experiment at the CERN LHC in 2016-2018, with an integrated luminosity of 138 fb$^{-1}$. Events are separated into single-lepton, same-sign charge dilepton, and multilepton channels. In the analysis of the single-lepton channel a multilayer neural network and jet identification techniques are employed to select signal events, while the same-sign dilepton and multilepton channels rely on the high-energy signature of the signal to distinguish it from standard model backgrounds. The data are consistent with standard model background predictions, and the production of vector-like quark pairs is excluded at 95% confidence level for T quark masses up to 1.54 TeV and B quark masses up to 1.56 TeV, depending on the branching fractions assumed, with maximal sensitivity to decay modes that include multiple top quarks. The limits obtained in this search are the strongest limits to date for $\mathrm{T\overline{T}}$ production, excluding masses below 1.48 TeV for all decays to third generation quarks, and are the strongest limits to date for $\mathrm{B\overline{B}}$ production with B quark decays to tW.

A search has been performed for heavy resonances decaying to ZZ or ZW and for axion-like particles (ALPs) mediating nonresonant ZZ or ZH production, in final states with two charged leptons ($\ell$ = e, $\mu$) produced by the decay of a Z boson, and two quarks produced by the decay of a Z, W, or Higgs boson H. The analysis is sensitive to resonances with masses in the range 450 to 2000 GeV. Two categories are defined corresponding to the merged or resolved reconstruction of the hadronically decaying boson. The search is based on data collected during 2016-2018 by the CMS experiment at the LHC in proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. No significant excess is observed in the data above the standard model background expectation. Upper limits on the production cross section of heavy, narrow spin-2 and spin-1 resonances are derived as functions of the resonance mass, and exclusion limits on the production of bulk graviton particles and W$'$ bosons are calculated in the framework of the warped extra dimensions and heavy vector triplet models, respectively. In addition, upper limits on the ALP-mediated diboson production cross section and ALP couplings to standard model particles are obtained in the framework of linear and chiral effective field theories. These are the first limits on nonresonant ALP-mediated ZZ and ZH production obtained by the LHC experiments.

A search is presented for a right-handed W boson (W$_\mathrm{R}$) and a heavy neutrino (N), in a final state consisting of two same-flavor leptons (ee or $\mu\mu$) and two quarks. The search is performed with the CMS experiment at the CERN LHC using a data sample of proton-proton collisions at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$. The search covers two regions of phase space, one where the decay products of the heavy neutrino are merged into a single large-area jet, and one where the decay products are well separated. The expected signal is characterized by an excess in the invariant mass distribution of the final-state objects. No significant excess over the standard model background expectations is observed. The observations are interpreted as upper limits on the product of W$_\mathrm{R}$ production cross sections and branching fractions assuming that couplings are identical to those of the standard model W boson. For N masses $m_\mathrm{N}$ equal to half the W$_\mathrm{R}$ mass $m_\mathrm{W_R}$ ($m_\mathrm{N}$ = 0.2 TeV), $m_\mathrm{W_R}$ is excluded at 95% confidence level up to 4.7 (4.8) and 5.0 (5.4) TeV for the electron and muon channels, respectively. This analysis provides the most stringent limits on the W$_\mathrm{R}$ mass to date.

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.

Results are presented from a search for physics beyond the standard model in proton-proton collisions at $\sqrt{s} =$ 13 TeV in channels with two Higgs bosons, each decaying via the process H $\to$$\mathrm{b\bar{b}}$, and large missing transverse momentum. The search uses a data sample corresponding to an integrated luminosity of 137 fb$^{-1}$ collected by the CMS experiment at the CERN LHC. The search is motivated by models of supersymmetry that predict the production of neutralinos, the neutral partners of the electroweak gauge and Higgs bosons. The observed event yields in the signal regions are found to be consistent with the standard model background expectations. The results are interpreted using simplified models of supersymmetry. For the electroweak production of nearly mass-degenerate higgsinos, each of whose decay chains yields a neutralino ($\tilde{\chi}^0_1$) that in turn decays to a massless goldstino and a Higgs boson, $\tilde{\chi}^0_1$ masses in the range 175 to 1025 GeV are excluded at 95% confidence level. For the strong production of gluino pairs decaying via a slightly lighter $\tilde{\chi}^0_2$ to H and a light $\tilde{\chi}^0_1$, gluino masses below 2330 GeV are excluded.

A search for new long-lived particles decaying to leptons using proton-proton collision data produced by the CERN LHC at $\sqrt{s}$ = 13 TeV is presented. Events are selected with two leptons (an electron and a muon, two electrons, or two muons) that both have transverse impact parameter values between 0.01 and 10 cm and are not required to form a common vertex. Data used for the analysis were collected with the CMS detector in 2016, 2017, and 2018, and correspond to an integrated luminosity of 118 (113) fb$^{-1}$ in the ee channel (e$\mu$ and $\mu\mu$ channels). The search is designed to be sensitive to a wide range of models with displaced e$\mu$, ee, and $\mu\mu$ final states. The results constrain several well-motivated models involving new long-lived particles that decay to displaced leptons. For some areas of the available phase space, these are the most stringent constraints to date.

A measurement is presented for the electroweak production of a W boson, a photon ($\gamma$), and two jets (j) in proton-proton collisions. The leptonic decay of the W boson is selected by requiring one identified electron or muon and large missing transverse momentum. The two jets are required to have large invariant dijet mass and large separation in pseudorapidity. The measurement is performed with the data collected by the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The cross section for the electroweak W$\gamma$jj production is 23.5 $^{+4.9}_{-4.7}$ fb, whereas the total cross section for W$\gamma$jj production is 113 $\pm$ 13 fb. Differential cross sections are also measured with the distributions unfolded to the particle level. All results are in agreement with the standard model expectations. Constraints are placed on anomalous quartic gauge couplings (aQGCs) in terms of dimension-8 effective field theory operators. These are the most stringent limits to date on the aQGCs parameters $f_\mathrm{M,2-5}$$/$$\Lambda^4$ and $f_\mathrm{T,6-7}$$/$$\Lambda^4$.

A measurement of the inclusive jet production in proton-proton collisions at the LHC at $\sqrt{s}$ = 13 TeV is presented. The double-differential cross sections are measured as a function of the jet transverse momentum $p_\mathrm{T}$ and the absolute jet rapidity $\lvert y \rvert$. The anti-$k_\mathrm{T}$ clustering algorithm is used with distance parameter of 0.4 (0.7) in a phase space region with jet $p_\mathrm{T}$ from 97 GeV up to 3.1 TeV and $\lvert y \rvert\lt$ 2.0. Data collected with the CMS detector are used, corresponding to an integrated luminosity of 36.3 fb$^{-1}$ (33.5 fb$^{-1}$). The measurement is used in a comprehensive QCD analysis at next-to-next-to-leading order, which results in significant improvement in the accuracy of the parton distributions in the proton. Simultaneously, the value of the strong coupling constant at the Z boson mass is extracted as $\alpha_\mathrm{S}$(Z) = 0.1170 $\pm$ 0.0019. For the first time, these data are used in a standard model effective field theory analysis at next-to-leading order, where parton distributions and the QCD parameters are extracted simultaneously with imposed constraints on the Wilson coefficient $c_1$ of 4-quark contact interactions. Note added: in the Addendum to this paper, available as Appendix B in this document, an improved value of $\alpha_\mathrm{S}$(Z) = 0.1166 $\pm$ 0.0017 has been extracted. This result supersedes the number in the above abstract of the original publication.

The mass of the top quark is measured in 36.3 fb$^{-1}$ of LHC proton-proton collision data collected with the CMS detector at $\sqrt{s}$ = 13 TeV. The measurement uses a sample of top quark pair candidate events containing one isolated electron or muon and at least four jets in the final state. For each event, the mass is reconstructed from a kinematic fit of the decay products to a top quark pair hypothesis. A profile likelihood method is applied using up to four observables to extract the top quark mass. The top quark mass is measured to be 171.77 $\pm$ 0.37 GeV. This approach significantly improves the precision over previous measurements.

The second-order ($v_2$) and third-order ($v_3$) Fourier coefficients describing the azimuthal anisotropy of prompt and nonprompt (from b-hadron decays) J/$\psi$, as well as prompt $\psi$(2S) mesons are measured in lead-lead collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The analysis uses a data set corresponding to an integrated luminosity of 1.61 nb$^{-1}$ recorded with the CMS detector. The J/$\psi$ and $\psi$(2S) mesons are reconstructed using their dimuon decay channel. The $v_2$ and $v_3$ coefficients are extracted using the scalar product method and studied as functions of meson transverse momentum and collision centrality. The measured $v_2$ values for prompt J/$\psi$ mesons are found to be larger than those for nonprompt J/$\psi$ mesons. The prompt J/$\psi$$v_2$ values at high $p_\mathrm{T}$ are found to be underpredicted by a model incorporating only parton energy loss effects in a quark-gluon plasma medium. Prompt and nonprompt J/$\psi$ meson $v_3$ and prompt $\psi$(2S) $v_2$ and $v_3$ values are also reported for the first time, providing new information about heavy quark interactions in the hot and dense medium created in heavy ion collisions.

A systematic study of the factorization of long-range azimuthal two-particle correlations into a product of single-particle anisotropies is presented as a function of pt and eta of both particles, and as a function of the particle multiplicity in PbPb and pPb collisions. The data were taken with the CMS detector for PbPb collisions at sqrt(s[NN]) = 2.76 TeV and pPb collisions at sqrt(s[NN]) = 5.02 TeV, covering a very wide range of multiplicity. Factorization is observed to be broken as a function of both particle pt and eta. When measured with particles of different pt, the magnitude of the factorization breakdown for the second Fourier harmonic reaches 20% for very central PbPb collisions but decreases rapidly as the multiplicity decreases. The data are consistent with viscous hydrodynamic predictions, which suggest that the effect of factorization breaking is mainly sensitive to the initial-state conditions rather than to the transport properties (e.g., shear viscosity) of the medium. The factorization breakdown is also computed with particles of different eta. The effect is found to be weakest for mid-central PbPb events but becomes larger for more central or peripheral PbPb collisions, and also for very high-multiplicity pPb collisions. The eta-dependent factorization data provide new insights to the longitudinal evolution of the medium formed in heavy ion collisions.

A search for the standard model Higgs boson decaying to a charm quark-antiquark pair, H $\to$$\mathrm{c\bar{c}}$, produced in association with a leptonically decaying V (W or Z) boson is presented. The search is performed with proton-proton collisions at $\sqrt{s}$ = 13 TeV collected by the CMS experiment, corresponding to an integrated luminosity of 138 fb$^{-1}$. Novel charm jet identification and analysis methods using machine learning techniques are employed. The analysis is validated by searching for Z $\to$$\mathrm{c\bar{c}}$ in VZ events, leading to its first observation at a hadron collider with a significance of 5.7 standard deviations. The observed (expected) upper limit on $\sigma$(VH)$ \mathcal{B}$(H $\to$$\mathrm{c\bar{c}}$) is 0.94 (0.50 $^{+0.22}_{-0.15}$) pb at 95% confidence level (CL), corresponding to 14 (7.6 $^{+3.4}_{-2.3}$) times the standard model prediction. For the Higgs-charm Yukawa coupling modifier, $\kappa_\mathrm{c}$, the observed (expected) 95% CL interval is 1.1 $\lt$ $\vert\kappa_\mathrm{C}\vert$ $\lt$ 5.5 ($\vert\kappa_\mathrm{c}\vert$ $\lt$ 3.4), the most stringent constraint to date.

The production of a Z boson, decaying into two leptons and produced in association with one or more b jets, is studied using proton-proton collisions delivered by the LHC at a centre-of-mass energy of 7 TeV. The data were recorded in 2011 with the CMS detector and correspond to an integrated luminosity of 5 inverse femtobarns. The Z(ll) + b-jets cross sections (where ll = mu mu or ee) are measured separately for a Z boson produced with exactly one b jet and with at least two b jets. In addition, a cross section ratio is extracted for a Z boson produced with at least one b jet, relative to a Z boson produced with at least one jet. The measured cross sections are compared to various theoretical predictions, and the data favour the predictions in the five-flavour scheme, where b quarks are assumed massless. The kinematic properties of the reconstructed particles are compared with the predictions from the MADGRAPH event generator using the PYTHIA parton shower simulation.

The ratio of the top-quark branching fractions $R = B(t \to Wb)/B(t \to Wq)$, where the denominator includes the sum over all down-type quarks (q = b, s, d), is measured in the $t\bar{t}$ dilepton final state with proton-proton collision data at $\sqrt{s}$ = 8 TeV from an integrated luminosity of 19.7 inverse-femtobarns, collected with the CMS detector. In order to quantify the purity of the signal sample, the cross section is measured by fitting the observed jet multiplicity, thereby constraining the signal and background contributions. By counting the number of b jets per event, an unconstrained value of R = 1.014 $\pm$ 0.003 (stat) $\pm$ 0.032 (syst) is measured, in good agreement with the standard model prediction. A lower limit R greater than 0.955 at the 95% confidence level is obtained after requiring R lower than one, and a lower limit on the Cabibbo-Kobayashi-Maskawa matrix element |$V_tb$| greater than 0.975 is set at 95% confidence level. The result is combined with a previous CMS measurement of the t-channel single-top-quark cross section to determine the top-quark total decay width, $\Gamma_t$ = 1.36 $\pm$ 0.02 (stat)$^{+0.14}_{-0.11}$ (syst) GeV.

Measurements of $W\gamma$ and $Z\gamma$ production in proton-proton collisions at $\sqrt{s}$ = 7 TeV are used to extract limits on anomalous triple gauge couplings. The results are based on data recorded by the CMS experiment at the LHC that correspond to an integrated luminosity of 5.0 inverse femtobarns. The cross sections are measured for photon transverse momenta $p_T^{\gamma} \gt$ 15 GeV, and for separations between photons and final-state charged leptons in the pseudorapidity-azimuthal plane of $\Delta R(l, \gamma) \gt$ 0.7 in $l \nu \gamma$ and $ll \gamma$ final states, where l refers either to an electron or a muon. A dilepton invariant mass requirement of $m_{ll} \gt$ 50 GeV is imposed for the Z$\gamma$ process. No deviations are observed relative to predictions from the standard model, and limits are set on anomalous WW$\gamma$, ZZ$\gamma$, and Z$\gamma\gamma$ triple gauge couplings.

Characteristics of multi-particle production in proton-proton collisions at $\sqrt{s}$=7 TeV are studied as a function of the charged-particle multiplicity, $N_{ch}$. The produced particles are separated into two classes: those belonging to jets and those belonging to the underlying event. Charged particles are measured with pseudorapidity |η|<2.4 and transverse momentum $p_T$ > 0.25 GeV/c. Jets are reconstructed from charged-particles only and required to have $p_T$ > 5 GeV/c. The distributions of jet $p_T$, average $p_T$ of charged particles belonging to the underlying event or to jets, jet rates, and jet shapes are presented as functions of $N_{ch}$ and compared to the predictions of the PYTHIA and HERWIG event generators. Predictions without multi-parton interactions fail completely to describe the $N_{ch}$-dependence observed in the data. For increasing $N_{ch}$, PYTHIA systematically predicts higher jet rates and harder $p_T$ spectra than seen in the data, whereas HERWIG shows the opposite trends. At the highest multiplicity, the data–model agreement is worse for most observables, indicating the need for further tuning and/or new model ingredients.