Signal selection efficiency times acceptance for charged Higgs boson with masses between $60\,$GeV to $168\,$GeV. The efficiency times acceptance drops at high masses rapidly due to the low momentum of the $b$-quark produced together with the charged Higgs boson. For the $168\,$GeV mass point the efficiency times acceptance is larger than for the $160\,$GeV mass point due to having a higher probability that the top quark, which decays into the charged Higgs boson and a $b$-quark, is off-shell.

Best fit values of $\mu_{t\bar{t}}$ and $f_{\mathrm{HF}}$ parameters at charged Higgs masses between $60\,$GeV and $168\,$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$).

Expected and observed upper limits on $\tan\beta$ as a function of $m_{H^{+}}$ in the hMSSM scenario of the MSSM. Limits are shown for $\tan\beta$ values in the range of 0.5-60, where predictions are available from both scenarios. The bands surrounding the expected limits show the 68% and 95% confidence intervals. The limits are based on the combination of the $\ell+$jets and $\ell\ell$ final states. The production cross-section of $t\bar{t}H$ and $tH$, as well as the branching ratios of the $H$, are fixed to their SM values at each point in the plane. Uncertainties on the predicted $H^{+}$ cross-sections or branching ratios are not considered.

Expected and observed lower limits on $\tan\beta$ as a function of $m_{H^{+}}$ in the hMSSM scenario of the MSSM. Limits are shown for $\tan\beta$ values in the range of 0.5-60, where predictions are available from both scenarios. The bands surrounding the expected limits show the 68% and 95% confidence intervals. The limits are based on the combination of the $\ell+$jets and $\ell\ell$ final states. The production cross-section of $t\bar{t}H$ and $tH$, as well as the branching ratios of the $H$, are fixed to their SM values at each point in the plane. Uncertainties on the predicted $H^{+}$ cross-sections or branching ratios are not considered.