<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $p_{T,jet}^{lead}$ region of $2.5 \;\mathrm{GeV} < p_{T,jet}^{lead} < 7 \;\mathrm{GeV}$ for $N_{jet} \geq 3$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $p_{T,jet}^{lead}$ region of $7 \;\mathrm{GeV} < p_{T,jet}^{lead} < 12 \;\mathrm{GeV}$ for $N_{jet} \geq 1$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $p_{T,jet}^{lead}$ region of $7 \;\mathrm{GeV} < p_{T,jet}^{lead} < 12 \;\mathrm{GeV}$ for $N_{jet} \geq 2$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $p_{T,jet}^{lead}$ region of $7 \;\mathrm{GeV} < p_{T,jet}^{lead} < 12 \;\mathrm{GeV}$ for $N_{jet} \geq 3$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $p_{T,jet}^{lead}$ region of $12 \;\mathrm{GeV} < p_{T,jet}^{lead} < 30 \;\mathrm{GeV}$ for $N_{jet} \geq 1$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $p_{T,jet}^{lead}$ region of $12 \;\mathrm{GeV} < p_{T,jet}^{lead} < 30 \;\mathrm{GeV}$ for $N_{jet} \geq 2$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $p_{T,jet}^{lead}$ region of $12 \;\mathrm{GeV} < p_{T,jet}^{lead} < 30 \;\mathrm{GeV}$ for $N_{jet} \geq 3$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $Q^{2}$ region of $10 \;\mathrm{GeV}^{2} < Q^{2} < 50 \;\mathrm{GeV}^{2}$ for $N_{jet} \geq 1$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $Q^{2}$ region of $10 \;\mathrm{GeV}^{2} < Q^{2} < 50 \;\mathrm{GeV}^{2}$ for $N_{jet} \geq 2$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $Q^{2}$ region of $10 \;\mathrm{GeV}^{2} < Q^{2} < 50 \;\mathrm{GeV}^{2}$ for $N_{jet} \geq 3$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $Q^{2}$ region of $50 \;\mathrm{GeV}^{2} < Q^{2} < 100 \;\mathrm{GeV}^{2}$ for $N_{jet} \geq 1$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $Q^{2}$ region of $50 \;\mathrm{GeV}^{2} < Q^{2} < 100 \;\mathrm{GeV}^{2}$ for $N_{jet} \geq 2$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $Q^{2}$ region of $50 \;\mathrm{GeV}^{2} < Q^{2} < 100 \;\mathrm{GeV}^{2}$ for $N_{jet} \geq 3$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $Q^{2}$ region of $100 \;\mathrm{GeV}^{2} < Q^{2} < 350 \;\mathrm{GeV}^{2}$ for $N_{jet} \geq 1$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $Q^{2}$ region of $100 \;\mathrm{GeV}^{2} < Q^{2} < 350 \;\mathrm{GeV}^{2}$ for $N_{jet} \geq 2$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

<b>Note: in the paper, uncertainties are given in relative terms. The HEPData table contains absolute numbers. The original data file, containing relative uncertainties as in the paper, is available via the 'Resources' button above.</b> Differential cross sections, $d\sigma/d\Delta\phi$, in the $Q^{2}$ region of $100 \;\mathrm{GeV}^{2} < Q^{2} < 350 \;\mathrm{GeV}^{2}$ for $N_{jet} \geq 3$, as obtained from the data, ARIADNE MC simulations, and perturbative calculations at $\mathcal{O}(\alpha_{s})$ and $\mathcal{O}(\alpha_{s}^{2})$ accuracy. Other details are as in the caption to Table 1.

QED correction factors for inclusive measurement, as estimated with RAPGAP.

QED correction factors for $2.5 \;\mathrm{GeV} < p_{T,jet}^{lead} < 7 \;\mathrm{GeV}$ and $N_{jet} \geq 1$, as estimated with RAPGAP.

QED correction factors for $2.5 \;\mathrm{GeV} < p_{T,jet}^{lead} < 7 \;\mathrm{GeV}$ and $N_{jet} \geq 2$, as estimated with RAPGAP.

QED correction factors for $2.5 \;\mathrm{GeV} < p_{T,jet}^{lead} < 7 \;\mathrm{GeV}$ and $N_{jet} \geq 3$, as estimated with RAPGAP.

QED correction factors for $7 \;\mathrm{GeV} < p_{T,jet}^{lead} < 12 \;\mathrm{GeV}$ and $N_{jet} \geq 1$, as estimated with RAPGAP.

QED correction factors for $7 \;\mathrm{GeV} < p_{T,jet}^{lead} < 12 \;\mathrm{GeV}$ and $N_{jet} \geq 2$, as estimated with RAPGAP.

QED correction factors for $7 \;\mathrm{GeV} < p_{T,jet}^{lead} < 12 \;\mathrm{GeV}$ and $N_{jet} \geq 3$, as estimated with RAPGAP.

QED correction factors for $12 \;\mathrm{GeV} < p_{T,jet}^{lead} < 30 \;\mathrm{GeV}$ and $N_{jet} \geq 1$, as estimated with RAPGAP.

QED correction factors for $12 \;\mathrm{GeV} < p_{T,jet}^{lead} < 30 \;\mathrm{GeV}$ and $N_{jet} \geq 2$, as estimated with RAPGAP.

QED correction factors for $12 \;\mathrm{GeV} < p_{T,jet}^{lead} < 30 \;\mathrm{GeV}$ and $N_{jet} \geq 3$, as estimated with RAPGAP.

QED correction factors for $10 \;\mathrm{GeV}^{2} < Q^{2} < 50 \;\mathrm{GeV}^{2}$ and $N_{jet} \geq 1$, as estimated with RAPGAP.

QED correction factors for $10 \;\mathrm{GeV}^{2} < Q^{2} < 50 \;\mathrm{GeV}^{2}$ and $N_{jet} \geq 2$, as estimated with RAPGAP.

QED correction factors for $10 \;\mathrm{GeV}^{2} < Q^{2} < 50 \;\mathrm{GeV}^{2}$ and $N_{jet} \geq 3$, as estimated with RAPGAP.

QED correction factors for $50 \;\mathrm{GeV}^{2} < Q^{2} < 100 \;\mathrm{GeV}^{2}$ and $N_{jet} \geq 1$, as estimated with RAPGAP.

QED correction factors for $50 \;\mathrm{GeV}^{2} < Q^{2} < 100 \;\mathrm{GeV}^{2}$ and $N_{jet} \geq 2$, as estimated with RAPGAP.

QED correction factors for $50 \;\mathrm{GeV}^{2} < Q^{2} < 100 \;\mathrm{GeV}^{2}$ and $N_{jet} \geq 3$, as estimated with RAPGAP.

QED correction factors for $100 \;\mathrm{GeV}^{2} < Q^{2} < 350 \;\mathrm{GeV}^{2}$ and $N_{jet} \geq 1$, as estimated with RAPGAP.

QED correction factors for $100 \;\mathrm{GeV}^{2} < Q^{2} < 350 \;\mathrm{GeV}^{2}$ and $N_{jet} \geq 2$, as estimated with RAPGAP.

QED correction factors for $100 \;\mathrm{GeV}^{2} < Q^{2} < 350 \;\mathrm{GeV}^{2}$ and $N_{jet} \geq 3$, as estimated with RAPGAP.

The correlation matrix for the inclusive measurement. The azimuthal correlation angle $\Delta\phi$ of each lepton--leading-jet pair was assigned a bin ID between 0 to 59 segmented uniformly from 0 to $\pi$. This bin ID was offset by multiples of 60 based on the jet multiplicity so that the pair from a single jet event, as suggested by the MC simulation, was assigned a value between 0 to 59, dijet events were given 60 to 119, trijet 120 to 179, and four or more jets 180 to 239.

The correlation matrix for $2.5 \;\mathrm{GeV} < p_{T,jet}^{lead} < 7 \;\mathrm{GeV}$. Other details are as in the caption to Fig. 13.

The correlation matrix for $7 \;\mathrm{GeV} < p_{T,jet}^{lead} < 12 \;\mathrm{GeV}$. Other details are as in the caption to Fig. 13.

The correlation matrix for $12 \;\mathrm{GeV} < p_{T,jet}^{lead} < 30 \;\mathrm{GeV}$. Other details are as in the caption to Fig. 13.

The correlation matrix for $10 \;\mathrm{GeV}^{2} < Q^{2} < 50 \;\mathrm{GeV}^{2}$. Other details are as in the caption to Fig. 13.

The correlation matrix for $50 \;\mathrm{GeV}^{2} < Q^{2} < 100 \;\mathrm{GeV}^{2}$. Other details are as in the caption to Fig. 13.

The correlation matrix for $100 \;\mathrm{GeV}^{2} < Q^{2} < 350 \;\mathrm{GeV}^{2}$. Other details are as in the caption to Fig. 13.

Expected and observed exclusion limits at 95% CL in the asymptotic approximation, on the product of the production cross section and branching fraction into two photons for an additional SM-like Higgs boson, for the VBF plus VH processes, from the analysis of the combined data from 2016, 2017, and 2018.

Expected and observed exclusion limits at 95% CL in the asymptotic approximation, on the product of the production cross section and branching fraction into two photons for an additional SM-like Higgs boson, assuming 100% production via the VBF processes, from the analysis of the combined data from 2016, 2017, and 2018.

Expected and observed exclusion limits at 95% CL in the asymptotic approximation, on the product of the production cross section and branching fraction into two photons for an additional SM-like Higgs boson, assuming 100% production via the VH processes, from the analysis of the combined data from 2016, 2017, and 2018.

Histogram of the Hbb AK8 jet softdrop mass for data, simulated signal and background events, and the background yields estimated from data in region A.

Histogram of the Hbb AK8 jet softdrop mass for data and simulated signal and background events.

Histogram of $p_{T}(H) + p_{T}(W)$ at the LHE level.

The signal yield estimated by simulation, background yield predicted from data, and observed yield in the signal region.

Observed and expected BDT output in the signal region used for the extraction of the SM-like signal strength. The yields and uncertainties of the expected contributions are shown after the fit to the data. The rightmost bin includes the overflow.

Limits on the signal strength of SM VBS WH production

Observed and expected significance of the measurement of SM VBS WH production

95% confidence limits on the signal strength in the SM analysis

A scan of $-2\Delta\ln L$, where $L$ is the likelihood, for the SM analysis plotted as a function of the signal strength $\mu$ for the expected (red) and the observed (black) data.

The data, fitted signal, and background distributions for the VBS dijet mass in the $t\overline{t}$ control region of the SM analysis.

The data, fitted signal, and background distributions for the W+jets BDT score in the W+jets control region of the SM analysis.

Fitted SME coefficients and their uncertainties at 1$\sigma$ and 2$sigma$, measured in fits of single coefficients while the coefficients corresponding to the three other directions are left floating, within the $c_L$, $c_R$, $c$, and $d$ families. The error bar includes statistical and systematic uncertainties.

Uncertainty breakdown for SME fits of single coefficients while the coefficients corresponding to the three other directions are left floating, by splitting according to the treatment of time dependence:\ flat across sidereal time (flat luminosity component, background normalization, theory), correlated in sidereal time bins (trigger, luminosity stability and linearity, pileup, MC statistical uncertainty, single top quark decay in the SME), systematics uncorrelated in sidereal time bins (other experimental uncertainties), and statistical uncertainty.

Figure 3 top left. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the RS graviton mass $m_{G}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 3 top right. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the heavy higgs mass $m_{S}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 3 mid left. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the RS graviton mass $m_{G}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 3 mid right. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the heavy higgs mass $m_{S}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 3 bottom left. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the RS graviton mass $m_{G}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 3 bottom right. Expected and observed 95% CL upper limits on the product of the production cross section and branching fraction as a function of the heavy higgs mass $m_{S}$ for the full Run 2 data set are shown. Expected $1\sigma$ and $2\sigma$ limit bands are shown in green and yellow, respectively

Figure 4 top left. Expected 95% CL upper limits on the product of the cross section and branching fraction as a function of the $m_{G}$ versus the resonance width for the full Run 2 data set.

Figure 4 top right. Observed 95% CL upper limits on the product of the cross section and branching fraction as a function of the $m_{G}$ versus the resonance width for the full Run 2 data set.

Figure 4 bottom left. Expected 95% CL upper limits on the product of the cross section and branching fraction as a function of the $m_{S}$ versus the resonance width for the full Run 2 data set.

Figure 4 bottom right. Observed 95% CL upper limits on the product of the cross section and branching fraction as a function of the $m_{S}$ versus the resonance width for the full Run 2 data set.

Figure 6. The $m_{\gamma\gamma}$ spectra and background prediction after nuisance paramter marginalization (post-fit) due to SM diphoton production ($\gamma\gamma$) and fake photon production (j$\gamma$, jj) for the EBEB, combining the 2016, 2017, and 2018 data sets. The prediction with a large extra dimensions signal (GRW convention with $M_{S}$ = 9 TeV) is also shown. The pull distributions, normalized to the statistical uncertainty only, are shown underneath the spectra. The last bin contains the overflow of events beyond $m_{\gamma\gamma}$ > 3.5 TeV.

Figure 6. The $m_{\gamma\gamma}$ spectra and background prediction after nuisance paramter marginalization (post-fit) due to SM diphoton production ($\gamma\gamma$) and fake photon production (j$\gamma$, jj) for the EBEE, combining the 2016, 2017, and 2018 data sets. The prediction with a large extra dimensions signal (GRW convention with $M_{S}$ = 9 TeV) is also shown. The pull distributions, normalized to the statistical uncertainty only, are shown underneath the spectra. The last bin contains the overflow of events beyond $m_{\gamma\gamma}$ > 3.5 TeV.

Figure 7. The exclusion limit for the clockwork framework over the k--M5 parameter space. The shaded region denotes where the theory becomes nonperturbative. The region below and to the left of the solid line constitutes the excluded region. Expected 1std and 2std limit bands are shown in green and yellow, respectively.

Weights from b tagging efficiency ratios as functions of the five-jet invariant mass in 2018 data for the high-mass selection, connecting the 2M1L and 3M regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the high-mass analysis only signals with mass above 800GeV are tested, so primarily the lower part of the distribution contributes to the final result.

Weights from b tagging efficiency ratios as functions of the five-jet invariant mass in 2018 data for the high-mass selection, connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the high-mass analysis only signals with mass above 800GeV are tested, so primarily the lower part of the distribution contributes to the final result.

The five-jet invariant mass distribution in the tH channel (black markers) after the high-mass selection in the QCD multijet 3T control region for the 2018 data set. The histograms are the corresponding reweighted 2M1L distributions. The background distribution is normalized to the number of entries in the data. The shaded area corresponds to the statistical uncertainties in the 2M1L control regions. A potential 900 GeV T' signal (red cross-hatched histogram) is added to the background histogram demonstrating a negligible contribution. Similar results are observed in the tZ channel, and for the other years, but with slightly larger statistical uncertainties.

The five-jet invariant mass distribution in the tH channel (black markers) after the high-mass selection in the ttbar 2T1L control region for the 2018 data set. The histograms are the corresponding reweighted 2M1L distributions. The background distribution is normalized to the number of entries in the data. The shaded area corresponds to the statistical uncertainties in the 2M1L control regions. A potential 900 GeV T' signal (red cross-hatched histogram) is added to the background histogram demonstrating a negligible contribution. Similar results are observed in the tZ channel, and for the other years, but with slightly larger statistical uncertainties.

The five-jet invariant mass distribution in the tH channel (black markers) after the high-mass selection in the 3M signal region for the 2018 data set. The histograms are the corresponding reweighted 2M1L distributions. The background distribution is normalized to the number of entries in the data. The shaded area corresponds to the statistical uncertainties in the 2M1L control regions. A potential 900 GeV T' signal (red cross-hatched histogram) is added to the background histogram demonstrating a negligible contribution. Similar results are observed in the tZ channel, and for the other years, but with slightly larger statistical uncertainties.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the high-mass selection. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the high-mass selection. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the high-mass selection. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the all-tH channel.

Observed p-values when considering the tH channel for each year and their combination.

Observed and expected 95% CL upper limits on the product of the production cross section for a single VLQ T production and the branching ration given as functions of $m_{T'}$ for a relative width of 1%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a single VLQ T production and the branching ration given as functions of $m_{T'}$ for a relative width of 1%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a single VLQ T production and the branching ration given as functions of $m_{T'}$ for a relative width of 1%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a single VLQ T production and the branching ration given as functions of $m_{T'}$ for a relative width of 1%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Weights from b tagging efficiency ratios as functions of the five-jet invariant mass in 2016 data for the low-mass selection, connecting the 2M1L and 3M regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the low-mass analysis only signals with mass below 800GeV are tested, so primarily the lower part of the distribution contributes to the final result.

Weights from b tagging efficiency ratios as functions of the five-jet invariant mass in 2017 data for the low-mass selection, connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the low-mass analysis only signals with mass below 800GeV are tested, so primarily the lower part of the distribution contributes to the final result.

Weights from b tagging efficiency ratios as functions of the five-jet invariant mass in 2016 data for the low-mass selection, connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the low-mass analysis only signals with mass below 800GeV are tested, so primarily the lower part of the distribution contributes to the final result.

Weights from b tagging efficiency ratios as functions of the five-jet invariant mass in 2017 data for the low-mass selection, connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the low-mass analysis only signals with mass below 800GeV are tested, so primarily the lower part of the distribution contributes to the final result.

Weights from b tagging efficiency ratios as functions of the five-jet invariant mass in 2016 data for the high-mass selection, connecting the 2M1L and 3M regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the high-mass analysis only signals with mass above 800GeV are tested, so primarily the lower part of the distribution contributes to the final result.

Weights from b tagging efficiency ratios as functions of the five-jet invariant mass in 2017 data for the high-mass selection, connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the high-mass analysis only signals with mass above 800GeV are tested, so primarily the lower part of the distribution contributes to the final result.

Weights from b tagging efficiency ratios as functions of the five-jet invariant mass in 2016 data for the high-mass selection, connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the high-mass analysis only signals with mass above 800GeV are tested, so primarily the lower part of the distribution contributes to the final result.

Weights from b tagging efficiency ratios as functions of the five-jet invariant mass in 2017 data for the high-mass selection, connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the high-mass analysis only signals with mass above 800GeV are tested, so primarily the lower part of the distribution contributes to the final result.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 2M1L regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3M regions for the high-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the low-mass selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Background: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the high-mass selection. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Data: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the high-mass selection. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.

Signal: Five-jet invariant mass distributions after a background-only fit (blue histogram) to the complete dataset (black markers) in the 3T regions for the high-mass selection. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 900 GeV T'. The fit is performed on the combined data from all 3 years in the tZ and tH channel.