Exclusive J/psi photoproduction cross section in gamma-p.

Dissociative J/psi photoproduction gamma-p.

Cross section ratio, dissociative J/psi photoproduction gamma-p over exclusive J/psi photoproduction gamma-p.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $0 < \Delta\Phi_{1,2} < 150^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $150 < \Delta\Phi_{1,2} < 170^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $170 < \Delta\Phi_{1,2} < 180^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $0 < \Delta\Phi_{1,2} < 150^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $150 < \Delta\Phi_{1,2} < 170^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $170 < \Delta\Phi_{1,2} < 180^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the leading $p_{T}$ jet ($p_{T1}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the subleading $p_{T}$ jet ($p_{T2}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the third leading $p_{T}$ jet ($p_{T3}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the fourth leading $p_{T}$ jet ($p_{T4}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity distribution in TUnfold binning ($N_{jets},\Delta\phi_{1,2},p_{T1}$) as indicated in the XML file provided as additional resource. The uncertainties follow the notation of Table 1.

Correlation matrix at particle level for the measured jet multiplicity in TUnfold binning ($N_{jets},\Delta\phi_{1,2},p_{T1}$) as indicated in the XML file provided as additional resource.

Jet $p_{T}$ distributions in TUnfold binning ($p_{T1},p_{T2},p_{T3},p_{T4}$) as indicated in the XML file provided as additional resource. The uncertainties follow the notation of Table 1.

Correlation matrix at particle level for the measured jet $p_{T}$ distributions in TUnfold binning ($p_{T1},p_{T2},p_{T3},p_{T4}$) as indicated in the XML file provided as additional resource.

Observed yields in the boost_0_ge2J signal region category of the hh channel.

Observed yields in the boost_1_1J signal region category of the hh channel.

Observed yields in the boost_1_1J signal region category of the hh channel.

Observed yields in the boost_1_ge2J signal region category of the hh channel.

Observed yields in the boost_1_ge2J signal region category of the hh channel.

Observed yields in the boost_2 signal region category of the hh channel.

Observed yields in the boost_2 signal region category of the hh channel.

Observed yields in the boost_3 signal region category of the hh channel.

Observed yields in the boost_3 signal region category of the hh channel.

Observed yields in the ttH_0 signal region category of the hh channel.

Observed yields in the ttH_0 signal region category of the hh channel.

Observed yields in the ttH_1 signal region category of the hh channel.

Observed yields in the ttH_1 signal region category of the hh channel.

Observed yields in the VH_0 signal region category of the hh channel.

Observed yields in the VH_0 signal region category of the hh channel.

Observed yields in the VH_1 signal region category of the hh channel.

Observed yields in the VH_1 signal region category of the hh channel.

Observed yields in the VBF_0 signal region category of the hh channel.

Observed yields in the VBF_0 signal region category of the hh channel.

Observed yields in the VBF_1 signal region category of the hh channel.

Observed yields in the VBF_1 signal region category of the hh channel.

Observed yields in the boost_0_1J signal region category of the lh channel.

Observed yields in the boost_0_1J signal region category of the lh channel.

Observed yields in the boost_0_ge2J signal region category of the lh channel.

Observed yields in the boost_0_ge2J signal region category of the lh channel.

Observed yields in the boost_1_1J signal region category of the lh channel.

Observed yields in the boost_1_1J signal region category of the lh channel.

Observed yields in the boost_1_ge2J signal region category of the lh channel.

Observed yields in the boost_1_ge2J signal region category of the lh channel.

Observed yields in the boost_2 signal region category of the lh channel.

Observed yields in the boost_2 signal region category of the lh channel.

Observed yields in the boost_3 signal region category of the lh channel.

Observed yields in the boost_3 signal region category of the lh channel.

Observed yields in the VH_0 signal region category of the lh channel.

Observed yields in the VH_0 signal region category of the lh channel.

Observed yields in the VH_1 signal region category of the lh channel.

Observed yields in the VH_1 signal region category of the lh channel.

Observed yields in the VBF_0 signal region category of the lh channel.

Observed yields in the VBF_0 signal region category of the lh channel.

Observed yields in the VBF_1 signal region category of the lh channel.

Observed yields in the VBF_1 signal region category of the lh channel.

Observed yields in the boost_0_1J signal region category of the ll channel.

Observed yields in the boost_0_1J signal region category of the ll channel.

Observed yields in the boost_0_ge2J signal region category of the ll channel.

Observed yields in the boost_0_ge2J signal region category of the ll channel.

Observed yields in the boost_1_1J signal region category of the ll channel.

Observed yields in the boost_1_1J signal region category of the ll channel.

Observed yields in the boost_1_ge2J signal region category of the ll channel.

Observed yields in the boost_1_ge2J signal region category of the ll channel.

Observed yields in the boost_2 signal region category of the ll channel.

Observed yields in the boost_2 signal region category of the ll channel.

Observed yields in the boost_3 signal region category of the ll channel.

Observed yields in the boost_3 signal region category of the ll channel.

Observed yields in the VH_0 signal region category of the ll channel.

Observed yields in the VH_0 signal region category of the ll channel.

Observed yields in the VH_1 signal region category of the ll channel.

Observed yields in the VH_1 signal region category of the ll channel.

Observed yields in the VBF_0 signal region category of the ll channel.

Observed yields in the VBF_0 signal region category of the ll channel.

Observed yields in the VBF_1 signal region category of the ll channel.

Observed yields in the VBF_1 signal region category of the ll channel.

Best-fit values and uncertainties for the pp $\rightarrow H\rightarrow\tau\tau$ cross-section measurement and the measurement in the four dominant production modes. All measurements include the branching ratio of $H\rightarrow\tau\tau$ and are performed with true Higgs boson rapidity $|y_{H}|<2.5$. The SM predictions for each region, computed using the inclusive cross-section calculations and the simulated event samples are also shown. The contributions to the total uncertainty in the measurements from statistical (Stat. unc.) or systematic uncertainties (Syst. unc.) are given separately.

Best-fit values and uncertainties for the pp $\rightarrow H\rightarrow\tau\tau$ cross-section measurement and the measurement in the four dominant production modes. All measurements include the branching ratio of $H\rightarrow\tau\tau$ and are performed with true Higgs boson rapidity $|y_{H}|<2.5$. The SM predictions for each region, computed using the inclusive cross-section calculations and the simulated event samples are also shown. The contributions to the total uncertainty in the measurements from statistical (Stat. unc.) or systematic uncertainties (Syst. unc.) are given separately.

The measured correlations between each parameter of interest in the measurement of the cross-sections per production mode

The measured correlations between each parameter of interest in the measurement of the cross-sections per production mode

Best-fit values and uncertainties for the $H\rightarrow\tau\tau$ cross-sections, in the reduced stage 1.2 STXS scheme. The EW production mode includes vector-boson fusion and qq$\rightarrow$V($\rightarrow$qq)H processes. All measurements include the branching ratio of $H\rightarrow\tau\tau$ and are performed with true Higgs boson rapidity $|y_{H}|<2.5$. The SM predictions for each region, computed using the inclusive cross-section calculations and the simulated event samples are also shown. The contributions to the total uncertainty in the measurements from statistical (Stat. unc.) or systematic uncertainties (Syst. unc.) are given separately. The asterisk symbol (*) indicates that the criteria for m$_{jj}$ only apply to events with at least two reconstructed jets.

Best-fit values and uncertainties for the $H\rightarrow\tau\tau$ cross-sections, in the reduced stage 1.2 STXS scheme. The EW production mode includes vector-boson fusion and qq$\rightarrow$V($\rightarrow$qq)H processes. All measurements include the branching ratio of $H\rightarrow\tau\tau$ and are performed with true Higgs boson rapidity $|y_{H}|<2.5$. The SM predictions for each region, computed using the inclusive cross-section calculations and the simulated event samples are also shown. The contributions to the total uncertainty in the measurements from statistical (Stat. unc.) or systematic uncertainties (Syst. unc.) are given separately. The asterisk symbol (*) indicates that the criteria for m$_{jj}$ only apply to events with at least two reconstructed jets.

The measured correlations between each pair of parameters of interest in the STXS measurement. The asterisk symbol (*) indicates that the criteria for mjj only apply to events with at least two reconstructed jets.

The measured correlations between each pair of parameters of interest in the STXS measurement. The asterisk symbol (*) indicates that the criteria for mjj only apply to events with at least two reconstructed jets.

The reconstructed m$_{tb}$ distributions in data and expected background in validation region for the data taking period of 2017. Yield in each bin is divided by the corresponding bin width.

The reconstructed m$_{tb}$ distributions in data and expected background in signal region for the data taking period of 2018. Yield in each bin is divided by the corresponding bin width.

The reconstructed m$_{tb}$ distributions in data and expected background in validation region for the data taking period of 2018. Yield in each bin is divided by the corresponding bin width.

Upper limits at 95% CL on the cross section for the production of right-handed W' boson using data and backgrounds in three years combined. Theory prediction is also shown.

Upper limits at 95% CL on the cross section for the production of left-handed W' boson using data and backgrounds in three years combined. Theory prediction is also shown.

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the rapidity of the pair. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $K^+K^-$ pairs as a function of the rapidity of the pair. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $K^+$, $K^-$ - $p_{\mathrm{T}} > 0.3~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(K^+), p_{\mathrm{T}}(K^-)) < 0.7~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $p\bar{p}$ pairs as a function of the rapidity of the pair. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $p$, $\bar{p}$ - $p_{\mathrm{T}} > 0.4~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(p), p_{\mathrm{T}}(\bar{p})) < 1.1~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the difference of azimuthal angles of the forward scattered protons. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $K^+K^-$ pairs as a function of the difference of azimuthal angles of the forward scattered protons. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $K^+$, $K^-$ - $p_{\mathrm{T}} > 0.3~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(K^+), p_{\mathrm{T}}(K^-)) < 0.7~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $p\bar{p}$ pairs as a function of the difference of azimuthal angles of the forward scattered protons. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $p$, $\bar{p}$ - $p_{\mathrm{T}} > 0.4~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(p), p_{\mathrm{T}}(\bar{p})) < 1.1~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the sum of the squares of the four-momenta losses in the proton vertices. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $K^+K^-$ pairs as a function of the sum of the squares of the four-momenta losses in the proton vertices. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $K^+$, $K^-$ - $p_{\mathrm{T}} > 0.3~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(K^+), p_{\mathrm{T}}(K^-)) < 0.7~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $p\bar{p}$ pairs as a function of the sum of the squares of the four-momenta losses in the proton vertices. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $p$, $\bar{p}$ - $p_{\mathrm{T}} > 0.4~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(p), p_{\mathrm{T}}(\bar{p})) < 1.1~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the invariant mass of the pair, for events with $\Delta\varphi < 90$ degrees. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $\Delta\varphi < 90$ degrees

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the invariant mass of the pair, for events with $\Delta\varphi > 90$ degrees. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $\Delta\varphi > 90$ degrees

Differential fiducial cross section for CEP of $K^+K^-$ pairs as a function of the invariant mass of the pair, for events with $\Delta\varphi < 90$ degrees. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $K^+$, $K^-$ - $p_{\mathrm{T}} > 0.3~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(K^+), p_{\mathrm{T}}(K^-)) < 0.7~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $\Delta\varphi < 90$ degrees

Differential fiducial cross section for CEP of $K^+K^-$ pairs as a function of the invariant mass of the pair, for events with $\Delta\varphi > 90$ degrees. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $K^+$, $K^-$ - $p_{\mathrm{T}} > 0.3~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(K^+), p_{\mathrm{T}}(K^-)) < 0.7~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $\Delta\varphi > 90$ degrees

Differential fiducial cross section for CEP of $p\bar{p}$ pairs as a function of the invariant mass of the pair, for events with $\Delta\varphi < 90$ degrees. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $p$, $\bar{p}$ - $p_{\mathrm{T}} > 0.4~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(p), p_{\mathrm{T}}(\bar{p})) < 1.1~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $\Delta\varphi < 90$ degrees

Differential fiducial cross section for CEP of $p\bar{p}$ pairs as a function of the invariant mass of the pair, for events with $\Delta\varphi > 90$ degrees. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $p$, $\bar{p}$ - $p_{\mathrm{T}} > 0.4~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(p), p_{\mathrm{T}}(\bar{p})) < 1.1~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $\Delta\varphi > 90$ degrees

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the invariant mass of the pair, for events with $\Delta\varphi < 90$ degrees and $\Delta p'_{\mathrm{T}} < 0.12~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $\Delta\varphi < 90$ degrees - $\Delta p'_{\mathrm{T}} < 0.12~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the invariant mass of the pair, for events with $\Delta\varphi < 90$ degrees and $\Delta p'_{\mathrm{T}} > 0.12~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $\Delta\varphi < 90$ degrees - $\Delta p'_{\mathrm{T}} > 0.12~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the $\cos{\theta^{CS}}$ of the positive charge pion in the Collins-Soper reference frame. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $K^+K^-$ pairs as a function of the $\cos{\theta^{CS}}$ of the positive charge kaon in the Collins-Soper reference frame. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $K^+$, $K^-$ - $p_{\mathrm{T}} > 0.3~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(K^+), p_{\mathrm{T}}(K^-)) < 0.7~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $p\bar{p}$ pairs as a function of the $\cos{\theta^{CS}}$ of the central state proton in the Collins-Soper reference frame. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $p$, $\bar{p}$ - $p_{\mathrm{T}} > 0.4~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(p), p_{\mathrm{T}}(\bar{p})) < 1.1~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the $\phi^{CS}$ of the positive charge pion in the Collins-Soper reference frame. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $K^+K^-$ pairs as a function of the $\phi^{CS}$ of the positive charge kaon in the Collins-Soper reference frame. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $K^+$, $K^-$ - $p_{\mathrm{T}} > 0.3~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(K^+), p_{\mathrm{T}}(K^-)) < 0.7~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $p\bar{p}$ pairs as a function of the $\phi^{CS}$ of the central state proton in the Collins-Soper reference frame. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $p$, $\bar{p}$ - $p_{\mathrm{T}} > 0.4~\mathrm{GeV}$ - $min(p_{\mathrm{T}}(p), p_{\mathrm{T}}(\bar{p})) < 1.1~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the rapidity of the pair, for invariant masses $m(\pi^+\pi^-) < 1~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $m(\pi^+\pi^-) < 1~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the difference of azimuthal angles of the forward scattered protons, for invariant masses $m(\pi^+\pi^-) < 1~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $m(\pi^+\pi^-) < 1~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the sum of the squares of the four-momenta losses in the proton vertices, for invariant masses $m(\pi^+\pi^-) < 1~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $m(\pi^+\pi^-) < 1~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the rapidity of the pair, for invariant masses $1 \mathrm{GeV} < m(\pi^+\pi^-) < 1.5~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $1 \mathrm{GeV} < m(\pi^+\pi^-) < 1.5~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the difference of azimuthal angles of the forward scattered protons, for invariant masses $1 \mathrm{GeV} < m(\pi^+\pi^-) < 1.5~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $1 \mathrm{GeV} < m(\pi^+\pi^-) < 1.5~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the sum of the squares of the four-momenta losses in the proton vertices, for invariant masses $1 \mathrm{GeV} < m(\pi^+\pi^-) < 1.5~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $1 \mathrm{GeV} < m(\pi^+\pi^-) < 1.5~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the rapidity of the pair, for invariant masses $m(\pi^+\pi^-) > 1.5~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $m(\pi^+\pi^-) > 1.5~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the difference of azimuthal angles of the forward scattered protons, for invariant masses $m(\pi^+\pi^-) > 1.5~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $m(\pi^+\pi^-) > 1.5~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the sum of the squares of the four-momenta losses in the proton vertices, for invariant masses $m(\pi^+\pi^-) > 1.5~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $m(\pi^+\pi^-) > 1.5~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the $\cos{\theta^{CS}}$ of the positive charge pion in the Collins-Soper reference frame, for invariant masses $m(\pi^+\pi^-) < 1~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $m(\pi^+\pi^-) < 1~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the $\cos{\theta^{CS}}$ of the positive charge pion in the Collins-Soper reference frame, for invariant masses $1 \mathrm{GeV} < m(\pi^+\pi^-) < 1.5~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $1 \mathrm{GeV} < m(\pi^+\pi^-) < 1.5~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the $\cos{\theta^{CS}}$ of the positive charge pion in the Collins-Soper reference frame, for invariant masses $m(\pi^+\pi^-) > 1.5~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $m(\pi^+\pi^-) > 1.5~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the $\phi^{CS}$ of the positive charge pion in the Collins-Soper reference frame, for invariant masses $m(\pi^+\pi^-) < 1~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $m(\pi^+\pi^-) < 1~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the $\phi^{CS}$ of the positive charge pion in the Collins-Soper reference frame, for invariant masses $1 \mathrm{GeV} < m(\pi^+\pi^-) < 1.5~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $1 \mathrm{GeV} < m(\pi^+\pi^-) < 1.5~\mathrm{GeV}$

Differential fiducial cross section for CEP of $\pi^+\pi^-$ pairs as a function of the $\phi^{CS}$ of the positive charge pion in the Collins-Soper reference frame, for invariant masses $m(\pi^+\pi^-) > 1.5~\mathrm{GeV}$. Systematic uncertainties assigned to data points are strongly correlated between bins and should be treated as allowed collective variation of all data points. There are two components of the total systematic uncertainty. The systematic uncertainty related to the experimental tools and analysis method is labeled "syst. (experimental)". The systematic uncertainty related to the integrated luminosity (fully correlated between all data points) is labeled "syst. (luminosity)". Fiducial region definition: * central state $\pi^+$, $\pi^-$ - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ - $|\eta| < 0.7$ * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$ * additional cuts - $m(\pi^+\pi^-) > 1.5~\mathrm{GeV}$

Integrated fiducial cross sections with statistical and systematic uncertainties for CEP of $\pi^+\pi^-$, $K^+K^-$ and $p\bar{p}$ pairs in two ranges of azimuthal angle difference, $\Delta\varphi$, between the two forward-scattered protons. Fiducial region definition: * central state - $p_{\mathrm{T}} > 0.2~\mathrm{GeV}$ ($\pi^+$, $\pi^-$) - $p_{\mathrm{T}} > 0.3~\mathrm{GeV}$, $min(p_{\mathrm{T}}^+, p_{\mathrm{T}}^-) < 0.7~\mathrm{GeV}$ ($K^+$, $K^-$) - $p_{\mathrm{T}} > 0.4~\mathrm{GeV}$, $min(p_{\mathrm{T}}^+, p_{\mathrm{T}}^-) < 1.1~\mathrm{GeV}$ ($p$, $\bar{p}$) - $|\eta| < 0.7$ (all species) * intact forward-scattered beam protons $p'$ - $p_x > -0.2~\mathrm{GeV}$ - $0.2~\mathrm{GeV} < |p_{y}| < 0.4~\mathrm{GeV}$ - $(p_x+0.3~\mathrm{GeV})^2 + p_y^2 < 0.25~\mathrm{GeV}^2$

Results of the fit to extrapolated $d\sigma/dm(\pi^+\pi^-)$ in two ranges of azimuthal angle difference $\Delta\varphi$ between forward-scattered protons. The fit describes the cross-section extrapolated to Lorentz-invariant phase-space defined below: - $|y(\pi^+\pi^-)| < 0.4$ - $0.05~\mathrm{GeV}^2 < -t_1, -t_2 < 0.16~\mathrm{GeV}^2$ The experimental systematic uncertainties "(syst.)" are calculated as the quadratic sum of the differences between the nominal fit result and the result of the fit to $d\sigma/dm(\pi^+\pi^-)$ with each systematic effect. The uncertainties related to the extrapolation "(model.)" are quoted as the largest deviation from the nominal fit result.

Results of the fit to extrapolated $d\sigma/dm(\pi^+\pi^-)$ in two ranges of azimuthal angle difference $\Delta\varphi$ between forward-scattered protons. The fit describes the cross-section extrapolated to Lorentz-invariant phase-space defined below: - $|y(\pi^+\pi^-)| < 0.4$ - $0.05~\mathrm{GeV}^2 < -t_1, -t_2 < 0.16~\mathrm{GeV}^2$ The experimental systematic uncertainties "(syst.)" are calculated as the quadratic sum of the differences between the nominal fit result and the result of the fit to $d\sigma/dm(\pi^+\pi^-)$ with each systematic effect. The uncertainties related to the extrapolation "(model.)" are quoted as the largest deviation from the nominal fit result.

Ratios of integrated cross sections of resonance production to cross sections for f2(1270) production, in the pi+pi- channel. The results are obtained for the Lorentz-invariant phase-space defined below: - $|y(\pi^+\pi^-)| < 0.4$ - $0.05~\mathrm{GeV}^2 < -t_1, -t_2 < 0.16~\mathrm{GeV}^2$ The experimental systematic uncertainties "(syst.)" are calculated as the quadratic sum of the differences between the nominal fit result and the result of the fit to $d\sigma/dm(\pi^+\pi^-)$ with each systematic effect. The uncertainties related to the extrapolation "(model.)" are quoted as the largest deviation from the nominal fit result.

Ratios of integrated cross sections of resonance production in two ranges of $\Delta\varphi$, in the pi+pi- channel. The results are obtained for the Lorentz-invariant phase-space defined below: - $|y(\pi^+\pi^-)| < 0.4$ - $0.05~\mathrm{GeV}^2 < -t_1, -t_2 < 0.16~\mathrm{GeV}^2$ The experimental systematic uncertainties "(syst.)" are calculated as the quadratic sum of the differences between the nominal fit result and the result of the fit to $d\sigma/dm(\pi^+\pi^-)$ with each systematic effect. The uncertainties related to the extrapolation "(model.)" are quoted as the largest deviation from the nominal fit result.

Exponential slope of the $t$-distribution in three ranges of $m(\pi^+\pi^-)$ and two ranges of $\Delta\varphi$. The results are obtained for the Lorentz-invariant phase-space defined below: - $|y(\pi^+\pi^-)| < 0.4$ - $0.05~\mathrm{GeV}^2 < -t_1, -t_2 < 0.16~\mathrm{GeV}^2$

Differential cross section as a function of the transverse momentum of the pion pair at 13 TeV, compared with generator-level simulations.

Differential cross section as a function of the rapidity of the pion pair at 5.02 TeV, compared with generator-level simulations.

Differential cross section as a function of the rapidity of the pion pair at 13 TeV, compared with generator-level simulations.

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.

Spline interpolation of the weight function (depending of $\cos \theta_{P_c}$) applied to $\Lambda_{b}^{0}$ candidates. The weight function is determined as the inverse of the density of $\Lambda_{b}^{0}$ candidates in the narrow pentaquark peak region while the interpolation is done between the bin centers.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Expected allowed ranges for 8 dimensional aQGC, no cutoff).

Inclusive Cross Section.

Data (bold boxes) and background estimates (colour fill) for m_beta vs. m_betagamma for the gluino R-hadron search (1000 GeV). The blue thin-line boxes illustrate the expected signal (on top of background) for the given R-hadron mass hypothesis. The black dashed vertical/horizontal lines at 500 GeV show the mass selection (signal region in the top-right). Two events pass this selection.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of bottom-squark R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Expected (dashed black line) and observed (solid red line) 95% CL upper limits on the cross section as a function of mass for the production of top-squark R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. The observed 8 TeV Run-1 limit and theory prediction [arXiv:1411.6795] are shown in dash-dotted and dotted lines, respectively.

Final selection requirements as a function of the simulated R-hadron mass.

Summary of all studied systematic uncertainties. Ranges indicate a dependency on the R-hadron mass hypothesis (from low to high masses).

Expected signal yield (Nsig) and efficiency (eff.), estimated background (Nbkg) and observed number of events in data (Nobs) for the full mass range after the final selection using 3.2/fb of data. The stated uncertainties include both the statistical and systematic contribution.

Distribution of the truth-level beta for gluino R-hadrons in exemplary signal MC samples and muons in a Zmumu MC sample. All distributions have been normalised to one. The last bin contains the overflow of the histograms. The distributions illustrate the good discriminating power of the variables.

Distribution of the truth-level betagamma for gluino R-hadrons in exemplary signal MC samples and muons in a Zmumu MC sample. All distributions have been normalised to one. The last bin contains the overflow of the histograms. The distributions illustrate the good discriminating power of the variables.

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 50 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 30 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Expected (dashed black line) and observed (solid red line) 95% confidence level upper limits on the cross section as a function of mass for the production of long-lived gluino R-hadrons. The theory prediction along with its +-1sigma uncertainty is show as a black line and a blue band, respectively. For meta-stable gluinos with a lifetime of 10 ns. (mass exclusion: about 1660 GeV expected, 1520 GeV observed).

Object-quality selection cut-flow with observed data and exemplary expected events (scaled to 3.2/fb for MC) in the gluino R-hadron search.

Object-quality selection cut-flow with observed data and exemplary expected events (scaled to 3.2/fb for MC) in the squark R-hadron search.

Expected signal yield (Nsig) and efficiency (eff.), estimated background (Nbkg) and observed number of events in data (Nobs) for the full mass range in the meta-stable gluino R-hadron search using 3.2/fb of data. The stated uncertainties include both the statistical and systematic contribution.

The distribution of the missing transverse momentum is shown in soft-lepton 2-jet W+jets control regions after normalising the ttbar and W+jets background processes in the simultaneous fit.

Expected background yields as obtained in the background-only fits in all hard-lepton and soft-lepton validation together with observed data are given. Uncertainties in the fitted background estimates combine statistical (in the simulated event yields) and systematic uncertainties.

Expected background yields as obtained in the background-only fits in all hard-lepton and soft-lepton signal together with observed data are given. Uncertainties in the fitted background estimates combine statistical (in the simulated event yields) and systematic uncertainties.

Distributions of mt for the hard-lepton 4-jet low-x signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of met/meff for the 4-jet high-x signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of mt for the hard-lepton 5-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of mt for the hard-lepton 6-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of met for the soft-lepton 2-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

Distributions of met for the soft-lepton 5-jet signal region. The requirement on the variable plotted is removed from the definitions of the signal regions, where the arrow indicates the position of the cut in the signal region. The lower panels of the plots show the ratio of the observed data to the total background prediction as derived in the background-only fit. The uncertainty bands plotted include all statistical and systematic uncertainties as discussed in Section 7. The component `Others' is the sum of Z+jets and ttbar+V. The last bin includes the overflow.

The observed combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The yellow band ($+ 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The yellow band ($- 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The observed combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The yellow band ($+ 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The yellow band ($- 1 \sigma$) of the combined 95% CL exclusion limits in the the gluino simplified models using for each model point the signal region with the best expected sensitivity. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$. The yellow band represents the $\pm 1 \sigma$ variation of the median expected limit due to the experimental and theoretical uncertainties.

The observed limits for the soft-lepton 2-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the soft-lepton 2-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The expected limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ models.

The observed limits for the soft-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the soft-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 4-jet low-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 4-jet high-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 4-jet high-x signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 5-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The observed limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

The expected limits for the hard-lepton 6-jet signal region. The limits are presented in the (gluino, x) plane for the chargino = 60 GeV models where $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

Number of generated events in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

Number of generated events in the (gluino, x) plane for the chargino = 60 GeV models.

Production cross-section in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

Production cross-section in the (gluino, x) plane for the chargino = 60 GeV models.

Acceptance times efficiency obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet low-x).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet high-x).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 5-jet).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 6-jet).

Acceptance times efficiency in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (soft-lepton 2-jet).

Acceptance times efficiency obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet low-x).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet high-x).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 5-jet).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 6-jet).

Acceptance times efficiency in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (soft-lepton 5-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet low-x).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet high-x).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 5-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 6-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (soft-lepton 2-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet low-x).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 4-jet high-x).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 5-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (hard-lepton 6-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$ (soft-lepton 2-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet low-x).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet high-x).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 5-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 6-jet).

The observed CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (soft-lepton 5-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet low-x).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 4-jet high-x).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 5-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (hard-lepton 6-jet).

The expected CLs values as obtained in the different signal regions in the (gluino, x) plane for the chargino = 60 GeV models (soft-lepton 5-jet).

The signal region yielding in the best expected limit is indicated for every signal point used in the the gluino simplified models for the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

The signal region yielding in the best expected limit is indicated for every signal point used in the the gluino simplified models for the (gluino, x) mass plane where for the chargino = 60 GeV and $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

Model-dependent 95% CL upper limits on the visible cross-section in addition to the observed and expected exclusion limits for the (gluino, chargino) mass plane for the scenario where the mass of the chargino is fixed to $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1)) = 1/2$.

Model-dependent 95% CL upper limits on the visible cross-section in addition to the observed and expected exclusion limits for the (gluino, x) mass plane where for the chargino = 60 GeV and $x=(m(\tilde\chi^\pm_1)-m(\tilde\chi^0_1))/(m(\tilde g) - m(\tilde\chi^0_1))$.

Simulated background event samples: the corresponding generator, parton shower, cross-section normalisation, PDF set and underlying-event tune are shown.

Overview of the selection criteria for the soft-lepton signal regions. The symbol $p_{T}^{l}$ refers to signal leptons.

Overview of the selection criteria for the hard-lepton signal regions. The symbol $p_{T}^{l}$ refers to signal leptons.

Background fit results for the hard-lepton and soft-lepton signal regions, for an integrated luminosity of 3.2 fb-1. Uncertainties in the fitted background estimates combine statistical (in the simulated event yields) and systematic uncertainties. The uncertainties in this table are symmetrised for propagation purposes but truncated at zero to remain within the physical boundaries.

Breakdown of upper limits. The columns show from left to right: the name of the respective signal region; the 95% confidence level (CL) upper limits on the visible cross-section and on the number of signal events the 95% CL upper limit on the number of signal events, given the expected number (and $\pm 1 \sigma$ variations on the expectation) of background events; the two-sided CLb value, i.e. the confidence level observed for the background-only hypothesis and the one-sided discovery p-value (p(s = 0)). The discovery p-values are capped to 0.5 in the case of observing less events than the fitted background estimates.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the mt distribution shown in figure 5 (top left). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the met/meff distribution shown in figure 5 (top right). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the mt distribution shown in figure 5 (middle left). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the mt distribution shown in figure 5 (middle right). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the met distribution shown in figure 5 (bottom left). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Table shows the data, fitted background and expected signal event counts for a benchmark signal point in each bin of the met distribution shown in figure 5 (bottom right). The fit results are shown for an integrated luminosity of 3.2 fb-1.

Cutflow table for the hard-lepton signal regions with representative target signal models. The weighted numbers are normalized to 3.2 fb-1 and rounded to the statistical error.

Cutflow table for the hard-lepton signal regions with representative target signal models. The weighted numbers are normalized to 3.2 fb-1 and rounded to the statistical error.

The leptonic asymmetry dependence on the invariant mass of top-antitop system ($m_{t\bar{t}}$) in the fiducial volume.

The leptonic asymmetry dependence on the longitudinal boost of the top-antitop system ($\beta_z^{t\bar{t}}$) in the fiducial volume.

The leptonic asymmetry dependence on the transverse momentum of the top-antitop system ($p_T^{t\bar{t}}$) in the fiducial volume.

The top-antitop inclusive asymmetry in the fiducial volume.

The top-antitop asymmetry dependence on the invariant mass of top-antitop system ($m_{t\bar{t}}$) in the fiducial volume.

The top-antitop asymmetry dependence on the longitudinal boost of the top-antitop system ($\beta_z^{t\bar{t}}$) in the fiducial volume.

The top-antitop asymmetry dependence on the transverse momentum of the top-antitop system ($p_T^{t\bar{t}}$) in the fiducial volume.

Excluded range of lifetimes as a function of gluino R-hadron mass. The central value of the lower-limit of the observed and expected excluded masses at the 95% CL is given.

Cross-section as a function of mass for stable gluino R-hadrons. The central value of the observed and expected 95% upper limit on excluded cross-section is given. The 1 sigma band on the expected UL is given as an uncertainty.

Measured fiducial cross section in fb as a function of Yll.

Measured fiducial cross section in fb/GeV as a function of the leading jet pT, PT(JET1).

Measured normalised fiducial cross section as a function of Njet. Jet PT>25 GeV for |eta|<2.4 and PT>30 GeV for 2.4<|eta|<4.5.

Measured normalised fiducial cross section in 1/GeV as a function of pTH.

Measured normalised fiducial cross section as a function of Yll.

Measured normalised fiducial cross section in 1/GeV as a function of the leading jet pT, PT(JET1).

Jet veto efficiency as a function of PT(JET1) pT threshold.

Full correlation matrix for the measured differential cross section as a function of Njet, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured differential cross section as a function of pTH, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured differential cross section as a function of Yll, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured differential cross section as a function of the leading jet pT, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured normalised differential cross section as a function of Njet, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured normalised differential cross section as a function of pTH, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured normalised differential cross section as a function of Yll, taking into account statistical and systematic uncertainties.

Full correlation matrix for the measured normalised differential cross section as a function of the leading jet pT, taking into account statistical and systematic uncertainties.

Fiducial cross section of H+0-jet events in femtobarn for different values of PT(JET1) pT threshold with the uncertainties for each bin given in absolute values and in percent. The asterisk for the 25 GeV column header indicates that the results are for a mixed jet pT threshold, which is raised from 25 GeV to 30 GeV for jets with 2.5 < |eta| < 4.5.

Correction factors from inclusive parton level to fiducial particle level for bins of Njet derived with POWHEG NNLOPS+Pythia8. Uncertainties due to the non-perturbative correction are dominant for the correction factors.

Correction factors from inclusive parton level to fiducial particle level for bins of the leading jet pT derived with POWHEG NNLOPS+Pythia8. Uncertainties due to the non-perturbative correction are dominant for the correction factors.

Correction factors from inclusive parton level to fiducial particle level for the jet-veto efficiency with different jet pT thresholds derived with POWHEG NNLOPS+Pythia8. The asterisk on the 25 GeV bin label indicates that the results are for a mixed pT threshold, which is raised from 25 GeV to 30 GeV for jets with 2.5 < |eta| < 4.5, corresponding to the selection used to define the signal regions for the analysis.

The measured fiducial cross sections are corrected to the dressed level. The ratio of $C_{WZ}^{\textrm{dressed}}/C_{WZ}^{\textrm{Born}}$ factors presented in this table can be used to correct fiducial integrated cross sections from dressed to Born level.

Cross section extrapolated to total phase space and all W and Z boson decays. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the theory uncertainty and the third is the luminosity uncertainty.

The total cross section is measured at dressed level and can also be corrected to Born level use this acceptance correction factor.

Ratio of fiducial cross sections measured for $W^{+} Z$ and $W^{-} Z$ production. The first systematic uncertainty is the combined systematic uncertainty excluding the luminosity uncertainty, the second is the luminosity uncertainty.

Observed and expected upper limits at $95\%$ CL in fb on the fiducial cross section $\sigma^{\mathrm{fid.}}_{W^{\pm} Zjj \textrm{-EW} \rightarrow \ell^{'} \nu \ell \ell}$, multiplied by the $W$ and $Z$ branching ratios in a single leptonic channel with muons or electrons in the VBS fiducial phase space. Values obtained with or without subtraction of the $tZj$ contribution are presented. The $1$ $\sigma$ and $2$ $\sigma$ uncertainty intervals around the expected limits are also indicated.

Observed and expected upper limits at $95\%$ CL in fb on the fiducial cross section $\sigma^{\mathrm{fid.}}_{W^{\pm} Zjj \textrm{-EW} \rightarrow \ell^{'} \nu \ell \ell}$, multiplied by the $W$ and $Z$ branching ratios in a single leptonic channel with muons or electrons in the aQGC fiducial phase space. Values obtained with or without subtraction of the $tZj$ contribution are presented. The $1$ $\sigma$ and $2$ $\sigma$ uncertainty intervals around the expected limits are also indicated.

Expected and observed one-dimensional 95% Confidence Level intervals on the anomalous couplings parameters. Anomalous couplings are introduced as form factors using a cut-off scale $\Lambda_{\mathrm{co}}$.

Expected and observed one-dimensional intervals at 95% Confidence Level on the Effective Field Theory (EFT) parameters.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Measured normalised fiducial cross section in all $\ell'^\pm \nu \ell^+ \ell^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding theory and luminosity uncertainties, the second is the luminosity uncertainty. The last bin is a cross section for all events above the lower end of the bin.

Ratio of $W^{+}Z$ over $W^{-}Z$ differential fiducial cross sections measured in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The last bin is a cross section for all events above the lower end of the bin.

Ratio of $W^{+}Z$ over $W^{-}Z$ differential fiducial cross sections measured in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The last bin is a cross section for all events above the lower end of the bin.

Ratio of $W^{+}Z$ over $W^{-}Z$ differential fiducial cross sections measured in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The last bin is a cross section for all events above the lower end of the bin.

Ratio of $W^{+}Z$ over $W^{-}Z$ differential fiducial cross sections measured in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The last bin is a cross section for all events above the lower end of the bin.

Ratio of $W^{+}Z$ over $W^{-}Z$ differential fiducial cross sections measured in all $\ell'^\pm \nu \ell^+ \ell'^-$ channels combined, where $\ell, \ell' = e, \mu$.

Extrapolation factors AWW and CWW for emu channel.

Extrapolation factors AWW and CWW for ee channel.

Extrapolation factors AWW and CWW for mumu channel.

Uncertainty on extrapolation factors AWW and CWW for emu channel.

Uncertainty on extrapolation factors AWW and CWW for ee channel.

Uncertainty on extrapolation factors AWW and CWW for mumu channel.

Measured differential cross sections as function of the leading lepton transverse momentum of WW production in the fiducial region for same different final states corresponding to different W decay channels: both W's decaying into electrons, both decaying to muon, and final states with one prompt electron and one prompt muon. The cross sections are defined for direct decays of the W bosons into prompt electrons or muons, intermediate decays into tau leptons are disregarded. The electrons are required to be contained within abs(eta)<2.47 and to lie outside of 1.37 < abs(eta) < 1.53, muons are required to lie within abs(eta)<2.4. The leading and subleading leptons in the events are required to have a transverse momentum above 25 and 20 GeV respectively. The transverse momentum of the vectorial sum of the neutrinos in the event should be larger than 20 GeV (PT(C=SUM(NU))). The transverse momentum of the vectorial sum of the neutrinos multiplied by the sin of azimuthal difference between lepton and the vectorial sum of the neutrinos in the event should be larger than 15 GeV. The invariant mass of the leptons should exceed 10 GeV. Particle-level jets are defined using the anti-kT algorithm with radius of 0.4. No jets above 25 GeV and within abs(eta)<4.5 are allowed in the event. Both, resonant and non-resonant WW production processes, are included in the cross sections. All differential distributions are obtained using this selections.

All differential distributions are obtained using the selection described for the transverse momentum distribution of leading lepton (PT(C=LEADING LEPTON)).

All differential distributions are obtained using the selection described for the transverse momentum distribution of leading lepton (PT(C=LEADING LEPTON)).

All differential distributions are obtained using the selection described for the transverse momentum distribution of leading lepton (PT(C=LEADING LEPTON)).

All differential distributions are obtained using the selection described for the transverse momentum distribution of leading lepton (PT(C=LEADING LEPTON)).

All differential distributions are obtained using the selection described for the transverse momentum distribution of leading lepton (PT(C=LEADING LEPTON)). Costheta is defined as: abs( tanh(DELTA-ETA(C=DILEPTON)/2 )).

Bin-to-bin correlation matrix for the differential cross section measurements for the pt(leading lepton) distribution.

Bin-to-bin correlation matrix for the normalised differential distribution for the pt(leading lepton) distribution.

Bin-to-bin correlation matrix for the differential cross section measurements for the transverse momentum distribution of the dilepton system.

Bin-to-bin correlation matrix for the normalised differential distribution for the transverse momentum distribution of the dilepton system.

Bin-to-bin correlation matrix for the differential cross section measurements for the invariant mass of the dilepton system.

Bin-to-bin correlation matrix for the normalised differential distribution for the invariant mass of the dilepton system.

Bin-to-bin correlation matrix for the differential cross section measurements for the difference in azimuthal angle between the two leptons.

Bin-to-bin correlation matrix for the normalised differential distribution for the difference in azimuthal angle between the two leptons.

Bin-to-bin correlation matrix for the differential cross section measurements for the absolute rapidity of the dilepton system.

Bin-to-bin correlation matrix for the normalised differential distribution for the absolute rapidity of the dilepton system.

Bin-to-bin correlation matrix for the differential cross section measurements for costheta distribution of the two leptons.

Bin-to-bin correlation matrix for the normalised differential distribution for the costheta distribution of the two leptons.

The expected and observed 95% confidence intervals for the anomalous coupling parameters defined in the no constraints scenario, LEP, HISZ and Equal Couplings scenarios. The results are shown with Lambda = infinity and Lamdba = 7 TeV.

The expected and observed 95% confidence intervals for the anomalous coupling parameters defined in the EFT frame work.

The average charged particle multiplicity for the more central jet and a charged particle threshold of 0.5 GeV as a function of the jet transverse momentum.

The average charged particle multiplicity for the more central jet and a charged particle threshold of 2 GeV as a function of the jet transverse momentum.

The average charged particle multiplicity for the more central jet and a charged particle threshold of 5 GeV as a function of the jet transverse momentum.

The average charged particle multiplicity differnce between the more forward and the more central jets and a charged particle threshold of 0.5 GeV as a function of the jet transverse momentum.

The average charged particle multiplicity differnce between the more forward and the more central jets and a charged particle threshold of 2 GeV as a function of the jet transverse momentum.

The average charged particle multiplicity differnce between the more forward and the more central jets and a charged particle threshold of 5 GeV as a function of the jet transverse momentum.

The average charged particle multiplicity for the sum of the more forward and the more central jets and a charged particle threshold of 0.5 GeV as a function of the jet transverse momentum.

The average charged particle multiplicity for the sum of the more forward and the more central jets and a charged particle threshold of 2 GeV as a function of the jet transverse momentum.

The average charged particle multiplicity for the sum of the more forward and the more central jets and a charged particle threshold of 5 GeV as a function of the jet transverse momentum.

three particle correlator, 3rd harmonic, difference neg-pos, as a function of centrality.

three particle correlator, 4th harmonic, as a function of centrality.

three particle correlator, 4th harmonic, difference neg-pos, as a function of centrality.

three particle correlator, 2nd harmonic, inferred slope parameter, as a function of centrality.

charge correlation, as a function of Delta eta.

three particle correlator, 2nd harmonic, as a function of Delta eta.

three particle correlator, 3rd harmonic, as a function of Delta eta.

three particle correlator, 4th harmonic, as a function of Delta eta.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} p_T(D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(\Upsilon(1S))}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(\Upsilon(1S))}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(D^0)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{\pi}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \left| \Delta \phi \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{\pi}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \left| \Delta \phi \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \left| \Delta y \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \left| \Delta y \right| }$ for $2<y(\Upsilon(1S))<4.5$, for $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} p_T(\Upsilon(1S)D^0) }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} p_T(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(\Upsilon(1S)D^0) }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \mathcal{A}_T }$, where $\mathcal{A}_T = \frac{p_T(\Upsilon(1S))-p_T(D^0)}{p_T(\Upsilon(1S))-p_T(D^0)}$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \mathcal{A}_T }$, where $\mathcal{A}_T = \frac{p_T(\Upsilon(1S))-p_T(D^0)}{p_T(\Upsilon(1S))-p_T(D^0)}$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} m(\Upsilon(1S)D^0) }$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} m(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Expected and observed measurements of the cross section and signal strength with 68% CL ranges and significances for t$\mathrm{\bar{t}}$W, in SS dilepton and 3$\ell$ channels.

Expected and observed measurements of the cross section and signal strength with 68% CL ranges and significances for t$\mathrm{\bar{t}}$Z, in OS dilepton, 3$\ell$, and 4$\ell$ channels.

Constraints from this t$\mathrm{\bar{t}}$Z and t$\mathrm{\bar{t}}$W measurement on selected dimension-six operators.

Measured differential four-jet cross section for R=0.4 jets, in bins of pT4, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of HT, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of m_4j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of min(m_2j)/m_4j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m_2j>500 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of min(m_2j)/m_4j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m_2j>1000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of min(m_2j)/m_4j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m_2j>1500 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of min(m_2j)/m_4j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as m_2j>2000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>700 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>1000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>700 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDphi_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>1000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>700 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>1000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>700 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of minDy_3j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>1000 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of maxDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of maxDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>250 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of maxDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of maxDy_2j, along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts, as well as pT1>550 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>1. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>1, as well as pT1>250 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>1, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>1, as well as pT1>550 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>2. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>2, as well as pT1>250 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>2, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>2, as well as pT1>550 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>3. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>3, as well as pT1>250 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>3, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>3, as well as pT1>550 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>4. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>4, as well as pT1>250 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>4, as well as pT1>400 GeV. All other details are as for pT1.

Measured differential four-jet cross section for R=0.4 jets, in bins of sum(pT), along with the uncertainties in the measurement. The events are selected using the inclusive analysis cuts and maxDy_2j>4, as well as pT1>550 GeV. All other details are as for pT1.

$p_{\rm T}$-differential yield of prompt ${\rm D}^{0}$ mesons in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 30-50% in the rapidity interval |y|<0.5. Branching ratio of ${\rm D}^{0}$->${\rm K}^{0}\pi^{+}$ : 0.0388. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

$p_{\rm T}$-differential yield of prompt ${\rm D}^{+}$ mesons in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 30-50% in the rapidity interval |y|<0.5. Branching ratio of ${\rm D}^{+}$->${\rm K}^{-}\pi^{+}\pi^{+}$ : 0.0913. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

$p_{\rm T}$-differential yield of prompt ${\rm D}^{*+}$ mesons in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 30-50% in the rapidity interval |y|<0.5. Branching ratio of ${\rm D}^{*+}$->${\rm D}^{0}\pi^{+}$->${\rm K}^{-}\pi^{+}\pi^{+}$ : 0.0388*0.677. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

Ratio of ${\rm D}^{+}$ to ${\rm D}^{0}$ meson yield in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 0-10% in |y| < 0.5 as a function of $p_{\rm T}$.

Ratio of ${\rm D}^{*+}$ to ${\rm D}^{0}$ meson yield in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 0-10% in |y| < 0.5 as a function of $p_{\rm T}$.

Nuclear modification factor $R_{\rm AA}$ of ${\rm D}^{0}$ mesons in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 0-10% in |y| < 0.5 as a function of $p_{\rm T}$.

Nuclear modification factor $R_{\rm AA}$ of ${\rm D}^{+}$ mesons in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 0-10% in |y| < 0.5 as a function of $p_{\rm T}$.

Nuclear modification factor $R_{\rm AA}$ of ${\rm D}^{*+}$ mesons in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 0-10% in |y| < 0.5 as a function of $p_{\rm T}$.

Nuclear modification factor $R_{\rm AA}$ of ${\rm D}^{0}$ mesons in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 30-50% in |y| < 0.5 as a function of $p_{\rm T}$.

Nuclear modification factor $R_{\rm AA}$ of ${\rm D}^{+}$ mesons in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 30-50% in |y| < 0.5 as a function of $p_{\rm T}$.

Nuclear modification factor $R_{\rm AA}$ of ${\rm D}^{*+}$ mesons in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 30-50% in |y| < 0.5 as a function of $p_{\rm T}$. Branching ratio of ${\rm D}^{*+}$->${\rm D}^{0}\pi^{+}$->${\rm K}^{-}\pi^{+}\pi^{+}$ : 0.0388*0.677.

Nuclear modification factor $R_{\rm AA}$ of Average ${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ mesons in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 0-10% in |y| < 0.5 as a function of $p_{\rm T}$.

Nuclear modification factor $R_{\rm AA}$ of Average ${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ mesons in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the centrality class 30-50% in |y| < 0.5 as a function of $p_{\rm T}$.

Ratio of prompt D meson (${\rm D}^{0}$, ${\rm D}^{+}$, ${\rm D}^{*+}$ average) and charged pion (charged particles for $p_{\rm T}>16~{\rm GeV/}c$) nuclear modification factors ($R_{\rm AA}$) as a function of $p_{\rm T}$ for D mesons with $|y|<0.5$ and charged pions (particles) with $|\eta|<0.8$, measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76~{\rm TeV}$ in the 0-10% centrality range.

The jet charge distribution standard deviation for kappa = 0.3 and the more central jet. Values are given in units of the positron charge, e.

The jet charge mean for kappa = 0.5 and the more forward jet. Values are given in units of the positron charge, e.

The jet charge distribution standard deviation for kappa = 0.5 and the more forward jet. Values are given in units of the positron charge, e.

The jet charge mean for kappa = 0.5 and the more central jet. Values are given in units of the positron charge, e.

The jet charge distribution standard deviation for kappa = 0.5 and the more central jet. Values are given in units of the positron charge, e.

The jet charge mean for kappa = 0.7 and the more forward jet. Values are given in units of the positron charge, e.

The jet charge distribution standard deviation for kappa = 0.7 and the more forward jet. Values are given in units of the positron charge, e.

The jet charge mean for kappa = 0.7 and the more central jet. Values are given in units of the positron charge, e.

The jet charge distribution standard deviation for kappa = 0.7 and the more central jet. Values are given in units of the positron charge, e.

Double differential cross-section for $J/\psi$-from-$b$ mesons as a function of $p_\perp$ in bins of $y$. The first uncertainties are statistical, the second are the correlated systematic uncertainties shared between bins and the last are the uncorrelated systematic uncertainties.

The fraction of $J/\psi$-from-$b$ mesons (in %) in bins of the $J/\psi$ $p_\perp$ and $y$. The uncertainties are statistical only. The systematic uncertainties are negligible.

The fraction of $J/\psi$-from-$b$ mesons (in %) in bins of the $J/\psi$ $p_\perp$ and $y$. The uncertainties are statistical only. The systematic uncertainties are negligible.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Acceptance times efficiency for the tight channel of the Razor analysis for the direct squark-squark model.

Acceptance for the 1tau channel of the Taus+jets+MET analysis for the gluino simplified model.

Efficiency for the 1tau channel of the Taus+jets+MET analysis for the gluino simplified model.

Acceptance for the 2tau channel of the Taus+jets+MET analysis for the gluino simplified model.

Efficiency for the 2tau channel of the Taus+jets+MET analysis for the gluino simplified model.

Acceptance for the 1tau channel of the Taus+jets+MET analysis for the squark simplified model.

Efficiency for the 1tau channel of the Taus+jets+MET analysis for the squark simplified model.

Acceptance for the 2tau channel of the Taus+jets+MET analysis for the squark simplified model.

Efficiency for the 2tau channel of the Taus+jets+MET analysis for the squark simplified model.

Observed limit on the cross-section times branching ratio in fb at 95% CL for the direct squark-squark model.

Observed limit on the cross-section times branching ratio in fb at 95% CL for the gluino simplified model, obtained by the Taus+jets+MET analysis.

Observed limit on the cross-section times branching ratio in fb at 95% CL for the squark simplified model, obtained by the Taus+jets+MET analysis.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M1=60 GeV and $m_{\tilde g}=$1.6 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M2=(M1+mqL)/2 and $m_{\tilde g}=$1.6 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M1=60 GeV with $m_{\tilde g}=$2.2 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M2=(M1+mqL)/2 and $m_{\tilde g}=$2.2 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M1=60 GeV and $m_{\tilde g}=$3.0 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed 95% CL cross section times branching ratio upper limit on pMSSM qLqL production with M2=(M1+mqL)/2 and $m_{\tilde g}=$3.0 TeV interpreted by the 0L analysis including its two dedicated signal regions.

Observed and expected limits on the cross-section times branching ratio in fb at 95% CL for the gluino-mediated stop model with off-shell stops, for mass points with four- or five-body gluino decays, obtained by the SS/3L and 0/1L3B analyses.

: Cut flow for signal regions 0L_4jt+ and 0L_5jt of baseline signal points normalized to 20.3 fb$^\mathrm{-1}$. The point used for the 0L_4jt+ cutflow is mqL = 1150 GeV, M1 = 60 GeV, M2 = 1000 GeV, m(gluino) = 2200 GeV. The point used for the 0L_5jt cutflow is mqL = 950 GeV, M1 = 100 GeV, M2 = 525 GeV, m(gluino) = 2200 GeV.

Cut flow for the 1L(H)_7-jet SR (electron channel) for the main background from the $t\bar{t}$ production. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for the 1L(H)_7-jet SR (muon channel) for the main background from the $t\bar{t}$ production. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for 1L(H)_7-jet SR (electron channel). The signal point is taken from the gluino simplified model with 1-step decay via charginos and parameters $m_{\tilde{g}}$ = 1145 GeV, $m_{{\tilde \chi}_{1}^{\pm}}$ = 785 GeV, $m_{\tilde{\chi}_{1}^{0}}$ = 425 GeV; 20000 events were generated for this point. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for 1L(H)_7-jet SR (muon channel). The signal point is taken from the gluino simplified model with 1-step decay via charginos and parameters $m_{\tilde{g}}$ = 1145 GeV, $m_{\tilde{\chi}_{1}^{\pm}}$ = 785 GeV, $m_{\tilde{\chi}_{1}^{0}}$ = 425 GeV; 20000 events were generated for this point. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for 1L(H)_7-jet SR (electron channel). The signal point is taken from the gluino simplified model with 1-step decay via charginos and parameters $m_{\tilde{g}}$ = 1025 GeV, $m_{\tilde{\chi}_{1}^{\pm}}$ = 545 GeV, $m_{\tilde{\chi}_{1}^{0}}$ = 65 GeV; 20000 events were generated for this point. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for 1L(H)_7-jet SR (muon channel). The signal point is taken from the gluino simplified model with 1-step decay via charginos and parameters $m_{\tilde{g}}$ = 1025 GeV, $m_{\tilde{\chi}_{1}^{\pm}}$ = 545 GeV, $m_{\tilde{\chi}_{1}^{0}}$ = 65 GeV; 20000 events were generated for this point. The events are normalized to 20.3 fb$^\mathrm{-1}$.

Cut flow for 0LRaz_SR$_{\rm loose}$ and 0LRaz_SR$_{\rm tight}$ signal regions of baseline signal point ($m(\tilde{q}, \tilde{\chi}_1^0)=(850, 100)$ GeV) normalized to 20.3 fb$^{-1}$. 5000 events were generated for ths point.

Detector response matrix (PROB) for the acoplanarity variable (ACO) for mu+ mu- channel (empty bins are not reported).

Acoplanarity (ACO) distributions unfolded for detector resolution, and lepton pair trigger, reconstruction and identification efficiencies for e+ e- channel (empty bins are not reported).

Acoplanarity (ACO) distributions unfolded for detector resolution, and lepton pair trigger, reconstruction and identification efficiencies for mu+ mu- channel (empty bins are not reported).

The lepton pair exclusive selection efficiency in the fiducial region for signal MC events as a function of M and Y for e+ e- channel.

The lepton pair exclusive selection efficiency in the fiducial region for signal MC events as a function of M and Y for mu+ mu- channel.

The lepton pair exclusive selection efficiency in the fiducial region for single-dissociative MC events as a function of M and Y for e+ e- channel.

The lepton pair exclusive selection efficiency in the fiducial region for single-dissociative MC events as a function of M and Y for mu+ mu- channel.

The observed 95\% C.L. limit for $pp\rightarrow G^{*}_{KK}\rightarrow hh\rightarrow b\bar{b}b\bar{b}$ in the bulk RS model with $k/\bar{M}_{Pl} = 2$, as a function of resonance mass.

The observed 95\% C.L. limit for $pp\rightarrow H\rightarrow hh\rightarrow b\bar{b}b\bar{b}$ with fixed $\Gamma_{H} = 1$ GeV, as a function of resonance mass.

The $\gamma$ differential transverse momentum cross-section in an inclusive $\gamma+\mathrm{jets}$, $N_{\mathrm{jets}} \geq2$ selection for central rapidities $\vert y_{\gamma} \vert > 1.4$.

The Z boson differential transverse momentum cross-section in an inclusive $Z/\gamma^{*}+\mathrm{jets}$, $N_{\mathrm{jets}} \geq2$ over $Z/\gamma^{*}+\mathrm{jets}$, $N_{\mathrm{jets}} \geq1$ selection.

The $\gamma$ boson differential transverse momentum cross-section in an inclusive $\gamma+\mathrm{jets}$, $N_{\mathrm{jets}} \geq2$ over $Z+\mathrm{jets}$, $N_{\mathrm{jets}} \geq1$ selection for central rapidities $\vert y_{\gamma} \vert < 1.4$.

The Z boson differential transverse momentum over $H_{\mathrm{T}}$ cross-section in an inclusive $Z/\gamma^{*}+\mathrm{jets}$, $N_{\mathrm{jets}} \geq2$ selection.

The Z boson differential transverse momentum over $H_{\mathrm{T}}$ cross-section in an inclusive $Z/\gamma^{*}+\mathrm{jets}$, $N_{\mathrm{jets}} \geq3$ selection.

The log (base 10) of the Z boson differential transverse momentum over the leading jet transverse momentum cross-section in an inclusive $Z/\gamma^{*}+\mathrm{jets}$, $N_{\mathrm{jets}} \geq2$ selection.

The log (base 10) of the Z boson differential transverse momentum over the leading jet transverse momentum cross-section in an inclusive $Z/\gamma^{*}+\mathrm{jets}$, $N_{\mathrm{jets}} \geq3$ selection.

Differential cross-section ratio of Z over $\gamma$ as a function of the total transverse momentum cross-section and for central bosons ($\vert y_{\mathrm{V}}\vert < 1.4$) in an inclusive $ N_{\mathrm{jets}}\geq 1 $ selection.

Differential cross-section ratio of Z over $\gamma$ as a function of the total transverse momentum cross-section and for central bosons ($\vert y_{\mathrm{V}}\vert < 1.4$) in a $ H_{\mathrm{T}}$% > 300 GeV, N_{\mathrm{jets}}\geq 1 $ selection.

Differential cross-section ratio of Z over $\gamma$ as a function of the total transverse momentum cross-section and for central bosons ($\vert y_{\mathrm{V}}\vert < 1.4$) in an inclusive $ N_{\mathrm{jets}}\geq 2 $ selection.

Differential cross-section ratio of Z over $\gamma$ as a function of the total transverse momentum cross-section and for central bosons ($\vert y_{\mathrm{V}}\vert < 1.4$) in an inclusive $ N_{\mathrm{jets}}\geq 3 $ selection.

Vertex-level efficiency if no re-tracking is performed, as a function of the vertex radial position for an RPV SUSY model of squark production with $\tilde{q}\to q[\tilde{\chi}_1^0\to\mu qq]$, where $m(\tilde{q}) = 700$ GeV, $m(\tilde{\chi}_1^0) = 494$ GeV and $c\tau(\tilde{\chi}_1^0)$ = 175 mm. The result with re-tracking is tabulated in http://hepdata.cedar.ac.uk/view/ins1362183/d1.

Upper limits (95% CL) on the number of DV+muon displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 7a) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark or gluino pair production, with masses in GeV as indicated, with either $\tilde{q}\to q[\tilde{\chi}_1^0\to\mu qq]$ or $\tilde{g}\to qq[\tilde{\chi}_1^0\to\mu qq]$ decays, respectively.

Upper limits (95% CL) on the number of DV+muon displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 7b) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark or gluino pair production, with masses in GeV as indicated, with either $\tilde{q}\to q\tilde{\chi}_1^0$ or $\tilde{g}\to qq\tilde{\chi}_1^0$ decays, respectively.

Upper limits (95% CL) on the number of DV+electron displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 7c) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark or gluino pair production, with masses in GeV as indicated, with either $\tilde{q}\to q[\tilde{\chi}_1^0\to e qq]$ or $\tilde{g}\to qq[\tilde{\chi}_1^0\to e qq]$ decays, respectively.

Upper limits (95% CL) on the number of DV+electron displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 7d) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider squark or gluino pair production, with masses in GeV as indicated, with either $\tilde{q}\to q\tilde{\chi}_1^0$ or $\tilde{g}\to qq\tilde{\chi}_1^0$ decays, respectively.

Upper limits (95% CL) on the number of dilepton $e^+e^-$ displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 5a) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to ee\nu]$ decays and masses in GeV as indicated. A dash indicates results omitted due to limited statistical precision.

Upper limits (95% CL) on the number of dilepton $e^{\pm}\mu^{\mp}$ displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 5b) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to e\mu\nu]$ decays and masses in GeV as indicated.

Upper limits (95% CL) on the number of $\mu^+\mu^-$ displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 5c) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to \mu\mu\nu]$ decays and masses in GeV as indicated. A dash indicates results omitted due to limited statistical precision.

Upper limits (95% CL) on the number of dilepton ($e^+e^-+e^{\pm}\mu^{\mp}+\mu^+\mu^-$) displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 5d) for two GGM SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to Z\tilde{G}]$ decays, $m(\tilde{g})$ = 1100 GeV and $m(\tilde{\chi}_1^0)$ = 400 GeV.

Upper limits (95% CL) on the number of dilepton ($e^+e^-+e^{\pm}\mu^{\mp}+\mu^+\mu^-$) displaced vertices in 20.3 fb$^{-1}$ and the corresponding vertex-level efficiencies (from Auxiliary Figure 5d) for two GGM SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to Z\tilde{G}]$ decays, $m(\tilde{g})$ = 1100 GeV and $m(\tilde{\chi}_1^0)$ = 1000 GeV.

Upper limits (95% CL) from the dilepton $e^+e^-+e^{\pm}\mu^{\mp}$ channels on the production cross-section and the corresponding event-level efficiencies (from Auxiliary Figure 6a) for four RPV SUSY models as a function of the $\tilde{\chi}_1^0$ proper decay distance $c\tau$. The models consider gluino pair production, with $\tilde{g}\to qq[\tilde{\chi}_1^0\to e\mu\nu/ee\nu]$ decays and masses in GeV as indicated. For comparison, the production cross-sections for $m(\tilde{g})$ = 600 GeV and 1300 GeV are $1280\pm230$ fb and $1.7\pm0.7$ fb, respectively.