<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.

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.

Isoscalar interaction coupling limit for Lagrangian 7

Isoscalar interaction coupling limit for Lagrangian 10

Isovector interaction coupling limit for Lagrangian 14

Isoscalar interaction coupling limit for Lagrangian 6

Isovector interaction coupling limit for Lagrangian 7

Isovector interaction coupling limit for Lagrangian 9

Isoscalar interaction coupling limit for Lagrangian 18

Isoscalar interaction coupling limit for Lagrangian 4

Isoscalar interaction coupling limit for Lagrangian 2

Isovector interaction coupling limit for Lagrangian 18

Isoscalar interaction coupling limit for Lagrangian 9

Isovector interaction coupling limit for Lagrangian 1

Isovector interaction coupling limit for Lagrangian 8

Isovector interaction coupling limit for Lagrangian 5

Isovector interaction coupling limit for Lagrangian 17

Isoscalar interaction coupling limit for Lagrangian 5

Isoscalar interaction coupling limit for Lagrangian 8

Isovector interaction coupling limit for Lagrangian 12

Isoscalar interaction coupling limit for Lagrangian 11

Isovector interaction coupling limit for Lagrangian 15

Isoscalar interaction coupling limit for Lagrangian 3

Isoscalar interaction coupling limit for Lagrangian 13

Isoscalar interaction coupling limit for Lagrangian 16

Isovector interaction coupling limit for Lagrangian 3

Isovector interaction coupling limit for Lagrangian 20

Isovector interaction coupling limit for Lagrangian 4

Isoscalar interaction coupling limit for Lagrangian 14

Isoscalar interaction coupling limit for Lagrangian 12

Isovector interaction coupling limit for Lagrangian 13

Isovector interaction coupling limit for Lagrangian 11

Isoscalar interaction coupling limit for Lagrangian 20

Isoscalar interaction coupling limit for Lagrangian 17

Isoscalar interaction coupling limit for Lagrangian 15

Isovector interaction coupling limit for Lagrangian 6

Isovector interaction coupling limit for Lagrangian 2

Isovector interaction coupling limit for Lagrangian 10

Isovector interaction coupling limit for Lagrangian 16

Electroweak corrections for the $R_{\Delta\phi}(p_\mathrm{T})$ distribution. The multiplicative combination is considered as the main electroweak correction factor and the additive combination as an uncertainty estimate.

Normalized distributions of EA in HT ($E_\mathrm{T}^\mathrm{trig}>4$ GeV) events selected by three ranges of leading jet $p_\mathrm{T}$. Events require $p_\mathrm{T,jet}^\mathrm{raw,lead}>4~\mathrm{GeV}/c$ and the presence of a recoil jet.

Jet spectra per trigger for trigger-side ($|\phi_\mathrm{jet}-\phi_\mathrm{trig}|<\pi/3$) and recoil-side ($|\phi_\mathrm{jet}-\phi_\mathrm{trig}|>2\pi/3$). Jets are of $R=0.4$, and the offline trigger threshold is $E_\mathrm{T}^\mathrm{trig}>4$ GeV. Spectra are shown for high- and low-EA events.

Ratio of the semi-inclusive jet spectra for trigger-side ($|\phi_\mathrm{jet}-\phi_\mathrm{trig}|<\pi/3$) and recoil-side ($|\phi_\mathrm{jet}-\phi_\mathrm{trig}|>2\pi/3$) in high-EA to low-EA events. Jets are of $R=0.4$, and the offline trigger threshold is $E_\mathrm{T}^\mathrm{trig}>4$ GeV.

Distributions of the azimuthal separation between the two highest-$p_\mathrm{T}$ jets for high- and low-EA events for two sets of minimum $p_\mathrm{T}^\mathrm{raw}$.

Ratios of distributions of the azimuthal separation between the two highest-$p_\mathrm{T}$ jets for two sets of minimum $p_\mathrm{T}^\mathrm{raw}$ for high-EA events relative to low-EA events.

Distributions of dijet $p_\mathrm{T}^\mathrm{raw}$ imbalance, $A_{J}$, for high- and low-EA events for two set of minimum $p_\mathrm{T}^\mathrm{raw}$ requirements.

Rastio of distributions of dijet $p_\mathrm{T}^\mathrm{raw}$ imbalance, $A_{J}$, for two sets of minimum $p_\mathrm{T}^\mathrm{raw}$ requirements for high-EA events relative to low-EA events.

Upper 95% CL limits on HH signal cross section scanned over the $\kappa_{\lambda}$ parameter while fixing the $\kappa_{2V}$ and $\kappa_{V}$ to their SM-predicted values.

Upper 95% CL limits on VHH signal cross section scanned over the $\kappa_{2V}$ parameter while fixing the $\kappa_{\lambda}$ and $\kappa_{V}$ to their SM-predicted values.

Upper 95% CL limits on HH signal cross section scanned over the $\kappa_{2V}$ parameter while fixing the $\kappa_{\lambda}$ and $\kappa_{V}$ to their SM-predicted values.

Upper 95% CL limits on VHH signal cross section scanned over the $\kappa_{V}$ parameter while fixing the $\kappa_{2V}$ and $\kappa_{\lambda}$ to their SM-predicted values.

Upper 95% CL limits on HH signal cross section scanned over the $\kappa_{V}$ parameter while fixing the $\kappa_{2V}$ and $\kappa_{\lambda}$ to their SM-predicted values.

The systematic covariance matrix for the $e^+e^- \to \pi^+ \pi^- \pi^0$ cross section measurement at the Belle II. The 212 x 212 matrix of the energy ranges from 0.62 to 3.50 GeV. This covariance matrix includes all systematic uncertainty given in Table I. The total covariance matrix for the measured cross section can be obtained by adding statistical and systematic covariance matrices.

This table provides the ``inverse'' of the total (stat.+ sys.) covariance matrix. This matrix can be used for a chi-square fit of the mass spectrum. We provide this inverse matrix since special care is needed to take the inverse of the given covariance matrix.

Differential cross sections normalized to the fiducial cross section as a function of the $p_T$ of the second-highest-$p_T$ jet

Differential cross sections normalized to the fiducial cross section as a function of the $|\eta|$ of the highest-$p_T$ jet

Differential cross sections normalized to the fiducial cross section as a function of the $|\eta|$ of the second-highest-$p_T$ jet

Differential cross sections normalized to the fiducial cross section as a function of the pseudorapidity difference between the two highest-$p_T$ jets in the event

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the two highest-$p_T$ jets in the event

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 0 jet, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 1 jet, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 2 jets, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with at least 3 jets, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 0 jet, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 1 jet, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 2 jets, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with at least 3 jets, in the full $m_{4l}$ range