Combined measurement of differential and total cross sections in the H → γ γ and the H → Z Z ∗ → 4 (cid:3) decay channels at √ s = 13 TeV with the ATLAS detector

A combined measurement of differential and inclusive total cross sections of Higgs boson production is performed using 36.1 fb − 1 of 13 TeV proton–proton collision data produced by the LHC and recorded by the ATLAS detector in 2015 and 2016. Cross sections are obtained from measured H → γγ and H → Z Z ∗ → 4 (cid:3) event yields, which are combined taking into account detector eﬃciencies, resolution, acceptances and branching fractions. The total Higgs boson production cross section is measured to be 57.0 + 6 . 0 − 5 . 9 (stat.) + 4 . 0 − 3 . 3 (syst.) pb, in agreement with the Standard Model prediction. Differential cross- section measurements are presented for the Higgs boson transverse momentum distribution, Higgs boson rapidity, number of jets produced together with the Higgs boson, and the transverse momentum of the leading jet. The results from the two decay channels are found to be compatible, and their combination agrees with the Standard Model predictions. 3 .

This Letter describes the combination of two fiducial crosssection measurements in the H → γ γ [11] and H → Z Z * → 4 [12] decay channels, which were obtained using 36.1 fb −1 of pp collision data produced by the Large Hadron Collider (LHC) in 2015 and 2016 with a centre-of-mass energy of 13 TeV and recorded by the ATLAS detector [13]. The combined cross section is extracted for the total phase space, increasing the degree of model dependence compared to the individual measurements, which were performed in a fiducial phase space close to the selection criteria for reconstructed events in the detector. Despite the additional systematic uncertainties assigned to the extrapolation to the total phase space, the combination significantly reduces the measurement uncertainty compared to the results in the individual decay channels.
The measured observables include the total production cross section, the Higgs boson's transverse momentum p H T , sensitive to perturbative QCD calculations, and the Higgs boson's rapidity |y H |, sensitive to the parton distribution functions (PDF). Furthermore the number of jets N jets is measured in events with a Higgs boson and jet transverse momentum above 30 GeV, as well as the leading jet's transverse momentum p j1 T . Both the N jets and p j1 T observables probe the theoretical modelling of high-p T QCD radiation in Higgs boson production. The N jets observable is also sensitive to the different Higgs boson production processes [14].
The cross sections are obtained from yields measured in the H → γ γ and H → Z Z * → 4 decay channels, which are combined taking into account detector efficiencies, resolution, acceptances and branching fractions. For each decay channel and each observable, the cross sections can be written as where i is the iterator over the bins of the observable of interest, σ i is the cross section in bin i, N sig i is the number of measured reconstructed signal events following the analysis selection, L is the integrated luminosity and B is the branching fraction. The term C i is the correction factor from the number of events reconstructed to the number of events at particle level produced in the respective fiducial phase space, and A i is the acceptance factor extrapolating  acceptance factors and uncertainties are calculated for the combination, as discussed in Section 3. Section 4 presents the combination methodology. The results are discussed in Section 5.

Higgs boson Monte Carlo samples, cross sections and branching fractions
Predictions of SM Higgs boson production are used in the calculation of the correction and acceptance factors, and are compared to the measured cross sections. The Monte Carlo (MC) event generators that were used to simulate gluon-gluon fusion (ggF), vector-boson fusion (VBF), associated Higgs boson production (V H, V = W , Z ), and Higgs boson production in association with a heavy-quark pair (tt H, bbH) are listed in Table 1. The accuracy of the calculations and the PDF sets used are also given, with the abbreviations NLO for next-to-leading order, NNLO for next-to-nextto-leading order, and NNLL for next-to-next-to-leading logarithm. For ggF, VBF, V H, bbH in both decay channels and tt H in the H → γ γ decay channel, Pythia8 [16,17] [34][35][36]. These predictions were obtained assuming a Higgs boson mass of 125.09 GeV [15] to calculate cross sections and branching ratios. Details are given in Table 2, including the accuracy of the calculations, and the composition of the production modes in the SM. N 3 LO is the abbreviation for next-tonext-to-next-to-leading order, and EW stands for electroweak.
In addition to the NNLOPS sample (see Table 1) scaled to the N 3 LO cross section with a K -factor of 1.1, further SM ggF predictions are compared with the measurements. If not mentioned otherwise, the cross sections predicted by the respective calculations are used. For the comparison with data, the non-ggF Higgs boson production processes are added using the samples and cross sections described above.

T distribution is compared with the predictions from
HRes [63,64], RaDISH + NNLOJET [65], and Madgraph5_aMC@ NLO. HRes includes resummation to NNLL and computes fixedorder cross sections for ggF Higgs boson production up to NNLO in QCD. It describes the p H T distribution at NLO. Finite t-, b-, and c-quark masses are included at NLO accuracy. The RaDISH + NNLOJET prediction includes resummation to NNLL and matching to the one-jet NNLO differential spectrum from NNLOJET [66,67]. It includes corrections from the finite t-and b-quark masses. The predictions from Madgraph5_aMC@NLO are scaled to the N 3 LO cross section with a K -factor of 1. 47. This generator provides NLO accuracy in QCD for zero, one, and two additional jets, merged with the FxFx scheme [68] and includes the finite top quark mass effects [30,69,70].
• The |y H | measurement is compared with predictions from Madgraph5_aMC@NLO merged with the FxFx scheme and SCETlib + MCFM8 [71,72], which achieves NNLO + NNLL ϕ accuracy 1 by applying a resummation of the virtual corrections to the gluon form factor. The underlying NNLO predictions are obtained using MCFM8 with zero-jettiness subtractions [73,74].
• Multiple predictions exist for different bins of the N jets distribution. Considered here are the STWZ-BLPTW prediction [14, 75,76], which includes NNLL +NNLO resummation for the p T of the leading jet, combined with a NLL +NLO resummation for the subleading jet, and the JVE-N 3 LO prediction [77], which includes NNLL resummation of the p T of the leading jet with small-R resummation and is matched to the N 3 LO total cross section. In addition, predictions from Madgraph5_aMC@NLO, are compared with the full N jets distribution.
For ggF, VBF and V H, the PDF4LHC set is varied according to its eigenvectors [25], and the envelope of the variations is used as the systematic uncertainty. The effect of PDF uncertainties on tt H and bbH is negligible and not included. The renormalization and fac- 1 The prime indicates that important parts of the N 3 LL (next-to-next-to-next-toleading logarithm) contribution are included along with the full NNLL corrections and the subscript ϕ indicates that resummation is applied to the gluon form factor.

Acceptance correction
The acceptance factors that extrapolate at particle-level from the H → γ γ and H → Z Z * → 4 fiducial phase space to the full phase space are estimated using the MC samples and cross sections described in Section 2. Their evaluation assumes SM Higgs boson production fractions and a Higgs boson mass of 125 GeV; the 90 MeV difference from 125.09 GeV has negligible impact on the Higgs boson kinematics and is covered by the systematic uncertainty from the Higgs boson mass measurement.
In the H → γ γ fiducial phase space [11], the selected events have two photons with pseudorapidity 3 |η| < 1.37 or 1.52 < |η| < 2.37 and p T the transverse momentum of the (sub)leading photon and m γ γ is the invariant mass of the two photons. The photons are required to be isolated: the p T of the system of charged generator-level particles within R = 0.2 of the photon is required to be less than 0.05 times the p T of the photon. In the H → Z Z * → 4 fiducial phase space [12], the selected events have four muons, four electrons, or two electrons and two muons. The three leading leptons are required to have p T > 20, 15, 10 GeV. The lowest-p T muon (electron) has to fulfil p T > 5 (7) GeV. The muons have to be within |η| < 2.7 and the electrons within |η| < 2.47. Following the selection of events in data, requirements are placed on the masses of the two same-flavour opposite-charge pairs, on the R of any two leptons, and the invariant mass of the four-lepton system, In the total 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. Simulated 3 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upwards. Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the z-axis. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2). Angular distance is measured in units of Theory uncertainties in the signal acceptance related to the PDF, higher-order corrections, and the parton shower are considered for the acceptance factors and are correlated between the two channels. Uncertainties due to the PDF and scales are estimated as described in Section 2. Uncertainties due to the parton shower are evaluated by comparing the ggF default showering Pythia8 with Herwig7. The uncertainty is derived from the full difference between the two cases. The Higgs boson mass is varied within the uncertainty of the ATLAS-CMS combined measurement [15]. To account for model dependence, the fractions of production modes are varied within the uncertainties from the dedicated measurements by the ATLAS and CMS collaborations [81]. For tt H, the 13 TeV ATLAS results are used [82]. The bbH cross section is varied within the uncertainties due to the PDF and higher-order corrections [14]. The total systematic uncertainties of the acceptance factors range between 0.4% and 5%, depending on the observable and bin. The parton shower uncertainty dominates.
The inclusive acceptance factors are 50% for the H → γ γ channel and 42% for the H → Z Z * → 4 channel (relative to the full phase space of H → Z Z * → 2 2 , where , = e or μ). The acceptance is lower for H → Z Z * → 4 than for H → γ γ since it is less likely for four leptons to fulfil the fiducial requirements. Fig. 1 shows the acceptance factors used for the differential observables and their systematic uncertainties. The fiducial acceptance falls off steeply as the Higgs boson rapidity increases, as both fiducial definitions include pseudorapidity requirements on the Higgs boson decay products. The fiducial acceptance in the H → γ γ channel as a function of p H T is shaped by the p T selection criteria on the photons.

Statistical procedure
The combined measurement is based on maximizing the profile-likelihood ratio [83]: Here σ are the parameters of interest, θ are the nuisance parameters, and L represents the likelihood function. The σ and θ terms denote the unconditional maximum-likelihood estimate of the parameters, while θ (σ ) is the conditional maximum-likelihood estimate for given parameter values.
The likelihood function L includes the signal extraction, the correction to particle level, and the extrapolation to the total phase space in each channel. Therefore, the total cross section as well as the cross sections in different bins for each observable can be derived directly as parameters of interest σ based on the combined data set from the H → γ γ and H → Z Z * → 4 channels.
The distribution shape and normalization systematic uncertainties of all components are included in the likelihood function as nuisance parameters θ with constraints from subsidiary measurements. This allows the uncertainties to be correlated between bins, decay channels, and correction and acceptance factors. The uncertainty components of the predicted branching ratios are correlated between the decay channels, as well as the uncertainties in the acceptance and correction factors due to production mode variations, PDF and higher-order corrections, and the parton shower. The uncertainty in the Higgs boson mass, including its effect on the predicted branching ratio, is also correlated between channels. Experimental uncertainties in the correction factors and the signal extraction in the H → Z Z * → 4 decay channel, like the energy scale and resolution of electrons, photons, and jets, and in the luminosity measurement and pileup modelling are also correlated. Where one bin in one of the measurements corresponds to two bins in the other, the wider bin size is used. The sum of the cross sections in the finer bins is considered as the parameter of interest in these cases, and an additional unconstrained nuisance parameter that floats in the fit describes the difference between the merged bins. The normalization and shape uncertainties of the H → γ γ background estimate [11] are fit to the data as nuisance parameters without any initial constraint.
The test statistic −2 ln is assumed to follow a χ 2 distribution for constructing confidence intervals [83]. This asymptotic assumption was tested with pseudo-experiments for bins with low numbers of events and found to be appropriate.
The level of agreement between the two channels in the total phase space is evaluated by using a profiled likelihood as a function of the difference of the cross sections in each bin i, σ i γ γ − σ i 4 . The number of degrees of freedom is the same as the number of bins in the tested distribution. The probability that a measured differential cross section is compatible with a theoretical prediction is found by computing a p-value based on the difference between the value of −2 ln at the unconditional maximum-likelihood estimate and the value obtained by fixing the cross sections in all bins to the ones predicted by the theory. The uncertainties in the theoretical predictions are ignored when calculating the p-values. Including these uncertainties would increase the p-values.

Results
The total cross section is measured to be 47. For all differential observables and bins, the measurement is dominated by statistical uncertainties, which vary between 20% and 30%. Significant uncertainties affecting all observables, includ- The level of agreement between the two channels in the total phase space is quantified by the corresponding p-values: 58% for p H T , 40% for |y H |, 53% for N jets and 67% for p j1 T . Table 3 shows the p-values indicating reasonable agreement between the probed SM predictions and the measurement. The relatively low p-value for HRes can be explained by the lower computed total cross section, as this prediction is at NNLO + NNLL accuracy. The lower p-values for p j1 T reflect the lower predictions compared to the measurement for high jet p T . Compatibility checks of individual bins indicate less than 3σ local discrepancy.

Conclusion
A combined measurement of the total and differential cross sections in the H → γ γ and H → Z Z * → 4 decay channels was performed, using 36.1 fb −1 of 13 TeV proton-proton collision data produced by the LHC and recorded by the ATLAS detector in 2015 and 2016. Good agreement is observed when comparing the re- Table 3 p-values in percent indicating the probabilities that the measured differential cross sections are compatible with various SM ggF predictions. The NNLOPS and Madgraph5_aMC@NLO predictions are scaled to the total N 3 LO cross section by the given K -factors. The non-ggF predictions are added, as discussed in Section 2. The uncertainties in the theoretical predictions are ignored when calculating the p-values.