Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Signal yields for two charge scenarios considered in the analysis, as well as their associated uncertainties.

Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Signal yields for two charge scenarios considered in the analysis, as well as their associated uncertainties.

Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the early 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 6%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text. Comparing with respect to the binomial fit starting at $N_{\text{hits}}^{\text{low dE/dx}} = 2$, and not $N_{\text{hits}}^{\text{low dE/dx}} = 1$, is needed to account for the fact that early 2016 data is more strongly affected low dE/dx by instrumental effects that widen the N hits distribution.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the early 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 6%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text. Comparing with respect to the binomial fit starting at $N_{\text{hits}}^{\text{low dE/dx}} = 2$, and not $N_{\text{hits}}^{\text{low dE/dx}} = 1$, is needed to account for the fact that early 2016 data is more strongly affected low dE/dx by instrumental effects that widen the N hits distribution.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the late 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 78%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the late 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 78%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2017 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 65%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2017 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 65%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2018 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 9%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2018 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 9%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Exclusion region (hatched) at 95% CL in the FCP charge-mass plane for the considered signal. The expected exclusion is shown with the associated 1 and 2 standard deviations $\sigma$ bands. Signal points at charges 0.9, 0.8, 2/3, 0.5, and 1/3 e are connected by straight lines to guide the eye. This is a conservative interpolation. Previous exclusions from CMS [Phys. Rev. D 87 (2013) 092008), JHEP 07 (2013) 122] as well as OPAL [Phys. Lett. B 572 (2003) 8] are given for comparison.

Exclusion region (hatched) at 95% CL in the FCP charge-mass plane for the considered signal. The expected exclusion is shown with the associated 1 and 2 standard deviations $\sigma$ bands. Signal points at charges 0.9, 0.8, 2/3, 0.5, and 1/3 e are connected by straight lines to guide the eye. This is a conservative interpolation. Previous exclusions from CMS [Phys. Rev. D 87 (2013) 092008), JHEP 07 (2013) 122] as well as OPAL [Phys. Lett. B 572 (2003) 8] are given for comparison.

Particle-level distributions of ungroomed AK4 pTD2 in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 thrust in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 thrust in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 width in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 width in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 LHA in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 LHA in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 LHA in 408 < PT < 1500 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 LHA in 1000 < PT < 4000 GeV in the central dijet region.

Particle-level distributions of ungroomed AK8 LHA in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK8 LHA in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 LHA (charged-only) in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 LHA (charged-only) in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of groomed AK4 LHA in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of groomed AK4 LHA in 120 < PT < 150 GeV in the central dijet region.

Mean of ungroomed LHA for AK4 jets as a function of PT in the Z+jet region.

Mean of ungroomed LHA for AK4 jets as a function of PT in the central dijet region.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of dimuon invariant mass. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Forward-backward ratios for $15<m_{\mu\mu}<60$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Forward-backward ratios for $60<m_{\mu\mu}<120$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of dimuon invariant mass.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV.

Reconstruction level differential cross section spectla, obtained for the central acceptance, $|\eta| < 3$, for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology compared to the to the EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The number of high purity tracksin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The transeverse momentum of high purity tracksin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The total energy of all Particle Flow candidatesin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

Unfolded diffraction enhanced differential cross section for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffraction enhanced differential cross section for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, compared to the to the hadron-level EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, compared to the to the hadron-level EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, compared to the to the hadron-level QGSJET II predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, compared to the to the hadron-level QGSJET II predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Reconstruction level diffraction enhanced differential cross section spectrum for Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap, compared with the spectrum obtained with events satisfying the ZDC veto requirement $E_{ZDC−} < 1 \mathrm{~TeV}$ which selects only the events without lead nuclear break up (without correction for HF undetectable energy). Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Cross section upper limit vs m(GLUINO) for SMS model T5ZZ.

Cross section upper limit vs m(GLUINO) for SMS model T5ZZ.

Cross section upper limit vs m(GLUINO) for SMS model T5ZZ.

Cross section vs m(GLUINO) for SMS model T5ZZ.

Cross section vs m(GLUINO) for SMS model T5ZZ.

Cross section vs m(GLUINO) for SMS model T5ZZ.

Number of events in the pTmiss CR, transfer factor, background prediction, and observed yield in each of the six pTmiss bins. The systematic uncertainties in the background prediction include the shape uncertainties in addition to the uncertainty in T. Also listed in the last column is the number of expected signal events and corresponding statistical uncertainties for one example mass point, m(GLUINO) = 1700 GeV. The last pTmiss bin includes overflows.

Pre-fit background covariance matrix $\sigma_{xy}$ for the 6 pTmiss bins.

Pre-fit background correlation matrix $\rho_{xy}$ for the 6 pTmiss bins.

Observed significance (in standard deviations) of the T5ZZ signal model as a function of the gluino mass.

Cut flow for the signal for two example mass points of the T5ZZ signal model. Here the hadronic baseline is defined by $N_{\mathrm{jets}}\ge 2$, $p_T^{\mathrm{miss}}$ > 300 GeV, $H_T$ > 300 GeV, $\Delta\phi$ cuts, and the isolated photon, lepton and track veto.

The observed DNN output distribution for data collected in 2016 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution for data collected in 2017 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution for data collected in 2018 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution in the VBF-SB regions for the combination of 2016, 2017, and 2018 data, compared to the post-fit prediction from SM processes. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution in the VBF-SR region for the combination of 2016, 2017, and 2018 data, compared to the post-fit prediction from SM processes. The post-fit distributions for the Higgs boson signal produced via ggH and VBF modes with mass of 125.38 GeV, along with the total post-fit and pre-fit uncertainties on the background prediction, are also reported.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ggH category, with dimuon mass between 110–150 GeV, are considered. The expected distributions for ggH, VBF, and other signal processes are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat3 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat4 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat5 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ttH hadronic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ttH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHhad-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHhad-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHhad-cat3 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHlep-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHlep-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the WH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ZH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

Comparison between the data and the total background extracted from a S+B fit performed in the WH-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the WH-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the WH-cat3 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ZH-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ZH-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Observed local p-values as a function of the Higgs boson mass hypothesis extracted from the combined fit across event categories as well as from each individual production category.

Expected local p-values as a function of the Higgs boson mass hypothesis extracted from the combined fit across event categories as well as from each individual production category. They are calculated using the background-only expectation obtained from the S+B fit and injecting a signal with mass of 125.38 GeV and strenght of one.

Signal strength modifiers measured for a Higgs boson with mass of 125.38 GeV in each production category as well as from the combined fit.

Scan of the profiled likelihood ratio as a function of $\mu_{ggH,t\overline{t}H}$ and $\mu_{VBF,VH}$.

The 68% CL contour for the profiled likelihood scan as a function of $\mu_{ggH,t\overline{t}H}$ and $\mu_{VBF,VH}$.

The 95% CL contour for the profiled likelihood scan as a function of $\mu_{ggH,t\overline{t}H}$ and $\mu_{VBF,VH}$.

The mass distribution for the weighted combination of VBF-SB and VBF-SR events. Each event is weighted proportionally to the S/(S+B) ratio, calculated as a function of the mass-decorrelated DNN output. The total post-fit and pre-fit uncertainties on the background prediction, along with the best fit signal contribution, are also reported.

Comparison between the data and the total background extracted for the S/(S+B) weighted combination of all event categories. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Observed and expected local p-values as a function of the Higgs boson mass hypothesis extracted from the combined fit performed on data recorded at centre-of-mass energies of 7, 8, and 13 TeV. The expected p-values are calculated using the background-only expectation obtained from the S+B fit and injecting a signal with mass of 125.38 GeV and strenght of one.

The observed profile likelihood ratio as a function of $\kappa_{\mu}$ for a Higgs boson mass of 125.38 GeV, obtained from a combined fit with Ref.[10] in the resolved $\kappa$-framework. The best fit value for $\kappa_{\mu}$ is 1.07 and the corresponding observed 68% CL interval is $0.85 < \kappa_{\mu} < 1.29$.

The best fit estimates for the reduced coupling modifiers extracted for fermions and weak bosons from the resolved $\kappa$-framework compared to their corresponding prediction from the SM. The error bars represent 68% CL intervals for the measured parameters.

The distribution of Top jet $ au_{3}/ au_{2}$ for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of W jet $ au_{2}/ au_{1}$ for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of m_{W} for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of Top jet DeepCSV score for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of m_{t} for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

Distributions of $m_{tW}$ in the signal region for an interval of $m_{t}$ (65-105 GeV).

Distributions of $m_{tW}$ in the signal region for an interval of $m_{t}$ (105-225 GeV).

Distributions of $m_{tW}$ in the signal region for an interval of $m_{t}$ (225-285 GeV).

Distributions of $m_{tW}$ in themultijet control region of the signal region for an interval of $m_{t}$ (65-105 GeV).

Distributions of $m_{tW}$ in themultijet control region of the signal region for an interval of $m_{t}$ (105-225 GeV).

Distributions of $m_{tW}$ in themultijet control region of the signal region for an interval of $m_{t}$ (225-285 GeV).

Distributions of $m_{t}$ in the ttbar measurement region for an interval of $m_{tar{t}}$ (1200-1300 GeV).

Distributions of $m_{t}$ in the ttbar measurement region for an interval of $m_{tar{t}}$ (1300-1800 GeV).

Distributions of $m_{t}$ in the ttbar measurement region for an interval of $m_{tar{t}}$ (1800-3000 GeV).

Distributions of $m_{t}$ in themultijet control region of the ttbar measurement region for an interval of $m_{tar{t}}$ (1200-1300 GeV).

Distributions of $m_{t}$ in themultijet control region of the ttbar measurement region for an interval of $m_{tar{t}}$ (1300-1800 GeV).

Distributions of $m_{t}$ in themultijet control region of the ttbar measurement region for an interval of $m_{tar{t}}$ (1800-3000 GeV).

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the HP VBF category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the HP ggF/DY category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the LP VBF category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the LP ggF/DY category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

Distributions of mT for high-purity ggF/DY CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for high-purity VBF CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity ggF/DY CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity VBF CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for high-purity ggF/DY SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for high-purity VBF SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity ggF/DY SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity VBF SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Expected and observed 95% CL upper limits on the product of the radion (R) production (gluon-gluon fusion) cross section and the R -> ZZ branching fraction versus the radion mass.

Expected and observed 95% CL upper limits on the product of the radion (R) production (vector boson fusion) cross section and the R -> ZZ branching fraction versus the radion mass.

Expected and observed 95% CL upper limits on the product of the W' (W') production (Drell-Yan) cross section and the W' -> WZ branching fraction versus the W' mass.

Expected and observed 95% CL upper limits on the product of the W' (W') production (vector boson fusion) cross section and the W' -> WZ branching fraction versus the W' mass.

Expected and observed 95% CL upper limits on the product of the graviton (G) production (gluon-gluon fusion) cross section and the G -> ZZ branching fraction versus the graviton mass.

Expected and observed 95% CL upper limits on the product of the graviton (G) production (vector boson fusion) cross section and the G -> ZZ branching fraction versus the graviton mass.