A measurement of inclusive four-jet production in proton-proton collisions at a center-of-mass energy of 13\TeV is presented. The transverse momenta of jets within $\lvert\eta\rvert \lt$ 4.7 reach down to 35, 30, 25, and 20 GeV for the first-, second-, third-, and fourth-leading jet, respectively. Differential cross sections are measured as functions of the jet transverse momentum, jet pseudorapidity, and several other observables that describe the angular correlations between the jets. The measured distributions show sensitivity to different aspects of the underlying event, parton shower, and matrix element calculations. In particular, the interplay between angular correlations caused by parton shower and double-parton scattering contributions is shown to be important. The double-parton scattering contribution is extracted by means of a template fit to the data, using distributions for single-parton scattering obtained from Monte Carlo event generators and a double-parton scattering distribution constructed from inclusive single-jet events in data. The effective double-parton scattering cross section is calculated and discussed in view of previous measurements and of its dependence on the models used to describe the single-parton scattering background.

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.

Collinear (small-angle) and large-angle, as well as soft and hard radiations are investigated in three-jet and Z + two-jet events collected in proton-proton collisions at the LHC. The normalized production cross sections are measured as a function of the ratio of transverse momenta of two jets and their angular separation. The measurements in the three-jet and Z + two-jet events are based on data collected at a center-of-mass energy of 8 TeV, corresponding to an integrated luminosity of 19.8 fb$^{-1}$. The Z + two-jet events are reconstructed in the dimuon decay channel of the Z boson. The three-jet measurement is extended to include $\sqrt{s} =$ 13 TeV data corresponding to an integrated luminosity of 2.3 fb$^{-1}$. The results are compared to predictions from event generators that include parton showers, multiple parton interactions, and hadronization. The collinear and soft regions are in general well described by parton showers, whereas the regions of large angular separation are often best described by calculations using higher-order matrix elements.

A measurement of the production cross section of top quark pairs in association with two b jets ($\mathrm{t\bar{t}}\mathrm{b\bar{b}}$) is presented using data collected in proton-proton collisions at $\sqrt{s} =$ 13 TeV by the CMS detector at the LHC corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The cross section is measured in the all-jet decay channel of the top quark pair by selecting events containing at least eight jets, of which at least two are identified as originating from the hadronization of b quarks. A combination of multivariate analysis techniques is used to reduce the large background from multijet events not containing a top quark pair, and to help discriminate between jets originating from top quark decays and other additional jets. The cross section is determined for the total phase space to be 5.5 $\pm$ 0.3 (stat)${}^{+1.6}_{-1.3}$ (syst) pb and also measured for two fiducial $\mathrm{t\bar{t}}\mathrm{b\bar{b}}$ definitions. The measured cross sections are found to be larger than theoretical predictions by a factor of 1.5-2.4, corresponding to 1-2 standard deviations.

Ratios of cross sections, $\sigma$(Z+c jets)/$\sigma$(Z+jets), $\sigma$(Z+b jets)/$\sigma$(Z+jets), and $\sigma$(Z+c jets)/$\sigma$(Z+b jets) in the associated production of a Z boson with at least one charm or bottom quark jet are measured in proton-proton collisions at $\sqrt{s}=$ 13 TeV. The data sample, collected by the CMS experiment at the CERN LHC, corresponds to an integrated luminosity of 35.9 fb$^{-1}$, with a fiducial volume of $p_\mathrm{T}>$ 30 GeV and $|\eta|<$ 2.4 for the jets, where $p_\mathrm{T}$ and $\eta$ represent transverse momentum and pseudorapidity, respectively. The Z boson candidates come from leptonic decays into electrons or muons with $p_\mathrm{T}>$ 25 GeV and $|\eta|<$ 2.4, and the dilepton mass satisfies 71 $\lt m_\mathrm{Z}\lt$ 111 GeV. The measured values are $\sigma$(Z+c jets)/$\sigma$(Z+jets) = 0.102 $\pm$ 0.002 $\pm$ 0.009, $\sigma$(Z+b jets)/$\sigma$(Z+jets) = 0.0633 $\pm$ 0.0004 $\pm$ 0.0015, and $\sigma$(Z+c jets)/$\sigma$(Z+b jets) = 1.62 $\pm$ 0.03 $\pm$ 0.15. Results on the inclusive and differential cross section ratios as functions of jet and Z boson transverse momentum are compared with predictions from leading and next-to-leading order perturbative quantum chromodynamics calculations. These are the first measurements of the cross section ratios at 13 TeV, and the measurement precision is better than that of the current theoretical predictions.

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.

Differential cross section measurements of W$^\pm\gamma$ production in proton-proton collisions at $\sqrt{s}$ = 13 TeV are presented. The data set used in this study was collected with the CMS detector at the CERN LHC in 2016-2018 with an integrated luminosity of 138 fb$^{-1}$. Candidate events containing an electron or muon, a photon, and missing transverse momentum are selected. The measurements are compared with standard model predictions computed at next-to-leading and next-to-next-to-leading orders in perturbative quantum chromodynamics. Constraints on the presence of TeV-scale new physics affecting the WW$\gamma$ vertex are determined within an effective field theory framework, focusing on the $\mathcal{O}_\mathrm{3W}$ operator. A simultaneous measurement of the photon transverse momentum and the azimuthal angle of the charged lepton in a special reference frame is performed. This two-dimensional approach provides up to a factor of ten more sensitivity to the interference between the standard model and the $\mathcal{O}_\mathrm{3W}$ contribution than using the transverse momentum alone.

The B$^0_s$ and B$^+$ production yields are measured in PbPb collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data sample, collected with the CMS detector at the LHC, corresponds to an integrated luminosity of 1.7 nb$^{-1}$. The mesons are reconstructed in the exclusive decay channels B$^0_s$$\to$ J/$\psi(\mu^+\mu^-)\phi($K$^+$K$^-)$ and B$^+$$\to$ J/$\psi(\mu^+\mu^-)$K$^+$, in the transverse momentum range 7-50 GeV/c and absolute rapidity 0-2.4. The B$^0_s$ meson is observed with a statistical significance in excess of five standard deviations for the first time in nucleus-nucleus collisions. The measurements are performed as functions of the transverse momentum of the B mesons and of the PbPb collision centrality. The ratio of production yields of B$^0_s$ and B$^+$ is measured and compared to theoretical models that include quark recombination effects.

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.