The ratio between the observed data and the simulation in the high $\Delta$m portion of the $\ell+\text{jets}$ control region as a function of $p_T^{miss}$ after scaling the simulation to match the total yield in data.

Comparison between data and simulation in the high $\Delta$m portion of the $\ell+\text{jets}$ control region as a function of $N_t$ after scaling the simulation to match the total yield in data. The hatched region indicates the total shape uncertainty in the simulation.

The ratio between the observed data and the simulation in the high $\Delta$m portion of the $\ell+\text{jets}$ control region as a function of $N_t$ after scaling the simulation to match the total yield in data.

Comparison between data and simulation in the high $\Delta$m portion of the $\ell+\text{jets}$ control region as a function of $N_W$ after scaling the simulation to match the total yield in data. The hatched region indicates the total shape uncertainty in the simulation.

The ratio between the observed data and the simulation in the high $\Delta$m portion of the $\ell+\text{jets}$ control region as a function of $N_W$ after scaling the simulation to match the total yield in data.

Comparison between data and simulation in the high $\Delta$m portion of the $\ell+\text{jets}$ control region as a function of $N_{\text{res}}$ after scaling the simulation to match the total yield in data. The hatched region indicates the total shape uncertainty in the simulation.

The ratio between the observed data and the simulation in the high $\Delta$m portion of the $\ell+\text{jets}$ control region as a function of $N_{\text{res}}$ after scaling the simulation to match the total yield in data.

Observed event yields in data (black points) and predicted SM background (filled histograms) for the low $\Delta$m search bins 0--52. The signal models are denoted in the legend with the masses in GeV of the SUSY particles in parentheses: $(m_{\tilde{t}}, m_{\tilde{\chi}^0_1})$ or $(m_{\tilde{g}}, m_{\tilde{\chi}^0_1})$ for the T2 or T1 signal models, respectively. The hatched bands correspond to the total uncertainty in the background prediction. The (unstacked) distributions for two example signal models are also shown.

The ratio of the data to the total background prediction for the low $\Delta$m search bins 0--52. The hatched bands correspond to the total uncertainty in the background prediction.

Observed event yields in data (black points) and predicted SM background (filled histograms) for the high $\Delta$m search bins 53--104. The signal models are denoted in the legend with the masses in GeV of the SUSY particles in parentheses: $(m_{\tilde{t}}, m_{\tilde{\chi}^0_1})$ or $(m_{\tilde{g}}, m_{\tilde{\chi}^0_1})$ for the T2 or T1 signal models, respectively. The hatched bands correspond to the total uncertainty in the background prediction. The (unstacked) distributions for two example signal models are also shown.

The ratio of the data to the total background prediction for the high $\Delta$m search bins 53--104. The hatched bands correspond to the total uncertainty in the background prediction.

Observed event yields in data (black points) and predicted SM background (filled histograms) for the high $\Delta$m search bins 105--152 with ${N_b = 2}$. The signal models are denoted in the legend with the masses in GeV of the SUSY particles in parentheses: $(m_{\tilde{t}}, m_{\tilde{\chi}^0_1})$ or $(m_{\tilde{g}}, m_{\tilde{\chi}^0_1})$ for the T2 or T1 signal models, respectively. The hatched bands correspond to the total uncertainty in the background prediction. The (unstacked) distributions for two example signal models are also shown.

The ratio of the data to the total background prediction for the high $\Delta$m search bins 105--152 with ${N_b = 2}$. The hatched bands correspond to the total uncertainty in the background prediction.

Observed event yields in data (black points) and predicted SM background (filled histograms) for the high $\Delta$m search bins 153--182 with ${N_b \geq 3}$. The signal models are denoted in the legend with the masses in GeV of the SUSY particles in parentheses: $(m_{\tilde{t}}, m_{\tilde{\chi}^0_1})$ or $(m_{\tilde{g}}, m_{\tilde{\chi}^0_1})$ for the T2 or T1 signal models, respectively. The hatched bands correspond to the total uncertainty in the background prediction. The (unstacked) distributions for two example signal models are also shown.

The ratio of the data to the total background prediction for the high $\Delta$m search bins 153--182 with ${N_b \geq 3}$. The hatched bands correspond to the total uncertainty in the background prediction.

The observed 95% CL upper limit on the production cross section of the T2tt simplified model as a function of the top squark and LSP masses. No interpretation is provided for signal models for which ${|{m_{\tilde{t}} - m_{\tilde{\chi}^0_1} - m_t}| < 25 GeV}$ and ${m_{\tilde{t}} < 275 GeV}$ as described in the text.

The expected 95% CL upper limit on the production cross section of the T2tt simplified model as a function of the top squark and LSP masses. No interpretation is provided for signal models for which ${|{m_{\tilde{t}} - m_{\tilde{\chi}^0_1} - m_t}| < 25 GeV}$ and ${m_{\tilde{t}} < 275 GeV}$ as described in the text.

The observed exclusion contour of the T2tt simplified model with respect to approximate NNLO+NNLL signal cross sections and the change in this contour due to variation of these cross sections within their theoretical uncertainties ($\sigma_{\text{theory}}$). No interpretation is provided for signal models for which ${|{m_{\tilde{t}} - m_{\tilde{\chi}^0_1} - m_t}| < 25 GeV}$ and ${m_{\tilde{t}} < 275 GeV}$ as described in the text.

The mean expected exclusion contour of the T2tt simplified model and the region containing 68 and 95\% ($\pm 1$ and $2\,\sigma_{\text{experiment}}$) of the distribution of expected exclusion limits under the background-only hypothesis. No interpretation is provided for signal models for which ${|{m_{\tilde{t}} - m_{\tilde{\chi}^0_1} - m_t}| < 25 GeV}$ and ${m_{\tilde{t}} < 275 GeV}$ as described in the text.

The observed 95% CL upper limit on the production cross section of the T2bW simplified model as a function of the top squark and LSP masses.

The expected 95% CL upper limit on the production cross section of the T2bW simplified model as a function of the top squark and LSP masses.

The observed exclusion contour of the T2bW simplified model with respect to approximate NNLO+NNLL signal cross sections and the change in this contour due to variation of these cross sections within their theoretical uncertainties ($\sigma_{\text{theory}}$).

The mean expected exclusion contour of the T2bW simplified model and the region containing 68 and 95\% ($\pm 1$ and $2\,\sigma_{\text{experiment}}$) of the distribution of expected exclusion limits under the background-only hypothesis.

The observed 95% CL upper limit on the production cross section of the T2tb simplified model as a function of the top squark and LSP masses.

The expected 95% CL upper limit on the production cross section of the T2tb simplified model as a function of the top squark and LSP masses.

The observed exclusion contour of the T2tb simplified model with respect to approximate NNLO+NNLL signal cross sections and the change in this contour due to variation of these cross sections within their theoretical uncertainties ($\sigma_{\text{theory}}$).

The mean expected exclusion contour of the T2tb simplified model and the region containing 68 and 95\% ($\pm 1$ and $2\,\sigma_{\text{experiment}}$) of the distribution of expected exclusion limits under the background-only hypothesis.

The observed 95% CL upper limit on the production cross section of the T2ttC simplified model as a function of the top squark mass and the difference between the top squark and LSP masses.

The expected 95% CL upper limit on the production cross section of the T2ttC simplified model as a function of the top squark mass and the difference between the top squark and LSP masses.

The observed exclusion contour of the T2ttC simplified model with respect to approximate NNLO+NNLL signal cross sections and the change in this contour due to variation of these cross sections within their theoretical uncertainties ($\sigma_{\text{theory}}$).

The mean expected exclusion contour of the T2ttC simplified model and the region containing 68\% ($\pm 1\,\sigma_{\text{experiment}}$) of the distribution of expected exclusion limits under the background-only hypothesis.

The observed 95% CL upper limit on the production cross section of the T2bWC simplified model as a function of the top squark mass and the difference between the top squark and LSP masses.

The expected 95% CL upper limit on the production cross section of the T2bWC simplified model as a function of the top squark mass and the difference between the top squark and LSP masses.

The observed exclusion contour of the T2bWC simplified model with respect to approximate NNLO+NNLL signal cross sections and the change in this contour due to variation of these cross sections within their theoretical uncertainties ($\sigma_{\text{theory}}$).

The mean expected exclusion contour of the T2bWC simplified model and the region containing 68\% ($\pm 1\,\sigma_{\text{experiment}}$) of the distribution of expected exclusion limits under the background-only hypothesis.

The observed 95% CL upper limit on the production cross section of the T2cc simplified model as a function of the top squark mass and the difference between the top squark and LSP masses.

The expected 95% CL upper limit on the production cross section of the T2cc simplified model as a function of the top squark mass and the difference between the top squark and LSP masses.

The observed exclusion contour of the T2cc simplified model with respect to approximate NNLO+NNLL signal cross sections and the change in this contour due to variation of these cross sections within their theoretical uncertainties ($\sigma_{\text{theory}}$).

The mean expected exclusion contour of the T2cc simplified model and the region containing 68\% ($\pm 1\,\sigma_{\text{experiment}}$) of the distribution of expected exclusion limits under the background-only hypothesis.

The observed 95% CL upper limit on the production cross section of the T1tttt simplified model as a function of the gluino and LSP masses.

The expected 95% CL upper limit on the production cross section of the T1tttt simplified model as a function of the gluino and LSP masses.

The observed exclusion contour of the T1tttt simplified model with respect to approximate NNLO+NNLL signal cross sections and the change in this contour due to variation of these cross sections within their theoretical uncertainties ($\sigma_{\text{theory}}$).

The mean expected exclusion contour of the T1tttt simplified model and the region containing 68 and 95\% ($\pm 1$ and $2\,\sigma_{\text{experiment}}$) of the distribution of expected exclusion limits under the background-only hypothesis.

The observed 95% CL upper limit on the production cross section of the T1ttbb simplified model as a function of the gluino and LSP masses.

The expected 95% CL upper limit on the production cross section of the T1ttbb simplified model as a function of the gluino and LSP masses.

The observed exclusion contour of the T1ttbb simplified model with respect to approximate NNLO+NNLL signal cross sections and the change in this contour due to variation of these cross sections within their theoretical uncertainties ($\sigma_{\text{theory}}$).

The mean expected exclusion contour of the T1ttbb simplified model and the region containing 68 and 95\% ($\pm 1$ and $2\,\sigma_{\text{experiment}}$) of the distribution of expected exclusion limits under the background-only hypothesis.

The observed 95% CL upper limit on the production cross section of the T5ttcc simplified model as a function of the gluino and LSP masses. The upper limits do not take into account contributions from direct top squark pair production; however, its effect is small for $m_{\tilde{\chi}^0_1} > 600 GeV$, which corresponds to the phase space beyond the exclusions based on direct top squark pair production. The excluded regions based on direct top squark pair production from this search and earlier searches are indicated by the hatched areas.

The expected 95% CL upper limit on the production cross section of the T5ttcc simplified model as a function of the gluino and LSP masses. The uppser limits do not take into account contributions from direct top squark pair production; however, its effect is small for $m_{\tilde{\chi}^0_1} > 600 GeV$, which corresponds to the phase space beyond the exclusions based on direct top squark pair production. The excluded regions based on direct top squark pair production from this search and earlier searches are indicated by the hatched areas.

The observed exclusion contour of the T5ttcc simplified model with respect to approximate NNLO+NNLL signal cross sections and the change in this contour due to variation of these cross sections within their theoretical uncertainties ($\sigma_{\text{theory}}$). The expected and observed upper limits do not take into account contributions from direct top squark pair production; however, its effect is small for $m_{\tilde{\chi}^0_1} > 600 GeV$, which corresponds to the phase space beyond the exclusions based on direct top squark pair production. The excluded regions based on direct top squark pair production from this search and earlier searches are indicated by the hatched areas.

The mean expected exclusion contour of the T5ttcc simplified model and the region containing 68% and 95% ($\pm 1$ and $2\,\sigma_{\text{experiment}}$) of the distribution of expected exclusion limits under the background-only hypothesis. The expected and observed upper limits do not take into account contributions from direct top squark pair production; however, its effect is small for $m_{\tilde{\chi}^0_1} > 600 GeV$, which corresponds to the phase space beyond the exclusions based on direct top squark pair production. The excluded regions based on direct top squark pair production from this search and earlier searches are indicated by the hatched areas.

The product of signal acceptance and efficiency in the 2l categories for the signal produced via vector boson fusion.

$m_{X}^{T}$ distribution in data in the 0l 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the non-VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis. The CMS search for a heavy resonance using 2016 data and the same final state [JHEP 11 (2018) 172] is shown as a comparison.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 2 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 3 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 4 TeV for the non-VBF signal.

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.

The $L_{T}$ distribution in the MisID-enriched region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 3L below-Z signal region. The last bin contains the overflow events.

The $M_{T}$ distribution in the 3L on-Z signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 3L above-Z signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 3L OSSF0 signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 4L OSSF0 signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 4L OSSF1 signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 4L OSSF2 signal region. The last bin contains the overflow events.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 4L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 4L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 4L(ee) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 4L(ee) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 4L(ee) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 4L(ee) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 4L($\mu\mu$) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 4L($\mu\mu$) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The 95% confidence level exclusion limits for the flavor-democratic scenario on the total production cross section of heavy fermion pairs.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(S) and $\mathcal{B}(\phi - {ee})$.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(PS) and $\mathcal{B}(\phi - {ee})$.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(S) and $\mathcal{B}(\phi - {\mu\mu})$.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(PS) and $\mathcal{B}(\phi - {\mu\mu})$.

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {ee})$ for tt$\phi$(S).

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {ee})$ for tt$\phi$(PS).

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$ for tt$\phi$(S).

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$ for tt$\phi$(PS).

Product of the fiducial acceptance and the event selection efficiency for the type-III Seesaw signal model at various signal mass hypotheses calculated after all analysis selection requirements.

Product of the fiducial acceptance and the event selection efficiency for the tt$\phi$ signal models at various signal mass hypotheses calculated after all analysis selection requirements.

Observed yields and pre-fit background predictions for Njets 8-9.

Observed yields and pre-fit background predictions for Njets >=10.

Observed yields and pre-fit background predictions for aggregated bins.

Pre-fit background covariance matrix $\sigma_{xy}$ for the 174 analysis bins.

Pre-fit background correlation matrix $\rho_{xy}$ for the 174 analysis bins.

Observed yields and post-fit background predictions for Njets 2-3.

Observed yields and post-fit background predictions for Njets 4-5.

Observed yields and post-fit background predictions for Njets 6-7.

Observed yields and post-fit background predictions for Njets 8-9.

Observed yields and post-fit background predictions for Njets >=10.

Simulated signal yields for representative gluino signal model light-LSP points in the aggregate search regions. The uncertainties shown are statistical.

Simulated signal yields for representative gluino signal model compressed points in the aggregate search regions. The uncertainties shown are statistical.

Simulated signal yields for representative squark signal model light-LSP points in the aggregate search regions. The uncertainties shown are statistical.

Simulated signal yields for representative squark signal model compressed points in the aggregate search regions. The uncertainties shown are statistical.

The 95% CL observed and expected upper limits on the production cross sections of the T1tttt signal model as a function of the gluino and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T1bbbb signal model as a function of the gluino and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T1qqqq signal model as a function of the gluino and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T5qqqqVV signal model as a function of the gluino and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T2tt signal model as a function of the top squark and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T2bb signal model as a function of the bottom squark and LSP masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

The 95% CL observed and expected upper limits on the production cross sections of the T2qq signal model as a function of the LSP and light squark masses.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (u d s c)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (single flavor)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (single flavor)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (single flavor)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (single flavor)

Exclusion limits assuming the approximate-NNLO+NNLL cross sections (single flavor)

The observed significance of the SMS models for T1tttt as function of the gluino and LSP masses.

The observed significance of the SMS models for T1bbbb as function of the gluino and LSP masses.

The observed significance of the SMS models for T1qqqq as function of the gluino and LSP masses.

The observed significance of the SMS models for T5qqqqVV as function of the gluino and LSP masses.

The observed significance of the SMS models for T2tt as function of the top squark and LSP masses.

The observed significance of the SMS models for T2bb as function of the bottom squark and LSP masses.

The observed significance of the SMS models for T2qq as function of the LSP and light squark masses.

The efficiency times acceptance of the SMS models for T1tttt as function of the gluino and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T1bbbb as function of the gluino and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T1qqqq as function of the gluino and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T5qqqqVV as function of the gluino and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T2tt as function of the top squark and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T2bb as function of the bottom squark and LSP masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

The efficiency times acceptance of the SMS models for T2qq as function of the LSP and light squark masses. The efficiency times acceptance is given for the total signal across the set of 174 search regions.

Absolute cumulative efficiencies in percent for each step of the event selection process, listed for four representative gluino pair production signal models with m(gluino) >> m(LSP). Only statistical uncertainties are shown.

Absolute cumulative efficiencies in percent for each step of the event selection process, listed for four representative gluino pair production signal models with m(gluino) ~ m(LSP). Only statistical uncertainties are shown.

Absolute cumulative efficiencies in percent for each step of the event selection process, listed for three representative squark pair production signal models with m(squark) >> m(LSP). Only statistical uncertainties are shown.

Absolute cumulative efficiencies in percent for each step of the event selection process, listed for three representative squark pair production signal models with m(squark) ~ m(LSP). Only statistical uncertainties are shown.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels in the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

The product of acceptance, efficiency, and branching fraction of the qbX signal with $\mathrm{X} \rightarrow \tau\tau$ and $m_{\mathrm{B}}=170~\mathrm{GeV}$ in $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels of the 1b1f and 1b1c event categories, as a function of the X boson mass. The selections are as described in the paper. The uncertainty refers to the statistical uncertainty only.

Observed $m_{\tau\tau}$ distribution in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1f category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1c category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mu\tau_\mathrm{h}$ channel of the 1b1f category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mu\tau_\mathrm{h}$ channel of the 1b1c category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 300 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 450 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1f category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1c category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1b1f category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1b1c category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Differential cross section dsigma(Z+c)/dpTZ times branching ratio and differential cross sections ratio dsigma(Z+c)/dpTZ / dsigma(Z+b)/dpTZ for three ranges of the transverse momentum of the Z boson and combining the dielectron and dimuon Z0 decay channels

Differential cross section dsigma(Z+c)/dpTjet times branching ratio and differential cross sections ratio dsigma(Z+c)/dpTZ / dsigma(Z+b)/dpTjet for three ranges of the transverse momentum of the jet and combining the dielectron and dimuon Z0 decay channels

Transverse momentum distribution of the selected muon-inside-a-jet for events with an identified muon among the jet constituents in the dielectron channel

Transverse momentum distribution of the selected muon-inside-a-jet for events with an identified muon among the jet constituents in the dimuon channel

The invariant mass distribution of three-prong secondary vertices for events selected in the D± mode, in the dielectron channel

The invariant mass distribution of three-prong secondary vertices for events selected in the D ± mode, in the dimuon channel

The invariant mass distribution of the three-track system composed of a two-prong secondary vertex and a primary particle for events selected in the D∗(2010)± mode, in the dielectron channel

The invariant mass distribution of the three-track system composed of a two-prong secondary vertex and a primary particle for events selected in the D ∗ (2010) ± mode, in the dimuon channel

Transverse momentum distribution of the c-tagged jet normalized to unity, in simulated W+c and Z+c samples and in W+c data events.

Number of reconstructed secondary vertices normalized to unity, in simulated W+c and Z+c samples and in W+c data events.

Distribution of the corrected secondary-vertex mass normalized to unity, in simulated W + c and Z + c samples, and in W + c data events.

Distribution of the JP discriminant for the D± mode normalized to unity, in simulated W + c and Z + c samples, and in W + c data events.

Distribution of the JP discriminant for the D∗(2010) ± mode normalized to unity, in simulated W + c and Z + c samples, and in W + c data events.

Distribution of the corrected secondary-vertex mass normalized to unity from simulated Z + b and eμ-tt data events

Corrected secondary-vertex mass distributions, after background subtraction in the dielectron channel for events selected in the semileptonic mode.

Corrected secondary-vertex mass distributions, after background subtraction in the dimuon channel for events selected in the semileptonic mode.

Background-subtracted distributions of the JP discriminant in the dielectron channel for Z + jets events with a D± → K∓ π± π± candidate.

Background-subtracted distributions of the JP discriminant in the dimuon channel for Z + jets events with a D± → K∓ π± π± candidate.

Background-subtracted distributions of the JP discriminant in the dielectron channel for Z +jets events with a D∗(2010)± → D0 π± → K∓ π± π± candidate.

Background-subtracted distributions of the JP discriminant in the dimuon channel for Z +jets events with a D∗(2010)± → D0 π± → K∓ π± π± candidate.

ontributions to the systematic uncertainty in the measured Z+c cross section and in the (Z+c)/(Z+b) cross section ratio.

Differential Z + c cross section as a function of the transverse momentum of the Z boson

(Z +c)/(Z +b) cross section ratio as a function of the transverse momentum of teh Z boson

Differential Z + c cross section as a function of the transverse momentum of the jet

(Z +c)/(Z +b) cross section ratio as a function of the transverse momentum of the jet.

Observed event yields and estimated backgrounds in the low- and high-mass control regions. The uncertainty in the background yield is the sum in quadrature of the statistical and systematic uncertainties.

Summary of the relative systematic uncertainties in heavy Majorana neutrino signal yields and the background from prompt same-sign leptons, both estimated from simulation. The relative systematic uncertainties assigned to the data-driven backgrounds for the misidentified lepton background and mismeasured charge background are also shown. The uncertainties are given for the low-mass selections.

Summary of the relative systematic uncertainties in heavy Majorana neutrino signal yields and the background from prompt same-sign leptons, both estimated from simulation. The relative systematic uncertainties assigned to the data-driven backgrounds for the misidentified lepton background and mismeasured charge background are also shown. The uncertainties are given for the high-mass selections.

Summary of contributions to the systematic uncertainty related to the prompt same-sign leptons background, misidentified lepton background, and mismeasured charge background on the total background uncertainty for the case of $m_{\rm N}$ = 100 and 500 GeV.

Observed event yields and estimated backgrounds after the application of all selection, except for the final optimization. The background predictions from prompt same-sign leptons (Prompt bkgd.), misidentified leptons (Misid. bkgd.), mismeasured charge (Charge mismeas. bkgd.) and the total background (Total bkgd.) are shown together withthe number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Dielectron channel results after final optimization. The background predictions from prompt same-sign leptons, misidentified leptons, and mismeasured charge are shown along with the total background estimate and the number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Electron-muon channel results after final optimization. The background predictions from prompt same-sign leptons and misidentified leptons are shown along with the total background estimate and the number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.