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$).

Cutflow for SR-$e\mu$ for 2 representative signal points. For the $\tilde{\tau}$ mass and lifetime points, the number of Monte Carlo events generated are: 30,000 for (200 GeV, 0.1 ns), and 104,000 for (300 GeV, 0.1 ns).

Cutflow for SR-$\mu\mu$ for 5 representative signal points. For the following $\tilde{\mu}$ mass and lifetime points, the number of Monte Carlo events generated are: 24,000 for (100 GeV, 0.01 ns), 16,000 for (300 GeV, 1 ns), and 12,000 for (500 GeV, 0.1 ns). For the $\tilde{\tau}$ mass and lifetime points, the number of Monte Carlo events generated are: 30,000 for (200 GeV, 0.1 ns), and 104,000 for (300 GeV, 0.1 ns).

Cutflow for SR-$\mu\mu$ for 5 representative signal points. For the following $\tilde{\mu}$ mass and lifetime points, the number of Monte Carlo events generated are: 24,000 for (100 GeV, 0.01 ns), 16,000 for (300 GeV, 1 ns), and 12,000 for (500 GeV, 0.1 ns). For the $\tilde{\tau}$ mass and lifetime points, the number of Monte Carlo events generated are: 30,000 for (200 GeV, 0.1 ns), and 104,000 for (300 GeV, 0.1 ns).

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, where all slepton flavors are mass degenerate (co-NLSP).

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, where all slepton flavors are mass degenerate (co-NLSP).

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the selectron signal model. Selectron ($\tilde{e}_{L, R}$) refers to the scalar superpartners of left- and right-handed electrons.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the selectron signal model. Selectron ($\tilde{e}_{L, R}$) refers to the scalar superpartners of left- and right-handed electrons.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the left-handed electrons, $\tilde{e}_L$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the left-handed electrons, $\tilde{e}_L$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the right-handed electrons, $\tilde{e}_R$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the right-handed electrons, $\tilde{e}_R$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the smuon signal model. Smuon ($\tilde{\mu}_{L, R}$) refers to the scalar superpartners of left- and right-handed muons.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the smuon signal model. Smuon ($\tilde{\mu}_{L, R}$) refers to the scalar superpartners of left- and right-handed muons.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the left-handed muons, $\tilde{\mu}_L$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the left-handed muons, $\tilde{\mu}_L$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the right-handed muons, $\tilde{\mu}_R$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the right-handed muons, $\tilde{\mu}_R$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the stau signal model. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin\theta_{\tilde\tau}=0.95$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the stau signal model. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin\theta_{\tilde\tau}=0.95$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for $\tilde{\tau}_L$ production, where $\tilde{\tau}_L$ is pure-state superpartner of the left-handed $\tau$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for $\tilde{\tau}_L$ production, where $\tilde{\tau}_L$ is pure-state superpartner of the left-handed $\tau$.

The expected and observed yields in the signal regions. Combined statistical and systematic uncertainties are presented. Estimates are truncated at 0 if the size of measured systematic uncertainties would yield a negative result.

The expected and observed yields in the signal regions. Combined statistical and systematic uncertainties are presented. Estimates are truncated at 0 if the size of measured systematic uncertainties would yield a negative result.

Reconstruction efficiency as a function of $|d_{0}|$ and $p_\text{T}$ for signal electrons. Monte Carlo samples with $\tilde{e}$ or $\tilde{\mu}$ with mass 400 GeV and 1 ns lifetime were used. Efficiency is defined as the number of leptons passing all signal requirements and matched to generator-level muons divided by the number of generator level leptons with $p_\text{T} > 65$ GeV, $|d_{0}| >$ 3 mm, and |$\eta$| $<$ 2.47 for electrons. No event-level selections are made. Reconstructed leptons must pass all quality criteria, including the cosmic veto. Electron selection scale factors are included on the reconstructed leptons. The overflow is not shown in these plots.

Reconstruction efficiency as a function of $|d_{0}|$ and $p_\text{T}$ for signal electrons. Monte Carlo samples with $\tilde{e}$ or $\tilde{\mu}$ with mass 400 GeV and 1 ns lifetime were used. Efficiency is defined as the number of leptons passing all signal requirements and matched to generator-level muons divided by the number of generator level leptons with $p_\text{T} > 65$ GeV, $|d_{0}| >$ 3 mm, and |$\eta$| $<$ 2.47 for electrons. No event-level selections are made. Reconstructed leptons must pass all quality criteria, including the cosmic veto. Electron selection scale factors are included on the reconstructed leptons. The overflow is not shown in these plots.

Reconstruction efficiency as a function of $|d_{0}|$ and $p_\text{T}$ for signal muons. Monte Carlo samples with $\tilde{e}$ or $\tilde{\mu}$ with mass 400 GeV and 1 ns lifetime were used. Efficiency is defined as the number of leptons passing all signal requirements and matched to generator-level muons divided by the number of generator level leptons with $p_\text{T} > 65$ GeV, $|d_{0}| >$ 3 mm, and |$\eta$| $<$ 2.5 for muons. No event-level selections are made. Reconstructed leptons must pass all quality criteria, including the cosmic veto. Muon selection scale factors are included on the reconstructed leptons. The overflow is not shown in these plots.

Reconstruction efficiency as a function of $|d_{0}|$ and $p_\text{T}$ for signal muons. Monte Carlo samples with $\tilde{e}$ or $\tilde{\mu}$ with mass 400 GeV and 1 ns lifetime were used. Efficiency is defined as the number of leptons passing all signal requirements and matched to generator-level muons divided by the number of generator level leptons with $p_\text{T} > 65$ GeV, $|d_{0}| >$ 3 mm, and |$\eta$| $<$ 2.5 for muons. No event-level selections are made. Reconstructed leptons must pass all quality criteria, including the cosmic veto. Muon selection scale factors are included on the reconstructed leptons. The overflow is not shown in these plots.

Acceptance for $\tilde{e}$ of various masses and lifetimes in SR-$ee$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{e}$ of various masses and lifetimes in SR-$ee$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\mu}$ of various masses and lifetimes in SR-$\mu\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\mu}$ of various masses and lifetimes in SR-$\mu\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Efficiency for $\tilde{e}$ of various masses and lifetimes in SR-$ee$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{e}$ of various masses and lifetimes in SR-$ee$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\mu}$ of various masses and lifetimes in SR-$\mu\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\mu}$ of various masses and lifetimes in SR-$\mu\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$ee$ from $\tilde{e}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47 for electrons, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$ee$ from $\tilde{e}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47 for electrons, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$\mu\mu$ from $\tilde{\mu}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 for muons, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$\mu\mu$ from $\tilde{\mu}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 for muons, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$ee$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$ee$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$e\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$e\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$\mu\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$\mu\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$ee$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$ee$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$e\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$e\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$\mu\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$\mu\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$e\mu$ from $\tilde{\tau}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$e\mu$ from $\tilde{\tau}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane where all slepton flavors and chiralities are mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane where all slepton flavors and chiralities are mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane where all slepton flavors and chiralities are mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane where all slepton flavors and chiralities are mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L, R}$) refer to the scalar superpartners of left- and right-handed electrons, which are assumed to be mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L, R}$) refer to the scalar superpartners of left- and right-handed electrons, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L, R}$) refer to the scalar superpartners of left- and right-handed electrons, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L, R}$) refer to the scalar superpartners of left- and right-handed electrons, which are assumed to be mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L}$) refer to the scalar superpartners of left-handed electrons.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L}$) refer to the scalar superpartners of left-handed electrons.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L}$) refer to the scalar superpartners of left-handed electrons.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L}$) refer to the scalar superpartners of left-handed electrons.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{R}$) refer to the scalar superpartners of right-handed electrons. Purple denotes the region excluded by LEP.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{R}$) refer to the scalar superpartners of right-handed electrons. Purple denotes the region excluded by LEP.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{R}$) refer to the scalar superpartners of right-handed electrons. Purple denotes the region excluded by LEP.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{R}$) refer to the scalar superpartners of right-handed electrons. Purple denotes the region excluded by LEP.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L, R}$) refer to the scalar superpartners of left- and right-handed muons, which are assumed to be mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L, R}$) refer to the scalar superpartners of left- and right-handed muons, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L, R}$) refer to the scalar superpartners of left- and right-handed muons, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L, R}$) refer to the scalar superpartners of left- and right-handed muons, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L}$) refer to the scalar superpartners of left-handed muons.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L}$) refer to the scalar superpartners of left-handed muons.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L}$) refer to the scalar superpartners of left-handed muons.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L}$) refer to the scalar superpartners of left-handed muons.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{R}$) refer to the scalar superpartners of right-handed muons. Purple denotes the region excluded by LEP.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{R}$) refer to the scalar superpartners of right-handed muons. Purple denotes the region excluded by LEP.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{R}$) refer to the scalar superpartners of right-handed muons. Purple denotes the region excluded by LEP.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{R}$) refer to the scalar superpartners of right-handed muons. Purple denotes the region excluded by LEP.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau})$ plane. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin(\theta_{\tilde\tau})=0.95$, which are assumed to be mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau})$ plane. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin(\theta_{\tilde\tau})=0.95$, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau})$ plane. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin(\theta_{\tilde\tau})=0.95$, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau})$ plane. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin(\theta_{\tilde\tau})=0.95$, which are assumed to be mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau}_L)$ plane, where $\tilde{\tau}_L$ is the pure-state super-partner of the left-handed $\tau$. Purple denotes the region excluded by LEP. This result does not present signficant sensitivity to the pure-state superpartner of the right-handed $\tau$.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau}_L)$ plane, where $\tilde{\tau}_L$ is the pure-state super-partner of the left-handed $\tau$. Purple denotes the region excluded by LEP. This result does not present signficant sensitivity to the pure-state superpartner of the right-handed $\tau$.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau}_L)$ plane, where $\tilde{\tau}_L$ is the pure-state super-partner of the left-handed $\tau$. Purple denotes the region excluded by LEP. This result does not present signficant sensitivity to the pure-state superpartner of the right-handed $\tau$.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau}_L)$ plane, where $\tilde{\tau}_L$ is the pure-state super-partner of the left-handed $\tau$. Purple denotes the region excluded by LEP. This result does not present signficant sensitivity to the pure-state superpartner of the right-handed $\tau$.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}(\tilde{\chi})$ scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}(\tilde{\tau})$ scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_h^{125}$(alignment) scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Observed and expected limits on $\tan\beta$ as a function of $m_{H^+}$ in the $M_{h_1}^{125}$(CPV) scenario. Limits are shown for $\tan\beta$ values in the range of 0.5-60 due to the availability of the model prediction. The bands surrounding the expected limits show the 68% and 95% confidence intervals. Uncertainties in the predicted $H^+$ cross-sections or branching ratios are not considered.

Event yields of the SM background processes and the 800 GeV $H^{+}$ sample in the four analysis regions before the fit to the data. Uncertainties include both statistical and systematic uncertainties. The yields of the $H^{+}$ signal sample correspond to a cross-section times branching fraction of 10 pb.

Event acceptance for the different $H^+$ mass signal samples.

Efficiency in the Single-bin SR as a function of the M(Sbottom) vs. M(N2) with the $\Delta$ M(N2,N1) $= 130$ GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Acceptance in the Multi-bin SR, $\min_{\Theta} < 0.5$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the acceptance is given in units of $10^{-4}$.

Efficiency in the Multi-bin SR, $\min_{\Theta} < 0.5$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Acceptance in the Multi-bin SR, $0.5 < \min_{\Theta} < 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the acceptance is given in units of $10^{-4}$.

Efficiency in the Multi-bin SR, $0.5 < \min_{\Theta} < 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Acceptance in the Multi-bin SR, $\min_{\Theta} > 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the acceptance is given in units of $10^{-4}$.

Efficiency in the Multi-bin SR, $\min_{\Theta} > 1.0$ bin as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV. Keep in mind that the efficiency is given in units of $10^{-2}$.

Observed upper limits on the signal cross section as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV.

Expected upper limits on the signal cross section as a function of the M(Sbottom) vs. M(N2) with the $\Delta M$(N2,N1) = 130 GeV.

Cutflows for the bechmarl signal point M(Sbottom) = 800 GeV, M(N2) = 180 GeV. Weighted event yields are reported starting with the "Preselection" line, normalized to an integrated luminosity of $139$ fb$^{−1}$.

Comparison of the expected and observed event yields in the signal regions. The top-quark and Z(mumu) background contributions are scaled with the normalization factors obtained from the background-only fit. The other contribution includes all the backgrounds not explicitly listed in the legend (V+jets except Z(mumu)+jets, di-/triboson, multijet). The hatched band indicates the total statistical and systematic uncertainties in the SM background. The contributions from three signal models to the signal regions are also displayed, where the masses M(Sbottom) and M(N2) are given in GeV in the legend. The lower panel shows the significance of the deviation of the observed yield from the expected background yield.

Dominant systematic uncertainties in the background prediction for the signal regions after the fit to the control regions. “Other” includes the uncertainties arising from muons, jet-vertex tagging, modeling of pile-up, the $E_{T}^{miss}$ computation, multijet background, and luminosity. The individual uncertainties can be correlated and do not necessarily add up quadratically to the total uncertainty.

Comparison between data and background prediction before the fit (Pre-Fit) for the mass of the FCNC top-quark candidate in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction before the fit (Pre-Fit) for the mass of the SM top-quark candidate in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction before the fit (Pre-Fit) for the transverse momentum of the Z boson candidate in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the mass sideband CR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{2}^{u}$ discriminant in the mass sideband CR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{2}^{u}$ discriminant in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the mass sideband CR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{2}^{c}$ discriminant in the mass sideband CR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{1}$ discriminant in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{2}^{c}$ discriminant in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the leading lepton $p_{T}$ in the $t\bar{t}$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the third lepton $p_{T}$ in the $t\bar{t}$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the leading lepton $p_{T}$ in the $t\bar{t}Z$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the $t\bar{t}Z$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Charged particle multiplicity distribution.

Average charged particle transverse momentum as a function of the number of charged particles in the event.

The unfolded differential charge asymmetry as a function of the transverse momentum of the top pair system. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NNLO in QCD and NLO in EW theory are listed. The scale uncertainty is obtained by varying renormalisation and factorisation scales independently by a factor of 2 or 0.5 around $\mu_0$ to calculate the maximum and minimum value of the asymmetry, respectively. The nominal value $\mu_0$ is chosen as $H_T/4$. The variations in which one scale is multiplied by 2 while the other scale is divided by 2 are excluded. Finally, the scale and MC integration uncertainties are added in quadrature.

The unfolded differential charge asymmetry as a function of the longitudinal boost of the top pair system. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NNLO in QCD and NLO in EW theory are listed. The scale uncertainty is obtained by varying renormalisation and factorisation scales independently by a factor of 2 or 0.5 around $\mu_0$ to calculate the maximum and minimum value of the asymmetry, respectively. The nominal value $\mu_0$ is chosen as $H_T/4$. The variations in which one scale is multiplied by 2 while the other scale is divided by 2 are excluded. Finally, the scale and MC integration uncertainties are added in quadrature.

The unfolded inclusive leptonic asymmetry. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the invariant mass of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the transverse momentum of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the longitudinal boost of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

Individual 68% and 95% CL bounds on the relevant Wilson coefficients of the SM Effective Field Theory in units of $\text{TeV}^{-2}$. The bounds are derived from the $A_C^{t\bar{t}}$ inclusive measurement. The experimental uncertainties are accounted for, in the form of the complete covariance matrix that keeps track of correlations between bins for the differential measurement. The theory uncertainty from the NNLO QCD + NLO EW calculation is included by explicitly varying the renormalization and factorization scales, or the parton density functions, in the calculation and registering the variations in the intervals.

Individual 68% and 95% CL bounds on the relevant Wilson coefficients of the SM Effective Field Theory in units of $\text{TeV}^{-2}$. The bounds are derived from the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. The experimental uncertainties are accounted for, in the form of the complete covariance matrix that keeps track of correlations between bins for the differential measurement. The theory uncertainty from the NNLO QCD + NLO EW calculation is included by explicitly varying the renormalization and factorization scales, or the parton density functions, in the calculation and registering the variations in the intervals.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ inclusive measurement. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0,0.3]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.3,0.6]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.6,0.8]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.8,1]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ < 500 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [500,750] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [750,1000] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [1000,1500] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ > 1500 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ < 30 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ $\in$ [30,120] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ > 120 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ inclusive measurement. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0,0.3]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.3,0.6]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.6,0.8]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.8,1]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ < 200 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ $\in$ [200,300] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ $\in$ [300,400] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ > 400 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ < 20 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ $\in$ [20, 70] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ > 70 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ inclusive measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ inclusive measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement.

Covariance matrix for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement.

Covariance matrix for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement.

Covariance matrix for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement.

Probability of finding at least one TAR jet, where the p_T-leading TAR jet passes the m_Wcand and D₂^&beta;=1 requirements, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=500 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1000 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1700 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied that do not pass the requirements of the merged category.

Observed exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=0.5 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1.7 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1.7 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=2.1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=2.1 TeV, shown in dashed blue.

Data overlaid on SM background yields stacked in each SR and CR category after the fit to data ('Post-fit'). The yields in the SR are broken down into their contributions to the individual bins. The maximum-likelihood estimators are set to the conditional values of the CR-only fit, and propagated to SR and CRs.

Dominant sources of uncertainty for three dark Higgs scenarios after the fit to data. The uncertainties are quantified in terms of their contribution to the fitted signal uncertainty that is expressed relative to the theory prediction. Three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01 and m_&chi;=200 GeV are considered, which are indicated using the (m_Z', m_s) format in units of GeV in the table columns.

Cumulative efficiencies in the merged category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Cumulative efficiencies in the resolved category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={300, 350, 400, 500, 750} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={1000, 1700} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={2100, 2500, 2900, 3300} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Measured fiducial cross section times leptonic branching ratio for Z/gamma* production in the Z/gamma* -> e+e- final state.

Measured fiducial cross section times leptonic branching ratio for W+ production in the W+ -> mu+ nu final state.

Measured fiducial cross section times leptonic branching ratio for W- production in the W- -> mu- nubar final state.

Measured fiducial cross section times leptonic branching ratio for W+/- production in the combined W+ -> mu+ nu and W- -> mu- nubar final state.

Measured fiducial cross section times leptonic branching ratio for Z/gamma* production in the Z/gamma* -> mu+ mu- final state.

Measured total cross section times leptonic branching ratio for W+ production in the W+ -> e+ nu final state.

Measured total cross section times leptonic branching ratio for W- production in the W- -> e- nubar final state.

Measured total cross section times leptonic branching ratio for W+/- production in the combined W+ -> e+ nu and W- -> e- nubar final state.

Measured total cross section times leptonic branching ratio for Z/gamma* production in the Z/gamma* -> e+ e- final state.

Measured total cross section times leptonic branching ratio for W+ production in the W+ -> mu+ nu final state.

Measured total cross section times leptonic branching ratio for W- production in the W- -> mu- nubar final state.

Measured total cross section times leptonic branching ratio for W+/- production in the combined W+ -> mu+ nu and W- -> mu- nubar final state.

Measured total cross section times leptonic branching ratio for Z/gamma* production in the Z/gamma* -> mu+ mu- final state.

Measured total cross section times leptonic branching ratio for W+ production in the W+ -> l+ nu (l = e, mu) final states.

Measured total cross section times leptonic branching ratio for W- production in the W- -> l- nubar (l = e, mu) final states.

Measured total cross section times leptonic branching ratio for W+/- production in the combined W+ -> l+ nu and W- -> l- nubar (l = e, mu) final states.

Measured total cross section times leptonic branching ratio for Z/gamma* production in the combined Z/gamma* -> l+ l- (l = e, mu) final states.

Measured total cross-section ratio R_{W+/Z} = sigma (W+ -> e+ nu) / sigma (Z/gamma^* -> e+ e-).

Measured total cross-section ratio R_{W-/Z} = sigma (W- -> e- nubar) / sigma (Z/gamma^* -> e+ e-).

Measured total cross-section ratio R_{W+-/Z} = sigma (W+/- -> e+/- nu/nubar) / sigma (Z/gamma^* -> e+ e-).

Measured total cross-section ratio R_{W+/Z} = sigma (W+ -> mu+ nu) / sigma (Z/gamma^* -> mu+ mu-).

Measured total cross-section ratio R_{W-/Z} = sigma (W- -> mu- nubar) / sigma (Z/gamma^* -> mu+ mu-).

Measured total cross-section ratio R_{W+-/Z} = sigma (W+/- -> mu+/- nu/nubar) / sigma (Z/gamma^* -> mu+ mu-).

Measured total cross-section ratio R_{W+/Z} = sigma (W+ -> l+ nu) / sigma (Z/gamma^* -> l+ l-) where (l = e, mu).

Measured total cross-section ratio R_{W-/Z} = sigma (W- -> l- nubar) / sigma (Z/gamma^* -> l+ l-) where (l = e, mu).

Measured total cross-section ratio R_{W+-/Z} = sigma (W+/- -> l+/- nu/nubar) / sigma (Z/gamma^* -> l+ l-) where (l = e, mu).

Measured lepton asymmetry from W+ -> e+ nu and W- -> e- nubar, binned as a function pseudorapidity.

Measured lepton asymmetry from W+ -> e+ nu and W- -> e- nubar.

Measured lepton asymmetry from W+ -> mu+ nu and W- -> mu- nubar, binned as a function pseudorapidity.

Measured lepton asymmetry from W+ -> mu+ nu and W- -> mu- nubar.

Measured lepton asymmetry from W+ -> l+ nu and W- -> l- nubar (l = e, mu), binned as a function pseudorapidity.

Measured lepton asymmetry from W+ -> l+ nu and W- -> l- nubar (l = e, mu).

CHI distribution for mass bin > 1200 GeV.

Centrality Ratio.

Inclusive jet double-differential cross sections in the |rapidity| range 1.2 to 2.1, using a jet resolution R value of 0.4. The three (sys) errors are respectively, the Absolute JES, the Unfolding and the Luminosity uncertainties.

Inclusive jet double-differential cross sections in the |rapidity| range 2.1 to 2.8, using a jet resolution R value of 0.4. The three (sys) errors are respectively, the Absolute JES, the Unfolding and the Luminosity uncertainties.

Inclusive jet double-differential cross sections in the |rapidity| range 0 to 0.3, using a jet resolution R value of 0.6. The three (sys) errors are respectively, the Absolute JES, the Unfolding and the Luminosity uncertainties.

Inclusive jet double-differential cross sections in the |rapidity| range 0.3 to 0.8, using a jet resolution R value of 0.6. The three (sys) errors are respectively, the Absolute JES, the Unfolding and the Luminosity uncertainties.

Inclusive jet double-differential cross sections in the |rapidity| range 0.8 to 1.2, using a jet resolution R value of 0.6. The three (sys) errors are respectively, the Absolute JES, the Unfolding and the Luminosity uncertainties.

Inclusive jet double-differential cross sections in the |rapidity| range 1.2 to 2.1, using a jet resolution R value of 0.6. The three (sys) errors are respectively, the Absolute JES, the Unfolding and the Luminosity uncertainties.

Inclusive jet double-differential cross sections in the |rapidity| range 2.1 to 2.8, using a jet resolution R value of 0.6. The three (sys) errors are respectively, the Absolute JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 0 to 0.3, using a jet resolution R value of 0.4. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 0.3 to 0.8, using a jet resolution R value of 0.4. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 0.8 to 1.2, using a jet resolution R value of 0.4. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 1.2 to 2.1, using a jet resolution R value of 0.4. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 2.1 to 2.8, using a jet resolution R value of 0.4. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 0 to 0.3, using a jet resolution R value of 0.6. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 0.3 to 0.8, using a jet resolution R value of 0.6. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 0.8 to 1.2, using a jet resolution R value of 0.6. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 1.2 to 2.1, using a jet resolution R value of 0.6. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 2.1 to 2.8, using a jet resolution R value of 0.6. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 340 to 520, using a jet resolution R value of 0.4. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 520 to 800, using a jet resolution R value of 0.4. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 800 to 1200, using a jet resolution R value of 0.4. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 340 to 520, using a jet resolution R value of 0.6. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 520 to 800, using a jet resolution R value of 0.6. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

Dijet double-differential cross sections in the |rapidity(max)| range 800 to 1200, using a jet resolution R value of 0.6. The four (sys) errors are respectively, the Absolute JES, the Relative JES, the Unfolding and the Luminosity uncertainties.

The performance of jet energy resolution (JER) for jets with |y| < 2.1 evaluated as a function of pT_truth in different centrality bins. Simulated hard scatter events were overlaid onto events from a dedicated sample of minimum-bias Xe+Xe data. The fit parameters are listed in a sperate table (Extras 1)

The relative magnitude of systematic uncertainties for per-pair normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for absolutely normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for rho distributions for leading jets in 0-10% Xe+Xe centrality

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

Parameter a,b, and c from JER fits in Figure 1b.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 0 to 2.8.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 0 to 2.8.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 0 to 2.8.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 0 to 2.8.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 0 to 2.8.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 310 to 400 GeV and absolute values of the jet rapidity from 0 to 2.8.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 400 to 500 GeV and absolute values of the jet rapidity from 0 to 2.8.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 500 to 600 GeV and absolute values of the jet rapidity from 0 to 2.8.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 30 to 40 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 30 to 40 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 30 to 40 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 30 to 40 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 30 to 40 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 40 to 60 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 40 to 60 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 40 to 60 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 40 to 60 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 40 to 60 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 60 to 80 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 60 to 80 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 60 to 80 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 60 to 80 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 60 to 80 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 310 to 400 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 310 to 400 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 310 to 400 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 310 to 400 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 400 to 500 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 400 to 500 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 400 to 500 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Differential Jet Shape RHO as a function of r for jet transverse momentum from 400 to 500 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Integrated Jet Shape 1-PSI as a function of jet transverse momentum for an r value of 0.3 and absolute values of the jet rapidity from 0 to 2.8.

Measured Integrated Jet Shape 1-PSI as a function of jet transverse momentum for an r value of 0.3 and absolute values of the jet rapidity from 0 to 0.3.

Measured Integrated Jet Shape 1-PSI as a function of jet transverse momentum for an r value of 0.3 and absolute values of the jet rapidity from 0.3 to 0.8.

Measured Integrated Jet Shape 1-PSI as a function of jet transverse momentum for an r value of 0.3 and absolute values of the jet rapidity from 0.8 to 1.2.

Measured Integrated Jet Shape 1-PSI as a function of jet transverse momentum for an r value of 0.3 and absolute values of the jet rapidity from 1.2 to 2.1.

Measured Integrated Jet Shape 1-PSI as a function of jet transverse momentum for an r value of 0.3 and absolute values of the jet rapidity from 2.1 to 2.8.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 30 to 40 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 30 to 40 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 30 to 40 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 30 to 40 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 30 to 40 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 30 to 40 GeV and absolute values of the jet rapidity from 0 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 40 to 60 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 40 to 60 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 40 to 60 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 40 to 60 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 40 to 60 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 40 to 60 GeV and absolute values of the jet rapidity from 0 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 60 to 80 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 60 to 80 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 60 to 80 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 60 to 80 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 60 to 80 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 60 to 80 GeV and absolute values of the jet rapidity from 0 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 80 to 110 GeV and absolute values of the jet rapidity from 0 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 110 to 160 GeV and absolute values of the jet rapidity from 0 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 160 to 210 GeV and absolute values of the jet rapidity from 0 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 210 to 260 GeV and absolute values of the jet rapidity from 0 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 2.1 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 260 to 310 GeV and absolute values of the jet rapidity from 0 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 310 to 400 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 310 to 400 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 310 to 400 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 310 to 400 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 310 to 400 GeV and absolute values of the jet rapidity from 0 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 400 to 500 GeV and absolute values of the jet rapidity from 0 to 0.3. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 400 to 500 GeV and absolute values of the jet rapidity from 0.3 to 0.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 400 to 500 GeV and absolute values of the jet rapidity from 0.8 to 1.2. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 400 to 500 GeV and absolute values of the jet rapidity from 1.2 to 2.1. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 400 to 500 GeV and absolute values of the jet rapidity from 0 to 2.8. This is additional data, not in the paper.

Measured Integrated Jet Shape PSI as a function of r for jet transverse momentum from 500 to 600 GeV and absolute values of the jet rapidity from 0 to 2.8. This is additional data, not in the paper.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of pseudorapidity for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of pseudorapidity for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of pseudorapidity for events with the number of charged particles >=6 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of pseudorapidity for events with the number of charged particles >=6 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 2360 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of transverse momentum for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of transverse momentum for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of transverse momentum for events with the number of charged particles >=6 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of transverse momentum for events with the number of charged particles >=6 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 900 GeV for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 2360 GeV for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 7000 GeV for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 900 GeV for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 7000 GeV for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 900 GeV for events with the number of charged particles >=6 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 7000 GeV for events with the number of charged particles >=6 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Average transverse momentum in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Average transverse momentum in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Average transverse momentum in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Average transverse momentum in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=2 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of pseudorapidity for events with the number of charged particles >=20 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of pseudorapidity for events with the number of charged particles >=20 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of pseudorapidity for events with the number of charged particles >=1 having transverse momentum >2500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of pseudorapidity for events with the number of charged particles >=1 having transverse momentum >2500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of transverse momentum for events with the number of charged particles >=20 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of transverse momentum for events with the number of charged particles >=20 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >2500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >2500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 900 GeV for events with the number of charged particles >=20 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 7000 GeV for events with the number of charged particles >=20 having transverse momentum >100 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 900 GeV for events with the number of charged particles >=1 having transverse momentum >2500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distributions in proton-proton collisions at a centre-of mass energy of 7000 GeV for events with the number of charged particles >=1 having transverse momentum >2500 MeV and absolute(pseudorapidity) <2.5.

Average transverse momentum in proton-proton collisions at a centre-of mass energy of 900 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >2500 MeV and absolute(pseudorapidity) <2.5.

Average transverse momentum in proton-proton collisions at a centre-of mass energy of 7000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >2500 MeV and absolute(pseudorapidity) <2.5.

The average charged-particle muliplicity per unit of rapidity for ETARAP=0 as a function of the centre-of-mass energy.

The average charged-particle muliplicity per unit of rapidity in the pseudorapidity region -2.5 to 2.5 for events with 2 or more charged particles as a function of the centre-of-mass energy.

The best-fit predicted and observed distributions of the combined NN output in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the high-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

Observed and expected upper limits on $\mathcal{B}(Z\rightarrow\ell\tau)$ with leptonically decaying $\tau$ at 95% confidence level.

Observed and expected upper limits on $\mathcal{B}(Z\rightarrow\ell\tau)$ at 95% confidence level, combination of hadronically and leptonically decaying $\tau$ channels.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the CRZ$\tau\tau$ for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the CRZ$\tau\tau$ for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the combined NN output in the CRZ$\tau\tau$ for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the combined NN output in the CRZ$\tau\tau$ for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the leading lepton transverse momentum $p_\text{T}(e)$ in the high-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the leading lepton transverse momentum $p_\text{T}(\mu)$ in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the missing transverse momentum $E^\text{miss}_\text{T}$ in the low-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the missing transverse momentum $E^\text{miss}_\text{T}$ in the low-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

IRing2 for HT2>=500 GeV, NJets>=5

IRing2 for HT2>=1000 GeV, NJets>=2

IRing2 for HT2>=1000 GeV, NJets>=3

IRing2 for HT2>=1000 GeV, NJets>=4

IRing2 for HT2>=1000 GeV, NJets>=5

IRing2 for HT2>=1500 GeV, NJets>=2

IRing2 for HT2>=1500 GeV, NJets>=3

IRing2 for HT2>=1500 GeV, NJets>=4

IRing2 for HT2>=1500 GeV, NJets>=5

IRing128 for HT2>=500 GeV, NJets>=2

IRing128 for HT2>=500 GeV, NJets>=3

IRing128 for HT2>=500 GeV, NJets>=4

IRing128 for HT2>=500 GeV, NJets>=5

IRing128 for HT2>=1000 GeV, NJets>=2

IRing128 for HT2>=1000 GeV, NJets>=3

IRing128 for HT2>=1000 GeV, NJets>=4

IRing128 for HT2>=1000 GeV, NJets>=5

IRing128 for HT2>=1500 GeV, NJets>=2

IRing128 for HT2>=1500 GeV, NJets>=3

IRing128 for HT2>=1500 GeV, NJets>=4

IRing128 for HT2>=1500 GeV, NJets>=5

ICyl16 for HT2>=500 GeV, NJets>=2

ICyl16 for HT2>=500 GeV, NJets>=3

ICyl16 for HT2>=500 GeV, NJets>=4

ICyl16 for HT2>=500 GeV, NJets>=5

ICyl16 for HT2>=1000 GeV, NJets>=2

ICyl16 for HT2>=1000 GeV, NJets>=3

ICyl16 for HT2>=1000 GeV, NJets>=4

ICyl16 for HT2>=1000 GeV, NJets>=5

ICyl16 for HT2>=1500 GeV, NJets>=2

ICyl16 for HT2>=1500 GeV, NJets>=3

ICyl16 for HT2>=1500 GeV, NJets>=4

ICyl16 for HT2>=1500 GeV, NJets>=5

IRing2 covariance for HT2>=500 GeV, NJets>=2 (Table 1)

IRing2 covariance for HT2>=500 GeV, NJets>=3 (Table 2)

IRing2 covariance for HT2>=500 GeV, NJets>=4 (Table 3)

IRing2 covariance for HT2>=500 GeV, NJets>=5 (Table 4)

IRing2 covariance for HT2>=1000 GeV, NJets>=2 (Table 5)

IRing2 covariance for HT2>=1000 GeV, NJets>=3 (Table 6)

IRing2 covariance for HT2>=1000 GeV, NJets>=4 (Table 7)

IRing2 covariance for HT2>=1000 GeV, NJets>=5 (Table 8)

IRing2 covariance for HT2>=1500 GeV, NJets>=2 (Table 9)

IRing2 covariance for HT2>=1500 GeV, NJets>=3 (Table 10)

IRing2 covariance for HT2>=1500 GeV, NJets>=4 (Table 11)

IRing2 covariance for HT2>=1500 GeV, NJets>=5 (Table 12)

IRing128 covariance for HT2>=500 GeV, NJets>=2 (Table 13)

IRing128 covariance for HT2>=500 GeV, NJets>=3 (Table 14)

IRing128 covariance for HT2>=500 GeV, NJets>=4 (Table 15)

IRing128 covariance for HT2>=500 GeV, NJets>=5 (Table 16)

IRing128 covariance for HT2>=1000 GeV, NJets>=2 (Table 17)

IRing128 covariance for HT2>=1000 GeV, NJets>=3 (Table 18)

IRing128 covariance for HT2>=1000 GeV, NJets>=4 (Table 19)

IRing128 covariance for HT2>=1000 GeV, NJets>=5 (Table 20)

IRing128 covariance for HT2>=1500 GeV, NJets>=2 (Table 21)

IRing128 covariance for HT2>=1500 GeV, NJets>=3 (Table 22)

IRing128 covariance for HT2>=1500 GeV, NJets>=4 (Table 23)

IRing128 covariance for HT2>=1500 GeV, NJets>=5 (Table 24)

ICyl16 covariance for HT2>=500 GeV, NJets>=2 (Table 25)

ICyl16 covariance for HT2>=500 GeV, NJets>=3 (Table 26)

ICyl16 covariance for HT2>=500 GeV, NJets>=4 (Table 27)

ICyl16 covariance for HT2>=500 GeV, NJets>=5 (Table 28)

ICyl16 covariance for HT2>=1000 GeV, NJets>=2 (Table 29)

ICyl16 covariance for HT2>=1000 GeV, NJets>=3 (Table 30)

ICyl16 covariance for HT2>=1000 GeV, NJets>=4 (Table 31)

ICyl16 covariance for HT2>=1000 GeV, NJets>=5 (Table 32)

ICyl16 covariance for HT2>=1500 GeV, NJets>=2 (Table 33)

ICyl16 covariance for HT2>=1500 GeV, NJets>=3 (Table 34)

ICyl16 covariance for HT2>=1500 GeV, NJets>=4 (Table 35)

ICyl16 covariance for HT2>=1500 GeV, NJets>=5 (Table 36)

IRing2 covariance, complete

1-IRing128 covariance, complete

1-ICyl16 covariance, complete

Expected 95% CL upper limits with 1-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 300 GeV.

Expected 95% CL upper limits with 2-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 300 GeV.

Observed 95% CL upper limits on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 300 GeV.

Expected 95% CL upper limits on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 600 GeV.

Expected 95% CL upper limits with 1-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 600 GeV.

Expected 95% CL upper limits with 2-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 600 GeV.

Observed 95% CL upper limits on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 600 GeV.

Expected 95% CL upper limits on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 900 GeV.

Expected 95% CL upper limits with 1-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 900 GeV.

Expected 95% CL upper limits with 2-sigma uncertainty on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 900 GeV.

Observed 95% CL upper limits on the production of a heavy Higgs boson as functions of fw and fww with mass equal to 900 GeV.

Observed (black solid curve) and expected (black dashed curve) 95% CL upper limits on the production of a heavy Higgs boson as functions of its mass with ( fW, fWW) fixed at (0, 6200). The green (inner) and yellow (outer) bands represent 1 and 2 sigma uncertainty in the expected limits. Theoretical predictions (red solid curve) as functions of the heavy Higgs boson mass are overlaid.

Observed (black solid curve) and expected (black dashed curve) 95% CL upper limits on the production of a heavy Higgs boson as functions of its mass with ( fW, fWW) fixed at (1350, 0). The green (inner) and yellow (outer) bands represent 1 and 2 sigma uncertainty in the expected limits. Theoretical predictions (red solid curve) as functions of the heavy Higgs boson mass are overlaid.

Summary of the total uncertainty in the background prediction for each SR of the tt0L-low, tt0L-high, tt1L and tt2L analysis channels in the statistical combination. Their dominant contributions are indicated by individual lines. Individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

$E_{\text{T}}^{\text{miss}}$ distribution in SR0X for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRWX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRTX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Representative fit distribution in the different flavour leptons signal region for the tt2L analysis: each bin of such distribution, starting from the red arrow, corresponds to a single SR included in the fit. 'FNP' includes the contribution from fake/non-prompt lepton background arising from jets (mainly $\pi/K$, heavy-flavour hadron decays and photon conversion) misidentified as leptons, estimated in a purely data-driven way. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The total uncertainty in the SM expectation is represented with hatched bands and the expected distributions for selected signal models are shown as dashed lines.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of selected reconstructed events divided by the acceptance.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Correlation matrix for the measurement of STXS parameters in each of the 28 regions considered.

Fitted values for kappa_g and kappa_y.

Correlation matrix for the measurement of kappa_g and kappa_y.

Summary of the 68% CL confidence intervals for individual measurements of SMEFT parameters observed in data. In each case, SMEFT parameters other than the one measured are fixed to 0.

Results of the EV parameter measurement in data, in the linear and linear+quadratic parameterizations of the SMEFT. All the EVs parameters are free to vary in the fits. The ranges correspond to 68% CL confidence intervals.

Observed linear correlation coefficients of the EVs parameters in the linear parameterization.

Observed linear correlation coefficients of the EVs parameters in the linear+quadratic parameterization.

Cross-sections times H->yy branching ratio for ggF +bbH, VBF, VH, ttH, and tH production. The values are obtained from a simultaneous fit to all categories. The theory uncertainties in the predictions include uncertainties due to missing higher-order terms in the perturbative QCD calculations and choices of parton distribution functions and value of alpha_s, as well as the H->yy branching ratio uncertainty.

Best-fit values and uncertainties for STXS parameters in each of the 28 regions considered. The values for the gg->H process also include the contributions from bbH production.

Correlation matrix for the measurement of STXS parameters in each of the 33 regions considered.

Cumulative efficiency for simulated signal events to pass each stage of the event selection in the hadhad category. The efficiencies are calculated with respect to HH to bbtautau decays in which both tau-leptons decay hadronically. The ‘Pre-selection’ consists of basic requirements, including that at least two tau-had-vis pass loose kinematic requirements, at least one of the tau-had-vis candidate passes identification requirement, and that events do not contain an electron or muon. The ‘Object selections’ require exactly two tau-had-vis candidates, and at least two jets with pT > 25 GeV and abs(eta) < 2.5. The ‘Loose tau ID’ requires that both tau-had-vis candidates pass identification requirements. The ‘DTT offline jet cuts’ are cuts placed on the pT of the reconstructed jet or tau-had-vis that are geometrically matched to the HLT objects, to ensure the efficiencies of the HLT objects reach the plateau region.

Cumulative efficiency for simulated signal events to pass each stage of the event selection in the lephad SLT category. The efficiencies are calculated with respect to HH to bbtautau decays in which one tau-lepton decays hadronically and one decays leptonically. The ‘Pre-selection’ consists of basic requirements, including that at least one tau-had-vis candidate and one lepton pass loose kinematic requirements. The ‘Object selections’ require exactly one tau-had-vis candidate, and at least two jets with pT > 25 GeV and abs(eta) < 2.5. The ‘Trigger specific offline pT cuts’ are cuts placed on the pT of the reconstructed jet or tau-had-vis that are geometrically matched to the HLT objects, to ensure the efficiencies of the HLT objects reach the plateau region.

Cumulative efficiency for simulated signal events to pass each stage of the event selection in the lephad LTT category. The efficiencies are calculated with respect to HH to bbtautau decays in which one tau-lepton decays hadronically and one decays leptonically. The ‘Pre-selection’ consists of basic requirements, including that at least one tau-had-vis candidate and one lepton pass loose kinematic requirements. The ‘Object selections’ require exactly one tau-had-vis candidate, and at least two jets with pT > 25 GeV and abs(eta) < 2.5. The ‘Trigger specific offline pT cuts’ are cuts placed on the pT of the reconstructed jet or tau-had-vis that are geometrically matched to the HLT objects, to ensure the efficiencies of the HLT objects reach the plateau region.

Post-fit expected number of signal and background events and observed number of data events after applying the selection criteria and requiring exactly 2 b-tagged jets and assuming a background-only hypothesis.

Observed and expected limits at 95% CL on the cross-section of HH production, for the non-resonant ggF+VBF HH search, and the resonant HH search for four values of the resonance mass mX.

Acceptance times efficiency for the full analysis selections as a function of the resonance mass mX in the hadhad, lephad SLT and lephad LTT trigger categories, and the combined lephad. The acceptance times efficiency is evaluated for HH to bbtautau decays, with respect to the targeted tau-lepton decay modes (lephad or hadhad).

Post-fit distribution of mHH in the hadhad channel.

Post-fit distribution of mHH in the lephad SLT channel.

Post-fit distribution of mHH in the lephad LTT channel.

Post-fit distribution of mtautau MMC in the hadhad channel.

Post-fit distribution of mtautau MMC in the hadhad channel.

Post-fit distribution of mtautau MMC in the hadhad channel.

Post-fit distribution of of di-b-jet mass in the hadhad channel.

Post-fit distribution of di-b-jet mass in the LTT lephad channel.

Post-fit distribution of di-b-jet mass in the LTT lephad channel.

BDT for SM HH in the hadhad channel

NN for SM HH in the lephad SLT channel

NN for SM HH in the lephad LTT channel

PNN for mX = 500 GeV resonant HH in the hadhad channel

PNN for mX = 500 GeV resonant HH in the lephad SLT channel

PNN for mX = 500 GeV resonant HH in the lephad LTT channel

PNN for mX = 1000 GeV resonant HH in the hadhad channel

PNN for mX = 1000 GeV resonant HH in the lephad SLT channel

PNN for mX = 1000 GeV resonant HH in the lephad LTT channel

Event yields as a function of log10(S/B) for data, background and non-resonant HH signal. Final discriminant bins from the hadhad, lephad SLT and lephad LTT categories are combined into bins of log10(S/B). The B is the fitted background yield assuming background-only hypothesis, and the signal S is scaled to the SM expected cross-section.

Observed and expected limits at 95% CL on the cross-section of the resonant HH production as a function of the scalar resonance mass mX.

Post-fit distribution of delta-R between the taus in the hadhad channel.

Post-fit distribution of delta-R between the b-tagged jets in the hadhad channel.

Post-fit distribution of delta-R between the taus in the lephad SLT channel.

Post-fit distribution of delta-R between the b-tagged jets in the lephad SLT channel.

Post-fit distribution of delta-pT between the tau and lepton in the lephad SLT channel.

Post-fit distribution of pT of the subleading b-tagged jet in the lephad SLT channel.

Post-fit distribution of MTW in the lephad SLT channel.

Post-fit distribution of missing transverse momentum in the lephad SLT channel.

Post-fit distribution of missing transverse momentum centrality in the lephad SLT channel.

Post-fit distribution of delta-phi between the Higgs boson candidates in the lephad SLT channel.

Post-fit distribution of delta-pT between the tau and lepton in the lephad LTT channel.

Post-fit distribution of delta-R between the taus in the lephad LTT channel.

Post-fit distribution of delta-phi between the lepton and the missing transverse momentum in the lephad LTT channel.

Post-fit distribution of delta-phi between the Higgs boson candidates in the lephad LTT channel.

Post-fit distribution of the total transverse momentum s in the lephad LTT channel.

PNN for mX = 300 GeV resonant HH in the hadhad channel

PNN for mX = 300 GeV resonant HH in the lephad SLT channel

PNN for mX = 300 GeV resonant HH in the lephad LTT channel

PNN for mX = 1600 GeV resonant HH in the hadhad channel

PNN for mX = 1600 GeV resonant HH in the lephad SLT channel

PNN for mX = 1600 GeV resonant HH in the lephad LTT channel

Local p-value of the background-only hypothesis as a function of the resonance mass.