Expected and observed 95% CL upper limits on the branching ratio of the Higgs boson decay to $Za$ times the branching ratio $a $->$ \gamma \gamma$ as a function of the $a$ particle mass in the merged ($m_a \le 2$ GeV) and the resolved ($m_a > 2$ GeV) categories. The figure uses $pp$ collision data at $\sqrt{s}=13$ TeV corresponding to 139 fb$^{-1}$.

ATLAS observed 95% CL exclusion contours limits in terms of the ALP mass and its effective coupling to photons, $|C_{\gamma\gamma}|/\Lambda$, for different values of the Higgs coupling to $Za$, $|C_{Zh}|/\Lambda$. The figure uses $pp$ collision data at $\sqrt{s}=13$ TeV corresponding to 139 fb$^{-1}$. The value 100 in the tables is presentational, denoting the upper edge of the plot.

Distribution of $m_{\mathrm{JJ}}$ in the SR (black points) compared to the background-only expectation, where the shape is derived from the CR and the normalisation is fitted in the SR.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model A, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model B, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d} = 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model C, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d}= 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model D, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Relative efficiency (in %) of each stage of the analysis selection with respect to the previous one for models A through D and $m_{Z^\prime}=2.5$ TeV. Approximately 90000 events were generated per signal mass hypothesis.

Electron $p_{T}$ distribution using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$p_{T}^{miss}$ distribution in the electron channel using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

Muon $p_T$ distribution using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$p_T^{miss}$ distribution in the muon channel using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$M_T$ distribution for the E+INVISIBLE system using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$M_T$ distribution for the MU+INVISIBLE system using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

Theoretical prediction for the production and decay cross-section of SSM WPRIME bosons (W --> LEPTON + NU) as a function of the WPRIME mass, for the range 200-6400 GeV and inclusive in phase-space (no cuts). The uncertainty shown covers PDF and ALPHAS uncertainties. These theoretical predictions are obtained with PYTHIA (LO) and corrected to NNLO in QCD using FEWZ program.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the electron decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the muon decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the combined electron and muon decay channels. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the electron decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the muon decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the combined electron and muon decay channels. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the coupling strength ratio, gWprime/gW as a function of the mass of the WPRIME boson, for the electron channel. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. The dotted line represents the case of SSM coupling, gWprime, being equal to gW, the SM coupling.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the coupling strength ratio, gWprime/gW, as a function of the mass of the WPRIME boson, for the muon channel. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. The dotted line represents the case of SSM coupling, gWprime , being equal to gW, the SM coupling.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the electron channel for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the muon channel for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the combined electron and muon channels for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the ${\lambda}_{231}$ coupling in the RPV SUSY model with a STAU mediator, as a function of the mass of STAU, for the electron channel and for the production coupling, ${\lambda}^{'}_{3ij}$=0.5 . The couplings ${\lambda}^{'}_{3ij}$, ${\lambda}_{231}$, and ${\lambda}_{132}$ are defined in Fig. 1. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. Limits for other values of ${\lambda}^{'}_{3ij}$ are obtained multiplying the ones shown by the ratio ${\lambda}^{'}_{3ij}$/0.5

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the ${\lambda}_{132}$ coupling in the RPV SUSY model with a stau mediator, as a function of the mass of the stau, for the muon channel and for the production coupling, ${\lambda}^{'}_{3ij}$ = 0.5 . The couplings ${\lambda}^{'}_{3ij}$, ${\lambda}_{231}$, and ${\lambda}_{132}$ are defined in Fig. 1. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. Limits for other values of ${\lambda}^{'}_{3ij}$ are obtained multiplying these ones by the ratio ${\lambda}^{'}_{3ij}$/0.5

Region in the oblique (W,Y) parameter phase space allowed by the current analysis at 95% CL. The combination of the electron and muon channel distributions from the 2017 and 2018 data sets at sqrt(s) = 13 TeV is used, having 101 fb-1 of luminosity.

Negative Log-Likelihood for the determination of the oblique W parameter, for the combination of the electron and muon channel and for 2017 and 2018 data sets (L=101 fb-1) at sqrt(s) = 13 TeV.

Total background estimate and observed data in all bins used in the limit calculation for the SSM WPRIME limits. The total event sample is split into the three years of data taking and the two lepton flavors.The distributions being shown are $M_T$ ones. The associated integrated luminosities for each year are of 35.9, 41.5, and 59.4 fb$̣^{-1}. The mass bins are numbered in ascending order.

Parametrization of the $A_{X}$ factor, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level, as function of $m_{X}$ obtained from the narrow width simulated signal samples produced in gluon fusion.

The correction factor, $C_{X}$, defined as the ratio between the number of reconstructed signal events passing the analysis cuts and the number of signal events at the particle level generated within the fiducial volume, and acceptance correction factor, $A_{X}$, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level. Both are computed for NWA spin-0 models as a function of $m_{X}$.

Effect of event selections on a scalar MC signal sample generated for $m_{X}$ = 15 GeV and on the data. For the MC sample, the efficiencies are shown after applying event weights and a truth level filter that requires two photons with $p^{\gamma\gamma}_{T}>40$ GeV; for the data, the absolute yields are shown. The initial yields for data include a trigger preselection that is the OR of a list of single photon and diphoton triggers. The "2 $loose$ photons" step includes the kinematic acceptance cuts.

Parameterization of the Double Sided Crystal Ball function parameters describing the scalar mass resolution model as a function of $m_{X}$ [GeV].

Expected and observed limits at 95% CL on the branching fraction B(hb --> PHI X) . B(PHI --> mu+mu-) as a function of signal mass for signal lifetime of 100 mm

Expected and observed limits at 95% CL on the branching fraction B(H --> ZD ZD) . B(ZD --> mu+mu-) as a function of signal mass for signal lifetime of 1 mm

Expected and observed limits at 95% CL on the branching fraction B(H --> ZD ZD) . B(ZD --> mu+mu-) as a function of signal mass for signal lifetime of 1 mm

Expected and observed limits at 95% CL on the branching fraction B(H --> ZD ZD) . B(ZD --> mu+mu-) as a function of signal mass for signal lifetime of 100 mm

Expected and observed limits at 95% CL on the branching fraction B(H --> ZD ZD) . B(ZD --> mu+mu-) as a function of signal mass for signal lifetime of 100 mm

Expected and observed limits at 95% CL on the branching fraction B(H --> ZD ZD) as a function of signal mass for signal lifetime of 1 mm

Expected and observed limits at 95% CL on the branching fraction B(H --> ZD ZD) as a function of signal mass for signal lifetime of 1 mm

Expected and observed limits at 95% CL on the branching fraction B(H --> ZD ZD) as a function of signal mass for signal lifetime of 100 mm

Expected and observed limits at 95% CL on the branching fraction B(H --> ZD ZD) as a function of signal mass for signal lifetime of 100 mm

Observed limits at 95% CL on the branching fraction B(hb --> PHI X) . B(PHI --> mu+mu-) in the parameter space containing signal mass and signal lifetime. The ROOT file containing the values are available as additional resource

Observed limits at 95% CL on the branching fraction B(hb --> PHI X) . B(PHI --> mu+mu-) in the parameter space containing signal mass and signal lifetime. The ROOT file containing the values is available as additional resource

Observed limits at 95% CL on the branching fraction B(H --> ZD ZD) in the parameter space containing signal mass and signal lifetime. The ROOT file containing the values are available as additional resource

Observed limits at 95% CL on the branching fraction B(H --> ZD ZD) in the parameter space containing signal mass and signal lifetime. The ROOT file containing the values is available as additional resource

Region excluded at 95% CL using this search for the inflaton model in the parameter space containing signal mass and square of the signal mixing angle

Region excluded at 95% CL using this search for the inflaton model in the parameter space containing signal mass and square of the signal mixing angle

Observed limits at 95% CL on the branching fraction B(H --> ZD ZD) in the parameter space containing signal mass and kinetic mixing parameter. The ROOT file containing the values are available as additional resource

Observed limits at 95% CL on the branching fraction B(H --> ZD ZD) in the parameter space containing signal mass and kinetic mixing parameter. The ROOT file containing the values is available as additional resource

Observed limits at 95% CL on the branching fraction B(hb --> PHI X) . B(PHI --> mu+mu-) as a function of signal lifetime for various representative signal mass hypotheses

Observed limits at 95% CL on the branching fraction B(H --> ZD ZD) as a function of signal lifetime for various representative signal mass hypotheses

CMS and LHCb observed limits at 95% CL on the branching fraction B(hb --> PHI X) . B(PHI --> mu+mu-) as a function of signal mass for signal lifetime of 1 mm

CMS and LHCb observed limits at 95% CL on the branching fraction B(hb --> PHI X) . B(PHI --> mu+mu-) as a function of signal mass for signal lifetime of 10 mm

CMS and LHCb observed limits at 95% CL on the branching fraction B(hb --> PHI X) . B(PHI --> mu+mu-) as a function of signal mass for signal lifetime of 100 mm

Observed limits at 95% CL on number of signal events in aggregate signal region 0.2 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 0 GeV and no requirement on muon isolation

Observed limits at 95% CL on number of signal events in aggregate signal region 0.2 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 0 GeV and two isolated muons

Observed limits at 95% CL on number of signal events in aggregate signal region 0.2 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 25 GeV and no requirement on muon isolation

Observed limits at 95% CL on number of signal events in aggregate signal region 0.2 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 25 GeV and two isolated muons

Observed limits at 95% CL on number of signal events in aggregate signal region 1.0 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 0 GeV and no requirement on muon isolation

Observed limits at 95% CL on number of signal events in aggregate signal region 1.0 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 0 GeV and two isolated muons

Observed limits at 95% CL on number of signal events in aggregate signal region 1.0 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 25 GeV and no requirement on muon isolation

Observed limits at 95% CL on number of signal events in aggregate signal region 1.0 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 25 GeV and two isolated muons

Observed limits at 95% CL on number of signal events in aggregate signal region 2.4 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 0 GeV and no requirement on muon isolation

Observed limits at 95% CL on number of signal events in aggregate signal region 2.4 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 0 GeV and two isolated muons

Observed limits at 95% CL on number of signal events in aggregate signal region 2.4 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 25 GeV and no requirement on muon isolation

Observed limits at 95% CL on number of signal events in aggregate signal region 2.4 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 25 GeV and two isolated muons

Observed limits at 95% CL on number of signal events in aggregate signal region 3.1 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 0 GeV and no requirement on muon isolation

Observed limits at 95% CL on number of signal events in aggregate signal region 3.1 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 0 GeV and two isolated muons

Observed limits at 95% CL on number of signal events in aggregate signal region 3.1 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 25 GeV and no requirement on muon isolation

Observed limits at 95% CL on number of signal events in aggregate signal region 3.1 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 25 GeV and two isolated muons

Observed limits at 95% CL on number of signal events in aggregate signal region 7.0 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 0 GeV and no requirement on muon isolation

Observed limits at 95% CL on number of signal events in aggregate signal region 7.0 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 0 GeV and two isolated muons

Observed limits at 95% CL on number of signal events in aggregate signal region 7.0 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 25 GeV and no requirement on muon isolation

Observed limits at 95% CL on number of signal events in aggregate signal region 7.0 cm $\leq$ lxy < 11.0 cm, dimuon pt $\geq$ 25 GeV and two isolated muons

Measured DPS effective cross section

pp -> J/psi J/psi J/psi X cross section

$pp \rightarrow J/\psi J/\psi J/\psi X~$ fiducial cross section

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Expected limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}$( B -> bH ) and $\mathcal{B}$( B -> bZ)

Observed limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}$( B -> bH ) and $\mathcal{B}$( B -> bZ)

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.