Measurement of the production cross section for a Higgs boson in association with a vector boson in the $H \rightarrow WW^{\ast} \rightarrow \ell\nu\ell\nu$ channel in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

TL;DR

This paper reports a measurement of the Higgs boson production cross sections in association with W and Z bosons, using the decay channel to WW* in proton-proton collisions at 13 TeV with the ATLAS detector, confirming Standard Model predictions.

Contribution

First measurement of associated Higgs production cross sections in the WW* decay channel at 13 TeV with ATLAS, providing data consistent with the Standard Model.

Findings

Measured cross sections agree with Standard Model predictions.

Quantified uncertainties in the production cross sections.

Provided the first such measurement at 13 TeV.

Abstract

A measurement of the Higgs boson production cross sections via associated $WH$ and $ZH$ production using $H \rightarrow WW^{\ast} \rightarrow \ell\nu\ell\nu$ decays, where $\ell$ stands for either an electron or a muon, is presented. Results for combined $WH$ and $ZH$ production are also presented. The analysis uses events produced in proton-proton collisions collected with the ATLAS detector at the Large Hadron Collider in 2015 and 2016. The data correspond to an integrated luminosity of 36.1 fb$ ^{-1}$ recorded at a centre-of-mass energy of 13 TeV. The products of the $H \rightarrow WW^{\ast}$ branching fraction times the $WH$ and $ZH$ cross sections are measured to be $0.67^{+0.31}_{-0.27}$(stat.)$^{+0.18}_{-0.14}$(syst.) pb and $0.54^{+0.31}_{-0.24}$(stat.)$^{+0.15}_{-0.07}$(syst.) pb respectively, in agreement with the Standard Model predictions.

Figures

Tables

Table 1. Table 1: Monte Carlo generators used to model the signal and background processes. Alternative generators, underlying event and parton-showering models, used to estimate systematic uncertainties, are shown in parentheses. In the last column the prediction order for the total cross section is shown. “P ythia 6” refers to version 6.428, “P ythia 8” refers to versions 8.210 or 8.186.

Process	Generator	UEPS model	Prediction order
	(alternative)	(alternative)	for total cross section
$q\overline{q}\to WH$	Powheg-Box v2 MiNLO	Pythia8	NNLO QCD + NLO EW [16, 17, 18]
		(Herwig 7)
$q\overline{q}\to ZH$	Powheg-Box v2 MiNLO	Pythia8	NNLO QCD + NLO EW [16, 17, 18]
		(Herwig 7)
$gg\to ZH$	Powheg-Box v2	Pythia8	NLO + NLL [19]
		(Herwig 7)
ggF $H$	Powheg-Box v2 NNLOPS	Pythia8	NNNLO QCD + NLO EW [20]
VBF $H$	Powheg-Box v2	Pythia8	NNLO QCD + NLO EW [21]
$t\overline{t}$	Powheg-Box v2	Pythia8	NNLO+NNLL [22]
		(Herwig 7)
	(Sherpa 2.2.1)	(Sherpa 2.2.1)
$Wt$	Powheg-Box v1	Pythia6	NLO [23]
$t\overline{t}W/Z$	MG5_aMC@LO	Pythia8	NLO [24, 25]
$tZ$	MG5_aMC@LO	Pythia6	LO [26]
$q\overline{q}/qg\to \ell \nu \ell \ell$	Sherpa 2.2.2	Sherpa 2.2.2	NLO [27]
	(Powheg-Box v2)	(Herwig++)
$q\overline{q}/qg\to \ell \ell \ell \ell$	Sherpa 2.1/2.2.2	Sherpa 2.1/2.2.2	NLO [27]
	(Powheg-Box v2)	(Herwig++)
$gg\to \ell \ell \ell \ell$	Sherpa 2.1.1	Sherpa 2.1.1	NLO [28]
$VVV$	Sherpa 2.2.2	Sherpa 2.2.2	NLO [29]
	(MG5_aMC@NLO)	(Pythia8)

Table 2. Table 2: Event selection criteria used to define the signal regions in the W ​ H 𝑊 𝐻 WH and Z ​ H 𝑍 𝐻 ZH analyses. The symbols are defined in the text.

	$WH$		$ZH$
Preselection	3 isolated leptons ($p_{\text{T}}$ $>$ 15  GeV)		4 isolated leptons ($p_{\text{T}}$ $>$ 10  GeV)
	total lepton charge $\pm 1$		total lepton charge 0
Category	$Z$-dominated	$Z$-depleted	2-SFOS	1-SFOS
Number of SFOS	2 or 1	0	2	1
Number of jets	$\le 1$	—	$\le 1$	$\le 2$
Number of $b$-jets	0	0	0	0
$E_{\text{T}}^{\text{miss}}$ [GeV]	$>30$	—	$>45$	—
$p_{\text{T}}^{4\ell }$ [GeV]	—	—	$>45$	—
$m_{\ell \ell }$ [GeV]	$>12$ (min. SFOS)	—	$>10$	$>10$
$\left|m_{\ell \ell }-m_Z\right|$ [GeV]	$>25$ (SFOS)	—
$m_{{\ell }_0{\ell }_1}$ [GeV]	—	—
$\Delta {\varphi }_{{\ell }_0{\ell }_1}^{\mathrm{boost}}$	—	—
$m_{\tau \tau }$ [GeV]	—	—	—
$\Delta {\varphi }_{{\ell }_0{\ell }_1,{\overrightarrow{p}}_{\text{T}}^{\text{ miss}}}$ [rad]	—	—	—	$>0.4$
$m_{4\ell }$ [GeV]	—	—	$>140$	—
BDT	BDT${}_{\text{Zdom}}{}^{}>0.3$	${\text{BDT}}_{t\overline{t}}>0.2$ ${\text{BDT}}_{WZ}>0.15$	—	—

Table 3. Table 3: Summary of the Z 𝑍 Z -depleted fit regions.

Bin 1	Bin 2	Bin 3	Bin 4	Bin 5	Bin 6
				$0.45\le {\mathrm{BDT}}_{t\overline{t}}$
	$0.35\le {\mathrm{BDT}}_{WZ}$		$0.35\le {\mathrm{BDT}}_{WZ}$		$0.35\le {\mathrm{BDT}}_{WZ}$

Table 4. Table 4: Definition of control regions in the W ​ H 𝑊 𝐻 WH and Z ​ H 𝑍 𝐻 ZH analyses. Selections indicated in boldface font are designed to ensure the control region (CR) is orthogonal to the relevant SR. In the last line normalisation factors which scale the corresponding yields in the signal region are shown with their uncertainties, including both the statistic and the systematic uncertainties.

Channel (Category)	$WH$ ($Z$-dominated and $Z$-depleted)		$ZH$ (2-SFOS)
CR	$WZ$	Top-quark	$ZZ$
Number of leptons	3	3	4
Total lepton charge	$\pm 1$	$\pm 1$	0
Number of SFOS	2 or 1	$\le 2$	2
Number of jets	$\le 1$	$\ge 1$	$\le 1$
Number of $b$-jets	0	$\ge \mathrm{𝟏}$	0
$E_{\text{T}}^{\text{miss}}$ [GeV]	$>50$	$>50$ (2 or 1 SFOS)	—
$\left|m_{\ell \ell }-m_Z\right|$ [GeV]		$>25$ (2 or 1 SFOS)
$m_{{\ell }_0{\ell }_1}$ [GeV]	—	—	$>\mathrm{𝟓𝟓}$
$m_{\ell \ell }$ (min. SFOS) [GeV]	$>12$	$>12$ (2 or 1 SFOS)	—
Normalisation factors	0.99$\pm$0.05	0.97$\pm$0.08	1.13$\pm$0.06

Table 5. Table 5: Post-fit predictions and data yields in the four SRs. The uncertainties include those from the sample statistics, and the theoretical and experimental systematic uncertainties. The sum of the single contributions may differ from the total value due to rounding. Moreover, the total uncertainty differs from the sum in quadrature of the single-process uncertainties due to the correlations. “Other Higgs” contains Higgs production mechanisms and decay processes different from V ​ H 𝑉 𝐻 VH and H → W ​ W ∗ → 𝐻 𝑊 superscript 𝑊 ∗ H\to WW^{\ast} respectively, except for the H → Z ​ Z ∗ → 𝐻 𝑍 superscript 𝑍 ∗ H\to ZZ^{\ast} contribution which is included in the “ Z ​ Z ∗ 𝑍 superscript 𝑍 ∗ ZZ^{\ast} ” row. The corresponding pre-fit predictions for the background processes differ only by the normalization factors listed in Table  4 .

Process	$WH$				$ZH$
	$Z$-dominated		$Z$-depleted		1-SFOS		2-SFOS
$WH$	$11$$\pm$	$6$	$5.8$$\pm$	$2.8$	—
$ZH$	$1.1$$\pm$	$0.6$	$0.61$$\pm$	$0.34$	$3.3$$\pm$	$1.7$	$1.8$$\pm$	$0.9$
$WZ/W{\gamma }^{\ast }$	$40.1$$\pm$	$2.8$	$1.7$$\pm$	$0.5$	—
$Z{Z}^{\ast }$	$2.4$$\pm$	$1.1$	$0.27$$\pm$	$0.09$	$0.14$$\pm$	$0.14$	$1.2$$\pm$	$0.3$
$VVV$	$1.5$$\pm$	$0.1$	$0.71$$\pm$	$0.11$	$0.32$$\pm$	$0.05$	$0.20$$\pm$	$0.03$
$tV/t\overline{t}V$	$0.14$$\pm$	$0.03$	$0.13$$\pm$	$0.03$	$0.04$$\pm$	$0.02$	$0.03$$\pm$	$0.01$
Other top-quark	$8.4$$\pm$	$2.6$	$1.9$$\pm$	$0.8$	—
Other Higgs	$0.31$$\pm$	$0.03$	$0.06$$\pm$	$0.01$	<0.01		$0.04$$\pm$	$0.01$
Misid. leptons	$9.7$$\pm$	$3.4$	<0.1		$0.19$$\pm$	$0.08$	$0.36$$\pm$	$0.12$
Total background	$62$$\pm$	$5$	$4.7$$\pm$	$1.0$	$0.65$$\pm$	$0.17$	$1.8$$\pm$	$0.3$
Observed	76		10		5		2

Table 6. Table 6: Breakdown of the main contributions to the total uncertainty in σ W ​ H ⋅ ℬ H → W ​ W ∗ ⋅ subscript 𝜎 𝑊 𝐻 subscript ℬ → 𝐻 𝑊 superscript 𝑊 \sigma_{WH}\cdot\mathcal{B}_{H\rightarrow WW^{*}} (left) and σ Z ​ H ⋅ ℬ H → W ​ W ∗ ⋅ subscript 𝜎 𝑍 𝐻 subscript ℬ → 𝐻 𝑊 superscript 𝑊 \sigma_{ZH}\cdot\mathcal{B}_{H\rightarrow WW^{*}} (right). The individual sources of systematic uncertainties are grouped together. The sum in quadrature of the individual components differs from the total uncertainty due to correlations between the components. Systematic uncertainties that affect the shape of the fitted distribution are indicated by an asterix.

Source	$\Delta {\sigma }_{WH}\cdot {ℬ}_{H\to W{W}^{\ast }}$ $[\%]$
Data statistics in SR	$43$
Data statistics in CR	$6$
Theoretical uncertainties	$16$
$WZ/W{\gamma }^{\ast }$(*)	$12$
Top-quark(*)	$8$
$WH$ signal(*)	$4$
$Z{Z}^{\ast }$	$2$
Experimental uncertainties	$12$
Impact parameter mismodelling	$8$
Misidentified leptons	$8$
$b$-tagging	$1$
MC statistics	$9$
Luminosity	$3$
TOTAL	$49$

Equations

$${\mu }_{WH}$$

$${\mu }_{ZH}$$

$${\mu }_{VH}$$

$${\sigma }_{WH}\cdot {\mathcal{B}}_{H\to W{W}^{\ast }}$$

$${\sigma }_{ZH}\cdot {\mathcal{B}}_{H\to W{W}^{\ast }}$$

Full text

\LEcontact

Pat Ward cpw1@hep.phy.cam.ac.uk

\AtlasTitleMeasurement of the production cross section for a Higgs boson in association with a vector boson in the $H\rightarrow WW^{\ast}\to\ell\nu\ell\nu$ channel in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector \AtlasAbstract A measurement of the Higgs boson production cross sections via associated $WH$ and $ZH$ production using $H\rightarrow WW^{\ast}\rightarrow\ell\nu\ell\nu$ decays, where $\ell$ stands for either an electron or a muon, is presented. Results for combined $WH$ and $ZH$ production are also presented. The analysis uses events produced in proton–proton collisions collected with the ATLAS detector at the Large Hadron Collider in 2015 and 2016. The data correspond to an integrated luminosity of $36.1\,\mathrm{fb}^{-1}$ recorded at a centre-of-mass energy of 13 TeV. The products of the $H\rightarrow WW^{\ast}$ branching fraction times the $WH$ and $ZH$ cross sections are measured to be $0.67^{+0.31}_{-0.27}$(stat.)${}^{+0.18}_{-0.14}$(syst.) pb and $0.54^{+0.31}_{-0.24}$(stat.)${}^{+0.15}_{-0.07}$(syst.) pb respectively, in agreement with the Standard Model predictions.

\AtlasRefCodeHIGG-2017-14 \PreprintIdNumberCERN-EP-2019-038 \AtlasJournalPhys. Lett. B. \AtlasJournalRefPhys. Lett. B 798 (2019) 134949 \AtlasDOI10.1016/j.physletb.2019.134949 \AtlasCoverSupportingNoteSupporting Notehttps://cds.cern.ch/record/2295764 \AtlasCoverCommentsDeadline??? \AtlasCoverAnalysisTeamAntonio Baroncelli, Kathrin Becker, Claudia Bertella, Michela Biglietti, Carsten Burgard, Dominik Duda, Ada Farilla, Frank Fiedler, Frank Filthaut, Ralf Gugel, Yesenia Hernandez Jimenez, Sana Ketabchi, Karsten Koeneke, Peter Krieger, Aaron Liblong, Bruce Mellado, Sophio Pataraia, Christian Schmitt, Marco Sessa, Jonas Strandberg, Marcel Weirich, Natalie Wieseotte \AtlasCoverEdBoardMemberZhiqing Zhang (chair), Aaron Armbruster, Andrea Bocci \AtlasCoverEgroupEditorsatlas-HIGG-2017-14-editors@cern.ch \AtlasCoverEgroupEdBoardatlas-HIGG-2017-14-editorial-board@cern.ch

1 Introduction

Higgs boson production in association with a $W$ or $Z$ boson, which is respectively denoted by $WH$ and $ZH$, and collectively referred to as $VH$ associated production in the following, provides direct access to the Higgs boson couplings to weak bosons. In particular, in the $WH$ mode with subsequent $H\rightarrow WW^{*}$ decay, the Higgs boson couples only to W bosons, at both the production and decay vertices.

This paper presents a measurement of the corresponding production cross sections through the decay $H\rightarrow WW^{*}\rightarrow\ell\nu\ell\nu$, using proton–proton collisions at a centre-of-mass energy of $\sqrt{s}=$ 13 TeV. The data correspond to an integrated luminosity of $36.1$ fb*-1* and were recorded by the ATLAS detector at the Large Hadron Collider (LHC). Previous measurements at $\sqrt{s}=8$ TeV were performed by the ATLAS [1] and CMS [2] Collaborations and recently at $\sqrt{s}=13$ TeV with $35.9\,\mathrm{fb}^{-1}$ of data by the CMS Collaboration [3]. Recent results at $\sqrt{s}=13$ TeV on $VH$ production in other decay modes can be found in Refs. [4, 5, 6, 7, 8, 9].

The analysis is performed using events with three ($3\ell$) or four ($4\ell$) charged leptons (electrons or muons) in the final state, targeting the $WH$ and $ZH$ channels respectively. Leptonic decays of $\tau$ leptons, from $H\rightarrow WW^{*}\rightarrow\tau\nu\tau\nu$ or $H\rightarrow WW^{*}\rightarrow\tau\nu\ell\nu$ or from the associated vector bosons, are considered as signal, while no specific selection is performed for events with hadronically decaying $\tau$ leptons in the final state. Events from $VH$ production with $H\rightarrow\tau\tau$ are considered as background. The leading-order Feynman diagrams for the $WH$ and $ZH$ production processes are depicted in Figure 1.

In the $WH$ channel, multivariate discriminants are used to maximise the sensitivity to the Higgs boson signal, while in the $ZH$ channel the analysis is performed through selection requirements. The distribution of these $WH$ discriminants, together with event counts in background control regions and the signal regions in the $ZH$ channel, are combined in a binned maximum-likelihood fit to extract the signal yield and the background normalisations. The maximum-likelihood fit provides results for the $WH$ and the $ZH$ channels separately and for their combination $VH$, assuming the Standard Model (SM) prediction for the relative cross sections of the two production processes.

2 ATLAS detector

The ATLAS experiment [10, 11, 12] is a multi-purpose particle detector with a forward–backward symmetric cylindrical geometry and a near $4\pi$ coverage in solid angle.111 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the $z$-axis along the beam pipe. The $x$-axis points from the IP to the centre of the LHC ring, and the $y$-axis points upwards. Cylindrical coordinates $(r,\phi)$ are used in the transverse plane, $\phi$ being the azimuthal angle around the $z$-axis. The pseudorapidity is defined in terms of the polar angle $\theta$ as $\eta=-\ln\tan(\theta/2)$. Angular distance is measured in units of $\Delta R\equiv\sqrt{(\Delta\eta)^{2}+(\Delta\phi)^{2}}$. Transverse momentum and energy are defined as $p_{\text{T}}=p\sin\theta$ and $E_{\text{T}}=E\sin\theta$ respectively. It consists of an inner tracking detector (ID) surrounded by a thin superconducting solenoid providing a $2\text{\,}\mathrm{T}$ axial magnetic field, electromagnetic (EM) and hadronic calorimeters, and a muon spectrometer (MS). The inner tracking detector covers the pseudorapidity range $|\eta|<2.5$. It consists of silicon pixel, silicon micro-strip, and transition-radiation tracking detectors. Lead/liquid-argon (LAr) sampling calorimeters provide electromagnetic energy measurements with high granularity. A hadronic (steel/scintillator-tile) calorimeter covers the central pseudorapidity range ($|\eta|<1.7$). The endcap and forward regions are instrumented with LAr calorimeters for both EM and hadronic energy measurements up to $|\eta|=4.9$. The muon spectrometer surrounds the calorimeters and is based on three large air-core toroidal superconducting magnet systems that provide a field integral between $2.0$ and $6.0\text{\,}\mathrm{T}\text{\,}\mathrm{m}$ across most of the detector. The muon spectrometer includes a system of precision tracking chambers covering the region $|\eta|<2.7$ and fast detectors for triggering within the range $|\eta|<2.4$. A two-level trigger system is used to select events [13].

3 Signal and background Monte Carlo simulation

Monte Carlo (MC) event generators are used to model signal and background processes. All signal samples were generated with a Higgs boson mass of 125 GeV [14, 15]. For most processes, separate programs were used to generate the hard scattering process and to model the underlying event and the parton showering (UEPS). A description of the MC samples is given in Table 1. They are normalised to cross-section predictions calculated with the QCD and electroweak (EW) orders specified in the last column of Table 1.

The $q\bar{q}\rightarrow WH$ and $q\bar{q}\rightarrow ZH$ processes were generated with Powheg-Box v2 [30] MiNLO interfaced to Pythia8 [31], with the AZNLO set of tuned parameters (tune) [32] for underlying event, showering and hadronisation. The $gg\to ZH$ process was simulated with Powheg-Box v2 + Pythia8 with the AZNLO tune for underlying event, showering and hadronisation. For the $VH$ samples, the PDF4LHC15 parton distribution function (PDF) set [33] was used for the hard scattering process in Powheg-Box v2 and the CTEQ6L1 PDF set [34] was used for the parton showering in Pythia8. Herwig 7  [35], with the MMHT2014lo68cl PDF set [36], was used as an alternative parton-showering model for $VH$. The uncertainty due to the PDF choice is smaller than the uncertainty obtained by using Herwig as an alternative parton shower model (Section 7).

The gluon–gluon fusion (ggF) events were generated with Powheg-Box v2 NNLOPS [37] interfaced to Pythia8 with the AZNLO tune. The vector-boson fusion (VBF) events were generated with Powheg-Box v2, interfaced to Pythia8. For the ggF and VBF samples, the PDF4LHC15 PDF set was used for the hard scattering process in Powheg-Box v2 and the CTEQ6L1 PDF set was used for the parton showering in Pythia8. The contribution from the $t\bar{t}H$ and $tH$ production modes is negligible.

The top-quark pair production ($t\bar{t}$) was simulated with Powheg-Box v2 [38] using the NNPDF 3.0 NNLO PDF set [39] and interfaced to Pythia8 using the NNPDF 2.3 PDF set [40] for parton showering, with the A14 tune [41]. For $t\bar{t}$ production, Sherpa [42] 2.2.1, with the NNPDF 3.0 PDF set, was used as an alternative generator while Herwig 7, with the MMHT2014lo68cl PDF set, was used as an alternative UEPS model. The single-top-quark production $Wt$ was generated with Powheg-Box v1 [23] interfaced to Pythia6 [43] for parton showering with the Perugia2012 tune [44]. EvtGen 1.2.0 [45] was used for the simulation of b-quark and c-quark decays. The $t\bar{t}W/Z$ and $tZ$ processes were generated at leading order (LO) with MG5_aMC@LO [25] version 2.2.2 ($t\bar{t}W/Z$) and 2.2.1 ($tZ$) interfaced to Pythia8 ($t\bar{t}W/Z$) and Pythia6 ($tZ$), using the NNPDF2.3 LO PDF set.

The $q\bar{q}/qg\to VV^{\ast}$ samples with final states $\ell\nu\ell\ell$ and $\ell\ell\ell\ell$ [46] were generated with Sherpa 2.2.2, with the exception of the $ZZ^{\ast}$ sample in the $WH$ analysis for which Sherpa 2.1 was used; the CT10 PDF set [47] and the NNPDF 3.0 PDF set were used for versions 2.1 and 2.2.2, respectively. Powheg-Box v2 [48] was used as an alternative generator for $VV^{\ast}$, with Herwig++, using the CTEQ6L1 PDF set, for parton showering. Among the loop-induced $gg$-initiated diboson processes, the only relevant process in this analysis is $gg\to ZZ^{\ast}$, for which a $K$-factor of 1.55 was used [28]. This process was simulated with Sherpa 2.1.1, using the CT10 PDF set.

The triboson $VVV$ samples were generated with Sherpa 2.2.2 and the NNPDF 3.0 PDF set. MG5_aMC@NLO was used as an alternative generator for $VVV$, with Pythia8, using the NNPDF2.3 LO PDF set. The same PDF sets were used for the hard scattering and the parton showering in all the Sherpa samples described above.

All simulated samples include the effect of pile-up from multiple interactions in the same and neighbouring bunch crossings. This was achieved by overlaying minimum-bias events, simulated using Pythia8 with the A2 tune [49] and MSTW2008LO PDF set [50]. All samples were processed through the Geant 4 [51] ATLAS detector simulation [52].

4 Event reconstruction

Candidate signal events are selected using triggers that require a single isolated lepton with minimum transverse momentum ($p_{\text{T}}$) thresholds between 24 GeV and 26 GeV for electrons and between 20 GeV and 26 GeV for muons, depending on the data-taking period. At least one of the leptons reconstructed offline is required to have triggered the event and to have a $p_{\text{T}}$ higher than the nominal trigger threshold by at least 1 GeV. The single-lepton trigger efficiencies on the plateau are approximately 70% for single muons with $|\eta|<1.05$, 90% for single muons in the range $1.05<|\eta|<2.40$ and greater than 90% for single electrons in the range $|\eta|<2.47$. The trigger efficiency for the signal events, estimated after the preselection, is 94% for $WH$ and 98.5% for $ZH$.

Selected events are required to have at least one primary vertex reconstructed from at least two associated tracks, each with transverse momentum $p_{\text{T}}>400\leavevmode\nobreak\ \text{Me\kern-1.00006ptV}$, as described in Ref. [53]. If an event has more than one reconstructed primary vertex, the vertex with the largest track $\sum p_{\text{T}}^{2}$ is selected for the analysis.

Electrons are reconstructed from clusters of energy deposits in the EM calorimeter matched to ID tracks, and are identified using criteria based on the calorimeter shower shape, the quality of the match between the track and the cluster and the amount of transition radiation emitted in the ID, as described in Ref. [54]. Electrons are required to satisfy $|\eta|<2.47$, excluding $1.37<|\eta|<1.52$, which corresponds to the transition region between the barrel and the endcap EM calorimeters. The efficiency for electron identification ranges from 88% to 94%, depending on electron $p_{\text{T}}$ and $\eta$. Muons are reconstructed by combining ID and MS tracks with consistent trajectories and curvatures. An overall fit of hits from the ID track, energy loss in the calorimeter and the hits of the track in the muon system is used to form muon candidates, as described in Ref. [55]. The efficiency for muon identification is close to 95% over the full instrumented $\eta$ range. To suppress particles misidentified as leptons, several identification requirements as well as impact parameter, calorimeter and track isolation criteria [54, 55] are applied.

Jets are reconstructed using the anti-$k_{t}$ algorithm with radius parameter $R=0.4$ [56, 57]. The four-momenta of jets are corrected for the effects of calorimeter non-compensation, energy loss in non-instrumented regions, and contributions from pile-up [58]. Jets are required to have $|\eta|<4.5$, with $p_{\text{T}}>25$ GeV for the region $|\eta|<2.5$ and $p_{\text{T}}>30$ GeV for the region $2.5<|\eta|<4.5$. A multivariate selection [59] is used to suppress jets with $p_{\text{T}}<60$ GeV and $|\eta|<2.4$ originated from pile-up. Furthermore, to suppress pile-up jets in the forward region, jet shapes and topological jet correlations in pile-up interactions are exploited [60]. Jets with $p_{\text{T}}$ > 20 GeV  and $|\eta|<2.5$ containing $b$-hadrons ($b$-jets) are identified using a multivariate technique [61] with an efficiency of 85%, estimated from simulated $t\bar{t}$ events. The multivariate technique gives rejection factors against jets originating from a light quark or gluon and jets containing $c$-hadrons of 33 and 3, respectively.

The missing transverse momentum $\vec{p}_{\text{T}}^{\text{\,miss}}$ with magnitude $E_{\text{T}}^{\text{miss}}$ in each event is calculated from the negative vectorial sum of the transverse momenta of electrons, muons, and jets. It uses both track-based and calorimeter-based measurements [62].

5 Event selection

In the $WH$ channel, exactly three isolated leptons with $p_{\text{T}}>15\leavevmode\nobreak\ \text{Ge\kern-1.00006ptV}$ are required with a total charge of $\pm 1$. The lepton with unique charge is labelled $\ell_{0}$, the lepton closest to $\ell_{0}$ in angular distance $\Delta R$ is labelled $\ell_{1}$, and the remaining lepton is labelled $\ell_{2}$. In signal events leptons $\ell_{0}$ and $\ell_{1}$ are most likely to originate from the $H\rightarrow WW^{*}$ decay, with probabilities of $99\%$ and $85\%$ respectively.

The most prominent background processes to the $WH$ channel are $WZ/W\gamma^{\ast}$ production and top-quark processes with either three prompt leptons, e.g. $t\bar{t}V$, or two prompt leptons and one non-prompt lepton from a $b$-hadron decay, e.g. $t\bar{t}$. Other important background processes are $ZZ^{\ast}$ (including $Z\gamma^{\ast}$), $Z\gamma$ and $Z$+jets production; they may satisfy the signal selection requirements if a lepton is undetected, in the case of $ZZ^{\ast}$, or if they contain a misidentified or non-prompt lepton, in the case of $Z\gamma$ and $Z$+jets production. Processes with three prompt leptons in the final state such as tribosons, in particular $WWW$, also contribute to the background. Contributions from background processes that include more than one misidentified lepton, such as $W$+jets production and inclusive $b{\bar{b}}$ pair production, are negligible. The background from top-quark production is suppressed by vetoing events if they contain any $b$-tagged jet.

The analysis of the $WH$ channel separates events with at least one same-flavour opposite-sign charge (SFOS) lepton pair from events with zero SFOS lepton pairs, which have different signal-to-background ratios. Due to the presence of $Z\rightarrow\ell\ell$ decays as a dominant background, the former is hereafter referred to as the $Z$-dominated category, while the latter is referred to as the $Z$-depleted category. In the $Z$-dominated category, the major background processes are those involving $Z$ bosons. Events are vetoed if they contain more than one jet. This requirement further suppresses top-quark events with an additional non-prompt lepton from $b$-hadron decays. In order to select final states with neutrinos, $E_{\text{T}}^{\text{miss}}$ is required to be above 30 GeV. The invariant masses $m_{\ell\ell}$ of all SFOS pairs are required to satisfy a $Z$-veto selection: $|m_{\ell\ell}-m_{Z}|>25\leavevmode\nobreak\ \text{Ge\kern-1.00006ptV}$. The last two requirements suppress $WZ/W\gamma^{\ast}$ and $ZZ^{\ast}$ events, and improve the $Z$+jets rejection. In order to suppress background events from heavy-flavour quarkonia, the smallest invariant mass of SFOS pairs is required to be greater than 12 GeV. The signal efficiency of this selection with respect to the preselection is $34.6\%$. A discriminant based on a boosted decision tree (BDT) [63, 64] is used to achieve a further separation between signal and background processes. The main purpose of the multivariate classifier, named BDT${}_{\textrm{Zdom}}$, is to distinguish between signal and the dominant $WZ/W\gamma^{\ast}$ and $ZZ^{\ast}$ background processes, and hence it is trained against these two background processes. The BDT${}_{\textrm{Zdom}}$ uses seven input variables. They are the magnitude of the vectorial sum of lepton transverse momenta ($|\Sigma{p_{\text{T}}^{\ell_{i}}}|$), the invariant masses of the first lepton pair ($m_{\ell_{0}\ell_{1}}$) and of the three leptons ($m_{\ell\ell\ell}$), the angular distance $\Delta R_{\ell_{0}\ell_{1}}$, $E_{\text{T}}^{\text{miss}}$ and the pseudorapidity separation between the leptons with the same charge ($\Delta\eta_{\ell_{1}\ell_{2}}$). Moreover the BDT uses the transverse mass of the $W$ boson, $m_{\textrm{T}}^{W}=\sqrt{2\,p_{\text{T}}^{\ell_{W}}\cdot E_{\text{T}}^{\text{miss}}\cdot\left(1-\cos{\Delta\phi\left({\ell_{W}},\vec{p}_{\text{T}}^{\text{\,miss}}\right)}\right)}$, built from the $\vec{p}_{\text{T}}^{\text{\,miss}}$ and the lepton $\ell_{W}$ which is the lepton not belonging to the SFOS lepton pair with invariant mass closer to the $Z$ boson mass, and could be either $\ell_{1}$ or $\ell_{2}$. Figure 2 shows the distribution of $\Delta R_{\ell_{0}\ell_{1}}$, which is the most powerful variable in the BDT${}_{\textrm{Zdom}}$ training, and the BDT${}_{\textrm{Zdom}}$ distribution in the $Z$-dominated category, before applying any selection requirement on the BDT${}_{\textrm{Zdom}}$ score.

The $Z$-dominated signal region (SR), defined as the events with high-ranking BDT${}_{\textrm{Zdom}}$ score (BDT{}_{\textrm{Zdom}}$$>0.3), is divided into three bins with increasing sensitivity: $0.3\leq\mathrm{BDT_{\textrm{Zdom}}}<0.5$, $0.5\leq\mathrm{BDT_{\textrm{Zdom}}}<0.7$ and $0.7\leq\mathrm{BDT_{\textrm{Zdom}}}<1.0$. This binning is the result of an optimisation to find the minimal number of BDT bins that gives the highest sensitivity. The expected signal-to-background ratio in these bins is about 0.05, 0.09 and 0.19, respectively. The full $Z$-dominated event selection is summarised in Table 2.

In the $Z$-depleted category, the two major background processes are $WZ/W\gamma^{\ast}$ with $Z/\gamma^{\ast}\to\tau\bar{\tau}$ and $t\bar{t}$, where $WZ/W\gamma^{\ast}$ has the same signature of the signal, namely three prompt leptons, while $t\bar{t}$ contains a misidentified or non-prompt lepton from a jet. Two separate BDTs, named $\text{BDT}_{WZ}$ and \text{BDT}_{\text{t\bar{t}}}, are trained to allow an optimal background rejection. The $\text{BDT}_{WZ}$ uses 11 input variables, of which three common to the BDT${}_{\textrm{Zdom}}$ are $m_{\ell_{0}\ell_{1}}$, $E_{\text{T}}^{\text{miss}}$ and $\Delta\eta_{\ell_{1}\ell_{2}}$; the other variables are the transverse momenta of the three leptons ($p_{\text{T}}^{\ell_{0}}$, $p_{\text{T}}^{\ell_{1}}$, $p_{\text{T}}^{\ell_{2}}$), the transverse mass ($m^{\textrm{T}}_{\ell_{0}\ell 1}$) built from $\ell_{0}$, $\ell_{1}$ and $\vec{p}_{\text{T}}^{\text{\,miss}}$, the invariant mass of the electrons with same-sign charge ($m_{ee}$), the transverse impact parameter significances of the lepton with lowest $p_{\text{T}}$ ($|d_{\textrm{0,sig,min}}|$), the transverse impact parameter significances of the lepton with second-lowest $p_{\text{T}}$ and opposite charge with respect to the lepton with lowest $p_{\text{T}}$ ($|d_{\textrm{0,sig,mid}}|$) and the compatibility of the event with the $WZ$ hypothesis $F_{\alpha}$.222Given the reconstructed charged-lepton momenta and the $\vec{p}_{\text{T}}^{\text{\,miss}}$, the event kinematics can be calculated under the $WZ$ with $Z\to\tau\tau$ hypothesis and using the collinear approximation for the $\tau$ decays with one remaining unknown, e.g. the ratio of one $\tau$ energy to the energy of the lepton from this $\tau$ decay. This unknown is varied and the number of physical kinematic solutions is taken as a measure of the compatibility with the $WZ$ hypothesis.

A definition of the most likely lepton from heavy-flavour decays ($\ell_{\textrm{HFL}}$) is crucial for an optimal performance of the \text{BDT}_{\text{t\bar{t}}}. For this purpose, a $\text{BDT}_{\text{HFL}}$ is trained purely on data using track and calorimeter isolation as well as impact-parameter variables as input. The lepton with the minimal $\text{BDT}_{\text{HFL}}$ output is selected as $\ell_{\textrm{HFL}}$. The \text{BDT}_{\text{t\bar{t}}} uses nine input variables, of which two common to the BDT${}_{\textrm{Zdom}}$ and $\text{BDT}_{WZ}$ are $m_{\ell_{0}\ell_{1}}$ and $\Delta\eta_{\ell_{1}\ell_{2}}$, one common to the BDT${}_{\textrm{Zdom}}$ is $\Delta R_{\ell_{0}\ell_{1}}$, and the other input variables are the number of jets ($N_{\textrm{jet}}$), the transverse momentum of the leading jet ($p_{\text{T}}^{j_{\textrm{lead}}}$), the invariant mass of the leptons with same-sign charge ($m_{\ell_{1}\ell_{2}}$), and three $\ell_{\textrm{HFL}}$-related variables: its $\text{BDT}_{\text{HFL}}$ output, its transverse momentum, $p_{\text{T}}^{\ell_{\textrm{HFL}}}$, and the invariant mass built from it together with the closest opposite-charge lepton ($m_{\ell_{\textrm{HFL}}\ell_{\textrm{cloc}}}$). Figure 3 shows the outputs of $\text{BDT}_{WZ}$ and \text{BDT}_{\text{t\bar{t}}} in the $Z$-depleted category, before applying any selection requirements on the BDT scores.

The choice of input variables for the different BDTs was the result of an optimisation study where several thousand different BDTs using different set of input variables have been compared and the best performing BDTs have been selected for the final analysis.

The full $Z$-depleted event selection is summarised in Table 2. The events with high-ranking BDT scores (\text{BDT}_{\text{t\bar{t}}} > 0.2 and $\text{BDT}_{WZ}$ > 0.15) are used to define the $Z$-depleted SR. In this region, the BDT scores are used as discriminant variables in the statistical analysis, with three bins in \text{BDT}_{\text{t\bar{t}}}, each being subdivided into two bins in $\text{BDT}_{WZ}$ as shown in Table 3. The expected signal-to-background ratio in these bins ranges from about $0.07$ in the first bin up to about $3.6$ in the last bin. The signal efficiency with respect to the preselection, before any cut on the BDT scores, is $23.5\%$.

The $ZH$ channel requires events to contain four isolated leptons with $p_{\textrm{T}}>10$ GeV and total electric charge of zero. Events that contain a SFOS lepton pair with $m_{\ell\ell}<10$ GeV are rejected to suppress the contamination from heavy-flavour quarkonia. Following this preselection, events are classified according to the number of SFOS lepton pairs: 1-SFOS and 2-SFOS. Events with no SFOS lepton pairs are not considered.

The reconstruction of the $ZH$ process proceeds through the identification of the leptons from the $Z$ boson, called $\ell_{2}$ and $\ell_{3}$, as the SFOS lepton pair with invariant mass closest to the $Z$ boson mass, $m_{Z}$. Then, the remaining two leptons, labelled $\ell_{0}$ and $\ell_{1}$, are candidates for originating from the Higgs boson decay. The background in the $ZH$ channel is almost exclusively due to $ZZ^{\ast}$ production. This constitutes $\sim$92% of the total background after the preselection is applied. Processes with four prompt leptons in the final state such as triboson production, in particular $ZWW$, which has the same signature as the signal, and $t\bar{t}Z$ also contribute to the background. Other background processes such as $WZ/W\gamma^{\ast}$, $Z$+jets and $tV$ may contribute when at least one jet, hadron or a converted photon is misidentified as a lepton.

In order to suppress the $t\bar{t}Z$ process, events containing $b$-tagged jets are rejected and at most one and two jets are allowed in 2-SFOS and 1-SFOS classes, respectively. To reduce the $ZZ^{\ast}$ background process in events with two SFOS lepton pairs, a threshold of 45 GeV is applied to the $E_{\textrm{T}}^{\textrm{miss}}$ and to the vector sum of the lepton transverse momenta, $p_{\textrm{T}}^{4\ell}$. The invariant mass of $\ell_{2}$ and $\ell_{3}$, $m_{\ell_{2}\ell_{3}}$, is required to satisfy $|m_{\ell_{2}\ell_{3}}-m_{Z}|<10$ GeV, and the invariant mass of $\ell_{0}$ and $\ell_{1}$, $m_{\ell_{0}\ell_{1}}$, is required to be between 10 GeV and 60 GeV (55 GeV) in 1-SFOS (2-SFOS) events. The variable $\Delta\phi_{\ell_{0}\ell_{1}}^{\mathrm{boost}}$ denotes the difference in azimuthal angle between the leptons from the Higgs boson candidate in the frame where the Higgs boson $p_{\textrm{T}}$ is zero. The Higgs boson transverse momentum is approximated with $\vec{p}_{\mathrm{T}}^{H}\approx-\vec{p}_{\mathrm{T}}^{Z}$, or $\vec{p}_{\mathrm{T}}^{H}\approx-\vec{p}_{\mathrm{T}}^{Z}-\sum\vec{p}_{\mathrm{T}}^{\mathrm{jet}}$ if at least one jet is present in the event. Events are required to satisfy $\Delta\phi_{\ell_{0}\ell_{1}}^{\mathrm{boost}}<1.9\leavevmode\nobreak\ (2.3)$ rad in the 1-SFOS (2-SFOS) class. In 1-SFOS events the $ZZ^{\ast}$ process contributes through the $Z\rightarrow\tau\tau$ decay, therefore the reconstructed mass of the $\tau$ pair, $m_{\tau\tau}$ is required to be below 50 GeV; $m_{\tau\tau}$ is computed using the collinear approximation method [65]. In addition, the azimuthal separation between the Higgs-candidate lepton pair and the $\vec{p}_{\text{T}}^{\text{\,miss}}$, $\Delta\phi_{\ell_{0}\ell_{1},\,\vec{p}_{\text{T}}^{\text{\,miss}}}$, is required to be above 0.4 rad. The final selection of these variables is optimised in order to maximise the signal significance. Figure 4(a) illustrates the discriminating power of the $\Delta\phi_{\ell_{0}\ell_{1},\,\vec{p}_{\text{T}}^{\text{\,miss}}}$ variable for 1-SFOS events after the preselection. This variable reduces the ZZ* contribution which becomes the dominant background source after the $b$-tagged jet veto and $m_{\ell_{0}\ell_{1}}$ selection are applied. Figure 4(b) shows the $E_{\textrm{T}}^{\textrm{miss}}$ distribution for 2-SFOS events after the preselection. In order to be orthogonal to the $H\rightarrow ZZ^{\ast}$ analysis of Ref. [66], 2-SFOS events are required to have an invariant mass of the four-lepton system, $m_{4\ell}$, above 140 GeV. The full event selection for $ZH$ is summarised in Table 2.

The total efficiency times acceptance of this selection for the process $WH$ with subsequent $H\rightarrow WW^{*}$ decay is about $0.073\%$ while for $ZH$ it is about $0.026\%$. These numbers are given with respect to all $W$ and $Z$ decays.

6 Background estimation

The main background contamination originates from processes with the same final state as the signal, namely diboson production ($WZ/W\gamma^{\ast}$, $ZZ^{\ast}$), top-quark processes with three or four prompt leptons such as $t\bar{t}V$, and triboson production. Other relevant background contributions arise from processes, such as $t\bar{t}$ or $Z+$jets, where the reconstructed leptons either originate from non-prompt decays of heavy-flavour hadrons or from jets misidentified as leptons.

Two dedicated regions, hereafter named control regions (CRs), are used to estimate the normalisation factors (NFs) of the main prompt background processes by fitting the expected yields from simulation to data: $WZ/W\gamma^{\ast}$ for the $WH$ channel and $ZZ^{\ast}$ for the $ZH$ channel in the 2-SFOS SR. In the 1-SFOS SR, $ZZ^{\ast}$ is estimated purely from simulation.

The CRs are made orthogonal to the corresponding SRs by inverting various selection criteria with respect to the SR definitions. The $WZ$ CR is defined by reversing the $Z$-veto in the $Z$-dominated $WH$ signal region. To improve the purity of the $WZ$ CR, the minimum $E_{\text{T}}^{\text{miss}}$ is increased from 30 GeV to 50 GeV. The $ZZ$ CR is defined by inverting the $m_{\ell_{0}\ell_{1}}$ requirement defined in the $ZH$ 2-SFOS SR. In order to increase the number of events, the $E_{\text{T}}^{\text{miss}}$, $\Delta\phi_{\ell_{0}\ell_{1}}^{\mathrm{boost}}$ , $p_{\text{T}}^{4\ell}$ , and $m_{4\ell}$ requirements are removed. Table 4 summarises the event selection for the $WZ$ and $ZZ$ CRs and the NFs for the background processes, obtained from the fit described in Section 8. The NFs are hereby completely dominated by the CR statistics and their values do not change when only the control regions are used in the fit.

Figures 5(a) and 5(b) show the distributions of $m_{\ell_{0}\ell_{1}}$ in the $WZ$ CR and the invariant mass of the four leptons in the $ZZ$ CR, as obtained from the final fit in the statistical analysis (post-fit) described in Section 8.

The background contributions with misidentified leptons are estimated using different techniques. The top-quark background in the $WH$ channels is normalised using a CR (top-quark CR) defined by requiring at least one $b$-tagged jet. To improve the purity of the top-quark CR, the minimum $E_{\text{T}}^{\text{miss}}$ is increased from 30 GeV to 50 GeV if at least one SFOS pair is present in the event. Processes with one misidentified lepton ($t\bar{t}$ and $Wt$) constitute 94$\%$ of the top-quark CR, the remaining events contain three prompt leptons from $t\bar{t}W$ decays. The full selection requirements applied to define the top-quark CR and the measured NF are also summarised in Table 4. Figure 5(c) shows the $m_{\ell_{0}\ell_{1}}$ distribution in the top-quark CR as obtained from the final fit.

A data-driven technique is used to estimate the $Z$+jets and $Z\gamma$ contribution in both the $WH$ and $ZH$ channels and the contributions from $WZ/W\gamma^{\ast}$ and top-quark processes in the $ZH$ channel. These contributions typically have one misidentified lepton ($Z$+jets and $Z\gamma$ in the $WH$ channel, and $WZ/W\gamma^{\ast}$ and $tZ/t\bar{t}W$ in the $ZH$ channel) or two misidentified leptons ($Z$+jets and $tW/t\bar{t}$ in the $ZH$ channel). A control sample where one or two of the lepton candidates fail to meet the nominal identification or isolation criteria but satisfy looser identification criteria, referred to as anti-identified leptons, is used. The contribution from misidentified leptons in the SR is then obtained by scaling the number of events in the control sample by extrapolation factors measured in a data sample enriched in $Z$+jets events. The latter is obtained by selecting events with two prompt leptons from a $Z$ boson decay and a loosely identified lepton considered to be the misidentified lepton candidate. The extrapolation factors are defined as the ratio of the number of misidentified lepton candidates that pass the nominal identification criteria to the number that pass the anti-identification criteria. In both the control sample and the data samples enriched in $Z$+jets events, the contribution from background events not estimated with this method is subtracted using MC expectations. Details of this method can be found in Ref. [67]. The uncertainty in the data-driven background processes described in this section includes the statistical uncertainty in the $Z$+jets enriched sample, the uncertainty from $Z$+jets MC modelling, and the theory uncertainty from the subtraction of other processes. For the $ZH$ channel, which can have events with two prompt leptons and two misidentified leptons, the uncertainty in the extrapolation from the control sample to the SRs is also included.

7 Systematic uncertainties

The systematic uncertainties can be categorised into those arising from experimental sources and those from theoretical sources. The dominant experimental uncertainties come from the misidentification of leptons (see Section 6), the mismodelling of the impact-parameter significance, the b-tagging efficiency [61] and the jet energy scale and resolution [58]. Other sources of uncertainty are due to the modelling of pile-up, the calibration of the missing transverse momentum measurement [62], and the luminosity measurement. The uncertainty in the combined 2015+2016 integrated luminosity is 2.1%. It is derived from the calibration of the luminosity scale using $x$-$y$ beam-separation scans, following a methodology similar to that detailed in Ref. [68], and using the LUCID-2 detector for the baseline luminosity measurements [69].

The impact of the uncertainties on lepton energy (momentum) scale and resolution, and identification and isolation criteria [70, 54, 55] is negligible. The experimental uncertainties are varied in a correlated way across all background processes and all signal- and control-region bins, so that uncertainties in the extrapolation from control to signal regions are correctly propagated. The luminosity uncertainty is only applied to background processes that are normalised to theoretical predictions, and to the Higgs boson signal. The theoretical uncertainties are evaluated by comparing nominal and alternative event generators and UEPS models as described in Section 3 and by varying PDF sets and the QCD renormalisation and factorisation scales. The uncertainty due to the PDF choice for the signal process is $1\%$ while the uncertainty obtained by using Herwig as an alternative parton shower model is $3-10\%$, depending on the signal region. All uncertainties are propagated through the full analysis chain and treated as being bin-dependent and region-dependent, i.e. potentially modifying not only the normalisation but also the shape of the BDT output distributions. Whenever the influence on the shape is found to be negligible, as in the case of the PDF and scale variations, only the normalisation uncertainties are used. A list of the systematic uncertainty sources and their impact on the cross-section measurement are shown later in Section 8.

8 Results

A binned likelihood function is constructed as a product of Poisson probability terms over the eleven bins of the different SRs defined in Section 5. The function has two independent scaling parameters: the signal strength parameter $\mu$, defined as the ratio of the measured signal yield to that predicted by the SM, for each of the $WH$ and the $ZH$ processes. Additionally, one Poisson probability term is added for each CR to determine simultaneously the normalisation of the corresponding background processes. Systematic uncertainties enter as nuisance parameters in the likelihood function and their correlations are taken into account. The final results are obtained using the profile likelihood method [71]. The resulting post-fit prediction and data yields in the four SRs are shown in Table 5.

Figure 6 shows the post-fit distribution of BDT scores in the $WH$ $Z$-dominated and $Z$-depleted SRs. The post-fit event yields in the $ZH$ 1-SFOS and 2-SFOS SRs are shown in Figure 7.

The measured signal strengths for the $WH$ and $ZH$ production modes in the $H\rightarrow WW^{*}$ decay channel are simultaneously determined to be

[TABLE]

compared to the expected results of $\mu_{WH}=1.0^{+0.86}_{-0.69}(\mathrm{stat.})^{+0.40}_{-0.30}(\mathrm{exp\leavevmode\nobreak\ syst.})^{+0.32}_{-0.26}(\mathrm{theo\leavevmode\nobreak\ syst.})=1.0^{+1.0}_{-0.8}$ and $\mu_{ZH}=1.0^{+1.3}_{-0.9}(\mathrm{stat.})^{+0.41}_{-0.20}(\mathrm{exp\leavevmode\nobreak\ syst.})^{+0.32}_{-0.16}(\mathrm{theo\leavevmode\nobreak\ syst.})=1.0^{+1.4}_{-0.9}$. The observed (expected) significances of $WH$ and $ZH$ production modes are 2.6 (1.3) standard deviations and 2.8 (1.2) standard deviations above the SM background, including other Higgs-boson processes, respectively. When determining the significance for $WH$ production, the $ZH$ signal-strength parameter is left floating in the fit, and vice versa. The combination of the $WH$ and $ZH$ channels leads to an observed (expected) significance for the combined $VH$ production mode of 4.1 (1.9) standard deviations above the SM background, including other Higgs-boson processes. The $p$-value with respect to the value predicted by the SM corresponds to about two standard deviations. The validity of the asymptotic approximation used in deriving these results was tested using pseudo-experiments. The combined signal strength for $VH$ is measured to be

[TABLE]

with an expected signal strength of $\mu_{VH}=1.0^{+0.66}_{-0.54}(\mathrm{stat.})^{+0.26}_{-0.17}(\mathrm{exp\leavevmode\nobreak\ syst.})^{+0.22}_{-0.17}(\mathrm{theo\leavevmode\nobreak\ syst.})=1.0^{+0.7}_{-0.6}$.

The cross-section times branching-fraction values, $\sigma_{WH}\cdot\mathcal{B}_{H\rightarrow WW^{*}}$ and $\sigma_{ZH}\cdot\mathcal{B}_{H\rightarrow WW^{*}}$, are simultaneously determined to be:

[TABLE]

The main contributions to the uncertainties in $\sigma_{WH}\cdot\mathcal{B}_{H\rightarrow WW^{*}}$ and $\sigma_{ZH}\cdot\mathcal{B}_{H\rightarrow WW^{*}}$ are summarised in Table 6. The predicted cross-section times branching-fraction values are $0.293\pm 0.007$ pb and $0.189\pm 0.007$ pb for $WH$ and $ZH$ [72], respectively. The 68% and 95% confidence level two-dimensional contours of $\sigma_{WH}\cdot\mathcal{B}_{H\to WW^{\ast}}$ and $\sigma_{ZH}\cdot\mathcal{B}_{H\to WW^{\ast}}$ are shown in Figure 8. The corresponding one-dimensional results, where the other parameter is left floating in the fit, are consistent with the SM predictions within $1.3\sigma$ for $WH$ and $1.5\sigma$ for $ZH$.

9 Conclusion

Measurements of the inclusive production of a Higgs boson in association with a $W$ or $Z$ boson are presented. Results of the $WH$ and $ZH$ production cross section times branching fraction for $H\rightarrow WW^{\ast}$ are obtained selecting the leptonic decays $W\rightarrow\ell\nu$ and $Z\rightarrow\ell\ell$, with $\ell=e,\ \mu$, in $\sqrt{s}=13\,\mathrm{TeV}$ proton–proton collisions recorded with the ATLAS detector at the LHC in a data sample corresponding to an integrated luminosity of $36.1\,\mathrm{fb}^{-1}$. The products of the $H\rightarrow WW^{*}$ branching fraction times the $WH$ and $ZH$ cross sections are measured to be $0.67^{+0.31}_{-0.27}$(stat.)${}^{+0.18}_{-0.14}$(syst.) pb and $0.54^{+0.31}_{-0.24}$(stat.)${}^{+0.15}_{-0.07}$(syst.) pb, respectively, compatible with the Standard Model predictions of $0.293\pm 0.007$ pb for $WH$ and $0.189\pm 0.007$ pb for $ZH$ within up to $1.5\sigma$. The $WH$ channel with $H\rightarrow WW^{\ast}$ is purely sensitive to the $HW$ coupling and this measurement is to date the most precise measurement of this channel.

Acknowledgements

We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently.

We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia; BMWFW and FWF, Austria; ANAS, Azerbaijan; SSTC, Belarus; CNPq and FAPESP, Brazil; NSERC, NRC and CFI, Canada; CERN; CONICYT, Chile; CAS, MOST and NSFC, China; COLCIENCIAS, Colombia; MSMT CR, MPO CR and VSC CR, Czech Republic; DNRF and DNSRC, Denmark; IN2P3-CNRS, CEA-DRF/IRFU, France; SRNSFG, Georgia; BMBF, HGF, and MPG, Germany; GSRT, Greece; RGC, Hong Kong SAR, China; ISF and Benoziyo Center, Israel; INFN, Italy; MEXT and JSPS, Japan; CNRST, Morocco; NWO, Netherlands; RCN, Norway; MNiSW and NCN, Poland; FCT, Portugal; MNE/IFA, Romania; MES of Russia and NRC KI, Russian Federation; JINR; MESTD, Serbia; MSSR, Slovakia; ARRS and MIZŠ, Slovenia; DST/NRF, South Africa; MINECO, Spain; SRC and Wallenberg Foundation, Sweden; SERI, SNSF and Cantons of Bern and Geneva, Switzerland; MOST, Taiwan; TAEK, Turkey; STFC, United Kingdom; DOE and NSF, United States of America. In addition, individual groups and members have received support from BCKDF, CANARIE, CRC and Compute Canada, Canada; COST, ERC, ERDF, Horizon 2020, and Marie Skłodowska-Curie Actions, European Union; Investissements d’ Avenir Labex and Idex, ANR, France; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by EU-ESF and the Greek NSRF, Greece; BSF-NSF and GIF, Israel; CERCA Programme Generalitat de Catalunya, Spain; The Royal Society and Leverhulme Trust, United Kingdom.

The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN, the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA), the Tier-2 facilities worldwide and large non-WLCG resource providers. Major contributors of computing resources are listed in Ref. [73].