Expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRGGWZ-H.

N-1 distribution for $E_{\mathrm{T}}^{\mathrm{miss}}$of observed data and expected background in SRGGSlep-M.

N-1 distribution for $\sum{p_{\mathrm{T}}^\mathrm{jet}}$of observed data and expected background in SRUDD-ge2b.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRLQD.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSWZ-H.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSSlep-H(loose).

Signal acceptance for SRGGWZ-H signal region from Fig 10(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-H signal region from Fig 15(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-M signal region from Fig 10(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-M signal region from Fig 15(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-L signal region from Fig 10(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-L signal region from Fig 15(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-L signal region from Fig 12(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-L signal region from Fig 17(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-M signal region from Fig 12(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-M signal region from Fig 17(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-H signal region from Fig 12(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-H signal region from Fig 17(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRUDD-1b signal region from Fig 14(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-1b signal region from Fig 19(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-2b signal region from Fig 14(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-2b signal region from Fig 19(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge2b signal region from Fig 14(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge2b signal region from Fig 19(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge3b signal region from Fig 14(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge3b signal region from Fig 19(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRLQD signal region from Fig 14(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal efficiency for SRLQD signal region from Fig 19(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal acceptance for SRSSWZ-L signal region from Fig 11(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-L signal region from Fig 16(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-ML signal region from Fig 11(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-ML signal region from Fig 16(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-MH signal region from Fig 11(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-MH signal region from Fig 16(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-H signal region from Fig 11(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-H signal region from Fig 16(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H signal region from Fig 13(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H signal region from Fig 18(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-MH signal region from Fig 13(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-MH signal region from Fig 18(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-L signal region from Fig 13(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-L signal region from Fig 18(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-ML signal region from Fig 13(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-ML signal region from Fig 18(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H(loose) signal region from Fig 13(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H(loose) signal region from Fig 18(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-H in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-M in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-L in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-L in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-M in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-H in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-1b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge3b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRLQD in a susy scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2200 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1870 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-L in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-ML in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-MH in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-H in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-MH in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-L in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-ML in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H(loose) in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Cross-section upper limits at 95% CL from Fig1(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Cross-section upper limits at 95% CL from Fig1(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Cross-section upper limits at 95% CL from Fig1(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Observed values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs and double-Higgs analyses combination, with all other coupling modifiers fixed to unity.

Observed values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the double-Higgs analyses, with all other coupling modifiers fixed to unity.

Observed values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs analyses, with all other coupling modifiers fixed to unity.

Observed values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs and double-Higgs combination for the generic model (free floating $\kappa_t$, $\kappa_b$, $\kappa_V$ and $\kappa_\tau$). The observed best-fit value of $\kappa_\lambda$ for the generic model is shifted slightly relative to the other models because of its correlation with the best-fit values of the $\kappa_b$, $\kappa_t$ and $\kappa_\tau$ parameters, which are slightly below, but compatible with unity.

Expected values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs and double-Higgs analyses combination derived from the combined single-Higgs and double-Higgs analyses, with all other coupling modifiers fixed to unity.

Expected values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the double-Higgs analyses.

Expected values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs analyses, with all other coupling modifiers fixed to unity.

Expected values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs and double-Higgs analyses for the generic model (free floating $\kappa_t$, $\kappa_b$, $\kappa_V$ and $\kappa_\tau$).

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs and double-Higgs combination. The solid lines show the 68% CL contours. The observed constraint for the single- and double-Higgs combination for $\kappa_t$ values below unity is slightly less stringent than that for the single-Higgs fit alone due to the slightly higher best-fit value for this coupling modifier.

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs and double-Higgs combination. The dashed lines show the 95% CL contours. The observed constraint for the single- and double-Higgs combination for $\kappa_t$ values below unity is slightly less stringent than that for the single-Higgs fit alone due to the slightly higher best-fit value for this coupling modifier.

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from double-Higgs analysis. The solid lines show the 68% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from double-Higgs analysis. The dashed lines show the 95% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs analysis. The solid lines show the 68% CL contours.

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs analysis. The dashed lines show the 95% CL contours.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs and double-Higgs combination. The solid lines show the 68% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs and double-Higgs combination. The dashed lines show the 95% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from double-Higgs analyses. The solid lines show the 68% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from double-Higgs analyses. The dashed lines show the 95% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs analyses. The solid lines show the 68% CL contours.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs analyses. The dashed lines show the 95% CL contours.

Observed and expected 95% CL upper limits on the sum of the ggF HH and VBF HH production cross-section from the bbbb, bb$\tau\tau$ and bb$\gamma\gamma$ decay channels, and their statistical combination. The value $m_H$=125.09 GeV is assumed when deriving the predicted SM cross section. The expected limit and the corresponding error bands are derived assuming the absence of the HH process with all nuisance parameters profiled to the observed data. The SM prediction together with its theoretical uncertainty is also shown (red vertical band).

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the HH to bbbb analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the HH to bb$\tau\tau$ analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the HH to bb$\gamma\gamma$ analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the double-Higgs combination. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for HH to bbbb analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for HH to bb$\tau\tau$ analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for HH to bb$\gamma\gamma$ analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for double-Higgs combination. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bbbb analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bb$\tau\tau$ analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bb$\gamma\gamma$ analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the double-Higgs combination. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bbbb analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bb$\tau\tau$ analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bb$\gamma\gamma$ analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the double-Higgs combination. All other coupling modifiers are fixed to their SM value.

Observed constraints in the $\kappa_{2V}$–$\kappa_{V}$ plane from double-Higgs combination. The solid lines show the 68% (95%) CL contours.

Observed constraints in the $\kappa_{2V}$–$\kappa_{V}$ plane from double-Higgs combination. The dashed lines show the 68% (95%) CL contours.

Expected constraints in the $\kappa_{2V}$-$\kappa_{V}$ plane from double-Higgs combination. The solid lines show the 68% CL contours.

Expected constraints in the $\kappa_{2V}$-$\kappa_{V}$ plane from double-Higgs combination. The dashed lines show the 95% CL contours.

Distribution of the diphoton invariant mass in validation region VR2. The solid histograms are stacked to show the SM expectations after the 2×2D background estimation technique is applied. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb̄h, all subdominant in this signature. Statistical and systematic uncertainties are indicated by the shaded area. The lower panel of each plot shows the ratio of the data to the SM prediction for the respective bin. The first and last bins include the underflows and overflows respectively.

Distribution of the missing transverse momentum in validation region VR2. The solid histograms are stacked to show the SM expectations after the 2×2D background estimation technique is applied. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb̄h, all subdominant in this signature. Statistical and systematic uncertainties are indicated by the shaded area. The lower panel of each plot shows the ratio of the data to the SM prediction for the respective bin. The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{γγ} itself, for signal region SR1h. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying χ̃⁰₁ mass are also overlaid (not stacked) assuming B(χ̃⁰₁ → hG̃ ) to equal 100%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb̄h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{γγ}-125 GeV|<5 GeV). The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{γγ} itself, for signal region SR1Z. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying χ̃⁰₁ mass are also overlaid (not stacked) assuming B(χ̃⁰₁ → hG̃ ) to equal 50%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb̄h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{γγ}-125 GeV|<5 GeV). The first and last bins include the underflows and overflows respectively.

Distribution of the diphoton invariant mass with all selections of the signal regions applied, except on m_{γγ} itself, for signal region SR2. The background estimation techniques described in the text are applied. The different backgrounds are stacked to add up to the total SM prediction in each bin. The predicted yields for signal benchmark models of varying χ̃⁰₁ mass are also overlaid (not stacked) assuming B(χ̃⁰₁ → hG̃ ) to equal 100%. Background and signal predictions are normalised to the luminosity. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb̄h, all subdominant in this signature. The sizes of the statistical and systematic uncertainties are indicated by the shaded areas. The lower panels show the ratio of the data to the SM prediction. Arrows indicate the borders of the signal region (|m_{γγ}-125 GeV|<5 GeV). The first and last bins include the underflows and overflows respectively.

Observed and expected limits on the pure higgsino cross-section at 95% CL assuming B(χ̃⁰₁ → hG̃ )=100% for different χ̃⁰₁ masses, obtained by a statistical combination of the three signal regions SR1h, SR1Z and SR2. The inner and outer bands indicate the 1σ and 2σ variation on the expected limit respectively.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Observed and expected 95% CL limits on the pure-higgsino branching fraction to B(χ̃⁰₁ → hG̃ ) as a function of the higgsino mass m(χ̃⁰₁) assuming it decays via either χ̃⁰₁→ hG̃ or χ̃⁰₁→ ZG̃. Limits are obtained by performing a statistical combination of the three signal regions SR1h, SR1Z and SR2. The ± 1σ variation on the expected limit is shown. The dotted lines indicate the observed limit obtained by a variation of theoretical prediction for the neutralino production cross-section by ±1 σ. Values of B(χ̃⁰₁ → hG̃ ) larger than the observed 95% CL limit are excluded, as indicated by the hatched area.

Numbers of signal and background events in the signal regions. The respective background estimation techniques are applied. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The different backgrounds are stacked to add up to the total Standard Model prediction in each bin. The predicted yields for signal benchmark models of varying &chi;&#771;⁰₁ mass are also plotted (not stacked), assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% and a &chi;&#771;⁰₁ mass of 130 or 200 GeV. The statistical and systematic uncertainties are indicated by the shaded areas in the top plot. The bottom panel shows the statistical significance&nbsp;<a href="https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PUBNOTES/ATL-PHYS-PUB-2020-025/">[Ref]</a> of the difference between the SM prediction and the observed data in each region.

Observed 95% CL limits in pb on the pure higgsino cross-section, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Limits are obtained by a statistical combination of the three signal regions SR1h, SR1Z and SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR1h, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR1Z, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Acceptances for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Acceptances are provided separately for signal region SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR1h, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR1Z, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Efficiencies for all signal model points considered in the analysis, shown in the m(&chi;&#771;⁰₁)-B(&chi;&#771;⁰₁ &rarr; hG&#771; ) plane. Efficiencies are provided separately for signal region SR2, assuming the neutralino to decay via either &chi;&#771;⁰₁&rarr; hG&#771; or &chi;&#771;⁰₁&rarr; ZG&#771;.

Observed and expected numbers of events in the three signal regions. The background category "h (other)" includes events originating from VBF, Vh, ggF, thq, thW and bb&#772;h, all subdominant in this signature. The table also includes model-independent 95% CL upper limits on the visible number of BSM events (S⁹⁵_obs), the number of BSM events given the expected number of background events (S⁹⁵_exp) and the visible BSM cross-section (&lang;&epsilon; &sigma;&rang;_obs⁹⁵), all calculated from pseudo-experiments. The discovery p-value (p₀) is also shown and its value is capped at 0.5 if the observed number of events is below the expected number of events.

Cutflows of two benchmark signal points assuming B(&chi;&#771;⁰₁ &rarr; hG&#771; )=100% for all three discovery signal regions. The initial selection includes the leptons veto. Only statistical uncertainties are included. Expected yields are normalised to a luminosity of 139 fb^-1.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-0 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-2 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-0.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-2.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-0 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-2 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Comparison of the background model (stacked histograms) with data (dots) in the $2b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $3b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $4b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-0 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-2 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The theoretical prediction for the bulk RS model with $k/\bar{M}_{\text{Pl}} = 1$ (solid red line) is shown; the decrease below 350 GeV is due to a sharp reduction in the $G^{*}_{\text{KK}} \rightarrow HH$ branching ratio. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

The impact of the eight most important systematic uncertainties on \(R_t\), in %.

The post-fit yields in the two SRs. All uncertainties applied in the analysis are included.

Limits on EFT parameters and Vtb extracted from the analysis. The upper and lower limits correspond to 95% CL.

Results of the main analysis.

Pre-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR plus.

Post-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR plus.

Post-fit distribution of the \(\Delta\phi(\vec{p}_{T}^{miss},\mu)\) variable in the muon-plus CR.

Pre-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR minus.

Post-fit distribution of the NN discriminant \(D_{\text{nn}}\) in SR minus.

Post-fit distribution of the \(\Delta\phi(\vec{p}_{T}^{miss},\mu)\) variable in the muon-minus CR.

Pre-fit distribution of the invariant mass of the untagged jet and the b-tagged jet, \(m(jb)\), in SR plus.

Pre-fit distribution of the absolute value of the pseudorapidity of the untagged jet, \(|\eta(j)|\), in SR plus.

Pre-fit distribution of the absolute value of the difference in \(p_{\text{T}}\) between the reconstructed W boson and the jet pair, \(|\Delta p_{\text{T}}(W, jb)|\), in SR plus.

Pre-fit distribution of the difference in azimuth angle between the reconstructed W boson and the jet pair, \(|\Delta\phi(W, jb)|\), in SR plus.

Pre-fit distribution of the invariant mass of the reconstructed top quark, \(m(t)\), in SR plus.

Pre-fit distribution of the absolute value of the difference in pseudorapidity between the charged lepton and the untagged jet, \(|\Delta\eta(\ell ,j)|\), SR plus.

Pre-fit distribution of the angular distance of the charged lepton and the untagged jet, \(\Delta R(\ell ,j)\), in SR plus.

Pre-fit distribution of the absolute value of the difference in pseudorapidity between the b-tagged jet and the charged lepton, \(|\Delta\eta(b,\ell )|\), in SR plus.

The symmetric cumulant $sc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The symmetric cumulant $sc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The symmetric cumulant $sc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The asymmetric cumulant $ac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The asymmetric cumulant $ac_{2}\{3\}$results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The asymmetric cumulant $ac_{2}\{3\}$results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The asymmetric cumulant $ac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The asymmetric cumulant $ac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The asymmetric cumulant $ac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 13 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized symmetric cumulant $nsc_{2\,3}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized symmetric cumulant $nsc_{2\,4}\{4\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The normalized asymmetric cumulant $nac_{2}\{3\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pp collisions at $\sqrt{s_{NN}}$ = 13 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in pPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The $v_{2}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{3}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The $v_{4}\{2\}$ results as a function of multiplicity ($N_{ch}$) in PbPb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV

The symmetric cumulant $ac_{2}\{3\}$ in Pb+Pb from different methods

The symmetric cumulant $ac_{2}\{3\}$ in Pb+Pb from different methods

The symmetric cumulant $ac_{2}\{3\}$ in p+Pb from different methods

The symmetric cumulant $ac_{2}\{3\}$ in p+Pb from different methods

The symmetric cumulant $ac_{2}\{3\}$ in pp from different methods

The symmetric cumulant $ac_{2}\{3\}$ in pp from different methods

Pre-fit background composition of the SR$3\ell$ Prod. The table shows the event yields as opposed to just the percentages of the relevant background processes.

Post-fit plot of $H_\text{T}(\text{jets})$ in the SR$2\ell$ Dec from a signal-blinded background-only fit.

Post-fit plot of $m(t_\text{SM}, b\text{-jet}_0)$ in the SR$2\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Dec from a signal-blinded background-only fit.

Post-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Prod from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Prod from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$e$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$\mu$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ $t\bar{t}V$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$e$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$\mu$ from the full fit to data.

Post-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}W$ from the full fit to data.

Post-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}Z$ from the full fit to data.

Observed and expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hu)$ for different analyses and their statistical combination.

Observed and expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hc)$ for different analyses and their statistical combination.

Post-fit normalisation factors of free-floating background processes and the signal normalisation.

Post-fit predicted and observed yields in all $2\ell$SS signal and control regions. Pre-fit signal contributions for a signal cross section equivalent to $\mathcal{B}(t\to Hq)=0.1\,\%$ are given as well.

Post-fit predicted and observed yields in all $3\ell$ signal and control regions. Pre-fit signal contributions for a signal cross section equivalent to $\mathcal{B}(t\to Hq)=0.1\,\%$ are given as well.

Expected upper limits on $\mathcal{B}(t\to Hq)$ for the nominal fit and alternative fit configurations. One contains the full phase space but only considers statistical uncertainties. Two other configurations consider the full set of systematic uncertainties, but only encompass one final state.

Expected and observed upper limits on $\mathcal{B}(t\to Hq)$ and $|C_{u\phi}^{qt,tq}|$ for the full fit containing all systematic uncertainties.

Pre-fit plot of $H_\text{T}(\text{jets})$ in the SR$2\ell$ Dec from a signal-blinded background-only fit.

Pre-fit plot of $m(t_\text{SM}, b\text{-jet}_0)$ in the SR$2\ell$ Prod from a signal-blinded background-only fit.

Pre-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Dec from a signal-blinded background-only fit.

Pre-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Prod from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$e$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$\mu$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ $t\bar{t}V$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$e$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$\mu$ from the full fit to data.

Pre-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}W$ from the full fit to data.

Pre-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}Z$ from the full fit to data.

Ranking of fit nuisance parameters according to their impact on the post-fit $tHu$ signal normalisation when fixed to $\pm1\sigma$

Ranking of fit nuisance parameters according to their impact on the post-fit $tHc$ signal normalisation when fixed to $\pm1\sigma$

Expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hu)$ for each individual final state and the full analysis.

Expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hc)$ for each individual final state and the full analysis.