Distributions in SR-low of the data and post-fit background prediction for H^boost. The SR-low event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-low of the data and post-fit background prediction for m_eff^3ℓ/H^boost. The SR-low event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-low of the data and post-fit background prediction for m_eff^3ℓ/H^boost. The SR-low event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-low of the data and post-fit background prediction for p_T^soft/(p_T^soft + m_eff^3ℓ). The SR-low event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-low of the data and post-fit background prediction for p_T^soft/(p_T^soft + m_eff^3ℓ). The SR-low event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for m_T. The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for m_T. The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for R(E_T^miss,jets). The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for R(E_T^miss,jets). The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for p_T^soft. The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for p_T^soft. The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for p_T^jets. The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Distributions in SR-ISR of the data and post-fit background prediction for p_T^jets. The SR-ISR event selections are applied for each distribution except for the variable shown, where the selection is indicated by a red arrow. The normalization factor for the WZ background is derived from the background-only estimation described in Section 7. The expected distribution for a benchmark signal model is included for comparison. The first (last) bin includes underflow (overflow). The "Top-quark like" category contains the tt̄, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. The bottom panel shows the ratio of the data to the post-fit background prediction. The hatched bands indicate the combined theoretical, experimental, and MC statistical uncertainties.

Observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Plus 1$\sigma$ uncertainty, varying the signal cross section within its uncertainty, on the observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Plus 1$\sigma$ uncertainty, varying the signal cross section within its uncertainty, on the observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Minus 1$\sigma$ uncertainty, varying the signal cross section within its uncertainty, on the observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Minus 1$\sigma$ uncertainty, varying the signal cross section within its uncertainty, on the observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Plus 1$\sigma$ uncertainty, due to uncertainties in the background prediction and experimental uncertainties affecting the signal, on the expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Plus 1$\sigma$ uncertainty, due to uncertainties in the background prediction and experimental uncertainties affecting the signal, on the expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Minus 1$\sigma$ uncertainty, due to uncertainties in the background prediction and experimental uncertainties affecting the signal, on the expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Minus 1$\sigma$ uncertainty, due to uncertainties in the background prediction and experimental uncertainties affecting the signal, on the expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Upper limits on observed wino-bino simplified model signal cross section $\sigma_\text{obs}^\text{95}$.

Upper limits on observed wino-bino simplified model signal cross section $\sigma_\text{obs}^\text{95}$.

Upper limits on expected wino-bino simplified model signal cross section $\sigma_\text{exp}^\text{95}$.

Upper limits on expected wino-bino simplified model signal cross section $\sigma_\text{exp}^\text{95}$.

Signal acceptance in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal acceptance in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal efficiency in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal efficiency in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal acceptance in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal acceptance in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal efficiency in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal efficiency in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal acceptance in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal acceptance in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal efficiency in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal efficiency in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal acceptance in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal acceptance in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal efficiency in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal efficiency in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

The observed and expected yields after the background-only fit in the SRs. The normalization factors of the $WZ$ sample for the low-mass and ISR regions are different and are treated separately in the combined fit. \The "Top-quark like" category contains the tt&#772;, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. Combined statistical and systematic uncertainties are presented. The individual uncertainties can be correlated and do not necessarily add in quadrature to equal the total background uncertainty.

The observed and expected yields after the background-only fit in the SRs. The normalization factors of the $WZ$ sample for the low-mass and ISR regions are different and are treated separately in the combined fit. \The "Top-quark like" category contains the tt&#772;, Wt, and WW processes while the "Others" category contains backgrounds from triboson production and processes that include a Higgs boson, 3 or more tops, and tops produced in association with W or Z bosons. Combined statistical and systematic uncertainties are presented. The individual uncertainties can be correlated and do not necessarily add in quadrature to equal the total background uncertainty.

Summary of the expected background and data yields in $\text{SR-low}$ and $\text{SR-ISR}$. The second and third columns show the data and total expected background with systematic uncertainties. The fourth column gives the model-independent upper limits at 95\% CL on the visible cross section ($\sigma_\text{vis}$). The fifth and sixth columns give the visible number of observed ($S^{95}_\text{obs}$) and expected ($S^{95}_\text{exp}$) events of a generic beyond-the-SM process, where uncertainties on $S^{95}_\text{exp}$ reflect the $\pm 1 \sigma$ uncertainties on the background estimation. The last column shows the discovery $p$-value and Gaussian significance $Z$ assuming no signal.

Summary of the expected background and data yields in $\text{SR-low}$ and $\text{SR-ISR}$. The second and third columns show the data and total expected background with systematic uncertainties. The fourth column gives the model-independent upper limits at 95\% CL on the visible cross section ($\sigma_\text{vis}$). The fifth and sixth columns give the visible number of observed ($S^{95}_\text{obs}$) and expected ($S^{95}_\text{exp}$) events of a generic beyond-the-SM process, where uncertainties on $S^{95}_\text{exp}$ reflect the $\pm 1 \sigma$ uncertainties on the background estimation. The last column shows the discovery $p$-value and Gaussian significance $Z$ assuming no signal.

Upper limits on observed (expected) wino-bino simplified model signal cross section $\sigma_\text{obs(exp)}^\text{95}$.

Upper limits on observed (expected) wino-bino simplified model signal cross section $\sigma_\text{obs(exp)}^\text{95}$.

Full list of event selections and MC generator-weighted yields and in $\text{SR-ISR}$ for the main $WZ$ background and a representative $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ signal point of mass 200 GeV and mass splitting $\Delta m = 100$ GeV normalized to 139 fb$^{-1}$. 40000 events were generated.

Full list of event selections and MC generator-weighted yields and in $\text{SR-ISR}$ for the main $WZ$ background and a representative $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ signal point of mass 200 GeV and mass splitting $\Delta m = 100$ GeV normalized to 139 fb$^{-1}$. 40000 events were generated.

Full list of event selections and MC generator-weighted yields and in $\text{SR-low}$ for the main $WZ$ background and a representative $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ signal point of mass 200 GeV and mass splitting $\Delta m = 100$ GeV normalized to 139 fb$^{-1}$. 40000 events were generated.

Full list of event selections and MC generator-weighted yields and in $\text{SR-low}$ for the main $WZ$ background and a representative $\tilde{\chi}^{\pm}_{1}\tilde{\chi}^{0}_{2}$ signal point of mass 200 GeV and mass splitting $\Delta m = 100$ GeV normalized to 139 fb$^{-1}$. 40000 events were generated.

The observed upper limits on the model cross-section in units of pb for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production. Three points at $M({\tilde{\chi}}^{0}_{1})=200GeV$ were removed from the plot but kept in the table because they overlapped with the plot's legend and are far from the exclusion contour.

The observed 95\% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with combined ${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production.

The observed 95\% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with combined ${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production.

The expected 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with combined ${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production.

The expected 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with combined ${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production.

The observed 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production.

The observed 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production.

The expected 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production.

The expected 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production.

Observed 95% CL exclusion limits for simplified models with direct stau pair production in SR-lowMass.

Observed 95% CL exclusion limits for simplified models with direct stau pair production in SR-lowMass.

Expected 95% CL exclusion limits for simplified models with direct stau pair production in SR-lowMass.

Expected 95% CL exclusion limits for simplified models with direct stau pair production in SR-lowMass.

Observed 95% CL exclusion limits for simplified models with direct stau pair production in SR-highMass.

Observed 95% CL exclusion limits for simplified models with direct stau pair production in SR-highMass.

Expected 95% CL exclusion limits for simplified models with direct stau pair production in SR-highMass.

Expected 95% CL exclusion limits for simplified models with direct stau pair production in SR-highMass.

Signal acceptance in SR highMass for combined stau final states

Signal acceptance in SR highMass for combined stau final states

Signal acceptance in SR lowMass for combined stau final states

Signal acceptance in SR lowMass for combined stau final states

Signal efficiency in SR highMass for combined stau final states

Signal efficiency in SR highMass for combined stau final states

Signal efficiency in SR lowMass for combined stau final states

Signal efficiency in SR lowMass for combined stau final states

Signal acceptance*efficiency in SR highMass for combined stau final states

Signal acceptance*efficiency in SR highMass for combined stau final states

Signal acceptance*efficiency in SR lowMass for combined stau final states

Signal acceptance*efficiency in SR lowMass for combined stau final states

Cutflow for two reference points (${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production) in SR. The column labelled $N_{weighted}$ shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$, while $N_{raw}$ in brackets shows the results for the generated number of events. The quoted uncertainties are statistical only. The "Generator filter" includes the requirements that two $\tau$ in the event have ${p}_{T} > 15$ GeV and $|\eta| <$ 2.6. The "Baseline Cut" includes the requirement of two baseline $\tau$ with a minimum value at 0.01 of the boosted decision tree discriminant (JetBDTSigTransMin $>$ 0.01) and ${p}_{T, \tau_{1}} > 50$ GeV and ${p}_{T, \tau_{2}} > 40$ GeV. At the step "Trigger & offline cuts", the following requirements are applied: the event is recorded using the asymmetric di-$\tau$ trigger (di-$\tau$ $E_{T}^{miss}$ trigger) in SR-lowMass (SR-highMass), and the lepton $p_{T}$ and $E_{T}^{miss}$ are required at plateau.

Cutflow for two reference points (${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production) in SR. The column labelled $N_{weighted}$ shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$, while $N_{raw}$ in brackets shows the results for the generated number of events. The quoted uncertainties are statistical only. The "Generator filter" includes the requirements that two $\tau$ in the event have ${p}_{T} > 15$ GeV and $|\eta| <$ 2.6. The "Baseline Cut" includes the requirement of two baseline $\tau$ with a minimum value at 0.01 of the boosted decision tree discriminant (JetBDTSigTransMin $>$ 0.01) and ${p}_{T, \tau_{1}} > 50$ GeV and ${p}_{T, \tau_{2}} > 40$ GeV. At the step "Trigger & offline cuts", the following requirements are applied: the event is recorded using the asymmetric di-$\tau$ trigger (di-$\tau$ $E_{T}^{miss}$ trigger) in SR-lowMass (SR-highMass), and the lepton $p_{T}$ and $E_{T}^{miss}$ are required at plateau.

Observed and expected numbers of events in the control and signal regions where all control and signal region bins are included as constraints in the likelihood. The expected event yields of SM processes are given after the background-only fit. The entries marked as "--" are negligible. The uncertainties correspond to the sum in quadrature of statistical and systematic uncertainties. The correlation of systematic uncertainties among control regions and among background processes is fully taken into account.

Observed and expected numbers of events in the control and signal regions where all control and signal region bins are included as constraints in the likelihood. The expected event yields of SM processes are given after the background-only fit. The entries marked as "--" are negligible. The uncertainties correspond to the sum in quadrature of statistical and systematic uncertainties. The correlation of systematic uncertainties among control regions and among background processes is fully taken into account.

The post-fit $m_{T2}$ distribution for SR-lowMass. The stacked histograms show the expected SM backgrounds. The multi-jet contribution is estimated from data using the ABCD method. The contributions of multi-jet and $W$+jets events are scaled with the corresponding normalization factors derived from the background-only fit. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distributions from the SUSY reference points are also shown as dashed lines. The last bin includes the overflow events.

The post-fit $m_{T2}$ distribution for SR-lowMass. The stacked histograms show the expected SM backgrounds. The multi-jet contribution is estimated from data using the ABCD method. The contributions of multi-jet and $W$+jets events are scaled with the corresponding normalization factors derived from the background-only fit. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distributions from the SUSY reference points are also shown as dashed lines. The last bin includes the overflow events.

The post-fit $m_{T2}$ distribution for SR-highMass. The stacked histograms show the expected SM backgrounds. The multi-jet contribution is estimated from data using the ABCD method. The contributions of multi-jet and $W$+jets events are scaled with the corresponding normalization factors derived from the background-only fit. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distributions from the SUSY reference points are also shown as dashed lines. The last bin includes the overflow events.

The post-fit $m_{T2}$ distribution for SR-highMass. The stacked histograms show the expected SM backgrounds. The multi-jet contribution is estimated from data using the ABCD method. The contributions of multi-jet and $W$+jets events are scaled with the corresponding normalization factors derived from the background-only fit. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distributions from the SUSY reference points are also shown as dashed lines. The last bin includes the overflow events.

The $m_{T2}$ post-fit distributions in the multi-jet background validation region VR-F (lowMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $m_{T2}$ post-fit distributions in the multi-jet background validation region VR-F (lowMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $m_{T2}$ post-fit distributions in the multi-jet background validation VR-F (highMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $m_{T2}$ post-fit distributions in the multi-jet background validation VR-F (highMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $E_{T}^{miss}$ post-fit distributions in the multi-jet background validation region VR-F (lowMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $E_{T}^{miss}$ post-fit distributions in the multi-jet background validation region VR-F (lowMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $E_{T}^{miss}$ post-fit distributions in the multi-jet background validation region VR-F (highMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $E_{T}^{miss}$ post-fit distributions in the multi-jet background validation region VR-F (highMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The pre-fit $m_{T2}$ distribution in the $WCR$. The SM backgrounds other than multi-jet production are estimated from MC simulation. The multi-jet contribution is estimated from data using the OS--SS method. The hatched bands represent the combined statistical and systematic uncertainties of the total SM background. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The pre-fit $m_{T2}$ distribution in the $WCR$. The SM backgrounds other than multi-jet production are estimated from MC simulation. The multi-jet contribution is estimated from data using the OS--SS method. The hatched bands represent the combined statistical and systematic uncertainties of the total SM background. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The post-fit yields in the $WVR$, $TVRs$, $ZVRs$ and $VVVRs$. The SM backgrounds other than multi-jet production are estimated from MC simulation. The multi-jet contribution is negligible and is estimated from data using the ABCD method, using CRs obtained with the same technique used for the SRs. The hatched bands represent the combined statistical and systematic uncertainties of the total SM background. For illustration, the distributions of the SUSY reference points are also shown as dashed lines. The lower panels show the ratio of data to the SM background estimate.

The post-fit yields in the $WVR$, $TVRs$, $ZVRs$ and $VVVRs$. The SM backgrounds other than multi-jet production are estimated from MC simulation. The multi-jet contribution is negligible and is estimated from data using the ABCD method, using CRs obtained with the same technique used for the SRs. The hatched bands represent the combined statistical and systematic uncertainties of the total SM background. For illustration, the distributions of the SUSY reference points are also shown as dashed lines. The lower panels show the ratio of data to the SM background estimate.

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-veto category of the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table.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 predicted mTtot distribution in the b-tag category of the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. 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 predicted mTtot distribution in the b-tag category of the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. 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 predicted mTtot distribution in the b-tag category of the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. 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 predicted mTtot distribution in the b-tag category of the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. 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 predicted mTtot distribution in the b-veto 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 predicted mTtot distribution in the b-veto 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 predicted mTtot distribution in the b-veto 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 predicted mTtot distribution in the b-veto 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 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 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 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 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 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 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 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.

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

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

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

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

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

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

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

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

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

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

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

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

The expected 95% CL upper limits with 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.

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

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

The expected 95% CL upper limits with 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.

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

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

The expected 95% CL upper limits with 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.

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

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

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

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

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

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

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

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

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

The expected 95% CL upper limits with 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.

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

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

The expected 95% CL upper limits with 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.

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

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

The expected 95% CL upper limits with 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.

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 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 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 gluon-gluon fusion as a function of the scalar boson mass.

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

Acceptance times efficiency for a scalar boson produced by b-associated production 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.

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.

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.

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.

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.

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.

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.

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

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

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 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 200 GeV signal mass point is shown in the HEPData table.

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

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 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 250 GeV signal mass point is shown in the HEPData table.

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

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 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 300 GeV signal mass point is shown in the HEPData table.

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

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 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 350 GeV signal mass point is shown in the HEPData table.

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

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 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 400 GeV signal mass point is shown in the HEPData table.

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

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 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 500 GeV signal mass point is shown in the HEPData table.

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

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 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 600 GeV signal mass point is shown in the HEPData table.

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

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 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 700 GeV signal mass point is shown in the HEPData table.

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

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 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 800 GeV signal mass point is shown in the HEPData table.

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

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 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 1000 GeV signal mass point is shown in the HEPData table.

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

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 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 1200 GeV signal mass point is shown in the HEPData table.

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

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 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 1500 GeV signal mass point is shown in the HEPData table.

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

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 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 2000 GeV signal mass point is shown in the HEPData table.

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

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.

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.

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

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.

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.

Summary of results for bottom muon v2 as a function of pT for different centrality. Uncertainties are statistical and systematic, respectively.

Summary of results for charm muon v3 as a function of pT for different centrality. Uncertainties are statistical and systematic, respectively.

Summary of results for bottom muon v3 as a function of pT for different centrality. Uncertainties are statistical and systematic, respectively.

Expected and observed 95% CL lower limits on the leptoquark mass as a function of the branching ratio $\mathcal{B}$(LQ->$t\mu$).

The best-fit expected and observed distributions of the combined NN output in the VRSS for the $e\tau$ channel for events with 3-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $e\tau$ channel for events with 1-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $e\tau$ channel for events with 3-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $\mu\tau$ channel for events with 1-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

The best-fit expected and observed distributions of the combined NN output in the SR for the $\mu\tau$ channel for events with 3-prong $\tau_\text{had-vis}$ candidates. The last bin in each plot includes overflow events.

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

The absolute differential cross-section measured in the fiducial phase-space as a function of the minimum $\Delta R$ between the photon and the leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute differential cross-section measured in the fiducial phase-space as a function of the $\Delta\phi$ between the two leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute differential cross-section measured in the fiducial phase-space as a function of the $|\Delta\eta|$ between the two leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised differential cross-section measured in the fiducial phase-space as a function of the photon pT in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised differential cross-section measured in the fiducial phase-space as a function of the photon $|\eta|$ in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised differential cross-section measured in the fiducial phase-space as a function of the minimum $\Delta R$ between the photon and the leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised differential cross-section measured in the fiducial phase-space as a function of the $\Delta\phi$ between the two leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised differential cross-section measured in the fiducial phase-space as a function of the $|\Delta\eta|$ between the two leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The total correlation matrix of the absolute differential cross-section measured in the fiducial phase-space as a function of the photon pT in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the absolute differential cross-section measured in the fiducial phase-space as a function of the photon $|\eta|$ in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the absolute differential cross-section measured in the fiducial phase-space as a function of the minimum $\Delta R$ between the photon and the leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the absolute differential cross-section measured in the fiducial phase-space as a function of the $\Delta\phi$ between the two leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the absolute differential cross-section measured in the fiducial phase-space as a function of the $|\Delta\eta|$ between the two leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the normalised differential cross-section measured in the fiducial phase-space as a function of the photon pT in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the normalised differential cross-section measured in the fiducial phase-space as a function of the photon $|\eta|$ in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the normalised differential cross-section measured in the fiducial phase-space as a function of the minimum $\Delta R$ between the photon and the leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the normalised differential cross-section measured in the fiducial phase-space as a function of the $\Delta\phi$ between the two leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the normalised differential cross-section measured in the fiducial phase-space as a function of the $|\Delta\eta|$ between the two leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The statistical correlation matrix of all the absolute differential cross-sections measured in the fiducial phase-space in the electron-muon channel.

The statistical correlation matrix of all the normalised differential cross-sections measured in the fiducial phase-space in the electron-muon channel.

Fiducial region definition.

Summary of results for charm and bottom muon v2 as a function of pT. Uncertainties are statistical and systematic, respectively.