Distribution of the resonance mass in the c-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the untagged Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the c-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The signal selection efficiency x acceptance summed over all signal regions, for all masses and LQ decay channels considered.

The observed and expected limits for all masses and LQ decay channels considered.

Cutflow Table in the electron channel, considering signal samples with LQ mass of 1 TeV.

Cutflow Table in the muon channel, considering signal samples with LQ mass of 1 TeV.

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

A comparison of the measured cell $E_T^{miss}$ distribution with that predicted by the two-component model for two pile-up scenarios. The magenta points extend the measured distribution using L1 $E_T^{miss} > 30\,$GeV and L1 $E_T^{miss} > 50\,$GeV data. The red curve is the prediction from the calorimeter-resolution part of the model. The green curve is the high $E_T^{miss}$ tail's probability distribution for the mean number of $pp$ interactions $\mu$ in each figure. The blue curve is the full model prediction computed by combining the $E_T^{miss}$ from these two individual sources shown in red and green, each calculated for $\mu=\langle\mu\rangle$. The black points show the unbiased $E_T^{miss}$ distribution measured in data. Corresponds to a prediction for $\langle\mu\rangle = 25$.

A comparison of the measured cell $E_T^{miss}$ distribution with that predicted by the two-component model for two pile-up scenarios. The magenta points extend the measured distribution using L1 $E_T^{miss} > 30\,$GeV and L1 $E_T^{miss} > 50\,$GeV data. The red curve is the prediction from the calorimeter-resolution part of the model. The green curve is the high $E_T^{miss}$ tail's probability distribution for the mean number of $pp$ interactions $\mu$ in each figure. The blue curve is the full model prediction computed by combining the $E_T^{miss}$ from these two individual sources shown in red and green, each calculated for $\mu=\langle\mu\rangle$. The black points show the unbiased $E_T^{miss}$ distribution measured in data. Corresponds to a prediction for $\langle\mu\rangle = 25$.

A comparison of the measured cell $E_T^{miss}$ distribution with that predicted by the two-component model for two pile-up scenarios. The magenta points extend the measured distribution using L1 $E_T^{miss} > 30\,$GeV and L1 $E_T^{miss} > 50\,$GeV data. The red curve is the prediction from the calorimeter-resolution part of the model. The green curve is the high $E_T^{miss}$ tail's probability distribution for the mean number of $pp$ interactions $\mu$ in each figure. The blue curve is the full model prediction computed by combining the $E_T^{miss}$ from these two individual sources shown in red and green, each calculated for $\mu=\langle\mu\rangle$. The black points show the unbiased $E_T^{miss}$ distribution measured in data. Corresponds to a prediction for $\langle\mu\rangle = 55$.

A comparison of the measured cell $E_T^{miss}$ distribution with that predicted by the two-component model for two pile-up scenarios. The magenta points extend the measured distribution using L1 $E_T^{miss} > 30\,$GeV and L1 $E_T^{miss} > 50\,$GeV data. The red curve is the prediction from the calorimeter-resolution part of the model. The green curve is the high $E_T^{miss}$ tail's probability distribution for the mean number of $pp$ interactions $\mu$ in each figure. The blue curve is the full model prediction computed by combining the $E_T^{miss}$ from these two individual sources shown in red and green, each calculated for $\mu=\langle\mu\rangle$. The black points show the unbiased $E_T^{miss}$ distribution measured in data. Corresponds to a prediction for $\langle\mu\rangle = 55$.

A comparison of the measured cell $E_T^{miss}$ distribution with that predicted by the two-component model for two pile-up scenarios. The magenta points extend the measured distribution using L1 $E_T^{miss} > 30\,$GeV and L1 $E_T^{miss} > 50\,$GeV data. The red curve is the prediction from the calorimeter-resolution part of the model. The green curve is the high $E_T^{miss}$ tail's probability distribution for the mean number of $pp$ interactions $\mu$ in each figure. The blue curve is the full model prediction computed by combining the $E_T^{miss}$ from these two individual sources shown in red and green, each calculated for $\mu=\langle\mu\rangle$. The black points show the unbiased $E_T^{miss}$ distribution measured in data. Corresponds to a prediction for $\langle\mu\rangle = 55$.

A comparison of the measured cell $E_T^{miss}$ distribution with that predicted by the two-component model for two pile-up scenarios. The magenta points extend the measured distribution using L1 $E_T^{miss} > 30\,$GeV and L1 $E_T^{miss} > 50\,$GeV data. The red curve is the prediction from the calorimeter-resolution part of the model. The green curve is the high $E_T^{miss}$ tail's probability distribution for the mean number of $pp$ interactions $\mu$ in each figure. The blue curve is the full model prediction computed by combining the $E_T^{miss}$ from these two individual sources shown in red and green, each calculated for $\mu=\langle\mu\rangle$. The black points show the unbiased $E_T^{miss}$ distribution measured in data. Corresponds to a prediction for $\langle\mu\rangle = 55$.

A comparison of the measured cell $E_T^{miss}$ distribution with that predicted by the two-component model for two pile-up scenarios. The magenta points extend the measured distribution using L1 $E_T^{miss} > 30\,$GeV and L1 $E_T^{miss} > 50\,$GeV data. The red curve is the prediction from the calorimeter-resolution part of the model. The green curve is the high $E_T^{miss}$ tail's probability distribution for the mean number of $pp$ interactions $\mu$ in each figure. The blue curve is the full model prediction computed by combining the $E_T^{miss}$ from these two individual sources shown in red and green, each calculated for $\mu=\langle\mu\rangle$. The black points show the unbiased $E_T^{miss}$ distribution measured in data. Corresponds to a prediction for $\langle\mu\rangle = 55$.

The $E_T^{miss}$ model predicted trigger rate as a function of $\mu$ for the cell $E_T^{miss}$ algorithm with a threshold of $80\,$GeV and $120\,$GeV, assuming no additional pile-up mitigation.

The L1 $E_T^{miss}$ trigger efficiency, shown as a function of $p_T(\mu\mu)$ in $Z\rightarrow\mu\mu$ events.

The efficiencies in the plot are shown for events satisfying a $Z\rightarrow\mu\mu$ selection and with $p_T(\mu\mu)$ larger than $150\,$GeV vs pile-up for each of the four years of data taking.

The L1 $E_T^{miss}$ trigger rate as a function of $\langle\mu\rangle$ for runs in three different periods ($\alpha$, $\beta$, $\gamma$) in the year 2017.

The L1 $E_T^{miss}$ trigger rate as a function of $\langle\mu\rangle$ for runs in three different periods ($\alpha$, $\beta$, $\gamma$) in the year 2017.

The L1 $E_T^{miss}$ trigger rate as a function of $\langle\mu\rangle$ for runs in three different periods ($\alpha$, $\beta$, $\gamma$) in the year 2017.

The L1 $E_T^{miss}$ trigger efficiency is shown as a function of mean pile-up for events satisfying a $Z\rightarrow\mu\mu$ selection and with $p_T(\mu\mu)$ larger than $150\,$GeV in three periods during the year 2017.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $0\leq\langle\mu\rangle < 20$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $0\leq\langle\mu\rangle < 20$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $0\leq\langle\mu\rangle < 20$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $0\leq\langle\mu\rangle < 20$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $0\leq\langle\mu\rangle < 20$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $20\leq\langle\mu\rangle < 30$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $20\leq\langle\mu\rangle < 30$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $20\leq\langle\mu\rangle < 30$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $20\leq\langle\mu\rangle < 30$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $20\leq\langle\mu\rangle < 30$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $30\leq\langle\mu\rangle < 40$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $30\leq\langle\mu\rangle < 40$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $30\leq\langle\mu\rangle < 40$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $30\leq\langle\mu\rangle < 40$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $30\leq\langle\mu\rangle < 40$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $40\leq\langle\mu\rangle$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $40\leq\langle\mu\rangle$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $40\leq\langle\mu\rangle$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $40\leq\langle\mu\rangle$.

Background acceptance vs signal efficiency for each of four individual HLT $E_T^{miss}$ algorithms for a $Z\rightarrow\mu\mu$ selection with $p_T(\mu\mu) > 175\,$GeV for data recorded in the year 2017. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger. $40\leq\langle\mu\rangle$.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Relative background acceptance fraction vs. relative efficiency for two different $p_T(\mu\mu)>150\,$GeV threshold for data recorded in the year 2018. Two of the curves show the performance of the stand-alone cell algorithm and the stand-alone pufit algorithm. The other two show combined algorithms each formed by requiring that the event satisfy both a fixed threshold (either 65 GeV or 70 GeV as shown in the legend) for the cell algorithm and a pufit threshold which varies along the curve. In each plot the background acceptance fractions and the efficiencies are relative to those of the pufit $E_T^{miss} > 110\,$GeV trigger and thus can be greater than one. The diamond indicates the performance of the pufit $E_T^{miss}>110\,$GeV trigger while the cross indicates the performance of the combined (pufit $E_T^{miss}>110\,$GeV and cell $E_T^{miss}>50\,$GeV) trigger.

Turn-on efficiency curves are shown for $Z\rightarrow\mu\mu$ events for three algorithms: the cell algorithm alone, the pufit algorithm alone and the combined cell+pufit algorithm. The thresholds are set such that the algorithms have equal rates, and the data were recorded in the year 2018. Here the trigger efficiency is shown with respect to $p_T(\mu\mu)$.

Turn-on efficiency curves are shown for $Z\rightarrow\mu\mu$ events for three algorithms: the cell algorithm alone, the pufit algorithm alone and the combined cell+pufit algorithm. The thresholds are set such that the algorithms have equal rates, and the data were recorded in the year 2018. Here the trigger efficiency is shown with respect to the offline $E_T^{miss}$ calculation with muons treated as being invisible.

Efficiencies for $Z\rightarrow\mu\mu$ events are shown for the L1 $E_T^{miss}>50\,$GeV trigger and for the complete L1+HLT trigger chain that also requires pufit > 110 GeV. Each is shown as a function of $\langle\mu\rangle$ for a $p_T(\mu\mu)$ threshold of 150 GeV.

Efficiencies for $Z\rightarrow\mu\mu$ events are shown for the L1 $E_T^{miss}>50\,$GeV trigger and for the complete L1+HLT trigger chain that also requires pufit > 110 GeV. Each is shown as a function of $\langle\mu\rangle$ for an offline $E_T^{miss}$ threshold of 150 GeV.

Efficiencies for $Z\rightarrow\mu\mu$ events are shown for the L1 $E_T^{miss}>50\,$GeV trigger and for the complete L1+HLT trigger chain that also requires pufit > 110 GeV. Each is shown as a function of $\langle\mu\rangle$ for an $p_T(\mu\mu)$ threshold of 175 GeV.

Efficiencies for $Z\rightarrow\mu\mu$ events are shown for the L1 $E_T^{miss}>50\,$GeV trigger and for the complete L1+HLT trigger chain that also requires pufit > 110 GeV. Each is shown as a function of $\langle\mu\rangle$ for an offline $E_T^{miss}$ threshold of 175 GeV.

HLT_xe70_mht trigger output rate as a function of $\langle\mu\rangle$ shown for an example run in year 2015.

HLT_xe90_mht trigger output rate as a function of $\langle\mu\rangle$ shown for an example run in year 2016.

HLT_xe110_mht trigger output rate as a function of $\langle\mu\rangle$ shown for an example run in year 2016.

HLT_xe110_pufit trigger output rate as a function of $\langle\mu\rangle$ shown for an example run in year 2017. This trigger included an implicit requirement of cell $E_T^{miss}>50\,$GeV.

HLT_xe110_pufit_xe65 trigger output rate as a function of $\langle\mu\rangle$ shown for an example run in year 2018.

HLT_xe110_pufit_xe70 trigger output rate as a function of $\langle\mu\rangle$ shown for an example run in year 2018.

Full-chain trigger efficiencies for each year as a function of $p_T(\mu\mu)$. The efficiency corresponds to that of the lowest unprescaled trigger that is adjusted throughout each year.

Full-chain trigger efficiencies for 2015, 2016, 2017 and 2018 as a function of $\langle\mu\rangle$ for $p_T(\mu\mu)>150$ GeV. The efficiency corresponds to that of the lowest unprescaled trigger that is adjusted throughout the year.

Efficiencies for the first-level trigger L1XE50 and the combined L1+HLT trigger chain HLT_xe110_pufit_xe65_L1XE50 in data recorded in the year 2018 are shown as a function of $\langle\mu\rangle$ for $W\rightarrow e\nu$ and $Z\rightarrow\mu\mu$ selections with offline $E_T^{miss}>150\,$GeV.

Efficiencies for the first-level trigger L1XE50 and the combined L1+HLT trigger chain HLT_xe110_pufit_xe65_L1XE50 in data recorded in the year 2018 are shown as a function of $\langle\mu\rangle$ for $t\bar t$ with offline $E_T^{miss}>150\,$GeV.

Efficiencies for the first-level trigger L1XE50 and the combined L1+HLT trigger chain HLT_xe110_pufit_xe65_L1XE50 in data recorded in the year 2018 are shown as a function of $\langle\mu\rangle$ for $W\rightarrow e\nu$ and $Z\rightarrow\mu\mu$ selections with offline $E_T^{miss}>175\,$GeV.

Efficiencies for the first-level trigger L1XE50 and the combined L1+HLT trigger chain HLT_xe110_pufit_xe65_L1XE50 in data recorded in the year 2018 are shown as a function of $\langle\mu\rangle$ for $t\bar t$ with offline $E_T^{miss}>175\,$GeV.

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. No theoretical uncertainty is considered when computing these limits.

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. No theoretical uncertainty is considered when computing these limits.

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. No theoretical uncertainty is considered when computing these limits.

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. No theoretical uncertainty is considered when computing these limits.

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. No theoretical uncertainty is considered when computing these limits.

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. No theoretical uncertainty is considered when computing these limits.

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. No theoretical uncertainty is considered when computing these limits.

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. No theoretical uncertainty is considered when computing these limits.

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. No theoretical uncertainty is considered when computing these limits.

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. No theoretical uncertainty is considered when computing these limits.

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. No theoretical uncertainty is considered when computing these limits.

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. No theoretical uncertainty is considered when computing these limits.

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.

Observed 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.

Observed 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.

Observed 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.

Observed 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.

Observed 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.

Observed 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.

Observed 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.

Observed 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.

Observed 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.

Observed 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.

Observed 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.

Observed 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.

Observed 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.

Observed 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.

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 theoretical uncertainty of signal cross section is considered.

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 theoretical uncertainty of signal cross section is considered.

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 theoretical uncertainty of signal cross section is considered.

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 theoretical uncertainty of signal cross section is considered.

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 theoretical uncertainty of signal cross section is considered.

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. The theoretical uncertainty of signal cross section is considered.

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 by 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 theoretical uncertainty of signal cross section is considered.

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 by 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 theoretical uncertainty of signal cross section is considered.

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 by 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 theoretical uncertainty of signal cross section is considered.

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 by 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 theoretical uncertainty of signal cross section is considered.

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 by 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 theoretical uncertainty of signal cross section is considered.

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 by 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 theoretical uncertainty of signal cross section is considered.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ 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 theoretical uncertainty of signal cross section is considered.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ 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 theoretical uncertainty of signal cross section is considered.

The observed 95% CL upper limits with one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. 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 theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. 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 theoretical uncertainty of signal cross section is considered.

Observed $CL_S$ limits on $\sigma \times BR$ for $m_{\Phi} = 200-400$ GeV.

Observed $CL_S$ limits on $\sigma \times BR$ for $m_{\Phi} = 600-1000$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 15$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 25$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 40$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 55$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{\Phi} = 200$ GeV, $m_s = 25$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{\Phi} = 200$ GeV, $m_s = 50$ GeV.

Comparison of the IDVx reconstruction and selection efficiency for a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV, using only standard tracking versus using standard and large radius tracking for all ID vertices in the signal MC sample.

Comparison of the IDVx reconstruction and selection efficiency for a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV, using only standard tracking versus using standard and large radius tracking for ID vertices passing the full IDVx selection criteria.

Comparison of the IDVx reconstruction and selection efficiency for a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV, using only standard tracking versus using standard and large radius tracking for all ID vertices in the signal MC sample.

Comparison of the IDVx reconstruction and selection efficiency for a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV, using only standard tracking versus using standard and large radius tracking for ID vertices passing the full IDVx selection criteria.

The IDVx selection efficiency as a function of long-lived particle decay $z$ position for MC signal samples with a 125 GeV Higgs boson decaying to LLPs with masses of 8, 25, and 55 GeV. The efficiency for the 55 GeV LLP sample is shown both with and without the material veto applied.

The IDVx selection efficiency as a function of long-lived particle decay $z$ position for MC signal samples with mediators of masses 200, 400, and 600 GeV decaying to LLPs with a mass of 50 GeV. The efficiency for the 200 GeV mediator sample is shown both with and without the material veto applied.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 400 GeV decaying to an LLP with a mass of 50 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 400 GeV decaying to an LLP with a mass of 50 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 400 GeV decaying to an LLP with a mass of 50 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 8$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 15$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 25$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 40$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 55$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 200$ GeV, $m_s = 8$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 200$ GeV, $m_s = 25$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 200$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 400$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 400$ GeV, $m_s = 100$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 600$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 600$ GeV, $m_s = 150$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 1000$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 1000$ GeV, $m_s = 150$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 1000$ GeV, $m_s = 400$ GeV.

The $L_{T}$ distribution in the MisID-enriched region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 3L below-Z signal region. The last bin contains the overflow events.

The $M_{T}$ distribution in the 3L on-Z signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 3L above-Z signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 3L OSSF0 signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 4L OSSF0 signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 4L OSSF1 signal region. The last bin contains the overflow events.

The $L_{T}+p^{miss}_{T}$ distribution in the 4L OSSF2 signal region. The last bin contains the overflow events.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 3L(ee) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 3L(ee) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 4L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 4L(ee) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 4L(ee) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 4L(ee) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{20}$ distribution in the 4L(ee) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dielectron $M_{OSSF}^{300}$ distribution in the 4L(ee) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {ee})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 0B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 0B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 1B, 400<$S_{T}$<800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 3L($\mu\mu$) 1B, $S_{T}$>800 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$<400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 4L($\mu\mu$) 0B, $S_{T}$>400 GeV signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{20}$ distribution in the 4L($\mu\mu$) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The dimuon $M_{OSSF}^{300}$ distribution in the 4L($\mu\mu$) 1B signal region. The last bin does not contain the overflow events. The signal is shown with $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$=0.05.

The 95% confidence level exclusion limits for the flavor-democratic scenario on the total production cross section of heavy fermion pairs.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(S) and $\mathcal{B}(\phi - {ee})$.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(PS) and $\mathcal{B}(\phi - {ee})$.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(S) and $\mathcal{B}(\phi - {\mu\mu})$.

The 95% confidence level exclusion limits on the product of the total production cross section of tt$\phi$(PS) and $\mathcal{B}(\phi - {\mu\mu})$.

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {ee})$ for tt$\phi$(S).

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {ee})$ for tt$\phi$(PS).

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$ for tt$\phi$(S).

The 95% confidence level exclusion limits on $g_{t}^2\mathcal{B}(\phi - {\mu\mu})$ for tt$\phi$(PS).

Product of the fiducial acceptance and the event selection efficiency for the type-III Seesaw signal model at various signal mass hypotheses calculated after all analysis selection requirements.

Product of the fiducial acceptance and the event selection efficiency for the tt$\phi$ signal models at various signal mass hypotheses calculated after all analysis selection requirements.

Dependence of the overall signal efficiency on the pT of the long-lived Z' for Z' -> ee. The error bars indicate the total uncertainties.

Dependence of the overall signal efficiency on the transverse decay radius Rxy of the long-lived Z' for Z' -> emu. The error bars indicate the total uncertainties.

Dependence of the overall signal efficiency on the pT of the long-lived Z' for Z' -> emu. The error bars indicate the total uncertainties.

Dependence of the overall signal efficiency on the transverse decay radius Rxy of the long-lived Z' for Z' -> mumu. The error bars indicate the total uncertainties.

Dependence of the overall signal efficiency on the pT of the long-lived Z' for Z' -> mumu. The error bars indicate the total uncertainties.

Overall signal efficiency as a function of the mean proper lifetime (ctau) of the neutralino for the lambda121 scenario of the RPV SUSY model. The error bars indicate the total uncertainties.

Overall signal efficiency as a function of the mean proper lifetime (ctau) of the neutralino for the lambda122 scenario of the RPV SUSY model. The error bars indicate the total uncertainties.

95% CL upper limits on the squark-antisquark production cross-section as a function of the mean proper lifetime (ctau) of the neutralino for the lambda121 scenario of the RPV SUSY model and a 700 GeV squark. The uncertainties of the expected limit indicate the +-1sigma variations.

95% CL upper limits on the squark-antisquark production cross-section as a function of the mean proper lifetime (ctau) of the neutralino for the lambda122 scenario of the RPV SUSY model and a 700 GeV squark. The uncertainties of the expected limit indicate the +-1sigma variations.

95% CL upper limits on the squark-antisquark production cross-section as a function of the mean proper lifetime (ctau) of the neutralino for the lambda121 scenario of the RPV SUSY model and a 1600 GeV squark. The uncertainties of the expected limit indicate the +-1sigma variations.

95% CL upper limits on the squark-antisquark production cross-section as a function of the mean proper lifetime (ctau) of the neutralino for the lambda122 scenario of the RPV SUSY model and a 1600 GeV squark. The uncertainties of the expected limit indicate the +-1sigma variations.

Fraction of detector volume covered by the material veto as a function of z and Rxy of the displaced dilepton vertex.

Fraction of detector volume covered by the disabled pixel modules veto veto as a function of z and Rxy of the displaced dilepton vertex.

Observed Rxy distribution of vertices composed of two non-leptonic tracks in a control sample in the data and the predicted distribution obtained from the event mixing. The error bars indicate the statistical uncertainties.

Rxy distributions of Kshort vertices from the large radius tracking in the data and background MC samples. Data has been normalised such that the total number of Kshort from the standard tracking in the data agrees with the total number of Kshort from the standard tracking in the MC. The error bars indicate the statistical uncertainties.

Acceptance per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> eenu. The error bars indicate the statistical uncertainties.

Acceptance per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> emunu. The error bars indicate the statistical uncertainties.

Acceptance per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> mumunu. The error bars indicate the statistical uncertainties.

Detection efficiency per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> eenu. The error bars indicate the total uncertainties.

Detection efficiency per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> emunu. The error bars indicate the total uncertainties.

Detection efficiency per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> mumunu. The error bars indicate the total uncertainties.

Overall signal efficiency (acceptance times efficiency) per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> eenu. The error bars indicate the total uncertainties.

Overall signal efficiency (acceptance times efficiency) per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> emunu. The error bars indicate the total uncertainties.

Overall signal efficiency (acceptance times efficiency) per decay as a function of the mean proper lifetime (ctau) of the neutralino for neutralino -> mumunu. The error bars indicate the total uncertainties.

Detection efficiency per decay for Rxy < 22 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for 22 <= Rxy < 38 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for 38 <= Rxy < 73 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for 73 <= Rxy < 111 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for 111 <= Rxy < 145 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for 145 <= Rxy < 300 mm as a function of the invariant mass and pT of the electron pair in LLP -> eeX.

Detection efficiency per decay for Rxy < 22 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for 22 <= Rxy < 38 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for 38 <= Rxy < 73 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for 73 <= Rxy < 111 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for 111 <= Rxy < 145 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for 145 <= Rxy < 300 mm as a function of the invariant mass and pT of the electron and muon pair in LLP -> emuX.

Detection efficiency per decay for Rxy < 22 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay for 22 <= Rxy < 38 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay for 38 <= Rxy < 73 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay for 73 <= Rxy < 111 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay for 111 <= Rxy < 145 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay for 145 <= Rxy < 300 mm as a function of the invariant mass and pT of the muon pair in LLP -> mumuX.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 100 GeV and LLP -> ee.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 100 GeV and LLP -> emu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 100 GeV and LLP -> mumu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 250 GeV and LLP -> ee.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 250 GeV and LLP -> emu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 250 GeV and LLP -> mumu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 500 GeV and LLP -> ee.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 500 GeV and LLP -> emu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 500 GeV and LLP -> mumu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 750 GeV and LLP -> ee.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 750 GeV and LLP -> emu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 750 GeV and LLP -> mumu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 1000 GeV and LLP -> ee.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 1000 GeV and LLP -> emu.

Detection efficiency per decay as a function of the transverse decay radius Rxy and the dilepton pT for a LLP mass of 1000 GeV and LLP -> mumu.

Acceptance, detection efficiency, and overall signal efficiency in the Z' toy model for mZ' = 100 GeV, three mean proper lifetimes (ctau) and the three decay modes of the Z'.

Acceptance, detection efficiency, and overall signal efficiency in the Z' toy model for mZ' = 250 GeV, three mean proper lifetimes (ctau) and the three decay modes of the Z'.

Acceptance, detection efficiency, and overall signal efficiency in the Z' toy model for mZ' = 500 GeV, three mean proper lifetimes (ctau) and the three decay modes of the Z'.

Acceptance, detection efficiency, and overall signal efficiency in the Z' toy model for mZ' = 750 GeV, three mean proper lifetimes (ctau) and the three decay modes of the Z'.

Acceptance, detection efficiency, and overall signal efficiency in the Z' toy model for mZ' = 1000 GeV, three mean proper lifetimes (ctau) and the three decay modes of the Z'.

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

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

Observed upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the electron+neutrino channel as a function of the $W^\prime$ pole mass.

Observed upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the muon+neutrino channel as a function of the $W^\prime$ pole mass.

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

Observed upper limits at the 95% CL on the visible cross section in the electron+neutrino channel as a function of the transverse mass threshold.

Observed upper limits at the 95% CL on the visible cross section in the muon+neutrino channel as a function of the transverse mass threshold.

Product of acceptance and efficiency for the electron and muon selections as a function of the $W^\prime$ pole mass.

Expected upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the electron+neutrino channel as a function of the $W^\prime$ pole mass.

Expected upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the muon+neutrino channel as a function of the $W^\prime$ pole mass.

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

Observed 95% confidence-level upper limits on the cross section for Drell-Yan spin-1/2 HECO production as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 4000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 200 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 1000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 1500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 2000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 2500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 3000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 4000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 200 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 1000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 1500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 2000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 2500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 3000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 4000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 200 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 1000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 1500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 2000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 2500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 3000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 4000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 200 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 1000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 1500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 2000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 2500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 3000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 4000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 200 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 1000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 1500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 2000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 2500 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 3000 GeV.

Total selection efficiency (i.e., the fraction of MC HECOs surviving the trigger and offline selection criteria) as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 4000 GeV.

Displaced HNL search observed 95% confidence level exclusion contour in $|U_{\mu}|^2$ as a function of HNL mass (LNV case).

Prompt HNL search observed and expected 95% confidence level exclusion upper limits in $|U_{e}|^2$ as a function of HNL mass.

Prompt HNL search observed and expected 95% confidence level exclusion upper limits in $|U_{\mu}|^2$ as a function of HNL mass.

Displaced HNL search expected 95% confidence level exclusion contour in $|U_{\mu}|^2$ as a function of HNL mass (LNC case).

Displaced HNL search expected 95% confidence level exclusion contour in $|U_{\mu}|^2$ as a function of HNL mass (LNV case).