Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the muon channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{T}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the invisible channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{T}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the invisible channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the electron channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the electron channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Expected upper limits at 95% CL on the product of the production cross section and the branching fraction as a function of the Z' boson mass.Each line corresponds to the three different final states and their combined result.

Expected upper limits at 95% CL on the product of the production cross section and the branching fraction as a function of the Z' boson mass.Each line corresponds to the three different final states and their combined result.

Expected and observed upper limits at 95% CL on the product of the production cross section and the branching fraction as functions of the Z' boson mass from the combination of all final states. The limits are compared with the predictions of the HVT model and the expected limits (shown by red curves) from a previous analysis.

Expected and observed upper limits at 95% CL on the product of the production cross section and the branching fraction as functions of the Z' boson mass from the combination of all final states. The limits are compared with the predictions of the HVT model and the expected limits (shown by red curves) from a previous analysis.

Prediction from the HVT model.

Prediction from the HVT model.

Observed upper limits at 95% CL on gF for different Z' boson masses as functions of the product of $g_{H}$ with the sign of $g_{F}$. The two benchmark scenarios of the HVT model are shown by the black markers.

Observed upper limits at 95% CL on gF for different Z' boson masses as functions of the product of $g_{H}$ with the sign of $g_{F}$. The two benchmark scenarios of the HVT model are shown by the black markers.

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 95% CL upper limits on the branching fraction $\mathcal{B}(\mathrm{H \to SS})$ with $\mathrm{S}\to\tau\tau$ decay, for different LLP masses $m_{\mathrm{S}}$ and proper decay lengths $c\tau_{0}$.

The 95% CL limits on the LLP mass $m_{\mathrm{S}}$ for different proper decay lengths $c\tau_{0}$ assuming a branching fraction of 1% for the $\mathrm{H \to SS}$ decay, and with subsequent $\mathrm{S \to b\overline{b}}$ decay.

The 95% CL limits on the LLP mass $m_{\mathrm{S}}$ for different proper decay lengths $c\tau_{0}$ assuming a branching fraction of 1% for the $\mathrm{H \to SS}$ decay, and with subsequent $\mathrm{S \to d \overline{d}}$ decay.

The 95% CL limits on the dark-sector top quark partner mass $m_{T}$ for different hidden glueball masses $m_{0}$ in the fraternal twin Higgs model.

The 95% CL limits on the dark-sector top quark partner mass $m_{T}$ for different hidden glueball masses $m_{0}$ in the folded SUSY model.

The 95% CL upper limits on $\sigma(\mathrm{pp}\to\mathrm{Z}')\mathcal{B}(\mathrm{Z}'\to\mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark pion mass $m_{\pi_{\mathrm{dark}}}$ is $20~\mathrm{GeV}$. For the dark sector parton showering and hadronization $N_{f}=7$, $N_{c}=3$, $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2$, and $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=4$ are assumed.

The 95% CL upper limits on $\sigma(\mathrm{pp}\to\mathrm{Z}')\mathcal{B}(\mathrm{Z}'\to\mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark pion mass $m_{\pi_{\mathrm{dark}}}$ is $15~\mathrm{GeV}$. For the dark sector parton showering and hadronization $N_{f}=7$, $N_{c}=3$, $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2$, and $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=4$ are assumed.

The 95% CL upper limits on $\sigma(\mathrm{pp}\to\mathrm{Z}')\mathcal{B}(\mathrm{Z}'\to\mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark pion mass $m_{\pi_{\mathrm{dark}}}$ is $10~\mathrm{GeV}$. For the dark sector parton showering and hadronization $N_{f}=7$, $N_{c}=3$, $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2$, and $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=4$ are assumed.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to SM particles through an offshell Higgs boson. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2.5$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=2.5$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to SM particles through an offshell Higgs boson. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=1$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=2.5$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to SM particles through an offshell Higgs boson. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=1$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=1$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to a pair of gluons. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2.5$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=2.5$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to a pair of gluons. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=1$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=2.5$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=1$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to a pair of gluons. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=1$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=1$.

The 95% CL upper limits on the branching fraction for the $125~\mathrm{GeV}$ Higgs boson to decay to a pair of dark quarks $\mathcal{B}(\mathrm{H}\to \mathrm{Q_{dark}\overline{Q}_{dark}})$ when the dark sector $N_{f}=7$, $N_{c}=3$. The dark pion $\pi_{\mathrm{dark}}$ is assumed to decay to a pair of down quarks. The upper limits are shown for different dark pion masses $m_{\pi_{\mathrm{dark}}}$ and proper decay lengths $c\tau_{\pi_{\mathrm{dark}}}$, with $\Lambda_{\mathrm{dark}}/m_{\pi_{\mathrm{dark}}}=2$, $m_{\rho_{\mathrm{dark}}}/m_{\pi_{\mathrm{dark}}}=4$.

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.

Post-fit distribution of the $M(ll)$ variable for the ultrahigh-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $M(ll)$ variable for the low-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '3l soft' signal region of the the '2/3l soft' analysis.

Post-fit distribution of the $M(ll)$ variable for the medium-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '3l soft' signal region of the the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for low-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for medium-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for high-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

Post-fit distribution of the $m_{\mathrm{T2}}(ll)$ variable for ultrahigh-$p_{\mathrm{T}}^{\mathrm{miss}}$ bins in the '2l soft' signal region of the '2/3l soft' analysis.

2SS $\ell/{\geq}\,3\ell$ search: observed and expected yields across the SRs in category A, events with three light leptons of which at least two form an OSSF pair, after the requirement that the leading-lepton $p_{\mathrm{T}}$ be greater than 30 GeV is applied.

2SS $\ell/{\geq}\,3\ell$ search: observed and expected yields across the SRs of the '${\geq}\ 3\ell$' search in category B, events with three light leptons and no OSSF pair, after the requirement that the leading-lepton $p_{\mathrm{T}}$ be greater than 30 GeV is applied.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production in the WZ topology for the full parameter space.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production in the WZ topology for the compressed space.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production in the WH topology for the full parameter space.

Wino-bino model: cross section limits in the model parameter space, for wino-like chargino-neutralino production with mixed topology with equal branching fraction to WZ and WH.

Wino-bino model: exclusion contours from the individual and combined analyses targeting WZ topology for the full parameter space. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Wino-bino model: exclusion contours from the individual and combined analyses targeting the corresponding compressed region. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Wino-bino model: exclusion contours from the individual and combined analyses targeting the WH topology for the full parameter space. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Wino-bino model: exclusion contours from the individual and combined analyses targeting combined contours for these two topologies. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

GMSB model: expected and observed cross section limits for the neutralino-neutralino production for the ZZ topology.

GMSB model: expected and observed cross section limits for the neutralino-neutralino production for the HH topology.

GMSB model: expected and observed cross section limits for the neutralino-neutralino production for the mixed topology with equal branching fraction to H and Z.

GMSB model: cross section limits for neutralino-neutralino production as a function of the NSLP mass and the branching fraction to the H boson for the combination of the searches.

GMSB model: exclusion limit for neutralino-neutralino production as a function of the NSLP mass and the branching fraction to the H boson for the combination of the searches along with the input searches. For visualization of the exclusion contours, linear interpolation is employed to account for the limited granularity of the available signal samples.

Cross section upper limit(s) in the mass plane of NLSP and LSP masses for the higgsino-bino model.

Mass plane cross section upper limit for direct slepton pair production, with observed and expected exclusion limits in the full mass plane from the combination.

Mass plane cross section upper limit for direct slepton pair production, with observed and expected exclusion limits in the compressed region from '2/3l' soft search.

Correlation matrix for the differential branching fraction extraction between different $q^2$ bins in the simultaneous fit.

The product of acceptance and efficiency ($A\epsilon$) of the $\mathcal{B}(\mathrm{B}^+\to\mathrm{K}^+\mu^+\mu^-)$ channel, as a function of the muon pair $q^2$, as measured in simulated signal events, after all the corrections applied. In the figure, regions corresponding to resonances are displayed with red markers.

Same as Table 5, but for the bins containing the J$\psi$ and $\psi$(2S) resonances.

Integrated branching fraction $\mathcal{B}(\mathrm{B}^\pm \to \mathrm{K}^\pm\mu^+\mu^-)$ in the low-$q^2$ region.

Ratio of branching fractions $\mathcal{B}(\text{B}^\pm \to \text{K}^\pm\mu^+\mu^-)$ to $\mathcal{B}(\text{B}^\pm \to \text{K}^\pm\text{e}^+\text{e}^-)$, determined from the measured double ratio $R$(K) of these decays to the respective branching fractions of the decay chains $\text{B}^\pm\to\text{J/}\psi\text{K}^\pm$ with $\text{J/}\psi\to\mu^+\mu^-$ and $\text{e}^+\text{e}^-$.

Negative log likelihood function from the fit profiled as a function of the inverse ratio $R(\mathrm{K})^{-1}$.

Post-fit distributon of BDT discriminants for $\rho_{tc}=1.0$ with $m_A$ = 350 GeV interfered with H.($m_A - m_H$ = 50 GeV)

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tu}$ = 0.1, 0.4, 1.0.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tc}$ = 0.1, 0.4, 1.0.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tu}$ = 0.1, 0.4, 1.0 with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tc}$ = 0.1, 0.4, 1.0 with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model.

Expected exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model

Observed exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model.

Expected exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model

Observed exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Expected exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Expected exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV