The collective behavior of K$^0_\mathrm{S}$ and $\Lambda/\bar{\Lambda}$ strange hadrons is studied by measuring the elliptic azimuthal anisotropy ($v_2$) using the scalar-product and multiparticle correlation methods. Proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV collected by the CMS experiment at the LHC are investigated. Nonflow effects in the pPb collisions are studied by using a subevent cumulant analysis and by excluding events where a jet with transverse momentum greater than 20\GeV is present. The strange hadron $v_2$ values extracted in \pPb collisions via the four- and six-particle correlation method are found to be nearly identical, suggesting the collective behavior. Comparisons of the pPb and PbPb results for both strange hadrons and charged particles illustrate how event-by-event flow fluctuations depend on the system size.

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.

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

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.

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.

Multiparticle azimuthal correlations of prompt D$^0$ mesons are measured in PbPb collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. For the first time, a four-particle cumulant method is used to extract the second Fourier coefficient of the azimuthal distribution ($v_2$) of D$^0$ mesons as a function of event centrality and the D$^0$ transverse momentum. The ratios of the four-particle $v_2$ values to previously measured two-particle cumulant results provide direct experimental access to event-by-event fluctuations of charm quark azimuthal anisotropies. These ratios are also found to be comparable to those of inclusive charged particles in the event. However, hints of deviations are seen in the most central and peripheral collisions. To investigate the origin of flow fluctuations in the charm sector, these measurements are compared with models implementing fluctuations of charm quark energy loss via collisional or radiative processes in the quark-gluon plasma. These models cannot quantitatively describe the data over the full transverse momentum and centrality ranges, although the calculations with collisional energy loss provide a better description of the data.

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.

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 combination of searches for top squark pair production using proton-proton collision data at a center-of-mass energy of 13 TeV at the CERN LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$ collected by the CMS experiment, is presented. Signatures with at least 2 jets and large missing transverse momentum are categorized into events with 0, 1, or 2 leptons. New results for regions of parameter space where the kinematical properties of top squark pair production and top quark pair production are very similar are presented. Depending on the model, the combined result excludes a top squark mass up to 1325 GeV for a massless neutralino, and a neutralino mass up to 700 GeV for a top squark mass of 1150 GeV. Top squarks with masses from 145 to 295 GeV, for neutralino masses from 0 to 100 GeV, with a mass difference between the top squark and the neutralino in a window of 30 GeV around the mass of the top quark, are excluded for the first time with CMS data. The results of theses searches are also interpreted in an alternative signal model of dark matter production via a spin-0 mediator in association with a top quark pair. Upper limits are set on the cross section for mediator particle masses of up to 420 GeV.

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.

The $B_\mathrm{c}^+$ meson is observed for the first time in heavy ion collisions. Data from the CMS detector are used to study the production of the $B_\mathrm{c}^+$ meson in lead-lead (PbPb) and proton-proton (pp) collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV, via the $B_\mathrm{c}^+ \to (J/\psi\to\mu^+\mu^-)\mu^+\nu_\mu$ decay. The $B_\mathrm{c}^+$ nuclear modification factor, derived from the PbPb-to-pp ratio of production cross sections, is measured in two bins of the trimuon transverse momentum and of the PbPb collision centrality. The B$_\mathrm{c}^+$ meson is shown to be less suppressed than quarkonia and most of the open heavy-flavor mesons, suggesting that effects of the hot and dense nuclear matter created in heavy ion collisions contribute to its production. This measurement sets forth a promising new probe of the interplay of suppression and enhancement mechanisms in the production of heavy-flavor mesons in the quark-gluon plasma.

The top quark pair production cross section is measured in proton-proton collisions at a center-of-mass energy of 5.02 TeV. The data were collected in a special LHC low-energy and low-intensity run in 2017, and correspond to an integrated luminosity of 302 pb$^{-1}$. The measurement is performed using events with one electron and one muon of opposite charge, and at least two jets. The measured cross section is 60.7 $\pm$ 5.0 (stat) $\pm$ 2.8 (syst) $\pm$ 1.1 (lumi) pb. To reduce the statistical uncertainty, a combination with the result in the single lepton + jets channel, based on data collected in 2015 at the same center-of-mass energy and corresponding to an integrated luminosity of 27.4 pb$^{-1}$, is then performed. The resulting measured value is 63.0 $\pm$ 4.1 (stat) $\pm$ 3.0 (syst+lumi) pb, in agreement with the standard model prediction of 66.8 $^{+2.9}_{-3.1}$ pb.

Results are presented from a search for charged-lepton flavor violating (CLFV) interactions in top quark production and decay in pp collisions at a center-of-mass energy of 13 TeV. The events are required to contain one oppositely charged electron-muon pair in the final state, along with at least one jet identified as originating from a bottom quark. The data correspond to an integrated luminosity of 138 fb$^{-1}$, collected by the CMS experiment at the LHC. This analysis includes both the production (q $\to$ e$\mu$t) and decay (t $\to$ e$\mu$q) modes of the top quark through CLFV interactions, with q referring to a u or c quark. These interactions are parametrized using an effective field theory approach. With no significant excess over the standard model expectation, the results are interpreted in terms of vector-, scalar-, and tensor-like CLFV four-fermion effective interactions. Finally, observed exclusion limits are set at 95% confidence levels on the respective branching fractions of a top quark to an e$\mu$ pair and an up (charm) quark of 0.13 $\times$ 10$^{-6}$ (1.31 $\times$ 10$^{-6}$), 0.07 $\times$ 10$^{-6}$ (0.89 $\times$ 10$^{-6}$), and 0.25 $\times$ 10$^{-6}$ (2.59 $\times$ 10$^{-6}$) for vector, scalar, and tensor CLFV interactions, respectively.

The first measurements of diboson production cross sections in proton-proton interactions at a center-of-mass energy of 5.02 TeV are reported. They are based on data collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 302 pb$^{-1}$. Events with two, three, or four charged light leptons (electrons or muons) in the final state are analyzed. The WW, WZ, and ZZ total cross sections are measured as $\sigma_\mathrm{WW} =$ 37.0 $^{+5.5}_{-5.2}$ (stat) $^{+2.7}_{-2.6}$ (syst) pb, $\sigma_\mathrm{WZ} =$ 6.4 $^{+2.5}_{-2.1}$ (stat) $^{+0.5}_{-0.3}$ (syst) pb, and $\sigma_\mathrm{ZZ} =$ 5.3 $^{+2.5}_{-2.1}$ (stat) $^{+0.5}_{-0.4}$ (syst) pb. All measurements are in good agreement with theoretical calculations at combined next-to-next-to-leading order quantum chromodynamics and next-to-leading order electroweak accuracy.

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.

Multijet events at large transverse momentum ($p_\mathrm{T}$) are measured at $\sqrt{s}$ = 13 TeV using data recorded with the CMS detector at the LHC, corresponding to an integrated luminosity of 36.3 fb$^{-1}$. The multiplicity of jets with $p_\mathrm{T}$$>$ 50 GeV that are produced in association with a high-$p_\mathrm{T}$ dijet system is measured in various ranges of the $p_\mathrm{T}$ of the jet with the highest transverse momentum and as a function of the azimuthal angle difference $\Delta\phi_{1,2}$ between the two highest $p_\mathrm{T}$ jets in the dijet system. The differential production cross sections are measured as a function of the transverse momenta of the four highest $p_\mathrm{T}$ jets. The measurements are compared with leading and next-to-leading order matrix element calculations supplemented with simulations of parton shower, hadronization, and multiparton interactions. In addition, the measurements are compared with next-to-leading order matrix element calculations combined with transverse-momentum dependent parton densities and transverse-momentum dependent parton shower.

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.

A search for long-lived particles (LLPs) produced in decays of standard model (SM) Higgs bosons is presented. The data sample consists of 137 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} =$ 13 TeV, recorded at the LHC in 2016-2018. A novel technique is employed to reconstruct decays of LLPs in the endcap muon detectors. The search is sensitive to a broad range of LLP decay modes and to masses as low as a few GeV. No excess of events above the SM background is observed. The most stringent limits to date on the branching fraction of the Higgs boson to LLPs subsequently decaying to quarks and $\tau^+\tau^-$ are found for proper decay lengths greater than 6, 20, and 40 m, for LLP masses of 7, 15, and 40 GeV, respectively.

A search for long-lived particles (LLPs) produced in association with a Z boson is presented. The study is performed using data from proton-proton collisions with a center-of-mass energy of 13 TeV recorded by the CMS experiment during 2016-2018, corresponding to an integrated luminosity of 117 fb$^{-1}$. The LLPs are assumed to decay to a pair of standard model quarks that are identified as displaced jets within the CMS tracker system. Triggers and selections based on Z boson decays to electron or muon pairs improve the sensitivity to light LLPs (down to 15 GeV). This search provides sensitivity to beyond the standard model scenarios which predict LLPs produced in association with a Z boson. In particular, the results are interpreted in the context of exotic decays of the Higgs boson to a pair of scalar LLPs (H $\to$ SS). The Higgs boson decay branching fraction is constrained to values less than 6% for proper decay lengths of 10-100 mm and for LLP masses between 40 and 55 GeV. In the case of low-mass ($\approx$ 15 GeV) scalar particles that subsequently decay to a pair of b quarks, the search is sensitive to branching fractions $\mathcal{B}$(H $\to$ SS) $\lt$ 20% for proper decay lengths of 10-50 mm. The use of associated production with a Z boson increases the sensitivity to low-mass LLPs of this analysis with respect to gluon fusion searches. In the case of 15 GeV scalar LLPs, the improvement corresponds to a factor of 2 at a proper decay length of 30 mm.

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 search for physics beyond the standard model (SM) in final states with an electron or muon and missing transverse momentum is presented. The analysis uses data from proton-proton collisions at a centre-of-mass energy of 13 TeV, collected with the CMS detector at the LHC in 2016–2018 and corresponding to an integrated luminosity of 138 fb−1. No significant deviation from the SM prediction is observed. Model-independent limits are set on the production cross section of W’ bosons decaying into lepton-plus-neutrino final states. Within the framework of the sequential standard model, with the combined results from the electron and muon decay channels a W’ boson with mass less than 5.7 TeV is excluded at 95% confidence level. Results on a SM precision test, the determination of the oblique electroweak W parameter, are presented using LHC data for the first time. These results together with those from the direct W’ resonance search are used to extend existing constraints on composite Higgs scenarios. This is the first experimental exclusion on compositeness parameters using results from LHC data other than Higgs boson measurements.

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.