The distributions of the variables used in the simultaneous fit for the SR. The black points with error bars represent the data and their statistical uncertainties, whereas the shaded band represents the predicted uncertainties. The bottom panel in each figure shows the ratio of the number of events observed in data to that of the total SM prediction. The last bin of each plot has been extended to include the overflow contribution.

Expected and observed 95% upper limits on the product of the cross section and branching fraction$\sigma(\mathrm{pp} \ ightarrow W\mathrm{a})\mathcal{B}(W ightarrow \ell^{+} v_{\ell})\mathcal{B}(\mathrm{a} \rightarrow Z\gamma)\mathcal{B}(Z \rightarrow \ell^{+} \ell^{-})$ s a function of the ALP mass from 110 to 400 GeV

Expected and observed 95% upper limits on the photophobic ALP model parameter $1/f_{\mathrm{a}}$ as a function of ALP mass reinterpreted from $1/f_{\mathrm{a}}=2~\mathrm{TeV}^{-1}$

Comparison of the shapes extrapolated with a linear or quadratic extrapolation of the QCD contribution in the W+ signal region.

Comparison of the shapes extrapolated with a linear or quadratic extrapolation of the QCD contribution in the W- signal region.

Post-fit distributions of the transverse mass in the W+ signal region.

Post-fit distributions of the transverse mass in the W- signal region.

Post-fit distributions of the dimuon mass in the Z boson signal region.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W+ signal region in the transverse mass bin 12 < mT < 18 GeV. The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W- signal region in the transverse mass bin 12 < mT < 18 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W+ signal region in the transverse mass bin 102 < mT < 108 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (larger value) and quadratic (smaller value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W- signal region in the transverse mass bin 102 < mT < 108 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

The table compares theoretical and measured values of the product of the fiducial cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value. The relative uncertainties on the measured values are given separately with and without the contribution from the luminosity uncertainty.

The table compares theoretical and measured values of the product of the total cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value. The relative uncertainties on the measured values are given separately with and without the contribution from the luminosity uncertainty.

The table compares theoretical and measured ratios of the product of the fiducial cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value.

The table compares theoretical and measured ratios of the product of the total cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value.

Comparison of the measured and theoretical cross sections for Z, W+, W-, and W production at different center-of-mass energies. The theoretical predictions are generated with DYTurbo in combination with the NNPDF 3.1 PDF set.

Expected and observed upper limits for the b-quark-associated Higgs boson production cross section times branching fraction of the decay into a b quark pair at 95% CL as functions of $m_\phi$ for the 2018 FH category. The vertical dashed lines indicate the boundaries of usage of the different fit ranges, as reflected in the rightmost column of Table 2.

Expected and observed upper limits for the b-quark-associated Higgs boson production cross section times branching fraction of the decay into a b quark pair at 95% CL as functions of $m_\phi$, corresponding to the Run 2 combination. The vertical dashed line separates the mass range where only the 2017 SL category contributes on its left, from the region where also the 2017 FH and 2018 FH categories contribute on its right. The expected limits from the 2017 SL, 2017 FH, and 2018 FH datasets as well as from the previously published result based on the 2016 dataset are also shown as colored lines.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = +1$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = -1$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = -2$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = -3$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $m_h^{mod+}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the hMSSM scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as a function of $m_{A}$ for cos(beta-alpha) = 0.1, for the 2HDM $\it Type-II$ scenario.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for masses of $m_{A} = m_{H} = 300$ GeV, for the 2HDM $\it Type-II$ scenario.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as a function of $m_{A}$ for cos(beta-alpha) = 0.1, for the 2HDM $\it Flipped$ scenario.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for masses of $m_{A} = m_{H} = 300$ GeV, for the 2HDM $\it Flipped$ scenario.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 140$ GeV.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 600$ GeV.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 900$ GeV.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 1200$ GeV.

Simulated signal shapes for three representative values of the Higgs boson mass $m_\phi$ in the 2017 SL channel.

Simulated signal shapes for three representative values of the Higgs boson mass $m_\phi$ in the 2017 FH channel.

Simulated signal shapes for three representative values of the Higgs boson mass $m_\phi$ in the 2018 FH channel.

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the SL category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the SL category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the SL category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

Distributions for data and expected signal and background contributions of the MVA score for the µ + jets channel in the 3j1b category, before the maximum likelihood fit. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty in the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-prediction ratio. The first and last bins in each distribution include underflow and overflow events, respectively.

Observed and predicted number of events in each of the eight categories of the signal region, before the maximum likelihood fit. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty in the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-prediction ratio.

Distributions for the e + jets final state before the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization before the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the e + jets final state after the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization after the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the µ + jets final state before the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization before the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the µ + jets final state after the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization after the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Covariance matrix for the lepton+jets ML fit.

Covariance matrix for the combination ML fit.

The $G_{\text{i}}^{\text{Strips}}$ distribution in the PASS region for events passing the event selection and with $55 < p_{\mathrm{T}} < 200$ GeV.

The $G_{\text{i}}^{\text{Strips}}$ distribution in the FAIL region for events passing the event selection and with $p_{\mathrm{T}} > 200$ GeV.

The $G_{\text{i}}^{\text{Strips}}$ distribution in the PASS region for events passing the event selection and with $p_{\mathrm{T}} > 200$ GeV.

Mass spectrum predicted in the signal region defined by $G_{\text{i}}^{\text{Strips}} > 0.22$ and $p_{\mathrm{T}} > 70$ GeV.

Cross section limits for gluino and supersymmetric top R-hadrons for both background prediction methods.

Cross section limits for supersymmetric tau models for both background prediction methods.

Cross section limits for DY-produced tau prime models (single and multicharged) for both background prediction methods.

Cross section limits for Z prime to multicharged tau prime model for both background prediction methods. All Z prime models assume a narrow width, and a 100% branching fraction to 600 GeV tau primes.

2D exclusion showing the observed cross section limit as a function of the multicharged tau prime mass and Z prime mass for the ionization method.

2D exclusion showing the observed cross section limit as a function of the multicharged tau prime mass and Z prime mass for the mass method.

Cumulative selection efficiency for the data and for two signal hypotheses.

Expected and observed mass limits obtained using 2017-2018 data for various HSCP candidate models,for the two background estimate methods.

Mass windows used in the mass method as a function of the signal target mass for the signal samples assuming a charge of $1e$ and $2e$, as appropriate to the signal model.

Mass spectrum predicted in the validation region defined by $0.018 < G_{\text{i}}^{\text{Strips}} < 0.057$ and $p_{\mathrm{T}} > 70$ GeV. The thresholds used in the $G_{\text{i}}^{\text{Strips}}$ requirement represent the 50% and 90% quantile of the distribution.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $\text{abs(}\eta\text{)}<0.3$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.3<\text{abs(}\eta\text{)}<0.6$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.6<\text{abs(}\eta\text{)}<0.9$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.9<\text{abs(}\eta\text{)}<1.2$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $1.2<\text{abs(}\eta\text{)}<2.1$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $2.1<\text{abs(}\eta\text{)}<2.4$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $\text{abs(}\eta\text{)}<0.3$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.3<\text{abs(}\eta\text{)}<0.6$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.6<\text{abs(}\eta\text{)}<0.9$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.9<\text{abs(}\eta\text{)}<1.2$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $1.2<\text{abs(}\eta\text{)}<2.1$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $2.1<\text{abs(}\eta\text{)}<2.4$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency as a function of $\beta$, for the $\tilde{g}$ signals.

Trigger efficiency as a function of $\beta$, for the $\tilde{\tau}$ signals.

Cross section limits for $\tilde{g}$ R-hadrons obtained with the ionization method.

Cross section limits for $\tilde{g}$ R-hadrons obtained with the mass method.

Cross section limits for $\tilde{t}$ R-hadrons obtained with the ionization method.

Cross section limits for $\tilde{t}$ R-hadrons obtained with the mass method.

Cross section limits for $\tilde{\tau}$ obtained with the ionization method.

Cross section limits for $\tilde{\tau}$ obtained with the mass method.

Cross section limits for $\tilde{\tau}$ production within the GMSB SPS7 modelobtained with the ionization method.

Cross section limits for $\tilde{\tau}$ production within the GMSB SPS7 model obtained with the mass method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 1e$ with the ionization method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 1e$ with the mass method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 2e$ with the ionization method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 2e$ with the mass method.

Cross section limits for the production of $Z'$ boson decaying into a pair of $\tau'$ fermions of charge $2e$ (with a branching fraction equal to 1 and a fixed $\tau'$ mass of 600 GeV), obtained with the ionization method.

Cross section limits for the production of $Z'$ boson decaying into a pair of $\tau'$ fermions of charge $2e$ (with a branching fraction equal to 1 and a fixed $\tau'$ mass of 600 GeV), obtained with the mass method.

Two-dimensional exclusion showing the observed cross section limit as a function of the masses of the $\tau'$ (on the $x$ axis) and of the $Z'$ boson (on the $y$ axis), for the ionization method.

Two-dimensional exclusion showing the observed cross section limit as a function of the masses of the $\tau'$ (on the $x$ axis) and of the $Z'$ boson (on the $y$ axis), for the mass method.

Expected background yield, expected signal yield, and observed data for the mass method for 2017.

Expected background yield, expected signal yield, and observed data for the mass method for 2018.

Selection efficieny for GMSB supersymmetric $\tau$.

Selection efficiency for gluino R-hardon samples.

Selection efficiency for supersymmetric top R-hadrons.

Selection efficiency for pair-produced supersymmetric $\tau$.

Selection efficiency for $\tau'$ with $|Q| = 1e$.

Selection efficiecny for $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ (mass = 3 TeV) going to $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ (mass 5 TeV) goes to $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ going to $\tau'$ (mass 600 GeV) with $|Q| = 2e$.

Selection efficiency for $Z'$ going to $\tau'$ (mass 1 TeV) with charge $|Q| = 2e$.

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.

Observed limits at the 95% CL on $C_{\mathrm{LL}}^{2332}$ vs. $C_{\mathrm{LR}}^{2332}$ (dark blue) showing the full range.

Observed limits at the 95% CL on $C_{\mathrm{LR}}^{2233}$ vs. $C_{\mathrm{LL}}^{2233}$ (light blue) showing the full range.

Observed limits at the 95% CL on $C_{\mathrm{LL}}^{2233}$ vs. $C_{\mathrm{LR}}^{2332}$ (red) showing an enlarged view of the region around values of 0.

Observed limits at the 95% CL on $C_{\mathrm{LR}}^{2233}$ vs. $C_{\mathrm{LL}}^{2332}$ (orange) showing an enlarged view of the region around values of 0.

Observed limits at the 95% CL on $C_{\mathrm{LL}}^{2332}$ vs. $C_{\mathrm{LR}}^{2332}$ (dark blue) showing an enlarged view of the region around values of 0.

The observed exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The expected exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The expected exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The expected upper $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The expected upper $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The expected lower $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The expected lower $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The $95\%$ CL observed upper limits on the cross section shown as a functions of $T_D$ and the pseudoscalar mass $m_{\phi}$ for $m_S = 300$ GeV, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The $95\%$ CL observed upper limits on the cross section shown as a functions of $T_D$ and the pseudoscalar mass $m_{\phi}$ for $m_S = 300$ GeV, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The $95\%$ CL observed upper limits on the cross section shown as a functions of $T_D$ and the pseudoscalar mass $m_{\phi}$ for $m_S = 600$ GeV, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The $95\%$ CL observed upper limits on the cross section shown as a functions of $T_D$ and the pseudoscalar mass $m_{\phi}$ for $m_S = 600$ GeV, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The $95\%$ CL observed upper limits on the cross section shown as a functions of $T_D$ and the pseudoscalar mass $m_{\phi}$ for $m_S = 1000$ GeV, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The $95\%$ CL observed upper limits on the cross section shown as a functions of $T_D$ and the pseudoscalar mass $m_{\phi}$ for $m_S = 1000$ GeV, for the case $m_{A'}=1.0$ GeV ($A' \rightarrow \pi^+\pi^-$ with $\mathcal{BR}=100\%$).

The observed exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The observed exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The expected exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The expected exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The expected upper $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The expected upper $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The expected lower $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The expected lower $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.5$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 40, 40, 20\%$).

The observed exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

The observed exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

The expected exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

The expected exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

The expected upper $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

The expected upper $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

The expected lower $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

The expected lower $68\%$ exclusions for the nominal $\sigma$ cross section in the plane of $m_{\phi}$ and $T_D$, for various $m_S$ values, for the case $m_{A'} = 0.7$ GeV ($A' \rightarrow e^+e^-, \mu^+\mu^-, \pi^+\pi^-$ with $\mathcal{BR} = 15, 15, 70\%$).

For a selection of signal samples in era 2018 with $m_{\phi} = T_D = 3$ GeV, the relative efficiencies and the total efficiency with respect to the selection on the generator-level $H_T$,$H^{gen}_T > 1000$ GeV.

For a selection of signal samples in era 2018 with $m_{\phi} = T_D = 3$ GeV, the relative efficiencies and the total efficiency with respect to the selection on the generator-level $H_T$,$H^{gen}_T > 1000$ GeV.

The $95\%$ CL cross section limit for all values of $m_{A'}$, $T_D$, $m_{\phi}$, and $m_S$ used in the scan.

The $95\%$ CL cross section limit for all values of $m_{A'}$, $T_D$, $m_{\phi}$, and $m_S$ used in the scan.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bHbH channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bZbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bHtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 5-jet bZtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bHbH channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bZbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bHtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the hadronic 6-jet bZtW channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 3-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 3-jet bZbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 4-jet bHbZ channel.

Distributions of reconstructed VLQ mass for expected postfit background (blue histogram), signal plus background (colored lines), and observed data (black points) for events in the semileptonic 4-jet bZbZ channel.

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 0% $\mathcal{B}(B \to bH)$, 100% $\mathcal{B}(B \to bH)$, and 0% $\mathcal{B}(B \to bH)$

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 25% $\mathcal{B}(B \to bH)$, 25% $\mathcal{B}(B \to bH)$, and 50% $\mathcal{B}(B \to bH)$

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 50% $\mathcal{B}(B \to bH)$, 50% $\mathcal{B}(B \to bH)$, and 0% $\mathcal{B}(B \to bH)$

The limit at 95% CL on the cross section for VLQ pair production for the branching fraction hypothesis 100% $\mathcal{B}(B \to bH)$, 0% $\mathcal{B}(B \to bH)$, and 0% $\mathcal{B}(B \to bH)$

Median expected exclusion limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}(B \to bH)$ and $\mathcal{B}(B \to tW)$, with $\mathcal{B}(B \to tW) = 1 - \mathcal{B}(B \to bH) - \mathcal{B}(B \to bZ)$. The grey area corresponds to the region where the exclusion limit is less than 1000 GeV.

Median observed exclusion limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}(B \to bH)$ and $\mathcal{B}(B \to tW)$, with $\mathcal{B}(B \to tW) = 1 - \mathcal{B}(B \to bH) - \mathcal{B}(B \to bZ)$. The grey area corresponds to the region where the exclusion limit is less than 1000 GeV.