Measurements of properties of the Higgs boson decaying to a W boson pair in pp collisions at $\sqrt{s} =$ 13 TeV

Measurements of the production of the standard model Higgs boson decaying to a W boson pair are reported. The W$^+$W$^-$ candidates are selected in events with an oppositely charged lepton pair, large missing transverse momentum, and various numbers of jets. To select Higgs bosons produced via vector boson fusion and associated production with a W or Z boson, events with two jets or three or four leptons are also selected. The event sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$, collected in pp collisions at $\sqrt{s} =$ 13 TeV by the CMS detector at the LHC during 2016. Combining all channels, the observed cross section times branching fraction is 1.28 $^{+0.18}_{-0.17}$ times the standard model prediction for the Higgs boson with a mass of 125.09 GeV. This is the first observation of the Higgs boson decay to W boson pairs by the CMS experiment.


Introduction
In the standard model (SM) of particle physics, the origin of the masses of the W and Z bosons is based on the spontaneous breaking of the electroweak symmetry. This symmetry breaking is achieved through the introduction of a complex doublet scalar field [1][2][3][4][5][6], leading to the prediction of the existence of one physical neutral scalar particle, commonly known as the Higgs boson (H). The observation of a new particle at a mass of approximately 125 GeV with Higgs boson-like properties was reported by the ATLAS [7] and CMS [8,9] Collaborations during the first running period of the CERN LHC in proton-proton (pp) collisions at center-of-mass energies of 7 and 8 TeV. Subsequent publications from both Collaborations, based on the 7 and 8 TeV data sets [10][11][12][13], established that all measured properties of the new particle, including its spin, parity, and coupling strengths to SM particles, are consistent within the uncertainties with those expected for the SM Higgs boson. A combination of the ATLAS and CMS results [14,15] further confirmed these observations and resulted in determining the boson mass to be m H = 125.09 ± 0.21 (stat) ± 0.11 (syst) GeV.
The Higgs boson decay to a pair of W bosons was studied by the ATLAS and CMS Collaborations using the 7 and 8 TeV data sets in leptonic final states, exploring several production mechanisms [16,17]. The probability of observing a signal at least as large as the one seen, under the background-only hypothesis, corresponded to a significance of 6.1 and 4.3 standard deviations (s.d.) for ATLAS and CMS respectively, while the expected significance for a SM Higgs boson was 5.8 s.d. for both Collaborations. A later CMS combination [12], that includes Higgs boson production in association with a top quark pair, reported an observed significance of 4.7 s.d. for this decay. The same decay channel was used by the ATLAS and CMS Collaborations to search for the Higgs boson off-shell production [18,19] and to perform fiducial and differential cross section measurements [20,21].
In 2015, the LHC restarted at √ s = 13 TeV, delivering high luminosity pp collisions. The new data are used to further constrain the properties of the Higgs boson: any significant deviation from the SM predictions would be a clear sign of new physics. This paper presents the analysis of the H → WW decay at 13 TeV, using a data sample corresponding to a total integrated luminosity of 35.9 fb −1 , collected during 2016.
Gluon fusion (ggH) is the dominant production mode for a Higgs boson with a mass of 125 GeV in pp collisions at √ s = 13 TeV. The large Higgs boson branching fraction to a W boson pair makes this channel one of the most suitable for a precision measurement of the Higgs boson production cross section, and also allows studies of subleading production channels, such as Higgs boson production via vector boson fusion (VBF) and associated production with a vector boson (VH). These channels are also targeted in this paper, contributing to the precision in the measurement of the Higgs boson couplings.
The leptonic decays of the two W bosons provide the cleanest decay channel, despite the presence of neutrinos in the final state that prevents the full reconstruction of the Higgs boson mass. The different-flavor (DF) leptonic decay mode eµ has the largest branching fraction, is the least affected by background processes, and therefore is the most sensitive channel of the analysis. The same-flavor (SF) e + e − and µ + µ − final states are also considered, although their sensitivity is limited by the contamination from the Drell-Yan (DY) background with missing transverse momentum due to instrumental effects.
Events with a pair of oppositely charged leptons (electrons and/or muons) and missing transverse momentum, due to the presence of neutrinos in the final state, are selected. This signature is common to other SM processes that contribute to the background in this analysis. The main contribution comes from nonresonant production of W boson pairs (WW), an irreducible background that shares the same final state and can only be separated from the signal using kinematic distributions. Backgrounds coming from top quark events (tt and tW) are also important, followed by other processes, such as W+jets and other diboson and triboson production processes. The DY process is the dominant source of background in the dielectron and dimuon final states, while it is subdominant in the electron-muon final state, where its contribution arises from the leptonic decays of the τ leptons emerging from Z/γ * → τ + τ − .
The events are categorized by jet multiplicity to better handle the tt background. In addition, dedicated categories are designed to enhance the sensitivity to the VBF and VH production mechanisms.

The CMS detector
The CMS detector is a multipurpose apparatus designed to study high transverse momentum (p T ) physics processes in proton-proton and heavy ion collisions, and is described in detail in Ref. [22] together with a definition of the coordinate system used. A superconducting solenoid occupies its central region, providing a magnetic field of 3.8 T parallel to the beam direction. Charged particle trajectories are measured by the silicon pixel and strip trackers, which cover a pseudorapidity region of |η| < 2.5. A lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter surround the tracking volume and cover |η| < 3. The steel and quartz fiber Cherenkov hadron forward calorimeter extends the coverage to |η| < 5. The muon system consists of gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid, and covers |η| < 2.4. The first level of the CMS trigger system, composed of custom hardware processors, is designed to select the most interesting events in less than 4 µs, using information from the calorimeters and muon detectors. The high-level trigger processor farm further reduces the event rate to about 1000 Hz, before data storage.

Data samples
The events used in this analysis are selected by high-level trigger algorithms that require the presence of one or two high-p T electrons or muons passing loose identification and isolation requirements. In single-lepton triggers, relatively tight lepton identification criteria are applied. The p T threshold is 25 GeV in the central region (|η| < 2.1) and 27 GeV for 2.1 < |η| < 2.5 for electrons, while it is 24 GeV for muons. In the dielectron trigger, the minimum required p T is 23 GeV for the leading and 12 GeV for the subleading electron. In the dimuon trigger, the minimum p T is 17 GeV for the leading and 8 GeV for the subleading muon. In the two dilepton eµ triggers used in the analysis, the minimum p T requirements are either 8 GeV for the muon and 23 GeV for the electron, or 23 GeV for the muon and 12 GeV for the electron. The combination of single-lepton and dilepton triggers provides an overall trigger efficiency in excess of 98% for signal events.

Event simulation
Several event generators are used to optimize the analysis and estimate the expected yields of signal and backgrounds, as well as their associated systematic uncertainties. Different Higgs boson production mechanisms are simulated. Both ggH and VBF are generated with POWHEG v2 [23][24][25][26], which describes the full next-to-leading order (NLO) perturbative quantum chromodynamics (QCD) properties of these processes. In addition, the ggH process is reweighted to match the calculation of POWHEG NNLOPS [27], which provides a next-to-next-to-leading order (NNLO) description for the inclusive Higgs boson production, NLO for the exclusive H + 1 jet production, and the leading order (LO) for the exclusive H + 2 jets production. The reweighting is performed by computing the ratio of the Higgs boson p T distribution from the NNLOPS generator to that from the POWHEG generator in each jet multiplicity bin, and applying this ratio to the ggH POWHEG simulation. The MINLO HVJ [28] extension of POWHEG is used to simulate the associated production of the Higgs boson with vector bosons (W + H, W − H, ZH), which simulates the VH + 0 and 1 jet processes with NLO accuracy. Higgs boson production in association with top or bottom quarks, such as ttH and bbH production mechanisms, are considered as well, although they only contribute to a minor extent in the phase space selected by this analysis. For the simulation of ttH production the POWHEG generator is used, while the MADGRAPH5 aMC@NLO v2.2.2 generator [29] is used to simulate the bbH production. The Higgs boson is generated with a mass of 125.09 GeV and is let to decay into a pair of W bosons, considering only leptonic W boson decays (e, µ, or τ). For Higgs bosons produced via ggH [30] and VBF [31] processes, their decay into two W bosons and subsequently into leptons is simulated using JHUGEN v5.2.5 [32,33]. For the associated production mechanisms, including gluon fusion produced ZH, the Higgs boson decay is simulated by PYTHIA 8.212 [34]. The simulated signal samples are normalized using cross sections [35] and decay rates [36] computed by the LHC Higgs Cross Section Working Group. In particular the most recent next-to-next-to-next-toleading order calculations for the inclusive gluon fusion production are used [35]. Additional simulated samples, where the Higgs boson decays into a pair of τ leptons, are also produced for each of the aforementioned production mechanisms. The H → ττ events are considered as signal events, although the expected contribution in the signal phase space is small compared to H → W + W − .
The various background processes in this study are simulated as follows: POWHEG v2 [37] is used for qq → WW production, whereas gg → WW production is generated using MCFM v7.0 [38]. A WW simulation with two additional jets is generated with MADGRAPH5 aMC@NLO at LO accuracy via diagrams with six electroweak (EW) vertices, referred to as WW EW production. In order to suppress the top quark background processes, the analysis is performed defining event categories with different number of high-p T jets (p T > 30 GeV). The classification of the events in bins of jet multiplicity spoils the convergence of fixed-order calculations of the qq → WW process and requires the use of dedicated resummation techniques for an accurate prediction of the differential distributions [39,40]. The simulated qq → WW events are therefore reweighted to reproduce the p T W W distribution from the p T -resummed calculation.
The LO cross section for the gg → WW process is obtained directly from MCFM. For this process, the difference between LO and NLO cross sections is significant; a K factor of 1.4 is calculated [41] and applied to the gg → WW simulation. Given the low sensitivity of the analysis to the invariant mass of the WW system, m WW , and the theoretical uncertainties in the K factor, an m WW -independent calculation is used.
Single top quark and tt processes are generated using POWHEG v2. The cross sections of the different single top quark processes are estimated at NLO accuracy [42], while the tt cross section is computed at NNLO accuracy, with next-to-next-to-leading-logarithmic soft-gluon resummation [43].
The DY production of Z/γ * is generated using MADGRAPH5 aMC@NLO at NLO accuracy using the FxFx jet matching and merging scheme with a merging scale µ Q = 30 GeV [44], and the Z/γ * p T distribution has been reweighted to match the distribution observed in data in dimuon events. Other multiboson processes, such as WZ, ZZ, and VVV (V = W, Z), are also simulated with MADGRAPH5 aMC@NLO at NLO accuracy.
All processes are generated using the NNPDF 3.0 [45,46] parton distribution functions (PDFs), with the accuracy matching that of the matrix element calculations. All the event generators are interfaced to PYTHIA for the showering of partons and hadronization, as well as the simulation of the underlying event (UE) and multiple-parton interactions based on the CUET8PM1 tune [47].
To estimate the systematic uncertainties related to the choice of the UE and multiple-parton interactions tune, the signal processes and the WW background are also generated with alternative tunes, which are representative of the uncertainties in the CUET8PM1 tuning parameters. The systematic uncertainty associated with showering and hadronization is estimated by interfacing the same samples with the HERWIG++ 2.7 generator [48,49], using the UE-EE-5C tune for the simulation of UE and multiple-parton interactions [47].
For all processes, the detector response is simulated using a detailed description of the CMS detector, based on the GEANT4 package [50]. Additional simulated minimum bias pp interactions from PYTHIA are overlapped with the event of interest in each collision to reproduce the number of interactions per bunch crossing (pileup) measured in data. The average number of pileup interactions is about 27 per event for the 2016 data set used in this analysis.

Event reconstruction
A particle-flow (PF) algorithm [51] is used to reconstruct the observable particles in the event.
Energy deposits (clusters) measured by the calorimeters, charged particle tracks identified in the central tracking system and the muon detectors are combined to reconstruct individual particles.
Among the vertices reconstructed in the event, the one with the largest value of summed physics-object p 2 T is taken to be the primary pp interaction vertex. The physics objects are the jets, clustered using the jet finding algorithm [52,53] with the tracks assigned to the vertex as inputs, and the associated missing transverse momentum, taken as the negative vector sum of the p T of those jets.
Electrons are reconstructed by matching clusters in the ECAL to tracks in the silicon tracker [54]. In this analysis, electron candidates are required to have |η| < 2.5. Additional requirements are applied to reject electrons originating from photon conversions in the tracker material or jets misreconstructed as electrons. Electron identification criteria rely on observables sensitive to the bremsstrahlung along the electron trajectory, the geometrical and momentum-energy matching between the electron track and the associated energy cluster in the ECAL, as well as ECAL shower shape observables and association with the primary vertex.
Muon candidates are reconstructed in the geometrical acceptance |η| < 2.4 by combining information from the silicon tracker and the muon system. Identification criteria based on the number of measurements in the tracker and in the muon system, the fit quality of the muon track, and its consistency with its origin from the primary vertex are imposed on the muon candidates to reduce the misidentification rate.
Prompt leptons coming from EW interactions are usually isolated, whereas misidentified lep-tons and leptons coming from jets are often accompanied by charged or neutral particles, and can arise from a secondary vertex. Hence charged leptons are required to satisfy the isolation criterion that the p T sum over charged PF candidates associated with the primary vertex, exclusive of the lepton itself, and neutral PF particles in a cone of a radius ∆R = where φ is the azimuthal angle in radians, centered on the muon (electron) direction is below a threshold of 15 (6)% relative to the muon (electron) p T . To mitigate the effect of the pileup on this isolation variable, a correction based on the average energy density in the event [55] is applied. Additional requirements on the transverse (|d xy |) and longitudinal (|d z |) impact parameters with respect to the primary vertex are included. Electrons detected by the ECAL barrel are required to have |d z | < 0.10 cm and |d xy | < 0.05 cm, while electrons in the ECAL endcap must satisfy |d z | < 0.20 cm and |d xy | < 0.10 cm. For muons, the |d z | parameter is required to be less than 0.10 cm, while |d xy | is required to be less than 0.01 cm for muons with p T < 20 GeV and less than 0.02 cm for p T > 20 GeV.
The jet reconstruction starts from all PF candidates, and removes the charged ones that are not associated with the primary vertex to mitigate the pileup impact. The remaining charged PF candidates and all neutral candidates are clustered by the anti-k T algorithm [52] with a distance parameter of 0.4. To reduce further the residual pileup contamination from neutral PF candidates, a correction based on the jet area [55] is applied. The jet energy is calibrated using both simulation and data following the technique described in Ref. [56]. To identify jets coming from b quarks (b jets), a multivariate (MVA) b tagging algorithm is used [57]. In this analysis, the chosen working point corresponds to about 80% efficiency for genuine b jets, and to a mistagging rate of about 10% for light-quark or gluon jets and of 35 to 50% for c jets. A per-jet scale factor is computed and applied to account for b tagging efficiency and mistagging rate differences between data and simulation.
The missing transverse momentum vector ( p miss T ), whose magnitude is denoted as p miss T , is reconstructed as the negative vectorial sum in the transverse plane of all PF particle candidate momenta. Since the presence of pileup induces a degradation of the p miss T measurement, affecting mostly backgrounds with no genuine p miss T , such as DY production, another p miss T that is constructed from only the charged particles (track p miss T ) is used in events with an SF lepton pair (ee or µµ). To suppress the remaining off-peak DY contribution in categories containing events with an SF lepton pair, a dedicated MVA selection based on a boosted decision tree algorithm (BDT) is used, combining variables related to lepton kinematics and p miss T . The BDT is trained on simulated samples separately for different jet multiplicity categories, and the output discriminator is used to define a phase space enriched in signal events and reduced DY background contamination.

Event preselection and categorization
Events are required to pass the single-lepton or dilepton triggers. For each event, this analysis requires at least two high-p T lepton candidates with opposite sign, originating from the primary vertex, categorized as dielectron, dimuon, or eµ pairs. Only jets with p T > 30 GeV (20 GeV for b jets) and |η| < 4.7 (|η| < 2.4 for b jets) are considered in the analysis. Jets are ignored if they overlap with an isolated lepton within a distance of ∆R = 0.3. In addition, the following kinematic selection is applied in the eµ final state: one electron and one muon are required to be reconstructed in the event with a minimum p T of 13 GeV for the electron and 10 GeV for the muon, the higher p T threshold for the electron resulting from the trigger definition. One of the two leptons should also have a p T greater than 25 GeV. In the case of SF e + e − and µ + µ − final states, the leading lepton is required to have p T greater than 25 GeV when it is an electron, or 20 GeV when it is a muon. The subleading electron is required to have p T greater than 13 GeV, while for the muon a minimum p T of 10 GeV is required. Both leptons are required to be well identified, isolated, and prompt.
Given the large background contribution from tt production in both DF and SF final states, events are further categorized based on the number of jets in the event, with the 0-jet category driving the sensitivity of the analysis. A categorization of the selected events is performed, targeting different production mechanisms and different flavor compositions of the WW decay products.

Different-flavor ggH categories
The categories described in this section target the ggH production mechanism and select the DF eµ final state. The main background processes are the nonresonant WW, top quark (both single and pair production), DY to τ lepton pairs, and W+jets when a jet is misidentified as a lepton. Smaller background contributions come from WZ, ZZ, Vγ, Vγ * , and triboson production. The WW background process can be distinguished from the signal by the different kinematic properties of the lepton system, since it is dominated by the on-shell W boson pairs that do not arise from a scalar resonance decay. The top quark background process is diluted by defining different categories that depend on the number of jets in the event, and reduced by vetoing any b-tagged jet with p T > 20 GeV.
The W+jets contribution (also referred to as nonprompt lepton background), where one jet mimics the signature of an isolated prompt lepton, is an important background process especially in the 0-and 1-jet ggH-tagged DF categories. This background is reduced by taking advantage of the charge symmetry of the signal, and the charge asymmetry of the W+jets process, in which the production of W + is favored over W − . Also, the fact that the probabilities for a jet to mimic an electron or a muon are different, and the fact that the misidentification rate is larger for lower-p T leptons, are exploited. Following these physics motivations the 0-and 1-jet ggH-tagged DF categories are further split into four categories according to the lepton flavor and charge: e + µ − , e − µ + , µ + e − , and µ − e + , where the first lepton is the one with the higher p T . In addition, the four categories are divided according to whether the subleading lepton p T (p T2 ) is above or below 20 GeV. This eight-fold partitioning of the 0-and 1-jet ggH-tagged categories provides an improvement in terms of the expected significance of about 15% with respect to the merged 0-and 1-jet categories.
To suppress background processes with three or more leptons in the final state, no additional identified and isolated leptons with p T > 10 GeV are allowed in the events for the dilepton categories. The dilepton invariant mass (m ) is required to be higher than 12 GeV, to reject low-mass resonances and background that comes from events with multiple jets that all arise through the strong interaction (referred to the multijet background). To suppress the background arising from DY events decaying to a τ lepton pair, which subsequently decays to the eµ final state, and to suppress processes without genuine missing transverse momentum, a minimum p miss T of 20 GeV is required. In the two lepton categories, the DY background is further reduced by requiring the dilepton p T (p T ) to be higher than 30 GeV, as on average eµ lepton pairs from Z → τ + τ − decays have lower p T than the ones from H → WW decays. These selection criteria also reduce contributions from H → WW → τντν and H → τ + τ − . Finally, to further suppress contributions from Z → τ + τ − and W+jets events, where the subleading lepton does not arise from a W boson decay, the transverse mass built with p miss T and the subleading lepton, defined as: is required to be greater than 30 GeV. Here ∆φ( 2, p miss T ) is the azimuthal angle between the subleading lepton momentum and p miss T .
Although the invariant mass of the Higgs boson cannot be reconstructed because of the undetected neutrinos, the expected kinematic properties of the Higgs boson production and decay can be exploited. The spin-0 nature of the SM Higgs boson results in the preferential emission of the two charged leptons in the same hemisphere. Moreover, the invariant mass of the two leptons in the signal is relatively small with respect to the one expected for a lepton pair arising from other processes, such as nonresonant WW and top quark production. On the other hand, several of the smaller remaining background processes, such as nonprompt leptons, DY→ τ + τ − , and Vγ populate the same m phase space as the Higgs boson signal. These can be partially disentangled from the signal by reconstructing the Higgs boson transverse mass as: where ∆φ( , p miss T ) is the azimuthal angle between the dilepton momentum and p miss T . These additional background processes populate different regions of the two-dimensional plane in m and m T . A shape analysis based on a two-dimensional binned template fit of m versus m T is performed to extract the Higgs boson signal in the DF ggH categories.
The observed events as a function of m and m T are shown in Figs. 1, 2, and 3, after the template fit to the (m , m T ) distribution. The 0-and 1-jet categories are split into p T2 < 20 GeV and p T2 > 20 GeV subcategories, to show the different purity of the two regions. In these figures the postfit number of events is shown, i.e. each signal and background process is normalized to the result of a simultaneous fit to all categories, assuming that the relative proportions for the different Higgs boson production mechanisms are those predicted by the SM. The events in each bin of one of the two variables are obtained by integrating over the other, and weighted using the ratio of fitted signal (S) to the sum of signal and background (S + B). As an example, this means that to obtain the weighted m distribution, this variable is integrated to compute the S/(S + B) ratio in each m T bin. This ratio is then used to perform a weighted sum of the m distributions in each m T bin. A similar weighting procedure is applied when merging the distributions of a given variable in different categories. The weighting procedure is used only for visualization purposes, and is not used for signal extraction.
The full list of DF ggH categories and their selection requirements is shown in Table 1.

Different-flavor VBF category
The VBF process is the second largest Higgs boson production mechanism at the LHC. This mode involves the production of a Higgs boson in association with two jets with large rapidity separations. After the common preselection, the VBF analysis requires events with exactly two jets with p T > 30 GeV, a pseudorapidity separation (|∆η jj |) between the two jets larger than 3.5, and an invariant mass (m jj ) greater than 400 GeV. The VBF analysis is based on the shape of the m distribution, and is split into two signal regions, one with 400 < m jj < 700 GeV and the other with m jj > 700 GeV, to profit from the higher purity of the m jj > 700 GeV region. The observed signal and background events as functions of m are shown in Fig. 4, for the two m jj regions separately. The list of event requirements applied in this category is presented in Table 2.     Preselection -m > 12 GeV, p T1 > 25 GeV, p T2 > 13 (10) GeV for e (µ), p miss T > 20 GeV, p T > 30 GeV no additional leptons with p T > 10 GeV electron and muon with opposite charges 0-jet ggH-tagged no jets with p T > 30 GeV no b-tagged jets with p T between 20 and 30 GeV no jets with p T > 30 GeV no b-tagged jets with p T between 20 and 30 GeV

Different-flavor VH with two jets category
The VH process involves the production of a Higgs boson in association with a W or Z boson. The 2-jet VH-tagged category targets final states where one vector boson (W or Z) decays into two resolved jets. This category with hadronically decaying vector bosons is affected by large backgrounds compared to the leptonic decays, but profits from a higher branching fraction. The 2-jet VH-tagged analysis reverses the pseudorapidity separation requirement of the VBF selection (|∆η| < 3.5) and requires m jj to be between 65 and 105 GeV. In addition, the two leading jets are required to be central (|η| < 2.5) to profit from more stringent b jet veto requirements, given that b tagging can only be performed for central jets.
The analysis is based on the shape of the m discriminant distribution, presented in Fig. 5. The list of event requirements applied is presented in Table 3.

Same-flavor ggH categories
Similarly to the DF ggH-tagged analysis described in Section 5.1, an analysis targeting ggH in the SF e + e − and µ + µ − channels is performed. The main challenge in this final state is the large DY background contribution. In order to control it, a BDT is trained to build a discriminator, called DYMVA, to identify DY events. Table 2: Analysis categorization and event requirements for the 2-jet VBF-tagged category, in the DF dilepton final state. The phase spaces defined by the 2-jet VBF-tagged requirements correspond to the events shown in Fig. 4.

Category
Subcategory Requirements Preselection -m > 12 GeV, p T1 > 25 GeV, p T2 > 13 (10) GeV for e(µ), p miss T > 20 GeV, p T > 30 GeV no additional leptons with p T > 10 GeV electron and muon with opposite charges 2-jet VBF-tagged eµ low m jj exactly two jets with p T > 30 GeV 60 < m T < 125 GeV leptons η between the two leading jets 400 < m jj < 700 GeV and |∆η jj | > 3.5 no b-tagged jets with p T > 20 GeV eµ high m jj exactly two jets with p T > 30 GeV 60 < m T < 125 GeV leptons η between the two leading jets m jj > 700 GeV and |∆η jj | > 3.5 no b-tagged jets with p T > 20 GeV Table 3: Analysis categorization and event requirements for the 2-jet VH-tagged category, in the DF dilepton final state. The phase space defined by the 2-jet VH-tagged requirements corresponds to the events shown in Fig. 5.

Category
Subcategory Requirements Preselection -m > 12 GeV, p T1 > 25 GeV, p T2 > 13 (10) GeV for e (µ) p miss T > 20 GeV, p T > 30 GeV no additional leptons with p T > 10 GeV electron and muon with opposite charges 2-jet VH-tagged eµ at least two jets with p T > 30 GeV two leading jets with |η| < 2.5 60 < m T < 125 GeV and ∆R < 2 no b-tagged jets with p T > 20 GeV 65 < m jj < 105 GeV and |∆η jj | < 3.5 A categorization based on the p T of the subleading lepton is introduced to better control the nonprompt lepton background, and a categorization in the number of jets is used to control the top quark backgrounds. The full list of event requirements is shown in Table 4.
This is an event counting analysis. The DY background estimations in these channels are based exclusively on control samples in data, as described in Section 6.

Associated WH production with three leptons in the final state
The three-lepton WH-tagged analysis selects events that have the leading lepton with p T1 > 25 GeV, the subleading with p T2 > 20 GeV, and the trailing with p T3 > 15 GeV. Events with a fourth lepton with p T > 10 GeV are discarded. A veto is applied to events with SF lepton pairs of opposite charge that are compatible with coming from the decay of a Z boson. Events containing jets with p T > 30 GeV or b-tagged jets with p T > 20 GeV are also vetoed. The azimuthal angle between p miss T and the three-lepton system p T , ∆φ( , p miss T ), is used to reduce the contamination of nonprompt lepton backgrounds. The rest of the three-lepton WH-tagged selection is in common with the other categories. These requirements are summarized in Table 5.
The events are further divided into two categories: same-sign SF (SSSF) lepton pairs, µ ± µ ± e ∓ /e ± e ± µ ∓ , and opposite-sign SF (OSSF) lepton pairs, µ ∓ µ ± e ∓ /e ∓ e ± µ ∓ . The two selections have different signal-over-background ratios, with the SSSF being the purest of the two. The main background contribution in both cases is the contamination from nonprompt leptons. In the OSSF category, events are required to have p miss T > 50 GeV to reduce the DY background.
The analysis is based on the minimum ∆R between oppositely charged leptons. The distribution of this variable is presented in Fig. 6, separately for the SSSF and OSSF categories. > 20 GeV, p T > 30 GeV no additional leptons with p T > 10 GeV two electrons or two muons with opposite charges 0-jet ggH-tagged DYMVA > 0.991, m < 55 GeV, m T > 50 GeV, e + e − p T2 < 20 GeV p T2 < 20 GeV, ∆φ < 1.7 µ + µ − p T2 < 20 GeV no jets with p T > 30 GeV no b-tagged jets with p T > 20 GeV exactly one jet with p T > 30 GeV no b-tagged jets with p T > 20 GeV Table 5: Analysis categorization and event requirements for the WH-tagged category, in the three-lepton final state. Here, min-m + − is the minimum m between the oppositely charged leptons. For the Z boson veto, the pair with the m closest to the Z boson mass is considered. Events that fulfill the three-lepton WH-tagged requirements correspond to the signal phase space shown in Fig. 6.

Associated ZH production with four leptons in the final state
The ZH final state is targeted by requiring exactly four isolated leptons with tight identification criteria and zero total charge, and large p miss T from the undetected neutrinos. The major background processes are ZZ and ttZ production.
Among the four leptons, the pair of SF leptons with opposite charge, and with the invariant mass closest to the Z boson mass, is chosen as the Z boson candidate. The remaining dilepton system, denoted as X, can be either SF or DF. Events are therefore divided in two categories, distinguishing between the cases in which the X candidate contains two DF leptons (XDF) or two SF leptons (XSF), as shown in Table 6.
The signal fraction is equally distributed in the two regions. In the XSF region, ZZ, DY, and ttZ production are the major background sources, while in the XDF region, ttZ and ZZ backgrounds are dominant. Backgrounds with two Z bosons fall predominantly into the XSF region, and enter the XDF selection only through the leptonic decays of the τ leptons. This makes the XDF region much cleaner than the XSF one.
Given the low expected signal yield in the XDF and XSF categories, the result in this case is extracted from event counting. 6 Background estimation 6

.1 Nonprompt lepton background
Events in which a single W boson is produced in association with jets may populate the signal region when a jet is misidentified as a lepton. These events contain a genuine lepton and p miss T from the W boson decay as well as a second nonprompt lepton from a misidentified jet, likely arising from a B hadron decay. A similar background arises from semileptonic decays of top quark pairs, especially in the 1-and 2-jets categories. At a lower rate, multijet production Table 6: Analysis categorization and event requirements for the ZH-tagged category, in the four-lepton final state. Here, X is defined as the remaining lepton pair after the Z boson candidate is chosen. The component leptons of X can be either same-flavor (XSF) or different-flavor (XDF).

Category Subcategory Requirements
Preselection -four tight and isolated leptons, with zero total charge p T > 25 GeV for the leading lepton p T > 15 GeV for the second leading lepton p T > 10 GeV for the remaining two leptons no additional leptons with p T > 10 GeV Z dilepton mass > 4 GeV X dilepton mass > 4 GeV 4-lepton ZH-tagged XSF |m − m Z | < 15 GeV 10 < m X < 50 GeV 35 < p miss T < 100 GeV four-lepton invariant mass > 140 GeV XDF |m − m Z | < 15 GeV 10 < m X < 70 GeV p miss T > 20 GeV no b-tagged jets with p T > 20 GeV and fully hadronic top quark pair decays also contribute. These backgrounds are particularly important for events with low-p T leptons and low m , and hence in the signal region of the analysis.
The nonprompt lepton background is suppressed by the identification and isolation requirements imposed on the electrons and muons, while the remaining contribution is estimated directly from data. A control sample is defined using events in which one lepton passes the standard lepton identification and isolation criteria and another lepton candidate fails the criteria but passes a looser selection, resulting in a sample of "pass-fail" lepton pairs. The efficiency ( misID ) for a jet that satisfies this looser selection to pass the standard selection is estimated directly from data in an independent sample dominated by events with nonprompt leptons from multijet processes. The efficiency misID is parameterized as a function of the p T and η of the leptons, and is used to weight the events in the pass-fail sample by misID /(1 − misID ), to obtain the estimated contribution from this background in the signal region. The systematic uncertainties associated with the determination of misID dominate the overall uncertainty of this method. The systematic uncertainty has two sources: the dependence of misID on the sample composition and the uncertainty in the method. The first source is estimated by modifying the jet p T threshold in the QCD multijet sample, which modifies the jet sample composition. The uncertainty in the method is obtained from a closure test, where misID is derived from simulated QCD multijet events and applied to a W+jets simulated sample to predict the number of background events. The total uncertainty in misID , including the statistical precision of the control sample, is about 40%.

Top quark background
Background contamination from single top quark processes, in particular tW associated production, and from tt production, arises because of the inefficiency of b jet identification and the relatively large top quark cross sections at 13 TeV. The shapes of the top quark background distributions in the various categories are obtained from simulation, taking into account the measured b jet identification inefficiencies. The normalizations are obtained from control regions enriched in top quark events. The background estimation is obtained separately for the 0-, 1-and 2-jet ggH-tagged categories, the 2-jet VBF-and VH-tagged categories, and for DF and SF final states.
The control region for the 0-jet ggH-tagged category is defined the same way as the signal region, except for the requirement that at least one jet with 20 < p T < 30 GeV is identified as a b jet by means of the b tagging algorithm. For the 1-jet ggH-tagged top quark enriched region, exactly one jet with p T > 30 GeV identified as a b jet is required. In the 2-jet top quark enriched regions (either ggH-, VH-, or VBF-tagged), two jets with p T > 30 GeV must be present in the event and at least one has to be identified as a b jet. To reduce other backgrounds in the top quark control regions, the dilepton mass is required to be higher than 50 GeV. The derived scale factors are shown in Table 7. The normalization of the top quark background in the three-and four-lepton categories is instead taken from simulation with its NNLO cross section uncertainty. The top quark p T in tt events is reweighted in simulated samples in order to have a better description of the p T distribution observed in data, as described in previous CMS analyses [58]. The difference between applying this reweighting, or not, is taken as a systematic shape uncertainty. The theoretical uncertainty related to the single top quark and tt cross sections is also taken into account. It is evaluated by varying the ratio between the single top quark and tt cross section by its uncertainty, which is 8% at 13 TeV [17]. All the experimental uncertainties described in Section 7 are also included as uncertainties on the top quark background shape.

Drell-Yan background
The DY → τ + τ − background is relevant for DF categories and, like the signal, populates the low-m T and low-m phase space. The kinematic variables of this background are predicted by the simulation after reweighting the Z boson p T spectrum to match the distribution measured in the data. The normalization is estimated in data control regions by selecting events with m T < 60 GeV and 30 < m < 80 GeV. Normalization scale factors are extracted, separately for the 0-, 1-, 2-jet ggH-tagged, the 2-jet VBF-and VH-tagged categories, and are shown in Table 8.
The effect of missing higher-order corrections in the DY simulation is estimated by varying the renormalization and factorization scales by a factor of two up and down. This effect is treated as a shape uncertainty and amounts to 1-2% in the DY yield. All experimental uncertainties described in Section 7 are considered as shape uncertainties for this background process.
In the SF categories, a dominant source of background is DY → e + e − and DY → µ + µ − . The contribution of the DY background is estimated from Z/γ * → events outside the Z boson mass region in data (dubbed as the out region, which corresponds to the signal region of the analysis) by counting the number of events in the Z boson mass region in data (in region), subtracting the non-Z-boson contribution from it, and scaling the yield by a ratio R out/in . This ratio is defined as the fraction of events outside and inside the Z boson mass region in Monte Carlo (MC) simulation, R out/in = N MC out /N MC in . The Z boson mass region is defined as |m − m Z | < 7.5 GeV. Such a tight mass window is chosen to reduce the non-Z-boson background contributions, which can be split into two categories. The first one is composed of the background processes, such as top quark pair and W + W − production, with equal decay rates into the four lepton-flavor final states (ee, eµ, µe, and µµ). Their contributions to the Z boson mass region in data, N background|in , can be estimated from the number of events in the e ± µ ∓ final state, N in eµ , applying a correction factor that accounts for the differences in the detection efficiency between electrons and muons (k ee and k µµ ): where stands for ee or µµ. The factor of 1/2 comes from the relative branching fraction between the and eµ final states. The second category, neglected in Eq. 3, is composed of background processes, such as WZ and ZZ (denoted as VV) production, with subsequent decay mostly into SF final states via the on-shell Z boson, which are determined from simulation. The number of events arising from these background processes are denoted as N in eµ (VV) and N in (VV), respectively for DF and SF final states.
Finally, the number of DY events in the signal region is estimated from the number of events in the SF final state, N in , separately for electrons and muons according to the following formula: The R out/in ratio is cross-checked in data after subtracting non-DY contributions; these are estimated from DF data for most of the backgrounds, but the VV background is taken from simulation. The difference of the R out/in values from the data and simulation is taken as a systematic uncertainty, and amounts to 10-25%.

The WZ and Wγ * background
The Wγ * EW production is included in the simulation as part of the WZ production, and the two processes are separated using a 4 GeV threshold on the Z/γ * mass at the generator level. For the final states with two leptons, the WZ and Wγ * processes may contribute to the signal region whenever one of the three leptons is not identified. Therefore, it is important to observe the process in data to validate the simulation.
The yield of the WZ background is measured in data by selecting events with three isolated leptons, two electrons and one muon (eeµ), or two muons and one electron (µµe). The SF lepton pair is identified as the Z boson candidate, and its invariant mass is required to be within the Z boson mass window defined in Section 6.3. This phase space is used to derive a scale factor for the WZ simulation, which is found to be 1.14 ± 0.18, with the uncertainty coming from the event count in the control region and from the difference between the eeµ and µµe phase spaces.
A Wγ * -enriched control region is defined by selecting events with two muons with invariant mass below 4 GeV, likely arising from a γ * decay, and a third isolated electron or muon passing a tight identification requirement. The dimuon invariant mass region close to the J/ψ resonance mass is discarded. This control region is used to derive a scale factor for the Wγ * simulation, which is found to be 0.9 ± 0.2, with the uncertainty coming from the event count in the µµe and µµµ samples.
All experimental uncertainties described in Section 7 are considered as shape and yield uncertainties for the WZ and Wγ * background determination.

Nonresonant WW and other backgrounds
The nonresonant WW background populates the entire two-dimensional phase space in m and m T , while the Higgs boson signal is concentrated at low m values, and m T values around the Higgs boson mass. The yield of this background is hence estimated directly from the fit procedure, separately for each category. The derived scale factors are shown in Table 9. In the qq → WW process, the p WW T spectrum in simulation is reweighted to match the resummed calculation [39,40]. The associated shape uncertainties related to missing higher orders are modeled by independently varying the factorization, renormalization, and resummation scales by a factor of two up and down. The cross section of the gluon-induced WW process is scaled to NLO accuracy and the uncertainty on this K factor is 15% [59]. In categories with at least two jets, the EW WW production is also taken into account. The theoretical uncertainty in the LO cross section of this process amounts to 11%, and is estimated by varying the renormalization and factorization scales by a factor of two up and down, including also the effect of PDF variations.
The WZ and Zγ * backgrounds in the three-lepton WH-tagged analysis are estimated using dedicated control regions from which the scale factors of 1.09 ± 0.06 and 1.61 ± 0.18, respectively, are derived. The ZZ background in the four-lepton ZH-tagged analysis is also estimated using a control region from which a scale factor of 0.96 ± 0.07 is derived.
All remaining backgrounds from diboson and triboson production are estimated according to their expected theoretical cross sections and the shape is taken from simulation.

Statistical procedure and systematic uncertainties
The statistical methodology used to interpret subsets of data selected for the H → WW analysis and to combine the results from the independent categories has been developed by the ATLAS and CMS Collaborations in the context of the LHC Higgs Combination Group. A general description of the methodology can be found in Ref. [60].
The number of events in each category and in each bin of the discriminant distributions used to extract the signal is modeled as a Poisson random variable, with a mean value that is the sum of the contributions from the processes under consideration. Systematic uncertainties are represented by individual nuisance parameters with log-normal distributions. The uncertainties affect the overall normalizations of the signal and backgrounds, as well as the shapes of the predictions across the distributions of the observables. Correlations between systematic uncertainties in different categories are taken into account.
The various control regions described in Section 6 used to constrain individual backgrounds are included in the fit in the form of single bins, representing the number of events in each of the control regions. Their normalizations are left unconstrained in the fit.
The remaining sources of systematic uncertainties of experimental and theoretical nature are described below. Effects due to the experimental uncertainties are estimated by scaling or smearing the targeted variable in the simulation and recalculating the analysis results. All experimental sources of systematic uncertainty, except for the integrated luminosity, have both a normalization and a shape component. The following experimental uncertainties are considered: • The uncertainty in the measured luminosity, which is 2.5% [61].
• The trigger efficiency uncertainty associated with the combination of single-lepton and dilepton triggers, which is 2% [62].
• The uncertainties in the lepton reconstruction and identification efficiencies, which vary within 2-5% for electrons [54] and 1-2% for muons [63], depending on p T and η. • The muon momentum and electron energy scale and resolution uncertainties, which amount to 0.6-1.0% for electrons and 0.2% for muons.
• The jet energy scale uncertainties, which vary in the range 1-13%, depending on the p T and η of the jet [64].
• The p miss T resolution uncertainty, which is taken into account by propagating lepton and jet energy scale and resolution uncertainties to p miss T . • The scale factors correcting the b tagging efficiency and mistagging rates, which are varied within their uncertainties. This systematic uncertainty, which varies between 0.5-1.0% [57], affects, in an anticorrelated way, the top quark control regions and the signal ones.
The uncertainties in the signal and background production rates due to the limited knowledge of the processes under study include several components, which are assumed to be independent: the choices of PDFs and the strong coupling constant α S , the UE and parton shower model, and the effects of missing higher-order corrections via variations of the renormalization and factorization scales. As most of the backgrounds are estimated from control regions in data, these theoretical uncertainties mostly affect the Higgs boson signal and they are implemented as normalization-only uncertainties unless stated otherwise.
The PDFs and α S uncertainties are further split between the cross section normalization uncertainties computed by the LHC Higgs Cross Section Working Group [36] for the Higgs boson signal and their effect on the acceptance [65]. The signal cross section normalization uncertainties amount to 3% for the ggH and 2% for the VBF Higgs boson production mechanism, between 1.6% and 1.9% for VH processes, and 3.6% for ttH production. The acceptance uncertainties are less than 1% for all production mechanisms.
The effect of missing higher order QCD corrections on the ggH production mechanism is split into several components that are propagated in a correlated way across the various jet multiplicity categories. Nine individual components have been identified [35] and the overall effect is about 10%. The effect of missing higher-order corrections in the VBF and VH simulations is less than 1%, while it amounts to about 8% for the ttH simulation.
The UE uncertainty is estimated by comparing two different MC simulation tunes, while the dependence on the parton shower is estimated by comparing samples interfaced with different parton shower models, as described in Section 3. The effect on the expected ggH signal yields is about 5% for the UE tuning and about 7% for the parton shower description.
Finally, the uncertainties arising from the limited number of events in the simulated samples are included independently for each bin of the discriminant distributions in each category.

Results
The signal strength modifier (µ), defined as the ratio between the measured signal cross section and the SM expectation in the H → WW → 2 2ν decay channel, is measured by performing a binned maximum likelihood fit using simulated binned templates for signal and background processes.
The combined results obtained using all the individual analysis categories are described in this section. A summary of the expected fraction of different signal production modes in each category is shown in Fig. 7, together with the total number of expected H → WW events. The chosen categorization proves effective in tackling the different production mechanisms, especially ggH, VBF, and VH. The measurements assume a Higgs boson mass of m H = 125.09 GeV, as reported in the ATLAS and CMS combined Higgs boson mass measurement [14]. The results reported below show a very weak dependence on the Higgs boson mass hypothesis.
The number of expected signal and background events, and the number of observed events in data, in each category after the full event selection are shown in Tables 10 and 11. Postfit event yields are also shown in parentheses, and correspond to the result of a simultaneous fit to all categories, assuming that the relative proportions of the different production mechanisms are those predicted by the SM.

Signal strength modifiers
The signal strength modifier is extracted by performing a simultaneous fit to all categories assuming that the relative proportions of the different production mechanisms are as the SM ones. As such, the value of µ provides an insight into the compatibility between this measurement and the SM. The combined observed signal strength modifier is: µ = 1.28 +0.18 −0.17 = 1.28 ± 0.10 (stat) ± 0.11 (syst) +0.10 −0.07 (theo), where the statistical, systematic, and theoretical uncertainties are reported separately. The statistical component is estimated by fixing all the nuisance parameters to their best fit values and   recomputing the likelihood profile. The breakdown of a given group of uncertainties (systematic or theoretical) is obtained by fixing all the nuisance parameters in the group to their best fit values, and recomputing the likelihood profile. The corresponding uncertainty is then taken as the difference in quadrature between the total uncertainty and the one obtained fixing the group of nuisance parameters. The expected and observed likelihood profiles as functions of the signal strength modifier are shown in Fig. 8, with the 68% and 95% confidence level (CL) indicated. The observed significance in the asymptotic approximation [66] of the Higgs boson production for the combination of all categories is 9.1 s.d., to be compared with the expected value of 7.1 s.d.. As such, this is the first observation of the Higgs boson decay to W boson pairs with the CMS experiment.
In order to assess the compatibility of the observed signal with the SM predictions in each category of the analysis and to ascertain the compatibility between the different categories, a simultaneous fit in which the signal strength modifier is allowed to float independently in each category is performed. The observed signal strength modifier for each category used in the combination is reported in Fig. 9 (left). Results are generally consistent with unity, with the largest deviation showing up in the 2-jet VH-tagged category (i.e., the category targeting the associated production of a Higgs boson with a vector boson decaying hadronically). The level of compatibility of the signal strength modifiers in each category with the combined signal strength modifier corresponds to an asymptotic p-value of 0.34.
Given the sensitivity of the analysis to various production mechanisms, a fit is performed in which a signal strength modifier is assigned to each production mechanism, i.e., µ ggH , µ VBF , µ WH , and µ ZH . A simultaneous fit to all categories is performed, and results are shown in Fig. 9 (right). The biggest deviation from unity is observed for the WH production mechanism, which is probed mainly by the 2-jet VH-tagged and 3-lepton WH-tagged categories. The level of compatibility of the signal strength modifiers associated with different production mechanisms with the combined signal strength modifier corresponds to an asymptotic p-value of 0.70.
A similar simultaneous fit has been performed to measure the cross section ratios corresponding to five Higgs boson production mechanisms, using a simplified fiducial phase space, as specified in the "stage-0" simplified template cross section framework [35]. The aim of this framework is to define simple generator-level requirements that define specific phase space regions, allowing one to reduce the dependence of the measurement on the theoretical uncertainty in the SM predictions, which would otherwise apply in the case of the extrapolation to the full phase space. The cross sections corresponding to five Higgs boson production processes (σ ggH , σ VBF , σ WH lep. , σ ZH lep. , σ VH had. ) are measured requiring the generator-level Higgs boson rapidity to be |y H | < 2.5. This analysis has a negligible acceptance for Higgs boson production above |y H | = 2.5. The measured cross sections, divided by the SM predictions, for the production channels in which the analysis has sensitivity are shown in Fig. 10. The observed deviation of the σ VH had. process with respect to the SM prediction corresponds to an asymptotic p-value of 0.02, and is driven by the excess of events already observed for µ WH . Compared to the µ WH fit, in this case the signal strength modifier for the hadronic decay of the associated W boson is fitted separately from the leptonic one, and is driven away from the SM prediction by the excess observed in the 2-jet VH-tagged category.

Higgs boson couplings
Given its large cross section times branching fraction, the H → WW channel has a significant potential for constraining the Higgs boson couplings to fermions and vector bosons. A fit is performed to probe these couplings. One signal strength modifier (µ F ) is used to scale fermioninduced production mechanisms, i.e., ggH, ttH, and bbH, and another one (µ V ) scales the production mechanisms associated with vector bosons, i.e., VBF and VH. The two-dimensional likelihood profile is shown in Fig. 11 (left), where the 68% and 95% CL contours in the (µ F , µ V ) plane are displayed. The best fit values for each signal strength modifier are µ F = 1.37 +0.21 −0.20 and µ V = 0.78 +0. 60 −0.57 . The determination of the Higgs boson coupling constants is a way to verify the theoretical predictions and to search for deviations with respect to the SM expectations. These couplings can be parametrized using two coupling modifiers associated either with fermion or vector   Figure 10: Observed cross sections for the main Higgs boson production modes, normalized to the SM predictions. Cross section ratios are measured in a simplified fiducial phase space defined requiring y H < 2.5, as specified in the "stage-0" simplified template cross section framework. The vertical line and band correspond to the SM prediction and associated theoretical uncertainty.
boson vertices, using the so-called κ-framework [35]. The two coupling modifiers are used to scale the expected product of cross section and branching fraction to match the observed signal yields in the data, according to the following formula: where κ H = κ H (κ F , κ V ) is the Higgs boson total width modifier, defined as a function of the two fit parameters κ F and κ V . The κ i coupling modifier is equal to κ F for the ggH, ttH, and bbH production modes, and to κ V for the VBF and VH production modes. No processes other than SM ones are considered to contribute to the total width modifier. The two-dimensional likelihood profile obtained using this approach, and the corresponding 68% and 95% CL contours, are shown in Fig. 11 (right). The best fit values for each coupling modifier are κ F = 1.52 +0.48  Figure 11: Two-dimensional likelihood profile as a function of (left) the signal strength modifiers associated with either fermion (µ F ) or vector boson (µ V ) couplings, and (right) the coupling modifiers associated with either fermion (κ F ) or vector boson (κ V ) vertices, using the κ-framework parametrization. The 68% and 95% CL contours are shown as continuous and dashed lines, respectively. The red circle represents the best fit value, while the black triangle corresponds to the SM prediction.

Summary
Measurements of the properties of the SM Higgs boson decaying to a W boson pair at the LHC have been reported. The data samples used in the analysis correspond to an integrated luminosity of 35.9 fb −1 collected by the CMS detector in proton-proton collisions at √ s = 13 TeV.
The W + W − candidates are selected in events with large missing transverse momentum and exactly two, three, or four leptons. In the case of events with two leptons, different categories are defined according to the lepton pair flavor, eµ, ee, or µµ. The analysis has specific categories for gluon fusion production, vector boson fusion, and vector boson associated production, with up to two jets in the final state.
The probability of observing a signal at least as large as the one seen by combining all channels, under the background-only hypothesis, corresponds to an observed significance of 9.1 standard deviations for m H = 125.09 GeV, to be compared with the expected value of 7.1 standard deviations. The observed global signal strength modifier is σ/σ SM = µ = 1.28 +0.18 −0.17 = 1.28 ± 0.10 (stat) ± 0.11 (syst) +0.10 −0.07 (theo). Measurements of the signal strength modifiers associated with the main Higgs boson production mechanisms are also performed, as well as measurements of the Higgs boson couplings to fermions and vector bosons. The measured Higgs boson production and decay properties are found to be consistent, within their uncertainties, with the SM expectation.

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 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 and the CMS detector provided by the following funding agencies:  [10] ATLAS Collaboration, "Measurements of the Higgs boson production and decay rates and coupling strengths using pp collision data at √ s = 7 and 8 TeV in the ATLAS experiment", Eur. Phys. J. C 76 (2016) 6, doi:10.1140/epjc/s10052-015-3769-y, arXiv:1507.04548.