A portrait of the Higgs boson by the CMS experiment ten years after the discovery

In July 2012, the ATLAS and CMS collaborations at the CERN Large Hadron Collider announced the observation of a Higgs boson at a mass of around 125 gigaelectronvolts. Ten years later, and with the data corresponding to the production of a 30-times larger number of Higgs bosons, we have learnt much more about the properties of the Higgs boson. The CMS experiment has observed the Higgs boson in numerous fermionic and bosonic decay channels, established its spin–parity quantum numbers, determined its mass and measured its production cross-sections in various modes. Here the CMS Collaboration reports the most up-to-date combination of results on the properties of the Higgs boson, including the most stringent limit on the cross-section for the production of a pair of Higgs bosons, on the basis of data from proton–proton collisions at a centre-of-mass energy of 13 teraelectronvolts. Within the uncertainties, all these observations are compatible with the predictions of the standard model of elementary particle physics. Much evidence points to the fact that the standard model is a low-energy approximation of a more comprehensive theory. Several of the standard model issues originate in the sector of Higgs boson physics. An order of magnitude larger number of Higgs bosons, expected to be examined over the next 15 years, will help deepen our understanding of this crucial sector.


1
The established theory of elementary particle physics, commonly referred to as the standard model (SM), provides a complete description of the electromagnetic (EM), weak, and strong interactions of matter particles, which are spin-1/2 fermions, through three different sets of mediators, which are spin-1 bosons. (In quantum mechanics, spin is an intrinsic form of angular momentum carried by elementary particles.) These vector bosons are the massless photons (gluons) for the EM (strong) interaction, and the heavy W and Z bosons for the weak interaction. The SM has been very successful in providing accurate predictions for essentially all particle physics experiments carried out so far. In 2012, the final missing particle of the SM, the Higgs boson, was observed by the ATLAS [1] and CMS [2,3] Collaborations at CERN. The Higgs boson is a prediction of a mechanism that took place in the early universe, less than a picosecond after the Big Bang, which led to the EM and the weak interactions becoming distinct in their actions. In the SM, this mechanism, labelled as the Brout-Englert-Higgs (BEH) mechanism, introduces a complex scalar (spin-0) field that permeates the entire universe. Its quantum manifestation is known as the SM Higgs boson. Scalar fields are described only by a number at every point in space that is invariant under Lorentz transformations. An analogy can be drawn of a map of an area where temperature is shown at various positions mimicking a scalar field. The same map, where instead the wind speed and direction are shown, would correspond to a vector field.

The long road to the Higgs boson
The BEH mechanism was first proposed in 1964 in the works of Brout and Englert [4], Higgs [5,6], and Guralnik, Hagen, and Kibble [7]. Further details of the mechanism were presented in 1966 by Higgs [8] and in 1967 by Kibble [9]. In 1967, Weinberg [10] and Salam [11], extending the 1961 work of Glashow [12], proposed the use of the BEH mechanism for a theory of the unification of the EM and weak interactions, labelled as the electroweak (EW) interaction. The key element in this work was the conjecture that nature possesses an EW symmetry, mathematically described by the Lagrangian of the theory, which is spontaneously broken, granting mass to the W and Z bosons. An additional feature of this model is that it provides a mechanism for granting masses to fermions as well, through the so-called Yukawa interactions [10,13]. Thus, the elementary particles interacting with the BEH field acquire mass. The impact is far reaching: for example, electrons become massive, allowing atoms to form, and endowing our universe with the observed complexity.
Salam and Weinberg had further conjectured that the model they put forward might be renormalizable (that is, give finite answers). In 1971, 't Hooft and Veltman [14,15] showed how indeed this theory could be renormalized. This development put the Glashow-Salam-Weinberg model on a firm basis deserving serious experimental scrutiny.
After the W and Z bosons were discovered by the UA1 and UA2 experiments at CERN in 1983 [16][17][18][19], the search for the Higgs boson became a central thrust in particle physics and an important motivation for the CERN Large Hadron Collider (LHC) [20], and the ATLAS and CMS experiments.
Finding the Higgs boson has been demanding. This is a consequence of its large mass, which puts it beyond the reach of previous electron-positron colliders, such as LEP [21] at CERN, and low cross section modes coupled with unfavourable decay channels in the range of mass in which it was eventually found, which made it challenging to observe at previous hadron colliders, such as the Tevatron [22] at Fermilab. In the SM, the Higgs boson is an elementary scalar particle, a type that had never been observed before. Fundamental scalar particles are subject to quantum corrections that can be as large as the scale of the physics beyond the SM (BSM). As this scale can be many orders of magnitude larger than the EW scale, which is ∼ 100 GeV, the measured mass of the Higgs boson is puzzlingly small. How to resolve this puzzle is part of the motivation for future work and accelerators.
The BEH mechanism does not predict the mass of the Higgs boson, but once the mass is fixed, all its other properties are precisely defined. The Higgs boson, once produced, decays directly to the heaviest allowed elementary particles. However, decays to massless particles can also occur through quantum loops. At the LHC, the production of Higgs bosons is dominated by gluon-gluon fusion (ggH) proceeding via a virtual top quark loop. The mass of a real particle is defined as m 2 = E 2 − p 2 , where E is the energy and p is the momentum vector of the particle. For a virtual particle this equation is not valid and thus a virtual particle does not have a defined value of the mass. A virtual particle is denoted by an asterisk, e.g. W * denoting a virtual W boson. Henceforth the distinction between real and virtual particles will be dropped, unless mentioned otherwise. At a mass of around 125 GeV the Higgs boson decays dominantly into a b quark and its antiquark. Henceforth, the distinction between a particle and its antiparticle will be dropped.
From the accurate observation and measurement of the products of the Higgs boson decays and of those associated with its production, experiments are able to infer its properties, including the strength of its self-interaction (λ) [23] and, potentially, decays into beyond-the-SM (BSM) particles.
This paper presents the combination of results from single Higgs boson production and decay, and its pair production, using data sets corresponding to an integrated luminosity (L) up to 138 fb −1 [24], collected by CMS in 2016-2018. An integrated luminosity of 1 fb −1 corresponds to ≈100 trillion proton-proton collisions at a centre-of-mass energy of 13 TeV.
In addition, a few projections are made for an assumed data sample corresponding to L = 3000 fb −1 , recorded at √ s = 14 TeV, expected to be accumulated by the end of the next decade during the high-luminosity operation of the LHC accelerator (HL-LHC).

The CMS experiment and data sets
The CMS apparatus [25], illustrated in Extended Data Fig. B.1, is a multipurpose, nearly hermetic detector, designed to trigger on [26,27] and identify electrons (e), muons (µ), photons (γ), and (charged and neutral) hadrons [28][29][30]. A trigger is a filter that selects interesting events, where event refers to the result of the selected interaction in a beam crossing, as observed in the detector. A global event reconstruction algorithm [31] combines the information provided by the all-silicon inner tracker, crystal electromagnetic calorimeter, and brass and scintillator hadron calorimeters, operating inside a 3.8 T superconducting solenoid, with data from gas-ionization muon detectors embedded in the solenoid flux-return yoke, to build electrons, muons, tau (τ) leptons, photons, hadronic jets, missing transverse momentum (p miss T ), and other physics objects [32][33][34]. Collimated streams of particles arising from the fragmentation of quarks or gluons are called jets. These jets are identified, and their energies measured, by specialized reconstruction algorithms [31,33]. The missing transverse momentum vector is measured with respect to the incoming proton beams, and it is computed as the negative vector sum of transverse momenta of all particles in an event.
By July 2012, CMS had collected L = 5.1 fb −1 at a proton-proton (pp) collision centre-of-mass energy √ s = 7 TeV (in 2011) and L = 5.3 fb −1 at 8 TeV (in the first half of 2012), with which the Higgs boson was discovered. By the end of 2012 (Run 1), CMS had collected data corresponding to L = 19.7 fb −1 at 8 TeV [35].
In LHC Run 2 (2015-2018), the accelerator delivered collisions at √ s = 13 TeV. At this larger energy the cross section for Higgs boson production increases by a factor of 2.2-4.0, depending on the production mode [36][37][38][39]. Physics analyses presented here are based on 2016-2018 data, corresponding to an L of up to 138 fb −1 (the additional ∼2 fb −1 recorded in 2015 are not used in this combination). This enabled not only a reduction of statistical, but also of systematic uncertainties, as well as a more precise calibration of the calorimeters and alignment of the tracking detectors. During Run 2, approximately 8 million Higgs bosons have been produced. Many more final states could be studied, since it was possible to separate the events by production mode and decay channel, as well as by kinematic properties; and differential distributions could be measured. Furthermore, improved analysis methods were deployed.
To enable comparison with the more precise experimental results, theoretical calculations have been carried out with commensurate improvements in accuracy [36][37][38][39], involving higher orders in perturbation theory.
The statistical procedure was developed in preparation for the search/discovery of the Higgs boson and hasn't changed much since then. It is based on building a combined likelihood from the various input channels (cf. Methods Section A.5). Parameter estimation and limit setting are performed using a profile likelihood technique with asymptotic approximation [40], taking into account the full correlation of the systematic uncertainties between individual channels and the years of data taking. The different channels included in the combination correlate nuisance parameters related to the same underlying effect, such as the uncertainty in the theoretical prediction or the energy scale uncertainty of the final-state objects. The inclusive signal strength (µ) combination has a total of O(10 4 ) nuisance parameters. The references to the individual analyses presented in the next section each contain more details of the statistical procedure employed for combining the several categories used, created according to various criteria, such as signal-to-background ratios, mass resolutions, and multiplicities of physics objects.

Portrait of the Higgs boson
The portrait of the Higgs boson is defined by its production modes, via cross sections, and its decay channels, via branching fractions. For the value of mass measured by CMS m H = 125.38 ± 0.14 GeV [41], these are given in Extended Data

Production
The rate of production of Higgs bosons is given by the product of the instantaneous luminosity, measured in units of cm −2 s −1 , and the cross section, measured in units of cm 2 . For m H = 125.38 GeV, the total cross section for the production of the SM Higgs boson at √ s = 13 TeV is 55.4 ± 2.6 pb [39]. (A cross section of 1 pb (picobarn) corresponds to an area of 10 −36 cm 2 .) This results in the production of one Higgs boson every second at an instantaneous luminosity of 2 × 10 34 cm −2 s −1 . The dominant production mode in the SM is ggH, where a pair of gluons, one from each of the incident protons, fuses, predominantly via a virtual top quark quantum loop. This is depicted in Fig. 1a and represents 87% of the total cross section. The next most important production mode is vector boson fusion (VBF) depicted in Fig. 1b, where a quark from each of the protons radiates a virtual vector boson (W or Z), which then fuse together to make a Higgs Higgs boson production modes a)

Higgs boson decay channels
Higgs boson pair production k)

Figure 1: Feynman diagrams for the leading Higgs boson interactions
Higgs boson production in (a) gluon-gluon fusion (ggH), (b) vector boson fusion (VBF), (c) associated production with a W or Z (V) boson (VH), (d) associated production with a top or bottom quark pair (ttH or bbH), (e, f) associated production with a single top quark (tH); with Higgs boson decays into (g) heavy vector boson pairs, (h) fermion-antifermion pairs, and (i, j) photon pairs or Zγ; Higgs boson pair production: (k, l) via gluon-gluon fusion, and (m, n, o) via vector boson fusion. The different Higgs boson interactions are labelled with the coupling modifiers κ, and highlighted in different colours for Higgs-fermion interactions (red), Higgsgauge-boson interactions (blue), and multiple Higgs boson interactions (green). The distinction between a particle and its antiparticle is dropped.
boson. Other processes, with smaller cross sections, are: production in association with a vector boson or Higgsstrahlung (VH) depicted in Fig. 1c, and production in association with top (tH, ttH) or bottom (bbH) quarks, depicted in Figs. 1d-f. The bbH mode has not been studied in the context of the SM Higgs boson because of limited sensitivity.
Events are categorized according to the signatures particular to each production mechanism. For example, they are categorized as VBF produced if there are two high transverse momentum (p T ) jets, or as VH produced if there are additional charged leptons ( ) and/or p miss T , or ttH, tH produced if there are jets identified as coming from b quarks, or otherwise ggH produced. (The top quark predominantly decays into a W boson and a b quark jet.)

Decays
In the SM, particle masses arise from spontaneous breaking of the gauge symmetry, through gauge couplings to the Higgs field in the case of vector bosons, and Yukawa couplings in the case of fermions. The SM Higgs boson couples to vector bosons, with an amplitude proportional to the gauge boson mass squared m 2 V , and to fermions with an amplitude proportional to the fermion mass m f . Hence, e.g. the coupling is stronger for the 3rd generation of quarks and leptons than for those in the 2nd generation. The observation of many Higgs boson decays to SM particles and the measurement of their branching fractions are a crucial test of the validity of the theory. Any sizeable deviation from the predictions could indicate the presence of BSM physics.
The Higgs boson, once produced, rapidly decays into a pair of fermions or a pair of bosons. In the SM its lifetime is τ H ≈ 1.6 × 10 −22 s, and its inverse, the natural width, is Γ = /τ H = 4.14 ± 0.02 MeV [39], where is the reduced Planck's constant. The natural width is the sum of all the partial widths, and the ratios of the partial widths to the total width are called branching fractions and represent the probabilities for that decay channel to occur. The Higgs boson does not couple directly to massless particles (e.g. the gluon or the photon), but can do so via quantum loops (e.g. Figs. 1a, 1i, and 1j).
By design, the event selections do not overlap among analyses targeting different final states. Where the final states are similar, the overlap has been checked and found to be negligible.
Detailed information on the analyses included in the new combination along with improvements, and the online and offline criteria used to select events for the analyses can be found in the Methods Section, in Extended Data Tables B.2

Higgs boson pair production
The measurement of the pair production of Higgs bosons can probe its self-interaction λ. The pair production modes are shown in Figs. 1k-o.
In the gluon-gluon fusion mode, there are two leading contributions: in the first (Fig. 1l), two Higgs bosons emerge from a top/bottom quark loop; in the second (Fig. 1k), a single virtual Higgs boson, H * , emerges from the top/bottom quark loop and then decays to two Higgs bosons (gg → H * → HH). Explicit establishment of the latter contribution, a direct manifestation of the Higgs boson's self-interaction, would elucidate the strikingly unusual potential of the BEH field.
In the VBF mode, there are three subprocesses that can lead to production of a pair of Higgs bosons: • via a virtual Higgs boson (Fig. 1m); • via a four-point interaction: VV → HH (Fig. 1n); and • via the exchange of a vector boson (Fig. 1o).
The scattering amplitudes of the processes giving rise to Higgs boson pair production via gluon-gluon fusion (Figs. 1k, 1l) are similar in magnitude, but have opposite signs and interfere destructively. This makes the overall Higgs boson pair production rate small, rendering its experimental observation challenging. The SM Higgs boson pair production cross section is calculated for m H = 125 GeV to be 32.76 +1.95 −6.83 fb [54-56], three orders of magnitude smaller than the single Higgs boson cross section.
The search for Higgs boson pair production is performed by combining Higgs bosons candidates reconstructed from different final states [57-62]. All final states analyzed are defined to be mutually exclusive so that they could be combined as statistically independent observations.

Measurement of the properties of the Higgs boson
At the time of the Higgs boson discovery [2,3], the combination of CMS data gave an observed (obs.) statistical significance of 5 Using the angular distributions of the leptons in the bosonic decay channels, the spin (J) and parity (P, a parity transformation that effectively turns a phenomenon into its mirror image) were also found to be compatible with the SM prediction (J P = 0 + ) with a large number of alternative spin-parity hypotheses ruled out at >99.9% confidence level (CL) [64, 65]. The total cross section, combining all of the different decay channels, was measured to be in agreement with the SM, with an uncertainty of 14% [35]. Each of the VBF, VH, and ttH production modes was measured at a level of 3 s.d. [35].
With Run 2 data, CMS has observed the Higgs boson decaying into a pair of tau leptons with a significance of 5.9 s.d.
[66], a pair of bottom quarks with a significance of 5.6 s.d. [48], and the ttH production mode at 5.2 s.d. [67]. The Higgs boson has also been seen in its decays into muons with a significance of 3 s.d. [52]. The mass of the Higgs boson has been measured to be 125.38 ± 0.14 GeV using the decay channels H → γγ and H → ZZ → 4 [41]. The natural width of the Higgs boson has been extracted and is found to be Γ H = 3.2 +2.4 −1.7 MeV by using offmass-shell and on-mass-shell Higgs boson production [68]. On-mass-shell refers to a particle with its physical mass, while off-mass-shell refers to a virtual particle.

The µ-framework for signal strengths.
The agreement between the observed signal yields and the SM expectations can be quantified by fitting the data with a model that introduces signal strength parameters. These are generically labelled µ, and scale the observed yields with respect to those predicted by the SM, without altering the shape of the distributions. The specific meaning of µ varies depending on the analysis. For given initial (i) and final (f) states, i → H → f , the signal strengths for individual production channels, µ i , and decay modes, µ f , are defined as where σ is the production cross section and B is the branching fraction. Perfect agreement with SM expectations would yield all µ equal to one.
A first test of compatibility is performed by fitting all data from production modes and decay channels with a common signal strength parameter, µ. At the time of discovery, the common µ was found to be 0.87 ± 0.23. The new combination of all Run 2 data yields µ = 1.002 ± 0.057, in excellent agreement with the SM expectation. The uncertainties in the new measurement correspond to an improvement by a factor of 4.5 in precision compared to what was achieved at the time of discovery. At present, the theoretical uncertainties in the signal prediction, the experimental statistical, and the systematic uncertainties separately contribute at a similar level, and they are 0.036, 0.029, and 0.033, respectively.
Relaxing the assumption of a common signal strength parameter, and introducing different µ i    and µ f , our measurements are shown in Fig. 2. The production modes ggH, VBF, WH, ZH, and ttH are all observed with a significance of 5 s.d. or larger.

The κ-framework for coupling modifiers.
Beyond-the-SM physics is expected to affect the production modes and decay channels in a correlated way if they are governed by similar interactions. Any modification in the interaction between the Higgs boson and, e.g. the W bosons and top quarks not only would affect the H → WW (Fig. 1g) or H → γγ ( Fig. 1i and 1j) decay rates, but also the production cross section for the ggH (Fig. 1a), WH (Fig. 1c), and VBF ( Fig. 1b) modes. To probe such deviations from the predictions of the SM, the κ-framework [38] is used. The quantities, such as σ i , Γ f , and Γ H , computed from the corresponding SM predictions, are scaled by κ 2 i , as indicated by the vertex labels in Fig. 1. As an example, for the decay H → γγ proceeding via the loop processes of Fig. 1i or 1j, the branching fraction is proportional to κ 2 γ or (1.26κ W − 0.26κ t ) 2 . In the SM all κ values are equal to one.
A first such fit to Higgs boson couplings introduces two parameters, κ V and κ f , scaling the Higgs boson couplings to massive gauge bosons and to fermions, respectively. With the limited data set available at the time of discovery, such a fit provided first indications for the existence of both kinds of couplings. The sensitivity with present data is much improved, and both coupling modifiers are measured to be in agreement, within an uncertainty of 10%, with the predictions from the SM, as shown in Fig. 3    (left) Constraints on the Higgs boson coupling modifiers to fermions (κ f ) and heavy gauge bosons (κ V ), in different data sets: discovery (red), the full LHC Run 1 (blue), and the data presented here (black). The SM prediction corresponds to κ V = κ f = 1 (diamond marker).
(right) The measured coupling modifiers of the Higgs boson to fermions and heavy gauge bosons, as functions of fermion or gauge boson mass, where υ is the vacuum expectation value of the BEH field (cf. Methods section A.7). For gauge bosons, the square root of the coupling modifier is plotted, to keep a linear proportionality to the mass, as predicted in the SM. The p-value with respect to the SM prediction for the right plot is 37.5%.
A second fit is performed to extract the coupling modifiers κ for the heavy gauge bosons (κ W , κ Z ) and the fermions probed in the present analyses (κ t , κ b , κ τ , κ µ ). Predictions for processes that in the SM occur via loops of intermediate virtual particles, e.g. Higgs boson production via ggH, or Higgs boson decay to a pair of gluons, photons, or Zγ, are computed in terms of the κ i above. The result is shown in Fig. 3 (right), as a function of the mass of the probed particles. The remarkable agreement with the predictions of the BEH mechanism over three orders of magnitude of mass is a powerful test of the validity of the underlying physics. Statistical and systematic uncertainties contribute at the same level to all measurements, except for κ µ , which still is dominated by the statistical uncertainty.
In extensions of the SM with new particles, the loop-induced processes may receive additional contributions. A more general fit for deviations in the Higgs boson couplings can then be defined by introducing additional modifiers for the effective coupling of the Higgs boson to gluons (κ g ), photons (κ γ ), and Zγ (κ Zγ ). Results for this fit are shown in Fig. 4 (left). Coupling modifiers are probed at a level of uncertainty of 10%, except for κ b and κ µ (≈20%), and κ Zγ (≈40%), and all measured values are compatible with the SM expectations, to within 1.5 s.d. These measurements correspond to an increase in precision by a factor of ≈5 compared to what was possible with the discovery data set. Parameter value  Figure 4: Coupling modifiers measurements and their evolution in time.
(left) Coupling modifiers resulting from the fit. The p-value with respect to the SM prediction is 28%. (right) Observed and projected values resulting from the fit in the κ-framework in different data sets: at the time of the Higgs boson discovery, using the full data from LHC Run 1, in the data set used in this paper, and the expected 1 s.d. uncertainty at the HL-LHC for L = 3000 fb −1 . The H → µµ and κ t measurements were not available for earlier data sets due to the lack of sensitivity.

Results from the search for Higgs boson pair production
The cross section for Higgs boson pair production in the SM is extremely small, thus escaping detection at the LHC so far. The results of the search are therefore expressed as an upper limit on the production cross section. Figure 5 (left) shows the expected and observed limits on Higgs boson pair production, expressed as ratios with respect to the SM expectation, in searches  using the different final states and their combination. With the current data set, and combining data from all currently studied modes and channels, the Higgs boson pair production cross section is found to be less than 3.4 times the SM expectation at 95% CL. Figure 5 (right) shows the evolution of the limits from the three most sensitive modes and the overall combination for: the first comprehensive set of measurements using early LHC Run 2 data (35.9 fb −1 ) [73], the present measurements using the full LHC Run 2 data (138 fb −1 ), and the projections for the HL-LHC (3000 fb −1 ) [69]. The HL-LHC projections are also expressed as limits, assuming that there is no Higgs boson pair production. The fact that the combined limit is expected to be below unity shows that the sensitivity is sufficient to establish the existence of the SM HH production. Figure 6 presents the expected and observed experimental limits on the HH production cross section as functions of the Higgs boson self-interaction coupling modifier κ λ and the quartic VVHH coupling modifier κ 2V . Cross section values above the solid black lines are experimentally excluded at 95% CL. The red lines show the predicted cross sections as functions of κ λ or κ 2V , which exhibit a characteristic dip in the vicinity of the SM values (κ = 1) due to the destructive interference of the contributing production amplitudes, as highlighted in Section 4. The experimental limits on the Higgs boson pair production cross section (black lines) also show a strong dependence on the assumed values of κ. This is because the interference between different subprocesses, besides changing the expected cross sections, also changes the differential kinematic properties of the two Higgs bosons, which in turn affects strongly the efficiency for detecting signal events. With the current data set we can ascertain at 95% CL that the Higgs boson self-interaction coupling modifier κ λ is in the range −1.24 to 6.49, while the quartic κ 2V coupling modifier is in the range 0.67 to 1.38. Figure 6 (right) shows that κ 2V = 0 is excluded, with a significance of 6.6 s.d., establishing the existence of the quartic coupling VVHH depicted in Fig. 1n. Combined expected and observed 95% CL upper limits on the HH production cross section for different values of κ λ (left) and κ 2V (right), assuming the SM values for the modifiers of Higgs boson couplings to top quarks and vector bosons. The green and yellow bands represent, respectively, the 1 and 2 s.d. extensions beyond the expected limit; the red solid line (band) shows the theoretical prediction for the HH production cross section (its 1 s.d. uncertainty). The areas to the left and to the right of the hatched regions are excluded at 95% CL.

Current knowledge and future prospects
The discovery of the Higgs boson in 2012 completed the particle content of the standard model (SM) of elementary particle physics, a theory that explains visible matter and its interactions in exquisite detail. The completion of the SM spanned 60 years of theoretical and experimental work. In the ten years following the discovery, great progress has been made in painting a clearer portrait of the Higgs boson.
In this paper, the CMS Collaboration reports the most up-to-date combination of results on the properties of the Higgs boson, based on data corresponding to an L of up to 138 fb −1 , recorded at 13 TeV. Many of its properties have been determined with accuracies better than 10%. All measurements made so far are found to be consistent with the expectations of the SM. In particular, the overall signal strength parameter has been measured to be µ = 1.002 ± 0.057. It has been shown that the Higgs boson directly couples to bottom quarks, tau leptons, and muons, which had not been observed at the time of the discovery, and also proven that it is indeed a scalar particle. The CMS experiment is approaching the sensitivity necessary to probe Higgs boson couplings to charm quarks [74]. The observed (expected) 95% CL value for κ c is found to be 1.1 < |κ c | < 5.5 (|κ c | < 3.40), the most stringent result to date. Moreover, the recent progress in searches for the pair production of Higgs bosons has allowed the setting of tight constraints on the Higgs boson self-interaction strength, and the setting of limits on the Higgs boson pair production cross section not much above twice the expected SM value.
Much evidence points to the fact that the SM is a low-energy approximation of a more comprehensive theory. In connection with the mechanism of spontaneous symmetry breaking, several puzzles appear: the so-called naturalness, a technical issue related to the fact that the Higgs boson mass is close to the electroweak scale; in relation with cosmology, the metastability of the vacuum state of the SM and the conjectured period of inflation in the early universe; the dynamics of electroweak phase transition and its connection to the matter-antimatter asymmetry of our universe. These issues motivate attempts at obtaining a deeper understanding of the physics of the Higgs boson. The impressive progress made over the last decade is foreseen to continue into the next one. The current data set is expected to be doubled in size by the middle of this decade, enabling the establishment of rare decays channels such as H → µµ and H → Zγ. Operation with the high-luminosity LHC is expected during the next decade and should yield ten times more data then originally foreseen. This should allow the ATLAS and CMS experiments to establish the SM Higgs boson pair production with a significance of 4 s.d., as well as the Higgs boson coupling to charm quarks, and to search for any exotic decays. Improvements in experimental techniques and theoretical calculations are also anticipated to continue.
The CMS experiment is entering the era of precision Higgs physics that will shed light on the physics beyond the SM.

Acknowledgments
We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centres and personnel of the Worldwide LHC Computing Grid and other centres for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC, the CMS detector, and the supporting computing infrastructure provided by the following funding agencies:  [23] G. P. Salam, L. T. Wang, and G. Zanderighi, "Higgs turns 10". In this volume [reference to be supplied by the publisher].  [32] CMS Collaboration, "Identification of hadronic tau lepton decays using a deep neural network", 1, 2022. arXiv:2201.08458. Submitted to JINST.   [57] CMS Collaboration, "Search for Higgs boson pair production in the four b quark final state in proton-proton collisions at √ s = 13 TeV", 2, 2022. arXiv:2202.09617. Submitted to PRL.
[68] CMS Collaboration, "First evidence for off-shell production of the Higgs boson and measurement of its width", 2022. arXiv:2202.06923. Submitted to Nature Phys.

A.1 Large Hadron Collider project and the Higgs boson
The primary goals of the LHC and its two general-purpose experiments, ATLAS and CMS, are to: • elucidate the mechanism of EW symmetry breaking and find the associated particle, which in the SM of particle physics is the Higgs boson [4][5][6], and, • search for BSM physics.
The necessity to study the wide range of processes in Fig. 1

A.2.1 Design criteria and the SM Higgs boson
In the early 1990s, during the design phase of the Compact Muon Solenoid (CMS) experiment, considerable emphasis was placed on the identification and measurement of high-energy electrons, photons, and muons, as these particles were expected to play an important role in the search for the SM Higgs boson and in the search for BSM physics.
As the rate of production of energetic muons at high-luminosity hadron colliders is very large, the online selection of events using muons is a particularly formidable task. The muon momentum has to be measured in real time and a momentum threshold placed to limit the rate. This requires a high bending power (high magnetic field) and an adequately precise and robust measurement of the trajectory of muons. This consideration determined the starting point of the design of CMS, and by implication the choice, size, and the power of the analyzing magnet. The next design priority was driven by the search for the Higgs boson via its decay H → γγ, requiring an excellent EM calorimeter (ECAL). The muon system and the ECAL were to be complemented by a precision inner tracking system, immersed in a high magnetic field, giving good momentum resolution and a hadron calorimeter (HCAL) that provided an almost full calorimetric coverage (e.g. for the search for the Higgs boson if its mass turned out to be larger than 500 GeV).

A.2.2 The CMS detector
The longitudinal cut-away view of the CMS detector is shown in Extended Data Fig. B.1. The CMS detector comprises four principal layers: the inner tracker, the ECAL, the HCAL, and the muon system. The various types of detecting elements and their channel counts are also indicated. Physics objects (e.g. electrons, photons, muons, quark or gluon jets, etc.) are identified by different combinations of the patterns of energy deposits and/or traces in these four layers.
The defining choice and the central element of the CMS detector is the long (13 m), large inner diameter (≈6 m), state-of-the-art high-field (3.8 T) superconducting solenoid, generating the magnetic field both for the inner tracker and the muon system. The large size of the solenoid allows the inner tracker and almost all the calorimetry to be installed inside the solenoid.

A.2.2.1 Inner tracking
Particles emerge from the interaction region into the inner tracker, housed in a cylindrical volume with a length of 5.8 m and a diameter of 2.5 m. The particles first encounter the pixel detector, configured in three (four) cylindrical layers of silicon sensors in the barrel region, and two (three) disks in the endcap region before (after) 2017. The pixel detector is surrounded by 10 concentric layers of silicon sensors in the barrel region, with 10 cm long or 20 cm long silicon microstrips, and 12 vertical planes in each endcap region. Points are measured with an accuracy of ≈15 µm in the bending plane. The geometric coverage extends down to angles of 9 • from the beamline.

A.2.2.2 Electromagnetic and hadron calorimeters
The ECAL employs dense lead tungstate scintillating crystals. Each crystal has a length of ≈23 cm that is sufficient to contain the full energy of high-energy electron and photon showers. The amount of generated or collected light is proportional to the energy of the incident particle. The fine transverse size of the crystals means that the energy of an EM shower is distributed over a cluster of crystals ranging from 9 (3×3) to 25 (5×5) crystals. The geometric coverage of the ECAL goes down to about 6 • from the beamline.
The HCAL, comprising ≈7000 channels, is a sandwich of ≈5 cm thick brass absorber plates and ≈4 mm thick scintillator plates. The charged particles in the shower, generated in the absorber plates, traverse the scintillator plates and produce light that is collected and guided by fibres to the photodetectors. The geometric coverage of the HCAL goes down to about 6 • from the beamline. This coverage is augmented by the very forward calorimeter, comprising an iron absorber with quartz fibres embedded in a matrix arrangement. The relativistic charged particles in the showers traverse the fibres and generate Cherenkov light, a part of which is guided by the fibres to the photodetectors. This calorimeter extends the calorimetric coverage down to an angle of ≈0.75 • from the beamline. The thickness of hadron calorimetry is sufficient to absorb almost all of the energy of high energy hadrons.

A.2.2.3 Muon system
Muons (and neutrinos) are the only particles that normally reach the muon system. All other particles deposit almost all of their energies in the calorimeters, and hence are said to have been absorbed. In addition to the measurements inside the inner tracker, the momentum of muons is measured a second time in gas-ionization chambers. These chambers are organized in four stations that measure several points, to a precision of ≈150 µm, and generate track segments whose direction is measured online with an angular precision of ≈5 mrad. An independent set of gas-ionization chambers, provide a signal timing resolution of ≈3 ns aiding the triggering process. The instrumented geometric coverage of the muon system goes down to an angle of 10 • from the beamline.

A.2.2.4 Event selection
As the resources needed to record data for later use from all ≈32 million beam crossings per second would be prohibitively costly, specific filters (known as triggers) are used to select the most interesting ones. An online two-tiered trigger system [26,27] is deployed, with the first tier (Level-1) being hardware-based and the second one (high-level or HLT) being software-based. The Level-1 uses custom hardware that processes coarse information from the calorimeters or the muon chambers to select around 100 thousand crossings of interest per second, corresponding to a reduction of a factor of ≈400. Crossings of interest are selected if the energy deposits in the calorimeters or the momentum of muons, are above predefined thresholds. Upon the issuance of a Level-1 trigger, and after a fixed latency of just under 4 µs, all data from the triggered crossing are off-loaded from the pipeline memories in the ≈100 million on-detector electronics channels. These data, after suitable treatment in electronics housed in the underground services cavern, are sent up 100 m to the surface as fragments on ≈1000 optical fibres and fed into a commercial telecommunication switch. The switch takes the individual fragments, puts them together, i.e. builds the event, and feeds the event into the next available CPU core, in a computer farm of some 50 000 CPU cores. There, in real time, fullevent physics-grade software algorithms, optimized for fast processing, reconstruct physics objects and select for permanent storage some 1000 events or crossings per second, based on topological and kinematic information (see Extended Table B.3).

A.2.2.5 Event reconstruction
The CMS experiment generates a large amount of collision and simulated data. To handle, store, and analyze all these data required the development of the worldwide LHC distributed computing grid (wLCG), providing universal access to data for all CMS Collaboration members.
The data from the stored events are transferred to the Tier-0 centre housed on CERN's main site, where a first processing stage is performed. The result of this stage is then distributed to seven other major centres worldwide, labelled Tier-1 centres, for offline analysis. The Tier-1s are designed to carry out tasks of further reconstruction of the collision data with improved calibration and alignment of the various CMS subdetectors, while the generation and reconstruction of Monte Carlo event samples is carried out both at the Tier-1 centres and smaller university-based locations, labelled Tier-2 centres.
The particle-flow (PF) algorithm [31] reconstructs and identifies each individual particle in an event, with an optimized combination of information from the various elements of the CMS detector. The energy of photons is obtained from the measurements in the ECAL. The energy of electrons is determined from a combination of the electron momentum at the primary interaction vertex as determined by the tracker, and the energy in the corresponding cluster of crystals, including the energy sum of all bremsstrahlung photons spatially compatible with originating from the electron track. The momentum of muons is derived from the curvature of the corresponding track. The energy of charged hadrons is determined from a combination of their momentum measured in the tracker, and the matching ECAL and HCAL energy deposits. The energy of neutral hadrons is obtained from the corresponding corrected ECAL and HCAL energies.
Hadronic jets, arising from quarks or gluons, are created from all the particles reconstructed by the PF algorithm within a cone of half-angle of ≈25 • , centred on the axis determined by the vectorial sum of the momenta of all particles in the jet.

A.2.2.6 Improvements of the CMS detector
Several improvements have been introduced into the CMS experiment since the discovery of the Higgs boson in 2012. These included: • The replacement, in late 2016, of the silicon pixel detector, with a new one comprising four concentric layers in the barrel region, at radii of 29, 68, 109, and 160 mm, and six endcap disks placed at ±34, ±41, and ±51 mm from the interaction point, along the beamline. The new configuration leads to an improvement in the reconstruction of the secondary vertices and in the quality of tagging of b quarks. The sensitivity of H → bb analysis is found to be improved by a factor of 2. • The replacement of photodetectors in HCAL (hybrid photodiodes replaced by sili-con photomultipliers) and implementation of more precise timing, allowing a reduction of accidental or instrumental backgrounds, e.g. stray or out-of-time particles.
• The installation in 2013 and 2014 of chambers in the fourth endcap muon station that were left out for Run 1.
• The upgrade of the Level-1 trigger hardware prior to LHC Run 2 to improve the selection of physics events of interest. The trigger rate from background processes is reduced and the trigger efficiency improved for a wide variety of physics signals. In the muon system, new trigger processor boards deploy powerful commercial field-programmable gate arrays (FPGAs). A time-multiplexed architecture was introduced that enabled data from all the calorimetry in each crossing to be pushed into a single FPGA of the type used in the muon trigger system. The FPGAs allow sophisticated and innovative algorithms to be implemented and evolved as conditions change.
• In the data acquisition system a new switch was installed and the CPU power of the computer farm increased. The whole fabric of the distributed computing systems was upgraded to allow more events to be stored (at least 1000 events/s instead of the initially foreseen 100 events/s).

A.3 Offline event analysis
The principal physics objects are required to have transverse momenta or energies above a set threshold. The thresholds are lowered for the second, or any further, objects. Typical values of these thresholds are listed in Extended Table B.3.
Lepton and photons resulting from the decays of Higgs bosons are expected to be unaccompanied by other particles; they are said to be isolated. Isolation criteria are imposed by requiring no additional energetic particles within a cone of ≈20 • opening angle around the object's direction. Particles, other than from decays of b and c quarks or τ leptons, are expected to emerge directly from the primary interaction vertex, defined as the vertex corresponding to the pp collision identified by the online selection.
Increased use of regression and classification algorithms implemented using machine learning (ML) methods, such as deep neural networks (DNNs) and boosted decision trees (BDTs), led to a simultaneous increase in purity and in efficiencies of identification and reconstruction of physics objects (electrons, muons, photons, b quarks, tau leptons, jets, and p miss T ), and improvements in the calibration of related kinematic observables.
All analyses make extensive use of Monte Carlo simulation of the signal and background processes. The CMS detector is precisely described in software code that is used to generate Monte Carlo event samples. Multiple interactions are included, that match the distribution of the number of pileup interactions observed in data. All the simulated event samples are then processed through the same chain of software programs and procedures as are collision data. Simulated samples are used to evaluate or determine geometric acceptances, energy/momentum/mass resolutions, as well as for online and offline particle identification and reconstruction efficiencies, and for training for the many BDT algorithms and DNNs.

A.4.1.1 H → γγ
The signal is extracted by measuring the narrow signal peak over a smoothly falling background distribution [42]. Despite its small branching fraction (0.23%), this mode is a sensitive one due to the excellent precision in the measurement of the energies of photons. The diphoton invariant mass resolution is σ m γ γ /m H ≈ 1%. All the principal production modes can be studied (ggH, VBF, VH, ttH, and tH). The background largely consists of an irreducible one from quantum chromodynamics (QCD) production of two photons. There is also a reducible background where one or more of the reconstructed photon candidates originate from misidentification of jet fragments, that is dominated by QCD Compton scattering from quarks.

A.4.1.2 H → ZZ → 4
The study of this decay channel uses the distinctive decay of the Z bosons to charged leptons ( ) leading to a final state with 4e, or 4µ, or 2e2µ [43]. The signal appears as a narrow peak on top of a smooth and small background. The momentum (energy) measurement of muons (electrons) is precise enough to give an invariant mass resolution with σ m 4 /m H ≈ 1%. The background comprises an irreducible part arising from the nonresonant production of two Z bosons or Zγ * , and a reducible part from the production of Z+jets and top pair events, where the jets originate from heavy quarks, and thus could contain charged leptons, or are misidentified as charged leptons. The event yield for this process is tiny due to the small branching fractions of H → ZZ (2.71%) and subsequent Z → (3.37% per lepton type) decays. To enhance the signal over background and to categorize events, discriminants exploiting the production and decay kinematics expected for the signal and background events based on a matrix element likelihood approach are used together with the invariant mass of the particle.
Extended Data Fig. B.2 (upper) shows a display of a candidate H → ZZ → eeµµ event produced in pp collisions at a centre-of-mass energy √ s = 13 TeV and recorded in the CMS detector.

A.4.1.3 H → WW → ν ν
Two high-p T and large p miss T characterize this final state [44] and benefit from the H → WW decay having one of the largest branching fractions (∼22%). Due to the presence of two neutrinos, the computation of the WW invariant mass is not possible. However, an associated variable, the transverse mass, m T , can be computed from the p T of the charged leptons and the p miss T . The square of transverse mass for a collection of parti- The dominant background arises from irreducible nonresonant WW production and is estimated from data. The channel has a good sensitivity to the ggH and VBF production processes. In the analysis, 3 and 4 categories are also included, sensitive to production of the Higgs boson in association with a leptonically decaying vector boson. The analysis does not target the ttH and tH production modes, which are covered by a dedicated analysis discussed in Section A.4.3.

A.4.1.4 H → Zγ
The signal is sought as a peak over a smoothly falling background distribution [45]. This analysis targets decays of the Z boson into 2e or 2µ. To increase the sensitivity to the signal, the events are divided into different categories based on the production mode. Multivariate analysis (MVA) techniques are used to further categorize regions with high and low signal-to-background ratios. The dominant background arises from Drell-Yan dilepton production in association with an initial-state photon.

A.4.2.1 H → ττ
Four different ditau final states are studied [46]: eµ, eτ h , µτ h , and τ h τ h , where τ h refers to a hadronically decaying tau lepton. The analysis of this decay channel targets the ggH, VBF, and VH production modes. The identification of τ h candidates uses DNN discriminants to reject quark and gluon jets misidentified as τ h . In order to separate the H → ττ signal events from the sizeable contribution of irreducible Z → ττ events, the likelihood estimate of the reconstructed mass of the ττ system is used. This analysis does not target ttH production, which is covered by the dedicated analysis discussed in Section A.4.3.

A.4.2.2 H → bb
The H → bb decay channel has by far the largest branching fraction of all the decay channels considered, with around 60% of Higgs bosons decaying in this way. The background from QCD production of pairs of b jets is very large, hence final states with special characteristics have been chosen to enhance the signal-to-background ratio [47-51].
To select jets most likely to originate from b quarks, a DNN algorithm is used [75,76]. It provides a continuous discriminant score, which combines information typical of b quark jets, such as the presence of tracks displaced from the primary vertex, identified secondary vertices, and the presence of low-p T leptons in the jet. The threshold on the discriminant score is set such that the misidentification rate for light (u, d, s) quarks or gluons is low. For example, setting this misidentification rate at 0.1% gives a 50% efficiency for b quark jet identification when applied to jets in top quark-antiquark events.
The VH production mode uses the presence of one or more leptons from the decay of the vector boson, or large p miss T . In the signal-sensitive region, DNNs are used to separate the signal from the background dominated by QCD multijet production.
The ttH and tH production modes are included in the combination and MVA techniques are used to separate the signal from the large multijet backgrounds. This analysis uses the 2016 data set.
Lastly, an inclusive analysis is included that targets Higgs bosons produced with large p T [51]. In this kinematic region the signal to background ratio is larger. The two b jets from decays of highly Lorentz-boosted Higgs bosons are close in space and appear in the detector as a single broad jet with distinctive internal structure.
Extended Data Fig. B.2 (lower) shows a display of a candidate H → bb event produced in pp collisions at a centre-of-mass energy √ s = 13 TeV and recorded in the CMS detector.

A.4.2.3 H → µµ
The signal is searched for as a peak in the dimuon mass distribution, over a smoothly falling background [52]. The dimuon invariant mass resolution is σ m µ µ /m H ≈ 1%. The analysis of this decay channel targets the ggH, VBF, VH, and ttH production modes, and is most sensitive in the first two modes. The largest background in this decay channel comes from Drell-Yan dimuon production in which an off-shell Z * boson decays to a pair of muons. Events are split into production modes based on their kinematical properties. To improve the sensitivity of the analysis, MVA techniques are used in each of these different categories.

A.4.2.4 H → cc
The analysis of this final state in the VH production mode (Fig. 1c) has recently been presented [74] but has not been included in the present combination. The analysis yields σ(VH)B(H → cc) < 0.94 pb at 95% CL. The observed 95% CL interval (expected upper limit) for κ c is found to be 1.1 < |κ c | < 5.5 (|κ c | < 3.4), the most stringent to date. A search for Z → cc in VZ events is used to validate the analysis strategy and yields a first observation of this decay channel, at a hadron collider, with a significance of 5.7 s.d.

A.4.3 ttH/tH with multileptons
The ttH (Fig. 1d) and tH (Figs. 1e, 1f) production channels, which probe the coupling of the Higgs boson to the top quarks, are studied in the case where the Higgs boson and the top quarks subsequently decay into final states with several leptons [53], supplementing dedicated studies of the H → γγ, H → ZZ → 4 , and H → bb decay modes.
This analysis uses a categorization based on the number of leptons and/or τ h candidates to target both the different Higgs boson final states and tt decay channels. Categories with at least two leptons, or one lepton and two τ h candidates, target cases where at least one top quark decays via a leptonically decaying W boson. Categories with one lepton and one τ h , or with no leptons and two τ h candidates are used to target events in which both top quarks decay via hadronically decaying W bosons. This analysis is sensitive to the H → WW, H → ττ, and H → ZZ decay channels. Several MVA techniques are employed to better separate the ttH and tH production modes.

A.4.4 Higgs boson decays beyond the SM
In addition to the invisible Higgs boson decays discussed in Section 5.2, other BSM decays are possible, into undetected particles. That is, these particles may or may not leave a trace in the CMS detector, but we do not have dedicated searches looking for these signatures. Nevertheless, the presence of undetected decays can be inferred indirectly from a reduction in the branching fraction for SM decays (or by an increase in the total Higgs boson width). In this interpretation, the total width becomes Γ H = ∑ Γ f (κ)/(1 − B Inv. − B Undet. ), where B Undet. is the branching fraction to undected particles.
To probe invisible or undetected decays of the Higgs boson, another fit can be performed, including B Inv. and B Undet. as additional floating parameters, while imposing as an upper bound on κ W and κ Z their SM values, also valid in most proposed extensions of the SM [77, 78]. As can be seen from Extended Data Fig. B.8 (right), B Inv. and B Undet. are found to be consistent with zero. The 95% CL upper limit on B Undet. is found to be <0. 16, with only small changes to the other κ i fitted values, as shown in Extended Data Fig. B.8 (right). The measurement of the width [68] of the Higgs boson will be used in the future to constrain these quantities without imposing bounds on κ W and κ Z .

A.5 Statistical analysis
The statistical framework used to build the combination of all the channels is based on an established combined likelihood method ([40] and references therein), and briefly detailed in this section.
Given the enormous number of pp collisions produced at the LHC and the relatively small probability that one of those collisions will produce a signal-like event, the observations in data are described by Poisson probability functions, P (k|λ) = e −λ λ k /k!, where k is the observed number of events, and the parameter λ is the expected number of events in a particular bin or region of one or more of the discriminating distributions used for each channel entering the combination. The combined likelihood is obtained from the product of the individual Poisson probability functions, accounting for the observed data and expected number of events across all channels.
The parameters λ are functions of the model parameters of interest: µ, which represent the Higgs boson couplings or signal strengths, and nuisance parameters θ, which model the effect of systematic uncertainties on the predicted signal and background contributions. Additional terms are included in the combined likelihood to represent constraints on the nuisance parameters due to external measurements, such as energy and momentum scale calibrations or an integrated luminosity determination. The measurements reported in this paper are determined using the profile likelihood ratio q(µ) = −2 ln L(µ,θ µ )/L(μ,θ) whereμ andθ are the values of the parameters of interest and nuisance parameters that maximize the likelihood L(µ, θ), and θ µ are the values of the nuisance parameters that maximize the likelihood for a fixed value of µ. The compatibility between a given set of measurements and their corresponding SM predictions is reported as a p-value, derived from the difference between q SM and q(μ). Expected intervals are derived from the Asimov data set, in which the nuisance parameters are set to their maximum likelihood estimator values.
All the reported confidence intervals, confidence regions, and p-values are obtained assuming various asymptotic approximations for the distributions of the (modified) likelihood ratio test statistic [79]. The validity of the asymptotic assumptions has been routinely checked in the context of individual analyses whenever the event yields are small or particular validity conditions are not met.

A.7 Notes on self-interaction strength
The potential energy of the BEH field (φ) is given by The first term accounts for the mass of the Higgs boson m H . The second term represents the Higgs boson self interaction, of strength λ. In the SM, λ = m 2 H /(2υ 2 ) (where the vacuum expectation value of the BEH field, corresponding to its minimum, is υ = 246 GeV) and it can be measured via the study of Higgs boson pair production. The third term represents the interaction of four Higgs bosons at a point, a process that is even rarer than its pair production. Knowledge of the exact shape of the potential V is crucial for understanding the phase transition that occurred in the early universe and its consequences [80].
The search for Higgs boson pair production is performed by combining Higgs boson pairs, each with differing decay modes. The decay modes that have been used are bb, ττ, and WW [57-60], benefitting from the large branching fractions, and γγ [61] and ZZ → 4 [62], benefitting from the presence of narrow mass peaks, thus improving the signal-to-background ratio. All final states analyzed are defined to be mutually exclusive so that they could be properly combined as statistically independent observations. Measurements of Higgs boson pair production are used to constrain the Higgs boson selfinteraction strength λ. Several combinations of individual Higgs boson decay modes are employed in this search. The highest rate for Higgs boson pair production and decays occurs when both Higgs bosons decay to b quark pairs, HH → bbbb, corresponding to ≈35% of all the possible HH decays in the SM.
The search in the 4b decay mode [57, 58] is performed separately under the assumptions that m H * 2m H or not. In the case m H * 2m H , each Higgs boson is energetic (and hence said to be boosted), such that its decay products, e.g. b quark jets, merge and appear as one broad jet, but with a distinctive internal structure. In the latter case, all four b quark jets rarely overlap, and hence said to be resolved. Measurements of single Higgs boson production and decay can also be used to constrain κ λ as quantum corrections to the SM Higgs boson production modes and decay channels depend on κ λ [81, 82]. These corrections have been derived [83] for the different production and decay modes entering the combination, as shown in Extended Data Table B

A.8 Upgrade of the CMS experiment for HL-LHC running
To exploit the full potential of the LHC, the accelerator and its experiments will be upgraded. The HL-LHC will operate at an instantaneous luminosity of 5 × 10 34 cm −2 s −1 . The intention is to collect ten times more data than the 300 fb −1 foreseen in the initial LHC phase. This means that the integrated radiation levels will be correspondingly larger.
The physics to be studied drives the technical choices for the upgrade. The physics goals are: • precise measurements of the properties of the Higgs boson and its self-coupling, to elucidate further the physics of EW symmetry breaking; • search for BSM physics; and • selected precision SM measurements.
The translation of these physics goals into experimental design goals requires: • The construction of a new higher-granularity, more radiation-hard silicon tracker. The design of the new front-end electronics will allow information from the inner tracker to participate in the Level-1 trigger. The size of the individual detecting el-ements will be decreased leading to about ten times larger number of electronics channels. All components inside the tracker (silicon sensors, front-end electronics, 10 Gb/s data links, etc.) will have to withstand integrated doses of up to 500 Mrad and fluences of 10 16 (1 MeV equivalent neutrons)/ cm 2 . The geometric coverage of the inner tracker will be increased, extending it down an angle of 2 • from the beamline.
• The replacement of other components affected by radiation. Principally these are the endcap calorimeters and the ECAL front-end electronics. The endcap calorimeters will be replaced with a new high-granularity "imaging" calorimeter with precision timing. It will be based on 600 m 2 of silicon sensors with detecting cells of sizes of 0.5 to 1.0 cm 2 . Regions in this calorimeter will reach integrated doses of up to 500 Mrad and fluences of 10 16 (1 MeV equivalent neutrons)/cm 2 . The new front-end electronics for the ECAL barrel will allow data from each crystal to be sent to the calorimeter Level-1 trigger processor, instead of the sum of 25 crystals today, and which will allow better measurement of the timing of the impact of electrons or photons.
• Higher-bandwidth Level-1 and high-level triggers. Information from the inner trackers will be used at Level-1. The Level-1 trigger latency will be increased from 4 to over 12 µs, requiring corresponding changes in the front-end electronics, allowing more processing time leading to a purer selection of events. The output rate from the Level-1 processors will be increased from 100 to 750 kHz and correspondingly the number of events stored for later analysis will be increased from 1 to 10 kHz.
• The introduction of precision timing detectors. A new set of detectors will be installed in the barrel and endcap regions, covering a region down to an angle of 9 • from the beamline. The precision timing of photons (in the barrel region) and charged tracks will greatly improve the localization of the correct interaction vertex. At HL-LHC, on the average, some 140 pairs of protons are expected to interact in each crossing, spread over a time characterized by σ ≈ 200 ps. Furthermore, suppression of energy can be carried out that is not consistent in time with the interaction of interest.
The upgraded CMS experiment at HL-LHC will be more powerful than the current one. Uncertainties in many measurements of the properties of the Higgs boson are expected to approach the percent level, benefitting from the anticipated larger event samples, reduced experimental systematic uncertainties, and more accurate theoretical calculations.

A.9 Theoretical references
The theoretical works used in our analyses can be found in the LHC

A.10 Data availability
Tabulated results are provided in the HEPData record for this analysis. Release and preservation of data used by the CMS collaboration as the basis for publications is guided by the CMS data preservation, re-use and open acess policy.

A.12 Ethics declarations
Author contributions All authors have contributed to the publication, being variously involved in the design and the construction of the detectors, in writing software, calibrating subsystems, operating the detectors and acquiring data, and finally analysing the processed data. The CMS Collaboration members discussed and approved the scientific results. The manuscript was prepared by a subgroup of authors appointed by the collaboration and subject to an internal collaboration-wide review process. All authors reviewed and approved the final version of the manuscript.      (upper) An event display of a candidate H → ZZ → eeµµ. (lower) An event display of an H → bb candidate produced in association with a Z boson decaying into an electron-positron pair, in pp collisions at √ s = 13 TeV recorded by CMS. The charged-particle tracks, as reconstructed in the inner tracker, are shown in yellow; the electrons are shown in green, the energy deposited by the electrons in the ECAL is shown as large green towers, the size of which is proportional of the amount of energy deposited; the blue towers are indicative of the energy deposits in the HCAL, while the red boxes are the muon chambers crossed by the muons (red tracks); the yellow cones represent the reconstructed jets. (lower, inset) The zoom into the collision region shows the displaced secondary vertices (in red) of the two b quarks decaying away from the primary vertex (in yellow). One of the bottom hadrons decays into a charm hadron that moves away from the secondary vertex before decaying (b → c → X; vertex in cyan).      Signal strength parameters per individual production mode and decay channel µ f i , and combined per production mode µ i and decay channel µ f . In this fit, ttH and tH are considered together and the µ i results are slightly different from those of Fig. 2 (left). The dashed vertical lines at 1 represent the SM value. Light grey shading indicates that µ is contained to be positive. Dark grey shading indicates the absence of measurement. The p-value with respect to the SM prediction is 5.8%.   Figure B.7: Time evolution of the signal strength measurements and their precision. Comparison of the signal strength parameter µ fit results in different data sets; in each panel, from left to right: at the time of the Higgs boson discovery, using the full data from LHC Run 1, in the data set analyzed for this paper, and the expected 1 s.d. uncertainty for HL-LHC for L = 3000 fb −1 . The H → µµ measurements were not available for the earlier data sets due to the lack of sensitivity.