Constraints on off-shell Higgs boson production and the Higgs boson total width in ZZ -> 4l and ZZ -> 2l2v final states with the ATLAS detector

Collaboration A measurement of off-shell Higgs boson production in the Z Z → 4 (cid:2) and Z Z → 2 (cid:2) 2 ν decay channels, where (cid:2) stands for either an electron or a muon, is performed using data from proton–proton collisions at a centre-of-mass energy of √ s = 13 TeV. The data were collected by the ATLAS experiment in 2015 and 2016 at the Large Hadron Collider, and they correspond to an integrated luminosity of 36 . 1 fb − 1 . An observed (expected) upper limit on the off-shell Higgs signal strength, deﬁned as the event yield normalised to the Standard Model prediction, of 3.8 (3.4) is obtained at 95% conﬁdence level (CL). Assuming the ratio of the Higgs boson couplings to the Standard Model predictions is independent of the momentum transfer of the Higgs production mechanism considered in the analysis, a combination with the on-shell signal-strength measurements yields an observed (expected) 95% CL upper limit on the Higgs boson total width of 14.4 (15.2) MeV. © 2018 The Author. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP 3 .


Introduction
The observation of the Higgs boson by the ATLAS and CMS experiments [1,2] at the Large Hadron Collider (LHC) marks a milestone towards the understanding of the mechanism of electroweak (EW) symmetry breaking [3][4][5]. Further studies of the spin, parity and couplings of the new particle have shown no significant deviation from the predictions for the Standard Model (SM) Higgs boson [6][7][8][9][10]. Efforts to measure the properties of the Higgs boson are primarily focused on on-shell production. For a Higgs boson at a mass of 125 GeV [10,11], the expected natural width of the SM Higgs boson is SM H ∼ 4.1 MeV [12]. However, above 125 GeV offshell production of the Higgs boson has a substantial cross section at the LHC [13][14][15][16], due to the increased phase space as the vector bosons (V = W , Z ) and top quark decay products become on-shell with the increasing energy scale. This provides an opportunity to study the Higgs boson properties at higher energy scales. Off-shell production can provide sensitivity to new physics that alters the interactions between the Higgs boson and other fundamental particles in the high-mass region [17][18][19][20][21][22][23][24].
The measured off-shell event yield from gluon-gluon fusion (ggF) production normalised to the SM prediction, where this E-mail address: atlas .publications @cern .ch. ratio is referred to as the signal strength μ off-shell , can be expressed as μ off-shell = σ gg→H * →Z Z off-shell σ gg→H * →Z Z off-shell,SM = κ 2 g,off-shell · κ 2 Z ,off-shell , where σ gg→H * →Z Z off-shell is the cross section of the off-shell Higgs boson production via ggF with subsequent decay into a Z Z pair, and κ g,off-shell and κ Z ,off-shell are the off-shell coupling modifiers relative to the SM predictions associated with the gg → H * production and the H * → Z Z decay, respectively. The off-shell Higgs boson signal cannot be treated independently of the gg → Z Z background, as sizeable negative interference effects appear [13]. The interference term is assumed to be proportional to √ μ off-shell = κ g,off-shell · κ Z ,off-shell . Similarly, μ on-shell for the on-shell Higgs boson production via ggF is given by: μ on-shell = σ gg→H→ Z Z * on-shell σ gg→H→ Z Z * on-shell,SM = κ 2 g,on-shell · κ 2 fiers [15,25]. The above formalism describing the ratio of off-shell to on-shell cross sections also applies to the vector-boson fusion (VBF) production mode. As in the previous measurement [26], for a measurement of H it is necessary to assume that the onshell and off-shell coupling modifiers are the same, and for an upper limit that the on-shell coupling modifiers are not larger than the off-shell couplings. It is also assumed that any new physics which modifies the off-shell signal strength and the off-shell couplings does not modify the relative phase of the interfering signal and background processes. Further, it is assumed that there are neither sizeable kinematic modifications to the off-shell signal nor new sizeable signals in the search region of this analysis unrelated to an enhanced off-shell signal strength.
The ATLAS and CMS experiments have presented studies of the off-shell production of the Higgs boson using Run-1 proton-proton (pp) collisions data [26][27][28][29]. ATLAS obtained an observed (expected) upper limit on the off-shell Higgs boson signal strength (μ off-shell ) in the range of 5.1-8.6 (6.7-11.0) [26], using the Z Z and W W channels. This range is determined by the assumption that the gg → Z Z and gg → W W background K -factors, corresponding to the ratio of the next-to-leading-order (NLO) QCD predictions to the leading-order (LO) predictions, lie between one-half and twice the value of the gg → H * → Z Z(W W ) signal K -factor.
An observed (expected) 95% confidence level (CL) upper limit of H < 23(33) MeV was obtained, assuming the gg → Z Z(W W ) background K -factor is equal to the gg → H * → Z Z(W W ) signal K -factor. CMS presented a similar study in the Z Z and W W channels, with observed (expected) 95% CL upper limit of H < 13(26) MeV [29]. By comparison, the precision of H from direct on-shell Higgs boson mass measurements alone is approximately 1 GeV [9, 30,31], limited by measurement resolution.
This Letter presents an analysis of off-shell Higgs boson production in the Z Z → 4 and Z Z → 2 2ν final states ( = e, μ), using 36.1 fb −1 of data collected by the ATLAS detector in pp collisions at √ s = 13 TeV. The off-shell region is defined by requiring the invariant mass of the Z Z system (m Z Z ) to be above the on-shell Z Z production threshold, hence well above the Higgs boson mass, and the on-shell region is defined by a mass window around the 125 GeV resonance. This analysis adopts the same methodology used in the Run-1 analysis reported in Ref. [26]. The analysis for the Z Z → 4 final state closely follows the Higgs boson measurements and high-mass search in the same final state described in Refs. [32,33]. The off-shell Higgs signal strength is extracted using a matrix-element discriminant, defined in Section 4, in a mass region 220 GeV < m 4 < 2000 GeV. The on-shell signal strength was measured in the 118 GeV < m 4 < 129 GeV region in Ref. [32]. The analysis of the Z Z → 2 2ν channel, described in Section 5, follows a strategy similar to that used in the search for heavy Z Z resonances described in Ref. [33]. For this channel, the signal strength is extracted from the transverse mass distribution in the 250 to 2000 GeV range. For off-shell production of the Higgs boson, the dominant processes of ggF and VBF are considered. Next-to-nextto-leading-order (NNLO) QCD and NLO EW corrections are known for the off-shell signal process gg → H * → Z Z [25]. More recently, NLO QCD corrections have also become available for the gg → Z Z background and for the signal-background interference [34,35], for which additional details are given in Section 3. Given that the QCD corrections for the off-shell signal processes have only been calculated inclusively in the jet multiplicity, the analysis is performed inclusively in jet observables and the event selection is designed to minimise the dependence on the momentum of the Z Z system, which is sensitive to the jet multiplicity.

ATLAS detector
The ATLAS experiment is described in Ref. [36]. ATLAS is a multipurpose detector with a forward-backward symmetric cylindrical geometry and a solid-angle 1 coverage of nearly 4π . The inner tracking detector, covering the region |η| < 2.5, consists of a silicon pixel detector, a silicon microstrip detector and a straw-tube transition-radiation tracker. The innermost layer of the pixel detector, the insertable B-layer [37], was installed between Run 1 and Run 2 of the LHC. The inner detector is surrounded by a thin superconducting solenoid providing a 2 T magnetic field, and by a finely segmented lead/liquid-argon (LAr) electromagnetic calorimeter covering the region |η| < 3.2. A steel/scintillator-tile hadronic calorimeter provides coverage in the central region |η| < 1.7. The endcap and forward regions, covering the pseudorapidity range 1.5 < |η| < 4.9, are instrumented with electromagnetic and hadronic LAr calorimeters, with copper or tungsten as the absorber material. A muon spectrometer system incorporating large superconducting toroidal air-core magnets surrounds the calorimeters. Three layers of precision wire chambers provide muon tracking in the range |η| < 2.7, while dedicated fast chambers are used for triggering in the region |η| < 2.4. The trigger system is composed of two stages [38]. The first stage, implemented with custom hardware, uses information from calorimeters and muon chambers to reduce the event rate from about 40 MHz to a maximum of 100 kHz. The second stage, called the high-level trigger, reduces the data acquisition rate to about 1 kHz on average. The high-level trigger is software-based and runs reconstruction algorithms similar to those used in the offline reconstruction.
NLO ratio of the QCD predictions over the data region considered in the analysis. Using these scaled NLO K -factors, the cross section for the gg → (H * →)Z Z process with any off-shell Higgs boson signal strength μ off-shell can be obtained from a parameterisation of three SM MC samples: the gg → H * → Z Z signal (σ SM gg→H * →Z Z ), the gg → Z Z continuum background (σ SM gg→Z Z, cont ) and the full process with signal, background and interference gg → (H * →)Z Z (σ SM gg→(H * →)Z Z ), where the last sample is required to derive the interference sample: (2) The electroweak pp → V V + 2 j processes containing both the VBF-like events and events from associated Higgs production with vector bosons (V H), which includes on-shell Higgs boson production, were simulated using MadGraph5_aMC@NLO [45] with matrix elements calculated at LO. The QCD renormalisation and factorisation scales were set to m W following the recommendation in Ref. [46] and the NNPDF23LO PDF set [47] was used.
Pythia 8.186 [48] was used for parton showering and hadronisation, with the A14 set of tuned parameters for the underlying event [49]. Due to the different H dependence, the on-shell and off-shell Higgs boson production processes are separated when weighting MC events as in Eqs.
(1) by requiring that the gener- The qq → Z Z background was simulated with Sherpa v2.2.2, using the NNPDF30NNLO PDF set for the hard-scattering process. NLO QCD accuracy is achieved in the matrix-element calculation for 0-and 1-jet final states and LO accuracy for 2-and 3-jet final states. The merging with the Sherpa parton shower was performed using the MePs@NLO prescription. NLO EW corrections are applied as a function of the particle-level m Z Z [50,51]. The W W and W Z backgrounds were simulated at NLO in QCD using the Powheg-Box v2 event generator [52] with the CT10NLO PDF set [53] and Pythia 8.186 for parton showering and hadronisation. The non-perturbative effects were modelled with the AZNLO set of tuned parameters [54]. The interference between the qq → Z Z and qq → W W processes for the 2 2ν final state is found to be negligible and thus is not considered.
Events containing a single Z boson with associated jets ( Z + jets) were simulated using the Sherpav v2.2.1 event generator. Matrix elements were calculated for up to two partons at NLO and four partons at LO using the Comix [55] and OpenLoops [41] matrix-element generators and merged with the Sherpa parton shower [43] using the MePs@NLO prescription. The NNPDF30NNLO PDF set was used in conjunction with dedicated parton-shower tuning developed by the Sherpa authors. The Z + jets events are normalised using the NNLO cross sections [56].
The triboson backgrounds Z Z Z , W Z Z, and W W Z with fully leptonic decays and at least four prompt charged leptons were modelled using Sherpav v2.2.1. The contribution from triboson backgrounds with one W or Z boson decaying hadronically is not included in the simulation, but the impact on the analysis is found to be negligible. For the fully leptonic tt + Z background, with four prompt charged leptons originating from the decays of the top quarks and Z boson, MadGraph5_aMC@NLO was used. The tt background, as well as the single-top and W t production, were modelled using Powheg-Box v2 interfaced to Pythia 6.428 [57] with the Perugia 2012 [58] set of tuned parameters for parton showering, hadronisation and the underlying event, and to Evt-Gen v1.2.0 [59] for properties of the bottom and charm hadron decays.
The particle-level events produced by each MC event generator were processed through the ATLAS detector simulation within the Geant 4 framework [60,61] or the fast detector simulation package Atlfast-II [61]. Additional pp interactions in the same and nearby bunch crossings (pile-up) are included in the simulation. The pileup events were generated using Pythia 8 with the A2 set of tuned parameters [62] and the MSTW2008LO PDF set [63]. The simulation samples were weighted to reproduce the observed distribution of the mean number of interactions per bunch crossing in the data.

Z Z → 4 analysis
The analysis for the Z Z → 4 final state closely follows the on-shell Higgs boson measurements and high-mass search in the same final state described in Refs. [32,33], with the same event reconstruction, trigger and event selections, and background estimation methods. A matrix-element-based (ME-based) discriminant computed at LO is constructed to enhance the separation between the gg → H * → Z Z signal and the gg → Z Z and qq → Z Z backgrounds, and this discriminant is subsequently used in a binned maximum-likelihood fit for the final result. To minimise the dependence of the gg → Z Z kinematics on higher-order QCD effects, the analysis is performed inclusively, ignoring the number of jets in the events.
The analysis is split into three channels (4μ, 2e2μ, 4e). Each electron (muon) must have transverse momentum p T > 7 (5) GeV and be measured in the pseudorapidity range |η| < 2.47 (|η| < 2.7). The highest-p T lepton in the quadruplet must satisfy p T > 20 GeV, and the second (third) lepton in p T order is required to have p T > 15 GeV (p T > 10 GeV). Lepton pairs are formed from same-flavour opposite-charge leptons. For each channel, the quadruplet with a lepton pair whose mass is closest to the Z boson mass is kept. This pair is referred to as the leading dilepton pair and its invariant mass, m 12 , is required to be between 50 GeV and 106 GeV. The second (subleading) pair is chosen from the remaining leptons as the pair closest in mass to the Z boson and in the range 50 GeV < m 34 < 115 GeV. The off-shell region is defined as the range 220 GeV < m 4 < 2000 GeV, while the on-shell region is defined as 118 GeV < m 4 < 129 GeV.
The dominant background in the Z Z → 4 channel arises from qq → Z Z events. This is modelled using MC simulation, accurate to NLO QCD and NLO EW corrections as explained in Section 3. Other backgrounds, such as triboson production, tt V , Z + jets, and top quark production, constitute less than 2% of the total background in the off-shell signal region, and are either taken from simulation or from dedicated data control regions. Fig. 1(a) shows the observed and expected distributions of m 4 combining all lepton channels in the off-shell region. The data are in agreement with the SM predictions, with two small excesses at m 4 around 240 GeV and 700 GeV, each having a significance of about two standard deviations (2σ ), as evaluated by the high-mass resonance search reported in Ref. [33]. Table 1 shows the expected and observed numbers of events in the signal region and additionally in the 400 GeV < m 4 < 2000 GeV mass range, which is backgrounds and the Higgs boson with μ off-shell = 5. The ratio plot shows the observed data yield divided by the SM prediction (black points) as well as the total expected event yield with μ off-shell = 5 divided by the SM prediction (dashed line) in each bin.

Table 1
The expected and observed numbers of events in the signal region for both final states. For the Z Z → 4 analysis, numbers are given for both the signal region and a signal-enriched region which covers the mass range 400 GeV < m 4 < 2000 GeV. The other backgrounds in the Z Z → 4 final state include contributions from Z + jets and top quark processes, while in the Z Z → 2 2ν final state they include contributions from tri-boson production, the W + jets process, and top quark processes other than pair production. For the Z Z → 2 2ν analysis, the range 250 GeV < m Z Z T < 2000 GeV is considered. The upper part of the table contains the expected events for the gg → (H * →)Z Z and VBF (H * →)Z Z processes which include the Higgs boson signal, background and interference for the SM predictions. The SM estimates for the signal (S) and background (B) event yields without interference are given in parentheses. The lower part of the table contains the corresponding predictions for μ off-shell = 5. The uncertainties in the number of expected events include the statistical uncertainties from MC samples and systematic uncertainties, summed in quadrature. Empty entries correspond to contributions with event yields smaller than 0.1 events.
1.67 ± 0.08 1.14 ± 0.04 5.45 ± 0.30) Other signal hypothesis 11.0 ± 0.5 4 .85 ± 0.22 8.8 ± 0.4 enriched in signal. The latter mass region was chosen for this table since it is optimal for a counting experiment. The matrix-element kinematic discriminant fully exploits the event kinematics in the centre-of-mass frame of the 4 system. It is computed from eight kinematic observables: the three masses m 4 , m 12 and m 34 , and the leading Z boson production angle and four decay angles defined in Ref. [64]. These observables are used to calculate the matrix elements for the different processes with the MCFM program [15] at LO. The following matrix elements are calculated for each event in the mass range 220 GeV < m 4 < 2000 GeV: • P qq : the matrix element squared for the qq → Z Z → 4 process, • P gg : the matrix element squared for the gg → (H * →)Z Z → 4 process, which includes the Higgs boson with SM couplings, the continuum background and their interference, • P H : the matrix element squared for the gg → H * → Z Z → 4 process without continuum background or interference.
The ME-based discriminant is defined as in Ref. [15]: where c = 0.1 is a constant whose value is chosen to balance the overall cross sections of the qq → Z Z and gg → (H * →)Z Z processes. The value of c has a small effect on the analysis sensitivity. Fig. 1(b) shows the observed and expected distributions of D ME .
Events with a D ME value between −4.5 and 0.5 are used for the final result.

Z Z → 2 2ν analysis
The analysis in the Z Z → 2 2ν final state closely follows the one performed to search for Z Z resonances [33]. The reconstruction, identification and selection of electrons, muons, jets, b-jets and missing transverse momentum are identical while the event selection is optimised for the current analysis.
To discriminate the signal from the background and enhance the sensitivity to off-shell Higgs boson production, the transverse mass of the Z Z system (m Z Z T ) is used, defined as: where p T is the transverse momentum of the dilepton system, m Z is the mass of the Z boson fixed to m Z = 91.187 GeV [65] and E miss T is the magnitude of the missing transverse momentum E miss T . The latter is computed as the negative sum of transverse momenta of all the leptons and jets, as well as the tracks originating from the primary vertex but not associated with any of the leptons or jets, the so-called soft term.
The event selection is designed to minimise the dependence on the p T of the Z Z system, and thus is performed inclusively in number of jets. First, events with two opposite-charge leptons of the same flavour are selected with the requirement of p T > 30 (20) GeV for the leading (sub-leading) lepton. The dilepton invariant mass m is required to be in the range 76 GeV < m < 106 GeV. Events with additional, loosely identified leptons with p T > 7 GeV are rejected to reduce the amount of W Z background. The two Z bosons originating from the decay of an off-shell Higgs boson are boosted and tend to be back-toback in the transverse plane. A series of selection requirements are applied to reduce the Z + jets background: the two leptons are required to be produced with an angular separation of R ≡ ( η) 2 + ( φ) 2 < 1.8; E miss T is required to be larger than 175 GeV; the azimuthal angle between the transverse momentum of the dilepton system and the missing transverse momentum is required to be large, φ( p T , E miss T ) > 2.7; p T is required to be The dominant backgrounds in the Z Z → 2 2ν channel consist of qq → Z Z events, followed by W Z events. Other background processes with two genuine leptons not directly originating from a Z boson decay include W W , tt, W t and Z → τ τ . The remaining background comes from Z + jets with poorly reconstructed E miss T , W +jets events with at least one misidentified electron or muon, semileptonic top decays, and multi-jet events.
The qq → Z Z background is modelled in the same manner as for the Z Z → 4 channel. The W Z background is estimated with simulation using a normalisation correction factor extracted from a dedicated control region (CR). This W Z-enriched CR is defined by selecting Z → candidates with an additional electron or muon with p T > 20 GeV. Events with a b-jet are rejected to suppress leptonic tt decays and a m T (W ) > 60 GeV requirement is applied to reduce the Z +jets contamination. The correction factor is then calculated in the CR as the number of data events, after subtracting the non-W Z contributions, divided by the predicted W Z yield, and is found to be 1. 29. The statistical uncertainty of the W Z estimate is about 2%, while the systematic uncertainty is estimated to be 5% from theoretical and experimental uncertainties in the simulation-based transfer factor between the three-lepton control region and the two-lepton signal region.
The selection criteria from MC simulation of the non-resonant-background. The m Z Z T distributions for the non-resonant-background are derived from the data CR and extrapolated to the SR. The total uncertainty in the non-resonant-estimate is about 40%, including the statistical uncertainty of the data in the control region, the extrapolation and the method bias estimated from simulation. The m Z Z T distribution differences between data and simulation are taken as a shape uncertainty (∼10%).
The Z + jets background, expected to be ∼2% of the total background, is estimated from a combination of MC and data-driven techniques. A Z + jets enriched CR is defined by reversing the E miss T /H T selection. Additionally, the b-jet veto and the requirement on φ( p T , E miss T ) are removed to allow more data events.
The estimation is performed by extrapolating the number of events observed in the CR, after subtracting non-Z -boson backgrounds, to the SR with a correction factor based on simulation. The m Z Z T distribution for the Z + jets background is derived from simulation. The total uncertainty in the Z + jets estimate is about 50% (80%) for the ee (μμ) channel, including the statistical uncertainty of the data in the control region and the extrapolation factor. The shape difference in m Z Z T between Z + jets MC events in the SR and those in the CR is taken into account as a systematic uncertainty.
Other backgrounds, such as triboson production, tt V , W + jets, and top quark processes other than pair production, constitute only a tiny fraction of the total background in the off-shell signal region, < 1%, and are taken from simulation. The contribution from the on-shell Higgs production is negligible in the off-shell signal region.
The expected and observed numbers of events in the signal region for the Z Z → 2 2ν analysis are summarised in Table 1. Fig. 2 shows the observed and expected distributions of m Z Z T in both the ee and μμ channels in the off-shell region.

Systematic uncertainties
Systematic uncertainty sources impacting the analysis of both channels can be divided into two categories: uncertainties in the theoretical description of the signal and background processes and experimental uncertainties related to the detector or to the reconstruction algorithms. The largest systematic uncertainties arise from the theoretical uncertainties in the gg-initiated Z Z processes and the qq → Z Z background process. The uncertainties from experimental measurements are generally small compared to the theoretical uncertainties in this analysis. The theoretical uncertainties originate from the PDF choice, the missing higher-order corrections, and the parton-shower modelling. The PDF uncertainty corresponds to the 68% CL variations of the nominal PDF set NNPDF30NNLO for both qq → Z Z and gg → (H * →)Z Z , as well as the difference from alternative PDF sets. The alternative PDF sets used are CT10NNLO [68] and MMHT2014NNLO [69]. The uncertainty due to PDF is found to be about 3% in the high-mass region considered. The uncertainty due to higher-order QCD corrections (QCD scale uncertainty) is estimated by varying the renormalisation and factorisation scales independently, ranging from a factor of one-half to two. The uncertainty in the K -factors due to the NLO QCD scale uncertainty is 10-20% as a function of m Z Z for the gg-initiated Z Z processes in the probed high-mass region, and ranges from 5% to 10% as a function of m Z Z for the qq → Z Z background. The QCD scale uncertainties are treated as correlated among the gg-initiated Z Z processes, and uncorrelated with the qq-initiated Z Z process. There are a few additional normalisation uncertainties associated with the NLO K -factors discussed in section 3. In the region below 2m t , the higher-order corrections are computed with a maximum jet transverse momentum of 150 GeV to ensure a good description by the 1/m t expansion. The default scale uncertainty is therefore doubled for events which have a jet with p T > 150 GeV, corresponding to about 8% of the events in this region. The scale uncertainty is also increased by 50% around the 2m t threshold, with a Gaussiansmoothed transition decreasing to the default uncertainty within 50 GeV of the threshold. This is intended to allow for possible effects on the K -factor which have not been estimated as the top quark moves on-shell. It is assumed that the 10-20% NLO QCD scale uncertainty for the gg-initiated Z Z processes covers the assumption of massless loops above the 2m t threshold, and as well the uncertainties in the 1.2 scale factor estimated only for the NNLO/NLO signal correction but also applied to the background and interference components. These NLO QCD scale uncertainties are larger than those associated with the NNLO QCD signal uncertainties. The EW correction uncertainty for qq → Z Z is evaluated using the same method as in Ref. [26] and its impact is estimated to be about 1%. The parton-shower uncertainty is evaluated by varying parameters in the parton-shower tunes according to Refs. [49,54] and found to be 2-3% in normalisation.
The theoretical uncertainties due to the missing higher-order corrections and PDF variations are small for V H-like and VBF-like processes pp → Z Z + 2 j; therefore, they are not included in the analysis.
For the Z Z → 4 analysis, the same sources of experimental uncertainty as in Ref.
[32] are evaluated. The leading experimental systematic uncertainties are due to the electron and muon reconstruction and selection efficiency uncertainties, which are smaller than the uncertainties associated with the theoretical predictions.
Similarly, for the Z Z → 2 2ν channel, the same sources of experimental uncertainty as in Ref. [70] are evaluated. These experimental uncertainties affect the sensitivity of the μ off-shell measurement only at the percent level.
The uncertainty in the combined 2015 and 2016 integrated luminosity is 2.1%, derived following a methodology similar to that detailed in Ref.
[71], from a preliminary calibration of the luminosity scale using x-y beam-separation scans. This uncertainty is applied to the normalisation of the signal and also to background contributions whose normalisations are derived from MC simulations. A variation in the pile-up reweighting of MC events is included to cover the uncertainty in the ratio of the predicted and measured inelastic cross sections in Ref. [72].

Results
The results for the Z Z → 4 and Z Z → 2 2ν analyses are first translated into limits on the off-shell signal strength μ off-shell .
A single off-shell signal-strength parameter is applied for all production modes, assuming that the ratio of the off-shell production rates via the ggF process to those via the VBF process are as predicted in the SM, namely μ ggF off-shell /μ VBF off-shell = 1. In a second step, the off-shell analyses are combined with the on-shell Z Z * → −0. 19 . The combination with the onshell analysis is performed with two assumptions that correspond to different interpretations of the results. In the first combination, the parameter of interest is the ratio of off-shell to on-shell signal strengths, which can be interpreted as the Higgs boson width normalised to its SM prediction: μ off-shell /μ on-shell = H / SM H . This interpretation assumes that the off-and on-shell coupling modifiers are the same for both ggF and VBF production modes (i.e., κ g,on-shell = κ g,off-shell = κ V ,on-shell = κ V ,off-shell ). In the second combination, the parameter of interest is the ratio of off-shell to on-shell signal strengths for the ggF production only, R gg = μ ggF off-shell /μ ggF on-shell , which can be interpreted as the ratio of off-shell to on-shell gluon couplings: R gg = κ 2 g,off-shell /κ 2 g,on-shell . In this case the coupling scale factors κ V = κ V ,on-shell = κ V ,off-shell associated with on-and off-shell VBF production and the H ( * ) → Z Z decay are assumed to be the same and fitted to the data (profiled). This also assumes that the total width is equal to the SM prediction.
The statistical analysis is based on the framework described in Refs. [74][75][76]. A binned likelihood function is constructed as a product of Poisson probability terms over all bins of the fit templates considered. This function depends on the parameter of interest μ, corresponding to one of the different interpretations discussed above (μ off-shell , H / SM H and R gg ), and θ , a set of nuisance parameters that encode the effects of systematic uncertainties on the signal and expected backgrounds, as described in Section 6. The nuisance parameters are constrained using either Gaussian or log-normal terms.
In the Z Z → 4 channel, a binned maximum-likelihood fit to the D ME distribution is performed to extract the limits on μ. The The PDF uncertainties and uncertainties from higher-order QCD corrections applied to the qq → Z Z process are considered correlated between the on-shell and off-shell measurements. Given the different theoretical computations, the corresponding uncertainties are considered uncorrelated for the gg-initiated Z Z processes between the on-shell and off-shell measurements, and the impact of such a correlation effect is found to be small. In addition to the main theoretical uncertainties, the common experimental system-  Table 2. As a result of the small data excess observed in the off-shell region, the observed limits on μ off-shell are less stringent than the expected ones. The observed (expected) limit on H / SM H is 3.5 (3.7) at the 95% CL. Due to the fact that the measured on-shell signal strength μ on-shell is larger than one [32], the observed limit on H / SM H is smaller than the expected limit. The limit on H / SM H can be translated into a limit on the total width of the Higgs boson, leading to an observed (expected) 95% CL upper limit on the Higgs boson total width of 14.4 (15.2) MeV.
These results are significantly improved compared to the Run-1 publication [26], the expected limit being about a factor two better.
If instead of constraining the qq → Z Z background to the theoretical expectation, the normalisation is left as a free parameter in the profile likelihood fit, the upper limits on μ on-shell are about 4% worse in the Z Z → 4 channel. If only the NLO K -factor are applied to the SM prediction of the gg-initiated Z Z processes, without the additional NNLO/NLO K -factor of 1.2 (Section 3), the upper limits on μ off-shell and H / SM H are about 10% worse. The impact of the various systematic uncertainties on the expected limit in the μ off-shell fit are listed in Table 3. The values in this table were derived by fixing all the nuisance parameters associated with the systematic uncertainties to the values derived from the SM-conditional fit to the data, with the exception of the one under study. The uncertainties with the largest impact on the sensitivity of μ off-shell are the theoretical uncertainties of the gg-and qq-initiated Z Z processes.

Conclusion
A determination of the off-shell Higgs boson signal strength in the Z Z → 4 and Z Z → 2 2ν final states and their combination is  Table 3 The expected 95% CL upper limit on μ off-shell with a ranked listing of the impact of the leading systematic uncertainty individually, comparing with no systematic uncertainty or all systematic uncertainties. The upper limits are evaluated using the CL s method. Assuming that the total width of the Higgs boson is as expected in the SM, and the coupling scale factors associated with on-and off-shell VBF production and the H ( * ) → Z Z decay are the same, the same combination can be interpreted as a limit on the ratio of the off-shell to the on-shell couplings to gluons R gg = κ 2 g,off-shell /κ 2 g,on-shell . An observed (expected) limit of 4.3 (4.1) at 95% CL on R gg is obtained.