Combined measurements of Higgs boson couplings in proton-proton collisions at s = 13 Te .

Combined measurements of the production and decay rates of the Higgs boson, as well as its couplings to vector bosons and fermions, are presented. The analysis uses the LHC proton-proton collision data set recorded with the CMS detector in 2016 at s = 13 Te , corresponding to an integrated luminosity of 35.9 fb - 1 . The combination is based on analyses targeting the five main Higgs boson production mechanisms (gluon fusion, vector boson fusion, and associated production with a W or Z boson, or a top quark-antiquark pair) and the following decay modes: H → γ γ , Z Z , W W , τ τ , b b , and μ μ . Searches for invisible Higgs boson decays are also considered. The best-fit ratio of the signal yield to the standard model expectation is measured to be μ = 1.17 ± 0.10 , assuming a Higgs boson mass of 125.09 Ge . Additional results are given for various assumptions on the scaling behavior of the production and decay modes, including generic parametrizations based on ratios of cross sections and branching fractions or couplings. The results are compatible with the standard model predictions in all parametrizations considered. In addition, constraints are placed on various two Higgs doublet models.

Introduction

Understanding the mechanism behind electroweak symmetry breaking (EWSB) remains one of the main objectives of the physics program at the CERN LHC. In the standard model (SM) of particle physics [1-4], EWSB is realized through the addition of a complex scalar doublet field. A salient feature of this is the prediction of one physical, neutral, scalar particle, the Higgs boson () [5-10]. The Higgs scalar field can also account for the fermion masses through Yukawa interactions [2, 11]. The Higgs boson was discovered by the ATLAS and CMS Collaborations [12-14], and is the subject of much study. The Yukawa coupling strengths are free parameters in the SM and do not explain the observed pattern of fermion masses. Furthermore, it is not understood why the Higgs boson mass is near the electroweak scale, since it is not protected in the SM from large quantum corrections [15-19]. This has led to the development of many beyond the SM (BSM) theories that can alter the properties of the Higgs boson [20-24]. Precision measurements of the properties of the Higgs boson are therefore an important test of the SM.

This paper describes combined measurements of the Higgs boson production rates, decay rates, and couplings using analyses of proton–proton collision data recorded with the CMS detector in 2016. The data set corresponds to an integrated luminosity of 35.9. The following decay channels are included in the combination: , , , , , and , as shown in Fig. 1. Here and in what follows, we do not distinguish between particles and antiparticles in our notations of production and decay processes. Searches for invisible decays of the Higgs boson, which are predicted to be considerably enhanced by several BSM theories [25-28], are also considered for selected measurements. The data samples considered for each decay channel are ensured to have negligible overlap to avoid introducing nontrivial correlations.

Fig. 1

Examples of leading-order Feynman diagrams for Higgs boson decays in the , , and (upper left); and (upper right); and (lower) channels

The analyses included in this combination target production via gluon fusion (), vector boson fusion (), associated production with a vector boson (, V or ), and associated production with a pair of top quarks (). The prediction for production has advanced to next-to-next-to-next-to-leading order (NLO) in perturbative quantum chromodynamics (QCD) [29, 30] and next-to-leading order (NLO) for electroweak (EW) corrections, reducing its uncertainty from (next-to-NLO) to . The calculations of the and cross sections are performed at next-to-NLO QCD and NLO EW accuracy, while the calculation of the cross section is performed at NLO QCD and NLO EW accuracy. The updated theoretical predictions used for the various production and decay modes in this paper can be found in Refs. [29-52] and are summarized in Ref. [53]. Examples of leading-order (LO) Feynman diagrams for these production processes can be seen in Figs. 2 and 3. In addition to the five main production processes, the contributions due to Higgs boson production in association with a single top quark () and either a boson () or a quark (), as shown in Fig. 4, are included in the analyses that have some sensitivity to them.

Fig. 2

Examples of leading-order Feynman diagrams for the (upper left), (upper right), (lower left), and (lower right) production modes

Fig. 3

Examples of leading-order Feynman diagrams for the production mode

Fig. 4

Examples of leading-order Feynman diagrams for production via the (upper left and right) and (lower) modes

For certain measurements in this paper, such as production and decay, the interference between the diagrams that contribute to the process is considered. In addition, the cross section is small in the SM, being approximately 14% of the  cross section, due to the destructive interference between the diagrams shown in Fig. 4, which involve the coupling of the Higgs boson to bosons ( process) and top quarks ( process). This interference becomes constructive, however, when the relative sign between these couplings is negative, and so the  process is sensitive to the relative sign of the and couplings.

The ATLAS and CMS Collaborations have published combined measurements of Higgs boson production rates, decay rates, and couplings with the and LHC Run 1 data [54, 55]. A combination of the Run 1 ATLAS and CMS analyses has also been performed [56]. All results were found to be in agreement, within their uncertainties, with the predictions of the SM. In this paper, due to the larger integrated luminosity and increased signal cross section at , the measured precision for several parameters of interest has significantly increased with respect to Ref. [56]. In particular, the predicted cross sections for the dominant production mode and the production mode increase by factors of approximately 2.3 and 3.8, respectively, between and . In addition, some of the theoretical predictions have improved, as mentioned earlier.

This paper is organized as follows: A brief description of the CMS detector is given in Sects. 2,  3 provides a summary of the various analyses included in the combination, and Sect. 4 describes the modifications made to these analyses to ensure a common signal and uncertainty model. Section 5 outlines the statistical procedure used to derive the results, and Sect. 6 outlines the treatment of the systematic uncertainties. Section 7 reports the results of the signal parametrizations in terms of signal strength modifiers and fiducial cross sections, while Sect. 8 describes the results obtained from an alternative set of signal parametrizations in terms of Higgs boson couplings. Section 9 details interpretations in terms of various two Higgs doublet models. The paper is summarized in Sect. 10.

The CMS detector

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter, each composed of a barrel and two endcap sections. Forward calorimeters extend the pseudorapidity coverage provided by the barrel and endcap detectors. Muons are detected in gas-ionization chambers embedded in the steel flux-return yoke outside the solenoid. A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [57].

Analyses included in the combination

In this section, the individual analyses included in the combination are briefly described. More detailed information on each analysis can be found in the corresponding references. Many of the analyses split their primary data sample in multiple event categories with specific signatures that enhance the discrimination power between different Higgs boson production processes. This is achieved through selections that require the presence of additional leptons or jets, as expected in the decay of a or boson in the and modes, or in top quark decays in the mode, and that exploit the distinctive kinematic properties of the final state objects, such as the presence of two jets with a large separation in pseudorapidity , and a large invariant mass , in the topology. In some categories, the kinematic features of an event as a whole are used to select particular production processes. For example, requiring a large missing transverse momentum , defined as the magnitude of the negative vector sum over the transverse momenta of all particles reconstructed in an event, targets production in which the boson decays to neutrinos. The event categories within and amongst the individual analyses are constructed to ensure a negligible level of overlap (i.e. the same event entering more than one category). In many cases, this is accomplished by synchronizing the object (e.g. electron, muon, tau, or jet) identification definitions and imposing strict requirements on the number of reconstructed objects. In other cases, the orthogonality is ensured by imposing opposing requirements on higher level observables formed using multiple objects. For rare cases where potential overlap is not explicitly removed, the lists of selected data events were checked and found to contain a negligible number of duplications. In total, up to 265 event categories are considered, and there are over 5500 nuisance parameters corresponding to various sources of experimental and theoretical systematic uncertainty. A summary of the production and decay modes, which are described in more detail in the following sections, is shown in Table 1.

Table 1

Summary of the event categories in the analyses included in this combination. The first column indicates the decay channel and the second column indicates the production mechanism targeted by an analysis. The third column provides the total number of categories per production tag, excluding control regions. Notes on the expected fractions of different Higgs signal production and decay modes with respect to the total signal yield in the given category are given in the fourth column. Where the numbers do not sum to 100%, the remaining contributions are from other signal production and decay processes. Finally, where relevant, the fifth column specifies the approximate expected relative mass resolution for the SM Higgs boson

Decay tags	Production tags	Number of categories	Expected signal fractions	Mass resolution
H→γγ, Sect. 3.1
γγ	Untagged	4	74–91% ggH	≈1–2%
	VBF	3	51–80% VBF
	VH hadronic	1	25% WH, 15% ZH
	WH leptonic	2	64–83% WH
	ZH leptonic	1	98% ZH
	VH pTmiss	1	59% VH
	ttH	2	80–89% ttH, ≈8% tH
H→ZZ^(∗)→4ℓ, Sect. 3.2
4μ, 2e2μ/2μ2e, 4e	Untagged	3	≈95% ggH	≈1–2%
	VBF 1, 2-jet	6	≈11–47% VBF
	VH hadronic	3	≈13% WH, ≈10% ZH
	VH leptonic	3	≈46% WH
	VH pTmiss	3	≈56% ZH
	ttH	3	≈71% ttH
H→WW^(∗)→ℓνℓν, Sect. 3.3
eμ/μe	ggH 0, 1, 2-jet	17	≈55–92% ggH, up to ≈15% H→ττ	≈20%
	VBF 2-jet	2	≈47% VBF, up to ≈25% H→ττ
ee+μμ	ggH 0, 1-jet	6	≈84–94% ggH
eμ+jj	VH 2-jet	1	22% VH, 21% H→ττ
3ℓ	WH leptonic	2	≈80% WH, up to 19% H→ττ
4ℓ	ZH leptonic	2	85–90% ZH, up to 14% H→ττ
H→ττ, Sect. 3.4
eμ,eτ_h,μτ_h,τ_hτ_h	0-jet	4	≈70–98% ggH, 29% H→WW in eμ	≈10–20%
	VBF	4	≈35–60% VBF, 42% H→WW in eμ
	Boosted	4	≈48–83% ggH, 43% H→WW in eμ
VH production with H→bb, Sect. 3.5
Z(νν)H(bb)	ZH leptonic	1	≈100% VH, 85% ZH	≈10%
W(ℓν)H(bb)	WH leptonic	2	≈100% VH, ≈97% WH
Z(ℓℓ)H(bb)	Low-p_T(V) ZH leptonic	2	≈100% ZH, of which ≈20% ggZH
	High-p_T(V) ZH leptonic	2	≈100% ZH, of which ≈36% ggZH
Boosted H Production with H→bb, Sect. 3.6
bb	p_T(H) bins	6	≈72–79% ggH	≈10%
ttH production with H→leptons, Sect. 3.7.1
2ℓss	ttH	10	WW/ττ≈4.5, ≈5% tH
3ℓ		4	WW:ττ:ZZ≈15:4:1, ≈5% tH
4ℓ		1	WW:ττ:ZZ≈6:1:1, ≈3% tH
1ℓ +2τ_h		1	96% ttH with H→ττ, ≈6% tH
2ℓss+1τ_h		2	ττ:WW≈5:4, ≈5% tH
3ℓ +1τ_h		1	ττ:WW:ZZ≈11:7:1, ≈3% tH
ttH production with H→bb, Sect. 3.7.2
bb	tt¯→ jets	6	≈83–97% ttH with H→bb
	tt¯→1ℓ+jets	18	≈65–95% ttH with H→bb, up to 20% H→WW
	tt¯→2ℓ+jets	3	≈84–96% ttH with H→bb
Search for H→μμ, Sect. 3.8
μμ	S/B bins	15	56–96% ggH, 1–42% VBF	≈1–2%
Search for invisible H decays, Sect. 3.9
Invisible	VBF	1	52% VBF, 48% ggH
	ggH + ≥1 jet	1	80% ggH, 9% VBF
	VH hadronic	1	54% VH, 39% ggH
	ZH leptonic	1	≈100% ZH, of which 21% ggZH

The analysis [58] provides good sensitivity to nearly all Higgs boson production processes. Since the decay proceeds mainly through - and top-loop processes, interference effects make its branching fraction sensitive to the relative sign of the fermion and vector boson couplings. The analysis measures a narrow signal peak in the diphoton invariant mass () spectrum over a smoothly falling continuum background, originating mainly from prompt, nonresonant diphoton production, or from events where at least one jet is misidentified as an isolated photon.

Exclusive event categories are defined using dedicated selections based on additional reconstructed objects to separate the different Higgs boson production mechanisms. The presence of additional leptons, , or jets is used to classify events into one of the following categories: leptonic, hadronic, leptonic, leptonic, loose leptonic with low requirement, , , and hadronic. The category is divided into three subcategories of increasing purity against production. Finally, the remaining events are divided into four untagged categories with increasing signal purity.

In each event class, the background in the signal region (SR) is estimated from a fit to the observed distribution in data. The dominant experimental uncertainties in the measurement of the rate of Higgs boson production in the decay channel are related to the modeling of the electromagnetic shower shape observables used in the photon identification and the background shape parametrization.

Despite the ( or ) decay having the lowest branching fraction of the decay channels considered, it also has the lowest background contamination, resulting in very good sensitivity to production processes with large cross sections, such as . It is also the most important decay channel in constraining the coupling. The  [59] analysis measures a narrow four-lepton invariant mass peak over a small continuum background. The dominant irreducible background in this analysis is due to nonresonant production with both bosons decaying to a pair of charged leptons, and is estimated from simulation. The , , and / decay channels are treated separately to better model the different mass resolutions and background rates arising from jets misidentified as leptons.

To separate the different Higgs boson production mechanisms, the following categories are defined on the basis of the presence of jets, -tagged jets, leptons, , and various matrix element discriminants that make use of the information about the additional objects: (1- and 2-jet), hadronic, leptonic, , , and untagged categories.

In the analysis, the dominant experimental uncertainties are related to the lepton efficiencies and the determination of the +jets background from data.

The analysis [60] profits from the fact that the decay mode has one of the largest branching fractions and has a relatively low-background final state. As a result, this decay channel has very good sensitivity to most production processes, in particular and . Imposing tight lepton identification criteria and requiring the absence of -tagged jets helps to reduce the misidentified lepton and top quark backgrounds, respectively. Several event categories with varying signal-to-background ratios are defined to improve the sensitivity to the signal. Events are selected that contain two leptons, denoted , which may be of different or same flavor. The different-flavor decay channel dominates the sensitivity since it has the largest branching fraction and is the least contaminated by backgrounds. The same-flavor and final states are also considered, although their sensitivity is limited by the contamination from Drell–Yan (DY) background events with misreconstructed . Given the large background contribution from production in both the different-flavor and same-flavor final states, events are further categorized into categories with 0, 1, and 2 associated jets, with the 0-jet category dominating the overall sensitivity. In addition, events are further categorized on the basis of the of the subleading lepton, since the background from misidentified leptons is larger in the low- region. In the different-flavor final state, dedicated 2-jet categories are included to enhance the sensitivity to VBF and production mechanisms.

The analysis also includes categories that are sensitive to the associated production of the Higgs boson with a vector boson that decays leptonically. Two categories that are sensitive to production are defined by requiring the presence of a total of three leptons (electrons or muons). The two are distinguished by whether or not they contain a pair of leptons with the same flavor and opposite sign. Events with four charged leptons, in which one pair is consistent with a boson decay, are separated into two categories depending on whether the remaining pair consists of same-flavor leptons or not. These categories are sensitive to the production mode. The signal extraction method depends on the event category.

When measuring the rate of Higgs boson production in the decay channel, the dominant experimental uncertainties arise from the determination of the top quark pair, and DY backgrounds from data, and the uncertainties related to the and dependent lepton reconstruction and identification efficiencies.

The analysis [61] benefits from a relatively large branching fraction and a reasonable mass resolution of 10–20%, providing competitive sensitivity to both the and production processes. It also provides the best sensitivity for the direct measurement of a fermionic Higgs boson coupling. The analysis utilizes the four most sensitive final states: , , , and , where denotes a hadronically decaying lepton. In the analysis of each decay channel, events are divided into three categories labeled 0-jet, boosted, and VBF.

The VBF category requires the presence of two additional jets with large and , designed to increase the purity of events. The 0-jet category does not have much sensitivity to the signal, but is useful to constrain systematic uncertainties in the background model. The boosted category contains all remaining events, and is binned as a function of of the system to increase the sensitivity to production. There is a nonnegligible contribution from the process in some categories, and this is treated consistently as an signal in this combined measurement.

The and energy scale uncertainties are the dominant experimental uncertainties in the measurement of the Higgs boson production rate in the decay channel, followed by the uncertainties in the determination of the jets background from data.

production with

The decay has the largest expected branching fraction in the SM (58.1% for ) and a reasonable mass resolution of 15%. By requiring production it is possible to increase the signal purity with respect to the inclusive case for which the background from QCD multijet production is dominant. The analysis of the decay targeting production () [62] provides the best sensitivity to the and processes as well as to the coupling. Selected events are categorized based on the presence of two -tagged jets, and two (), one () or no () electrons or muons in the final state. The categories are subdivided into low-boost () and high boost () regions. Events selected in the category are further required to have .

The main backgrounds come from or boson production in association with light- and heavy-flavor (LF and HF) jets, as well as from top quark pair and diboson production. The dominant experimental uncertainties in in this analysis are related to the determination of these backgrounds, and uncertainties in the tagging discriminator shapes and efficiencies.

Boosted production with

The decay is also measured in an analysis that targets inclusive production of the Higgs boson [63], exploiting the higher signal to background ratio at high (the transverse momentum of the Higgs boson). The decay products of a high- system are reconstructed using the anti- algorithm [64, 65] with a distance parameter of 0.8 (AK8 jet), and the soft-drop algorithm [66, 67] is used to reconstruct the jet mass , which peaks at the Higgs boson mass for signal events. Events containing substantial , or identified and isolated electrons, muons or leptons are vetoed to reduce the background contributions from vector boson production and top quark processes.

The main background component, QCD multijet production, is estimated from a signal-depleted Control Region (CR). The selected events are divided according to the jet into six bins of increasing width from 450 GeV to 1.

The dominant experimental uncertainties in this analysis are the uncertainties related to the extrapolation of the QCD multijet and top quark pair backgrounds from the CRs.

production

Measurements of the rate of the production process provide a direct test of the Higgs boson’s coupling to top quarks. A recent measurement by CMS combining the , 8 and datasets was able to establish the first observation of the production process [68]. Dedicated analyses targeting the  [69] and  [70, 71] decay channels using data are described in this section.

production with

The analysis of production with  [69] is mainly sensitive to the Higgs boson decaying to , or with electrons, muons and/or in the final state. This analysis provides the best sensitivity to the production process. The main irreducible backgrounds come from V and diboson production. Reducible backgrounds containing misidentified leptons or leptons with misidentified charge are estimated from CRs in data. Events are categorized according to their lepton content. The light-lepton (/) categories are defined as:The categories, which require the presence of hadronically decaying leptons, are defined as:In the / and categories, the dominant experimental uncertainties on the measurement of the rate of Higgs boson production in the mode are related to the lepton reconstruction efficiencies, and the estimation of the reducible background contributions from data.

2ss: Events with two leptons having the same sign and at least four additional jets. A veto on the presence of hadronic tau decays is applied. Further categories based on lepton charge, flavor and the number of -tagged jets are defined within this class.

3: Events containing three leptons, with the sum of lepton charges equal to ±1, and at least two additional jets of which one or two are tagged.

4: Events with four leptons, with an explicit veto on events as these are selected by the analysis described in Sect. 3.2.

1+2: Events with two oppositely charged candidates and an additional /, at least three additional jets, and at least one -tagged jet.

2ss+1: Events containing three leptons, with sum of lepton charges equal to ±1, and at least two additional jets of which one or two are tagged. These events are further sorted into two subcategories based on whether or not all of the jets expected in the process are reconstructed.

3+1: Events with three light leptons, one and at least two additional jets and one -tagged jet.

There are two analyses that target the associated production of the Higgs boson with a pair of top quarks in the decay mode [70, 71]. The leptonic analysis requires at least one lepton to be present in the final state, from the decay system, while the hadronic analysis selects events in the all-hadronic final state. These analyses provide good sensitivity to the production process and improve the precision in the measurement of the coupling.

In the leptonic analysis, events are sorted into the or classes, depending on the presence of one or two well-identified leptons. Events are further categorized based on the number of reconstructed jets (Nj) and the number of jets that are tagged as jets (Nb) in each event. The largest background is due to top quark pair production with additional jets that contain heavy flavor hadrons. In the class, three categories are used: 4j 3, 5j 3, and 6j 3. In each category a multi-classification deep neural network (DNN) [72] is used to define six classes on the basis of the most probable event hypothesis for each event, yielding a total of 18 categories. In the class, there are two jet categories: 4j 3 and 4j 4. The 4j 4 category is further divided into two subcategories.

The all-hadronic final-state analysis selects events that contain at least seven jets, at least three of which are tagged as jets. These events are divided into seven categories: 7j 3, 7j 4, 8j 3, 8j 4, 9j 3, and 9j 4. Events containing electrons or muons are vetoed to maintain an orthogonal selection to the leptonic final state analysis. The dominant background is QCD multijet production, with other important backgrounds coming from +jets processes.

The dominant experimental uncertainties in the measurement of the rate of production with decay in the leptonic and all-hadronic final states are due to uncertainties in the determination of the backgrounds and tagging efficiencies. In the all-hadronic final state, the uncertainty in the determination of the QCD multijet background also has a significant contribution to the overall systematic uncertainty.

Search for

The search [73] is the only analysis included here that is sensitive to the coupling of the Higgs boson to second-generation fermions. The analysis searches for a narrow peak in the dimuon invariant mass () spectrum above a large continuum background from DY production of muon pairs. Events are categorized using variables that are uncorrelated with , in order to avoid introducing an irregular shape in the background spectrum. Variables that distinguish between the and signals, and the DY and backgrounds, are used to define event categories with varying signal-to-background ratios. The categories are further divided based on the momentum of the muon with the largest in the dimuon pair, to exploit the differences in the resolution. Since there are more variables associated with production that can be used to separate signal and background, the events with the highest BDT output value are most compatible with that process.

In each event category, the background is estimated from a fit to the observed distribution. As in the analysis, the parameters of the functions used to describe the background contribute to the statistical uncertainty in the measurements. This is the dominant source of uncertainty in constraining the rate of Higgs boson decay in the decay channel. The observed upper limit on the cross section times branching fraction of obtained in Ref. [73] is 2.93 times the SM value.

Search for invisible

The direct search for the Higgs boson decaying into particles that cannot be detected provides a constraint on the invisible Higgs boson branching fraction (), which is predicted to be enhanced in BSM scenarios [25–28, 74]. The search is performed using events with large , and containing additional particles consistent with Higgs boson production via the  [75], with  [76], with the or boson decaying hadronically, or modes [77].

Events selected in the category are required to contain two jets, with a large and a large . The hadronic and categories comprise events containing either a high- AK8 jet, consistent with a boosted, hadronically decaying vector boson, or a jet from initial-state radiation, reconstructed in the fiducial volume of the tracker. The dominant backgrounds in these categories are due to the + and + processes. These are estimated from dedicated lepton and photon CRs in data. In all three categories, the dominant uncertainties are related to the extrapolation of the lepton and photon CRs to determine the + and + backgrounds in the SR.

The leptonic category is defined by selecting events that contain a pair of oppositely charged electrons or muons consistent with the decay of a boson. The dominant backgrounds arise from and diboson production and are estimated using a combination of CRs in data containing additional leptons, and simulated events. The dominant uncertainty in this category is related to theoretical uncertainties in the higher-order corrections used in the simulation of these backgrounds.

The observed upper limit on the branching fraction of assuming SM Higgs production rates is 26% [75]. As described later in Sect. 8, the analyses are included only in models for which a nonzero invisible branching fraction of the Higgs boson is considered.

Modifications to the input analyses

This section describes the changes in each analysis, as implemented for this combination, compared to their respective publications.

Gluon fusion modeling

In order to consistently combine the various analyses, it is necessary to use the same theoretical predictions for the signal. The most significant difference between the input analyses is the modeling of the dominant production mode in the , , , and decay channels. The published results in these analyses used different generators with next-to-leading order matrix elements merged with parton showering (NLO+PS). In order to synchronize these analyses and take advantage of the most accurate simulation of available, a reweighting is applied. Gluon fusion events are generated using the powheg 2.0 [78-81], MadGraph 5_amc@nlo version 2.2.2 [82, 83], and nnlops [84, 85] generators. The nnlops simulation, which is the highest order parton shower matched simulation available, includes the effects of finite quark masses. Events are separated into 0, 1, 2, and 3 jet bins, where the jets used for counting are clustered from all stable particles, excluding the decay products of the Higgs boson or associated vector bosons, and have . The sums of weights in each sample are first normalized to the inclusive NLO cross section. The ratio of the distribution from the nnlops generator to that from the powheg or MadGraph 5_amc@nlo generators in each jet bin is applied to the signal samples. The reweighting procedure has been checked against fully simulated nnlops samples in the and decay channels and was found to give results compatible within the statistical uncertainty of the simulated samples. The and boosted analyses, which are much less sensitive to production than other decay channels, use the NLO + PS simulation.

Theoretical uncertainties in gluon fusion

The cross section uncertainty scheme for the and decay channels has been updated to the one proposed in Ref. [53], as already used in the and analyses. This uncertainty scheme includes 9 nuisance parameters accounting for uncertainties in the cross section prediction for exclusive jet bins (including the migration between the 0- and 1-jet, as well as between the 1- and 2-jet bins), the 2-jet and 3-jet phase spaces, different regions, and the uncertainty in the distribution due to missing higher-order finite top quark mass corrections. The boosted search, which is only sensitive to in the high tail, uses a dedicated prediction in this region, and hence the theoretical uncertainties assigned are assumed to be uncorrelated with the other analyses.

Statistical uncertainties in simulation

In the combination, many of the nuisance parameters originate from the use of a limited number of Monte Carlo events to determine SM signal and background expectations. Some of the input analyses have been modified to use the “Barlow-Beeston lite” approach, which assigns a single nuisance parameter per bin that scales the total bin yield [86, 87]. This differs from the previous implementation, which utilized separate nuisance parameters for each process per bin. With the Barlow-Beeston approach, the maximum likelihood estimator for each of these nuisance parameters is independent from the others, and can be solved for analytically. This has been found to provide a significant reduction in the minimization time, while reproducing the results obtained with the full treatment to within 1%.

Combination procedure

The overall statistical methodology used in this combination is the same as the one developed by the ATLAS and CMS Collaborations, and described in Ref. [56]. The procedures used in this paper are described in more detail in Refs. [14, 88, 89] and are based on the standard LHC data modeling and handling toolkits RooFit [90] and RooStats [91].

The parameters of interest (POI) for a particular model are estimated with their corresponding confidence intervals using a profile likelihood ratio test statistic  [92], in which experimental or theoretical uncertainties are incorporated via nuisance parameters (NP) :The likelihood functions in the numerator and denominator of Eq. (1) are constructed using products of probability density functions of signal and background for the various discriminating variables used in the input analyses, as well as constraint terms for certain NPs. The probability density functions are derived from simulation for the signal and from both data and simulation for the background. The quantities and denote the unconditional maximum likelihood estimates of the parameter values, while denotes the conditional maximum likelihood estimate for fixed values of the parameters of interest . The choice of the POIs, e.g., signal strengths (), couplings modifiers, production cross sections, branching fractions or related ratios of the above quantities, depends on the specific model under consideration, while the remaining parameters are treated as NPs. An individual NP represents a single source of systematic uncertainty, and its effect is therefore considered fully correlated between all of the input analyses included in the fit.

For each model considered, the maximum likelihood estimates are identified as the best fit parameter values. The and confidence level () intervals for one-dimensional measurements of each POI are determined as the union of intervals for which and , respectively, unless otherwise stated. In models with more than one POI, these intervals are determined treating the other POIs as NPs. The differences between the boundaries of the and intervals and the best fit value yield the and uncertainties on the measurement. In cases where a physical boundary restricts the interval, we report a truncated interval and determine the uncertainty from that interval. (See Fig. 6 and Table 3, for example). In these cases, the intervals are not expected to maintain coverage. In the case where the intervals are not contiguous, the interval that contains the best fit point is used to determine these uncertainties. The 2D and regions are determined from the set of parameter values for which and , respectively, unless otherwise stated.

Fig. 6

Summary plot of the fit to the production–decay signal strength products . The points indicate the best fit values while the horizontal bars indicate the intervals. The hatched areas indicate signal strengths that are restricted to nonnegative values as described in the text

Table 3

Best fit values and uncertainties for the parameters of the model with one signal strength parameter for each production and decay mode combination. The entries in the table represent the parameter , where i is indicated by the column and f by the row. The expected uncertainties are given in brackets. Some of the signal strengths are restricted to nonnegative values, as described in the text

Decay mode	Production process
	ggH				VBF				WH				ZH				ttH
	Best fit		Uncertainty		Best fit		Uncertainty		Best fit		Uncertainty		Best fit		Uncertainty		Best fit		Uncertainty
	value		stat	syst	value		stat	syst	value		stat	syst	value		stat	syst	value		stat	syst
H→bb	2.51	-2.01+2.43	-1.92+1.96	-0.59+1.44	—				1.73	-0.68+0.70	-0.51+0.53	-0.45+0.46	0.99	-0.45+0.47	-0.40+0.41	-0.20+0.23	0.91	-0.43+0.45	-0.23+0.24	-0.36+0.38
		(-1.86+2.06)	(-1.83+1.85)	(-0.33+0.89)	—					(-0.67+0.69)	(-0.51+0.52)	(-0.44+0.45)		(-0.44+0.46)	(-0.39+0.40)	(-0.20+0.23)		(-0.42+0.44)	(-0.23+0.23)	(-0.35+0.37)
H→ττ	1.05	-0.47+0.53	-0.25+0.25	-0.40+0.46	1.12	-0.43+0.45	-0.35+0.37	-0.25+0.25	—				—				0.23	-0.88+1.03	-0.71+0.80	-0.52+0.65
		(-0.41+0.45)	(-0.23+0.23)	(-0.34+0.38)		(-0.43+0.45)	(-0.35+0.37)	(-0.24+0.25)	—				—					(-0.87+0.98)	(-0.73+0.80)	(-0.47+0.56)
H→WW	1.35	-0.19+0.21	-0.12+0.12	-0.15+0.17	0.28	-0.60+0.64	-0.53+0.58	-0.28+0.28	3.91	-2.01+2.26	-1.72+1.89	-1.05+1.24	0.96	-1.46+1.81	-1.44+1.74	-0.22+0.50	1.60	-0.59+0.65	-0.39+0.40	-0.45+0.52
		(-0.16+0.17)	(-0.10+0.10)	(-0.12+0.13)		(-0.58+0.62)	(-0.53+0.57)	(-0.25+0.26)		(-1.19+1.47)	(-1.06+1.32)	(-0.54+0.64)		(-1.37+1.67)	(-1.35+1.61)	(-0.20+0.45)		(-0.53+0.56)	(-0.38+0.38)	(-0.37+0.41)
H→ZZ	1.22	-0.21+0.23	-0.19+0.20	-0.10+0.12	-0.09	-0.76+1.02	-0.72+1.00	-0.22+0.21	0.00	+0.00+2.33	-0.00+2.31	-0.00+0.30	0.00	+0.00+4.26	-0.00+4.19	-0.00+0.80	0.00	+0.00+1.50	-0.00+1.47	-0.00+0.30
		(-0.20+0.22)	(-0.19+0.20)	(-0.07+0.10)		(-0.99+1.27)	(-0.97+1.25)	(-0.21+0.23)		(-0.99+4.46)	(-0.99+4.42)	(-0.00+0.57)		(-1.00+7.57)	(-1.00+7.45)	(-0.00+1.33)		(-0.99+2.95)	(-0.99+2.89)	(-0.00+0.59)
H→γγ	1.16	-0.18+0.21	-0.15+0.17	-0.10+0.13	0.67	-0.46+0.59	-0.42+0.49	-0.18+0.32	3.76	-1.35+1.48	-1.33+1.45	-0.24+0.33	0.00	+0.00+1.14	-0.00+1.14	-0.00+0.09	2.18	-0.75+0.88	-0.74+0.82	-0.14+0.32
		(-0.16+0.17)	(-0.14+0.14)	(-0.08+0.11)		(-0.48+0.59)	(-0.43+0.48)	(-0.21+0.34)		(-1.16+1.28)	(-1.16+1.27)	(-0.06+0.13)		(-1.04+2.51)	(-1.04+2.50)	(-0.00+0.25)		(-0.63+0.74)	(-0.63+0.72)	(-0.06+0.16)
H→μμ	0.31	-1.79+1.80	-1.78+1.79	-0.19+0.19	2.72	-7.03+7.12	-7.04+7.12	-0.00+0.26	—				—				—
		(-1.65+1.69)	(-1.67+1.67)	(-0.00+0.28)		(-6.94+7.02)	(-6.93+7.01)	(-0.50+0.38)	—				—				—

The likelihood functions are constructed with respect to either the observed data or an Asimov data set [92] constructed using the expected values of the POIs for the SM, in order to obtain the observed or expected results, respectively. Because of fluctuations in the observed data the observed intervals may differ from the expected ones.

Finally, the SM predictions for the production and decay rates of the Higgs boson depend on the mass of the Higgs boson, . For all measurements in this paper, the mass is taken to be , determined from the ATLAS and CMS combined measurement, from the LHC Run 1 data, using the high-resolution and decay channels [93].

Systematic uncertainties

For many of the POIs, the systematic uncertainties in their determination are expected to be as large as, or larger than, the statistical uncertainties. The theoretical uncertainties affecting the signal are among the most important contributions to the systematic uncertainties. The uncertainties in the total cross section prediction for the signal processes arising from the parton distribution functions, the renormalization and factorization scales used in the calculations and the branching fraction predictions are correlated between all input analyses. Instead, theoretical uncertainties that affect kinematic distributions and cause migrations between event categories are largely uncorrelated between the input analyses. An exception is the set of theoretical uncertainties for the production mode, where the correlation scheme described in Sect. 4.2 is used to correlate both the normalization and shape uncertainties between input analyses. The theoretical uncertainties affecting the background predictions, including the parton distribution function uncertainties, are assumed to be uncorrelated with those affecting the signal predictions [88], with the exception of the uncertainties from the underlying event and parton shower model.

The majority of the systematic uncertainties arising from experimental sources are uncorrelated between the input analyses, with a few exceptions. The uncertainties in the integrated luminosity measurement [94], and in the modeling of additional collisions in the event (pileup), are correlated between all of the input input analyses. Certain analyses, namely the , , and analyses, are able to further constrain the jet energy scale uncertainties determined in auxiliary measurements. The jet energy scale uncertainty in these analyses is decomposed into several nuisance parameters corresponding to different sources of uncertainty (for example, different flavor composition and kinematic regions) that are correlated among these analyses but uncorrelated with the other analyses. An independent jet energy scale uncertainty is assumed to be correlated between the input analyses that are not sensitive to the different sources of uncertainty. The uncertainties in the tagging efficiency are correlated between the analyses, but are uncorrelated from the analysis, which is sensitive to different kinematic regions. A separate set of NPs is used to describe the uncertainty in the tagging efficiency in the , , and analyses. The uncertainty in the efficiency of the double--tagger algorithm described in Sect. 3.6 is taken to be uncorrelated from the single tagging uncertainties. Finally, the uncertainties in the lepton efficiency and misidentification rate in the - and - event classes are correlated, since the same reconstruction and identification algorithms were used. In other input analyses, different algorithms were used and therefore the uncertainties are assumed to be uncorrelated.

The free parameters describing the shapes and normalizations of the background models, and parameters that allow for the choice of the background parametrization in each of the analysis categories are fully determined by the data without any additional constraints, and are therefore assigned to the statistical uncertainty of a measurement. The remaining uncertainties are assigned to groups of systematic uncertainties.

Signal strength and cross section fits

The signal strength modifier , defined as the ratio between the measured Higgs boson yield and its SM expectation, has been used extensively to characterize the Higgs boson yields. However, the specific meaning of varies depending on the analysis. For a specific production and decay channel , the signal strengths for the production, , and for the decay, , are defined as:respectively. Here and are, respectively, the production cross section for and the branching fraction for . The subscript ”SM” refers to their respective SM predictions, so by definition, the SM corresponds to . Since and cannot be separately measured without additional assumptions, only the product of and can be extracted experimentally, leading to a signal strength  for the combined production and decay:This parametrization makes use of the narrow width approximation, and the reliability of this approximation was studied in Ref. [95] and found to be adequate for global fits.

In this section, results are presented for several signal strength parametrizations starting with a single global signal strength , which is the most restrictive in terms of the number of assumptions. Further parametrizations are defined by relaxing the constraint that all production and decay rates scale with a common signal strength modifier.

The combined measurement of the common signal strength modifier at is,where the total uncertainty has been decomposed into statistical, signal theoretical systematic, and other systematic components. The largest single source of uncertainty apart from the signal theoretical systematic uncertainties is the integrated luminosity (), which is correlated between all of the input analyses. In this measurement and others, however, the other systematic uncertainty component is mostly dominated by uncertainties that only affect a single input analysis.

Relaxing the assumption of a common production mode scaling, but still assuming the relative SM branching fractions, leads to a parametrization with five production signal strength modifiers: , , , , and . In this parametrization, as well as all subsequent parametrizations involving signal strengths or cross sections, the production is assumed to scale like . Conversely, relaxing the common decay mode scaling, but assuming the relative SM production cross sections, leads to one with the modifiers: , , , , , and . Results of the fits in these two parametrizations are summarized in Fig. 5. The numerical values, including the decomposition of the uncertainties into statistical and systematic components, and the corresponding expected uncertainties, are given in Table 2.

Fig. 5

Summary plot of the fit to the per-production mode (left) and per-decay mode (right) signal strength modifiers. The thick and thin horizontal bars indicate the and uncertainties, respectively. Also shown are the systematic components of the uncertainties. The last point in the per-production mode summary plot is taken from a separate fit and indicates the result of the combined overall signal strength

Table 2

Best fit values and uncertainties for the parametrizations with per-production mode and per-decay mode signal strength modifiers. The expected uncertainties are given in brackets

Production process	Best fit value		Uncertainty
			stat.	syst.
ggH	1.22	-0.12+0.14	-0.08+0.08	-0.10+0.12
		(-0.11+0.11)	(-0.07+0.07)	(-0.08+0.09)
VBF	0.73	-0.27+0.30	-0.23+0.24	-0.15+0.17
		(-0.27+0.29)	(-0.23+0.24)	(-0.15+0.16)
WH	2.18	-0.55+0.58	-0.45+0.46	-0.32+0.34
		(-0.51+0.53)	(-0.42+0.43)	(-0.29+0.30)
ZH	0.87	-0.42+0.44	-0.38+0.39	-0.18+0.20
		(-0.41+0.43)	(-0.37+0.38)	(-0.17+0.19)
ttH	1.18	-0.27+0.30	-0.16+0.16	-0.21+0.26
		(-0.25+0.28)	(-0.15+0.16)	(-0.20+0.23)

The improvement in the precision of the measurement of the production rate of 50% (from 20% to 10%) compared to Ref. [55] and 33% (from 15% to 10%) compared to Ref. [56], can be attributed to the combined effects of an increased production cross section, and a reduction in the associated theoretical uncertainties. The improvements in the precision are up to 20% for the and production rates compared to Ref. [55]. The uncertainty in the measurement of the production rate is reduced by around 50% compared to Ref. [56]. This is in part due to the increase in the cross section between 8 and 13, but also due to the inclusion of additional exclusive event categories for this production process.

The most generic signal strength parametrization has one signal strength parameter for each production and decay mode combination, . Given the five production and six decay modes listed above, this implies a model with 30 parameters of interest. However not all can be experimentally constrained in this combination. There is no dedicated analysis from CMS at targeting and production with decay, or production with decay, therefore these signal strength modifiers are fixed to the SM expectation and are not included in the maximum likelihood fit. Likewise, the , , and production rates with decay are fixed to the SM expectation. In the case of , , and production with decay, as well as production with decay, the background contamination is sufficiently low so that a negative signal strength can result in an overall negative event yield. Therefore, these signal strengths are restricted to nonnegative values. Figure 6 summarizes the results in this model along with the intervals. The numerical values, including the uncertainty decomposition into statistical and systematic parts, and the corresponding expected uncertainties, are given in Table 3.

Ratios of cross sections and branching fractions, relative to

Results are presented for a model based on the ratios of cross sections and branching fractions. These are given relative to a well-measured reference process, chosen to be (). Using ratios has the advantage that some systematic or theoretical uncertainties common to both the numerator and denominator cancel. The following ratios are used: , , ,, , , , , and . These results are summarized in Fig. 7, and the numerical values are given in Table 4. The uncertainties in the SM predictions are included in the measurements.

Fig. 7

Summary of the cross section and branching fraction ratio model. The thick and thin horizontal bars indicate the and uncertainties, respectively. Also shown are the systematic components of the uncertainties

Table 4

Best fit values and uncertainties for the parameters of the cross section and branching fraction ratio model. The expected uncertainties are given in brackets

Parameter	Best fit		Uncertainty		Parameter	Best fit		Uncertainty
			stat	syst				stat	syst
μggHZZ	1.07	-0.18+0.20	-0.15+0.16	-0.09+0.11	B^bb/B^ZZ	0.84	-0.27+0.37	-0.21+0.27	-0.17+0.26
		(-0.16+0.19)	(-0.14+0.15)	(-0.08+0.11)			(-0.37+0.56)	(-0.28+0.38)	(-0.25+0.40)
μ_VBF/μ_ggH	0.60	-0.24+0.30	-0.21+0.24	-0.13+0.17	B^ττ/B^ZZ	1.07	-0.30+0.37	-0.21+0.25	-0.22+0.27
		(-0.32+0.40)	(-0.27+0.31)	(-0.17+0.24)			(-0.28+0.35)	(-0.20+0.25)	(-0.19+0.25)
μ_WH/μ_ggH	2.19	-0.69+0.86	-0.56+0.68	-0.39+0.52	B^WW/B^ZZ	1.23	-0.22+0.27	-0.18+0.22	-0.13+0.16
		(-0.52+0.65)	(-0.44+0.53)	(-0.29+0.39)			(-0.19+0.24)	(-0.16+0.19)	(-0.11+0.14)
μ_ZH/μ_ggH	0.88	-0.44+0.57	-0.41+0.49	-0.17+0.30	B^γγ/B^ZZ	1.14	-0.20+0.28	-0.18+0.23	-0.09+0.16
		(-0.47+0.68)	(-0.41+0.53)	(-0.23+0.43)			(-0.18+0.23)	(-0.16+0.20)	(-0.08+0.11)
μ_ttH/μ_ggH	1.06	-0.27+0.34	-0.18+0.20	-0.20+0.27	B^μμ/B^ZZ	0.63	-1.21+1.24	-1.20+1.24	-0.11+0.15
		(-0.30+0.36)	(-0.21+0.23)	(-0.21+0.27)			(-1.19+1.26)	(-1.19+1.25)	(-0.03+0.20)

Stage 0 simplified template cross sections

Measurements of production cross sections, which are complementary to the signal strength parametrization, are made for seven processes defined according to the simplified template cross sections proposed in Ref. [53]. The results given here are for the stage 0 fiducial regions defined by the rapidity of the Higgs boson . All input analyses have a negligible acceptance for . Defining the fiducial region in this way reduces the theoretical uncertainty that would otherwise apply while extrapolating to the fully inclusive phase space. Subsequent stages involve splitting the fiducial regions into a number of smaller ones, for example based on ranges of the Higgs boson . The measured cross sections are defined as:In addition to the cross sections, the Higgs boson branching fractions are also included as POIs via ratios with respect to . A summary of the results in this model, normalized to the expected SM cross sections, is given in Fig. 8 and Table 5. Since cross sections are measured and not signal strength modifiers, the theoretical uncertainties in these cross sections do not enter as sources of uncertainty. In Fig. 8, the uncertainties in the SM predictions are indicated by gray bands.

Fig. 8

Summary of the stage 0 model, ratios of cross sections and branching fractions. The points indicate the best fit values, while the error bars show the and uncertainties. The uncertainties on the measurements considering only the contributions from the systematic uncertainties are also shown. The uncertainties in the SM predictions are indicated

Table 5

Best fit values and uncertainties for the parameters of the stage 0 simplified template cross section model. The values are all normalized to the SM predictions. The expected uncertainties are given in brackets

Parameter	Best fit		Uncertainty		Parameter	Best fit		Uncertainty
			stat	syst				stat	syst
σ_ggHB^ZZ	1.00	-0.16+0.19	-0.15+0.16	-0.07+0.09	B^bb/B^ZZ	0.96	-0.31+0.44	-0.24+0.32	-0.20+0.30
		(-0.16+0.18)	(-0.15+0.16)	(-0.07+0.09)			(-0.38+0.57)	(-0.29+0.40)	(-0.25+0.41)
σ_VBFB^ZZ	0.66	-0.26+0.32	-0.22+0.27	-0.13+0.17	B^ττ/B^ZZ	0.98	-0.28+0.35	-0.20+0.24	-0.20+0.25
		(-0.32+0.40)	(-0.27+0.33)	(-0.16+0.22)			(-0.29+0.36)	(-0.21+0.26)	(-0.19+0.25)
σ_H+V(qq)B^ZZ	3.93	-1.71+2.00	-1.53+1.77	-0.75+0.93	B^WW/B^ZZ	1.30	-0.24+0.29	-0.20+0.24	-0.13+0.17
		(-1.05+1.66)	(-1.05+1.49)	(-0.00+0.72)			(-0.20+0.24)	(-0.16+0.20)	(-0.11+0.14)
σ_H+W(ℓν)B^ZZ	1.95	-0.68+0.88	-0.57+0.72	-0.38+0.51	B^γγ/B^ZZ	1.14	-0.20+0.26	-0.18+0.23	-0.09+0.13
		(-0.52+0.69)	(-0.44+0.56)	(-0.29+0.40)			(-0.19+0.23)	(-0.17+0.21)	(-0.08+0.11)
σ_H+Z(ℓℓ/νν)B^ZZ	0.84	-0.43+0.57	-0.40+0.49	-0.17+0.29	B^μμ/B^ZZ	0.67	-1.36+1.40	-1.35+1.39	-0.13+0.18
		(-0.46+0.71)	(-0.41+0.56)	(-0.22+0.44)			(-1.28+1.35)	(-1.28+1.34)	(-0.05+0.17)
σ_ttHB^ZZ	1.08	-0.30+0.37	-0.22+0.26	-0.19+0.26
		(-0.31+0.38)	(-0.23+0.28)	(-0.20+0.26)

: gluon fusion and -associated production. While Ref. [53] proposes separate bins for these modes, they are merged here because of the current lack of sensitivity to the associated production with quarks.

: production.

: Associated production with a or boson, either quark or gluon initiated, in which the vector boson decays hadronically.

: Associated production with a boson, in which the boson decays leptonically. While Ref. [53] proposes separate bins for the quark- and gluon-initiated modes, they are merged here because they cannot easily be distinguished experimentally, and therefore, their measurements would be highly anticorrelated.

: Associated production with a boson, in which the decays leptonically.

: Associated production with a pair of top quarks or a single top quark. While Ref. [53] proposes separate bins for these modes, they are merged here because of the lack of a dedicated analysis targeting production in this combination.

Measurements of Higgs boson couplings

In the -framework [96], coupling modifiers are introduced in order to test for deviations in the couplings of the Higgs boson to other particles. In order to measure the individual Higgs couplings in this framework, some assumption must be made to constrain the total Higgs boson width since it cannot be directly measured at the LHC. Unless stated otherwise, it is assumed that there are no BSM contributions to the total Higgs boson width. With this assumption, the cross section times branching fraction for a production process i and decay f can be expressed as,where is the total width of the Higgs boson and is the partial width of the Higgs boson decay to the final state f. A set of coupling modifiers, , is introduced to parameterize potential deviations in the bosonic and fermionic couplings of the Higgs boson from the SM predictions. For a given production process or decay mode j, a coupling modifier is defined such that,In the SM, all values are positive and equal to unity. In this parametrization it is assumed that the higher-order accuracy of the QCD and electroweak corrections to the SM cross sections and branching fractions is preserved when the values of deviate from unity. While this does not hold in general, for the parameter ranges considered in this paper the dominant higher-order QCD corrections largely factorize from the rescaling of the couplings, therefore the approach is considered valid. Individual coupling modifiers, corresponding to tree-level Higgs boson couplings to the different particles, are introduced, as well as effective coupling modifiers and that describe production and decay. This is possible because the presence of any BSM particles in these loops is not expected to significantly change the corresponding kinematic properties of the processes. This approach is not possible for production, which occurs at leading order through box and triangular loop diagrams, because a tree-level contact interaction from BSM physics would likely exhibit a kinematic structure very different from the SM, and is expected to be highly suppressed [97]. Other possible BSM effects on the process are related to modifications of the and vertices, which are best taken into account, within the limitation of the framework, by resolving the loop in terms of the corresponding coupling modifiers, and . More details on the development of this framework as well as its theoretical and phenomenological foundations and extensions can be found, for example, in Refs. [98-112].

The normalization scaling effects of each of the parameters are given in Table 6. Loop processes such as and can be studied through either the effective coupling modifiers, thereby providing sensitivity to potential BSM physics in the loops, or the modifiers of the SM particles themselves. Interference between different diagrams, such as those that contribute to , provides some sensitivity to relative signs between Higgs boson couplings to different particles. Modifications to the kinematic distributions of the production are also expected when the relative sign of and is negative. These effects were studied and the distributions of the final observables were found to be insensitive with the present dataset to the relative sign of and .

Table 6

Normalization scaling factors for all relevant production cross sections and decay partial widths. For the parameters representing loop processes, the resolved scaling in terms of the fundamental SM couplings is also given

Effective
	Loops	Interference	scaling factor	Resolved scaling factor
Production
σ(ggH)	✓	g-t	κg2	1.04κt2+0.002κb2-0.038κ_tκ_b
σ(VBF)	—	—		0.73κW2+0.27κZ2
σ(WH)	—	—		κW2
σ(qq/qg→ZH)	—	—		κZ2
σ(gg→ZH)	✓	Z-t		2.46κZ2+0.47κt2-1.94κ_Zκ_t
σ(ttH)	—	—		κt2
σ(gb→WtH)	—	W-t		2.91κt2+2.31κW2-4.22κ_tκ_W
σ(qb→tHq)	—	W-t		2.63κt2+3.58κW2-5.21κ_tκ_W
σ(bbH)	—	—		κb2
Partial decay width
Γ^ZZ	—	—		κZ2
Γ^WW	—	—		κW2
Γ^γγ	✓	W-t	κγ2	1.59κW2+0.07κt2-0.67κ_Wκ_t
Γ^ττ	—	—		κτ2
Γ^bb	—	—		κb2
Γ^μμ	—	—		κμ2
Total width for B_BSM=0
				0.58κb2+0.22κW2+0.08κg2+
Γ_H	✓	—	κH2	+0.06κτ2+0.026κZ2+0.029κc2+
				+0.0023κγ2+0.0015κZγ2+
				+0.00025κs2+0.00022κμ2

Generic model within -framework assuming resolved loops

Under the assumption that there are no BSM particles contributing to the production or decay loops, these processes can be expressed in terms of the coupling modifiers to the SM particles as described previously. There are six free coupling parameters: , , , , , and . Without loss of generality, the value of is restricted to be positive, while both negative and positive values of , and are allowed. In this model, the rates of the and processes, which occur through loop diagrams at leading order, are resolved, meaning that they are described by the functions of , , , and given in Table 6. The results of the fits with this parametrization are given in Fig. 9 and Table 7.

Fig. 9

Summary of the -framework model assuming resolved loops and . The points indicate the best fit values while the thick and thin horizontal bars show the and intervals, respectively. In this model, the and loops are resolved in terms of the remaining coupling modifiers. For this model, both positive and negative values of , , and are considered. Negative values of in this model are disfavored by more than

Table 7

Best fit values and uncertainties for the parameters of the model in which the loop processes are resolved. The expected uncertainties are given in brackets

Parameter	Best fit value		Uncertainty
			stat.	syst.
κ_W	1.10	-0.17+0.12	-0.16+0.08	-0.06+0.08
		(-0.10+0.11)	(-0.08+0.08)	(-0.06+0.06)
κ_Z	0.99	-0.12+0.11	-0.10+0.09	-0.07+0.06
		(-0.11+0.11)	(-0.09+0.09)	(-0.06+0.06)
κ_t	1.11	-0.10+0.12	-0.07+0.07	-0.08+0.09
		(-0.12+0.11)	(-0.08+0.07)	(-0.09+0.09)
κ_b	-1.10	-0.23+0.33	-0.16+0.29	-0.17+0.15
		(-0.22+0.22)	(-0.15+0.15)	(-0.16+0.17)
κ_τ	1.01	-0.20+0.16	-0.17+0.11	-0.11+0.12
		(-0.15+0.17)	(-0.10+0.12)	(-0.11+0.12)
κ_μ	0.79	-0.79+0.58	-0.80+0.56	-0.00+0.14
		(-1.01+0.50)	(-1.01+0.50)	(-0.10+0.08)

The rate of the decay and production depend only on the absolute value of . The interference between the two diagrams shown in Fig. 3, however, allows contributions from the production mode to break the degeneracy between the signs, leading to a positive value of being preferred. As these contributions are typically small compared to other production modes, the 1 and 2 intervals also include negative values of . Although a negative value of is preferred in this model, the difference in q between the best fit point and the minimum in the region is smaller than 0.1.

An additional fit is performed using a phenomenological parametrization relating the masses of the fermions and vector bosons to the corresponding modifiers using two parameters, denoted M and  [113, 114]. In such a model one can relate the coupling modifiers to M and as for fermions and for vector bosons. Here, , is the SM Higgs boson vacuum expectation value [115]. The SM expectation, , is recovered when .

The lepton and vector boson mass values are taken from Ref. [115], while the top quark mass is taken to be 172.5 for consistency with theoretical calculations used in setting the SM predictions. The bottom quark mass is evaluated at the scale of the Higgs boson mass, .

The and regions in the fit are shown in Fig. 10 (left). The results of the fit using the six parameter model are plotted versus the particle masses in Fig. 10 (right), and the result of the fit is also shown for comparison. For the quark, since the best fit point for is negative, the absolute value of this coupling modifier is shown. In order to show both the Yukawa and vector boson couplings in the same plot, a “reduced” vector boson coupling is shown.

Fig. 10

Likelihood scan in the M- plane (left). The best fit point and the and regions are shown, along with the SM prediction. Result of the phenomenological fit overlayed with the resolved -framework model (right)

Generic model within -framework with effective loops

The results of the fits to the generic model where the and loops are scaled using the effective coupling modifiers and are given in Fig. 11 and Table 8. In this parametrization, additional contributions from BSM decays are allowed for by rewriting the total width of the Higgs boson, relative to its SM value, as,where is defined in Table 6.

Fig. 11

Summary plots for the -framework model in which the and loops are scaled with effective couplings. The points indicate the best fit values while the thick and thin horizontal bars show the and intervals, respectively. In the left figure the constraint is imposed, and both positive and negative values of and are considered. In the right figure a constraint is imposed (same sign of and ), while and are free parameters

Table 8

Best fit values and uncertainties for the parameters of the -framework model with effective loops. The expected uncertainties are given in brackets

B_BSM=0					B_BSM>0,  |κ_V|<1
Parameter	Best fit		Uncertainty		Parameter	Best fit		Uncertainty
			stat	syst				stat	syst
κ_Z	1.00	-0.11+0.11	-0.09+0.09	-0.07+0.06	κ_Z	-0.87	-0.08+0.08	-0.06+0.07	-0.04+0.04
		(-0.11+0.11)	(-0.09+0.09)	(-0.06+0.06)			(-0.12+0.00)	(-0.10+0.00)	(-0.06+0.00)
κ_W	-1.13	-0.13+0.16	-0.10+0.15	-0.08+0.06	κ_W	-1.00	-0.00+0.09	-0.00+0.07	-0.00+0.05
		(-0.12+0.12)	(-0.09+0.09)	(-0.07+0.07)			(-0.12+0.00)	(-0.09+0.00)	(-0.07+0.00)
κ_t	0.98	-0.14+0.14	-0.08+0.08	-0.11+0.12	κ_t	1.02	-0.15+0.19	-0.09+0.13	-0.13+0.13
		(-0.15+0.14)	(-0.09+0.08)	(-0.12+0.12)			(-0.15+0.18)	(-0.09+0.13)	(-0.12+0.13)
κ_τ	1.02	-0.17+0.17	-0.13+0.11	-0.10+0.12	κ_τ	0.93	-0.13+0.13	-0.09+0.08	-0.10+0.11
		(-0.15+0.16)	(-0.11+0.11)	(-0.11+0.12)			(-0.15+0.14)	(-0.10+0.09)	(-0.11+0.11)
κ_b	1.17	-0.31+0.27	-0.29+0.18	-0.10+0.20	κ_b	0.91	-0.16+0.17	-0.12+0.11	-0.11+0.13
		(-0.23+0.25)	(-0.17+0.18)	(-0.16+0.17)			(-0.22+0.19)	(-0.16+0.14)	(-0.15+0.13)
κ_g	1.18	-0.14+0.16	-0.09+0.10	-0.10+0.12	κ_g	1.16	-0.13+0.18	-0.09+0.14	-0.10+0.12
		(-0.12+0.14)	(-0.09+0.10)	(-0.09+0.10)			(-0.12+0.17)	(-0.09+0.13)	(-0.09+0.11)
κ_γ	1.07	-0.15+0.14	-0.14+0.10	-0.05+0.09	κ_γ	0.96	-0.09+0.09	-0.08+0.06	-0.06+0.06
		(-0.12+0.12)	(-0.09+0.10)	(-0.07+0.07)			(-0.11+0.09)	(-0.09+0.07)	(-0.07+0.05)
κ_μ	0.80	-0.80+0.59	-0.81+0.56	-0.00+0.17	κ_μ	0.72	-0.72+0.50	-0.71+0.50	-0.07+0.00
		(-1.01+0.51)	(-1.01+0.50)	(-0.09+0.09)			(-1.01+0.49)	(-1.00+0.48)	(-0.08+0.06)
					B_inv	0.07	-0.07+0.08	-0.03+0.03	-0.06+0.07
							(+0.00+0.09)	(-0.00+0.04)	(-0.00+0.08)
					B_undet	0.00	+0.00+0.17	-0.00+0.14	-0.00+0.09
							(+0.00+0.20)	(-0.00+0.17)	(-0.00+0.11)

Two different model assumptions are made concerning the BSM branching fraction. In the first parametrization, it is assumed that , whereas in the second, and are allowed to vary as POIs, and instead the constraint is imposed. This avoids a complete degeneracy in the total width where all of the coupling modifiers can be scaled equally to account for a non-zero . The parameter represents the total branching fraction to any final state that is not detected by the analyses included in this combined analysis. The likelihood scan for the parameter in this model, and the 2D likelihood scan of vs. are given in Fig. 12. The 68 and 95% regions for Fig. 12 (right) are determined as the regions for which and 5.99, respectively. The 95% upper limits of and are determined, corresponding to the value for which  [92]. The uncertainty in the measurement of is reduced by nearly 40% compared to Ref. [56]. This improvement is because of the improved sensitivity to the production mode as described in Section 7.

Fig. 12

Results within the generic -framework model with effective loops and with the constraint (same sign of and ), and with and as free parameters. Scan of the test statistic q as a function of (left), and 68 and 95% regions for vs. (right). The scan of the test statistic q as a function of expected assuming the SM is also shown in the left figure

In both of the generic models, the best fit point for is negative. The value of as a function of in the two cases is shown in Fig. 13. While different combinations of signs for and are shown, the minimum value of q across all combinations is used to determine the best fit point and the and regions.

Fig. 13

Scan of the test statistic q as a function of in the generic model assuming (left) and allowing and to float (right). The different colored lines indicate the value of q for different combinations of signs for and . The solid black line shows the minimum value of in each case and is used to determine the best fit point and the and regions. The scan in the right figure is truncated because of the constraints of and , which are imposed in this model

The preferred negative value of is due to the interference between some of the diagrams describing production, which contributes in several analyses entering the combination. In particular, the excess in the tagged categories of the analysis can be accounted for by a negative value of as this increases the contribution of production. In these models, the decay is treated as an effective coupling so that it has no dependence on . This means that a negative value of will not result in excesses in the other categories of the analysis.

Using Eq. (7), this model is also reinterpreted as a constraint on the total Higgs boson width, and the corresponding likelihood scan is shown in Fig. 14. Using this parametrization, the total Higgs boson width relative to the SM expectation is determined to be . The different behavior between the observed and expected likelihood scans for large is due to the preference in data for the relative sign combination.

Fig. 14

The scan of the test statistic q as a function of obtained by reinterpreting the model allowing for BSM decays of the Higgs boson. The expected scan of q as a function of assuming the SM is also shown

An additional fit is performed assuming that the only BSM contributions to the Higgs couplings appear in the loop-induced and processes. In this fit, and are the POIs, and are floated, and the other couplings are fixed to their SM predictions. The best fit point and the and regions in the - plane for this model are shown in Fig. 15.

Fig. 15

The and regions in the vs. parameter space for the model assuming the only BSM contributions to the Higgs boson couplings appear in the loop-induced processes or in BSM Higgs decays

Generic model with effective loops and coupling modifier ratios

An analogous parametrization to the ratios of cross sections and branching fractions described in the previous section can be derived in terms of ratios of the coupling modifiers (). In this parametrization a reference combined coupling modifier is defined that accounts for modifications to the total event yield of a specific production times decay process, thereby avoiding the need for assumptions on the total Higgs boson width. The reference coupling modifier is taken to be . The remaining parameters of interest are ratios of the form: , , , , , . A summary of the results in this model is given in Fig. 16, and the numerical values along with the  uncertainties are shown in Table 9.

Fig. 16

Summary of the model with coupling ratios and effective couplings for the and loops. The points indicate the best fit values while the thick and thin horizontal bars show the and intervals, respectively. For this model, both positive and negative values of and are considered

Table 9

Best fit values and uncertainties for the parameters of the coupling modifier ratio model. The expected uncertainties are given in brackets

Parameter	Best fit		Uncertainty		Parameter	Best fit		Uncertainty
			stat	syst				stat	syst
κ_gZ	1.03	-0.09+0.09	-0.07+0.07	-0.05+0.05	λ_γZ	1.07	-0.10+0.12	-0.08+0.10	-0.05+0.06
		(-0.09+0.09)	(-0.07+0.07)	(-0.05+0.05)			(-0.09+0.11)	(-0.08+0.09)	(-0.04+0.05)
λ_WZ	-1.13	-0.11+0.10	-0.09+0.08	-0.06+0.06	λ_bZ	1.17	-0.20+0.23	-0.14+0.16	-0.14+0.16
		(-0.09+0.11)	(-0.08+0.09)	(-0.05+0.06)			(-0.19+0.22)	(-0.14+0.16)	(-0.13+0.15)
λ_tg	0.83	-0.13+0.14	-0.08+0.08	-0.10+0.11	λ_τZ	1.02	-0.15+0.16	-0.10+0.11	-0.11+0.12
		(-0.16+0.17)	(-0.11+0.11)	(-0.12+0.12)			(-0.14+0.16)	(-0.10+0.11)	(-0.10+0.11)
λ_Zg	0.85	-0.13+0.14	-0.12+0.10	-0.05+0.09	λ_μZ	0.81	-0.81+0.57	-0.82+0.56	-0.00+0.11
		(-0.16+0.17)	(-0.13+0.13)	(-0.09+0.11)			(-1.01+0.50)	(-1.01+0.49)	(-0.07+0.07)

Fits of vector boson and fermion coupling modifiers

A more constrained version of the loop-resolved model is defined by assuming a common scaling of all vector boson and fermion couplings, respectively. Two models are defined: one in which all signal processes are scaled according to these two and parameters, and one in which separate and parameters are defined for each of the five decay processes. The best fit points and the and regions in the - plane for both models are shown in Fig. 17, and the results are summarized in Table 10. For large values of the likelihood becomes essentially flat, resulting in the best fit point for this parameter being beyond the scale of the axis shown. The 1D 68% region for can be expressed as .

Fig. 17

The and regions in the vs. parameter space for the model assuming a common scaling of all the vector boson and fermion couplings

Table 10

Best fit values and uncertainties for the parameters of the model. The expected uncertainties are given in brackets

Parameter	Best fit		Uncertainty		Parameter	Best fit		Uncertainty
			stat	syst				stat	syst
κVWW	1.10	-0.08+0.08	-0.06+0.06	-0.06+0.06	κFWW	1.49	-0.38+1.55	-0.33+1.04	-0.18+1.15
		(-0.08+0.08)	(-0.06+0.06)	(-0.05+0.05)			(-0.20+0.38)	(-0.17+0.31)	(-0.11+0.21)
κVZZ	0.96	-0.08+0.09	-0.07+0.08	-0.04+0.05	κFZZ	—		—	—
		(-0.11+0.12)	(-0.10+0.11)	(-0.05+0.05)			(-0.31+1.79)	(-0.30+1.55)	(-0.06+0.91)
κVbb	1.15	-0.22+0.23	-0.18+0.18	-0.13+0.15	κFbb	1.01	-0.18+0.16	-0.10+0.09	-0.15+0.13
		(-0.22+0.22)	(-0.18+0.17)	(-0.12+0.13)			(-0.18+0.16)	(-0.10+0.09)	(-0.15+0.13)
κVττ	1.10	-0.29+0.38	-0.24+0.26	-0.18+0.27	κFττ	0.98	-0.16+0.15	-0.09+0.08	-0.14+0.13
		(-0.29+0.32)	(-0.23+0.24)	(-0.17+0.20)			(-0.14+0.14)	(-0.08+0.07)	(-0.12+0.11)
κVγγ	1.10	-0.08+0.14	-0.07+0.11	-0.05+0.09	κFγγ	1.14	-0.29+0.67	-0.26+0.53	-0.14+0.41
		(-0.08+0.10)	(-0.06+0.08)	(-0.04+0.05)			(-0.25+0.47)	(-0.23+0.41)	(-0.10+0.23)

Benchmark models with resolved loops to test the symmetry of fermion couplings

Several BSM models predict the existence of an extended Higgs sector. In such scenarios, the couplings to up-and down-type fermions, or to leptons and quarks, can be separately modified. In order to probe such models, parametrizations are introduced in which the couplings of the Higgs boson to fermions are scaled either by separate common modifiers for up-type () and down-type () fermions or by separate common modifiers for quarks and leptons ().

Figure 18 shows the results of the fits where the ratio of the couplings to up- and down-type fermions is determined along with the ratio and . Also shown are the results of the fit where the ratio of the coupling to leptons and to quarks is determined along with the ratio and . The results of these two parametrizations are summarized in Table 11.

Fig. 18

Summary plots of the 3-parameter models comparing up- and down-type fermions, and floating the ratio of the vector coupling to the up-type coupling (left) and comparing lepton and quark couplings (right). The points indicate the best fit values while the thick and thin horizontal bars show the and intervals, respectively. Both positive and negative values of , , , and are considered

Table 11

Best fit values and uncertainties for the parameters of the two benchmark models with resolved loops to test the symmetry of fermion couplings. The expected uncertainties are given in brackets

λ_Vu				λ_du				κ_uu
Best fit		Uncertainty		Best fit		Uncertainty		Best fit		Uncertainty
value		stat	syst	value		stat	syst	value		stat	syst
0.97	-0.10+0.14	-0.08+0.11	-0.06+0.09	0.92	-0.12+0.12	-0.09+0.09	-0.08+0.08	1.14	-0.20+0.20	-0.16+0.13	-0.12+0.15
	(-0.11+0.15)	(-0.09+0.12)	(-0.07+0.09)		(-0.13+0.13)	(-0.10+0.10)	(-0.08+0.08)		(-0.19+0.16)	(-0.16+0.12)	(-0.11+0.11)

Compatibility of measurements with the SM

Table 12 shows a summary of the compatibility of the different models considered, as described in Sects. 7 and 8, with the SM predictions. For each model, the value of q at the values of the POIs for the SM expectation () is converted to a p-value with respect to the SM. This is done assuming q is distributed according to a function with the number of degrees of freedom equal to the number of POIs. This p-value is found to be greater than 5% for all parametrizations.

Table 12

Compatibility of the fit results with the SM prediction under various signal parametrizations. The value of q at the values of the POIs for which the SM expectation is obtained () is shown along with the corresponding p-value, with respect to the SM, assuming q is distributed according to a function with the specified number of degrees of freedom (DOF)

Parameterization	p-value (q_SM)	DOF	Parameters of interest
Global signal strength	6.28% (3.46)	1	μ
Production processes	9.87% (9.27)	5	μ_ggH, μ_VBF, μ_WH, μ_ZH, μ_ttH
Decay modes	53.8% (5.05)	6	μ^γγ, μ^ZZ, μ^WW, μ^ττ, μ^bb, μ^μμ
σ_iB^f products	61.2% (21.5)	24	σ_ggHB^bb, σ_ggHB^ττ, σ_ggHB^μμ, σ_ggHB^WW, σ_ggHB^ZZ, σ_ggHB^γγ, σ_VBFB^ττ, σ_VBFB^μμ, σ_VBFB^WW, σ_VBFB^ZZ, σ_VBFB^γγ, σ_WHB^bb, σ_WHB^WW, σ_WHB^ZZ, σ_WHB^γγ, σ_ZHB^bb, σ_ZHB^WW, σ_ZHB^ZZ, σ_ZHB^γγ, σ_ttHB^ττ, σ_ttHB^WW, σ_ttHB^ZZ, σ_ttHB^γγ, σ_ttHB^bb
Ratios of σ and B relative to gg→H→ZZ	32.3% (11.5)	10	μggHZZ, μ_VBF/μ_ggH, μ_WH/μ_ggH, μ_ZH/μ_ggH, μ_ttH/μ_ggH, μ^WW/μ^ZZ, μ^γγ/μ^ZZ, μ^ττ/μ^ZZ, μ^bb/μ^ZZ, μ^bb/μ^μμ
Simplified template cross sections with branching fractions relative to B^ZZ	21.2% (14.4)	11	σ_ggHB^ZZ, σ_VBFB^ZZ, σ_H+V(qq)B^ZZ, σ_H+W(ℓν)B^ZZ, σ_H+Z(ℓℓ/νν)B^ZZ, σ_ttHB^ZZ, B^bb/B^ZZ, B^ττ/B^ZZ, B^μμ/B^ZZ, B^WW/B^ZZ, B^γγ/B^ZZ
Couplings, SM loops	45.6% (5.71)	6	κ_Z, κ_W, κ_t, κ_τ, κ_b, κ_μ
Couplings vs. mass	16.8% (3.57)	2	M, ϵ
Couplings, BSM loops	18.5% (11.3)	8	κ_Z, κ_W, κ_t, κ_τ, κ_b, κ_μ, κ_γ, κ_g
Couplings, BSM loops and decays including H→invisible analyses	32.4% (11.5)	10	κ_Z, κ_W, κ_t, κ_τ, κ_b, κ_μ, κ_γ, κ_g, B_inv, B_undet
Ratios of coupling modifiers	18.1% (11.4)	8	κ_gZ, λ_WZ, λ_γZ, λ_tg, λ_bZ, λ_τZ, λ_μZ, λ_Zg
Fermion and vector couplings	16.9% (3.55)	2	κ_F, κ_V
Fermion and vector couplings, per decay mode	76.7% (8.2)	12	κFbb, κFττ, κFμμ, κFWW, κFZZ, κFγγ, κVbb, κVττ, κVμμ, κVWW, κVZZ, κVγγ
Up vs. down-type couplings	25.5% (4.06)	3	λ_Vu, λ_du, κ_uu
Lepton vs. quark couplings	27.2% (3.91)	3	λ_lq, λ_Vq, κ_qq

Constraints on benchmark two Higgs doublet models

The generic models described in Sect. 8.5 can also be interpreted in the context of explicit benchmark BSM models that contain a second Higgs doublet (2HDM) [116-118]. Only models with CP conservation and a discrete symmetry to prevent tree-level flavor changing neutral currents are considered. Under these assumptions, four 2HDM types are possible, referred to as Types I, II, III, and IV. Each of these 2HDMs contain seven free parameters. Under the additional assumption that the Higgs boson with a mass of 125.09 is the lightest CP-even, neutral Higgs boson in the extended Higgs sector, the predicted rates for its production and decay are sensitive at leading order to only two 2HDM parameters: the angles and that diagonalize the mass-squared matrices of the scalars and pseudoscalars. These two parameters are conventionally substituted by and , without loss of generality. In all of the 2HDMs, the coupling of the Higgs boson to vector bosons is modified by a factor . The 2HDM types differ in how the fermions couple to the Higgs doublets. In the Type I model, all fermions couple to just one of the Higgs doublets. In Type II, the up-type fermions couple to one of the Higgs doublets, while the down-type fermions and the right-handed leptons couple to the second. In Type III 2HDM, also referred to as “lepton-specific”, the quarks couple to one of the Higgs doublets and the right-handed leptons couple to the other. In the Type IV 2HDM, also referred to as “flipped”, the up-type fermions and right-handed leptons couple to one of the Higgs doublets, while the down-type quarks couple to the other. Table 13 shows the relation between the coupling modifiers to vector bosons, quarks and leptons and the 2HDM model parameters.

Table 13

Modifications to the couplings of the Higgs bosons to up-type () and down-type () fermions, and vector bosons (), with respect to the SM expectation, in 2HDM and for the hMSSM. The coupling modifications for the hMSSM are completed by the expressions for and , as given by Eqs. (8) and (9)

	2HDM				hMSSM
	Type I	Type II	Type III	Type IV
κ_V	sin(β-α)	sin(β-α)	sin(β-α)	sin(β-α)	s_d+s_utanβ1+tan^2β
κ_u	cos(α)/sin(β)	cos(α)/sin(β)	cos(α)/sin(β)	cos(α)/sin(β)	s_u1+tan^2βtanβ
κ_d	cos(α)/sin(β)	-sin(α)/cos(β)	cos(α)/sin(β)	-sin(α)/cos(β)	s_d1+tan^2β
κ_l	cos(α)/sin(β)	-sin(α)/cos(β)	-sin(α)/cos(β)	cos(α)/sin(β)	s_d1+tan^2β

The minimal supersymmetric standard model (MSSM) [119-122] is a specific example of a 2HDM of Type II that includes additional particle content compared to the SM. The additional strong constraints given by the nontrivial fermion-boson symmetry fix all mass relations between the Higgs bosons and the angle , at tree-level, leaving only two free parameters to fully constrain the MSSM Higgs sector, usually chosen to be and . The hMSSM scenario [123, 124], in particular, is an effective MSSM model, trading the precise knowledge of against unknown higher-order corrections such that across the , parameter space. Another requirement of the scenario is that be identified as the lightest of the two neutral scalar Higgs bosons. Furthermore, one also obtains relatively simple relations between , , and the Higgs boson coupling modifiers [125], which are shown in Table 13 and completed by Eqs. (8) and (9). Although many other MSSM benchmark models have also been defined [126], the lack of analytic expressions for the Higgs boson couplings renders these models technically more challenging to consider and they are therefore beyond the scope of this paper.

To set constraints on the 2HDM model parameters, 3-dimensional likelihood scans of the parametrizations described in Sect. 8.5 (with necessary modifications to the lepton coupling modifiers to describe the Type IV 2HDM) are performed. A test-statistic is then defined, for example in the Types I, II and hMSSM scenarios,where are the values of the POIs that maximize the likelihood. An interpolation scheme is used to determine the value of q as a function of and , or and , for the Types I and II, or hMSSM scenarios, respectively, using the relations in Table 13.

A second quantity is defined as,where is the maximum likelihood value attained in the planes of –, or –. The allowed regions are determined as the points in each plane for which the difference between q and () is less than 5.99. This value corresponds to the 95% confidence region assuming is distributed as a function with 2 degrees of freedom. A similar procedure is performed using the model with the parameters , and , to determine the allowed region for the Type III scenario.

Figure 19 shows the results of the fits for the different 2HDM benchmark scenarios. The lobe features that can be seen in the Types II, III, and IV constraints for are due to negative values of , , and , which are not excluded with the current sensitivity. In all of these 2HDM models, the Higgs boson couplings are the same as those predicted in the SM for .

Fig. 19

Constraints in the vs. plane for the Types I, II, III, and IV 2HDM, and constraints in the vs. plane for the hMSSM. The white regions, bounded by the solid black lines, in each plane represent the regions of the parameter space that are allowed at the 95% , given the data observed. The dashed lines indicate the boundaries of the allowed regions expected for the SM Higgs boson

The results for the hMSSM scenario are also shown in Fig. 19. The constraints observed are more stringent than those expected under the SM. This is due to the best fit value of being smaller than 1, while in the hMSSM for , is strictly greater than 1 and asymptotically approaches unity only at large . Therefore the observed data disfavors small values of , leading to the stronger constraint.

The constraints in the 2HDM and hMSSM scenarios are complementary to those obtained from direct searches for additional Higgs bosons [127-131].

Summary

A set of combined measurements of Higgs boson production and decay rates has been presented, along with the consequential constraints placed on its couplings to standard model (SM) particles, and on the parameter spaces of several beyond the standard model (BSM) scenarios. The combination is based on analyses targeting the gluon fusion and vector boson fusion production modes, and associated production with a vector boson or a pair of top quarks. The analyses included in the combination target Higgs boson production in the , , and decay channels, using proton–proton collision data collected in 2016 and corresponding to an integrated luminosity of 35.9. Additionally, searches for invisible Higgs boson decays are included to increase the sensitivity to potential interactions with BSM particles.

Measurements of the Higgs boson production cross section times branching fractions are presented, along with a generic parametrization in terms of ratios of production cross sections and branching fractions, which makes no assumptions about the Higgs boson total width. The combined signal yield relative to the SM prediction has been measured as at . An improvement in the measured precision of the gluon fusion production rate of around 50% is achieved compared to previous ATLAS and CMS measurements. Additionally, a set of fiducial Higgs boson cross sections, in the context of the simplified template cross section framework, is presented for the first time from a combination of six decay channels. Furthermore, interpretations are provided in the context of a leading-order coupling modifier framework, including variants for which effective couplings to the photon and gluon are introduced. All of the results presented are compatible with the SM prediction. The invisible (undetected) branching fraction of the Higgs boson is constrained to be less than 22 (38%) at 95% Confidence Level. The results are additionally interpreted in two BSM models, the minimal supersymmetric model and the generic two Higgs doublet model. The constraints placed on the parameter spaces of these models are complementary to those that can be obtained from direct searches for additional Higgs bosons.