Measured diferential cross sections as a function of the transverse momentum of the lepton from the W decay.

Measured diferential cross sections ratio $R=\sigma(W^+ + \overline{c}) / \sigma(W^- + c)$ as a function of the absolute value of the pseudorapidity of the lepton from the W decay.

Measured diferential cross sections ratio $R=\sigma(W^+ + \overline{c}) / \sigma(W^- + c)$ as a function of the transverse momentum of the lepton from the W decay.

Distribution of the three leading leptons flavour in the SR-WZ with uncertainties evaluated after the inclusive cross section fit

Efficiency, acceptance, and proportion of events with leptonic tau decays in WZ production

WZ fiducial cross section in the four flavour exclusive and the flavour inclusive channels

WZ total cross section extrapolated from the four flavour exclusive and the flavour inclusive channels

Distribution of the total lepton charge in the SR-WZ with uncertainties evaluated after the inclusive cross section fit

W$^{+}$Z fiducial cross section in the four flavour exclusive and the flavour inclusive channels

W$^{-}$Z fiducial cross section in the four flavour exclusive and the flavour inclusive channels

WZ charge asymmetry ratio measured on each of the four flavour exclusive and the flavour inclusive channels

Distribution of the cosine of the W polarization angle times total lepton charge in the SR-WZ with uncertainties evaluated after the W polarization fit

Distribution of the cosine of the Z polarization angle in the SR-WZ with uncertainties evaluated after the Z polarization fit

Best fits to the W and Z polarization fractions

2D confidence regions at the 68, 95, and 99% CL in the $f_O^W$-$f_{L}^W-f_R^W$ plane

2D confidence regions at the 68, 95, and 99% CL in the $f_O^Z$-$f_{L}^Z-f_R^Z$ plane

Distribution of the invariant mass of the WZ system in the SR-WZ with uncertainties evaluated after the inclusive cross section fit

Best fit values and one dimensional confidence regions in several EFT coefficients obtained from the EFT fit considering both the SM interferences and purely BSM (order $\Lambda^{-2}$ and $\Lambda^{-4}$) terms

2D confidence regions at the 68, 95, and 99% CL in the $c_{www}$-$c_{w}$ plane

2D confidence regions at the 68, 95, and 99% CL in the $c_{w}$-$c_{b}$ plane

2D confidence regions at the 68, 95, and 99% CL in the $c_{www}$-$c_{w}$ plane

Best fit values and one dimensional confidence regions in several EFT coefficients obtained from the EFT fit considering only the SM-EFT interference (order $\Lambda^{-2}$) terms

Evolution of the best fit and expected and observed 95% CI for the $c_{w}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $c_{b}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $c_{www}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $\tilde{c}_{www}$ parameter as a function of the cutoff scale

Evolution of the best fit and expected and observed 95% CI for the $\tilde{c}_{w}$ parameter as a function of the cutoff scale

Differential cross section with respect to the transverse momentum of the Z boson

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the Z boson

Response matrix for the $p_{T}$ of the Z boson obtained with POWHEG

Differential cross section with respect to the transverse momentum of the leading jet

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the leading jet

Response matrix for the $p_{T}$ of the leading jet obtained with POWHEG

Differential cross section with respect to the jet multiplicity

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the jet multiplicity

Response matrix for the jet multiplicity obtained with POWHEG

Differential cross section with respect to the invariant mass of the WZ system

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the invariant mass of the WZ system

Response matrix for the invariant mass of the WZ system obtained with POWHEG

Differential cross section with respect to the transverse momentum of the lepton associated to the W boson

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the lepton associated to the W boson

Response matrix for the $p_{T}$ of the lepton associated to the W boson obtained with POWHEG

Differential cross section with respect to the transverse momentum of the lepton associated to the W boson, W$^{+}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the lepton associated to the W boson, W$^{+}$Z only

Differential cross section with respect to the transverse momentum of the lepton associated to the W boson, W$^{-}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the $p_{T}$ of the lepton associated to the W boson, W$^{-}$Z only

Differential cross section with respect to the cosine of the W polarization angle times total lepton charge

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the W polarization angle times total lepton charge

Response matrix for the cosine of the W polarization angle times total lepton charge obtained with POWHEG

Differential cross section with respect to the cosine of the W polarization angle times total lepton charge, W$^{+}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the W polarization angle times total lepton charge, W$^{+}$Z only

Differential cross section with respect to the cosine of the W polarization angle times total lepton charge, W$^{-}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the W polarization angle times total lepton charge, W$^{-}$Z only

Differential cross section with respect to the cosine of the Z polarization angle

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the Z polarization angle

Response matrix for the cosine of the Z polarization angle obtained with POWHEG

Differential cross section with respect to the cosine of the Z polarization angle, W$^{+}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the Z polarization angle, W$^{+}$Z only

Differential cross section with respect to the cosine of the Z polarization angle, W$^{-}$Z only

Correlation matrix for the unfolded results obtained using NNLO bias, area con-straint, and no additional regularization for the cosine of the Z polarization angle, W$^{-}$Z only

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The efficiency of $N_{station} > 1$ requirement in bins of hadronic energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bin corresponds to LLPs that decayed leptonically with 0 hadronic energy. The efficiency is calculated with respect to clusters that pass all cluster-level cuts described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 10%.

The geometric signal acceptance as a function of $c\tau$. The acceptance is shown for four different LLP mass hypotheses: 7, 15, 40, and 55 GeV. The acceptance is defined by requiring at least 1 LLP to decay in the region defined as 400 cm $< |z| <$ 1100 cm, $r < $ 695.5 cm, and $|\eta| <$ 2.4.

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.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Differential signal yields for various signal hypotheses.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Background and data yields in the control and signal region bins. The prediction before ("prefit") and after the background only fit ("b-only") are given separately.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit to the control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit to the corresponding control regions, and serve as input to the simplified likelihood reinterpretation scheme.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Unweighted signal acceptance times efficiency at every cut stage. The requirements called "HCAL mitigation" refer to the requirements imposed in the 2018 data set in order to mitigate the localized failure of the HCAL detector.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with axial couplings.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the simplified model with vector couplings.

Upper limits on the coupling $g_{\chi}$ in the simplified model with a axial mediator.

Upper limits on the coupling $g_{q}$ in the simplified model with a axial mediator.

Upper limits on the coupling $g_{\chi}$ in the simplified model with a vector mediator.

Upper limits on the coupling $g_{q}$ in the simplified model with a vector mediator.

Exclusion limits on the signal strength in the simplified model with scalar couplings.

Exclusion limits on the signal strength in the simplified model with pseudoscalar couplings.

Exclusion limits on the fundamental Planck scale $M_{D}$ as a function of the number of extra dimensions $d$.

Median Expected exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Observed exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Expected plus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Expected minus 1 s.d. exclusion contour in the $m_{med}$-$m_{\chi}$ plane in the fermion portal model.

Tagging efficiency for AK8 jets. The efficiency includes the effect of the machine-learning based DeepAK8 tagger, as well as the application of the mass window requirement on the jet. The efficiency is split depending on the matching generator-level object. For reinterpretation purposes, this efficiency can directly be applied to any AK8 jet that is matched to a given type of generator-level object, with no prior selection on the jet mass. Other acceptance requirements on jet $\p_{T}$ and $\eta$ should still be applied.

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