Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of 4. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1800 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of $\geq$ 5. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1900 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of $\geq$ 5. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1800 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

The distribution of $n_{\mathrm{track}}$ in different $S_{\mathrm{ML}}$ regions for simulated background events. Events with 0 $ < S_{\mathrm{ML}} < $ 0.2 (blue), 0.2 $ < S_{\mathrm{ML}} < $ 0.6 (red), and 0.6 $ < S_{\mathrm{ML}} < $ 1.0 (green) are compared. All distributions are normalized to unity. The similar $n_{\mathrm{track}}$ distributions demonstrate that $n_{\mathrm{track}}$ and $S_{\mathrm{ML}}$ are decorrelated.

The distribution of $d_{\mathrm{BV}}$ in $K_{\mathrm{S}}^{\mathrm{0}}$ vertices in data (black) and simulation (purple). The lower panel shows the ratio between data and simulation.

The vertex reconstruction efficiency for artificially displaced vertices in data (black) and simulation (red). In this example, the artificially displaced vertices are corrected to mimic split-SUSY signal events with gluino mass of 2000 GeV and neutralino mass of 1800 GeV.

The ML tagging efficiency for artificially displaced vertices in data (black) and simulation (red). In this example, the artificially displaced vertices are corrected to mimic split-SUSY signal events with gluino mass of 2000 GeV and neutralino mass of 1800 GeV.

Number of predicted and observed events in the control, validation, and search regions. Predictions are calculated using Eqs. (2) and (3) and fitting the data under the background-only hypothesis. Regions are organized by $S_{\mathrm{ML}}$ and $n_{\mathrm{track}}$ values, and region names corresponding with Fig. 7 are given in parentheses. The predicted number of events that pass the $S_{\mathrm{ML}}$ selection and the observed number of events that pass or fail the $S_{\mathrm{ML}}$ selection are shown in seperate rows.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY signal model with a mass splitting of 100 GeV, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY signal model with a mass splitting of 100 GeV, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY model with a $c\tau$ of 10 mm, shown as a function of gluino mass and mass splitting. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY model with a $c\tau$ of 10 mm, shown as a function of gluino mass and mass splitting. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the GMSB SUSY signal model, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the GMSB SUSY signal model, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

Efficiency of signal region A for all signal samples interperted in this search. The cross sections of all signal models are taken to be 1 fb in the table.

Dilepton mass spectra for the high mass $W\phi($ee$)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $W\phi($ee$)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $W\phi($ee$)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $Z\phi($ee$)$ SR event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $Z\phi($ee$)$ SR event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi($ee$)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi($ee$)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi($ee$)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi($ee$)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi($ee$)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi($ee$)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $W\phi(\mu\mu)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $W\phi(\mu\mu)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $W\phi(\mu\mu)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $W\phi(\mu\mu)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $Z\phi(\mu\mu)$ SR event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $Z\phi(\mu\mu)$ SR event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi(\mu\mu)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi(\mu\mu)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi(\mu\mu)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi(\mu\mu)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi(\mu\mu)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi(\mu\mu)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $W\phi(\tau\tau)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $Z\phi(\tau\tau)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $W\phi(\tau\tau)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $Z\phi(\tau\tau)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $W\phi(\tau\tau)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $Z\phi(\tau\tau)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR4 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR5 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR6 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR7 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with scalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with pseudoscalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with scalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with pseudoscalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with scalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with pseudoscalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with scalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with pseudoscalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with scalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with pseudoscalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with scalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with pseudoscalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with H-like production in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with H-like production in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with H-like production in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with H-like production in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with H-like production in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with H-like production in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with scalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with pseudoscalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with scalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with pseudoscalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with scalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with pseudoscalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on $g^2_{tS}$ for the dilaton-like $t\bar{t} \phi$ signal model. Masses of the $\phi$ boson above 300 GeV are not probed for the dilaton-like signal model as the $\phi$ branching fraction into top quark-antiquark pairs becomes nonnegligible.

The 95% confidence level upper limits on $g^2_{tPS}$ for the axion-like $t\bar{t} \phi$ signal model. Masses of the $\phi$ boson above 300 GeV are not probed for the axion-like signal model as the $\phi$ branching fraction into top quark-antiquark pairs becomes nonnegligible.

The 95% confidence level upper limits on the product of $sin^2 \theta$ and branching fraction for the H-like production of X$\phi \rightarrow$ ee. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level upper limits on the product of $sin^2 \theta$ and branching fraction for the H-like production of X$\phi \rightarrow \mu\mu$. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level upper limits on $sin^2 \theta$ for the H-like production and decay of X$\phi$ signal model.

Cross section in units of pb for the W$\phi$, Z$\phi$, and $t\bar{t}\phi$ signals as a function of the $\phi$ boson mass in GeV. All cross sections are inclusive of all W, Z, $t\bar{t}$ and $\phi$ decay modes.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tS}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $t\bar{t} \phi$ signal with scalar couplings, where $g_{tS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tS}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $t\bar{t} \phi$ signal with scalar couplings, where $g_{tS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tS}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $t\bar{t} \phi$ signal with scalar couplings, where $g_{tS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tPS}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $t\bar{t} \phi$ signal with pseudoscalar couplings, where $g_{tPS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tPS}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $t\bar{t} \phi$ signal with pseudoscalar couplings, where $g_{tPS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tPS}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $t\bar{t} \phi$ signal with pseudoscalar couplings, where $g_{tPS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $t\bar{t} \phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $t\bar{t} \phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $t\bar{t} \phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $W\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $W\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $W\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $W\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $W\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $W\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $W\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $W\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $W\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $Z\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $Z\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $Z\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $Z\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $Z\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $Z\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $Z\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $Z\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $Z\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level observed upper limits on the product of $\sigma$($t \bar{t} \phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $t \bar{t} \phi$ signal with scalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $t \bar{t} \phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($t \bar{t} \phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $t \bar{t} \phi$ signal with pseudoscalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $t \bar{t} \phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($W\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with scalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $W\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($W\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with pseudoscalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $W\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($W\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with H-like production, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $W\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($Z\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with scalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $Z\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($Z\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with pseudoscalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $Z\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($Z\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with H-like production, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $Z\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $g^{2}_{tS}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $t \bar{t} \phi$ signal with scalar couplings, where $g_{tS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $g^{2}_{tPS}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $t \bar{t} \phi$ signal with pseudoscalar couplings, where $g_{tPS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $t \bar{t} \phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

Selected signal shapes of the $W\phi$(ee) signal for illustration purposes. All shape parametrizations for all coupling scenarios of the $X\phi$(ee) signal are provided in the SignalShapes_XPhiToEleEle.root file, and a README file with instructions is provided under Additional Resources.

Selected signal shapes of the $W\phi$$(\mu\mu)$ signal for illustration purposes. All shape parametrizations for all coupling scenarios of the $X\phi$$(\mu\mu)$ signal are provided in the SignalShapes_XPhiToMuMu.root file, and a README file with instructions is provided under Additional Resources.

Selected signal shapes of the $W\phi$$(\tau\tau)$ signal for illustration purposes. All shape parametrizations for all coupling scenarios of the $X\phi$$(\tau\tau)$ signal are provided in the SignalShapes_XPhiToTauTau.root file, and a README file with instructions is provided under Additional Resources.

The distribution of $\mathrm{p_{T}^{miss}}$ (GeV) in 2L1T MisID CR events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction.

The distribution of $\mathrm{M_{T}}$ (GeV) in 3L OnZ CR events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction.

The distribution of $\mathrm{H_{T}}$ (GeV) in 3L ttZ CR events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction.

Distribution of BDT score from the SS-M ($\mathrm{B_{e}=B_{\mu}=B_{\tau}}$) BDT for the 3L+2L1T CR events for the combined 2016-2018 data set. The 3L+2L1T CR consists of the 3L OnZ, 3L Z$\mathrm{\gamma}$, and 2L1T MisID CRs. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction.

The distribution of visible diboson $\mathrm{p_{T}}$ (GeV) in 4L ZZ CR events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction.

Distribution of BDT score from the SS-M ($\mathrm{B_{e}=B_{\mu}=B_{\tau}}$) BDT for the 4L ZZ CR events for the combined 2016-2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction.

The distribution of $\mathrm{L_{T}}$ in all seven multilepton channels for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton of $\mathrm{m_{\tau'}}$ = 1 TeV in the doublet scenario, before the fit, is also overlaid.

The distribution of $\mathrm{p_{T}^{miss}}$ in all seven multilepton channels for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermion of $\mathrm{m_{\Sigma}}$ = 1 TeV in the flavor-democratic scenario, before the fit, is also overlaid.

The distribution of $\mathrm{H_{T}}$ in all seven multilepton channels for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. For illustration, an example signal hypothesis for the production of the scalar leptoquark of $mathrm{m_{S}}$ = 1 TeV coupled to a top quark and a $\tau$ lepton, before the fit, is also overlaid.

The distribution of $\mathrm{M_{OSSF}}$ in channels with at least one light lepton pair (4L, 3L1T, 3L, 2L2T, and 2L1T) for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermion of $\mathrm{m_{\Sigma}}$ = 1 TeV in the flavor-democratic scenario, before the fit, is also overlaid.

The $\mathrm{N_{b}}$ distribution in 4L, 3L1T, 3L, 2L2T, 2L1T, 1L3T, and 1L2T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction.

The invariant mass distribution of the opposite-sign same-flavor ($\mathrm{M_{OSSF}}$) tau lepton pair distribution in 2L2T, 1L3T, and 1L2T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction.

The $\mathrm{M_{T}^{12}}$ distribution in 4L, 3L1T, 3L, 2L2T, 2L1T, 1L3T, and 1L2T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction.

The $\mathrm{N_{b}}$ distribution in 3L, 2L1T, and 1L2T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The $\mathrm{L_{T}}$ distribution in 3L, 2L1T, and 1L2T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The $\mathrm{p_{T}^{miss}}$ distribution in 3L, 2L1T, and 1L2T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The $\mathrm{H_{T}}$ distribution in 3L, 2L1T, and 1L2T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The $\mathrm{M_{OSSF}}$ distribution in 3L, and 2L1T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The invariant mass distribution of the opposite-sign different-flavor ($\mathrm{M_{OSDF}}$) light lepton pair distribution in 3L, and 2L1T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The invariant mass distribution of the opposite-sign same-flavor ($\mathrm{M_{OSSF}}$) tau lepton pair distribution in 1L2T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The invariant mass distribution of the opposite-sign different-flavor ($\mathrm{M_{OSDF}}$) light lepton and tau lepton pair distribution in 2L1T, and 1L2T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The $\mathrm{M_{T}^{1}}$ distribution in 3L, 2L1T, and 1L2T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The $\mathrm{M_{T}^{12}}$ distribution in 3L, 2L1T, and 1L2T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The model independent fundamental table categories for the combined 2016-2018 data set, as defined in Table 1. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The $\mathrm{N_{b}}$ distribution in 4L, 3L1T, 2L2T, and 1L3T events for the combined 2016-2018 data set. The rightmost bin contains the overflow events. The gray band represents the sum of statistical and systematic uncertainties on the SM background predictions.

The SR distributions of the fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table in 3L channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table in 2L1T and 1L2T channels for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table in 3L channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table in 2L1T and 1L2T channels for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table in 4L, 3L1T, 2L2T, and 1L3T channels for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table in 4L, 3L1T, 2L2T, and 1L3T channels for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the fundamental $\mathrm{S_{T}}$ table in 3L channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the fundamental $\mathrm{S_{T}}$ table in 2L1T and 1L2T channels for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the fundamental $\mathrm{S_{T}}$ table in 3L channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the fundamental $\mathrm{S_{T}}$ table in 2L1T and 1L2T channels for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the fundamental $\mathrm{S_{T}}$ table in 4L, 3L1T, 2L2T, and 1L3T channels for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the fundamental $\mathrm{S_{T}}$ table in 4L, 3L1T, 2L2T, and 1L3T channels for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 3L 0B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 3L 1B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 3L 2B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 2L1T 0B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 2L1T 1B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 2L1T 2B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 1L2T channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. An example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid. For this category, the signal yield is negligible and is not visible in the figure.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 3L 0B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 3L 1B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 3L 2B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 2L1T 0B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 2L1T 1B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 2L1T 2B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 1L2T channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. An example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid. For this category, the signal yield is negligible and is not visible in the figure.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 4L 0B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 4L 1B/2B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 3L1T, 2L2T, and 1L3T channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 4L 0B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 4L 1B/2B channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The SR distributions of the advanced $\mathrm{S_{T}}$ table in 3L1T, 2L2T, and 1L3T channel for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The VLL-L BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 200 GeV, before the fit, is also overlaid.

The VLL-L BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 200 GeV, before the fit, is also overlaid.

The VLL-L BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 200 GeV, before the fit, is also overlaid.

The VLL-M BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 400 GeV, before the fit, is also overlaid.

The VLL-M BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 400 GeV, before the fit, is also overlaid.

The VLL-M BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 400 GeV, before the fit, is also overlaid.

The VLL-H BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 900 GeV, before the fit, is also overlaid.

The VLL-H BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 900 GeV, before the fit, is also overlaid.

The VLL-H BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 900 GeV, before the fit, is also overlaid.

The VLL-L BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 200 GeV, before the fit, is also overlaid.

The VLL-L BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 200 GeV, before the fit, is also overlaid.

The VLL-L BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 200 GeV, before the fit, is also overlaid.

The VLL-M BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 400 GeV, before the fit, is also overlaid.

The VLL-M BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 400 GeV, before the fit, is also overlaid.

The VLL-M BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 400 GeV, before the fit, is also overlaid.

The VLL-H BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 900 GeV, before the fit, is also overlaid.

The VLL-H BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 900 GeV, before the fit, is also overlaid.

The VLL-H BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 900 GeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 850 GeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 850 GeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 850 GeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 100 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-L $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 300 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-M $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 550 GeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 850 GeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 850 GeV, before the fit, is also overlaid.

The SS-H $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the scenario with mixing exclusively to $\tau$ lepton for $\mathrm{m_{\Sigma}}$ = 850 GeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 1.2 TeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 1.2 TeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{\tau}=1}$ BDT regions for the 3-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 1.2 TeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 1.2 TeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 1.2 TeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{e}+B_{\mu}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and an electron for $\mathrm{m_{S}}$ = 1.2 TeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-VL $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 200 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-L $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 400 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-M $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 700 GeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2016 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 1.2 TeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2017 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 1.2 TeV, before the fit, is also overlaid.

The LQ-H $\mathrm{B_{\tau}=1}$ BDT regions for the 4-lepton channels for the 2018 data set. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background prediction. The expected SM background distributions and the uncertainties are shown after fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a $\tau lepton$ for $\mathrm{m_{S}}$ = 1.2 TeV, before the fit, is also overlaid.

Observed and expected upper limits at 95%% CL on the production cross section for the type-III seesaw fermions in the flavor-democratic scenario using the table schemes and the BDT regions of the SS-M and the SS-H $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ BDTs. To the left of the vertical dashed gray line, the limits are shown from the advanced $\mathrm{S_{T}}$ table, and to the right the limits are shown from the BDT regions.

Observed and expected upper limits at 95%% CL on the production cross section for the vector-like $\mathrm{\tau}$ leptons: doublet model. To the left of the vertical dashed gray line, the limits are shown from the advanced $\mathrm{S_{T}}$ table, and to the right the limits are shown from the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks: $\mathrm{B_{\tau}=1}$ and $\mathrm{\beta=1}$. To the left of the vertical dashed gray line, the limits are shown from the advanced $\mathrm{S_{T}}$ table, and to the right the limits are shown from the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks: $\mathrm{B_{e}=1}$ and $\mathrm{\beta=1}$. To the left of the vertical dashed gray line, the limits are shown from the advanced $\mathrm{S_{T}}$ table, and to the right the limits are shown from the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks: $\mathrm{B_{\mu}=1}$ and $\mathrm{\beta=1}$. To the left of the vertical dashed gray line, the limits are shown from the advanced $\mathrm{S_{T}}$ table, and to the right the limits are shown from the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the vector-like $\mathrm{\tau}$ leptons: singlet model. The limit is shown from the advanced $\mathrm{S_{T}}$ table for all masses.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ scenario using the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ scenario using the Fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ scenario using the Fundamental $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{e}=B_{\mu}=B_{\tau}}$ scenario using the Advanced $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{\mu}=1}$ scenario using the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{\mu}=1}$ scenario using the Fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{\mu}=1}$ scenario using the Fundamental $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{\mu}=1}$ scenario using the Advanced $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{e}=1}$ scenario using the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{e}=1}$ scenario using the Fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{e}=1}$ scenario using the Fundamental $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{e}=1}$ scenario using the Advanced $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{\tau}=1}$ scenario using the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{\tau}=1}$ scenario using the Fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{\tau}=1}$ scenario using the Fundamental $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the type-III seesaw fermions in the $\mathrm{B_{\tau}=1}$ scenario using the Advanced $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{\mu}=1}$ scenario using the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{\mu}=1}$ scenario using the Fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{\mu}=1}$ scenario using the Fundamental $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{\mu}=1}$ scenario using the Advanced $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{e}=1}$ scenario using the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{e}=1}$ scenario using the Fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{e}=1}$ scenario using the Fundamental $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{e}=1}$ scenario using the Advanced $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{\tau}=1}$ scenario using the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{\tau}=1}$ scenario using the Fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{\tau}=1}$ scenario using the Fundamental $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{\tau}=1}$ scenario using the Advanced $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the vector-like $\tau$ leptons in the doublet scenario using the BDT regions.

Observed and expected upper limits at 95% CL on the production cross section for the vector-like $\tau$ leptons in the doublet scenario using the Fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the vector-like $\tau$ leptons in the doublet scenario using the Fundamental $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the vector-like $\tau$ leptons in the doublet scenario using the Advanced $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the vector-like $\tau$ leptons in the singlet scenario using the Fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the vector-like $\tau$ leptons in the singlet scenario using the Fundamental $\mathrm{S_{T}}$ table.

Observed and expected upper limits at 95% CL on the production cross section for the vector-like $\tau$ leptons in the singlet scenario using the Advanced $\mathrm{S_{T}}$ table.

Observed lower limits at 95% CL on the mass of the type-III seesaw fermions in the plane defined by $\mathrm{B_{e}}$ and $\mathrm{B_{\tau}}$, with the constraint that $\mathrm{B_{e}+B_{\mu}+B_{\tau}=1}$. These limits arise from the SS-H $\mathrm{B_{\tau}=1}$ BDT when $\mathrm{B_{\tau}\geq0.9}$, and by the SS-H $\mathrm{B_{e}+B_{\mu}+B_{\tau}=1}$ BDT for the other decay branching fraction combinations.

Median Expected lower limits at 95% CL on the mass of the type-III seesaw fermions in the plane defined by $\mathrm{B_{e}}$ and $\mathrm{B_{\tau}}$, with the constraint that $\mathrm{B_{e}+B_{\mu}+B_{\tau}=1}$. These limits arise from the SS-H $\mathrm{B_{\tau}=1}$ BDT when $\mathrm{B_{\tau}\geq0.9}$, and by the SS-H $\mathrm{B_{e}+B_{\mu}+B_{\tau}=1}$ BDT for the other decay branching fraction combinations.

Acceptance times efficiency values for the major SM backgrounds WZ, ZZ, and ttZ in the signal regions of all seven multilepton channels. The product is defined as the ratio of the total reconstructed yield in a given channel (after all the corrections and scale factor implementation) to the product of luminosity and the production cross section of the given simulation sample. The statistical uncertainty on the acceptance times efficiency values is insignificant with respect to the quoted precision.

Acceptance times efficiency values with statistical uncertainty for the vector-like $\mathrm{\tau}$ lepton model in the doublet scenario in the signal regions of all seven multilepton channels. The product is defined as the ratio of the total reconstructed yield in a given channel (after all the corrections and scale factor implementation) to the product of luminosity and the production cross section of the given simulation sample.

Acceptance times efficiency values with statistical uncertainty for the vector-like $\mathrm{\tau}$ lepton model in the singlet scenario in the signal regions of all seven multilepton channels. The product is defined as the ratio of the total reconstructed yield in a given channel (after all the corrections and scale factor implementation) to the product of luminosity and the production cross section of the given simulation sample.

Acceptance times efficiency values with statistical uncertainty for the type-III seesaw fermions in the $\mathrm{(B_{e}=B_{\mu}=B_{\tau})}$ scenario in the signal regions of all seven multilepton channels. The product is defined as the ratio of the total reconstructed yield in a given channel (after all the corrections and scale factor implementation) to the product of luminosity and the production cross section of the given simulation sample.

Acceptance times efficiency values with statistical uncertainty for the type-III seesaw fermions in the $\mathrm{(B_{e}=1)}$ scenario in the signal regions of all seven multilepton channels. The product is defined as the ratio of the total reconstructed yield in a given channel (after all the corrections and scale factor implementation) to the product of luminosity and the production cross section of the given simulation sample.

Acceptance times efficiency values with statistical uncertainty for the type-III seesaw fermions in the $\mathrm{(B_{\mu}=1)}$ scenario in the signal regions of all seven multilepton channels. The product is defined as the ratio of the total reconstructed yield in a given channel (after all the corrections and scale factor implementation) to the product of luminosity and the production cross section of the given simulation sample.

Acceptance times efficiency values with statistical uncertainty for the type-III seesaw fermions in the $\mathrm{B_{\tau}=1)}$ scenario in the signal regions of all seven multilepton channels. The product is defined as the ratio of the total reconstructed yield in a given channel (after all the corrections and scale factor implementation) to the product of luminosity and the production cross section of the given simulation sample.

Acceptance times efficiency values with statistical uncertainty for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{\tau}=1}$ scenario in the signal regions of all seven multilepton channels. The product is defined as the ratio of the total reconstructed yield in a given channel (after all the corrections and scale factor implementation) to the product of luminosity and the production cross section of the given simulation sample.

Acceptance times efficiency values with statistical uncertainty for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{e}=1}$ scenario in the signal regions of all seven multilepton channels. The product is defined as the ratio of the total reconstructed yield in a given channel (after all the corrections and scale factor implementation) to the product of luminosity and the production cross section of the given simulation sample.

Acceptance times efficiency values with statistical uncertainty for the scalar leptoquarks with $\mathrm{\beta=1}$ in the $\mathrm{B_{\mu}=1}$ scenario in the signal regions of all seven multilepton channels. The product is defined as the ratio of the total reconstructed yield in a given channel (after all the corrections and scale factor implementation) to the product of luminosity and the production cross section of the given simulation sample.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of SM gauge bosons (W/Z/h) for 0.2<dRmin<0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of SM gauge bosons (W/Z/h) for dRmin>0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for 0.2<dRmin<0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for dRmin>0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for 0.2<dRmin<0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for dRmin>0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of $\tau$ leptons for 0.2<dRmin<0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of $\tau$ leptons for dRmin>0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of $\tau$ leptons for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of $\tau$ leptons for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of $\tau$ leptons for 0.2<dRmin<0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of $\tau$ leptons for dRmin>0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of $\tau$ leptons for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of $\tau$ leptons for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of $\tau$ leptons for 0.2<dRmin<0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of $\tau$ leptons for dRmin>0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of $\tau$ leptons for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for electrons in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of $\tau$ leptons for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|<1.2}$ region, arising from the decay of SM gauge bosons (W/Z/h) for 0.2<dRmin<0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|<1.2}$ region, arising from the decay of SM gauge bosons (W/Z/h) for dRmin>0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|<1.2}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|<1.2}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|>1.2}$ region, arising from the decay of SM gauge bosons (W/Z/h) for 0.2<dRmin<0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|>1.2}$ region, arising from the decay of SM gauge bosons (W/Z/h) for dRmin>0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|>1.2}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|>1.2}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|<1.2}$ region, arising from the decay of $\tau$ leptons for 0.2<dRmin<0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|<1.2}$ region, arising from the decay of $\tau$ leptons for dRmin>0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|<1.2}$ region, arising from the decay of $\tau$ leptons for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|<1.2}$ region, arising from the decay of $\tau$ leptons for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|>1.2}$ region, arising from the decay of $\tau$ leptons for 0.2<dRmin<0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|>1.2}$ region, arising from the decay of $\tau$ leptons for dRmin>0.4. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|>1.2}$ region, arising from the decay of $\tau$ leptons for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for muons in the $\mathrm{|\eta|>1.2}$ region, arising from the decay of $\tau$ leptons for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 1-prong $\tau_{h}$ in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of SM gauge bosons (W/Z/h) for dRmin>0.2. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 1-prong $\tau_{h}$ in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 1-prong $\tau_{h}$ in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 1-prong $\tau_{h}$ in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for dRmin>0.2. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 1-prong $\tau_{h}$ in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 1-prong $\tau_{h}$ in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 1-prong $\tau_{h}$ in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for dRmin>0.2. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 1-prong $\tau_{h}$ in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 1-prong $\tau_{h}$ in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 3-prong $\tau_{h}$ in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of SM gauge bosons (W/Z/h) for dRmin>0.2. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 3-prong $\tau_{h}$ in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 3-prong $\tau_{h}$ in the $\mathrm{|\eta|<1.1}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 3-prong $\tau_{h}$ in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for dRmin>0.2. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 3-prong $\tau_{h}$ in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 3-prong $\tau_{h}$ in the $\mathrm{|\eta|>1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 3-prong $\tau_{h}$ in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for dRmin>0.2. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 3-prong $\tau_{h}$ in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}<2}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

Reconstruction efficiency and associated uncertainty maps for 3-prong $\tau_{h}$ in the $\mathrm{1.1<|\eta|<1.6}$ region, arising from the decay of SM gauge bosons (W/Z/h) for $\mathrm{N_{j}>1}$. The lepton efficiency is estimated in a simulated event sample for the ZZ process. For a given input generator-level $\mathrm{p_{T}}$, the efficiency map provides the probability distribution of the reconstructed $\mathrm{p_{T}}$, accounting for reconstruction and identification efficiency, and the $\mathrm{p_{T}}$ resolution. The x-axis and the y-axis represent bins in the reconstructed and generated lepton $\mathrm{p_{T}}$, respectively.

The SR distributions of the Fundamental $\mathrm{L_{T}+p_{T}^{miss}}$ table for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the type-III seesaw heavy fermions in the flavor-democratic scenario for $\mathrm{m_{\Sigma}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the Fundamental $\mathrm{S_{T}}$ table for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the vector-like $\tau$ lepton in the doublet scenario for $\mathrm{m_{\tau'}}$ = 1 TeV, before the fit, is also overlaid.

The SR distributions of the Advanced $\mathrm{S_{T}}$ table for the combined 2016-2018 data set. The detailed description of the bin numbers can be found in Tables 3-6 in the paper. The lower panel shows the ratio of observed events to the total expected background prediction. The gray band on the ratio represents the sum of statistical and systematic uncertainties in the SM background predictions. The expected SM background distributions and the uncertainties are shown before fitting the data under the background-only hypothesis. For illustration, an example signal hypothesis for the production of the scalar leptoquark coupled to a top quark and a muon for $\mathrm{m_{S}}$ = 1.4 TeV, before the fit, is also overlaid.