$\gamma p \rightarrow J/\psi p$ differential cross sections 10.36–11.44 GeV beam energy range, average $t$ and beam energy in bins of $t$. The first cross section uncertainties are statistical, and the second are systematic. The overall average beam energy is 10.82 GeV. There is an additional fully correlated systematic uncertainty of 19.5% on the total cross section, not included here.

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.

This table contains thinned out samples of the Markov chains used in the parameter estimation of the SDME measurements for $-(t-t_\text{0}) = 0.400\pm0.056~\text{GeV}^2/c^2$, reported in the main article. One in about 250 steps in the chain, which results in 200 different sets of SDMEs, is provided. These values should be used instead of bootstrapping of the results, in order to estimate uncertainties of physics models fitted to this data. To assess how the uncertainties propagate to the model uncertainties, one should evaluate the model under scrutiny for each of the 200 different sets of SDMEs. Plotting all resulting lines in a single plot will create bands which reflect the influence of the uncertainties in the data on the model. This method has the great advantage that all correlations are accurately taken into account.

This table contains thinned out samples of the Markov chains used in the parameter estimation of the SDME measurements for $-(t-t_\text{0}) = 0.597\pm0.057~\text{GeV}^2/c^2$, reported in the main article. One in about 250 steps in the chain, which results in 200 different sets of SDMEs, is provided. These values should be used instead of bootstrapping of the results, in order to estimate uncertainties of physics models fitted to this data. To assess how the uncertainties propagate to the model uncertainties, one should evaluate the model under scrutiny for each of the 200 different sets of SDMEs. Plotting all resulting lines in a single plot will create bands which reflect the influence of the uncertainties in the data on the model. This method has the great advantage that all correlations are accurately taken into account.

This table contains thinned out samples of the Markov chains used in the parameter estimation of the SDME measurements for $-(t-t_\text{0}) = 0.793\pm0.057~\text{GeV}^2/c^2$, reported in the main article. One in about 250 steps in the chain, which results in 200 different sets of SDMEs, is provided. These values should be used instead of bootstrapping of the results, in order to estimate uncertainties of physics models fitted to this data. To assess how the uncertainties propagate to the model uncertainties, one should evaluate the model under scrutiny for each of the 200 different sets of SDMEs. Plotting all resulting lines in a single plot will create bands which reflect the influence of the uncertainties in the data on the model. This method has the great advantage that all correlations are accurately taken into account.

This table contains thinned out samples of the Markov chains used in the parameter estimation of the SDME measurements for $-(t-t_\text{0}) = 0.992\pm0.058~\text{GeV}^2/c^2$, reported in the main article. One in about 250 steps in the chain, which results in 200 different sets of SDMEs, is provided. These values should be used instead of bootstrapping of the results, in order to estimate uncertainties of physics models fitted to this data. To assess how the uncertainties propagate to the model uncertainties, one should evaluate the model under scrutiny for each of the 200 different sets of SDMEs. Plotting all resulting lines in a single plot will create bands which reflect the influence of the uncertainties in the data on the model. This method has the great advantage that all correlations are accurately taken into account.

This table contains thinned out samples of the Markov chains used in the parameter estimation of the SDME measurements for $-(t-t_\text{0}) = 1.189\pm0.058~\text{GeV}^2/c^2$, reported in the main article. One in about 250 steps in the chain, which results in 200 different sets of SDMEs, is provided. These values should be used instead of bootstrapping of the results, in order to estimate uncertainties of physics models fitted to this data. To assess how the uncertainties propagate to the model uncertainties, one should evaluate the model under scrutiny for each of the 200 different sets of SDMEs. Plotting all resulting lines in a single plot will create bands which reflect the influence of the uncertainties in the data on the model. This method has the great advantage that all correlations are accurately taken into account.

This table contains thinned out samples of the Markov chains used in the parameter estimation of the SDME measurements for $-(t-t_\text{0}) = 1.435\pm0.086~\text{GeV}^2/c^2$, reported in the main article. One in about 250 steps in the chain, which results in 200 different sets of SDMEs, is provided. These values should be used instead of bootstrapping of the results, in order to estimate uncertainties of physics models fitted to this data. To assess how the uncertainties propagate to the model uncertainties, one should evaluate the model under scrutiny for each of the 200 different sets of SDMEs. Plotting all resulting lines in a single plot will create bands which reflect the influence of the uncertainties in the data on the model. This method has the great advantage that all correlations are accurately taken into account.

This table contains thinned out samples of the Markov chains used in the parameter estimation of the SDME measurements for $-(t-t_\text{0}) = 1.761\pm0.111~\text{GeV}^2/c^2$, reported in the main article. One in about 250 steps in the chain, which results in 200 different sets of SDMEs, is provided. These values should be used instead of bootstrapping of the results, in order to estimate uncertainties of physics models fitted to this data. To assess how the uncertainties propagate to the model uncertainties, one should evaluate the model under scrutiny for each of the 200 different sets of SDMEs. Plotting all resulting lines in a single plot will create bands which reflect the influence of the uncertainties in the data on the model. This method has the great advantage that all correlations are accurately taken into account.

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$

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.

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.

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

Axis error includes +- 0.0/0.0 contribution (?////).

Bare cross section $\sigma^\mathrm{bare}(e^+e^-\to\pi^+\pi^-(\gamma_\mathrm{FSR}))$ of the process $e^+e^-\to\pi^+\pi^-$ measured using the initial state radiation method. The data is corrected concerning final state radiation and vacuum polarization effects. The final state radiation is added using the Schwinger term at born level.

Statistical covariance matrix of the bare cross section $\sigma^\mathrm{bare}(e^+e^-\to\pi^+\pi^-(\gamma_\mathrm{FSR}))$.

Statistical covariance matrix of the bare cross section $\sigma^\mathrm{bare}(e^+e^-\to\pi^+\pi^-(\gamma_\mathrm{FSR}))$.

Pion form factor $|F_\pi|^2$ measured using the initial state radiation method. The data is corrected concerning vacuum polarization effects.

Pion form factor $|F_\pi|^2$ measured using the initial state radiation method. The data is corrected concerning vacuum polarization effects.

Statistical covariance matrix of the pion form factor $|F_\pi|^2$.

Statistical covariance matrix of the pion form factor $|F_\pi|^2$.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

Dijet |k_{T,y}| distributions in p-Pb collisions in the 0-40% V0A multiplicity event class.

Dijet |k_{T,y}| distributions in p-Pb collisions in the 40-100% V0A multiplicity event class.

Dijet |k_{T,y}| distributions in p-Pb collisions in the 0-40% V0A multiplicity event class.

Dijet |k_{T,y}| distributions in p-Pb collisions in the 0-40% V0A multiplicity event class.

Dijet |k_{T,y}| distributions in p-Pb collisions in the 0-40% V0A multiplicity event class.

Mean dijet |k_{T,y}| as function of full jet pT in p-Pb collisions in the 0-40% V0A multiplicity event class.

Mean dijet |k_{T,y}| as function of associated charged jet pT in p-Pb collisions in the 0-40% V0A multiplicity event class.

The $V_2${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_2${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_2${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_2${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_3${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_3${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_3${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_3${2} and difference between the pairs at ($\eta_{\alpha}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, -$\eta_{\beta}$). The data are from 20-30% central Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $V_2^{1/2}$ difference between quadruplets at ($\eta_{\alpha}$, $\eta_{\alpha}$, $\eta_{\beta}$, $\eta_{\beta}$) and ($\eta_{\alpha}$, $\eta_{\alpha}$, -$\eta_{\beta}$, -$\eta_{\beta}$).

The two-particle cumulants, V2{2} as a function of $\eta$ for one particle while averaged over $\eta$ of the partner particle.

The two-particle cumulants, V2{2} as a function of $\eta$.

The two-particle cumulants, V2{2} as a function of $\eta$ for one particle while averaged over $\eta$ of the partner particle.

The difference in two-particle cumulants, V2{2} as a function of $\eta$ for one particle while averaged over $\eta$ of the partner particle.

No description provided.

V3{2} as a function of $\eta$.

<$V_3^2$> as a function of $\eta$.

The $\Delta\eta$-dependent component of the two-particle cumulant with $\Delta\eta$-gap, $D^-$ in Eq. (11), of the second harmonic is shown as a function of ??-gap |$\Delta\eta$| > x.

The $\Delta\eta$-dependent component of the two-particle cumulant with $\Delta\eta$-gap, $D^-$ in Eq. (11), of the third harmonic is shown as a function of ??-gap |$\Delta\eta$| > x.

The nonflow, $\sqrt{\bar{D}_2}$, $\sqrt{\delta_2}$, $\sqrt{\bar{D}_3}$ and flow, $\sqrt{v_2^2/2}$, $\sqrt{<v_2^2>}$ results are shown as a function of centrality percentile for the second harmonic.

The nonflow, $\sqrt{\bar{D}_2}$, $\sqrt{\delta_2}$, $\sqrt{\bar{D}_3}$ and flow, $\sqrt{v_2^2/2}$, $\sqrt{<v_2^2>}$ results are shown as a function of centrality percentile for the second harmonic.

The nonflow, $\sqrt{\bar{D}_2}$, $\sqrt{\delta_2}$, $sqrt{\bar{D}_3}$ and flow, $\sqrt{v_2^2/2}$, $\sqrt{<v_2^2>}$ results are shown as a function of centrality percentile for the second harmonic.

The nonflow, $\sqrt{\bar{D}_2}$, $\sqrt{\delta_2}$, $sqrt{\bar{D}_3}$ and flow, $\sqrt{v_2^2/2}$, $\sqrt{<v_2^2>}$ results are shown as a function of centrality percentile for the third harmonic.

The nonflow, $\sqrt{\bar{D}_2}$, $sqrt{\delta_2}$, $\sqrt{\bar{D}_3}$ and flow, $\sqrt{v_2^2/2}$, $\sqrt{<v_2^2>}$ results are shown as a function of centrality percentile for the third harmonic.

The relative elliptic flow fluctuation $\sigma_2/<v_2>$ centrality dependence in $\sqrt{s_{NN}}$ = 200 GeV Au+Au collisions.

Differential cross sections for inclusive jet production in charm events as a function of ETARAP(JET) for -1.6 < ETARAP(JET) < 2.2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Differential cross sections for inclusive jet production in beauty events as a function of Q**2 for 5 < Q**2 < 1000 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Differential cross sections for inclusive jet production in beauty events as a function of Q**2 for 5 < Q**2 < 1000 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Differential cross sections for inclusive jet production in beauty events as a function of XB. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Differential cross sections for inclusive jet production in charm events as a function of XB. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Double-differential cross sections for inclusive jet production in beauty events as a function of XB for 5 < Q**2 < 20 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Double-differential cross sections for inclusive jet production in beauty events as a function of XB for 20 < Q**2 < 60 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Double-differential cross sections for inclusive jet production in beauty events as a function of XB for 60 < Q**2 < 120 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Double-differential cross sections for inclusive jet production in beauty events as a function of XB for 120 < Q**2 < 400 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Double-differential cross sections for inclusive jet production in beauty events as a function of XB for 400 < Q**2 < 1000 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Double-differential cross sections for inclusive jet production in charm events as a function of XB for 5 < Q**2 < 20 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Double-differential cross sections for inclusive jet production in charm events as a function of XB for 20 < Q**2 < 60 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Double-differential cross sections for inclusive jet production in charm events as a function of XB for 60 < Q**2 < 120 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Double-differential cross sections for inclusive jet production in charm events as a function of XB for 120 < Q**2 < 400 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

Double-differential cross sections for inclusive jet production in charm events as a function of XB for 400 < Q**2 < 1000 GeV**2. The measurements are given together with their statistical and systematic uncertainties. Hadronisation and QED radiative corrections, CHAD and CRAD, respectively, are also shown.

The structure function F2(BOTTOM BOTTOMBAR) and the reduced beauty cross sections SIGR(BOTTOM BOTTOMBAR), both as functions of XB and Q**2. The measurements are given together with their statistical, systematic and extrapolation uncertainties.

The structure function F2(BOTTOM BOTTOMBAR) and the reduced beauty cross sections SIGR(BOTTOM BOTTOMBAR), both as functions of XB and Q**2. The measurements are given together with their statistical, systematic and extrapolation uncertainties.

The structure function F2(BOTTOM BOTTOMBAR) and the reduced beauty cross sections SIGR(BOTTOM BOTTOMBAR), both as functions of XB and Q**2. The measurements are given together with their statistical, systematic and extrapolation uncertainties.

The structure function F2(BOTTOM BOTTOMBAR) and the reduced beauty cross sections SIGR(BOTTOM BOTTOMBAR), both as functions of XB and Q**2. The measurements are given together with their statistical, systematic and extrapolation uncertainties.

The structure function F2(BOTTOM BOTTOMBAR) and the reduced beauty cross sections SIGR(BOTTOM BOTTOMBAR), both as functions of XB and Q**2. The measurements are given together with their statistical, systematic and extrapolation uncertainties.

The structure function F2(CHARM CHARMBAR) and the reduced charm cross sections SIGR(CHARM CHARMBAR), both as functions of XB and Q**2. The measurements are given together with their statistical, systematic and extrapolation uncertainties.

The structure function F2(CHARM CHARMBAR) and the reduced charm cross sections SIGR(CHARM CHARMBAR), both as functions of XB and Q**2. The measurements are given together with their statistical, systematic and extrapolation uncertainties.

The structure function F2(CHARM CHARMBAR) and the reduced charm cross sections SIGR(CHARM CHARMBAR), both as functions of XB and Q**2. The measurements are given together with their statistical, systematic and extrapolation uncertainties.

The structure function F2(CHARM CHARMBAR) and the reduced charm cross sections SIGR(CHARM CHARMBAR), both as functions of XB and Q**2. The measurements are given together with their statistical, systematic and extrapolation uncertainties.

The structure function F2(CHARM CHARMBAR) and the reduced charm cross sections SIGR(CHARM CHARMBAR), both as functions of XB and Q**2. The measurements are given together with their statistical, systematic and extrapolation uncertainties.

The three-point correlator, $\gamma$, as a function of centrality for Au+Au collisions at 19.6 GeV.

The three-point correlator, $\gamma$, as a function of centrality for Au+Au collisions at 11.5 GeV.

The three-point correlator, $\gamma$, as a function of centrality for Au+Au collisions at 7.7.

The two-particle correlation as a function of centrality for Au+Au collisions at 62.4 GeV.

The two-particle correlation as a function of centrality for Au+Au collisions at 39 GeV.

The two-particle correlation as a function of centrality for Au+Au collisions at 27 GeV.

The two-particle correlation as a function of centrality for Au+Au collisions at 19.6 GeV.

The two-particle correlation as a function of centrality for Au+Au collisions at 11.5 GeV.

The two-particle correlation as a function of centrality for Au+Au collisions at 7.7 GeV.

$H_{SS}-H{OS}$, as a function of beam energy for 60-80% centrality in Au+Au collisions.

$H_{SS}-H{OS}$, as a function of beam energy for 30-60% centrality in Au+Au collisions.

$H_{SS}-H{OS}$, as a function of beam energy for 10-30% centrality in Au+Au collisions.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 0.15 GeV, for centrality 30-50%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 0.15 GeV, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 10-30%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 30-50%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 10-30%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 30-50%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 0.15 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 0.15 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 10 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 10 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. Central collisions are defined to have centrality 0-10% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. Central collisions are defined to have centrality 10-30% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. Central collisions are defined to have centrality 30-50% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. Central collisions are defined to have centrality 0-10% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. Central collisions are defined to have centrality 10-30% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. Central collisions are defined to have centrality 30-50% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. Central collisions are defined to have centrality 0-10% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. Central collisions are defined to have centrality 10-30% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. Central collisions are defined to have centrality 30-50% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, as a function of the average number of participants in the collision (<Npart>), for charged jets with pT = 60-70 GeV and either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. The correspondence between <Npart> and the centrality classes is given in Table 1 of the paper, and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, as a function of the average number of participants in the collision (<Npart>), for charged jets with pT = 60-70 GeV and either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. The correspondence between <Npart> and the centrality classes is given in Table 1 of the paper, and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, as a function of the average number of participants in the collision (<Npart>), for charged jets with pT = 60-70 GeV and either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. The correspondence between <Npart> and the centrality classes is given in Table 1 of the paper, and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet pT-spectra with radius parameter R = 0.2 and 0.3 and a leading charged particle with pT > 5 GeV, for two centrality classes 0-10% and 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.