The $m_{e\mu}$ distribution of the observed data is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in red solid line, and the background component of the fit in red dotted line. Events and fit in each category are weighted by $S/(S+B)$. The one and two standard deviation uncertainty bands of the background component are shown in green and yellow. The lower panel shows the residuals after subtracting the background component of the fit from data.

Observed (expected) 95% confidence level upper limits on $\sigma(p p \to X(146) \to e \mu)$ for each individual analysis category (as shown in the left axis label) and for the combination of all analysis categories assuming the relative SM-like production cross sections of the ggH and VBF production modes.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 0 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 1 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 2 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 3 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category VBF cat 0 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category VBF cat 1 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

Best fit of $\sigma(p p \to X(146) \to e \mu)$ for each analysis category (black point) compared to the best fit for the combination of all analysis categories (blue line). The one standard deviation uncertainty of the per-category best fit is shown in red line and the one standard deviation uncertainty of the combined best fit is shown in green band.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.0.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.1.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.2.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.3.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.4.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.5.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.6.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.7.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.8.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.9.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 1.0.

Observed and expected 95% CL upper limits on ${{(\sigma_{\text{H}}/\sigma_{\text{H}}^{\mathrm{SM}}) \times {{\mathcal{B}(\text{H} \to \text{inv})}}}}$ for the combined VBF, ttH, VH and ggH channels using all available CMS data for 2016--2018, Run1, 2015, and their combination, assuming a SM Higgs boson with a mass of 125 GeV.

Observed and expected 95% CL upper limits on ${{(\sigma_{\text{H}}/\sigma_{\text{H}}^{\mathrm{SM}}) \times {{\mathcal{B}(\text{H} \to \text{inv})}}}}$ for the VBF, ttH(lep), ttH(had), V(qq)H, Z(ll)H, Z(bb)H and ggH channels using all available CMS data, and their combination, assuming a SM Higgs boson with a mass of 125 GeV.

The 90% CL upper limits on the spin-independent DM-nucleon scattering cross section in Higgs-portal models, assuming a scalar, fermionic or vectorial DM candidate (with a UV-complete vector model), compared to direct detection experimets CRESST-III, DarkSide-50, LUX_ZEPLIN and PandaX-4T (2022 results).

Postfit event yields in each category and hadronic recoil bin of the ttH category SR for the different background and signal processes summed across 2016, 2017 and 2018 samples. The background and signal yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only. The Postfit signal contributions for a Higgs boson, produced in ttH, VH (split into WH and ZH), ggH and VBF, decaying to invisible particles with a branching fraction of $\mathcal{B}(\text{H} \to \text{inv}) = 1$, and the observed event yields are reported.

Postfit event yields in each category and hadronic recoil bin of the VH category SR for the different background and signal processes summed across 2016, 2017 and 2018 samples. The background and signal yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only. The Postfit signal contributions for a Higgs boson, produced in ttH, VH (split into WH and ZH), ggH and VBF, decaying to invisible particles with a branching fraction of $\mathcal{B}(\text{H} \to \text{inv}) = 1$, and the observed event yields are reported.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the ttH category e$+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the VH category e$+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the ttH category $\mu+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the VH category $\mu+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the ttH category ee$+$jets, $\mu\mu+$jets, and $\ell\ell+$jets CRs for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the VH category ee$+$jets, $\mu\mu+$jets, and $\gamma+$jets CRs for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Feynman diagrams of the effective neutralino pair production in the GMSB simplified model in which the two neutralinos decay into two gravitinos ($\tilde{G}$) and two $Z$ bosons (left), a $Z$ and a Higgs boson ($H$) (center), or two Higgs bosons (right).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

TD jet tagger score distributions (left) for signal (red, lighter) and collision background (blue, darker). Identification probability for the signal versus the misidentification probability for the background (right) with the tagger working point (w.~p.) used in the analysis shown as a blue marker. Both are evaluated using an independent sample of testing events.

TD jet tagger score distributions (left) for signal (red, lighter) and collision background (blue, darker). Identification probability for the signal versus the misidentification probability for the background (right) with the tagger working point (w.~p.) used in the analysis shown as a blue marker. Both are evaluated using an independent sample of testing events.

The efficiency of the TD jet tagger working point used in the analysis is shown as a function of the lab frame transverse decay length for simulated signals with $\chi$ mass of 400 GeV. The uncertainties shown account for lifetime dependence and statistical uncertainty.

The TD jet tagger score distributions for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. The last bin contains jets with tagger scores greater than 0.996, the threshold used to tag signal jets. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The TD jet tagger misidentification probability measured using the nominal $W$+jets MR (black round markers) is shown along with the systematic uncertainty (gray band), quantifying the degree of process dependence measured from alternative MRs. The measurements in the alternative MRs are displayed as well ($Z$+jets MR as green round markers, $t\bar{t}$ MR as red squared markers, QCD MR as blue triangular markers) along with their respective statistical uncertainty. On the left, this probability is shown for the first 19.9 fb$^{-1}$ of data collected in 2016, while on the right it is shown for the last 16.4 fb$^{-1}$ of data collected in 2016combined with data collected in 2017-2018.

The TD jet tagger misidentification probability measured using the nominal $W$+jets MR (black round markers) is shown along with the systematic uncertainty (gray band), quantifying the degree of process dependence measured from alternative MRs. The measurements in the alternative MRs are displayed as well ($Z$+jets MR as green round markers, $t\bar{t}$ MR as red squared markers, QCD MR as blue triangular markers) along with their respective statistical uncertainty. On the left, this probability is shown for the first 19.9 fb$^{-1}$ of data collected in 2016, while on the right it is shown for the last 16.4 fb$^{-1}$ of data collected in 2016combined with data collected in 2017-2018.

Distribution of the number of TD tagged jets for the $m_{\chi} = 400$ GeVsimulated signal samples with $c\tau_{\chi} = 0.5$ m (solid red line) and $c\tau_{\chi} = 3.0$ m (dotted green line), estimated background (blue square markers), and data (black round markers). The signal distributions are normalized to the expected cross section limit. The blue shaded region indicates the systematic uncertainty in the background prediction. No background prediction is shown for the bin with zero TD tagged jets as it is the main control region used to predict the background for the other two bins. There are zero observed events in the bin with two or more TD tagged jets.

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $m_\chi$ in a scenario with $\mathcal{B}(\chi\to HG) = 0.5$ and $c\tau = 0.5$ m (left) or 3 m (right).

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $m_\chi$ in a scenario with $\mathcal{B}(\chi\to HG) = 0.5$ and $c\tau = 0.5$ m (left) or 3 m (right).

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.5$ and $m_{\chi} = 400$ GeV (left) or 1000 GeV (right).

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.5$ and $m_{\chi} = 400$ GeV (left) or 1000 GeV (right).

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.5$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The distribution of the jet charged hadron energy fraction, a variable used as input to the TD jet tagger score, for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The distribution of the jet neutral hadron energy fraction, a variable used as input to the TD jet tagger score, for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The distribution of the number of track constituents in the jet, a variable used as input to the TD jet tagger score, for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The $\eta$ distribution of TD-tagged jets in a background-enriched control region that comprises events with exactly one TD-tagged jet. Observed data (black round markers) and the corresponding prediction based on control samples in data (empty squared markers), measured using the nominal $W$+jets MR, are compared. The prediction uncertainty (gray band) includes the systematic uncertainty quantifying the degree of process dependence measured from alternative MRs. The predictions for the shape and the normalization of the $\eta$ distribution are consistent with the data.

Jet time distribution in a sample of b-tagged jets from dilepton $t \bar{t}$ events in 2017 data-taking period (black round markers) and simulation (filled histogram). A Gaussian smearing procedure is applied to the jet time in the $t \bar{t}$ sample (green line) to correct for effects that are difficult to simulate (timing drift caused by crystals transparency loss due to detector aging, electronics jitter).

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 1$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.75$, $\mathcal{B}(\chi\to Z\tilde{G}) = 0.25$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.25$, $\mathcal{B}(\chi\to Z\tilde{G}) = 0.75$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to Z\tilde{G}) = 1$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 127 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 127 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 127 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 150 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 150 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 150 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 175 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 175 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 175 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 200 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 200 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 200 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 300 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 300 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 300 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 400 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 400 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 400 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 600 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 600 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 600 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1000 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1000 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1000 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1500 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1500 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1500 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 127 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 127 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 150 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 150 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 175 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 175 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 200 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 200 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 250 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 250 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 300 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 300 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 400 GeV.

Expected and observed upper limits on the AQGC operators $a^Z_C/\Lambda^2$, with no unitarization. The $y$ axis shows the limit on the ratio of the observed cross section to the cross section predicted for each anomalous coupling value ($\sigma_\mathrm{AQGC}$).

Limits on LEP-like dimension-6 anomalous quartic gauge coupling parameters, with and without unitarization via a clipping procedure.

Conversion of limits on $a^W_0$ to dimension-8 $f_{M,i}$ operators, using the assumption of vanishing $WWZ\gamma$ couplings to eliminate some parameters. When quoting limits on one of the operators, the other is fixed to zero. The results for $|f_{M,0}/\Lambda^{4}|$ and $|f_{M,4}/\Lambda^{4}|$ are shown with and without clipping of the signal model at 1.4 TeV, when the other parameter is fixed to the SM value of zero.

Conversion of limits on $a^W_0$ and $a^W_C$ to dimension-8 $f_{M,i}$ operators, using the assumption that all $f_{M,i}$ except one are equal to zero. The results are shown with and without clipping of the signal model at 1.4 TeV.

{Expected and observed limits in the two-dimensional plane of $a^W_0/\Lambda^2$ vs. $a^W_C/\Lambda^2$. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^W$ and $a_C^W$ couplings, respectively.

{Expected and observed limits in the two-dimensional plane of $a^W_0/\Lambda^2$ vs. $a^W_C/\Lambda^2$ with unitarization imposed by clipping the signal model at 1.4 TeV. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^W$ and $a_C^W$ couplings, respectively.

{Expected and observed limits in the two-dimensional plane of $a^Z_0/\Lambda^2$ vs. $a^Z_C/\Lambda^2$. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^Z$ and $a_C^Z$ couplings, respectively.

Distribution of the MLP W+jets score in the CRs for the $T\overline{T}$ search. The observed data, predicted $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario, and the background are all shown. Statistical and systematic uncertainties in the background prediction before performing the fit to data are also shown.

Distribution of DeepAK8 jet tags in the DeepAK8 CR of the $T\overline{T}$ search. The observed data are shown using black markers, $T\overline{T}$ signals using solid lines, and background using filled histograms. Statistical and systematic uncertainties in the background prediction before performing the fit to data are shown by the hatched region.

Distribution of HT in the W+jets CR of the $T\overline{T}$ search. The observed data are shown using black markers, $T\overline{T}$ signals using solid lines, and background using filled histograms. Statistical and systematic uncertainties in the background prediction before performing the fit to data are shown by the hatched region.

Distribution of HT in the $T\overline{T}$ CR of the $T\overline{T}$ search. The observed data are shown using black markers, $T\overline{T}$ signals using solid lines, and background using filled histograms. Statistical and systematic uncertainties in the background prediction before performing the fit to data are shown by the hatched region.

Numbers of predicted and observed CR events in 2016--2018 data for the $T\overline{T}$ signal hypotheses in the single-lepton channel, after a background-only fit to data. Predicted numbers of signal events before the fit to data are included for comparison, using the singlet branching fraction scenario. Uncertainties include statistical and systematic components, and lepton flavor categories are combined with their uncertainties added in quadrature. Values in the DeepAK8 CR represent the number of large-radius jets rather than number of events.

Numbers of predicted and observed CR events in 2016--2018 data for the $B\overline{B}$ signal hypotheses in the single-lepton channel, after a background-only fit to data. Predicted numbers of signal events before the fit to data are included for comparison, using the singlet branching fraction scenario. Uncertainties include statistical and systematic components, and lepton flavor categories are combined with their uncertainties added in quadrature. Values in the DeepAK8 CR represent the number of large-radius jets rather than number of events.

Numbers of predicted and observed CR events in 2017--2018 data in the same-sign dilepton channel, after a background-only fit to data. Predicted numbers of signal events before the fit to data are included for comparison, using the singlet branching fraction scenario. Uncertainties include statistical and systematic components. Predictions for 2017 and 2018 are combined with their uncertainties added in quadrature.

Distribution of lepton $p_T$ in 2017 in the multilepton channel nonprompt lepton control region for all flavor categories, evaluated with the best-fit nonprompt rates. The uncertainty shown is the quadratic sum of statistical and systematics uncertainties.

Distribution of lepton $p_T$ in 2018 in the multilepton channel nonprompt lepton control region for all flavor categories, evaluated with the best-fit nonprompt rates. The uncertainty shown is the quadratic sum of statistical and systematics uncertainties.

Numbers of predicted and observed CR events in 2017--2018 data in the multilepton channel, after a background-only fit to data. Predicted numbers of signal events before the fit to data are included for comparison, using the singlet branching fraction scenario. Uncertainties include statistical and systematic components. Predictions for 2017 and 2018 are combined with their uncertainties added in quadrature.

Numbers of predicted and observed $T\overline{T}$ SR events in 2016--2018 data in the single-lepton channel, after a background-only fit to data. Predicted numbers of signal events before the fit to data are included for comparison, using the singlet branching fraction scenario. Uncertainties include statistical and systematic components, and lepton flavor categories are combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 1. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 2. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 3. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 4. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 5. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 6. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 7. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 8. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 9. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 10. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 11. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of the VLQ score in single-lepton category 12. Observed data is shown using black markers, $T\overline{T}$ signal with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region. Electron and muon categories have been combined with their uncertainties added in quadrature.

Distribution of $H_T^{\mathrm{lep}}$ in the same-sign dilepton signal region for $ee$ category. Data from 2017 and 2018 is shown using black markers, $T\overline{T}$ signals with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region.

Distribution of $H_T^{\mathrm{lep}}$ in the same-sign dilepton signal region for $e\mu$ category. Data from 2017 and 2018 is shown using black markers, $T\overline{T}$ signals with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region.

Distribution of $H_T^{\mathrm{lep}}$ in the same-sign dilepton signal region for $\mu\mu$ category. Data from 2017 and 2018 is shown using black markers, $T\overline{T}$ signals with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region.

Distribution of $S_T$ in the multilepton signal region for $eee$ category. Data from 2017 and 2018 is shown using black markers, $T\overline{T}$ signals with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region.

Distribution of $S_T$ in the multilepton signal region for $ee\mu$ category. Data from 2017 and 2018 is shown using black markers, $T\overline{T}$ signals with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region.

Distribution of $S_T$ in the multilepton signal region for $e\mu\mu$ category. Data from 2017 and 2018 is shown using black markers, $T\overline{T}$ signals with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region.

Distribution of $S_T$ in the multilepton signal region for $\mu\mu\mu$ category. Data from 2017 and 2018 is shown using black markers, $T\overline{T}$ signals with mass of 1.2 (1.5) TeV in the singlet scenario using solid (dashed) lines, and post-fit background using filled histograms. Statistical and systematic uncertainties in the background prediction after performing the fit to data are shown by the hatched region.

Numbers of predicted and observed $B\overline{B}$ SR events in 2016--2018 data in the single-lepton channel, after a background-only fit to data. Predicted numbers of signal events before the fit to data are included for comparison, using the singlet branching fraction scenario. Uncertainties include statistical and systematic components, and lepton flavor categories are combined with their uncertainties added in quadrature.

Numbers of predicted and observed SR events in 2017--2018 data in the same-sign dilepton channel, after a background-only fit to data. Predicted numbers of signal events before the fit to data are included for comparison, using the singlet branching fraction scenario. Uncertainties include statistical and systematic components. Predictions for 2017 and 2018 are combined with their uncertainties added in quadrature.

Numbers of predicted and observed SR events in 2017--2018 data in the multilepton channel, after a background-only fit to data. Predicted numbers of signal events before the fit to data are included for comparison, using the singlet branching fraction scenario. Uncertainties include statistical and systematic components. Predictions for 2017 and 2018 are combined with their uncertainties added in quadrature.

Expected and observed lower limits on the VLT quark masses.

Expected and observed lower limits on the VLB quark masses.

Expected 95% CL upper limits on $T\overline{T}$ production cross sections for the singlet combination.

Expected 95% CL upper limits on $T\overline{T}$ production cross sections for the doublet combination.

Expected 95% CL upper limits on $B\overline{B}$ production cross sections for the singlet combination.

Expected 95% CL upper limits on $B\overline{B}$ production cross sections for the doublet combination.

The 95% CL expected lower mass limits on pair-produced T quark masses as functions of their branching ratios to W and H bosons.

The 95% CL observed lower mass limits on pair-produced T quark masses as functions of their branching ratios to W and H bosons.

The 95% CL expected lower mass limits on pair-produced B quark masses as functions of their branching ratios to W and H bosons.

The 95% CL observed lower mass limits on pair-produced B quark masses as functions of their branching ratios to W and H bosons.

The unfolded differential charge asymmetry as a function of the transverse momentum of the top pair system. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NNLO in QCD and NLO in EW theory are listed. The scale uncertainty is obtained by varying renormalisation and factorisation scales independently by a factor of 2 or 0.5 around $\mu_0$ to calculate the maximum and minimum value of the asymmetry, respectively. The nominal value $\mu_0$ is chosen as $H_T/4$. The variations in which one scale is multiplied by 2 while the other scale is divided by 2 are excluded. Finally, the scale and MC integration uncertainties are added in quadrature.

The unfolded differential charge asymmetry as a function of the longitudinal boost of the top pair system. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NNLO in QCD and NLO in EW theory are listed. The scale uncertainty is obtained by varying renormalisation and factorisation scales independently by a factor of 2 or 0.5 around $\mu_0$ to calculate the maximum and minimum value of the asymmetry, respectively. The nominal value $\mu_0$ is chosen as $H_T/4$. The variations in which one scale is multiplied by 2 while the other scale is divided by 2 are excluded. Finally, the scale and MC integration uncertainties are added in quadrature.

The unfolded inclusive leptonic asymmetry. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the invariant mass of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the transverse momentum of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the longitudinal boost of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

Individual 68% and 95% CL bounds on the relevant Wilson coefficients of the SM Effective Field Theory in units of $\text{TeV}^{-2}$. The bounds are derived from the $A_C^{t\bar{t}}$ inclusive measurement. The experimental uncertainties are accounted for, in the form of the complete covariance matrix that keeps track of correlations between bins for the differential measurement. The theory uncertainty from the NNLO QCD + NLO EW calculation is included by explicitly varying the renormalization and factorization scales, or the parton density functions, in the calculation and registering the variations in the intervals.

Individual 68% and 95% CL bounds on the relevant Wilson coefficients of the SM Effective Field Theory in units of $\text{TeV}^{-2}$. The bounds are derived from the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. The experimental uncertainties are accounted for, in the form of the complete covariance matrix that keeps track of correlations between bins for the differential measurement. The theory uncertainty from the NNLO QCD + NLO EW calculation is included by explicitly varying the renormalization and factorization scales, or the parton density functions, in the calculation and registering the variations in the intervals.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ inclusive measurement. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0,0.3]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.3,0.6]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.6,0.8]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.8,1]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ < 500 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [500,750] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [750,1000] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [1000,1500] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ > 1500 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ < 30 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ $\in$ [30,120] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ > 120 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ inclusive measurement. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0,0.3]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.3,0.6]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.6,0.8]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.8,1]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ < 200 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ $\in$ [200,300] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ $\in$ [300,400] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ > 400 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ < 20 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ $\in$ [20, 70] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ > 70 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ inclusive measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ inclusive measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. The total (stat. + syst.) uncertainties are considered.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement. The total (stat. + syst.) uncertainties are considered.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement. The total (stat. + syst.) uncertainties are considered.

Covariance matrix for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement. The total (stat. + syst.) uncertainties are considered.

Covariance matrix for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement. The total (stat. + syst.) uncertainties are considered.

Covariance matrix for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement. The total (stat. + syst.) uncertainties are considered.

The observed and expected (in parentheses) 95% CL lower mass limits on RPV SUSY, Z′ ($\mathcal{B}=0.1$) , and QBH signals for the e$\mu$, e$\tau$, and $\mu\tau$ channels.

Expected and observed 95% CL upper limits on the product of cross section times branching fraction as a function of the $ au$ sneutrino mass in an RPV SUSY model for the e$\mu$ channel. The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits. The red and blue solid lines show the product of cross section times branching fraction as a function of the tau sneutrino mass for two different values of couplings.

Expected and observed 95% CL upper limits on the product of cross section times branching fraction as a function of the $ au$ sneutrino mass in an RPV SUSY model for the e$\tau$ channel. The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits. The red and blue solid lines show the product of cross section times branching fraction as a function of the tau sneutrino mass for two different values of couplings.

Expected and observed 95% CL upper limits on the product of cross section times branching fraction as a function of the $ au$ sneutrino mass in an RPV SUSY model for the $\mu\tau$ channel. The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits. The red and blue solid lines show the product of cross section times branching fraction as a function of the tau sneutrino mass for two different values of couplings.

Expected (black dashed line) and observed (black solid line) 95% CL upper limits on the product of cross section and branching fraction for a Z′ ($\mathcal{B}=0.1$) boson with LFV decays, in the e$\mu$ channel.The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits. The red solid lines show the product of cross section times branching fraction as a function of the Z′ mass.

Expected (black dashed line) and observed (black solid line) 95% CL upper limits on the product of cross section and branching fraction for a Z′ ($\mathcal{B}=0.1$) boson with LFV decays, in the e$\tau$ channel.The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits. The red solid lines show the product of cross section times branching fraction as a function of the Z′ mass.

Expected (black dashed line) and observed (black solid line) 95% CL upper limits on the product of cross section and branching fraction for a Z′ ($\mathcal{B}=0.1$) boson with LFV decays, in the $\mu\tau$ channel.The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits. The red solid lines show the product of cross section times branching fraction as a function of the Z′ mass.

Expected (black dashed line) and observed (black solid line) 95% CL upper limits on the product of cross section and branching fraction for quantum black hole production in an ADD model with $n=4$ extra dimensions, in the e$\mu$ channel. The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits. The red solid lines show the product of cross section times branching fraction as a function of the QBH threshold mass.

Expected (black dashed line) and observed (black solid line) 95% CL upper limits on the product of cross section and branching fraction for quantum black hole production in an ADD model with $n=4$ extra dimensions, in the e$\tau$ channel. The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits. The red solid lines show the product of cross section times branching fraction as a function of the QBH threshold mass.

Expected (black dashed line) and observed (black solid line) 95% CL upper limits on the product of cross section and branching fraction for quantum black hole production in an ADD model with $n=4$ extra dimensions, in the $\mu\tau$ channel. The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits. The red solid lines show the product of cross section times branching fraction as a function of the QBH threshold mass.

Upper limits at 95% CL on the RPV SUSY model in the plane of $\tau$ sneutrino mass and $\lambda'$ coupling, for four values of $\lambda$ couplings for the e$\mu$ channel. The regions to the left of and above the curves are excluded.

Upper limits at 95% CL on the RPV SUSY model in the plane of $\tau$ sneutrino mass and $\lambda'$ coupling, for four values of $\lambda$ couplings for the e$\tau$ channel. The regions to the left of and above the curves are excluded.

Upper limits at 95% CL on the RPV SUSY model in the plane of $\tau$ sneutrino mass and $\lambda'$ coupling, for four values of $\lambda$ couplings for the $\mu\tau$ channel. The regions to the left of and above the curves are excluded.

Model-independent upper limits at 95% CL on the product of cross section, branching fraction, and acceptance are shown. Observed (expected) limits are shown in black solid (dashed) lines for the e$\mu$ channel. The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits.

Model-independent upper limits at 95% CL on the product of cross section, branching fraction, and acceptance are shown. Observed (expected) limits are shown in black solid (dashed) lines for the e$\tau$ channel. The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits.

Model-independent upper limits at 95% CL on the product of cross section, branching fraction, and acceptance are shown. Observed (expected) limits are shown in black solid (dashed) lines for the $\mu\tau$ channel. The shaded bands represent the one and two standard deviation (s.d.) uncertainties in the expected limits.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "channel_year_binnumber", following the binning from Figure 2.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "channel_year_binnumber", following the binning from Figure 2.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "channel_year_binnumber", following the binning from Figure 2.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "channel_year_binnumber", following the binning used in Figure 2.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "channel_year_binnumber", following the binning used in Figure 2.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "channel_year_binnumber", following the binning used in Figure 2.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=1600\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2000\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2200\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2400\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2600\;\mathrm{GeV}$.

Upper limits at 95% CL on $\sigma\times\mathcal{B}^2(\mathrm{H}_1\rightarrow b\bar{b})$ as a function of $m_{\mathrm{H_1}}$ for $m_{\mathrm{SUSY}}=2800\;\mathrm{GeV}$.

Number of predicted signal events in the $\mu\mu X$ final states, for several benchmark signal hypotheses for Majorana HNL. The quoted uncertainties include statistical and systematic uncertainties.

The 95% CL limits on $|V_{Ne}|^2$ as a function of the HNL mass for a Majorana HNL. Values of $-1$ indicate that no limit is available for the mass point.

The 95% CL limits on $|V_{N\mu}|^2$ as a function of the HNL mass for a Majorana HNL. Values of $-1$ indicate that no limit is available for the mass point.

The 95% CL limits on $|V_{Ne}|^2$ as a function of the HNL mass for a Dirac HNL. Values of $-1$ indicate that no limit is available for the mass point.

The 95% CL limits on $|V_{N\mu}|^2$ as a function of the HNL mass for a Dirac HNL. Values of $-1$ indicate that no limit is available for the mass point.

Covariance matrix for the background prediction in the signal region bins for the $eeX$ final states. The bin labels specify the channel, the $m(\ell_2,\ell_3)$ bin, and the $\Delta_{\mathrm{2D}}$ bin.

Covariance matrix for the background prediction in the signal region bins for the $\mu\mu X$ final states. The bin labels specify the channel, the $m(\ell_2,\ell_3)$ bin, and the $\Delta_{\mathrm{2D}}$ bin.

Correlation matrix for the background prediction the signal region bins for the $eeX$ final states. The bin labels specify the channel, the $m(\ell_2,\ell_3)$ bin, and the $\Delta_{\mathrm{2D}}$ bin.

Correlation matrix for the background prediction the signal region bins for the $\mu\mu X$ final states. The bin labels specify the channel, the $m(\ell_2,\ell_3)$ bin, and the $\Delta_{\mathrm{2D}}$ bin.

Mean lifetime and production cross section for HNLs of Majorana nature, assuming the decay to three charged leptons and coupling to only electrons or muons, as evaluated from simulated event samples.