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.

Dijet signal efficiencies as a function of the signal mass and lifetime for events satisfying all event and vertex requirements, with corrections based on systematic differences in the vertex reconstruction efficiency between data and simulation.

The distribution of $d_{\mathrm{BV}}$ for $\geq$5-track one-vertex events in data and three simulated multijet signal samples each with a mass of 1600 GeV. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for the samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The last bin includes the overflow events. This bin includes one event in data with a vertex with large $d_{\mathrm{BV}}$ that appears to arise from tracks originating from separate pp interaction vertices, consistent with background.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 3-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 4-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized using $\geq$5-track one-vertex event information. The two vertical red dashed lines separate the regions used in the fit.

Predicted yields for the background-only normalized template, predicted yields for three simulated multijet signals each with a mass of 1600 GeV, and the observed yield in each $d_{\mathrm{VV}}$ bin. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The uncertainties in the signal yields and the systematic uncertainties in the background prediction reflect the systematic uncertainties given in the text.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals. For $m$ = 2400 GeV, the expected neutralino cross section is $\approx 8\times 10^{-5}$ fb and is not shown.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Data-to-simulation efficiency correction factors, shown for multijet and dijet signal topologies in several ranges of $c\tau$. Note that these correction factors account for the two long-lived particles in the simulated events, and are therefore the total correction factors used to scale event yields rather than the correction factors one would apply to individual vertices.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 3-track one-vertex events in data, and is normalized to the number of 3-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track and 3-track one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track one-vertex events in data, and is normalized to the number of 4-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from $\geq$5-track one-vertex events in data, and is normalized using one-vertex event information. No $\geq$5-track two-vertex data events pass the selection.