The post-fit distribution of the $M(\ell\ell)$ variable is shown for the high-MET bin for the $\text{t}\bar{\text{t}}$ CR. Uncertainties include both the statistical and systematic components.

The post-fit distribution of the $M(\ell\ell)$ variable is shown for the low-MET bin for the WZ-enriched region. Uncertainties include both the statistical and systematic components.

The post-fit distribution of the $M(\ell\ell)$ variable is shown for the high-MET bin for the WZ-enriched region. Uncertainties include both the statistical and systematic components.

The post-fit distribution of the $M(\ell\ell)$ variable is shown for the high-MET bin for the SS CR. Uncertainties include both the statistical and systematic components.

The 2$\ell$-Ewk SR: the post-fit distribution of the $M(\ell\ell)$ variable is shown for the low-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 2$\ell$-Ewk SR: the post-fit distribution of the $M(\ell\ell)$ variable is shown for the med-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 2$\ell$-Ewk SR: the post-fit distribution of the $M(\ell\ell)$ variable is shown for the high-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 2$\ell$-Ewk SR: the post-fit distribution of the $M(\ell\ell)$ variable is shown for the ultra-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 3$\ell$-Ewk search regions: the post-fit distribution of the $M^{\text{min}}_{\text{SFOS}}(\ell\ell)$ variable is shown for the low-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 3$\ell$-Ewk search regions: the post-fit distribution of the $M^{\text{min}}_{\text{SFOS}}(\ell\ell)$ variable is shown for the high-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the TChiWZ and the simplified higgsino models in the scenario where the product of $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ eigenvalues is positive and negative, respectively. The numbers after the model name in the legend indicate the mass of the NLSP and the mass splitting between the NLSP and LSP, in GeV.

The 2$\ell$-Stop SR: the post-fit distribution of the leading lepton $p_{T}$ variable is shown for the low-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the T2bff$\tilde{\chi}^0_1$ and the T2bW models. The numbers after the model name in the legend indicate the mass of the top squark and the mass splitting between the top squark and LSP, in GeV.

The 2$\ell$-Stop SR: the post-fit distribution of the leading lepton $p_{T}$ variable is shown for the med-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the T2bff$\tilde{\chi}^0_1$ and the T2bW models. The numbers after the model name in the legend indicate the mass of the top squark and the mass splitting between the top squark and LSP, in GeV.

The 2$\ell$-Stop SR: the post-fit distribution of the leading lepton $p_{T}$ variable is shown for the high-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the T2bff$\tilde{\chi}^0_1$ and the T2bW models. The numbers after the model name in the legend indicate the mass of the top squark and the mass splitting between the top squark and LSP, in GeV.

The 2$\ell$-Stop SR: the post-fit distribution of the leading lepton $p_{T}$ variable is shown for the ultra-MET bin. Uncertainties include both the statistical and systematic components. The signal distributions overlaid on the plot are from the T2bff$\tilde{\chi}^0_1$ and the T2bW models. The numbers after the model name in the legend indicate the mass of the top squark and the mass splitting between the top squark and LSP, in GeV.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the TChiWZ model. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. Results are reported for the $\widetilde{m}_{\tilde{\chi}^0_2} \widetilde{m}_{\tilde{\chi}^0_1} > 0$ $M(\ell\ell)$ spectrum reweighting scenario. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the TChiWZ model. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. Results are reported for the $\widetilde{m}_{\tilde{\chi}^0_2} \widetilde{m}_{\tilde{\chi}^0_1} < 0$ $M(\ell\ell)$ spectrum reweighting scenario. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the simplified higgsino model. The simplified model includes both neutralino pair and neutralino-chargino production modes. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. The results are reported for the $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ $M(\ell\ell)$ spectrum reweighting scenario. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the pMSSM higgsino model. The pMSSM one includes all possible production modes. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. The results are reported for the $\widetilde{m}_{\tilde{\chi}^0_2}\widetilde{m}_{\tilde{\chi}^0_1}$ $M(\ell\ell)$ spectrum reweighting scenario. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the T2bff$\tilde{\chi}^0_1$ simplified model. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

The observed 95% CL exclusion contours (black curves) assuming the NLO+NLL cross sections, with the variations (thin lines) corresponding to the uncertainty in the cross section for the T2bW simplified model. The red curves present the 95% CL expected limits with the band (thin lines) covering 68% of the limits in the absence of signal. The range of luminosities of the analysis regions included in the fit is indicated on the plot.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Comparison of predicted and observed event yields in each of the control and signal regions in the dilepton channel after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison of predicted and observed event yields in each of the control and signal regions in the single-lepton channels after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV (yield only). The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ lo}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ hi}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit yields of signal ($S$) and total background ($B$) as a function of $\log (S/B)$, compared with data. Final-discriminant bins in all dilepton and single-lepton analysis regions are combined into bins of $\log (S/B)$, with the signal normalised to the SM prediction used for the computation of $\log (S/B)$. The signal is then shown normalised to the best-fit value and the SM prediction. The lower frame reports the ratio of data to background, and this is compared with the expected ${t\bar {t}H}$-signal-plus-background yield divided by the background-only yield for the best-fit signal strength (solid red line) and the SM prediction (dashed orange line).

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta \eta $ between $b$-tagged jets, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the likelihood discriminant, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta R$ for all possible combinations of $b$-tagged jet pairs, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate mass for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Fitted values of the ${t\bar {t}H}$ signal strength parameter in the individual channels and in the inclusive signal-strength measurement.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the fit. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Pre-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Post-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal strength.

95% CL simplified template cross-section upper limits in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive limit. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown. The hatched uncertainty bands correspond to the theory uncertainty in the fiducial cross-section prediction in each bin.

The ratios $S/B$ (black solid line, referring to the vertical axis on the left) and $S/\sqrt{B}$ (red dashed line, referring to the vertical axis on the right) for each category in the inclusive analysis in the dilepton channel (left) and in the single-lepton channels (right), where $S$ ($B$) is the number of selected signal (background) events predicted by the simulation and normalised to a luminosity of 139 fb$^{-1}$ .

Comparison between data and prediction for the $\Delta R$ between the Higgs candidate and the ${t\bar {t}}$ candidate system, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the number of $b$-tagged jet pairs with an invariant mass within 30 GeV of 125 GeV, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the $\Delta R$ between the two highest ${p_{{T}}}$ $b$-tagged jets, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $0\le {\hat{p}_{{T}}^{H}}<120$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $120\le {\hat{p}_{{T}}^{H}}<200$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $200\le {\hat{p}_{{T}}^{H}}<300$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $300\le {\hat{p}_{{T}}^{H}}<450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for ${\hat{p}_{{T}}^{H}}\ge 450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

95% confidence level upper limits on signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal-strength limit, after the fit used to extract multiple signal-strength parameters. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown.

Post-fit correlation matrix (in percentages) between the $\mu $ values obtained in the STXS bins.

Performance of the Higgs boson reconstruction algorithms. For each row of `truth' ${\hat{p}_{{T}}^{H}}$, the matrix shows (in percentages) the fraction of Higgs boson candidates which are truth-matched to ${b\bar {b}}$ decays, with reconstructed $p_T^H$ in the various bins of the dilepton (left), single lepton resolved (middle) and boosted (right) channels.

Pre-fit event yields in the dilepton signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Post-fit event yields in the dilepton signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Pre-fit event yields in the single-lepton resolved and boosted signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Post-fit event yields in the single-lepton resolved and boosted signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Breakdown of the contributions to the uncertainties in $\mu$. The contributions from the different sources of uncertainty are evaluated after the fit. The $\Delta \mu $ values are obtained by repeating the fit after having fixed a certain set of nuisance parameters corresponding to a group of systematic uncertainties, and then evaluating $(\Delta \mu)^2$ by subtracting the resulting squared uncertainty of $\mu $ from its squared uncertainty found in the full fit. The same procedure is followed when quoting the effect of the ${t\bar {t}+{\geq }1b}$ normalisation. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters existing in the fit.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the dilepton channel.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the single-lepton channels.

Predicted SM ${t\bar {t}H}$ cross-section in each of the five STXS ${\hat{p}_{{T}}^{H}}$ bins and signal acceptance times efficiency (including all event selection criteria) in each STXS bin as well as for the inclusive ${\hat{p}_{{T}}^{H}}$ range.

Number of expected signal events before the fit, after each selection requirement applied to enter the dilepton channel $SR^{\geq 4j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel resolved $SR^{\geq 6j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel boosted $SR_{boosted}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

The JER for R = 0.2 jets in Pb+Pb collisions as a function of $2|\Psi_2-\phi^{reco}|$ for centrality selections of 0-5%, 5-10%, 10-20%, 20-40% and 40-60%.

The systematic uncertainties in v2 for 20-40% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for 5-10% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 20-40% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 5-10% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 20-40% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 5-10% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for pT = 71--398 GeV jets as a function of centrality.

The systematic uncertainties in v3 for pT = 71--398 GeV jets as a function of centrality.

The systematic uncertainties in v4 for pT = 71--398 GeV jets as a function of centrality.

Angular distribution of jets as a function of the observed psi2 plane for jets with 71 < pT < 79 GeV in the 10-20% centrality bin.

Angular distribution of jets as a function of the observed psi3 plane for jets with 71 < pT < 79 GeV in the 10-20% centrality bin.

Angular distribution of jets as a function of the observed psi4 plane for jets with 71 < pT < 79 GeV in the 10-20% centrality bin.

The v2 values for R = 0.2 jets as a function of centrality for jets in several pT ranges.

The v2 values for R = 0.2 jets as a function of pT for 0-5%, 5-10%, and 20-40% centrality collisions.

The v2, v3, and v4 as a function of centrality for jets with pT = 71-398 GeV.

The v3 values for R = 0.2 jets as a function of centrality for jets in several pT ranges.

The v4 values for R = 0.2 jets as a function of centrality for jets in several pT ranges.

The v2 as a function of pT for jets in 10-20% centrality collisions.

The v3 as a function of pT for jets in 10-20% centrality collisions.

The v2 as a function of pT for jets in 20-40% centrality collisions.

The v3 as a function of pT for jets in 20-40% centrality collisions.

The v2 for jets in 10-20% centrality collisions.

The v3 for jets in 10-20% centrality collisions.

R2max as a function of pT (filled circles). Also shown is 1 - 4v2/(1+2v2).

R3max as a function of pT (filled circles). Also shown is 1 - 4v3/(1+2v3).

The systematic uncertainties in v2 for 40-60% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 40-60% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 40-60% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for 10-20% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 10-20% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 10-20% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for 0-5% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 0-5% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 0-5% centrality Pb+Pb collisions as a function of pT.

The v2 as a function of pT for jets in 40-60% centrality collisions.

The v2 as a function of pT for jets in 5-10% centrality collisions.

The v2 as a function of pT for jets in 0-5% centrality collisions.

The v3 as a function of pT for jets in 40-60% centrality collisions.

The v3 as a function of pT for jets in 5-10% centrality collisions.

The v3 as a function of pT for jets in 0-5% centrality collisions.

The efficiency for simulated VBF events to pass each analysis cut.

The efficiency for simulated VH events to pass each analysis cut.

The efficiency for simulated ttH events to pass each analysis cut.

The fractional contribution of each production mode to a given analysis region around the Higgs boson peak, defined as 105 < mJ < 140 GeV. The fraction is given with respect to the total signal yield in the analysis region in question.

Signal acceptance times efficiency within the fiducial volume used in the fiducial region.

Signal acceptance times efficiency for the STXS volumes in the differential measurement. Along with the pTH requirements shown, |yH | < 2 is required. For events with pTH < 300 GeV, the acceptance times efficiency is less than 0.1 × 10−2.

Event yields and associated uncertainties after the global likelihood fit in the pT(H) > 450 GeV fiducial region. The total background yield can differ from the sum of the individual components due to rounding.

Event yields and associated uncertainties after the global likelihood fit of the differential regions. The total background yield can differ from the sum of the individual components due to rounding. The five STXS volumes are labeled pT0–pT4 corresponding to pTH < 300 GeV, 300 – 450 GeV, 450 – 650 GeV, 650 – 1000 GeV, and > 1000 GeV, respectively. The pT0 event yield is constrained to its SM value within the theoretical and experimental uncertainties and free parameters act independently on the remaining four volumes.

The observed exclusion contours as a function of (m(med), m($\chi$)), with Vector mediator)

The expected exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.35)

The observed exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.35)

The expected exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.7)

The observed exclusion contours as a function of (m(a), tan($\beta$)), with sin($\theta$) = 0.7)

The expected exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.35)

The observed exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.35)

The expected exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.7)

The observed exclusion contours as a function of (m(H), tan($\beta$)), with sin($\theta$) = 0.7)

The expected exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.35)

The observed exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.35)

The expected exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.7)

The observed exclusion contours as a function of (m(a), m(H)), with sin($\theta$) = 0.7)

Expected lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 200 GeV, m(H) = 600 GeV.

Observed lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 200 GeV, m(H) = 600 GeV.

Expected lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 350 GeV, m(H) = 1000 GeV.

Observed lower limit on signal strength at 95% CL as a function of sin($\theta$), with m(a) = 350 GeV, m(H) = 1000 GeV.

Observed lower limit on WIMP-nucleon cross section at 90% CL as a function of m(WIMP), assuming Higgs-portal scenario with Scalar WIMP.

Observed lower limit on WIMP-nucleon cross section at 90% CL as a function of m(WIMP), assuming Higgs-portal scenario with Majorana WIMP.

Observed lower limit on the spin-dependent WIMP–proton scattering cross-section.

Observed lower limit on the spin-independent WIMP–nucleon scattering cross-section.

Cutflow of unweighted and weighted events of $ZH$ signals in the electron channel.

Cutflow of unweighted and weighted events of $ZH$ signals in the muon channel.

Cutflow of unweighted and weighted events of DM signals (simplified DM axial-vector with m(med) = 900 GeV and m($\chi$) = 40 GeV, 2HDM+$a$ signal with m(A) = 600 GeV, m(a) = 400 GeV, tan($\beta$) = 1.0 and m($\chi$) = 10 GeV) in the electron channel.

Cutflow of unweighted and weighted events of DM signals (simplified DM axial-vector with m(med) = 900 GeV and m($\chi$) = 40 GeV, 2HDM+$a$ signal with m(A) = 600 GeV, m(a) = 400 GeV, tan($\beta$) = 1.0 and m($\chi$) = 10 GeV) in the muon channel.

The inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ measured in $1.5 < |y| < 2.0$ for jets clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.4$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ measured in $|y| < 0.5$ for jets clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$.

The inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ measured in $0.5 < |y| < 1.0$ for jets clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$.

The inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ measured in $1.0 < |y| < 1.5$ for jets clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$.

The inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ measured in $1.5 < |y| < 2.0$ for jets clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $|y| < 0.5$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $0.5 < |y| < 1.0$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.0 < |y| < 1.5$and $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $|y| < 0.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $0.5 < |y| < 1.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$and $1.0 < |y| < 1.5$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The correlation matrix at the particle level of inclusive jet production cross section as a function of the jet transverse momentum~$p_\mathrm{T}$ for jets in $1.5 < |y| < 2.0$ clustered using the anti-$k_\mathrm{t}$ algorithm with $R=0.7$. It contains contributions from the data and from the \textsc{PYTHIA}~8 sample used to perform the unfolding.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The corrections were calculated by Dittmaier, Huss, Speckner.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The values correspond to the average of the corrections obtained with \textsc{{PYTHIA}}~8 and with \textsc{{HERWIG}}++.

The correlation functions (corrected for nonuniform detector efficiency in $\phi$; not corrected for the absolute detection efficiency) vs. azimuthal angle difference between forward ($2.6<\eta<4.0$) $\pi^{0}$s in $p$+$p$ collisions at $\sqrt{s_{\mathrm{_{NN}}}}=200$ GeV at high $p_{T}$ ($p^{trig}_{T}$=2.5-3 GeV/c, $p^{asso}_{T}$=2-2.5 GeV/c)

The correlation functions (corrected for nonuniform detector efficiency in $\phi$; not corrected for the absolute detection efficiency) vs. azimuthal angle difference between forward ($2.6<\eta<4.0$) $\pi^{0}$s in $p$+Al collisions at $\sqrt{s_{\mathrm{_{NN}}}}=200$ GeV at high $p_{T}$ ($p^{trig}_{T}$=2.5-3 GeV/c, $p^{asso}_{T}$=2-2.5 GeV/c)

The correlation functions (corrected for nonuniform detector efficiency in $\phi$; not corrected for the absolute detection efficiency) vs. azimuthal angle difference between forward ($2.6<\eta<4.0$) $\pi^{0}$s in $p$+$Au collisions at $\sqrt{s_{\mathrm{_{NN}}}}=200$ GeV at high $p_{T}$ ($p^{trig}_{T}$=2.5-3 GeV/c, $p^{asso}_{T}$=2-2.5 GeV/c)

Relative area of MinBias $pAu/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative area of MinBias $pAl/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative area of $pAu/pp$ at b=0 of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from prediction

Relative width of MinBias $pAu/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative width of MinBias $pAl/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative pedestal of MinBias $pAu/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative pedestal of MinBias $pAl/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative area of MinBias $pp/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$), the ratio is unity, error is zero

Relative area of MinBias $pA/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

SR1 diboson invariant mass distribution in the resolved tagged category after fitting the signal and sideband regions using a signal (ALP linear ZZ) plus background model. The last bin includes events with diboson invariant masses up to 3000 GeV. The signal is normalized to the 95% CL cross section limit at $f_a$ = 3 TeV (the scale factor used in the original figure for better visibility is not applied here).

SR2 diboson invariant mass distribution in the boosted untagged category after fitting the signal and sideband regions using a signal (ALP chiral ZH) plus background model. The last bin includes events with diboson invariant masses up to 3000 GeV. The signal is normalized to the 95% CL cross section limit at $f_a$ = 3 TeV (the scale factor used in the original figure for better visibility is not applied here).

SR2 diboson invariant mass distribution in the boosted tagged category after fitting the signal and sideband regions using a signal (ALP chiral ZH) plus background model. The last bin includes events with diboson invariant masses up to 3000 GeV. The signal is normalized to the 95% CL cross section limit at $f_a$ = 3 TeV.

SR2 diboson invariant mass distribution in the resolved untagged category after fitting the signal and sideband regions using a signal (ALP chiral ZH) plus background model. The last bin includes events with diboson invariant masses up to 3000 GeV. The signal is normalized to the 95% CL cross section limit at $f_a$ = 3 TeV (the scale factor used in the original figure for better visibility is not applied here).

SR2 diboson invariant mass distribution in the resolved tagged category after fitting the signal and sideband regions using a signal (ALP chiral ZH) plus background model. The last bin includes events with diboson invariant masses up to 3000 GeV. The signal is normalized to the 95% CL cross section limit at $f_a$ = 3 TeV (the scale factor used in the original figure for better visibility is not applied here).

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(G -> ZZ) as a function of the resonance mass, taking into account all statistical and systematic uncertainties. The electron and muon channels and the various categories used in the analysis are combined. The inner green band and the outer yellow band represent the 68 and 95% coverage of the expected limits in the background-only hypothesis.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(W' -> ZW) as a function of the resonance mass, taking into account all statistical and systematic uncertainties. The electron and muon channels and the various categories used in the analysis are combined. The inner green band and the outer yellow band represent the 68 and 95% coverage of the expected limits in the background-only hypothesis.

Observed and expected 95% CL upper limit on the product of the ALP coupling to gluons and the coupling to Z bosons as a function of the mass scale $f_a$ for ALP masses smaller than 100 GeV.

Observed and expected 95% CL upper limit on the product of the ALP coupling to gluons and the chiral $a_\widetilde{2D}$ coupling as a function of the mass scale $f_a$ for ALP masses smaller than 100 GeV.

Relative uncertainties on the measured fractional differential $p_{\mathrm{T}}^{\gamma}$ cross section.

Correlation matrix for the measured absolute differential $p_{\mathrm{T}}^{\gamma}$ cross section.

Correlation matrix for the measured fractional differential $p_{\mathrm{T}}^{\gamma}$ cross section.

Response matrix from the fiducial to the reconstructed bins of the differential $p_{\mathrm{T}}^{\gamma}$ cross section measurement.

Measured absolute differential $\eta^{\gamma}$ cross section, compared to the MG5_aMC+PY8, GENEVA, MATRIX and MCFM predictions. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$.

Measured fractional differential $\eta^{\gamma}$ cross section, compared to the MG5_aMC+PY8, GENEVA, MATRIX and MCFM predictions. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$.

Relative uncertainties on the measured absolute differential $\eta^{\gamma}$ cross section.

Relative uncertainties on the measured fractional differential $\eta^{\gamma}$ cross section.

Correlation matrix for the measured absolute differential $\eta^{\gamma}$ cross section.

Correlation matrix for the measured fractional differential $\eta^{\gamma}$ cross section.

Response matrix from the fiducial to the reconstructed bins of the differential $\eta^{\gamma}$ cross section measurement.

Measured absolute differential $\Delta R(\ell,\gamma)$ cross section, compared to the MG5_aMC+PY8, GENEVA, MATRIX and MCFM predictions. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$.

Measured fractional differential $\Delta R(\ell,\gamma)$ cross section, compared to the MG5_aMC+PY8, GENEVA, MATRIX and MCFM predictions. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$.

Relative uncertainties on the measured absolute differential $\Delta R(\ell,\gamma)$ cross section.

Relative uncertainties on the measured fractional differential $\Delta R(\ell,\gamma)$ cross section.

Correlation matrix for the measured absolute differential $\Delta R(\ell,\gamma)$ cross section.

Correlation matrix for the measured fractional differential $\Delta R(\ell,\gamma)$ cross section.

Response matrix from the fiducial to the reconstructed bins of the differential $\Delta R(\ell,\gamma)$ cross section measurement.

Measured absolute differential $\Delta\eta(\ell,\gamma)$ cross section, compared to the MG5_aMC+PY8, GENEVA, MATRIX and MCFM predictions. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$.

Measured fractional differential $\Delta\eta(\ell,\gamma)$ cross section, compared to the MG5_aMC+PY8, GENEVA, MATRIX and MCFM predictions. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$.

Relative uncertainties on the measured absolute differential $\Delta \eta(\ell,\gamma)$ cross section.

Relative uncertainties on the measured fractional differential $\Delta \eta(\ell,\gamma)$ cross section.

Correlation matrix for the measured absolute differential $\Delta \eta(\ell,\gamma)$ cross section.

Correlation matrix for the measured fractional differential $\Delta \eta(\ell,\gamma)$ cross section.

Response matrix from the fiducial to the reconstructed bins of the differential $\Delta \eta(\ell,\gamma)$ cross section measurement.

Measured absolute differential $m_{\mathrm{T}}^{\mathrm{cluster}}$ cross section, compared to the MG5_aMC+PY8, GENEVA, MATRIX and MCFM predictions. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$.

Measured fractional differential $m_{\mathrm{T}}^{\mathrm{cluster}}$ cross section, compared to the MG5_aMC+PY8, GENEVA, MATRIX and MCFM predictions. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$.

Relative uncertainties on the measured absolute differential $m_{\mathrm{T}}^{\mathrm{cluster}}$ cross section.

Relative uncertainties on the measured fractional differential $m_{\mathrm{T}}^{\mathrm{cluster}}$ cross section.

Correlation matrix for the measured absolute differential $m_{\mathrm{T}}^{\mathrm{cluster}}$ cross section.

Correlation matrix for the measured fractional differential $m_{\mathrm{T}}^{\mathrm{cluster}}$ cross section.

Response matrix from the fiducial to the reconstructed bins of the differential $m_{\mathrm{T}}^{\mathrm{cluster}}$ cross section measurement.

Measured absolute differential $N_{\mathrm{jets}}$ cross section, compared to the MG5_aMC+PY8, GENEVA, MATRIX and MCFM predictions. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$.

Correlation matrix for the measured absolute differential $N_{\mathrm{jets}}$ cross section. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$.

Response matrix from the fiducial to the reconstructed bins of the differential $N_{\mathrm{jets}}$ cross section measurement.

Relative uncertainties on the measured absolute differential $N_{\mathrm{jets}}$ cross section.

Measured absolute differential $\Delta\eta(\ell,\gamma)$ cross section, compared to the MG5_aMC+PY8, GENEVA, MATRIX and MCFM predictions. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$. Requirements of $m_{\mathrm{T}}^{\mathrm{cluster}}>150\,\mathrm{GeV}$ and a veto on the presence of any jet with $p_{\mathrm{T}} > 30\,\mathrm{GeV}$ and $|\eta| < 2.5$ are applied in addition to the baseline selection.

Measured fractional differential $\Delta\eta(\ell,\gamma)$ cross section, compared to the MG5_aMC+PY8, GENEVA, MATRIX and MCFM predictions. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 30\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 30\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$. Requirements of $m_{\mathrm{T}}^{\mathrm{cluster}}>150\,\mathrm{GeV}$ and a veto on the presence of any jet with $p_{\mathrm{T}} > 30\,\mathrm{GeV}$ and $|\eta| < 2.5$ are applied in addition to the baseline selection.

Relative uncertainties on the measured absolute differential $\Delta \eta(\ell,\gamma)$ cross section. Requirements of $m_{\mathrm{T}}^{\mathrm{cluster}}>150\,\mathrm{GeV}$ and a veto on the presence of any jet with $p_{\mathrm{T}} > 30\,\mathrm{GeV}$ and $|\eta| < 2.5$ are applied in addition to the baseline selection.

Relative uncertainties on the measured fractional differential $\Delta \eta(\ell,\gamma)$ cross section. Requirements of $m_{\mathrm{T}}^{\mathrm{cluster}}>150\,\mathrm{GeV}$ and a veto on the presence of any jet with $p_{\mathrm{T}} > 30\,\mathrm{GeV}$ and $|\eta| < 2.5$ are applied in addition to the baseline selection.

Correlation matrix for the measured absolute differential $\Delta \eta(\ell,\gamma)$ cross section. Requirements of $m_{\mathrm{T}}^{\mathrm{cluster}}>150\,\mathrm{GeV}$ and a veto on the presence of any jet with $p_{\mathrm{T}} > 30\,\mathrm{GeV}$ and $|\eta| < 2.5$ are applied in addition to the baseline selection.

Correlation matrix for the measured fractional differential $\Delta \eta(\ell,\gamma)$ cross section. Requirements of $m_{\mathrm{T}}^{\mathrm{cluster}}>150\,\mathrm{GeV}$ and a veto on the presence of any jet with $p_{\mathrm{T}} > 30\,\mathrm{GeV}$ and $|\eta| < 2.5$ are applied in addition to the baseline selection.

Response matrix from the fiducial to the reconstructed bins of the differential $\Delta \eta(\ell,\gamma)$ cross section measurement. Requirements of $m_{\mathrm{T}}^{\mathrm{cluster}}>150\,\mathrm{GeV}$ and a veto on the presence of any jet with $p_{\mathrm{T}} > 30\,\mathrm{GeV}$ and $|\eta| < 2.5$ are applied in addition to the baseline selection.

Fiducial cross section scaling terms as a function of $C_{3W}$ in all $p_{\mathrm{T}}^{\gamma} \times |\phi_{f}|$ bins. Values are given relative to the SM prediction.

Scans of the profile likelihood test statistic $q$ as a function of $C_{3W}$, given with and without the pure BSM term.

Scans of the profile likelihood test statistic $q$ as a function of $C_{3W}$, given with and without the pure BSM term.

Best-fit values of $C_{3W}$ and corresponding 95% CL confidence intervals as a function of the maximum $p_{\mathrm{T}}^{\gamma}$ bin included in the fit, with and without the pure BSM term.

Best-fit values of $C_{3W}$ and corresponding 95% CL confidence intervals as a function of the maximum $p_{\mathrm{T}}^{\gamma}$ bin included in the fit, with and without the binning in $|\phi_{f}|$.

Response matrix from the fiducial to the reconstructed bins of the differential $p_{\mathrm{T}}^{\gamma} \times |\phi_{f}|$ cross section measurement.

Measured double-differential $p_{\mathrm{T}}^{\gamma} \times |\phi_{f}|$ cross section, compared to the MG5_aMC+PY8 NLO prediction. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 80\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 150\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$. A veto on the presence of any jet with $p_{\mathrm{T}} > 30\,\mathrm{GeV}$ and $|\eta| < 2.5$ is also applied for this result.

Measured double-differential $p_{\mathrm{T}}^{\gamma} \times |\phi_{f}|$ cross section, compared to the MG5_aMC+PY8 NLO prediction. This table gives the entries for the $(0,\pi/6)$ slice. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 80\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 150\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$. A veto on the presence of any jet with $p_{\mathrm{T}} > 30\,\mathrm{GeV}$ and $|\eta| < 2.5$ is also applied for this result.

Measured double-differential $p_{\mathrm{T}}^{\gamma} \times |\phi_{f}|$ cross section, compared to the MG5_aMC+PY8 NLO prediction. This table gives the entries for the $(\pi/6,\pi/3)$ slice. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 80\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 150\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$. A veto on the presence of any jet with $p_{\mathrm{T}} > 30\,\mathrm{GeV}$ and $|\eta| < 2.5$ is also applied for this result.

Measured double-differential $p_{\mathrm{T}}^{\gamma} \times |\phi_{f}|$ cross section, compared to the MG5_aMC+PY8 NLO prediction. This table gives the entries for the $(\pi/3,\pi/2)$ slice. The differential cross sections $\sigma_{j}(\mathrm{pp}\rightarrow\mathrm{W}^{\pm}\gamma\rightarrow\ell^{\pm}\nu\gamma)$, where $\ell$ denotes all three lepton flavors, are measured in the following fiducial region: $p_{\mathrm{T}}^{\ell} > 80\,\mathrm{GeV}$, $|\eta^{\ell}| < 2.5$, $p_{\mathrm{T}}^{\gamma} > 150\,\mathrm{GeV}$, $|\eta^{\gamma}| < 2.5$, $p_{\mathrm{T}}^{\mathrm{miss}} > 40\,\mathrm{GeV}$, and $\Delta R(\ell, \gamma) > 0.7$. The leptons are dressed by adding the four-momenta of any photons with $\Delta R(\ell, \gamma) < 0.1$ to the four-momentum of the lepton. A smooth-cone photon isolation is also applied, with parameters $\delta_{0}=0.4$, $\epsilon=1.0$, and $n=1$. A veto on the presence of any jet with $p_{\mathrm{T}} > 30\,\mathrm{GeV}$ and $|\eta| < 2.5$ is also applied for this result.

Correlation matrix for the measured absolute differential $p_{\mathrm{T}}^{\gamma} \times |\phi_{f}|$ cross section.