Measurement of the total and differential Higgs boson production cross-sections at $$ \sqrt{s} $$ = 13 TeV with the ATLAS detector by combining the H → ZZ* → 4ℓ and H → γγ decay channels

<jats:title>A<jats:sc>bstract</jats:sc>                      </jats:title><jats:p>The total and differential Higgs boson production cross-sections are measured through a combined statistical analysis of the <jats:italic>H</jats:italic> → <jats:italic>ZZ</jats:italic><jats:sup>*</jats:sup> → 4<jats:italic>ℓ</jats:italic> and <jats:italic>H</jats:italic> → <jats:italic>γγ</jats:italic> decay channels. The results are based on a dataset of 139 fb<jats:sup><jats:italic>−</jats:italic>1</jats:sup> of proton–proton collisions at a centre-of-mass energy of 13 TeV, recorded by the ATLAS detector at the Large Hadron Collider. The measured total Higgs boson production cross-section is <jats:inline-formula><jats:alternatives><jats:tex-math>$$ {55.5}_{-3.8}^{+4.0} $$</jats:tex-math><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML">                   <mml:msubsup>                     <mml:mn>55.5</mml:mn>                     <mml:mrow>                       <mml:mo>−</mml:mo>                       <mml:mn>3.8</mml:mn>                     </mml:mrow>                     <mml:mrow>                       <mml:mo>+</mml:mo>                       <mml:mn>4.0</mml:mn>                     </mml:mrow>                   </mml:msubsup>                 </mml:math></jats:alternatives></jats:inline-formula> pb, consistent with the Standard Model prediction of 55<jats:italic>.</jats:italic>6 ± 2<jats:italic>.</jats:italic>5 pb. All results from the two decay channels are compatible with each other, and their combination agrees with the Standard Model predictions. A combined statistical interpretation of the measured fiducial cross-sections as a function of the Higgs boson transverse momentum is performed in order to probe the Yukawa couplings to the bottom and charm quarks. A similar interpretation is performed by including also the constraints from the measurements of Higgs boson production in association with a <jats:italic>W</jats:italic> or <jats:italic>Z</jats:italic> boson in the <jats:italic>H →</jats:italic><jats:inline-formula><jats:alternatives><jats:tex-math>$$ b\overline{b} $$</jats:tex-math><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML">                   <mml:mi>b</mml:mi>                   <mml:mover>                     <mml:mi>b</mml:mi>                     <mml:mo>¯</mml:mo>                   </mml:mover>                 </mml:math></jats:alternatives></jats:inline-formula> and <jats:inline-formula><jats:alternatives><jats:tex-math>$$ c\overline{c} $$</jats:tex-math><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML">                   <mml:mi>c</mml:mi>                   <mml:mover>                     <mml:mi>c</mml:mi>                     <mml:mo>¯</mml:mo>                   </mml:mover>                 </mml:math></jats:alternatives></jats:inline-formula> decay channels.</jats:p>

section_content: Following the discovery of a Higgs boson (H) with a mass around 125 GeV ten years ago [1, 2], by the ATLAS and CMS collaborations at the Large Hadron Collider (LHC) at CERN, an intense programme to measure the properties of this particle and compare them with those of the Higgs boson predicted by the Standard Model (SM) of particle physics [3, 4] has been carried out. 

In particular, total and differential fiducial Higgs boson production cross-sections have been measured, probing the kinematic features of the Higgs boson and of the particles produced in association with it. Both the ATLAS and CMS collaborations have measured total and differential fiducial Higgs boson production cross-sections at a proton-proton (pp) centre-of-mass energy √ s = 13 TeV in the H → ZZ * → 4 (where = e, µ) [5, 6], H → γγ [7, 8], H → W W * → eνµν [9], and H → τ τ [10] decay channels. The collaborations have also performed combinations of some of the most sensitive results [11, 12]. The measurements are performed in fiducial phase spaces that closely match the selection requirements of the detector-level analysis after the event reconstruction. This approach significantly reduces the model dependence that would otherwise be introduced by relying on the acceptance factors predicted by the model under consideration to extrapolate the measured signal yields to the full phase space.
section_name: Introduction

section_content: The most recent measurements of these cross-sections published by the ATLAS collaboration, exploiting 139 fb -1 [13, 14] of 13 TeV proton-proton collisions produced during the whole second data-taking phase of the LHC (Run 2, 2015-2018) and recorded by the ATLAS detector [15], have been performed using the H → ZZ * → 4 [5] and H → γγ [7] final states. 

The results of these two publications are combined in this article. The measurements are extrapolated to the full phase space and the measured cross-sections are compared with SM predictions. Additional systematic uncertainties introduced by the extrapolation to the full phase space are counterbalanced by a significant reduction of the statistical uncertainty of the measurement, which is the main limitation to the precision of the measurements in the individual decay channels. 

The measurements include the total production cross-section and one and twodimensional differential production cross-sections as a function of the Higgs boson transverse momentum 1 p H T, sensitive to perturbative QCD calculations, and of the Higgs boson rapidity |y H |, sensitive to the parton distribution functions (PDF). Furthermore, differential cross-sections for jet multiplicity N jets and the transverse momentum of the highest-p T jet p lead. jet T are also measured. Both N jets and p lead. jet T observables probe the theoretical modelling of high-p T QCD radiation in Higgs boson production. These distributions are also sensitive to the different Higgs boson production processes. The measurements provide a stringent test of the SM predictions and any deviations from these predictions can indicate the presence of physics beyond the SM (BSM). 

This article also presents a combined statistical interpretation, in terms of the b-and c-quark Yukawa coupling strengths to the Higgs boson, of the fiducial differential crosssections measured as a function of p H T in the two decay channels. Another interpretation, also including the constraints on the b and c Yukawa coupling strengths obtained from the measurements of Higgs boson production in association with a W or Z boson, with the Higgs boson decaying to b-or c-quark pairs [16, 17], is presented. 

The results presented in this article update and supersede those of a previous publication [11] based on the same final states and a partial Run 2 dataset corresponding to an integrated luminosity of 36. 1 fb -1. With respect to the previous publication, both measurements included in this article use an improved jet reconstruction [18] and an improved unfolding procedure that is based on a detector response matrix included in the likelihood fit. Full descriptions of the measurements and the respective improvements in the H → ZZ * → 4 and H → γγ decay channels used in this article are given in refs. [5, 7]. In both decay channels, the cross-sections in the full phase space are obtained from these unfolded yields by taking into account the luminosity, detector effects, acceptance factors, and branching fractions. The SM values of the Higgs boson branching fractions are as-
section_name: JHEP05(2023)028

section_content: sumed, and the acceptance factors are based on SM predictions. The value of the Higgs boson mass is assumed to be 125. 09 GeV [19]. 

The paper is organised as follows. Section 2 describes the simulated Higgs boson event samples and inclusive theory cross-section calculations used to obtain the total and fiducial cross-section predictions. The signal acceptance factors for extrapolating the results to the full phase space are detailed in section 3. The statistical procedure for the combination of the two channels is illustrated in section 4, yielding the results summarised in section 5. The differential cross-sections measured as a function of p H T are then used to constrain the Yukawa couplings of the Higgs boson to the bottom and charm quarks in section 6.
section_name: JHEP05(2023)028

section_content: The Monte Carlo (MC) event generators used for the calculation of the acceptance factors and detector effects, and for the SM predictions, are described in detail in refs. [5, 7]. Their main features are summarised in this section. 

Gluon-gluon fusion (ggF) events are simulated using Powheg NNLOPS [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] with the PDF4LHC15 next-to-next-to-leading order (NNLO) set of parton distribution functions [31], while other production modes are simulated with Powheg [20] [21] [22] with the PDF4LHC15 next-to-leading order (NLO) set except for b bH and tH, which are simulated using MadGraph5_aMC@NLO [32, 33] with the NNPDF3. 0 NLO PDF set [34]. These samples are generated assuming a Higgs boson with a mass m H = 125 GeV and are normalised to cross-sections obtained from the best available predictions as provided by the LHC Higgs Working Group [35] for m H = 125. 09 GeV, which are 48. 5 ± 2. 4 pb, 3. 78 ± 0. 08 pb, 2. 25 ± 0. 06 pb, 0. 49 ± 0. 11 pb and 0. 59 ± 0. 05 pb for the ggF, VBF, V H, b bH and t tH + tH processes respectively. In the case of the ggF NNLOPS prediction, this corresponds to a rescaling to the fixed order N 3 LO cross-section by a global K-factor of 1. 1. The impact of the 90 MeV difference between the values of m H used in the simulation and in the analysis is negligible, as discussed in section 3. 

For all production mechanisms the Pythia 8. 2 generator [36] is used to model the H → ZZ * → 4 and H → γγ decays, as well as for the parton shower and the underlying event. The AZNLO set of tuned parameters [37] is used for ggF, VBF and V H production, while the A14 tune [38] is used for the other production modes. Alternative ggF, VBF, V H, tH (t tH) samples are produced by interfacing the nominal matrix element generator with Herwig 7. 1. 3 (Herwig 7. 0. 4) [39, 40], using the H7UE set of tuned parameters [40], in order to estimate uncertainties in the signal acceptance factors related to the modelling of the parton shower. 

The measurements are also compared with an alternative prediction obtained by summing the expected cross-sections of non-ggF Higgs boson production processes described previously and an alternative SM ggF prediction obtained using Mad-Graph5_aMC@NLO (MG5 FxFx). This matrix-element generator provides NLO accuracy in QCD for zero, one, and two additional jets, using the FxFx merging scheme [32, 41], and includes the top and bottom quark mass effects [42] [43] [44]. The events are generated using the NNPDF30 NLO PDF set. The generator is interfaced to Pythia 8 for the
section_name: Higgs boson simulation samples and theoretical predictions

section_content: The acceptance factors that extrapolate at particle-level from the respective H → ZZ * → 4 and H → γγ fiducial phase spaces to the full phase space are estimated using the simulated event samples and cross-sections described in section 2. The definitions of the fiducial phase spaces are summarised in table 1 and table 2, with more details provided in refs. [5, 7] respectively. The evaluation of the acceptance factors assumes SM Higgs boson production fractions and a Higgs boson mass of 125 GeV: the 90 MeV difference from the measured mass value of 125. 09 GeV has a negligible impact on the Higgs boson kinematics. 

In the full phase space, the quantities p H T and |y H | are computed directly from the simulated Higgs boson momentum instead of its decay products, as in the fiducial analyses. The acceptance factors implicitly include the correction for this difference. Simulated particle-level jets are built from all stable particles with cτ > 10 mm, including neutrinos, photons, and leptons from hadron decays or produced in the shower. All decay products from the Higgs boson decay and the leptonic decays of associated vector bosons are removed from the inputs to the jet algorithm. Jets are reconstructed using the anti-k t algorithm [45] with a radius parameter R = 0. 4, and are required to have p T > 30 GeV. 

Theory uncertainties related to the PDF, higher-order corrections, and the parton shower model are taken into account when evaluating acceptance factors. For each channel, the uncertainties in the acceptance factors are correlated with the impact of these theoretical sources on the detector response matrix used in the unfolding. Due to this procedure, compared with the results in ref. [7], the H → γγ results presented in this article have these additional theoretical uncertainties in the detector response matrix. Uncertainties due to the PDF and missing higher-order corrections are estimated as described in -4 -
section_name: Acceptance factors

section_name: JHEP05(2023)028

section_content: Photon kinematics p T threshold for two leading photons:

Table 2. Summary of the particle-level fiducial definitions in the H → γγ analysis [7]. E iso T (∆R, p T, charged) is the scalar sum of the transverse momenta of charged stable particles with a transverse momentum above the specified threshold within a ∆R cone centred on the photon direction. section 2. Uncertainties due to the parton shower model are evaluated by comparing the acceptance factors estimated using MC samples with the default Pythia8 showering with the acceptance factors computed using MC samples relying on the Herwig7 showering model. To account for the uncertainties in the SM Higgs boson production cross-sections when calculating the total acceptance factor from the sum of the various production modes, the fractions of production modes are independently varied within their measured uncertainties taken from ref. [46]. The total systematic uncertainties in the acceptance factors range between 0. 5% and 7%, depending on the observable and bin, with the parton shower uncertainty being the dominant source.
section_name: Event selection

section_content: The inclusive acceptance factors, relative to the full phase space, are about 50% for both the H → ZZ * → 4 and the H → γγ channels. Figure 1 shows the acceptance factors and their systematic uncertainties as a function of p H T and N jets. In the H → ZZ * → 4 channel, the acceptance factor drops in the highest p H T bin due to the lepton separation requirement, while the shape of the acceptance factor for the H → γγ channel as a function of p H T is due to the p T selection criteria on the photons.
section_name: JHEP05(2023)028

section_content: A likelihood combination of the two decay channels is performed, following the method described in ref. [11]. For some observables, such as p H T and p lead. jet T, the binning in the H → γγ analysis is finer than that in the H → ZZ * → 4 analysis. Where needed, the sum of the consecutive H → γγ sub-bins is combined with one H → ZZ * → 4 bin Table 3. Bin boundaries used in the combination of the cross-section of various differential observables. For the p lead. jet T distribution, the first bin contains all events with a leading jet with p T less than 30 GeV, and corresponds exactly with the 0-jet bin in the N jets differential distribution.
section_name: Statistical procedure

section_content: such that the measured bin boundaries match between the two results. A summary of the bin boundaries used in the combined results is presented in table 3. Higgs boson events that are outside of the fiducial region but are reconstructed within the signal region are accounted for in the likelihood function by a small correction (around 1-2%) of the signal normalisation, as described in refs. [5, 7]. 

Experimental and theoretical uncertainties that affect both channels are correlated via common nuisance parameters. The correlated experimental uncertainties include the uncertainties in the integrated luminosity, in the description of the pile-up in the simulation, in the jet reconstruction and calibration, in the common electron-photon energy scale, in the Higgs boson mass value, and in the contributions of the different Higgs boson production modes. Additionally, the common sources of theoretical uncertainty in the H → ZZ * → 4 and H → γγ branching ratios (strong coupling constant, b and c quark masses, and partial decay widths to the main decay channels, such as two vector bosons, two gluons, or a b b pair) are also correlated. Finally, the theoretical uncertainties in the acceptance factor and response matrix due to missing higher-order QCD effects, PDF variations, variations of the modelling of the parton shower, and signal composition uncertainties are also correlated across the Higgs boson decay channels. The asymptotic approximation [47] for the distribution of the profile likelihood ratio is assumed in the computation of uncertainties on all reported measurements. The validity of this approximation has been verified in previous analyses by performing pseudoexperiments.
section_name: JHEP05(2023)028

section_content: The total Higgs boson production cross-section at 13 TeV is measured to be 53. 0 +5. 3 -5. 1 pb ( +4. 9 -4. 8 (stat. ) +2. 0 -1. 7 (syst. )) using the H → ZZ * → 4 decay channel and 58. 1 +5. 7 -5. 4 pb (±4. 2(stat. ) +3. 9 -3. 5 (syst. )) using the H → γγ decay channel. The total cross-section obtained combining the two results is 55. 5 +4. 0 -3. 8 pb (±3. 2(stat. ) +2. 4 -2. 2 (syst. )). All three results are in agreement with the SM prediction of 55. 6 ± 2. 5 pb. The measurements in the two decay channels are compatible with each other with a p-value of 49%, and the compatibility of the combined result with the SM prediction has a p-value of 98%. All compatibility JHEP05(2023)028 Total pp → H + X cross-sections measured at centre-of-mass energies of 7, 8 and 13 TeV, compared with Standard Model predictions taken from ref. [35]. The measurements with the H → ZZ * → 4 channel (blue triangles), H → γγ channel (magenta inverted triangles) and their combination (black dots) are shown. The individual channel results are offset along the x-axis for display purposes. The black boxes around the combined measurements represent the systematic uncertainty, while the error bars show the total uncertainty. The light grey band shows the uncertainty in the prediction due to missing QCD higher-order corrections. The dark grey band indicates the total theoretical uncertainty, corresponding to the dominant QCD higher-ordercorrection uncertainty summed in quadrature with the sum of the PDF and α S uncertainties, and is partially correlated across values of the centre-of-mass energy. checks are performed using a likelihood ratio approach, based on the test statistic variation under different hypotheses in the asymptotic approximation. 

The total cross-section measured using the two channels, their combination, and the SM prediction for a Higgs boson mass of 125. 09 GeV are shown in figure 2. The figure also includes the results of the measurements using data collected at a pp centre-of-mass energies of √ s = 8 TeV and 7 TeV, and the corresponding theoretical expectations. The event samples, selections and the cross-section measurement techniques used for the 8 TeV measurements are described in refs. [48, 49] ; similar techniques are used to measure the cross-sections at 7 TeV as described in refs. [50, 51]. For both the 7 and 8 TeV results, the signal yields in the two decay channels are measured inclusively and corrected for acceptance and detector effects. The results at each centre-of-mass energy are then combined using a likelihood-based technique described in ref. [52]. The total Higgs boson production cross-section at 7 TeV is measured to be 33 +21 -16 pb using the H → ZZ * → 4 channel, 35 +13 -16 pb using the H → γγ decay channel, and 34 +11 -10 pb (±10(stat. ) +4 -2 (syst. )) from their combination. This is to be compared with the SM expectation of 19. 2 ± 0. 9 pb. At 8 TeV, the total Higgs boson production cross-section is measured to be 37 +9 -8 pb using the H → ZZ * → 4 channel, 30. 5 for the individual channels and their combination are shown in figures 3 and 4. The uncertainty band around the SM prediction includes PDF and α S uncertainties as well as those due to missing QCD higher-order corrections, obtained following the method described in ref. [53]. When compared with the results from the individual channels, the total uncertainty for the combined results is lower by 20%-40% and the impact of uncorrelated systematic uncertainties is reduced by approximately 40%. The observed correlation matrices among the cross-sections measured in different bins of the same observable are shown in appendix A. The correlations are small (< 10%) for the Higgs-related observables (p H T, |y H |), characterised by better experimental resolution, and larger (up to about 40%) for jet-related observables (N jets and p lead. jet T ) with worse resolution and larger migrations. All combined measurements are dominated by statistical uncertainties. Significant systematic uncertainties affecting the total and all differential cross-sections arise from the background modelling in the H → γγ signal extraction [7] (typical error of 2-5%) and the integrated luminosity (1. 7%). For the N jets and p lead. jet T differential cross-section measurements, the uncertainties in the reconstruction of the jet energy scale and resolution are important as well, with impacts on the results typically in the range of 2-9%. The dominant theoretical source of uncertainty is the parton shower modelling for ggF signal and has an impact of 2-6%. 

The p-values for the compatibility among the individual measurements are given in table 4. The p-values for the compatibility of the measurements with various theoretical predictions are given in table 5 for the differential cross-section results. For all observables, the measurements in the two channels are compatible with each other, with p-values ranging between 20% and 80%. The combined measurements are also in good agreement with the predictions, with p-values ranging between 20% and 98%. The prediction based on the NNLOPS simulation of gluon-gluon fusion events is lightly favoured over that based on the MG5 FxFx simulation.
section_name: Results

section_content: The observations of the Higgs boson decays to b b [16, 54] provided stringent constraints on the possible modification of the b-quark Yukawa coupling with respect to its SM prediction, whereas current searches for Higgs boson decays to charm final states [17, 55] still allow for a relatively large modification of the c-quark coupling. These measurements have been

Theory/Data (a) Table 5. p-values for the compatibility of the measured cross-sections with the SM predictions when the distributions for gluon-gluon fusion events obtained with either NNLOPS or MG5 FxFx, scaled to the fixed order N 3 LO total gluon-gluon fusion cross-section, are used. The uncertainties in the theoretical predictions are included when calculating the p-values. 

interpreted in terms of the Yukawa coupling modifiers for b-and c-quarks, κ b and κ c, defined as multipliers of the SM values of these couplings [35]. The measured value of κ b agrees with the SM prediction of one with a precision of about 10% [56] to 20% [57], whereas the constraints on κ c are significantly looser: |κ c | < 5. 7 [56] or 1. 1 < |κ c | < 5. 5 [58] at the 95% confidence level (CL). The Higgs boson p T distribution is sensitive to modifications of the Yukawa couplings of the Higgs boson to the b-and c-quarks [59]. This sensitivity is driven by the contributions of b-and c-quarks to the loop-induced ggF production and by the quark-initiated production of the Higgs boson. The former production mode includes an interference term between b-and c-quark loop-mediated amplitudes which is proportional to the product of the two couplings and is therefore sensitive to their relative sign. Modifications of the coupling strength to b-and c-quarks result in changes to both the overall cross-section and the shape of the p H T distribution. In addition, the branching ratio for the H → γγ decay would be affected by corresponding changes to its partial decay width, and both the H → ZZ * → 4 and H → γγ branching ratios would also be affected by the changes to the total Higgs boson decay width.
section_name: Constraints on the band c-quark Yukawa couplings

section_content: This section presents constraints on κ b and κ c, inferred from the measured p H T distributions. A combined interpretation is then performed in terms of κ b and κ c by including also the constraints from the measurement of Higgs bosons, produced in association with a vector boson, decaying to b b [16], and from the search for Higgs bosons produced in a similar way and decaying to cc [17]. All tree-level couplings of the Higgs boson to particles other than the b-or c-quarks are set to their SM values and loop-induced Higgs boson couplings are resolved to their SM expectation, with κ b and κ c as free parameters.
section_name: JHEP05(2023)028

section_content: The constraints on κ b and κ c from the observed p H T distributions in the H → ZZ * → 4 and H → γγ final states are derived in two scenarios: one in which only modification to the shape of the measured p H T distributions is considered ("shape-only"), and one in which the impact on the overall expected normalisation, through modifications of the total cross-sections, branching ratios and Higgs boson decay width, is also considered ("shape+normalisation"). 

The theoretical predictions used for these interpretations are detailed in ref. [7]. The predictions for κ b and κ c modifications of the ggF production are computed with SCETlib [60, 61]. For the H → γγ decay channel, the calculations are performed directly after applying the particle-level requirements defining the H → γγ fiducial phase space, while for the H → ZZ * → 4 decay channel, the calculations are performed in the full phase space and then extrapolated to the H → ZZ * → 4 fiducial phase space using the acceptance factors obtained from the NNLOPS ggF prediction. It has been verified that the dependence of the acceptance factors in each p H T bin on the b and c Yukawa coupling modifiers is negligible. The predictions for quark-initiated b b → H and cc → H production modes are computed with MadGraph5_aMC@NLO 2. 7. 3. The simulation of the Higgs boson decay, the parton shower, hadronisation and underlying event, is performed with Pythia 8 using a dedicated PDF set [62] and the A14 tune. The inclusive b b → H and cc → H cross-sections are then normalised to the state-of-the-art NNLO computations available in refs. [62, 63]. All the other Higgs boson production modes remain unchanged with κ b and κ c variations, and they are estimated as detailed in section 2. 

Theoretical uncertainties related to the QCD modelling and PDF uncertainties on the differential cross-sections are considered using the procedure detailed in ref. [7]. For the theoretical calculation of the ggF process, uncertainties related to numerical integration, fixed-order scale, hard resummation phase, resummation scheme, matching scale and nonperturbative scheme are implemented [60]. For the b-and c-quark initiated processes, the uncertainty related to missing higher order QCD effects are estimated by varying the renormalisation, factorisation and merging scales; the uncertainties related to the PDF set are estimated by varying the mass and scale associated with the b-quark for the b b → H process and by using the MC replicas of the nominal PDF set for the cc → H process; the uncertainties due to the choice of the FxFx merging scale are estimated by using alternative values of this scale. Theoretical uncertainties in the other production modes that do not depend on κ b or κ c, from higher-order QCD effects, PDF and α s, and the parton shower model, are estimated as described in section 2.
section_name: Constraints from the Higgs boson transverse momentum distributions

section_content: Channel Parameter Observed Expected 95% confidence interval 95% confidence interval Table 6. Observed and expected 95% confidence intervals for the Yukawa coupling modifiers when modifications to only the p H T shape are considered (shape-only), for the individual decay channels and their combination. The results for one coupling modifier are obtained while fixing the other one to the SM expectation (κ = 1). 

The statistical interpretation is performed by first parameterising the fiducial crosssections as a function of κ b and κ c for each decay channel. The two likelihood models are then jointly interpreted using the same procedure as detailed in section 4. 

The expected and observed 95% confidence intervals for κ b and κ c are shown in tables 6 and 7 for the shape-only and shape+normalisation scenarios, respectively. The limits on a given κ parameter are determined with the other one fixed to SM prediction (κ = 1). If κ b is unconstrained in the fit, the 95% confidence intervals for κ c are about 10% (twice) larger than if κ b is fixed to the SM value of one, in the shape-only (shape+normalisation) approach. 

In the shape-only approach, the combined observed limits are less stringent than the individual H → ZZ * → 4 result. This is due to the quadratic dependency of the crosssection and the differential distribution on the κ parameters leading to a double minimum in the profile likelihood ratio, and due to the combined best-fit value for the κ b parameter being further from the SM expectation when probing only the H → γγ decay channel. For κ c, the observed combined best-fit value is similar to the best-fit value in the H → ZZ * → 4 channel. However, due to the correlation between the κ b and κ c parameters, different best-fit κ b observations between the channels, as well as the data fluctuations in some of the p H T bins, the 68% CL observed combined limits on κ c are worse than the results from the H → ZZ * → 4 channel. The corresponding 95% CL limits are similar to those from the H → ZZ * → 4 channel. 

In the shape+normalisation scenario the constraints on the coupling modifiers are tighter, since a large fraction of the allowed ranges for κ b and κ c from the shape-only approach lead to values of the total width and thus of the H → ZZ * → 4 and H → γγ branching ratios and overall normalisation that are inconsistent with the data. Two-dimensional confidence regions on κ b and κ c are also derived for both scenarios, as shown in figure 5.
section_name: JHEP05(2023)028

section_content: Channel Parameter Observed Expected 95% confidence interval 95% confidence interval 7. Observed and expected 95% confidence intervals for the Yukawa coupling modifiers when modifications to both the p H T shape and normalisation are considered (shape+normalisation), for the individual decay channels and their combination. The results for one coupling modifier are obtained while fixing the other one to the SM expectation (κ = 1). T differential cross-section (shape-only) or (b) also its normalisation (shape+normalisation) is used to constrain the parameters for the combined and individual decay channels results. The SM predictions ( * ) and the observed best-fit values (+) are indicated on the plots.
section_name: JHEP05(2023)028

section_content: The measurement of Higgs boson decays to b b and the search for Higgs boson decays to cc in Higgsstrahlung events (V H) constrain the b-and c-quark coupling modifiers through the quadratic dependence on κ 2 b and κ 2 c of the partial widths of the Higgs boson to these two final states. This section describes the methodology and the results of a simultaneous fit to the Higgs boson transverse momentum distributions of the H → ZZ * → 4 and H → γγ fiducial cross section measurements and to the multivariate discriminant used to measure the V H(q q) (q = b, c) signal strength [16, 17]. 

Two scenarios are considered for this combination. The first scenario is the "shape+normalisation" scenario as described previously. In the second scenario, the JHEP05(2023)028 Higgs boson is also allowed to decay to BSM particles and the associated partial width is included in the total width. The partial width for BSM decays is parameterised as

where Γ is the Higgs boson total width, and B BSM is its branching ratio to BSM particles. The second scenario reduces the assumptions of the model, at the cost of reduced sensitivity. 

In the combination, most common experimental systematic uncertainties and signal theory uncertainties are modelled as correlated between the four channels (H → ZZ * → 4, H → γγ, V H(b b), V H(cc)). Jet energy calibration and flavour tagging efficiency uncertainties are not modelled as correlated between the channels due to the use of different jet clustering algorithms. 

The observed 68% and 95% CL contours in the 2D κ b vs κ c plane are shown in figure 6(a) for the shape+normalisation scenario where B BSM is fixed to zero and in figure 6(b) for the case where B BSM is a free parameter. The fit prefers a positive value of κ b, but negative values are not excluded at 68% CL, leading to two disconnected allowed regions, corresponding to positive or negative values of κ b. One-dimensional confidence intervals for κ c with κ b unconstrained in the fit are summarised in table 8. Excluding the V H(cc) channel would worsen the one-dimensional constraints on κ c by about 10% for the B BSM = 0 scenario, and by a factor two for the alternative scenario where B BSM is not fixed to zero.
section_name: Combination with the constraints from V H(b b) and V H(cc) production

section_content: A combined measurement of the total and differential Higgs production cross-sections in the H → γγ and H → ZZ * → 4 decay channels was performed using 139 fb -1 of 13 TeV proton-proton collision recorded by the ATLAS detector during the LHC Run 2. Good agreement is observed when comparing the results from the two channels, after -15 - Differential cross-sections are measured as a function of the Higgs boson transverse momentum and rapidity, the number of jets produced together with the Higgs boson and the transverse momentum of the leading jet. The larger data set and the combination of the two decay channels result in measurement uncertainties that are significantly smaller than in previous results. Notably, the differential cross-section as a function of the Higgs boson transverse momentum is measured with 20-30% precision up to 300 GeV and about 60% precision in the 300-650 GeV range. The combined differential distributions agree with the Standard Model predictions.
section_name: Conclusions

section_content: The measured fiducial differential cross-sections as a function of p H T are used to derive limits on the bottom and charm-quark Yukawa couplings modifiers, κ b and κ c, assuming SM values of the other tree-level Higgs boson couplings. Fixing the value of κ b to one, the 95% confidence interval for κ c is [-8. 6, 17. 3] using only the observed shape of the p H T distribution, and [-2. 27, 2. 27] when considering also the impact of these couplings on the normalisation of the measured p H T fiducial cross-sections. A combined fit with the ATLAS measurement of Higgs bosons produced in association with a W or Z boson and decaying to b-or c-quark pairs allows constraints to be set on the charm quark coupling modifier without any assumption on the bottom quark coupling. The 95% CL allowed range for κ c when the Higgs boson is assumed to decay only to SM particles is [-2. 47, 2. 53] while in a more generic scenario in which BSM Higgs boson decays are allowed, the constraint is loosened to [-4. 46, 4. 81 ]. These represent the most stringent constraints on κ c to date in these scenarios. The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN, the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA), the Tier-2 facilities worldwide and large non-WLCG resource providers. Major contributors of computing resources are listed in ref. [64].
section_name: JHEP05(2023)028

section_content: -17 - 3, with a higher label index corresponding to a higher bin for the given variable. For the correlation matrix for the Higgs boson transverse momentum vs Higgs boson rapidity, lower rapidity bins are labelled first with ascending bins in p H T.
section_name: JHEP05(2023)028

section_content: -18 -JHEP05(2023)028 3, with a higher label index corresponding to a higher bin for the given variable. 

Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4. 0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited. SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.
section_name: A Correlation matrices between the measured cross-sections

section_name: Acknowledgments

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

We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia ; BMWFW and FWF, Austria; ANAS, Azerbaijan; CNPq and FAPESP, Brazil ; NSERC, NRC and CFI, Canada ; CERN; ANID, Chile ; CAS, MOST and NSFC, China ; Minciencias, Colombia; MEYS CR, Czech Republic; DNRF and DNSRC, Denmark ; IN2P3-CNRS and CEA-DRF/IRFU, France; SRNSFG, Georgia; BMBF, HGF and MPG, Ger-
section_name: Acknowledgments