Event yields for a simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ sample with $m_a = 17.5\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. Cut 0 corresponds to the initial number of events. Cut 1 requires the single lepton trigger. Cut 2 requires 2 identified leptons. Cut 3 requires the Z-boson mass window. Cut 4 requires 2 reconstructed $a\to b\bar{b}$ candidates. Cut 5a requires 2 identified $a\to b\bar{b}$ candidates in the 1HPC1LPC region. Cut 6a requires the 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region to be inside the Higgs mass window. Cut 5b requires 2 identified $a\to b\bar{b}$ candidates in the 2HPC region. Cut 6b requires the 2 $a\to b\bar{b}$ candidates in the 2HPC region to be inside the Higgs mass window.

Background yield table for Z+jets, $t\bar{t}$, and rare sources. Observed data yield. Signal $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ yield with $m_a = 20\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$, with a branching ratio set to 1 for the $H \rightarrow aa$ decay, whereas the ATLAS figure attached to this entry instead uses the upper-limit branching ratio (smaller than 1). The table includes the yields in two signal regions with leptons consistent with an on-shell Z-boson decay, one with 2 $a\to b\bar{b}$ candidates in the 2HPC region and one with 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region. The table also includes the yields in four control regions, one with leptons consistent with an on-shell Z-boson decay and 2 $a\to b\bar{b}$ candidates in the Low Purity Category (LPC), and three others where the leptons are not consistent an on-shell Z-boson decay.

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the Higgs boson mass.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

Acceptance times efficiency for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance times efficiency for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Observed 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Expected 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

Distributions of $m_{T2}$ in VR-top-low for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Distributions of $m_{T2}$ in VR-top-high for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Distributions of $E_T^{miss}$ in VR-WW-0J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Distributions of $E_T^{miss}$ in VR-WW-1J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Distributions of $E_T^{miss}$ significance in VR-VZ for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Distributions of $E_T^{miss}$ significance in VR-top-WW for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow.

Breakdown of the dominant systematic uncertainties on background estimates in the inclusive SRs requiring $m_{T2}$>100 GeV after performing the profile likelihood fit. Note that the individual uncertainties can be correlated, and do not necessarily add up quadratically to the total background uncertainty. The percentages show the size of the uncertainty relative to the total expected background. "Top theoretical uncertainties" refers to $t\bar t$ theoretical uncertainties and the uncertainty associated to $Wt-t\bar t$ interference added quadratically.

Observed events and predicted post-fit background yields for the DF inclusive SRs. The model independent upper limits at 95% confidence level (CL) on the observed and expected number of beyond the SM events $S^{0.95}_{obs/exp}$ and the effective beyond the SM cross-section $\sigma^{0.95}_{obs}$ are also reported. The last row reports the $p_0$-value of the SM-only hypothesis. For SRs where the data yield is smaller than expected, the $p$-value is truncated at 0.50. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$+V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted post-fit background yields for the DF inclusive SRs. The model independent upper limits at 95% confidence level (CL) on the observed and expected number of beyond the SM events $S^{0.95}_{obs/exp}$ and the effective beyond the SM cross-section $\sigma^{0.95}_{obs}$ are also reported. The last row reports the $p_0$-value of the SM-only hypothesis. For SRs where the data yield is smaller than expected, the $p$-value is truncated at 0.50. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$+V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted post-fit background yields for the SF inclusive SRs. The model independent upper limits at 95% confidence level (CL) on the observed and expected number of beyond the SM events $S^{0.95}_{obs/exp}$ and the effective beyond the SM cross-section $\sigma^{0.95}_{obs}$ are also reported. The last row reports the $p_0$-value of the SM-only hypothesis. For SRs where the data yield is smaller than expected, the $p$-value is truncated at 0.50. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$+V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted post-fit background yields for the SF inclusive SRs. The model independent upper limits at 95% confidence level (CL) on the observed and expected number of beyond the SM events $S^{0.95}_{obs/exp}$ and the effective beyond the SM cross-section $\sigma^{0.95}_{obs}$ are also reported. The last row reports the $p_0$-value of the SM-only hypothesis. For SRs where the data yield is smaller than expected, the $p$-value is truncated at 0.50. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$+V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Distributions of $m_{T2}$ in SRSF-0J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow. Distributions for three benchmark signal points are overlaid for comparison.

Distributions of $m_{T2}$ in SRSF-1J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow. Distributions for three benchmark signal points are overlaid for comparison.

Distributions of $m_{T2}$ in SRDF-0J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow. Distributions for three benchmark signal points are overlaid for comparison.

Distributions of $m_{T2}$ in SRDF-1J for data and the estimated SM backgrounds. The normalisation factors extracted from the corresponding CRs are used to rescale the $t\bar t$, single top, WW, WZ and ZZ backgrounds. The fake and non-prompt leptons background (FNP) is calculated using the data-driven matrix method. The uncertainty band includes all sources of systematic and statistical errors and the last bin includes the overflow. Distributions for three benchmark signal points are overlaid for comparison.

Observed exclusion limits on SUSY simplified models for chargino-pair production with $W$ boson mediated decays. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for chargino-pair production with $W$ boson mediated decays. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for chargino-pair production with slepton/sneutrino mediated mediated decays. The mass relation $m(\tilde{l}_L)=\frac{1}{2}[m(\tilde{\chi}^{\pm}_1 + m(\tilde{\chi}^{0}_1)]$ is assumed. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for chargino-pair production with slepton/sneutrino mediated mediated decays. The mass relation $m(\tilde{l}_L)=\frac{1}{2}[m(\tilde{\chi}^{\pm}_1 + m(\tilde{\chi}^{0}_1)]$ is assumed. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for slepton-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for slepton-pair production. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for selectron-pair production, with left and right handed selectron production combined. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for selectron-pair production, with left and right handed selectron production combined. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for left-handed selectron-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for left-handed selectron-pair production. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for right-handed selectron-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for right-handed selectron-pair production. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for smuon-pair production, with left and right handed smuon production combined. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for smuon-pair production, with left and right handed smuon production combined. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for left-handed smuon-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for left-handed smuon-pair production. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for right-handed smuon-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for right-handed smuon-pair production. All limits are computed at 95% CL.

Observed events and predicted background yields from the fit in the binned DF SRs with $n_{non-b-tagged jets}=0$. For backgrounds whose normalisation is extracted from the fit in the CRs, the yield expected from the simulation before the fit is also reported. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$ +V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted background yields from the fit in the binned DF SRs with $n_{non-b-tagged jets}=1$. For backgrounds whose normalisation is extracted from the fit in the CRs, the yield expected from the simulation before the fit is also reported. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$ +V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted background yields from the fit in the binned SF SRs with $n_{non-b-tagged jets}=0$. For backgrounds whose normalisation is extracted from the fit in the CRs, the yield expected from the simulation before the fit is also reported. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$ +V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed events and predicted background yields from the fit in the binned SF SRs with $n_{non-b-tagged jets}=1$. For backgrounds whose normalisation is extracted from the fit in the CRs, the yield expected from the simulation before the fit is also reported. The background denoted as "Other" in the Table includes the non-dominant background sources for this analysis, i.e. Z+jets, $t\bar t$ +V, Higgs and Drell-Yan events. A "–" symbol indicates that the background contribution is negligible.

Observed exclusion limits on SUSY simplified models for left-handed slepton-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for left-handed slepton-pair production. All limits are computed at 95% CL.

Observed exclusion limits on SUSY simplified models for right-handed slepton-pair production. All limits are computed at 95% CL.

Expected exclusion limits on SUSY simplified models for right-handed slepton-pair production. All limits are computed at 95% CL.

Upper limits on signal cross-section (fb) for chargino-pair production with W -boson-mediated decays.

Upper limits on signal cross-section (fb) for chargino-pair production with slepton/sneutrino-mediated decays. The mass relation $m(\tilde{l}_L)=\frac{1}{2}[m(\tilde{\chi}^{\pm}_1 + m(\tilde{\chi}^{0}_1)]$ is assumed.

Upper limits on signal cross-section (fb) for slepton-pair production.

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[160,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[160,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[160,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[160,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[160,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[160,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[160,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[160,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[120,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[120,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[120,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[120,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[120,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[120,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[120,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[120,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,105).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[100,105).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,105).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[100,105).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,105).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[100,105).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,105).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[100,105).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[105,110).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[105,110).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[105,110).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[105,110).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[105,110).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[105,110).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[105,110).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[105,110).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[110,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[110,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[110,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[110,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[110,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[110,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[110,120).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[110,120).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[120,140).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[120,140).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[120,140).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[120,140).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[120,140).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[120,140).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[120,140).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[120,140).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[140,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[140,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[140,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[140,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[140,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[140,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[140,160).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[140,160).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[160,180).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[160,180).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[160,180).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[160,180).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[160,180).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[160,180).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[160,180).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[160,180).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[180,220).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[180,220).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[180,220).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[180,220).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[180,220).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[180,220).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[180,220).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[180,220).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[220,260).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[220,260).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[220,260).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[220,260).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[220,260).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[220,260).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[220,260).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[220,260).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[260,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-0J-[260,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[260,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-SF-1J-[260,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[260,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-0J-[260,inf).

Signal Acceptance for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[260,inf).

Signal Efficiency for direct chargino-pair production with W-boson mediated decays in SR-DF-1J-[260,inf).

Cutflow for supersymmetric model where $\tilde{\chi}_1^{\pm}\tilde{\chi}_1^{\mp}$ decay via $W^{\pm}W^{\mp}$. The masses of the two charginos are 300 GeV, while the mass of $\tilde{\chi}_1^{0}$ is 50 GeV. The numbers are normalised to the luminosity of 139~fb$^{-1}$.

Cutflow for supersymmetric model where $\tilde{\chi}_1^{\pm}\tilde{\chi}_1^{\mp}$ decay via slepton-neutrino/sneutrino-lepton pair. The masses of the two charginos are 600 GeV, while the mass of $\tilde{\chi}_1^{0}$ is 1 GeV. The slepton/sneutrino masses are 300 GeV. The numbers are normalised to the luminosity of 139~fb$^{-1}$.

Cutflow for supersymmetric model where $\tilde\ell\tilde\ell$ are produced. Only $\tilde{e}$ and $\tilde{\mu}$ are considered in this model. The masses of the two sleptons are 400 GeV, while the mass of $\tilde{\chi}_1^{0}$ is 200 GeV. The numbers are normalised to the luminosity of 139~fb$^{-1}$.

The observed upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed.

The observed upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed.

The observed and expected upper limits at $95\%$ CL on the gluino pair production cross section for a gluino GMSB model with $m_{\tilde{g}}=2400$ GeV. The one (two) standard deviation variation in the expected limit is shown in the inner green (outer yellow) band. The blue solid line shows the observed limit obtained by the CMS displaced jet search

The distribution (normalized to unity) of number of ECAL cells hit in the jet for jets in a background enriched data sample (satisfying $|\eta| < 1.48$, $PV_{\rm track}^{\rm fraction} > 1/12$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets satisfying signal region requirements (except those on $E_{\mathrm{ECAL}}$ and $N^{\mathrm{cell}}_{\mathrm{ECAL}}$).

The distribution (normalized to unity) of $\mathrm{HEF}$ for a data sample enriched in beam halo and noise jets (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30 \,\mathrm{GeV}$, $PV_{\rm track}^{\rm fraction} < 1/12$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$) and for signal jets passing signal region selections (except on $\mathrm{HEF}$).

The distribution (normalized to unity) of $E_{\mathrm{HCAL}}$ for a data sample enriched in beam halo and noise jets (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30 \,\mathrm{GeV}$, $PV_{\rm track}^{\rm fraction} < 1/12$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$) and for signal jets passing signal region selections (except on $E_{\mathrm{HCAL}}$).

The distribution (normalized to unity) of $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ for data sample enriched in jets from noise (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $PV_{\rm track}^{\rm fraction} > 1/12$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20 \,\mathrm{GeV}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets passing signal region selections (except on $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ and $t^{\mathrm{RMS}}_\mathrm{jet}$)

The distribution (normalized to unity) of $t^{\mathrm{RMS}}_\mathrm{jet}$ for data sample enriched in jets from noise (satisfying $|\eta| < 1.48$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $PV_{\rm track}^{\rm fraction} > 1/12$, $\mathrm{HEF} > 0.2$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20 \,\mathrm{GeV}$ and $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$) and for signal jets passing signal region selections (except on $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}}$ and $t^{\mathrm{RMS}}_\mathrm{jet}$)

The distribution (normalized to unity) of $PV_{\rm track}^{\rm fraction}$ for a data sample enriched in main bunch backgrounds (satisfying $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $|t_{\mathrm{jet}}| < 3\,\mathrm{ns}$ and $E_{\mathrm{ECAL}} >20\,\mathrm{GeV}$) and for signal jets passing signal selections (except on $PV_{\rm track}^{\rm fraction}$).

The distribution (normalized to unity) of $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}}$ for a data sample enriched in beam halo (satisfying $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $PV_{\rm track}^{\rm fraction} < 1/12$, $\mathrm{HEF} < 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} < -3\,\mathrm{ns}$ and $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$) and for signal jets passing signal selections (except on $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}}$)

The distribution (normalized to unity) of $\mathrm{max}(\Delta \phi_{\mathrm{DT}})$ for a data sample enriched in cosmic muons (satisfying $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $PV_{\rm track}^{\rm fraction} < 1/12$, $\mathrm{HEF} > 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} > 3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and failing the HCAL noise rejection quality filters) and for signal jets passing signal selections (except on $\mathrm{max}(\Delta \phi_{\mathrm{DT}})$).

The distribution (normalized to unity) of $\mathrm{max}(\Delta \phi_{\mathrm{RPC}})$ for a data sample enriched in cosmic muons (satisfying $E^{\mathrm{CSC}}_\mathrm{ECAL}/E_{\mathrm{ECAL}} < 0.8$, $p_\mathrm{T} > 30\,\mathrm{GeV}$, $|\eta| < 1.48$, $PV_{\rm track}^{\rm fraction} < 1/12$, $\mathrm{HEF} > 0.2$, $t^{\mathrm{RMS}}_\mathrm{jet}/t_{\mathrm{jet}} < 0.4$, $t_{\mathrm{jet}} > 3\,\mathrm{ns}$, $E_{\mathrm{ECAL}} > 20\,\mathrm{GeV}$ and failing the HCAL noise rejection quality filters) and for signal jets passing signal selections (except on $\mathrm{max}(\Delta \phi_{\mathrm{RPC}})$).

Distribution of $t_{\mathrm{jet}}$ for jets with the full Run 2 dataset with no cleaning selection applied (a) and after all jet cleaning selections are applied (b) in events satisfying the trigger requirements and satisfying $p_T^\mathrm{miss} > 300$. The jets are required to pass an inverted selection of $PV_{\rm track}^{\rm fraction} > 1/12$ to enrich the sample in those originating from main and satellite bunch backgrounds. The cleaning selections are shown to reduce the backgrounds by many orders of magnitude.

Distribution of $t_{\mathrm{jet}}$ for jets with the full Run 2 dataset in events satisfying the trigger requirements and satisfying $p_T^\mathrm{miss} < 300$. The jets are required to pass an inverted selection of $PV_{\rm track}^{\rm fraction} > 1/12$ to enrich the sample in those originating from main and satellite bunch backgrounds (all other jet cleaning selections are applied). Clear contributions from jets from satellite bunch collisions can be seen peaked around -5, 5 and 10 ns.

Contribution to the delayed time of jets from the $\beta$ of the gluino is plotted against the delay contribution from the difference (assuming straight line paths) between the path taken by the gluino and the particle forming the jet from the path length for a particle travelling directly to the same position on the ECAL barrel for gluino $c\tau_{0} = 10$ m and mass of 1000 Gev (top) and 3000 GeV (bottom). The dominant contribution is shown to be the gluino $\beta$.

Selection efficiencies for the GMSB model with $m_{\tilde{g}}=1000$ and various proper decay lengths

Selection efficiencies for the GMSB model with $m_{\tilde{g}}=2400$ and various proper decay lengths

Selection efficiencies for the GMSB model with $m_{\tilde{g}}=3000$ and various proper decay lengths

The observed upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected plus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected minus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section relative to the theoretical cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected plus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

The expected minus 1 $\sigma_{\mathrm{experiment}}$ upper limits at $95\%$ CL for the gluino pair production cross section in the GMSB model, shown in the plane of $m_{\tilde{g}}$ and $c\tau_{0}$. A branching fraction of $100\%$ for the gluino decay to a gluon and a gravitino is assumed. The area below the thick black curve represents the observed exclusion region, while the dashed red lines indicate the expected limits and their $\pm 1$ standard deviation ranges from experimental uncertainties. The thin black lines show the effect of the theoretical uncertainties on the signal cross section.

Expected and observed events in the 16 discovery regions along with the according control regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic R-hadron search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector R-hadron search.

p0-values and model-independent upper limits on cross-section x acceptance x efficiency for the 16 discovery regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic search for metastable gluino R-hadrons.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector direct-stau search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector chargino search.

Upper cross-section limit in gluino R-hadron search.

Upper cross-section limit in sbottom R-hadron search.

Upper cross-section limit in stop R-hadron search.

Upper cross-section limit in stau search.

Upper cross-section limit in chargino search.

Lower mass limit as function of gluino lifetime.

Acceptance x efficiency, acceptance and efficiency for the full range of simulated masses in the MS-agnostic R-hadron search.

Upper cross-section limit in meta-stable gluino R-hadron search.

Flavor composition of 800 GeV stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

ETmiss trigger efficiency as function of true ETmiss (EtmissTurnOn).

Single-muon trigger efficiency as function of $|\eta|$ and $\beta$ (SingleMuTurnOn).

Candidate reconstruction efficiency for ID+Calo selection (IDCaloEff).

Candidate reconstruction efficiency for loose selection (LooseEff).

Efficiency for a loose candidate to be promoted to a tight candidate (TightPromotionEff).

Resolution and average of reconstructed dE/dx mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for FullDet candidates.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels in the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for bbA prodcution and the branching fraction $\mathrm{A} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1 b tag category, as a function of the A boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

The product of acceptance, efficiency, and branching fraction of the qbX signal with $\mathrm{X} \rightarrow \tau\tau$ and $m_{\mathrm{B}}=170~\mathrm{GeV}$ in $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels of the 1b1f and 1b1c event categories, as a function of the X boson mass. The selections are as described in the paper. The uncertainty refers to the statistical uncertainty only.

Observed $m_{\tau\tau}$ distribution in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1f category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1c category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mu\tau_\mathrm{h}$ channel of the 1b1f category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed $m_{\tau\tau}$ distribution in the $\mu\tau_\mathrm{h}$ channel of the 1b1c category, compared to the expected SM background contributions. The distributions of the qbX signal with $m_{\mathrm{B}}=170~\mathrm{GeV}$, and X boson mass 40 and 60 GeV are overlaid to illustrate the sensitivity. They are normalized to the cross section times branching fraction of 20 pb. The uncertainty band represents the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panel shows the ratio between the observed and expected events in each bin.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 300 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 450 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained from the combination of the $\mathrm{e}\tau_\mathrm{h}$ and $\mu\tau_\mathrm{h}$ channels, and the 1b1f and 1b1c categories, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1f category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mathrm{e}\tau_\mathrm{h}$ channel of the 1b1c category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1b1f category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of the the cross section for the prodcution of the qbX signal with mB = 170 GeV and the branching fraction $\mathrm{X} \rightarrow \tau\tau$, obtained in the $\mu\tau_\mathrm{h}$ channel of the 1b1c category, as a function of the X boson mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits.

Measured misidentification probability distribution as a function of track multiplicity for the EMJ-1 criteria group defined in Table 2. The red up-pointing triangles are for b jets while the blue down-pointing triangles are for light-flavor jets. The horizontal lines on the data points indicate the variable bin width.

The $H_\mathrm{T}$ distribution for the observed data events (black points) and the predicted background estimation (blue) for selection set 8 (SM QCD-enhanced), requiring at least two jets tagged by loose emerging jet criteria.

The number of associated tracks distribution for the observed data events (black points) and the predicted background estimation (blue) for selection set 8 (SM QCD-enhanced), requiring at least two jets tagged by loose emerging jet criteria.

The $H_\mathrm{T}$ distribution for the observed data events (black points) and the predicted background estimation (blue) for selection set 9 (SM QCD-enhanced), requiring at least one jet tagged by loose emerging jet criteria and large $p_\mathrm{T}^\mathrm{miss}$.

The number of associated tracks distribution for the observed data events (black points) and the predicted background estimation (blue) for selection set 9 (SM QCD-enhanced), requiring at least one jet tagged by loose emerging jet criteria and large $p_\mathrm{T}^\mathrm{miss}$.

Upper limits at 95% CL on the signal cross section for models with dark pion mass $(m_{\pi_\mathrm{DK}})$ of 5 GeV in the proper decay length $(c\tau_{\pi_\mathrm{DK}})$ versus dark mediator mass $(m_{\mathrm{X_{DK}}})$ plane.

Dark mediator mass exclusion limits at 95% CL derived from theoretical cross sections for models with dark pion mass $(m_{\pi_\mathrm{DK}})$ of 5 GeV in the proper decay length $(c\tau_{\pi_\mathrm{DK}})$ versus dark mediator mass $(m_{\mathrm{X_{DK}}})$ plane.

2$\tau$ HighMass SR eff.

2$\tau$ multibin SR eff.

2$\tau$ GMSB SR eff.

1$\tau$ Compressed SR eff.

1$\tau$ MediumMass SR eff.

2$\tau$ Compressed SR eff.

2$\tau$ HighMass SR eff.

2$\tau$ multibin SR eff.

2$\tau$ GMSB SR eff.

1$\tau$ Compressed SR acceptance.

1$\tau$ MediumMass SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ GMSB SR acceptance.

1$\tau$ Compressed SR acceptance.

1$\tau$ MediumMass SR acceptance.

2$\tau$ Compressed SR acceptance.

2$\tau$ HighMass SR acceptance.

2$\tau$ multibin SR acceptance.

2$\tau$ GMSB SR acceptance.

Cutflow table of the $1\tau$ compressed SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $1\tau$ medium-mass SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ compressed SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ high-mass SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ multibin SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Cutflow table of the $2\tau$ GMSB SR for the four signal benchmark scenarios of low, medium, and high mass-splitting in the simplified model as well as the GMSB model.

Best performing fit setups entering the final combination as a function of the LSP mass and the gluino mass. 'S' marks the simultaneous fit of the four simplified model single-bin SRs, 'M' denotes the simultaneous fit of the two $1\tau$ SRs and the $2\tau$ multibin SR.

Observed exclusion contour at 95% CL as a function of tanBeta and the SUSY-breaking mass scale Lambda.

Expected exclusion contour at 95% CL as a function of tanBeta and the SUSY-breaking mass scale Lambda.

Observed exclusion contour at 95% CL as a function of the LSP mass and the gluino mass.

Expected exclusion contour at 95% CL as a function of the LSP mass and the gluino mass.

Observed upper limits on the production cross section at 95% CL in pb as a function of tanBeta and SUSY breaking mass scale Lambda.

Observed upper limits on the production cross section at 95% CL in pb as a function of the LSP mass and the gluino mass.

Yields of the expected background from the SM in the bins of the multibin SR of the $2\tau$ channel with all bins being simultaneously used to constrain the background prediction. Expectation is given with the scalings computed in the combined fit applied. Uncertainties are statistial plus systematrics. Only the subsamples contributing the respective region are considered.

$m_{\mathrm{T}}^{\tau}$ in the compressed $m_{\mathrm{T}}^{\tau}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$E_{\mathrm{T}}^{\mathrm{miss}}$ in the compressed $E_{\mathrm{T}}^{\mathrm{miss}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau}$ in the medium-mass $m_{\mathrm{T}}^{\tau}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$E_{\mathrm{T}}^{\mathrm{miss}}$ in the medium-mass $E_{\mathrm{T}}^{\mathrm{miss}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$H_{\mathrm{T}}$ in the medium-mass $H_{\mathrm{T}}$ VR of the $1\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau_1}$ + $m_{\mathrm{T}}^{\tau_2}$ in the top VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$H_{\mathrm{T}}$ in the $W$ VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau_1}$ + $m_{\mathrm{T}}^{\tau_2}$ in the $Z$ VR of the $2\tau$ channel, illustrating the background modeling after the fit. The last bin includes overflow events.

$m_{\mathrm{T}}^{\tau}$ in the compressed SR of the $1\tau$ channel before application of the $m_{\mathrm{T}}^{\tau}$ > 80 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the medium-mass SR of the $1\tau$ channel before application of the $H_{\mathrm{T}}$ > 1000 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$m_{\mathrm{T}}^{\mathrm{sum}}$ in the compressed SR of the $2\tau$ channel before application of the $m_{\mathrm{T}}^{\mathrm{sum}}$ > 1600 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the high-mass SR of the $2\tau$ channel before application of the $H_{\mathrm{T}}$ > 1100 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

mT(tau_1) + mT(tau_2) in the multibin SR of the 2T channel. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

$H_{\mathrm{T}}$ in the GMSB SR of the $2\tau$ channel before application of the $H_{\mathrm{T}}$ > 1900 GeV requirement. The last bin includes overflow events. Signal predictions corresponding to the simplified model scenarios of low (LM), medium (MM), and high mass-splitting (HM) as well as for the GMSB benchmark are given.

Distribution of $R_\text{ISR}$ for SRC1-5 after the likelihood fit. The stacked histograms show the SM expectation and the hatched uncertainty band around the SM expectation shows the MC statistical and detector-related systematic uncertainties. A representative signal point is shown for each distribution.

Distribution of $m_\text{T}^{b,\text{max}}$ for SRD-high after the likelihood fit. The stacked histograms show the SM expectation and the hatched uncertainty band around the SM expectation shows the MC statistical and detector-related systematic uncertainties. A representative signal point is shown for each distribution.

Distribution of $H_\text{T}$ for SRE after the likelihood fit. The stacked histograms show the SM expectation and the hatched uncertainty band around the SM expectation shows the MC statistical and detector-related systematic uncertainties. A representative signal point is shown for each distribution.

Expected (blue solid line) exclusion limits at 95% CL as a function of stop and LSP masses in the scenario where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Observed (red solid line) exclusion limits at 95% CL as a function of stop and LSP masses in the scenario where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses and branching fraction to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ in the Natural SUSY-inspired mixed grid scenario where $m_{\tilde{\chi^{\pm}_{1}}}=m_{\tilde{\chi^{0}_{1}}}$+1 GeV.

Expected exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a small tan$\beta$ assumption.

Observed exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a small tan$\beta$ assumption.

Expected exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a large tan$\beta$ assumption.

Observed exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a large tan$\beta$ assumption.

Expected exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a small right-handed top-squark mass parameter assumption.

Observed exclusion limits at 95% CL as a function of $m_{\tilde{t}}$ and $m_{\tilde{\chi^{0}_{1}}}$ for the pMSSM-inspired non-asymptotic Higgsino simplified model for a small right-handed top-squark mass parameter assumption.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for the Wino NLSP pMSSM model for a negative value of $\mu$.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for the Wino NLSP pMSSM model for a negative value of $\mu$.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for the Wino NLSP pMSSM model for a positive value of $\mu$.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for the Wino NLSP pMSSM model for a positive value of $\mu$.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for for the left-handed top-squark mass parameter scan in the well-tempered pMSSM model.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for for the left-handed top-squark mass parameter scan in the well-tempered pMSSM model.

Expected exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for for the right-handed top-squark mass parameter scan in the well-tempered pMSSM model.

Observed exclusion limits at 95% CL as a function of $\tilde{t}$ and $\tilde{\chi^{0}_{1}}$ masses for for the right-handed top-squark mass parameter scan in the well-tempered pMSSM model.

Expected exclusion limits at 95% CL exclusion as a function of $\tilde{g}$ and $\tilde{t}$ masses in the scenario where both gluinos decay via $\tilde{g}\to t \tilde{t}\to t\tilde{\chi^{0}_{1}}+$soft and $\Delta m(\tilde{t},\tilde{\chi^{0}_{1}})=5$ GeV.

Observed exclusion limits at 95% CL exclusion as a function of $\tilde{g}$ and $\tilde{t}$ masses in the scenario where both gluinos decay via $\tilde{g}\to t \tilde{t}\to t\tilde{\chi^{0}_{1}}+$soft and $\Delta m(\tilde{t},\tilde{\chi^{0}_{1}})=5$ GeV.

Results of the exclusion fits for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ from the combination of SRA, SRB, SRC, SRD and SRE, based on the best expected $CL_s$. The numbers centered on the grid points indicate which of the signal regions gave the best expected $CL_s$ (with 1, 2, 3, 4, 5, 6 corresponding to SRA, SRB, SRC, SRD-low,SR D-high, SRE respectively).

Results of the exclusion fits in the grid with two stop decay channels: $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV. The results are shown as a function of the branching ratio to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$: 0%. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRC, SRD-low and SRD-high, The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high, 5: SRC).

Results of the exclusion fits in the grid with two stop decay channels: $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV. The results are shown as a function of the branching ratio to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$: 25%. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRC, SRD-low and SRD-high, The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high, 5: SRC).

Results of the exclusion fits in the grid with two stop decay channels: $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV. The results are shown as a function of the branching ratio to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$: 50%. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRC, SRD-low and SRD-high, The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high, 5: SRC).

Results of the exclusion fits in the grid with two stop decay channels: $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV. The results are shown as a function of the branching ratio to $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$: 0 (top left), 25% (top right), 50% (middle left), 75% middle right) and 100% (bottom). The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRC, SRD-low and SRD-high, The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high, 5: SRC).

Results of the exclusion fits in the wino NLSP grid for negative values of $\mu$. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 4: SRD-low, 5: SRD-high).

Results of the exclusion fits in the wino NLSP grid for positive values of $\mu$. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 4: SRD-low, 5: SRD-high).

Results of the exclusion fits in the well-tempered neutralino grid for the $m_{q3L}$ scenario. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 4: SRD-low, 5: SRD-high).

Results of the exclusion fits in the well-tempered neutralino grid for the $m_{tR}$ scenario. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 4: SRD-low, 5: SRD-high).

Results of the exclusion fits in the non-asymptotic higgsino grid with $m(\tilde{\chi^{\pm}_{1}}) - m(\tilde{\chi^{0}_{1}}) = 5$ GeV. A scenarios with large tan$\beta$ (top left) is shown. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high).

Results of the exclusion fits in the non-asymptotic higgsino grid with $m(\tilde{\chi^{\pm}_{1}}) - m(\tilde{\chi^{0}_{1}}) = 5$ GeV. A scenarios with small tan$\beta$ are shown. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high).

Results of the exclusion fits in the non-asymptotic higgsino grid with $m(\tilde{\chi^{\pm}_{1}}) - m(\tilde{\chi^{0}_{1}}) = 5$ GeV. A scenario with a mostly right-handed top squark partner is shown. The results are based on taking the signal region with the best expected $CL_s$, using SRA, SRB, SRD-low and SRD-high, where SRA and SRB are the statistical combinations of their respective regions. The numbers centered on the grid points indicate the signal region used (1: SRA, 2: SRB, 3: SRD-low, 4: SRD-high).

Acceptance for SRA-TT for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRB-TT for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRE for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRA-TW for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRB-TW for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRA-T0 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRB-T0 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRC1 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRC2 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRC3 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRC4 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRC5 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Acceptance for SRD-low for the natural SUSY-inspired mixed grid in which two decay modes are considered, the $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV.

Acceptance for SRD-high for the natural SUSY-inspired mixed grid in which two decay modes are considered, the $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV.

Acceptance for SRE for gluino pair production in the case where both gluinos decay via $\tilde{g}\to t \tilde{t} \to t\tilde{\chi^{0}_{1}}+$soft and $\Delta m(\tilde{t},\tilde{\chi^{0}_{1}})=5$ GeV.

Efficiencies for SRA-TT for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRB-TT for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRE for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRA-TW for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRB-TW for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRA-T0 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRB-T0 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRC1 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRC2 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRC3 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRC4 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRC5 for top squark pair production in the case where both top squarks decay via $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$.

Efficiencies for SRD-low for the natural SUSY-inspired mixed grid in which two decay modes are considered, the $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV.

Efficiencies for SRD-high for the natural SUSY-inspired mixed grid in which two decay modes are considered, the $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ and $\tilde{t}\to b \tilde{\chi^{\pm}_{1}} \to b W^{(*)} \tilde{\chi^{0}_{1}}$, with $m(\tilde{\chi^{\pm}_{1}})-m(\tilde{\chi^{0}_{1}}) = 1$ GeV.

Efficiencies for SRE for gluino pair production in the case where both gluinos decay via $\tilde{g}\to t \tilde{t} \to t\tilde{\chi^{0}_{1}}+$soft and $\Delta m(\tilde{t},\tilde{\chi^{0}_{1}})=5$ GeV.

Upper limit cross-section, in femtobarn, for the $\tilde{t}\to t^{(*)} \tilde{\chi^{0}_{1}}$ grid.

Upper limit cross-section, in femtobarn, for the $\tilde{g}\to t \tilde{t} \to t\tilde{\chi^{0}_{1}}+$soft grid.

Cutflow for SRA for a signal model with top squark pair production in the case where both top squarks decay via $\tilde{t}_1\to t^{(*)} \tilde\chi^0_1$ with $m(\tilde{t}_1,\tilde\chi^0_1)=$ (800,1) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow SRB for a signal model with top squark pair production in the case where both top squarks decay via $\tilde{t}_1\to t^{(*)} \tilde\chi^0_1$ with $m(\tilde{t}_1,\tilde\chi^0_1)=$ (600,300) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow for SRC1 and SRC2 for a signal model with top squark pair production in the case where both top squarks decay via $\tilde{t}_1\to t^{(*)} \tilde\chi^0_1$ with $m(\tilde{t}_1,\tilde\chi^0_1)=$ (250,77) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow for SRC3, SRC4, and SRC5 for a signal model with top squark pair production in the case where both top squarks decay via $\tilde{t}_1\to t^{(*)} \tilde\chi^0_1$ with $m(\tilde{t}_1,\tilde\chi^0_1)=$ (500,327) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow for SRD-high for a signal model with bottom squark pair production in the case where both bottom squarks decay via $b\tilde\chi^{\pm}_1\to bW^{(*)} \tilde\chi^0_1$, with $m(\tilde\chi^{\pm}_1)-m(\tilde\chi^0_1) = 1$ GeV with $m(\tilde{t}_1,\tilde\chi^{\pm}_1)=$ (800,100) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow for SRD-low for a signal model with bottom squark pair production in the case where both bottom squarks decay via $b\tilde\chi^{\pm}_1\to bW^{(*)} \tilde\chi^0_1$, with $m(\tilde\chi^{\pm}_1)-m(\tilde\chi^0_1) = 1$ GeV witht $m(\tilde{t}_1,\tilde\chi^{\pm}_1)=$ (600,200) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Cutflow for SRE for a signal model with gluino pair production in the case where both gluinos decay via $\tilde{g}\to t\tilde{t}_1\to t\tilde\chi^0_1+$soft and $\Delta m(\tilde{t}_1, \tilde\chi^0_1)=5$ GeV with $m(\tilde{g},\tilde{t}_1)=$ (1700,400) GeV. An integrated luminosity of 36.1 $\mathrm{fb}^{-1}$ is assumed when calculating the weighted yields. For the derivation skim at least one of the following four criteria is required: $H_{\mathrm{T}}$ $>$ 150 GeV; at least one loose electron with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two loose electrons with $p_{\mathrm{T}}$ $>$ 20 GeV; at least one muon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two muons with $p_{\mathrm{T}}$ $>$ 20 GeV; or at least one photon with $p_{\mathrm{T}}$ $>$ 100 GeV or at least two photons with $p_{\mathrm{T}}$ $>$ 50 GeV.

Upper limits at the 95% CL on the cross section for SSM W' production and decay to the muon+neutrino channel as a function of the W' pole mass.

Combined (electron and muon channels) upper limits at the 95% CL on the cross section for SSM W' production and decay to a single lepton generation as a function of the W' pole mass.

Product of acceptance and efficiency for the electron and muon selections as a function of the W' pole mass.

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Despite listing this as an exclusive final state (as there must be at least one b-jets), there is no explicit selection on the presence of additional light-flavour jets. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 300, 500 and 800 GeV and $\tan\beta$ = 10 in the hMSSM scenario are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. In the paper, the first bin is cut off at 60 GeV for aesthetics but contains underflows down to 50 GeV as in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and predicted mTtot distribution for the b-inclusive selection in the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The prediction for a SSM Zprime with masses of 1500, 2000 and 2500 GeV are also provided.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed and expected 95% CL upper limits on the Drell Yan production cross section times ditau branching fraction as a function of the Zprime boson mass.

Observed and expected 95% CL upper limits on the Higgs boson production cross section times ditau branching fraction as a function of the boson mass and the relative strength of the b-associated production.

Ratio of the 95% CL upper limits on the production cross section times branching fraction for alternate Zprime models with respect to the SSM, both observed and expected are shown.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance, acceptance times efficiency and b-tag category fraction for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Acceptance and acceptance times efficiency for a heavy gauge boson produced by Drell Yan as a function of the gauge boson mass.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times braching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the Higgs boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. Vaules are provided for the fit to the observed data and to the expected data, which is the sum of Standard Model contributions not including the SM Higgs boson. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed event yields and estimated backgrounds in the low- and high-mass control regions. The uncertainty in the background yield is the sum in quadrature of the statistical and systematic uncertainties.

Summary of the relative systematic uncertainties in heavy Majorana neutrino signal yields and the background from prompt same-sign leptons, both estimated from simulation. The relative systematic uncertainties assigned to the data-driven backgrounds for the misidentified lepton background and mismeasured charge background are also shown. The uncertainties are given for the low-mass selections.

Summary of the relative systematic uncertainties in heavy Majorana neutrino signal yields and the background from prompt same-sign leptons, both estimated from simulation. The relative systematic uncertainties assigned to the data-driven backgrounds for the misidentified lepton background and mismeasured charge background are also shown. The uncertainties are given for the high-mass selections.

Summary of contributions to the systematic uncertainty related to the prompt same-sign leptons background, misidentified lepton background, and mismeasured charge background on the total background uncertainty for the case of $m_{\rm N}$ = 100 and 500 GeV.

Observed event yields and estimated backgrounds after the application of all selection, except for the final optimization. The background predictions from prompt same-sign leptons (Prompt bkgd.), misidentified leptons (Misid. bkgd.), mismeasured charge (Charge mismeas. bkgd.) and the total background (Total bkgd.) are shown together withthe number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Dielectron channel results after final optimization. The background predictions from prompt same-sign leptons, misidentified leptons, and mismeasured charge are shown along with the total background estimate and the number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

Electron-muon channel results after final optimization. The background predictions from prompt same-sign leptons and misidentified leptons are shown along with the total background estimate and the number of events observed in data. The uncertainties shown are the statistical and systematic components, respectively.

HT inclusive distribution in SRZ.

Jet multiplicity distribution in SRZ.

B-tagged jet multiplicity distribution in SRZ.

Dilepton invariant mass distribution in SRMedium.

Dilepton invariant mass distribution in SRMedium.

Dilepton invariant mass in SRHigh.

Expected 95% exclusion contour for the gluino simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Observed 95% exclusion contour for the gluino simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Expected 95% exclusion contour for the squark simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Observed 95% exclusion contour for the squark simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Expected 95% exclusion contour for the gluino simplified model with varying N1 mass and decays to quarks and Z bosons in SRZ.

Observed 95% exclusion contour for the gluino simplified model with varying N1 mass and decays to quarks and Z bosons in SRZ.

Expected 95% exclusion contour for the gluino simplified model with sleptons in the edge SRs.

Observed 95% exclusion contour for the gluino simplified model with sleptons in the edge SRs.

Expected 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in the edge SRs.

Observed 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in the edge SRs.

Expected 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in SRZ.

Observed 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in SRZ.

Cutflow table for three benchmark signal points from the $\tilde{g}$--$\chi_{2}^{0}$ on-shell grid in SRZ for the $ee$ and $\mu\mu$ channels separately. The number of generated events is N RAW, while N WEIGHTED corresponds to the number of events expected in 14.7/fb, after all MC simulation corrections have been applied.

Cutflow table for three benchmark signal points for the slepton model in SR-low, SR-medium and SR-high for the $ee$ and $\mu\mu$ channels separately. The number of generated events is N RAW, while N WEIGHTED corresponds to the number of events expected in 14.7/fb, after all MC simulation corrections have been applied. The event generation was filtered to require at least two leptons (electrons or muons) with $p_{T}>5$ GeV and $|\eta|<2.8$, with 20000 events required to pass the filter. The corresponding number of unfiltered events is provided in the ``Pre-filter events'' row of the table.

Signal acceptance for the gluino simplified model with decays to quarks and Z bosons in SRZ.

Signal efficiency for the gluino simplified model with decays to quarks and Z bosons in SRZ.

Signal acceptance for the gluino simplified model with decays via sleptons in SRLow.

Signal efficiency for the gluino simplified model with decays via sleptons in SRLow.

Signal acceptance for the gluino simplified model with decays via sleptons in SRMedium.

Signal efficiency for the gluino simplified model with decays via sleptons in SRMedium.

Signal acceptance for the gluino simplified model with decays via sleptons in SRHigh.

Signal efficiency for the gluino simplified model with decays via sleptons in SRHigh.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Upper limits on the signal cross-section at 95% CL for the squark simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with varying N1 mass and decays to quarks and Z bosons in SRZ.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with decays via sleptons.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with decays to on-shell or off-shell Z bosons.

Median expected and observed 95% CL upper limits on the cross-section times branching ratio (sigma*B) for W'_SSM production for the exclusive muon and electron channels, and for both channels combined.

Kinematic distributions for "extended SR selections" of the 2-tau channel, for mTsum in the Compressed SR selection without the mTsum>1400 GeV requirement. The last bin includes overflow events. Cited uncertainties are statistical uncertainties only. Signal predictions are overlaid for several benchmark models, normalised to their predicted cross sections. For the simplified model, "MM" refers to a medium mass splitting, with m(gluino)=1145 GeV and m(neutralino)=265 GeV; "HM" denotes a high mass splitting scenario, with m(gluino)=1305 GeV and m(neutralino)=105 GeV. The GMSB benchmark model corresponds to Lambda = 90 TeV and tanbeta = 40.

Kinematic distributions for "extended SR selections" of the 2-tau channel, for mTtau1+mTtau2 in the High-Mass SR selection without the mTtau1+mTtau2>350GeV requirement. The last bin includes overflow events. Cited uncertainties are statistical uncertainties only. Signal predictions are overlaid for several benchmark models, normalised to their predicted cross sections. For the simplified model, "MM" refers to a medium mass splitting, with m(gluino)=1145 GeV and m(neutralino)=265 GeV; "HM" denotes a high mass splitting scenario, with m(gluino)=1305 GeV and m(neutralino)=105 GeV. The GMSB benchmark model corresponds to Lambda = 90 TeV and tanbeta = 40.

Kinematic distributions for "extended SR selections" of the 2-tau channel, for HT in the GMSB SR selection without the HT > 1700 GeV requirement. The last bin includes overflow events. Cited uncertainties are statistical uncertainties only. Signal predictions are overlaid for several benchmark models, normalised to their predicted cross sections. For the simplified model, "MM" refers to a medium mass splitting, with m(gluino)=1145 GeV and m(neutralino)=265 GeV; "HM" denotes a high mass splitting scenario, with m(gluino)=1305 GeV and m(neutralino)=105 GeV. The GMSB benchmark model corresponds to Lambda = 90 TeV and tanbeta = 40.

Expected exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on the combined results from the 1tau and 2tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Observed exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on the combined results from the 1tau and 2tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Expected exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on results from the 2tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Expected exclusion contour at the 95% confidence level for the simplified model of gluino pair production, based on results from the 1tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data.

Observed exclusion contours at the 95% confidence level for the gauge-mediated supersymmetry-breaking model, based on results from the 2 tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data. Additional model parameters are M(mess) = 250 TeV, N5 = 3, mu>0 and Cgrav =1.

Expected exclusion contours at the 95% confidence level for the gauge-mediated supersymmetry-breaking model, based on results from the 2 tau channel. The result is obtained using 3.2 fb-1 of sqrt(s) = 13 TeV ATLAS data. Additional model parameters are M(mess) = 250 TeV, N5 = 3, mu>0 and Cgrav =1.

Observed upper cross section limits in pb for the simplified model of gluino pair production for the combination of all SRs.

Best expected signal region for the simplified model of gluino pair production. The respective SR has been used in the combination of the results.

Acceptance for the gluino production simplified model grid in the Compressed 1tau signal region.

Efficiency for the gluino production simplified model grid in the Compressed 1tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the Compressed 1tau signal region.

Acceptance for the gluino production simplified model grid in the medium mass 1tau signal region.

Efficiency for the gluino production simplified model grid in the medium mass 1tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the medium mass 1tau signal region.

Acceptance for the gluino production simplified model grid in the high mass 1tau signal region.

Efficiency for the gluino production simplified model grid in the high mass 1tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the high mass 1tau signal region.

Acceptance for the gluino production simplified model grid in the compressed 2tau signal region.

Efficiency for the gluino production simplified model grid in the compressed 2tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the compressed 2tau signal region.

Acceptance for the gluino production simplified model grid in the high mass 2tau signal region.

Efficiency for the gluino production simplified model grid in the high mass 2tau signal region.

Acceptance times Efficiency for the gluino production simplified model grid in the high mass 2tau signal region.

Acceptance for the GMSB model grid in the 2tau signal region.

Efficiency for the GMSB model grid in the 2tau signal region.

Acceptance times Efficiency for the GMSB model grid in the 2tau signal region.

Observed and expected background and signal effective mass distributions for SR2jl. For signal, a squark direct decay model with $m(\tilde q)=800$ GeV and $m(\tilde\chi^0_1)=400$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR2jm. For signal, a gluino direct decay model with $m(\tilde g)=750$ GeV and $m(\tilde\chi^0_1)=660$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR2jt. For signal, a squark direct decay model with $m(\tilde q)=1200$ GeV and $m(\tilde\chi^0_1)=0$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR4jt. For signal, a gluino direct decay model with $m(\tilde g)=1400$ GeV and $m(\tilde\chi^0_1)=0$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR5j. For signal, a gluino one-step decay model with $m(\tilde g)=1265$ GeV, $m(\tilde\chi^\pm_1)=945$ GeV and $m(\tilde\chi^0_1)=625$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR6jm. For signal, a gluino one-step decay model with $m(\tilde g)=1265$ GeV, $m(\tilde\chi^\pm_1)=945$ GeV and $m(\tilde\chi^0_1)=625$ GeV is shown.

Observed and expected background and signal effective mass distributions for SR6jt. For signal, a gluino one-step decay model with $m(\tilde g)=1385$ GeV, $m(\tilde\chi^\pm_1)=705$ GeV and $m(\tilde\chi^0_1)=25$ GeV is shown.

Expected limit at 95% CL for squark direct decay model grid.