Measurements of the Higgs boson inclusive and differential fiducial cross sections in the 4$\ell$ decay channel at $\sqrt{s}$ = 13 TeV

Inclusive and differential fiducial cross sections of the Higgs boson are measured in the $H \to ZZ^{*} \to 4\ell$ ($\ell = e,\mu$) decay channel. The results are based on proton$-$proton collision data produced at the Large Hadron Collider at a centre-of-mass energy of 13 TeV and recorded by the ATLAS detector from 2015 to 2018, equivalent to an integrated luminosity of 139 fb$^{-1}$. The inclusive fiducial cross section for the $H \to ZZ^{*} \to 4\ell$ process is measured to be $\sigma_\mathrm{fid} = 3.28 \pm 0.32$ fb, in agreement with the Standard Model prediction of $\sigma_\mathrm{fid, SM} = 3.41 \pm 0.18 $ fb. Differential fiducial cross sections are measured for a variety of observables which are sensitive to the production and decay of the Higgs boson. All measurements are in agreement with the Standard Model predictions. The results are used to constrain anomalous Higgs boson interactions with Standard Model particles.


Introduction
The observation of the Higgs boson by the ATLAS and CMS Collaborations [1,2] using data from protonproton (pp) collisions at the Large Hadron Collider (LHC) recorded in 2011 and 2012 at centre-of-mass energies of √ s = 7 TeV and 8 TeV, respectively, was a major step forward in the understanding of the electroweak (EW) symmetry breaking mechanism [3][4][5]. Studies of the spin and parity of the Higgs boson, its coupling structure to other particles, and measurements of fiducial and differential cross sections have been performed . These show no significant deviations from the Standard Model (SM) predictions for the Higgs boson with a mass of 125.09 ± 0. 24 GeV [29].
This paper presents updated inclusive and differential cross-section measurements of the Higgs boson in the H → Z Z * → 4 decay channel (where = e or µ). The full ATLAS Run 2 dataset, consisting of pp collision data at √ s = 13 TeV taken between 2015 and 2018, is used for this analysis. The total integrated luminosity after imposing data quality requirements is 139 fb −1 , with a data-taking efficiency of 91.5%.
All measurements are performed with the assumption that the mass of the Higgs boson is 125 GeV, and are compared with SM predictions. The signal is extracted from a binned likelihood fit to the four-lepton invariant mass, m 4 , distribution. All major background processes are estimated from data. In particular, the normalisation of the dominant non-resonant Z Z * background is now constrained from dedicated data sidebands rather than from simulation. Signal events are corrected for detector measurement inefficiency and resolution by unfolding using the detector response matrix in the likelihood fit, in place of a bin-by-bin correction. Compared with the previous published results [11], this paper also benefits from the full LHC Run 2 integrated luminosity, improved event and electron reconstruction [30,31], and improved lepton isolation to mitigate the impact of additional pp interactions in the same or neighbouring bunch crossing (pile-up). The fiducial phase-space definition has also been updated with respect to the previous publication to harmonise the selection of the leptons.
The paper is organised as follows. A brief introduction of the ATLAS detector is given in Section 2, while in Section 3, the data and simulated signal and background samples are described. The selection of the Higgs boson candidate events is detailed in Section 4. Section 5 outlines the fiducial phase-space definition and the observables that are unfolded, while the background modelling is described in Section 6. The unfolding strategy is described in Section 7. The experimental and theoretical systematic uncertainties, detailed in Section 8, are taken into account for the statistical interpretation of the data. The final results are presented in Section 9 and their interpretation to constrain possible beyond the SM (BSM) contact interactions or non-SM values of the band c-quark Yukawa couplings are shown in Section 10. Concluding remarks are given in Section 11. More information about general aspects of the analysis is contained in the concurrent Ref. [32], where, in particular, details of the event selection and background estimation can be found. the identification of electrons. Lead/liquid-argon (LAr) sampling calorimeters provide electromagnetic energy measurements with high granularity. A steel/scintillator-tile hadron calorimeter covers the central pseudorapidity range (|η| < 1.7). The endcap and forward regions are instrumented up to |η| = 4.9 with LAr calorimeters for both the EM and hadronic energy measurements. The calorimeters are surrounded by the muon spectrometer, which has three large air-core toroidal superconducting magnets with eight coils each. The field integral of the toroid magnets ranges between 2.0 and 6.0 T m across most of the detector. The muon spectrometer includes a system of precision tracking chambers and fast detectors for triggering with a coverage of |η| < 2.7. Events are selected using a first-level trigger implemented in custom electronics, which reduces the event rate to a maximum of 100 kHz using a subset of detector information. Software algorithms with access to the full detector information are then used in the high-level trigger to yield a recorded event rate of about 1 kHz [36].

Theoretical predictions and event simulation
The production of the SM Higgs boson via gluon-gluon fusion (ggF), via vector-boson fusion (VBF), with an associated vector boson (VH, where V is a W or Z boson), and with a top quark pair (ttH) was modelled with the P -B v2 Monte Carlo (MC) event generator [37][38][39][40][41][42][43][44]. Table 1 summarises the predicted SM production cross sections and branching ratios for the H → Z Z * → 4 decay for m H = 125 GeV together with their theoretical accuracy. Table 1: Predicted SM Higgs boson production cross sections (σ) for ggF, VBF and five associated production modes in pp collisions for m H = 125 GeV at √ s = 13 TeV . For bbH the accuracy of calculations in the 4-and 5-flavour schemes (FS) is reported. The quoted uncertainties correspond to the total theoretical systematic uncertainties calculated by adding in quadrature the uncertainties due to missing higher-order corrections and PDF+α s . The decay branching ratios (B) with the associated uncertainty for H → Z Z * and H → Z Z * → 4 , with = e, µ, are also given. For ggF, the PDF4LHC next-to-next-to-leading-order (NNLO) set of parton distribution functions (PDF) was used, while for all other production modes, the PDF4LHC next-to-leading-order (NLO) set was used [72]. The simulation of ggF Higgs boson production used the P method for merging the NLO Higgs + jet cross section with the parton shower and the M NLO method [76] to simultaneously achieve NLO accuracy for the inclusive Higgs boson production. In a second step, a reweighting procedure (NNLOPS) [77], exploiting the Higgs boson rapidity distribution, was applied using the HNNLO program [78,79] to achieve NNLO accuracy in the strong coupling constant α s .
The matrix elements of the VBF, qq → VH and ttH production mechanisms were calculated to NLO accuracy in QCD. For VH production, the M NLO method was used to merge 0-and 1-jet events [44,76]. The gg → Z H contribution was modelled at leading order (LO) in QCD.
The production of a Higgs boson in association with a bottom quark pair (bbH) was simulated at NLO with M G 5_ MC@NLO v2.3.3 [80], using the CT10 NLO PDF [81]. The production in association with a single top quark (tH+X where X is either jb or W, defined in the following as tH) was simulated at NLO with M G 5_ MC@NLO v2.6.0 using the NNPDF30 PDF set [75].
For all production mechanisms the P 8 [82] generator was used for the H → Z Z * → 4 decay as well as for the parton shower modelling. The AZNLO set of tuned parameters [83] was used, except for ttH, where, like for the tt samples, the A14 tune [84] was employed. The event generator was interfaced to E G v1.2.0 [85] for simulation of the bottom and charm hadron decays. All signal samples were simulated for a Higgs boson mass m H = 125 GeV.
For additional cross-checks, the ggF sample was also generated with M G 5_ MC@NLO. This simulation has NLO QCD accuracy for zero, one and two additional partons merged with the F F merging scheme [86,87], and top and bottom quark mass effects are taken into account [88][89][90]. Some final results are also compared with ggF predictions calculated with R ISH, which provides resummation at N 3 LL+NNLO accuracy [91][92][93][94][95], and uses MATRIX for the fixed-order calculation [96,97]. Similarly, ggF predictions are also obtained from NNLOJET for distributions of Higgs plus one-or two-jet events [98][99][100]. Neither of these two predictions are included for the case in which there are zero jets. Additionally, final results for several of the variables that probe the kinematics of the Higgs boson decay products include comparisons with H 4 and P . These two programs include the full NLO electroweak corrections to the Higgs boson decay into four charged leptons [69][70][71][101][102][103][104][105][106].
The samples are normalised to cross sections obtained from the best available predictions as provided in Refs. [45-47, 67, 68, 73-75, 107]. The SM branching ratio prediction, taken from P 4 [69,102], includes the full NLO EW corrections, and interference effects which result in a branching ratio that is 10% higher for same-flavour final states (4µ and 4e) than for different-flavour states (2e2µ and 2µ2e).
For the BSM interpretation, described in Section 10.1, deviations from the SM are studied using a ggF sample generated with M G 5 at LO using F R [108] and the NNPDF23 PDF set. The sample was interfaced to P 8 using the A14 parameter set [84]. For studies of the Yukawa couplings described in Section 10.2, the gluon-initiated component of the prediction was calculated using R DISH, while M G 5_ MC@NLO was used for the quark-initiated component with F F merging for 0and 1-jet final states.
The Z Z * continuum background from quark-antiquark annihilation was modelled using S 2.2.2 [109][110][111], which provides a matrix element calculation accurate to NLO in α s for 0-and 1-jet final states, and LO accuracy for 2-and 3-jet final states. The merging with the S parton shower [112] was performed using the ME+PS@NLO prescription [113]. The NLO EW corrections were applied as a function of the invariant mass of the Z Z * system m Z Z * [114,115]. This process was also simulated using two additional MC generators. The first is P -B v2 interfaced to P 8 for parton showering and hadronisation, with E G for the simulation of bottom and charm hadron decays. The second is M G 5_ MC@NLO with F F merging at NLO for 0-and 1-jet final states and interfaced to P 8 for parton showering.
The gluon-induced Z Z * production was modelled by S 2.2.2 [109][110][111] at LO in QCD for 0-and 1-jet final states. The higher-order QCD effects for the gg → Z Z * continuum production have been calculated for massless quark loops [116][117][118] in the heavy top-quark approximation [119], including the gg → H * → Z Z processes [120,121]. The gg → Z Z simulation cross section is scaled by a K-factor of 1.7±1.0, defined as the ratio of the higher-order to leading-order cross section predictions. Production of Z Z * via vector-boson scattering was simulated at LO in QCD with the S 2.2.2 generator.
The WZ background was modelled using P -B v2 interfaced to P 8 and E G v1.2.0 for the simulation of bottom and charm hadron decays. The triboson backgrounds ZZZ, WZZ, and WWZ with four or more prompt leptons (denoted by VVV hereafter) were modelled using S 2.2.2. The simulation of tt + Z events with both top quarks decaying semileptonically and the Z boson decaying leptonically was performed with M G 5_ MC@NLO interfaced to P 8. The total cross section is normalised to the prediction of Ref. [63], which includes the two dominant terms at both the LO and the NLO in a mixed perturbative expansion in the QCD and EW couplings. For modelling comparisons, S 2.2.1 was used to simulate tt + Z events at LO. The smaller tW Z, ttW + W − , ttt, tttt and t Z background processes were simulated with M G 5_ MC@NLO interfaced to P 8.
The modelling of events containing Z bosons with associated jets (Z + jets) was performed using the S 2.2.1 generator. Matrix elements were calculated for up to two partons at NLO and four partons at LO using C [110] and O L [111], and merged with the S parton shower [112] using the ME+PS@NLO prescription [113]. The NNPDF3.0 NNLO PDF set was used in conjunction with a dedicated set of tuned parton shower parameters.
The tt background was modelled using P -B v2 interfaced to P 8 for parton showering, hadronisation, and the underlying event, and to E G v1.2.0 for heavy-flavour hadron decays. For this sample, the A14 tune was used [122]. Simulated Z + jets and tt background samples are normalised to the data-driven estimates described in Section 6. Generated events were processed through the ATLAS detector simulation [123] within the G 4 framework [124] and reconstructed in the same way as collision data. Additional pp interactions in the same and nearby bunch crossings are included in the simulation. The pile-up was modelled by overlaying the original hard-scattering event with simulated inelastic pp events generated with P 8.186 [82] using the NNPDF2.3LO set of PDFs [125] and the A3 tune [126].

Event reconstruction and selection
The details of the selection and reconstruction of Higgs boson candidate events are provided in Ref.
[32], while a brief description is provided here. Single-lepton, dilepton, and trilepton triggers are employed and ensure a signal selection efficiency above 98%. Data events are subjected to quality requirements and are required to have at least one vertex with two associated ID tracks with transverse momentum p T > 500 MeV. The primary interaction vertex is selected as the one with the largest p 2 T of all associated tracks.
The lepton identification requirements follow the inclusive event selection described in Ref. [32]. All muons are required to satisfy p T > 5 GeV and |η| < 2.7, except those that are reconstructed with ID tracks matched to energy deposits in the calorimeter (calorimeter-tagged), which must satisfy p T > 15 GeV and |η| < 0.1. No more than one calorimeter-tagged or stand-alone muon is allowed per event, where stand-alone muons have not been matched to an ID track. Electrons are required to satisfy E T > 7 GeV and |η| < 2.47. Jets are reconstructed using the anti-k t algorithm with a radius parameter R = 0.4 and applied to Particle Flow objects [127]. Jets are required to have p T > 30 GeV and |η| < 4.5. Jets within |η| < 2.5 are identified as containing a b-hadron using the MV2c10 b-tagging algorithm at the 70% efficiency working point [128,129]. If a jet overlaps geometrically with a reconstructed a muon (electron) within a cone of radial size ∆R = 0.1(0.2), the jet is removed.
Same-flavour opposite-charge (SFOC) lepton pairs are selected to form Higgs boson candidates. The SFOC lepton pair with mass m 12 closest to the Z boson mass is called the leading pair, while the other becomes the subleading pair, with mass m 34 . If multiple combinations of SFOC pairs exist, the Higgs boson candidate with m 12 closest to the Z boson mass is chosen. The three leading leptons of each Higgs boson candidate are required to satisfy p T > 20, 15, 10 GeV. Higgs boson candidate events are subjected to further selection requirements on the dilepton masses, lepton separation, J/ψ veto, impact parameter significance (d 0 /σ(d 0 )), and vertex quality, as outlined in Table 2. In addition, isolation requirements are imposed on the leptons to suppress the tt and Z + jets reducible backgrounds. If an extra prompt lepton with p T > 12 GeV passing all identification and isolation requirements detailed previously is present in the event, the final Higgs boson candidate is chosen using a method based on the matrix element (ME). The matrix element is calculated at LO using M G 5_ MC@NLO and the quadruplet with the highest ME value is chosen. This increases the probability of selecting the correct Higgs boson candidate in cases where the extra lepton comes from the decay of a vector boson or top quark in VH-leptonic or ttH/tH production. The four-lepton mass resolution is improved by accounting for reconstructed final-state radiation (FSR) photons in the Z boson decay. After selection criteria are applied, events are divided into bins for each variable of interest for the differential cross-section measurements. Finally, all measurements presented in this paper are performed within a four-lepton mass window of 105 < m 4 < 160 GeV. The signal selection efficiency is about 31%, 21%, 17%, and 16% for the 4µ, 2e2µ, 2µ2e, and 4e final states, respectively. Here, the first lepton pair refers to the lepton pair with an invariant mass closest to the Z boson mass.

Fiducial phase space and unfolded observables
The fiducial cross sections are defined using simulation at particle level and the selection requirements outlined in Table 3. In order to minimise model-dependent acceptance extrapolations, these are chosen to closely match the selection requirements of the detector-level analysis after the event reconstruction.
The fiducial selection is applied to final-state electrons and muons that do not originate from hadrons or τ-lepton decays, after 'dressing' them, i.e., the four-momenta of photons within a cone of size ∆R = 0.1 around the lepton are added to the lepton's four-momentum. The photons which originate from hadron decays are excluded. Particle-level jets are reconstructed from final-state neutral and charged particles using the anti-k t algorithm with radius parameter R = 0.4. Electrons, muons, neutrinos (if they are not from hadron decays) and photons from Higgs decays as well as those used to dress leptons are excluded from the jet clustering. A jet is labelled as a b-jet if there is a b-hadron with p T > 5 GeV within a cone of size ∆R = 0.3 around the jet axis. Jets are removed if they are within a cone of size ∆R = 0.1 around a selected lepton.
Quadruplet selection using the selected dressed leptons follows the same procedure as for reconstructed events. In the case of VH or ttH production, additional leptons not originating from a Higgs boson decay can induce a 'lepton mispairing' when assigning them to the leading and subleading Z bosons. To improve the lepton pairing efficiency, the matrix-element-based pairing method as described in Section 4 is employed. The variables used in the differential cross-section measurement are calculated using the dressed leptons of the quadruplets.
The acceptance of the fiducial selection, defined as the ratio of the number of events passing the particle-level selection to the number of events generated in a given bin or final state (with respect to the full phase space of H → Z Z * → 2 2 , where , = e or µ), is about 49% for each final state for a SM Higgs boson with m H = 125 GeV. The ratio of the number of events passing the selection after detector simulation and event reconstruction to those passing the particle-level selection is about 45%. Due to resolution effects, about 1.6% of the events which pass the detector-level selection fail the particle-level selection.
GeV for all SFOC lepton pairs Mass window 105 GeV< m 4 < 160 GeV If extra lepton with p T > 12 GeV Quadruplet with largest matrix element value Within the fiducial phase space defined above, differential cross sections are measured for variables which are sensitive to both the production and decay of the Higgs boson. For example, the transverse momentum distribution of the Higgs boson provides a test of perturbative QCD calculations, is sensitive to the structure of the Higgs boson interactions and is sensitive to charm and bottom Yukawa couplings. The rapidity of the Higgs boson is sensitive to the choice of parton distribution functions for the colliding protons, and is also influenced by QCD radiative corrections. The invariant masses of the leading and subleading lepton pair are sensitive to higher-order electroweak corrections to the Higgs boson decay, and are sensitive to BSM contributions. These two variables and the angular variables of the Higgs boson decay are also of interest due to their sensitivity to the spin and parity of the Higgs boson, as well as to same-flavour pair final-state interference and EW corrections. Variables related to jets probe QCD radiation effects and the Higgs boson production. The jet multiplicity is sensitive to different production mechanisms and provides sensitivity to the theoretical modelling of high-p T quark and gluon emission. The transverse momentum of the jets directly probes the quark and gluon radiation. The invariant mass of the two leading jets is also sensitive to the production mechanisms of the Higgs boson, while the signed angle in the transverse plane of the two leading jets is a test of the spin and parity of the Higgs boson. Jet-related variables, in particular double differential variables, also probe the effects of QCD resummation. Additional variables which combine the properties related to the kinematics of the Higgs boson and the jets are also considered. A summary of all the variables and their descriptions is given in Table 4. Table 4: Definitions of observables for which differential cross sections are measured. The angular variables are defined as in Ref. [130]. In addition to the single observables listed, the following double differential observables are built using variables defined below: m 12 vs. m 34  vs. |y lead. jet | (where |y lead. jet || is the rapidity of the leading jet). Jet-related variables are inclusive, while for the jet multiplicity the results are provided in both the inclusive and exclusive jet bins. ∆φ j j is defined as φ lead. jet − φ sublead. jet if η lead. jet > η sublead. jet or as φ sublead. jet − φ lead. jet if η sublead. jet > η lead. jet . If ∆φ j j < 0, 2π is added to the value.

Higgs boson kinematic-related variables
T , |y 4 | Transverse momentum and rapidity of the four-lepton system m 12 , m 34 Invariant mass of the leading and subleading lepton pair | cos θ * | Magnitude of the cosine of the decay angle of the leading lepton pair in the four-lepton rest frame relative to the beam axis cos θ 1 , cos θ 2 Production angles of the anti-leptons from the two Z bosons, where the angle is relative to the Z vector. φ, φ 1 Two azimuthal angles between the three planes constructed from the Z bosons and leptons in the Higgs boson rest frame.

Jet-related variables N jets , N b-jets
Jet and b-jet multiplicity p lead. jet T , p sublead. jet T Transverse momentum of the leading and subleading jet, for events with at least one and two jets, respectively. Here, the leading jet refers to the jet with the highest p T in the event, while subleading refers to the jet with the second-highest p T . m j j , |∆η j j |, ∆φ j j Invariant mass, difference in pseudorapidity, and signed difference in φ of the leading and subleading jets for events with at least two jets Higgs boson and jet-related variables p Transverse momentum and invariant mass of the four-lepton system and leading and subleading jets, for events with at least two jets

Background estimation
Non-resonant SM (Z ( * ) /γ * )(Z ( * ) /γ * ) production via qq annihilation and gluon-gluon fusion, referred to as Z Z * , can result in four prompt leptons in the final state and constitutes the largest background for this analysis. While for previous analyses [11,12] both the shape and the normalisation of this background were exclusively estimated with simulation, in this paper the normalisation is constrained with a data-driven technique. The systematic uncertainty is reduced because both the theoretical and luminosity uncertainties no longer contribute to the normalisation uncertainty. The normalisation of the non-resonant Z Z * component, which dominates outside the Higgs boson peak region, is obtained from data by extending the mass interval considered from 115-130 GeV to 105-160 GeV. The increased mass interval allows an estimation of this process with minimal impact on the expected sensitivity for the signal process. This contribution is determined as part of the 4 mass fit (discussed in Section 7) in the full four-lepton mass region, with the shape of the background taken from simulation.
The Z Z * normalisation is estimated separately in each bin of each differential observable, where a different Z Z * scaling factor is used for each observable bin. In phase-space regions where the Z Z * component in the m 4 sidebands is too low to provide a reliable estimate of its contribution, the estimate is evaluated simultaneously for several differential bins.2 Other background processes, such as Z + jets, tt, and W Z, contain at least one jet, photon or lepton from a hadron decay that is misidentified as a prompt lepton. These reducible backgrounds are significantly smaller than the non-resonant Z Z * background and are estimated using data where possible, following slightly different approaches for the + µµ and + ee final states [11,12,32].
In the + µµ final states, the normalisations for the Z + jets and tt backgrounds are determined by performing fits to the invariant mass of the leading lepton pair in dedicated independent control regions which target each background process for each bin of the differential observables. Depending on the background process being targeted, the control regions are formed by relaxing the χ 2 requirement on the four-lepton vertex fit, and by inverting or relaxing isolation and/or impact-parameter requirements on the subleading muon pair. Additional control regions (eµµµ and + µ ± µ ± ) are used to improve the background estimate by reducing the statistical uncertainty of the fitted normalisation. Transfer factors to extrapolate from the control regions to the signal region are obtained separately for tt and Z + jets using simulation. This method is performed in each differential bin. The m 4 shape for both processes in each bin is obtained from simulation.
The + ee control-region selection requires the electrons in the subleading lepton pair to have the same charge, and relaxes the identification, impact parameter and isolation requirements on the electron candidate with the lowest transverse energy. This electron candidate, denoted by X, can be a light-flavour jet, an electron from photon conversion or an electron from heavy-flavour hadron decay. The heavy-flavour background is completely determined from simulation, whereas the light-flavour and photon conversion background is obtained with the sPlot method [131]. This is based on a fit to the number of hits in the innermost ID layer in the data control region. Transfer factors to extrapolate from the + ee control region to the signal region for the light-flavour jets and converted photons, obtained from simulated samples, are corrected using a Z + X data control region. The corrected transfer factors are then used to extrapolate the extracted yields to the signal region. Both the extraction of the global yield in the control region and the extrapolation to the signal mass region are performed in bins of the transverse momentum of the electron candidate and the jet multiplicity. In order to extract the shape of the backgrounds from light-flavour jets and photon conversions for each observable, a similar method is used, except that the extraction and extrapolation is performed only as a function of the transverse momentum of the electron candidate, ignoring the binning in jet multiplicity.
Additional contributions from rare processes, such as t X X (tt Z, ttW, tW Z and other rare top-associated processes) and VVV are estimated from simulation.

Signal extraction and unfolding
To extract the number of signal events in each bin of a differential distribution (or for each decay final state for the inclusive fiducial cross section), invariant mass templates for the Higgs boson signal and the background processes are fitted to the m 4 distribution in data. Compared to the previous analysis [11], the non-resonant Z Z * background is fitted simultaneously with the signal and constrained by extending the m 4 fit range from 115-130 GeV to 105-160 GeV. For the total and fiducial cross sections in different final states, the same normalisation factor is used for the Z Z * contribution. For the differential cross-section measurements, multiple Z Z * normalisation factors are introduced in the model, as described in Section 6. The reducible background, composed of Z + jets, tt, and W Z processes, is estimated from dedicated control regions as described in Section 6 and its overall normalisation and shape can vary within the associated systematic uncertainties. Finally, for the differential distributions, no splitting into decay final states is performed, and the SM Z Z * → 4 decay fractions are assumed.
The number of expected events N i in each observable reconstruction bin i, expressed as a function of m 4 , is given by where A j is the acceptance in the fiducial phase space and σ j the total cross section in fiducial bin j, L is the integrated luminosity, B is the branching ratio and N bkg i (m 4 ) is the background contribution. The index j runs over all observable bins in the fiducial phase space. The term P i (m 4 ) is the m 4 signal shape containing the fraction of events as a function of m 4 expected in each reconstruction bin, taken from MC simulation. The term r i j represents the detector response matrix, created with simulated signal samples and averaged across the different production modes. These factors correspond to the probability that an event generated within the fiducial volume in the observable bin j is reconstructed in bin i.
The normalisation, f nonfid i , represents the fraction of events which are outside of the fiducial region but are reconstructed within the signal region. This ranges from 1.1% to 1.7% depending on the bin of the unfolded observable or final state.
The detector response matrix accounts for bin-to-bin migrations in the unfolding. It was chosen over the bin-by-bin correction factor technique used in the previous analyses [11,12] due to its lower model dependence. Biases introduced via the unfolding method are minimised when using the response matrix; however, matrix unfolding can amplify small fluctuations in data when the response matrix is characterised by a large condition number.3 The binning choice made for all observables ensures a statistical significance of more than 2σ for the signal process. The binning is also chosen to minimise migrations between bins. In general, the bin width is more than twice the experimental resolution. As a result, the response matrices for all the variables considered are well-conditioned, with a condition number less than 2.5. The fluctuations of the unfolded distribution can be further reduced using regularisation techniques. Unfolding tests done with toy data sets indicate that while regularisation provides a modest reduction of the statistical uncertainty, this reduction is counterbalanced by the bias introduced by this technique. Therefore, no regularisation of the unfolding was applied.   particle-level events in that bin. For the N jets observable, the migrations are more relevant due to the relatively worse jet energy resolution and the presence of pile-up jets in the reconstructed events. This brings the purity for the for N jets ≥ 3 bin down to 68%. The p lead. jet T migrations are similarly larger, with the lowest purity value of 67% occurring in the lowest p lead. jet T bin. The m 12 vs. m 34 observable, like p 4 T , has a higher purity. All bins have a purity of around 90% except the first bin, which has a purity of 78%.

Systematic uncertainties
The systematic uncertainties include experimental uncertainties, such as those in object reconstruction, identification, isolation, resolution, and trigger efficiencies, as well as theoretical uncertainties related to the modelling of the signal and background processes. More detail is provided in Ref.
[32], while a brief overview of the dominant sources of uncertainty is provided here. The impacts of the experimental and theoretical uncertainties on the measurements are summarised in Table 5.

Experimental uncertainties
The uncertainty in the predicted yields due to pile-up modelling ranges between 1% and 2%. The uncertainty in the integrated luminosity is 1.7% and affects the signal yields and simulated background estimates when not constrained by the sidebands.
The electron (muon) reconstruction and identification efficiency uncertainties are approximately 1.0-2.0% (< 1.0%). The uncertainty in the expected yields due to the muon and electron isolation efficiencies is also considered, and is approximately 1%. Lepton energy momentum scale and resolution uncertainties have negligible impacts on the presented results.
The impact of uncertainties in the jet energy scale and resolution (of between 1% and 3%) is only relevant for the jet-related differential cross-section measurements, where their impact is typically between 3% and 5%, and is negligible in the other measurements. The uncertainty in the performance of the b-tagging algorithm is at the level of a few percent over most of the jet p T range [129].
The impact of the precision of the Higgs boson mass measurement, m H = 125.09 ± 0.24 GeV [15], on the signal acceptance due to the signal region mass-window requirement is negligible.
For the data-driven measurement of the reducible background, three sources of uncertainty are considered: statistical uncertainty, overall systematic uncertainty for each of + µµ and + ee, and a shape systematic uncertainty which varies with the differential variable. Impacts from these sources of uncertainty range from less than 1% to a maximum of around 3%. The inclusive reducible background estimate has a relatively small (3%) statistical uncertainty, which has minimal impact on the cross section.

Theoretical uncertainties
Sources of theoretical uncertainty include missing higher-order corrections, parton shower and underlying event modelling, and PDF+α s uncertainties, and these all affect modelling of the signal and background processes. For measurements of the cross section, the impact of these theory systematic uncertainties on the signal comes from their effects on the response matrix. Table 5: Fractional uncertainties for the inclusive fiducial and total cross sections, and ranges of systematic uncertainties for the differential measurements. The columns 'e/µ' and 'Jets' represent the experimental uncertainties in lepton and jet reconstruction and identification, respectively. The Z + jets, tt, t X X (Other Bkg.) column includes uncertainties related to the estimation of these background sources. The Z Z * theory (Z Z * th.) uncertainties include the PDF and scale variations. Signal theory (Sig th.) uncertainties include PDF choice, QCD scale, and shower modelling of the signal. Finally, the column labelled 'Comp.' contains uncertainties related to production mode composition and unfolding bias which affect the response matrices. The uncertainties have been rounded to the nearest 0.5%, except for the luminosity uncertainty, which has been measured to be 1.7%.  The prediction of the ggF process in different N jets categories and migration effects on the N jets ggF cross sections are large sources of theoretical uncertainty, which are accounted for using the approach detailed in Ref. [107]. The QCD scale uncertainty from the factorisation and renormalisation scales, resummation scales, and migrations between N-jet phase-space bins are considered [53,[132][133][134][135]. The impact of QCD scale variations on the Higgs boson p T distribution as well as the uncertainty of the p T distribution in the 0-jet bins are also taken into account. Higher-order impacts on the p T distribution predictions due to treating the top quark mass as infinite in the heavy-quark loop are accounted for by comparing these predictions with finite-mass calculations. For the VBF production mode, the uncertainty due to missing higher orders in QCD are considered, including migration effects in number of jets, transverse momentum of the Higgs boson, transverse momentum of the Higgs boson and leading dijet system, and the invariant mass of the two leading jets as outlined in the scheme presented in Ref. [136].
For production modes other than ggF and VBF, the effects of QCD scale uncertainties are estimated by considering all configurations of renormalisation and factorisation scales varied by a factor of two. In each experimental bin, the largest difference between all the variations and the nominal configuration is assigned as uncertainty.
The effects of parton shower and multiple-parton interaction modelling uncertainties on the acceptance are estimated using tune eigenvector variations as well as comparisons between acceptances calculated with P 8 and H 7 parton showering algorithms.
PDF uncertainty impacts are estimated using the eigenvector variations of the PDF4LHC_NLO_30 Hessian PDF set, following the PDF4LHC recommendations [72].
For the cross sections extrapolated to the full phase space, an additional uncertainty (2.2%) related to the H → Z Z * branching ratio [69,70] is included in the measurement.
Since the Z Z * process normalisation is constrained by performing a simultaneous fit of sideband regions enriched in this contribution together with the signal region, most of the theoretical uncertainty in the normalisation for this background vanishes.4 The uncertainties due to missing higher-order effects in QCD are estimated by varying the factorisation and renormalisation QCD scales by a factor of two; the impact of the PDF uncertainty is estimated using the MC replicas of the NNPDF 3.0 PDF set. Uncertainties due to the parton shower modelling for the Z Z * process are considered as well. The impact of these uncertainties is below 2% for all the fiducial differential cross sections. In addition, the m 4 shape obtained from S is compared with that obtained from P and M G 5_ MC@NLO and the difference is taken as an additional source of systematic uncertainty. In each m 4 bin, the largest difference between S and P or M G 5_ MC@NLO is used, and the systematic uncertainty is determined by interpolating between these shapes.
The uncertainty in the gluon-induced Z Z * process is taken into account as well by changing the relative composition between the quark-initiated and gluon-initiated Z Z * components according to the theoretical uncertainty in the predicted cross sections.
Finally, unfolding-related uncertainties arise from uncertainties in the production mode composition that affect the response matrices, as well as from uncertainties in the bias introduced by the unfolding method. For the former, an uncertainty is assessed by varying the production cross sections within their measured uncertainties taken from Ref. [12], and has an impact of less than 1%. In the latter case, the uncertainty in the bias is obtained independently per bin by comparing the unfolded cross section from simulation with that expected when varying the underlying true cross sections of the simulated data sample within the expected statistical error. The impact of this uncertainty is typically negligible in distributions such as p 4 T , where the response matrix is largely diagonal, but can be of the order of 10% in distributions with larger bin migrations, such as N jets .

Results
Results are presented for the full set of inclusive and differential variables outlined in Section 5. Section 9.1 presents the data yields from the full Run 2 data set. Section 9.2 provides details of the statistical procedure used for the extraction of the measurements. Cross-section results, and comparisons with SM predictions, are provided in Sections 9.3 and 9.4.

Measured data yields
The observed number of events in each of the four decay final states, and the expected signal and background yields before fitting to data (pre-fit), are presented in Table 6. These events have passed the event selection and fall in a narrow window around the Higgs boson mass peak (115< m 4 <130 GeV).       In distribution (c), the first bin contains events with zero jets, while in distribution (d), the first bin contains events with fewer than two jets. In both (c) and (d), all bins except the first are divided by the bin width. The uncertainty in the prediction is shown by the hatched band, which includes the theoretical uncertainties of the SM cross section for the signal and the Z Z * background.

Statistical analysis
The inclusive fiducial and differential cross sections are measured using a binned profile-likelihood-ratio fit [137], taking into account all bins of a given distribution. The likelihood function includes the shape and normalisation uncertainties of the signal and background predictions as nuisance parameters, as outlined in Section 8. The cross sections are extracted by minimising two times the negative logarithm of the profile likelihood ratio, −2 ln Λ. In the asymptotic approximation, i.e. the large sample limit, −2 ln Λ behaves as a χ 2 distribution with one degree of freedom. The compatibility of a measured cross section and its theoretical prediction is tested by computing a p-value based on the difference between the value of −2 ln Λ at the best-fit value and the value obtained by fixing the cross section in each bin to that predicted by theory. These p-values do not include the uncertainties in the theoretical predictions. For all measured observables the asymptotic approximation is validated with pseudo-experiments, and where the number of observed events is less than three, the uncertainties are corrected to the values obtained with the pseudo-experiments.
For the fiducial and differential cross-section measurements, the fitted m 4 distribution in each final state or differential bin is used to extract the measured cross section following Eq. (1). The fiducial cross sections of the four final states can either be summed to obtain an inclusive fiducial cross section, or they can be combined assuming the SM Z Z * → 4 relative branching ratios. The latter combination is more model dependent, but benefits from a smaller statistical uncertainty.

Inclusive fiducial cross-section measurements
The fiducial production cross sections of the H → Z Z * → 4 process are presented in Table 7 and Figure 18. The left panel in Figure 18(a) shows the fiducial cross sections for the four individual decay final states: 4µ, 4e decays (hereafter referred to as same flavour), and 2µ2e, 2e2µ decays (hereafter referred to as different flavour). The middle panel shows the cross sections for same-and different-flavour decays, which can provide a probe of same-flavour interference effects, as well as the inclusive fiducial cross sections obtained by either summing all 4 decay final states or combining them assuming relative SM branching ratios.
The data are compared with the SM prediction after accounting for the fiducial acceptance as determined from the SM Higgs boson simulated samples (see Section 3).
The combined inclusive fiducial cross section is extrapolated to the full phase space, as shown in the right panel of Figure 18, using the fiducial acceptance as well as the branching ratios, with the uncertainties described in Section 8. The total cross section is also compared with the cross sections predicted by NNLOPS, M G 5_ MC@NLO-FxFx (MG5-FxFx) and H 2.3 [52,138] for ggF, while for all other production modes the predictions described in Section 3 are used. For ggF, all generators predict cross sections that are lower than the N 3 LO calculation. The p-values, calculated as described in Section 9.2, are shown in Table 7. The probability of compatibility of the measured fiducial cross section (σ comb ) and the Standard Model expectation is at the level of 67%.  Figure 18: (a) The fiducial cross sections (left two panels) and total cross section (right panel) of Higgs boson production measured in the 4 final state. The fiducial cross sections are shown separately for each decay final state, and for same-and different-flavour decays. The inclusive fiducial cross section is measured as the sum of all final states, as well as by combining the per-final-state measurements assuming SM Z Z * → 4 relative branching ratios. The total SM prediction is accurate to N 3 LO in QCD for the ggF process. The cross sections for all other Higgs boson production modes X H are added. For the fiducial cross-section predictions, the SM cross sections are multiplied by the acceptances determined using the NNLOPS sample for ggF and the samples discussed in Section 3 for the other production modes. For the total cross section, the predictions by the generators NNLOPS , H , and M G 5_ MC@NLO-FxFx are also shown. The error bars on the data points show the total uncertainties, while the systematic uncertainties are indicated by the boxes. The shaded bands around the theoretical predictions indicate the PDF and scale uncertainties, calculated as described in Section 8.2. (b) The correlation between the fiducial cross sections for the four individual decay final states and the Z Z * normalisation factor. Table 7: The fiducial and total cross sections of Higgs boson production measured in the 4 final state. The fiducial cross sections are given separately for each decay final state, and for same-and different-flavour decays. The inclusive fiducial cross section is measured as the sum of all final states (σ sum ), as well as by combining the per-final-state measurements assuming SM Z Z * → 4 relative branching ratios (σ comb ). For the total cross section (σ tot ), the Higgs boson branching ratio at m H 125 GeV is assumed. The total SM prediction is accurate to N 3 LO in QCD and NLO EW for the ggF process. For the fiducial cross-section predictions, the SM cross sections are multiplied by the acceptances determined using the NNLOPS sample for ggF. For all the other production modes, the cross sections from the samples discussed in Section 3 are added. The p-values indicating the probability of compatibility of the measurement and the SM prediction are shown as well. They do not include the systematic uncertainty in the theoretical predictions.

Differential cross-section measurements
The measured differential production cross sections for the transverse momentum p 4 T of the Higgs boson are shown in Figure 19, while the measured differential cross sections with respect to the masses of the leading and subleading Z bosons resulting from the Higgs boson decay, m 12 and m 34 , are provided in Figure 20. Figures 21, 22, and 23 show the measured differential production cross sections with respect to angular variables, |y 4 |, |cos θ * |, cos θ 1 , cos θ 2 , φ, and φ 1 , that probe the kinematics of the Higgs boson decay products.
Differential production cross-section measurements with respect to variables that probe the jet activity in reconstructed Higgs boson events follow in Figures 24-28. These include the exclusive and inclusive jet multiplicities, N jets , the b-jet multiplicity, N b−jets , variables measuring the transverse momentum of the jets, p lead. jet T and p sublead. jet T , as well as variables that probe the kinematics of pairs of jets in events with at least two jets, m jj , ∆η jj , and ∆φ jj .
In addition, differential cross-section measurements are provided for observables aimed at studying the relationship between the reconstructed Higgs boson and accompanying jets. These are presented in Figures 29-30.    The central panel of (a) shows the ratio of different predictions to the data, and the grey area represents the total uncertainty of the measurement. The bottom panel of (a) shows the ratios of the fitted values of the Z Z * normalisation factors to the predictions from MC simulation discussed in Section 3. As indicated by the horizontal error bars, the Z Z * normalisation is estimated in each of the first three p 4 T bins separately, while the next two bins share a common estimation factor, as do the last five bins.                 Figure 24: Differential fiducial cross sections for (a) the jet multiplicity, N jets , in the selected events, and (c), the inclusive jet multiplicity. In the N jets distribution in (a), the first three bins are exclusive in number of jets, while the fourth is inclusive. The corresponding correlation matrix between the measured cross sections and the Z Z * background normalisation factors is also shown in (b). In the N jets distribution in (c), all bins are inclusive, with the first bin including all events, the second including all events with at least one jet, and so on. The first bin is filled with events which do not have any jets, the second is filled with events with at least one jet but no b-tagged jets, while the third includes all events with at least one b-tagged jet. The corresponding correlation matrices between the measured cross sections and the Z Z * background normalisation factors are also shown in (b)).   , in events with at least two jets. Leading and subleading jets refer to the jets with the highest and second-highest transverse momenta. The first bin contains events which do not pass the jet requirements. The corresponding correlation matrices between the measured cross sections and the Z Z * background normalisation factors are also shown ((b) and (d)).    Figure 28: Differential fiducial cross sections for (a) the distance between these two jets in pseudorapidity, ∆η jj , and (c) the distance between the two jets in φ, ∆φ jj . The first bin contains events with fewer than two jets that pass the jet selection requirements. Finally, the corresponding correlation matrices between the measured cross sections and the Z Z * background normalisation factors are provided ((b) and (d)).                      Figure 39: (a) Double differential fiducial cross section for the transverse momentum of leading vs. subleading jet, p lead. jet T vs. p sublead. jet T , and (b) the corresponding correlation matrix between the measured cross sections and the Z Z * background normalisation factor. The bin boundaries are defined in Figure 16.
Overall, there is good agreement between measured cross sections and predictions. Small differences between measurement and prediction occur in several of the angular observables, as well as in bins of m 4 jj , and several of the double differential measurements. For example, the p-value for the double differential distribution p lead. jet T vs. |y lead. jet | in Figure 38 is particularly low due to the downward fluctuation in bin 2. However, when considering the size of the uncertainties these differences are not significant. Since no events are observed in the highest bin for p 4 T in Figure 19, an upper limit of 27 ab at 95% confidence level (CL) is set on the cross section using CLs [139]. Similarly, a limit of σ < 38 ab at 95% CL is also set in the last bin of p 10 Interpretation of differential distributions The measured differential fiducial cross sections can be used to probe possible effects of physics beyond the SM. Two possible interpretations of the results are presented. In Section 10.1, the m 12 vs. m 34 double differential cross section is used to probe several BSM scenarios within the framework of pseudoobservables [140], while in Section 10.2, the p 4 T differential cross section is used to constrain the Yukawa couplings of the Higgs boson with the band c-quarks [141].

Constraints on BSM effects within the pseudo-observables framework
In this interpretation, the couplings related to the BSM contact interactions of the Higgs boson decaying into four leptons are considered. As defined in Ref. [142], the pseudo-observables framework introduces modified contact terms between the Higgs boson, the Z boson, and the left-or right-handed leptons Z, L and Z, R . In order to reduce the number of independent parameters considered in the pseudo-observables framework for the H → 4 decay amplitudes, specific symmetries are imposed [142]. In all the scenarios considered, the parameters associated with other pseudo-observables affecting the angular distributions, such as (C P) Z Z , (C P) Zγ and (C P) γγ , are set to zero. Thus, the contact terms considered have the same Lorentz structure as the SM term and only affect the dilepton invariant mass distributions.
Four scenarios are investigated [142]. In the first scenario, referred to as the flavour-universal contact terms, the parameters of interest are the Z, L and the Z, R couplings, where the interactions described by these contact terms have the same strength for electrons and muons. The second scenario considered is linear EFT-inspired, where lepton-flavour universality is again imposed and the Higgs boson is assumed to be part of a SU(2) L doublet. This is reflected in the condition R = 0.48 L [142]. The parameters of interest are L and the coupling strength of the Higgs boson to the Z boson, κ Z Z . In the following two scenarios, lepton-flavour universality can be violated. For the third scenario, referred to as flavour non-universal vector contact terms, the helicity structure of the couplings is fixed to be vector ( Z,e L = Z,e R , Z,µ L = Z,µ R ) and the independent parameters are the couplings to electrons Z,e R and muons Z,µ R . Finally, a fourth scenario with flavour non-universal axial-vector contact terms is considered. In this case the helicity structure of the couplings is fixed to be axial-vector, with the parameters of interest being the couplings to electrons Z,e R and muons Z,µ R and the condition Z, L = − Z, R is imposed. Using the m 12 vs. m 34 double differential cross sections for these interpretations provides sensitivity to distinguish between potential contributions from the contact terms and those from changes to the coupling strength of the Higgs boson to the Z boson.
The variation of the fiducial cross section as a function of the BSM couplings is computed relative to the SM by M G 5_ MC@NLO in each of the bins of the measured m 12 vs. m 34 differential cross section. This is done for a grid of points in the BSM parameter space in each scenario. These relative variations are then fit to a two-dimensional quadratic function. The parameterisation, which also includes any changes in the acceptance, is then encoded into the likelihood and corresponding limits are set for each scenario. Figure 40 shows the limits on BSM interactions of the Higgs boson for the four considered cases. The corresponding 95% confidence intervals for each of the parameters are listed in Table 8.  Table 8: Confidence intervals for the scenarios considered in the pseudo-observables framework. Based on the observed 2D exclusion contours, 1D exclusion intervals are provided for the EFT-inspired, flavour non-universal vector, and flavour non-universal axial-vector scenarios. The observed limits are calculated while profiling the other parameters of interest. For the EFT-inspired interpretation, the limits are derived assuming κ Z Z ≥ 0. This constraint has no impact on the limit as the analysis is not sensitive to the sign of this parameter.

Interpretation
Parameter best-fit value 95% confidence interval

Constraints on Yukawa couplings
Although the couplings of the Higgs boson to the top and bottom quarks have been established recently, obtaining evidence for the coupling of the Higgs boson to the charm quark is more challenging. Direct methods are limited either by low branching fraction (H → J/ψγ → µ + µ − γ) or by large backgrounds (H → cc). Nevertheless, it has been shown recently that it is possible to indirectly constrain the Yukawa coupling to quarks by analysing the p H T spectrum [19,141]. In particular, the effects of BSM contributions to the coupling modifiers for the Higgs boson to charm quarks, κ c , and for the Higgs boson to bottom quarks, κ b , are investigated.
The fiducial cross section is parameterised as a function of the κ c and κ b values in each measured bin of p 4 T . Both the gluon-initiated and quark-initiated components of the prediction show a larger variation, different in size and shape, of the cross section especially at p 4 T < 10 GeV. The theoretical uncertainties of these predictions are calculated separately for the gluon-initiated and quark-initiated components by varying the normalisation and factorisation scales by factors of two. The configuration with largest uncertainty across all the p 4 T bins across κ c ∈ [−10, 10] and κ b ∈ [−2, 2] ranges is used to define the systematic uncertainty for the predictions. The impact of this uncertainty is about 20% on the expected limits.
Three different scenarios are considered, with an increasing level of model dependency. In the first case, the modified fiducial cross sections in each bin due to value of the band c-quark Yukawa couplings are fit to the data together with a global normalisation factor. The corresponding observed limits on κ c and κ b are shown in Figure 41(a). The sensitivity in this case comes mainly from the modification of the shape induced by κ c and κ b , while possible overall normalisation effects are factorised out. In a second scenario, no additional normalisation factor is introduced in the likelihood and the obtained limits for the Yukawa couplings are shown in Figure 41(b). Finally, in a third scenario, a modification to the total width, and correspondingly to the branching ratio as function of the modified Yukawa couplings, is also encoded in the likelihood and the corresponding limits are shown in Figure 41(c). The 95% confidence intervals for the first and second scenarios are also listed in Table 9. These are comparable to results from direct searches in V H, H → cc [143,144]. Constraining κ b to the results from Ref. [145] leads to a less than 5% improvement in the observed limits for κ c for the scenarios considered.  Figure 41: Observed limits at 95% CL on Yukawa couplings κ c and κ b for the three scenarios considered: (a) only the p 4 T shape is used to constrain κ c and κ b ; (b) the predicted p 4 T differential cross section is used; (c) both the prediction of the p 4 T differential cross section and the modification to the branching ratio due to the κ c and κ b values are used. The p-values shown represent the probability of compatibility between the data and the p 4 T prediction corresponding to the best-fit values of κ c and κ b . The SM predictions ( * ) and the observed best-fit values (+) are indicated on the plots.

Summary
Fiducial inclusive and differential cross-section measurements of the Higgs boson in the H → Z Z * → 4 decay channel are presented. They are based on 139 fb −1 of √ s = 13 TeV proton-proton collisions recorded by the ATLAS detector at the LHC in 2015−2018. The inclusive fiducial cross section in the H → Z Z * → 4 decay channel is measured to be σ fid = 3.28 ± 0.30 (stat.) ± 0.11 (syst.) fb, in agreement with the Standard Model prediction σ fid,SM = 3.41 ± 0.18 fb. The measurement is about 40% more precise than the previous ATLAS result. The inclusive fiducial cross section is also extrapolated to the full phase space. Differential cross sections defined in a fiducial region close to the reconstructed event selection are measured for several variables sensitive to the Higgs boson production and decay such as the transverse momentum of the Higgs boson, the number of jets produced in association with the Higgs boson, the leading and subleading invariant masses of the lepton pairs. The measured cross sections are compared with different Standard Model predictions and in general good agreement is found. The results are also used to set new and more stringent constraints on BSM scenarios where contact term interactions in the H → 4 amplitudes are introduced. In addition, the p 4 T spectrum is used to constrain the band c-quark Yukawa couplings of the Higgs boson. In the scenario with minimal assumptions, values of κ c outside the range κ c ∈ [−12, +11] are excluded at 95% CL.  Figure 42: (a), (c) Differential fiducial cross sections as a function of the jet multiplicity, N jets , and leading jet p T in events with at least one jet, using a regularised matrix unfolding with the τ parameters set to 0.6 and 0.7 respectively. The corresponding correlation matrix between the measured cross sections and the Z Z * background normalisation factors are also shown in (b) and (d). Figure 43 presents results for the invariant mass of the leading lepton pair in same-flavour and differentflavour final states.  [7] CMS Collaboration, Constraints on the spin-parity and anomalous HVV couplings of the Higgs boson    [32] ATLAS Collaboration, Higgs boson production cross section measurements and their EFT interpretation in the 4 decay channel at √ s =13 TeV with the ATLAS detector, CERN-EP-2020-034, to appear soon.