Cross sections (fb) at large Z $p_{T}$ values in the Z -> $\ell\ell$ and Z -> $\nu\nu$ channels, and their combination. The theoretical predictions from Madgraph at NLO in QCD and corrected to NLO in electroweak using the NNPDF 3.0 are also reported. With the exception of the largest Z $p_{T}$ bin, the statistical uncertainties in the measurements are much smaller than the systematic uncertainties. Both measurements and predictions correspond to $\sigma \mathcal{B}(Z -> \ell\ell)$, where $\sigma$ is the total fiducial cross section, $\mathcal{B}$ is the branching fraction, and $\ell$ is a charged lepton. The $Z -> \nu\nu$ measurement corresponds to $\sigma \mathcal{B}(Z -> \ell\ell)/\mathcal{B}(Z -> \nu\nu)$.

Cross sections normalized to the total cross section measurements at high Z $p_{T}$ values in the Z -> $\ell\ell$ and Z -> $\nu\nu$ channels, and in their combination. The uncertainty includes both statistical and systematic uncertainties.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on charged leptons.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on charged leptons.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on neutrinos.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on neutrinos.

The relative statistical and systematic uncertainties from various sources for the absolute cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

The relative statistical and systematic uncertainties from various sources for the relative cross section measurements in bins of Z $p_{T}$ on neutrinos and charged leptons.

p LambdaBar plus pbar Lambda correlation function for centrality 10-20% from Pb-Pb collisions at 2.76 TeV

Lambda LambdaBar plus pbar Lambda correlation function for centrality 10-20% from Pb-Pb collisions at 5.02 TeV

Lambda LambdaBar plus pbar Lambda correlation function for centrality 10-20% from Pb-Pb collisions at 2.76 TeV

Differential fiducial cross section of the Z boson $|y|$ in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the leading b-jet $|y|$ in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta \phi$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta y$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta R$ between Z boson and leading $b$-jet in events with $Z(\rightarrow ll)\ge+1$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta \phi$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta y$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $\Delta R$ between the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the invariant mass of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the Z boson $p_{\text{T}}$ in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the $p_{\text{T}}$ of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Differential fiducial cross section of the ratio between the $p_{\text{T}}$ and the invariant mass of the first two leading $b$-jets in events with $Z(\rightarrow ll)\ge+2$ b-jets. The statistical uncertainties and the individual components of systematic uncertainty are given in each bin. Statistical uncertainties are bin-to-bin uncorrelated.

Event yields for a simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ sample with $m_a = 17.5\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. Cut 0 corresponds to the initial number of events. Cut 1 requires the single lepton trigger. Cut 2 requires 2 identified leptons. Cut 3 requires the Z-boson mass window. Cut 4 requires 2 reconstructed $a\to b\bar{b}$ candidates. Cut 5a requires 2 identified $a\to b\bar{b}$ candidates in the 1HPC1LPC region. Cut 6a requires the 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region to be inside the Higgs mass window. Cut 5b requires 2 identified $a\to b\bar{b}$ candidates in the 2HPC region. Cut 6b requires the 2 $a\to b\bar{b}$ candidates in the 2HPC region to be inside the Higgs mass window.

Background yield table for Z+jets, $t\bar{t}$, and rare sources. Observed data yield. Signal $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ yield with $m_a = 20\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$, with a branching ratio set to 1 for the $H \rightarrow aa$ decay, whereas the ATLAS figure attached to this entry instead uses the upper-limit branching ratio (smaller than 1). The table includes the yields in two signal regions with leptons consistent with an on-shell Z-boson decay, one with 2 $a\to b\bar{b}$ candidates in the 2HPC region and one with 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region. The table also includes the yields in four control regions, one with leptons consistent with an on-shell Z-boson decay and 2 $a\to b\bar{b}$ candidates in the Low Purity Category (LPC), and three others where the leptons are not consistent an on-shell Z-boson decay.

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.

Detailed breakdown of systematic uncertainties for the combined measurement of muon, electron central and electron central-forward channels. Common systematic uncertainty on the luminosity measurment of 1.8% is not included. The source 'sys,uncor' represents bin-to-bin uncorrelated systematic uncertainty. The cross sections are given at the Born QED level.

Detailed breakdown of systematic uncertainties for the combined measurement, integerated in cos theta_CS (differential in y, Mll) Common systematic uncertainty on the luminosity measurment of 1.8% is not included. The source 'sys,uncor' represents bin-to-bin uncorrelated systematic uncertainty. The cross sections are given at the Born QED level.

Detailed breakdown of systematic uncertainties for the combined measurement, integerated in cos theta_CS and y (differential in Mll) Common systematic uncertainty on the luminosity measurment of 1.8% is not included. The source 'sys,uncor' represents bin-to-bin uncorrelated systematic uncertainty. The cross sections are given at the Born QED level.

Powheg based prediction for AFB in the central-central fiducial phase space, as reported in Fig 16 of the paper. Powheg prediction is corrected to NNLO QCD and NLO EWK, as described in the paper. PDF uncertainties are computed using CT10 PDF set scaled to 68%.

Powheg based prediction for AFB in the central-fiducial fiducial phase space, as reported in Fig 17 of the paper. Powheg prediction is corrected to NNLO QCD and NLO EWK, as described in the paper. PDF uncertainties are computed using CT10 PDF set scaled to 68%.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, where all slepton flavors are mass degenerate (co-NLSP).

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the selectron signal model. Selectron ($\tilde{e}_{L, R}$) refers to the scalar superpartners of left- and right-handed electrons.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the left-handed electrons, $\tilde{e}_L$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the right-handed electrons, $\tilde{e}_R$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the smuon signal model. Smuon ($\tilde{\mu}_{L, R}$) refers to the scalar superpartners of left- and right-handed muons.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the left-handed muons, $\tilde{\mu}_L$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the superpartners of the right-handed muons, $\tilde{\mu}_R$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for the stau signal model. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin\theta_{\tilde\tau}=0.95$.

Upper limits on observed signal cross section, $\sigma_\text{obs}^\text{95}$, for $\tilde{\tau}_L$ production, where $\tilde{\tau}_L$ is pure-state superpartner of the left-handed $\tau$.

The expected and observed yields in the signal regions. Combined statistical and systematic uncertainties are presented. Estimates are truncated at 0 if the size of measured systematic uncertainties would yield a negative result.

Reconstruction efficiency as a function of $|d_{0}|$ and $p_\text{T}$ for signal electrons. Monte Carlo samples with $\tilde{e}$ or $\tilde{\mu}$ with mass 400 GeV and 1 ns lifetime were used. Efficiency is defined as the number of leptons passing all signal requirements and matched to generator-level muons divided by the number of generator level leptons with $p_\text{T} > 65$ GeV, $|d_{0}| >$ 3 mm, and |$\eta$| $<$ 2.47 for electrons. No event-level selections are made. Reconstructed leptons must pass all quality criteria, including the cosmic veto. Electron selection scale factors are included on the reconstructed leptons. The overflow is not shown in these plots.

Reconstruction efficiency as a function of $|d_{0}|$ and $p_\text{T}$ for signal muons. Monte Carlo samples with $\tilde{e}$ or $\tilde{\mu}$ with mass 400 GeV and 1 ns lifetime were used. Efficiency is defined as the number of leptons passing all signal requirements and matched to generator-level muons divided by the number of generator level leptons with $p_\text{T} > 65$ GeV, $|d_{0}| >$ 3 mm, and |$\eta$| $<$ 2.5 for muons. No event-level selections are made. Reconstructed leptons must pass all quality criteria, including the cosmic veto. Muon selection scale factors are included on the reconstructed leptons. The overflow is not shown in these plots.

Acceptance for $\tilde{e}$ of various masses and lifetimes in SR-$ee$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\mu}$ of various masses and lifetimes in SR-$\mu\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Efficiency for $\tilde{e}$ of various masses and lifetimes in SR-$ee$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\mu}$ of various masses and lifetimes in SR-$\mu\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$ee$ from $\tilde{e}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47 for electrons, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$\mu\mu$ from $\tilde{\mu}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 for muons, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$ee$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$e\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Acceptance for $\tilde{\tau}$ of various masses and lifetimes in SR-$\mu\mu$. Acceptance is defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$ee$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth electrons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.47, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$e\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Efficiency for $\tilde{\tau}$ of various masses and lifetimes in SR-$\mu\mu$. Efficiency is defined as the number of selected reconstruced events divided by the acceptance. To be accepted, events are required to have at least 2 truth muons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5, $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. To be selected, events must satisfy all signal region requirements.

Acceptance as a function of the generator-level $p_\text{T}$ of the leading and subleading lepton in SR-$e\mu$ from $\tilde{\tau}$ decays. The plot is made from signal Monte Carlo events with $\tilde{\ell}$ with mass of 400 GeV and lifetime of 1 ns. To be accepted, events are required to have at least 2 truth leptons with $p_\text{T} > 65$ GeV, |$\eta$| $<$ 2.5 (2.47) for muons (electrons), $|d_{0}| > 3$ mm, and $\Delta R_{\ell\ell} > 0.2$. Events are also required to fall into one of the acceptance regions of the triggers used. At generator level, events must have one of the following: one electron with $p_\text{T} >$ 160 GeV, 2 electrons each with $p_\text{T} >$ 60 GeV, or 1 muon with $p_\text{T} >$ 60 GeV and $|\eta| < 1.07$. The overflow is not shown in these plots.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane where all slepton flavors and chiralities are mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane where all slepton flavors and chiralities are mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L, R}$) refer to the scalar superpartners of left- and right-handed electrons, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L, R}$) refer to the scalar superpartners of left- and right-handed electrons, which are assumed to be mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L}$) refer to the scalar superpartners of left-handed electrons.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{L}$) refer to the scalar superpartners of left-handed electrons.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{R}$) refer to the scalar superpartners of right-handed electrons. Purple denotes the region excluded by LEP.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$ee$ targeting selectron production. Selectrons ($\tilde{e}_{R}$) refer to the scalar superpartners of right-handed electrons. Purple denotes the region excluded by LEP.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L, R}$) refer to the scalar superpartners of left- and right-handed muons, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L, R}$) refer to the scalar superpartners of left- and right-handed muons, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L}$) refer to the scalar superpartners of left-handed muons.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{L}$) refer to the scalar superpartners of left-handed muons.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{R}$) refer to the scalar superpartners of right-handed muons. Purple denotes the region excluded by LEP.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\ell})$ plane in SR-$\mu\mu$ targeting smuon production. Smuons ($\tilde{\mu}_{R}$) refer to the scalar superpartners of right-handed muons. Purple denotes the region excluded by LEP.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau})$ plane. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin(\theta_{\tilde\tau})=0.95$, which are assumed to be mass degenerate.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau})$ plane. Staus, $\tilde{\tau}_{1,2}$ are the mixed states of the superpartners of the left- and right-handed $\tau$ leptons, with mixing angle $\sin(\theta_{\tilde\tau})=0.95$, which are assumed to be mass degenerate.

Observed 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau}_L)$ plane, where $\tilde{\tau}_L$ is the pure-state super-partner of the left-handed $\tau$. Purple denotes the region excluded by LEP. This result does not present signficant sensitivity to the pure-state superpartner of the right-handed $\tau$.

Expected 95% CL exclusion sensitivity. The limit is displayed in the lifetime vs. $m(\tilde{\tau}_L)$ plane, where $\tilde{\tau}_L$ is the pure-state super-partner of the left-handed $\tau$. Purple denotes the region excluded by LEP. This result does not present signficant sensitivity to the pure-state superpartner of the right-handed $\tau$.

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

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Reconstructed VLQ mass distributions for simulated signal events with a generated VLQ mass $m_{B} = 1200\GeV$. A moderate requirement of $\chi^{2}$/ndf < 2$ is applied to the events. Mass distributions for 4-jet (left), 5-jet (center), and 6-jet (right) events are shown for the three decay modes: bHbH (upper row), bHbZ (middle row), and bZbZ (lower row).

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distribution of $\chi^{2}$/ndf for the best jet combination for simulated 1200\GeV VLQ events (red histogram) and data (black points), for 4-jet events (left), 5-jet events (center), and 6-jet events (right). The simulated signal events and data events are normalized to the same value within the displayed $\chi^{2}$/ndf range.

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Distributions of the average reconstructed mass of VLQ candidates for the jet combination with the least $\chi^{2}$ in 4-jet (left), 5-jet (center), and 6-jet (right) multiplicity events. The red lines show the exponential fit in the range 1000-2000 GeV. The lower panels show the fractional difference, (data-fit)/fit

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on the average reconstructed VLQ mass in the control region, 12 < $\chi^{2}$/ndf < 48 for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

Dependence of the BJTF on $\chi^{2}$/ndf in the low-mass VLQ region, for 4-jet (left column), 5-jet, (center column) and 6-jet (right column) multiplicities and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes. The data is shown in black points, and the linear fit and its uncertainty are shown as the red line and the pale red band.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in data events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

The reduction factor in $m_B$ = 1200 GeV VLQ signal events, for 4-jet (left column), 5-jet (center column), and 6-jet (right column) multiplicities, and for the bHbH (upper row), bHbZ (middle row), and bZbZ (lower row) event modes.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Data (black points), expected background (solid blue histogram), and expected background plus a VLQ signal for different VLQ masses (colored lines). Left column: 4-jet events, center column: 5-jet events, right column: 6-jet events; upper row: bHbH, middle row: bHbZ,lower row: bZbZ. For the signal, a 100% B -> bH branching fraction is assumed. The hatched regions for the background and background plus signal distributions indicate the systematic uncertainties. All three data-taking years are combined.

Expected limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}$( B -> bH ) and $\mathcal{B}$( B -> bZ)

Observed limits on the VLQ mass at 95% CL as a function of the branching fractions $\mathcal{B}$( B -> bH ) and $\mathcal{B}$( B -> bZ)

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.

The 95% confidence limit on the cross section for VLQ pair production as a function of VLQ mass for three branching fraction hypotheses: B( B -> bH ) = 100% (upper left), B( B -> bZ ) = 100% (upper right), and and B( B -> bH ) = B( B -> bZ ) = 50% (lower). The solid black line indicates the observed limit and the dashed line indicates the expected limit with 1 sigma (green band) and 2 sigma (yellow band) uncertainties. The theoretical cross section and its uncertainty are shown as the red line and pale red band.

Transverse momentum spectrum of $K^{*0} + \overline{K^{*0}}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $\phi$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $p + \overline{p}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $\Lambda + \overline{\Lambda}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $\Xi^{-} + \overline{\Xi}^{+}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $\Omega^{-} + \overline{\Omega}^{+}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 13 TeV. The normalization uncertainty of $\pm2.6\%$ is excluded.

Transverse momentum spectrum of $K^{0}_{S}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 7 TeV. The normalization uncertainty of $^{+7.3}_{-3.5}\%$ is excluded.

Transverse momentum spectrum of $\Lambda + \overline{\Lambda}$ measured at midrapidity ($|y|<0.5$) in inelastic pp collisions at $\sqrt{s}$ = 7 TeV. The normalization uncertainty of $^{+7.3}_{-3.5}\%$ is excluded.

Ratios of the transverse-momentum spectra of $\pi^{+} + \pi^{-}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $K^{+} + K^{-}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $K^{0}_{S}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $K^{*0} + \overline{K^{*0}}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $\phi$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $p + \overline{p}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $\Lambda + \overline{\Lambda}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $\Xi^{-} + \overline{\Xi}^{+}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

Ratios of the transverse-momentum spectra of $\Omega^{-} + \overline{\Omega}^{+}$ in inelastic pp collisions at $\sqrt{s}=13$ TeV to those for $\sqrt{s}={7}$ TeV. The normalization uncertainty of $^{+10.8}_{-6.3}\%$ is excluded.

$(K^{+} + K^{-})/(\pi^{+}+\pi^{-})$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$(K^{*0} + \overline{K^{*0}})/(\pi^{+}+\pi^{-})$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$K^{0}_{S}/(\pi^{+}+\pi^{-})$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$\phi/(\pi^{+}+\pi^{-})$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$\phi/(\pi^{+}+\pi^{-})$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 7 TeV.

Average transverse momentum of $\pi^{+}+\pi^{-}$ as a function of the center-of-mass energy.

Average transverse momentum of $K^{+} + K^{-}$ as a function of the center-of-mass energy.

Average transverse momentum of $K^{0}_{S}$ as a function of the center-of-mass energy.

Average transverse momentum of $K^{*0} + \overline{K^{*0}}$ as a function of the center-of-mass energy.

Average transverse momentum of $\phi$ as a function of the center-of-mass energy.

Average transverse momentum of $p + \overline{p}$ as a function of the center-of-mass energy.

Average transverse momentum of $\Lambda+\overline{\Lambda}$ as a function of the center-of-mass energy.

Average transverse momentum of $\Xi^{-} + \overline{\Xi}^{+}$ as a function of the center-of-mass energy.

Average transverse momentum of $\Omega^{-} + \overline{\Omega}^{+}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $K^{+} + K^{-}$ to $\pi^{+} + \pi^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $p + \overline{p}$ to $\pi^{+} + \pi^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda+\overline{\Lambda}$ to $\pi^{+} + \pi^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $\Xi^{-} + \overline{\Xi}^{+}$ to $\pi^{+} + \pi^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $\Omega^{-} + \overline{\Omega}^{+}$ to $\pi^{+} + \pi^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $\phi$ to $K^{+} + K^{-}$ as a function of the center-of-mass energy.

Ratio of $p_{\rm T}$-integrated yields of $K^{*0} + \overline{K^{*0}}$ to $K^{+} + K^{-}$ as a function of the center-of-mass energy.

$(p + \overline{p})$/($\pi^{+}+\pi^{-}$) particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$(\Lambda + \overline{\Lambda})/K^{0}_{S}$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$(\Omega^{-} + \overline{\Omega}^{+})/\phi$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$(\Xi^{-} + \overline{\Xi}^{+})/\phi$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$(\Omega^{-} + \overline{\Omega}^{+})/\phi$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 7 TeV.

$(\Xi^{-} + \overline{\Xi}^{+})/\phi$ particle ratio as a function of $p_{\rm T}$ measured in pp collisions at $\sqrt{s}$ = 7 TeV.

The balance function in terms of $\Delta \eta$ for all charged particles with $0.2 < p_{T} < 2.0$ GeV/$c$ from central Au+Au collisions (0-5%) for $\sqrt{s_{NN}}=27$ GeV. The data are the measured balance functions corrected by subtracting balance functions calculated using mixed events. Also shown are balance functions calculated using shuffled events.

The balance function in terms of $\Delta \eta$ for all charged particles with $0.2 < p_{T} < 2.0$ GeV/$c$ from central Au+Au collisions (0-5%) for $\sqrt{s_{NN}}=39$ GeV. The data are the measured balance functions corrected by subtracting balance functions calculated using mixed events. Also shown are balance functions calculated using shuffled events.

The balance function in terms of $\Delta \eta$ for all charged particles with $0.2 < p_{T} < 2.0$ GeV/$c$ from central Au+Au collisions (0-5%) for $\sqrt{s_{NN}}=62.4$ GeV. The data are the measured balance functions corrected by subtracting balance functions calculated using mixed events. Also shown are balance functions calculated using shuffled events.

The balance function in terms of $\Delta \eta$ for all charged particles with $0.2 < p_{T} < 2.0$ GeV/$c$ from central Au+Au collisions (0-5%) for $\sqrt{s_{NN}}=200$ GeV. The data are the measured balance functions corrected by subtracting balance functions calculated using mixed events. Also shown are balance functions calculated using shuffled events.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Energy dependence of the balance function widths compared with the widths of the balance functions calculated using shuffled events. Also shown are the balance function widths calculated using UrQMD. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution. Error bars represent the statistical error and the shaded bands represent the systematic error.

Balance function widths for the most central events ($0-5\%$) compared with balance function widths calculated using shuffled events. Also shown are balance function widths calculated using UrQMD and shuffled UrQMD events. The dashed line represents the width of the balance function calculated using shuffled events for a constant $dN/d\eta$ distribution.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

Acceptance-corrected balance function widths for Au+Au measured over the range $0.1 < \Delta \eta < 1.6$ normalized to the most peripheral centrality bin compared with similar results from Pb+Pb collisions from ALICE. Only statistical errors are shown. Lines represent fits of the form $a + b(N_{part})^{0.01}$.

The $m_{W}^{T}$ distribution from data (points) and simulation (shaded histograms) in the 2j1t (left) and 3j1t (right) categories for the muon (upper) and electron (lower) channels. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 2j1t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and b jet momenta system (right), shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distributionfrom the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 2j1t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and b jet momenta system (right), shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distributionfrom the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 2j1t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and b jet momenta system (right), shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distributionfrom the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 2j1t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and b jet momenta system (right), shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distributionfrom the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j1t category: the $p_{T}^{miss}$ in the transverse plane (left) and the value of the MVA b tagger discriminator when applied to the extra jet (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points andthe hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j1t category: the $p_{T}^{miss}$ in the transverse plane (left) and the value of the MVA b tagger discriminator when applied to the extra jet (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points andthe hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j1t category: the $p_{T}^{miss}$ in the transverse plane (left) and the value of the MVA b tagger discriminator when applied to the extra jet (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points andthe hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j1t category: the $p_{T}^{miss}$ in the transverse plane (left) and the value of the MVA b tagger discriminator when applied to the extra jet (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points andthe hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j2t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and non-b-tagged jet system (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j2t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and non-b-tagged jet system (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j2t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and non-b-tagged jet system (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distributions of the two most discriminating variables from data (points) and simulation (shaded histograms) in the 3j2t category: the $|\eta|$ of the non-b-tagged jet $|\eta_{j^{'}}|$ (left) and the invariant mass of lepton and non-b-tagged jet system (right) are shown for the muon (upper) and electron (lower) channels, respectively. The vertical lines on the points and the hatched bands show the experimental and MC statistical uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the MC prediction.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Distribution of the multivariate discriminators, comparing data to simulation normalised after the fit procedure, for the muon channel on the left and for the electron channel on the right, for 2j1t (upper), 3j1t (middle), and 3j2t (lower). The vertical lines on the points and the hatched bands show the experimental and fit uncertainties, respectively. The expected distribution from the STq,b+STb,q processes (multiplied by a factor of 1000) is shown by the solid blue line. The lower panels show the ratio of the data to the fit.

Values of the third-row elements of the CKM matrix inferred from low-energy measurements with the respective values of the top quark decay branching fractions. The q in bs{V_{tq}} and \mathcal{B}(t o Wq) in the first column refers to b, s, and d quarks, according to the quark label shown in the header row.

For each of the production and decay vertices, the cross section times branching fraction for the corresponding signal process from simulation. The uncertainties shown includethose from the factorisation and renormalisation scales, the PDFs, and any experimental uncer-tainties, where appropriate.

The sources and relative values in percent of the systematic uncertainty in the measurement of the STb,b cross section. The uncertainties are broken up into profiled and nonprofiled sources.

The signal strengths obtained under SM assumptions.

Results obtained under SM assumptions.

Results obtained under first BSM scenario.

Results obtained under second BSM scenario.

Confidence limits on associated yield with thrust coordinates as a function of avg N_trk^corr.

The correlation matrix for the njet measurement. The last bin is inclusive.

The fiducial integrated signal strength and cross section The uncertainty breakdown is given in terms of statistical (stat), experimental (exp), theoretical uncertainties on the background (bkg) and on the signal (sig), and the luminosity uncertainty (lumi). The fiducial cross section and its full uncertainty in each bin are also given.

Invariant $p_{T}$-differential spectra of $\Omega- + \bar{\Omega+}$ in p+p and p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in -1.8 < $y_{CM}$ < 0

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0 < $y_{CM}$ < 1.8

Invariant $p_{T}$-differential spectra of ${K_{0}}^{S}$ in p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

Invariant $p_{T}$-differential spectra of $\Lambda + \bar{\Lambda}$ in p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.3 < $y_{CM}$ < 0.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.8 < $y_{CM}$ < 1.3

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 1.3 < $y_{CM}$ < 1.8

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.

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

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.

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

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

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

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

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

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

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

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.

Polarization of $\Lambda$ and $\bar\Lambda$ as a function of observed charge asymmetry $A_{ch}$ normalized with its RMS $\sigma_{A_{ch}}$ for the $20\%–60\%$ centrality bin in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. Open boxes and vertical lines show systematic and statistical uncertainties. Solid and dashed lines show linear fitting functions.

The invariant yield versus transverse momentum for |y| < 1 in 0-60% centrality in Au+Au collisions (solid circles). The results are compared to high-$p_T$ (3 < $p_T$ < 10 GeV/c) results from STAR [9] (solid squares) and PHENIX data [8] (open squares). Also shown is the yield in 0-60% centrality Cu+Cu collisions for low $p_T$ (solid triangles) and high $p_T$ [48] (open triangles). The models are described in the text [49, 50]. The $J/\psi$ cross section in p+p collisions is also shown at STAR (stars) and PHENIX [51] (diamonds), and the scale is indicated on the right axis.

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 0-60% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 0-20% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 20-40% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 40-60% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ invariant yield and nuclear modification factor as a function of centrality for |y| < 1 in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $J/\psi$ invariant yield $Bd^2N/(2\pi p_Tdydp_T) [(GeV/c)^{-2}]$ and nuclear modification factor as a function of transverse momentum for |y| < 1 in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $J/\psi$ invariant yield and nuclear modification factor for |y| < 1 in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV.

v2 as a function of pT for Pb-Pb collisions at \sqrt{s_NN} = 5.02 TeV and centrality 20-30%.

v2 as a function of pT for Pb-Pb collisions at \sqrt{s_NN} = 5.02 TeV and centrality 30-40%.

v2 as a function of pT for Pb-Pb collisions at \sqrt{s_NN} = 5.02 TeV and centrality 40-50%.

v2 as a function of pT for Pb-Pb collisions at \sqrt{s_NN} = 5.02 TeV and centrality 50-60%.

v2 as a function of pT for Pb-Pb collisions at \sqrt{s_NN} = 5.02 TeV and centrality 60-70%.

v3 as a function of pT for Pb-Pb collisions at \sqrt{s_NN} = 5.02 TeV and centrality 0-20%.

v3 as a function of pT for Pb-Pb collisions at \sqrt{s_NN} = 5.02 TeV and centrality 20-40%.

v3 as a function of pT for Pb-Pb collisions at \sqrt{s_NN} = 5.02 TeV and centrality 40-60%.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

The overall flow harmonic $v_4\{\Psi_{4}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

The overall flow harmonic $v_5\{\Psi_{5}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

The overall flow harmonic $v_6\{\Psi_{6}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

The overall flow harmonic $v_6\{\Psi_{6}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

The overall flow harmonic $v_7\{\Psi_{7}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

The overall flow harmonic $v_4\{\Psi_{4}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

The overall flow harmonic $v_5\{\Psi_{5}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

The overall flow harmonic $v_6\{\Psi_{6}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

The overall flow harmonic $v_6\{\Psi_{6}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

The overall flow harmonic $v_7\{\Psi_{7}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 5.02 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 5.02 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 2.76 TeV as a function of PT in the 0-20% centrality range.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 2.76 TeV as a function of PT in the 20-60% centrality range.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_4\{\Psi_{22}\}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_5\{\Psi_{23}\}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_6\{\Psi_{222}\}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_6\{\Psi_{33}\}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Mixed higher-order flow harmonic $v_7\{\Psi_{223}\}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 5.02 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{422}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{523}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{6222}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{633}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Nonlinear response coefficient $\chi_{7223}$ from the scalar-product method at 2.76 TeV as a function of centrality.

Distribution of $E^{miss}_T$ in data and for the expected SM background in the Two-Electron CR after performing the 'simplified shape fit'. The $E^{miss}_$ calculation in this CR does not include electron contribution. The error bars are statistical, and the dashed band includes statistical and systematic uncertainties determined by the fit. The lower panel shows the ratio of data to expected background event yields.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.