Contributions from individual sources of systematic uncertainty for the ($p^{\mu}_{T} > 35$, $E_T^{missing} > 35$) GeV kinematic region. All uncertainty values are multiplied by 100. The columns (1-7) correspond to: 1.0 = Electro-Weak background 2.0 = Multi-Jet background 3.0 = Charge mis-identification 4.0 = Relative charge efficiency 5.0 = Magnet polarity weighting 6.0 = Momentum/$E_T^{missing}$ resolution 7.0 = Trigger isolation.

The normalised differential cross section of the azmiuthal decorrelation variable DELTA for the GAMMA+2JET sample for the PT of the second jet in the range 25 TO 30 GeV.

Distribution of $S$, which is one of the discriminating variables used to define the $2\ell^{sc}$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $E_{T}^{miss}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $\Delta R_{\ell^{\pm}\ell^{\pm}}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $m_{x\ell}$ ($x$=3), which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $p_{T}^{leading jet}$, which is one of the discriminating variables used to define the $3\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $E_{T}^{miss}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $\Delta R_{\ell^{\pm}\ell^{\pm}}^{min}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $m_{x\ell}$ ($x$=4), which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Distribution of $p_{T}^{\ell_{1}}$, which is one of the discriminating variables used to define the $4\ell$ SRs. The events are selected with the preselection requirements listed in Table 4 in the paper. The data (dots) are compared with the expected contributions from the relevant background sources (histograms). The expected signal distributions for $m_{H^{\pm\pm}} = 300~GeV$ are also shown, scaled to the observed number of events. The last bin includes overflows.

Contributions from different categories of uncertainties relative to the expected background yields in the defined SRs, as obtained after performing the likelihood ratio test discussed in Section 9 in the paper. The uncertainties are shown for the combination of the individual channels of the $2\ell^{sc}$, $3\ell$ and $4\ell$ SRs. The SRs are indicated along the horizontal axis. In the HEPData entry, the x-axis is simplified for easier visualisation. The first number indicates the sub channel (2:$2\ell^{sc}$, 3:$3\ell$, 4:$4\ell$), while the second number indicates the mass point (2:200, 3:300, 4:400, 5:500).

Data event yields compared with the expected contributions from relevant background sources, for the combination of the individual channels of the $2\ell^{sc}$, $3\ell$ and $4\ell$ SRs. The total uncertainties in the expected event yields are shown as the hatched bands. The SRs are indicated along the horizontal axis. In the HEPData entry, the x-axis is simplified for easier visualisation. The first number indicates the sub channel (2:$2\ell^{sc}$, 3:$3\ell$, 4:$4\ell$), while the second number indicates the mass point (2:200, 3:300, 4:400, 5:500).

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $2\ell^{sc}$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $3\ell$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.

The $E_{T}^{miss}$ distribution for the SRs of the $m_{H^{\pm\pm}} = 300~GeV$ signal mass hypothesis, where the selection requirement on $E_{T}^{miss}$ has been removed. In the attached plot, the signals are stacked on top of the backgrounds while individuals contributions of the $4\ell$ channel are shown in HEPData. The last bin, isolated by a vertical red dashed line, is inclusive and corresponds to the SR.

Observed and expected upper limits of the $H^{\pm\pm}$ pair production cross section times branching fraction at 95% CL obtained from the combination of 2$\ell^{sc}$, 3$\ell$ and 4$\ell$ channels. The region above the observed limit is excluded by the measurement. The bands represent the expected exclusion curves within one and two standard deviations.

The theoretical prediction of Figure 9(a) in the paper.

Observed and expected upper limits of the $H^{\pm\pm}$ and $H^{\pm}$ production cross section times branching fraction at 95% CL obtained from the combination of 2$\ell^{sc}$, 3$\ell$ and 4$\ell$ channels. The region above the observed limit is excluded by the measurement. The bands represent the expected exclusion curves within one and two standard deviations.

The theoretical prediction of Figure 9(b) in the paper.

Data event yields compared with the estimated background in the $m_{H^{\pm\pm}} = 200~GeV$ or $m_{H^{\pm\pm}} = 220~GeV$ SRs. SFOC 0 and SFOC 1,2 refer to the number of same-flavour opposite charge lepton pairs. The total uncertainties in the estimated background yields are shown as the hashed bands. In the HEPData entry, the x-axis is simplified for easier visualisation (1:$e^{\pm}e^{\pm}$, 2:$e^{\pm}\mu^{\pm}$, 3:$\mu^{\pm}\mu^{\pm}$, 4:SFOC 0, 5:SFOC 1,2, 6:$4\ell$).

Data event yields compared with the estimated background in the $m_{H^{\pm\pm}} = 300~GeV$ or $m_{H^{\pm\pm}} = 350~GeV$ SRs. SFOC 0 and SFOC 1,2 refer to the number of same-flavour opposite charge lepton pairs. The total uncertainties in the estimated background yields are shown as the hashed bands. In the HEPData entry, the x-axis is simplified for easier visualisation (1:$e^{\pm}e^{\pm}$, 2:$e^{\pm}\mu^{\pm}$, 3:$\mu^{\pm}\mu^{\pm}$, 4:SFOC 0, 5:SFOC 1,2, 6:$4\ell$).

Data event yields compared with the estimated background in the $m_{H^{\pm\pm}} = 400~GeV$ or $m_{H^{\pm\pm}} = 450~GeV$ SRs. SFOC 0 and SFOC 1,2 refer to the number of same-flavour opposite charge lepton pairs. The total uncertainties in the estimated background yields are shown as the hashed bands. In the HEPData entry, the x-axis is simplified for easier visualisation (1:$e^{\pm}e^{\pm}$, 2:$e^{\pm}\mu^{\pm}$, 3:$\mu^{\pm}\mu^{\pm}$, 4:SFOC 0, 5:SFOC 1,2, 6:$4\ell$).

Data event yields compared with the estimated background in the $m_{H^{\pm\pm}} = 500~GeV$ or $m_{H^{\pm\pm}} = 550~GeV$ or $m_{H^{\pm\pm}} = 600~GeV$ SRs. SFOC 0 and SFOC 1,2 refer to the number of same-flavour opposite charge lepton pairs. The total uncertainties in the estimated background yields are shown as the hashed bands. In the HEPData entry, the x-axis is simplified for easier visualisation (1:$e^{\pm}e^{\pm}$, 2:$e^{\pm}\mu^{\pm}$, 3:$\mu^{\pm}\mu^{\pm}$, 4:SFOC 0, 5:SFOC 1,2, 6:$4\ell$).

Differential fiducial cross section of the Z boson $|y|$ in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the leading b-jet $|y|$ in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta \phi$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta y$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta R$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta \phi$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta y$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta R$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the invariant mass of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the Z boson $p_{\text{T}}$ in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $p_{\text{T}}$ of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the ratio between the $p_{\text{T}}$ and the invariant mass of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Expected and observed 95% CL upper limits on $\sigma(pp\to H)B(H\to Z(\eta_c~\text{or}~J/\psi))/$pb. The smaller (larger) quoted ranges around the expected limits represent $\pm 1\sigma$ ($\pm 2\sigma$) fluctuations.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 0.33 < ln(R/#DeltaR) < 0.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 0.67 < ln(R/#DeltaR) < 1.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 1.00 < ln(R/#DeltaR) < 1.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 1.33 < ln(R/#DeltaR) < 1.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 1.67 < ln(R/#DeltaR) < 2.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 2.00 < ln(R/#DeltaR) < 2.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 2.33 < ln(R/#DeltaR) < 2.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 2.67 < ln(R/#DeltaR) < 3.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 3.00 < ln(R/#DeltaR) < 3.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 3.33 < ln(R/#DeltaR) < 3.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 3.67 < ln(R/#DeltaR) < 4.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 4.00 < ln(R/#DeltaR) < 4.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 0.69 < ln(1/z) < 0.97.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 0.97 < ln(1/z) < 1.25.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 1.25 < ln(1/z) < 1.52.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 1.52 < ln(1/z) < 1.80.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 1.80 < ln(1/z) < 2.08.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.08 < ln(1/z) < 2.36.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.36 < ln(1/z) < 2.63.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.63 < ln(1/z) < 2.91.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.91 < ln(1/z) < 3.19.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 3.19 < ln(1/z) < 3.47.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 3.47 < ln(1/z) < 3.74.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 3.74 < ln(1/z) < 4.02.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.02 < ln(1/z) < 4.30.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.30 < ln(1/z) < 4.57.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.57 < ln(1/z) < 4.85.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.85 < ln(1/z) < 5.13.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 5.13 < ln(1/z) < 5.41.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 5.41 < ln(1/z) < 5.68.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 5.68 < ln(1/z) < 5.96.

The summed covariance matrix of all systematic and statistical uncertainties associated with the measurement in bins of $\ln{(1/z)} \times \ln{(R/\Delta R)}$.

The summed covariance matrix of all statistical uncertainties associated with the measurement in bins of $\ln{(1/z)} \times \ln{(R/\Delta R)}$.

The measured fiducial production cross-sections times branching ratio for $W^+\rightarrow\mu^+\nu$ and $W^-\rightarrow\mu^-\bar{\nu}$, their sum, and their ratio for data

The measured fiducial production cross-sections times branching ratio for $W^+\rightarrow\mu^+ u$ and $W^-\rightarrow\mu^-\bar{\nu}$, their sum, and their ratio for the predictions from DYNNLO (CT14 NNLO PDF set)

Size of the $W^{+}$ the cross-section (differential in $\eta_{\mu}$, as a function of $|\eta_{\mu}|$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. gThe uncertainties are given in percent.

Size of the $W^{+}$ the cross-section (differential in $\eta_{\mu}$, as a function of $|\eta_{\mu}|$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. gThe uncertainties are given in percent.

Size of the asymmetry as a function of $|\eta_{\mu}|$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The uncertainties are given as an absolute difference from the nominal.

Post-fit BDT distributions after the VBF event selection for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. The VBF signal is shown for $\mu = 0.73$ and $\tilde d = -0.01$. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

The observed $\Delta\mathrm{NLL}$ curve as a function of $\tilde d$ values. For comparison, expected $\Delta\mathrm{NLL}$ curves are also shown. The constraints on the nuisance parameters and normalization factors are first determined in a CR-only fit, and then a fit to pseudo-data corresponding to these nuisance parameters, normalization factors, and to $\tilde d=0, \mu = 1$ or $\tilde d =0, \mu = 0.73$ is performed to obtain these $\Delta\mathrm{NLL}$ curves. Moreover, a pre-fit expected $\Delta\mathrm{NLL}$ is shown, using pseudo-data corresponding to $\tilde d =0$ and $\mu = 1$ in the signal and control regions.

The expected $\Delta\mathrm{NLL}$ curves comparing the sensitivity of the fit with and without systematic uncertainties. For comparison, other curves are shown which remove the effect of jet-based systematic uncertainties, $\tau$-based systematic uncertainties, and MC statistical uncertainties.

The observed $\Delta\mathrm{NLL}$ curves for each analysis channel as a function of $\tilde d$, compared to the combined result. For the individual analysis channel $\Delta\mathrm{NLL}$ curves, only event yield information in the other SRs is used, ensuring that the Optimal Observable distributions in the other SRs do not influence the preferred value of $\tilde d$. The signal strength is constrained to be positive in these individual channel $\Delta\mathrm{NLL}$ curves. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit BDT distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit BDT distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit BDT distributions in the $Z\to \ell\ell$ CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the $Z\to \ell\ell$ CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit distribution of weighted event yields as a function of the Optimal Observable for all four SRs combined. The contributions of the different SRs are weighted by a factor of ln(1 + S/B), where S and B are the post-fit expected numbers of signal and background events in that region, respectively. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Dijet invariant mass distribution for the category with at least one b-tagged jet.

Dijet invariant mass distribution for the category with both jets b-tagged.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of q* signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of QBH signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of W' signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of W* signal.

The upper limits on the DM mediator Z' signal at 95% CL from the inclusive category. The 95% CL upper limits are set on the universal quark coupling $g_q$ as a function of the Z' mass.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the mass of the b* signal.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the DM mediator Z' with $g_q$ = 0.25 model.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the SSM Z' model.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the graviton with $k/\overline{M}_{PL}=0.2$ model.

The 95% CL upper limit on the cross-section times kinematic acceptance times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the inclusive category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The 95% CL upper limit on the cross-section times kinematic acceptance times b-tagging efficiency times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the 1 b category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The 95% CL upper limit on the cross-section times kinematic acceptance times b-tagging efficiency times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the 2 b category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The expected 95% CL upper limits on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the DM mediator Z' mass for the current and previous iterations of the analysis. The upper limit of the previous result was obtained with the Bayesian method and is also shown scaled to the 139 fb$^{-1}$ integrated luminosity of the current result to illustrate the effect of the analysis improvements.

Acceptance for the QBH, Q* and W' benchmark signal models in the inclusive category, as a function of the simulated mass.

Acceptance for the DM Z' benchmark signal model for various $g_q$ coupling parameters in the inclusive category, as a function of the simulated mass.

Acceptance for the W* benchmark signal model in the inclusive category, as a function of the simulated mass.

Acceptance times b-tagging efficiency for the b* benchmark signal model in the 1b category, as a function of the simulated mass.

Acceptance times b-tagging efficiency for the DM Z' with $g_q$=0.25, SSM Z' and graviton with $k/\overline{M}_{PL}=0.2$ benchmark signal models in the 2b category, as a function of the simulated mass.

The expected 95% CL upper limits on the cross-section times branching ratio as a function of the DM mediator Z' mass for the current and previous iterations of the analysis. The upper limit of the previous result was obtained with the Bayesian method and is also shown scaled to the 139 fb$^{-1}$ integrated luminosity of the current result to illustrate the effect of the analysis improvements. The current b-tagging requirement is tighter than the previous one for high-$p_T$ jets, resulting in a data sample with limited size for mass above 4 TeV. The background rejection, instead, has improved significantly across the entire mass spectrum inspected by the analysis.