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.