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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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

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_{\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_{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_{\chi}$ in the simplified model with a vector mediator.

Upper limits on the coupling $g_{q}$ 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 scalar couplings.

Exclusion limits on the signal strength in the simplified model with pseudoscalar 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$.

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.

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.

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

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.

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.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 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 40 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 40 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 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 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 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 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 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 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 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 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 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 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 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 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.

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.

The timing distribution of the background sources predicted to contribute to the signal region, compared to those for a representative signal model. The time is defined by the jet in the event with the largest $t_{\mathrm{jet}}$ passing the relevant selection. The distributions for the major backgrounds are taken from control regions and normalized to the predictions. The observed data is shown by the black points. No events are observed in data for $t_{\mathrm{jet}} > 3\,$ns (indicated with a vertical black line).

The product of the acceptance and efficiency in the $c\tau_{0}$ vs. $m_{\tilde{g}}$ plane for the GMSB model, after all requirements.

The product of the acceptance and efficiency in the $c\tau_{0}$ vs. $m_{\tilde{g}}$ plane for the GMSB model, after all requirements.

The observed upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed.

The observed upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The observed upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed.

The observed upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The observed and expected upper limits at $95\%$ CL on the gluino pair production cross section for a gluino GMSB model with $m_{\tilde{g}}=2400$ GeV. The one (two) standard deviation variation in the expected limit is shown in the inner green (outer yellow) band. The blue solid line shows the observed limit obtained by the CMS displaced jet search

The observed and expected upper limits at $95\%$ CL on the gluino pair production cross section for a gluino GMSB model with $m_{\tilde{g}}=2400$ GeV. The one (two) standard deviation variation in the expected limit is shown in the inner green (outer yellow) band. The blue solid line shows the observed limit obtained by the CMS displaced jet search

The distribution (normalized to unity) of number of ECAL cells hit in the jet for jets in a background enriched data sample (satisfying $|\eta| < 1.48$, $PV_{\rm track}^{\rm fraction} > 1/12$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets satisfying signal region requirements (except those on $E_{\mathrm{ECAL}}$ and $N^{\mathrm{cell}}_{\mathrm{ECAL}}$).

The distribution (normalized to unity) of number of ECAL cells hit in the jet for jets in a background enriched data sample (satisfying $|\eta| < 1.48$, $p_\mathrm{T}tf > 0.08$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets satisfying signal region requirements (except those on $E_{\mathrm{ECAL}}$ and $N^{\mathrm{cell}}_{\mathrm{ECAL}}$).

The distribution (normalized to unity) of $\mathrm{HEF}$ for a data sample enriched in beam halo and noise jets (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30 \,\mathrm{GeV}$, $PV_{\rm track}^{\rm fraction} < 1/12$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$) and for signal jets passing signal region selections (except on $\mathrm{HEF}$).

The distribution (normalized to unity) of $\mathrm{HEF}$ for a data sample enriched in beam halo and noise jets (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30 \,\mathrm{GeV}$, $p_\mathrm{T}tf < 0.08$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$) and for signal jets passing signal region selections (except on $\mathrm{HEF}$).

The distribution (normalized to unity) of $E_{\mathrm{HCAL}}$ for a data sample enriched in beam halo and noise jets (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30 \,\mathrm{GeV}$, $PV_{\rm track}^{\rm fraction} < 1/12$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$) and for signal jets passing signal region selections (except on $E_{\mathrm{HCAL}}$).

The distribution (normalized to unity) of $E_{\mathrm{HCAL}}$ for a data sample enriched in beam halo and noise jets (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30 \,\mathrm{GeV}$, $p_\mathrm{T}tf < 0.08$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$) and for signal jets passing signal region selections (except on $E_{\mathrm{HCAL}}$).

The distribution (normalized to unity) of $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ for data sample enriched in jets from noise (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $PV_{\rm track}^{\rm fraction} > 1/12$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20 \,\mathrm{GeV}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets passing signal region selections (except on $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ and $t^{\mathrm{RMS}}_\mathrm{jet}$)

The distribution (normalized to unity) of $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ for data sample enriched in jets from noise (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $p_\mathrm{T}tf > 0.08$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20 \,\mathrm{GeV}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets passing signal region selections (except on $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ and $t^{\mathrm{RMS}}_\mathrm{jet}$)

The distribution (normalized to unity) of $t^{\mathrm{RMS}}_\mathrm{jet}$ for data sample enriched in jets from noise (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $PV_{\rm track}^{\rm fraction} > 1/12$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20 \,\mathrm{GeV}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets passing signal region selections (except on $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ and $t^{\mathrm{RMS}}_\mathrm{jet}$)

The distribution (normalized to unity) of $t^{\mathrm{RMS}}_\mathrm{jet}$ for data sample enriched in jets from noise (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $p_\mathrm{T}tf > 0.08$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20 \,\mathrm{GeV}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets passing signal region selections (except on $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ and $t^{\mathrm{RMS}}_\mathrm{jet}$)

The distribution (normalized to unity) of $PV_{\rm track}^{\rm fraction}$ for a data sample enriched in main bunch backgrounds (satisfying $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $|t_{\mathrm{jet}}| < 3\,\mathrm{ns}$ and $E_{\mathrm{ECAL}} >20\,\mathrm{GeV}$) and for signal jets passing signal selections (except on $PV_{\rm track}^{\rm fraction}$).

The distribution (normalized to unity) of $p_\mathrm{T}tf$ for a data sample enriched in main bunch backgrounds (satisfying $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $|t_{\mathrm{jet}}| < 3\,\mathrm{ns}$ and $E_{\mathrm{ECAL}} >20\,\mathrm{GeV}$) and for signal jets passing signal selections (except on $p_\mathrm{T}tf$).

The distribution (normalized to unity) of $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}}$ for a data sample enriched in beam halo (satisfying $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $PV_{\rm track}^{\rm fraction} < 1/12$, $\mathrm{HEF} < 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$ and $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$) and for signal jets passing signal selections (except on $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}}$)

The distribution (normalized to unity) of $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}}$ for a data sample enriched in beam halo (satisfying $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $p_\mathrm{T}tf < 0.08$, $\mathrm{HEF} < 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$ and $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$) and for signal jets passing signal selections (except on $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}}$)

The distribution (normalized to unity) of $\mathrm{max}(\Delta \phi_{\mathrm{DT}})$ for a data sample enriched in cosmic muons (satisfying $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $PV_{\rm track}^{\rm fraction} < 1/12$, $\mathrm{HEF} > 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} > 3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and failing the HCAL noise rejection quality filters) and for signal jets passing signal selections (except on $\mathrm{max}(\Delta \phi_{\mathrm{DT}})$).

The distribution (normalized to unity) of $\mathrm{max}(\Delta \phi_{\mathrm{DT}})$ for a data sample enriched in cosmic muons (satisfying $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $p_\mathrm{T}tf < 0.08$, $\mathrm{HEF} > 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} > 3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and failing the HCAL noise rejection quality filters) and for signal jets passing signal selections (except on $\mathrm{max}(\Delta \phi_{\mathrm{DT}})$).

The distribution (normalized to unity) of $\mathrm{max}(\Delta \phi_{\mathrm{RPC}})$ for a data sample enriched in cosmic muons (satisfying $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $PV_{\rm track}^{\rm fraction} < 1/12$, $\mathrm{HEF} > 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} > 3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and failing the HCAL noise rejection quality filters) and for signal jets passing signal selections (except on $\mathrm{max}(\Delta \phi_{\mathrm{RPC}})$).

The distribution (normalized to unity) of $\mathrm{max}(\Delta \phi_{\mathrm{RPC}})$ for a data sample enriched in cosmic muons (satisfying $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $p_\mathrm{T}tf < 0.08$, $\mathrm{HEF} > 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} > 3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and failing the HCAL noise rejection quality filters) and for signal jets passing signal selections (except on $\mathrm{max}(\Delta \phi_{\mathrm{RPC}})$).

Distribution of $t_{\mathrm{jet}}$ for jets with the full Run 2 dataset with no cleaning selection applied (a) and after all jet cleaning selections are applied (b) in events satisfying the trigger requirements and satisfying $p_T^\mathrm{miss} > 300$. The jets are required to pass an inverted selection of $PV_{\rm track}^{\rm fraction} > 1/12$ to enrich the sample in those originating from main and satellite bunch backgrounds. The cleaning selections are shown to reduce the backgrounds by many orders of magnitude.

Distribution of $t_{\mathrm{jet}}$ for jets with the full Run 2 dataset with no cleaning selection applied (a) and after all jet cleaning selections are applied (b) in events satisfying the trigger requirements and satisfying $p_T^\mathrm{miss} > 300$. The jets are required to pass an inverted selection of $PV_{\rm track}^{\rm fraction} > 0.08$ to enrich the sample in those originating from main and satellite bunch backgrounds. The cleaning selections are shown to reduce the backgrounds by many orders of magnitude.

Distribution of $t_{\mathrm{jet}}$ for jets with the full Run 2 dataset in events satisfying the trigger requirements and satisfying $p_T^\mathrm{miss} < 300$. The jets are required to pass an inverted selection of $PV_{\rm track}^{\rm fraction} > 1/12$ to enrich the sample in those originating from main and satellite bunch backgrounds (all other jet cleaning selections are applied). Clear contributions from jets from satellite bunch collisions can be seen peaked around -5, 5 and 10 ns.

Distribution of $t_{\mathrm{jet}}$ for jets with the full Run 2 dataset in events satisfying the trigger requirements and satisfying $p_T^\mathrm{miss} < 300$. The jets are required to pass an inverted selection of $PV_{\rm track}^{\rm fraction} > 0.08$ to enrich the sample in those originating from main and satellite bunch backgrounds (all other jet cleaning selections are applied). Clear contributions from jets from satellite bunch collisions can be seen peaked around -5, 5 and 10 ns.

Contribution to the delayed time of jets from the $\beta$ of the gluino is plotted against the delay contribution from the difference (assuming straight line paths) between the path taken by the gluino and the particle forming the jet from the path length for a particle travelling directly to the same position on the ECAL barrel for gluino $c\tau_{0} = 10$ m and mass of 1000 Gev (top) and 3000 GeV (bottom). The dominant contribution is shown to be the gluino $\beta$.

Contribution to the delayed time of jets from the $\beta$ of the gluino is plotted against the delay contribution from the difference (assuming straight line paths) between the path taken by the gluino and the particle forming the jet from the path length for a particle travelling directly to the same position on the ECAL barrel for gluino $c\tau_{0} = 10$ m and mass of 1000 Gev (top) and 3000 GeV (bottom). The dominant contribution is shown to be the gluino $\beta$.

Selection efficiencies for the GMSB model with $m_{\tilde{g}}=1000$ and various proper decay lengths

Selection efficiencies for the GMSB model with $m_{\tilde{g}}=1000$ and various proper decay lengths

Selection efficiencies for the GMSB model with $m_{\tilde{g}}=2400$ and various proper decay lengths

Selection efficiencies for the GMSB model with $m_{\tilde{g}}=2400$ and various proper decay lengths

Selection efficiencies for the GMSB model with $m_{\tilde{g}}=3000$ and various proper decay lengths

Selection efficiencies for the GMSB model with $m_{\tilde{g}}=3000$ and various proper decay lengths

The observed upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The observed upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected plus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected plus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected minus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected minus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected plus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected plus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected minus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected minus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

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.

Measured misidentification probability distribution as a function of track multiplicity for the EMJ-1 criteria group defined in Table 2. The red up-pointing triangles are for b jets while the blue down-pointing triangles are for light-flavor jets. The horizontal lines on the data points indicate the variable bin width.

The $H_\mathrm{T}$ distribution for the observed data events (black points) and the predicted background estimation (blue) for selection set 8 (SM QCD-enhanced), requiring at least two jets tagged by loose emerging jet criteria.

The number of associated tracks distribution for the observed data events (black points) and the predicted background estimation (blue) for selection set 8 (SM QCD-enhanced), requiring at least two jets tagged by loose emerging jet criteria.

The $H_\mathrm{T}$ distribution for the observed data events (black points) and the predicted background estimation (blue) for selection set 9 (SM QCD-enhanced), requiring at least one jet tagged by loose emerging jet criteria and large $p_\mathrm{T}^\mathrm{miss}$.

The number of associated tracks distribution for the observed data events (black points) and the predicted background estimation (blue) for selection set 9 (SM QCD-enhanced), requiring at least one jet tagged by loose emerging jet criteria and large $p_\mathrm{T}^\mathrm{miss}$.

Upper limits at 95% CL on the signal cross section for models with dark pion mass $(m_{\pi_\mathrm{DK}})$ of 5 GeV in the proper decay length $(c\tau_{\pi_\mathrm{DK}})$ versus dark mediator mass $(m_{\mathrm{X_{DK}}})$ plane.

Dark mediator mass exclusion limits at 95% CL derived from theoretical cross sections for models with dark pion mass $(m_{\pi_\mathrm{DK}})$ of 5 GeV in the proper decay length $(c\tau_{\pi_\mathrm{DK}})$ versus dark mediator mass $(m_{\mathrm{X_{DK}}})$ plane.