The fraction of photon-induced background as compared with the total amount of DY signal plus photon-induced events ($N_{\gamma\gamma}/(N_{\gamma\gamma} + N_\mathrm{DY})$) in different dilepton mass bins. These numbers are averaged across the different years and channels.

Exclusion limits at 95% CL on the coupling K_L for a Z' in the sequential standard model as a function of the Z' mass.

Covariance matrix of the Absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $p_T^{t,h}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $p_T^{t,l}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $p_T^{t,l}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $m^{t\bar{t}}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $m^{t\bar{t}}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $y^{t,h}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $y^{t,h}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $y^{t,l}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $y^{t,l}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $y^{t\bar{t}}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $y^{t\bar{t}}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $H_T^{t\bar{t}}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $H_T^{t\bar{t}}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $p_T^{t\bar{t}}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $p_T^{t\bar{t}}$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $p_T^{t\bar{t}}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $H_T^{t\bar{t}+jets}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $H_T^{t\bar{t}+jets}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $N^j$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $N^j$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $N^j$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $p_T^{j,1}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $p_T^{j,1}$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $p_T^{j,1}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $m(j_1, t_h)$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $m(j_1, t_h)$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $m(j_1, t_h)$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of $p_T^{j,2}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of $p_T^{j,2}$ at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative differential cross-section as a function of $p_T^{j,2}$ at particle level in the boosted topology. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of $p_T^{j,1}$ vs $N^j$ at particle level in the boosted topology in $N^j$ = 1. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of $p_T^{j,1}$ vs $N^j$ at particle level in the boosted topology in $N^j$ = 2. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of $p_T^{j,1}$ vs $N^j$ at particle level in the boosted topology in $N^j$ $\geq$ 3. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $p_T^{j,1}$ vs $N^j$ at particle level in the boosted topology in $N^j$ = 1. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $p_T^{j,1}$ vs $N^j$ at particle level in the boosted topology in $N^j$ = 2. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $p_T^{j,1}$ vs $N^j$ at particle level in the boosted topology in $N^j$ $\geq$ 3. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ = 1 and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ = 1 at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ = 2 and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ = 1 at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ = 2 and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ = 2 at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ $\geq$ 3 and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ = 1 at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ $\geq$ 3 and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ = 2 at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ $\geq$ 3 and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $N^j$ in $N^j$ $\geq$ 3 at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative double-differential cross-section as a function of $p_T^{j,1}$ vs $p_T^{t,h}$ at particle level in the boosted topology in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of $p_T^{j,1}$ vs $p_T^{t,h}$ at particle level in the boosted topology in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of $p_T^{j,1}$ vs $p_T^{t,h}$ at particle level in the boosted topology in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $p_T^{j,1}$ vs $p_T^{t,h}$ at particle level in the boosted topology in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $p_T^{j,1}$ vs $p_T^{t,h}$ at particle level in the boosted topology in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $p_T^{j,1}$ vs $p_T^{t,h}$ at particle level in the boosted topology in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of $p_T^{j,1}$ vs $p_T^{t,h}$ in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ at particle level in the boosted topology in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ at particle level in the boosted topology in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ at particle level in the boosted topology in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ at particle level in the boosted topology in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ at particle level in the boosted topology in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ at particle level in the boosted topology in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 355.0 GeV < $p_T^{t,h}$ < 398.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 398.0 GeV < $p_T^{t,h}$ < 496.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $p_T^{t,h}$ in 496.0 GeV < $p_T^{t,h}$ < 2000.0 GeV at particle level in the boosted topology, accounting for the statistical uncertainties.

Relative double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ at particle level in the boosted topology in $N^j$ = 1. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ at particle level in the boosted topology in $N^j$ = 2. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ at particle level in the boosted topology in $N^j$ $\geq$ 3. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ at particle level in the boosted topology in $N^j$ = 1. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ at particle level in the boosted topology in $N^j$ = 2. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ at particle level in the boosted topology in $N^j$ $\geq$ 3. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ = 1 and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ = 1 at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ = 2 and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ = 1 at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ = 2 and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ = 2 at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ $\geq$ 3 and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ = 1 at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ $\geq$ 3 and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ = 2 at particle level in the boosted topology, accounting for the statistical uncertainties.

Covariance matrix between the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ $\geq$ 3 and the Absolute double-differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ vs $N^j$ in $N^j$ $\geq$ 3 at particle level in the boosted topology, accounting for the statistical uncertainties.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t,h}$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t,l}$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t,l}$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m^{t\bar{t}}$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m^{t\bar{t}}$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m^{t\bar{t}}$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t,h}$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t,h}$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t,h}$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t,h}$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t,l}$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t,l}$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t,l}$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t,l}$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t,l}$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t\bar{t}}$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t\bar{t}}$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t\bar{t}}$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t\bar{t}}$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t\bar{t}}$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $y^{t\bar{t}}$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}}$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}}$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}}$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}}$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}}$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}}$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}}$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ and the absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t\bar{t}}$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t\bar{t}}$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t\bar{t}}$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t\bar{t}}$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t\bar{t}}$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t\bar{t}}$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t\bar{t}}$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t\bar{t}}$ and the absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{t\bar{t}}$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ and the absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ and the absolute differential cross-section as function of $p_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ and the absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ and the absolute differential cross-section as function of $p_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of $p_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $N^j$ and the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $p_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,1}$ and the absolute differential cross-section as function of $N^j$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $p_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $N^j$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $m(j_1, t_h)$ and the absolute differential cross-section as function of $p_T^{j,1}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $p_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $N^j$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $p_T^{j,1}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ and the absolute differential cross-section as function of $m(j_1, t_h)$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $p_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $N^j$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $p_T^{j,1}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $m(j_1, t_h)$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $p_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $N^j$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $p_T^{j,1}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $m(j_1, t_h)$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of Inclusive $\sigma$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $p_T^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $p_T^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $m^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $y^{t,h}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $y^{t,l}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $y^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $H_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(b_{l},t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $p_T^{t\bar{t}}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(t_{h},t_{l})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $H_T^{t\bar{t}+jets}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $N^j$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $p_T^{j,1}$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $m(j_1, t_h)$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_1, t_h)$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{2}, t_{h})$ at particle level in the boosted topology.

Statistical covariance matrix between the absolute differential cross-section as function of $p_T^{j,2}$ and the absolute differential cross-section as function of $\frac{\Delta\phi}{\pi}(j_{1}, j_{2})$ at particle level in the boosted topology.

Differential cross section times dimuon branching fraction of Y(1S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Nuclear modification factor of Y(1S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

RFB of Y(1S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(2S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(3S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(1S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(2S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(3S) versus $E^{|\eta_{lab}|>4}_{T}$.

Nuclear modification factor of Y(1S), Y(2S), and Y(3S) integrated over pT and $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

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.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-0 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-2 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-0.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-2.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-0 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-2 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Comparison of the background model (stacked histograms) with data (dots) in the $2b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $3b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $4b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-0 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-2 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The theoretical prediction for the bulk RS model with $k/\bar{M}_{\text{Pl}} = 1$ (solid red line) is shown; the decrease below 350 GeV is due to a sharp reduction in the $G^{*}_{\text{KK}} \rightarrow HH$ branching ratio. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

GGsm vs -2dNLL (f_{#Lambda1} (u) observed)

GGsm vs -2dNLL (SM-like (f_{ai}=0) expected)

GGsm vs -2dNLL (f_{a2} (u) expected)

GGsm vs -2dNLL (f_{a3} (u) expected)

GGsm vs -2dNLL (f_{#Lambda1} (u) expected)

GGsm vs -2dNLL (2l2#nu off-shell + 4l on-shell observed)

GGsm vs -2dNLL (4l off-shell + 4l on-shell observed)

GGsm vs -2dNLL (2l2#nu off-shell + 4l on-shell expected)

GGsm vs -2dNLL (4l off-shell + 4l on-shell expected)

r_offshell vs -2dNLL (R_{V,F}^{off-shell}=1 observed)

r_offshell vs -2dNLL (R_{V,F}^{off-shell} (u) observed)

r_offshell vs -2dNLL (R_{V,F}^{off-shell}=1 (2l2#nu) observed)

r_offshell vs -2dNLL (R_{V,F}^{off-shell} (u, 2l2#nu) observed)

r_offshell vs -2dNLL (R_{V,F}^{off-shell}=1 expected)

r_offshell vs -2dNLL (R_{V,F}^{off-shell} (u) expected)

r_offshell vs -2dNLL (R_{V,F}^{off-shell}=1 (2l2#nu) expected)

r_offshell vs -2dNLL (R_{V,F}^{off-shell} (u, 2l2#nu) expected)

rf_offshell vs -2dNLL (2l2#nu+4l observed)

rf_offshell vs -2dNLL (Only 2l2#nu observed)

rf_offshell vs -2dNLL (2l2#nu+4l expected)

rf_offshell vs -2dNLL (Only 2l2#nu expected)

rv_offshell vs -2dNLL (2l2#nu+4l observed)

rv_offshell vs -2dNLL (Only 2l2#nu observed)

rv_offshell vs -2dNLL (2l2#nu+4l expected)

rv_offshell vs -2dNLL (Only 2l2#nu expected)

rf_offshell vs rv_offshell observed

fa2 vs -2dNLL (#Gamma_{H}=4.1 MeV observed)

fa2 vs -2dNLL (#Gamma_{H} (u) observed)

fa2 vs -2dNLL (On-shell 4l observed)

fa2 vs -2dNLL (#Gamma_{H}=4.1 MeV expected)

fa2 vs -2dNLL (#Gamma_{H} (u) expected)

fa2 vs -2dNLL (On-shell 4l expected)

fa3 vs -2dNLL (#Gamma_{H}=4.1 MeV observed)

fa3 vs -2dNLL (#Gamma_{H} (u) observed)

fa3 vs -2dNLL (On-shell 4l observed)

fa3 vs -2dNLL (#Gamma_{H}=4.1 MeV expected)

fa3 vs -2dNLL (#Gamma_{H} (u) expected)

fa3 vs -2dNLL (On-shell 4l expected)

fL1 vs -2dNLL (#Gamma_{H}=4.1 MeV observed)

fL1 vs -2dNLL (#Gamma_{H} (u) observed)

fL1 vs -2dNLL (On-shell 4l observed)

fL1 vs -2dNLL (#Gamma_{H}=4.1 MeV expected)

fL1 vs -2dNLL (#Gamma_{H} (u) expected)

fL1 vs -2dNLL (On-shell 4l expected)

J/PSI(1S)-->MU+MU- nuclear modification in p+Al collisions as a function of rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

PSI(2S)-->MU+MU- nuclear modification in p+Au collisions as a function of rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/PSI(1S)-->MU+MU- nuclear modification in p+Au collisions as a function of rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

PSI(2S)-->MU+MU- nuclear modification in p+Au collisions as a function of centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

PSI(2S)-->MU+MU- nuclear modification in p+Au collisions as a function of centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/PSI(1S)-->MU+MU- nuclear modification in p+Au collisions as a function of centrality. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

Invariant yield ratio of PSI(2S)-->MU+MU- with respect to J/PSI(1S)-->MU+MU- in p+p collisions at rapidity 1.2 < |y| < 2.2. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. The global systematic uncertainty cancels when taking the ratio.

Invariant yield ratio of PSI(2S)-->MU+MU- with respect to J/PSI(1S)-->MU+MU- in p+Au collisions as a function of centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. The global systematic uncertainty cancels when taking the ratio.

Invariant yield ratio of PSI(2S)-->MU+MU- with respect to J/PSI(1S)-->MU+MU- in p+Au collisions as a function of centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. The global systematic uncertainty cancels when taking the ratio.