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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The observed 95% CL upper limits with one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The fiducial cross section measured in bins of the azimuthal angle difference between the mesons for events in the fiducial region with 2 Y(1S) with absolute rapidity less than 2.0.

The expected and observed 95% CL limits on the product of the cross section and branching fraction for a tetraquark. The branching fraction is the product of the branching fraction for the decay of the resonance to a Y(1S) meson and two muons, and the Y(1S) meson dimuon branching fraction.

The expected and observed 95% CL limits on the product of the cross section and branching fraction for a scalar state. The branching fraction is the product of the branching fraction for the decay of the resonance to a Y(1S) meson and two muons, and the Y(1S) meson dimuon branching fraction.

The expected and observed 95% CL limits on the product of the cross section and branching fraction for a pseudoscalar state. The branching fraction is the product of the branching fraction for the decay of the resonance to a Y(1S) meson and two muons, and the Y(1S) meson dimuon branching fraction.

The expected and observed 95% CL limits on the product of the cross section and branching fraction for a spin-2 state. The branching fraction is the product of the branching fraction for the decay of the resonance to a Y(1S) meson and two muons, and the Y(1S) meson dimuon branching fraction.

Kinematic properties and charges of the four muons selected in data in the signal region for the measurement of the YY production cross section. MUON1 and MUON2 form the dimuon pair associated to the dimuon object required at trigger-level. The invariant mass of both muon pairs is between 8.6 and 11 GeV. All events pass the trigger requirements and have four muons with a total charge of 0 and with vertex probability above 5%. The muons have pT above 3.5 GeV if they have an absolute pseudorapidity below 0.9 and 2.5 GeV otherwise. The muons pass the medium muon identification criterion.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle level.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle level.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid}\times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}p_{T}^{ll}$ measured on born level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid}\times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}\phi_{\eta}^{*}$ measured on born level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid}\times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}p_{T}^{ll}$ measured on born level for the $Z\rightarrow ee$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) and Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. for the electron momentum scale and resolution uncertainties; Elec. (Reco), Elec. (ID), Isolation, Trigger and Charge-ID denote the correlated uncertainties of the data/MC scale-factors for the electron reconstruction, identification, isolation, trigger and charge-identification efficiencies; The uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid} \times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}\phi_{\eta}^{*}$ measured on born level for the $Z\rightarrow ee$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) and Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. for the electron momentum scale and resolution uncertainties; Elec. (Reco), Elec. (ID), Isolation, Trigger and Charge-ID denote the correlated uncertainties of the data/MC scale-factors for the electron reconstruction, identification, isolation, trigger and charge-identification efficiencies; The uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid} \times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}p_{T}^{ll}$ measured on bare level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid} \times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}\phi_{\eta}^{*}$ measured on bare level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Measured combined normalized differential cross-section in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle.

Measured combined normalized differential cross-section in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle.

The distribution of events passing the selection requirements in the electron channel as a function of dilepton invariant mass $m_{ll}$ , the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of dilepton invariant mass $m_{ll}$, the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the electron channel as a function of lepton pseudorapidity $\eta$, the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of lepton pseudorapidity $\eta$, the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the electron channel as a function of dilepton transverse momentum. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of dilepton transverse momentum. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the electron channel as a function of $\phi_{\eta}^{*}$. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of $\phi_{\eta}^{*}$. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The measured normalized cross section as a function of $p_{ll}$ for the electron and muon channels and the combined result as well as their ratio together with the total uncertainties, shown as a blue band. The pull distribution between the electron and muon channels, defined as the difference between the two channels divided by the combined uncorrelated uncertainty, is also shown. The $p_{ll}$ distribution is split into linear and logarithmic scales at 30 GeV.

The measured normalized cross section as a function of $\phi_{\eta}^{*}$ for the electron and muon channels and the combined result as well as their ratio together with the total uncertainties, shown as a blue band. The pull distribution between the electron and muon channels, defined as the difference between the two channels divided by the combined uncorrelated uncertainty, is also shown.

Comparison of the normalized $p_{ll}$ distributions predicted by different computations: Pythia8 with the AZ tune, Powheg+Pythia8 with the AZNLO tune, Sherpa v2.2.1 and RadISH with the Born level combined measurement. The uncertainties of the measurement are shown as vertical bars and uncertainties of the Sherpa and RadISH predictions are indicated by the coloured bands.

Comparison of the normalized $\phi_{\eta}^{*}$ distributions predicted by different computations: Pythia8 with the AZ tune, Powheg+Pythia8 with the AZNLO tune, Sherpa v2.2.1 and RadISH with the Born level combined measurement. The uncertainties of the measurement are shown as vertical bars and uncertainties of the Sherpa and RadISH predictions are indicated by the coloured bands.

Comparison of the normalized $p_{ll}$ distribution in the range $p_{ll}$ > 10 GeV. The Born level combined measurement is compared with predictions by Sherpa v2.2.1, fixed-order NNLOjet and NNLOjet supplied with NLO electroweak corrections. The uncertainties in the measurement are shown as vertical bars and the uncertainties in the predictions are indicated by the coloured bands.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at dressed level.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at dressed level.

The measured fiducial cross section vs $m(\ell\ell\gamma)$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

The measured fiducial cross section vs $p_{\mathrm{T}}^{\ell\ell\gamma}$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

The measured fiducial cross section vs $\Delta\phi(ll,\gamma)$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

The measured fiducial cross section vs $p_{\mathrm{T}}^{\ell\ell\gamma}/m(\ell\ell\gamma)$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Fig. A1: Number density $N_{ch}$ (left) and $\\Sigma p_{T}$ (right) distributions as a function of pleadingT and the comparisons to MC predictions in Toward (top), Transverse (middle), and Away (bottom) regions for $p_{T}^{track} >$ 0.5 GeV/$c$. The shaded areas in the upper panels represent the systematic uncertainties and vertical error bars indicate statistical uncertainties. In the lower panels, the shaded areas are the sum in quadrature of statistical and systematic uncertainties from the upper panels. No uncertainties are given for the MC calculations.

Fig. A2: Number density $N_{ch}$ (left) and $\\Sigma p_{T}$ (right) distributions as a function of pleadingT and the comparisons to MC predictions in Toward (top), Transverse (middle), and Away (bottom) regions for $p_{T}^{track} >$ 1.0 GeV/$c$. The shaded areas in the upper panels represent the systematic uncertainties and vertical error bars indicate statistical uncertainties. In the lower panels, the shaded areas are the sum in quadrature of statistical and systematic uncertainties from the upper panels. No uncertainties are given for the MC calculations.

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class IV+V, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class VI, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class VII, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class VIII, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class IX, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class X, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class I

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class II

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class III

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class IV+V

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VI

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VII

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VIII

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class IX

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class X

$\phi$ transverse momentum spectrum - V0M multiplicity class I

$\phi$ transverse momentum spectrum - V0M multiplicity class II

$\phi$ transverse momentum spectrum - V0M multiplicity class III

$\phi$ transverse momentum spectrum - V0M multiplicity class IV

$\phi$ transverse momentum spectrum - V0M multiplicity class V

$\phi$ transverse momentum spectrum - V0M multiplicity class VI

$\phi$ transverse momentum spectrum - V0M multiplicity class VII

$\phi$ transverse momentum spectrum - V0M multiplicity class VIII

$\phi$ transverse momentum spectrum - V0M multiplicity class IX

$\phi$ transverse momentum spectrum - V0M multiplicity class X

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class I

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class II

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class III

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class IV

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class V

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VI

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VII

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VIII

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class IX

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class X

K$^{*0}$ average transverse momentum vs multiplicity - V0M multiplicity classes

$\phi$ average transverse momentum vs multiplicity - V0M multiplicity classes

K$^{*0}$ transverse momentum integrated yield vs multiplicity - V0M multiplicity classes, average of particle and antiparticle

$\phi$ transverse momentum integrated yield vs multiplicity - V0M multiplicity classes

K$^{*0}$/K$_{\mathrm{S}}^{0}$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes, average of particle and antiparticle

$\phi$/K$_{\mathrm{S}}^{0}$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

$(\Xi^{-} + \bar{\Xi}^{+})$/$\phi$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

$\phi$/$(\pi^{+} + \pi^{-})$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

K$^{*0}$/K$_{\mathrm{S}}^{0}$ yield ratio vs transverse momentum - V0M multiplicity class II, average of particle and antiparticle

K$^{*0}$/K$_{\mathrm{S}}^{0}$ yield ratio vs transverse momentum - V0M multiplicity class X, average of particle and antiparticle

$\phi$/K$_{\mathrm{S}}^{0}$ yield ratio vs transverse momentum - V0M multiplicity class II

$\phi$/K$_{\mathrm{S}}^{0}$ yield ratio vs transverse momentum - V0M multiplicity class X

K$^{*0}$/K$_{\mathrm{S}}^{0}$ double yield ratio vs transverse momentum - V0M multiplicity class II / V0M multiplicity class X, average of particle and antiparticle

$\phi$/K$_{\mathrm{S}}^{0}$ double yield ratio vs transverse momentum - V0M multiplicity class II / V0M multiplicity class X

Significance of K$^{*0}$/K$_{\mathrm{S}}^{0}$ double yield ratio vs transverse momentum - V0M multiplicity class II / V0M multiplicity class X

Significance of $\phi$/K$_{\mathrm{S}}^{0}$ double yield ratio vs transverse momentum - V0M multiplicity class II / V0M multiplicity class X

($\rm{K}^{*0}$+$\bar{\rm{K}}^{*0}$)/(K$^{+}$+K$^{-}$) transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

($\rm{K}^{*0}$+$\bar{\rm{K}}^{*0}$)/($\pi^{+}$+$\pi^{-}$) transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

2$\phi$/(K$^{+}$+K$^{-}$) transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

$\phi$/($\rm{P}$ + $\bar{\rm{P}}$) transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

$\phi$/$(\Lambda~+~\bar{\Lambda})$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

$(\Omega^{-}~+~\bar{\Omega}^{+})$/$\phi$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes