The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 9 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 11 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 12 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 13 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for 14 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $1\ell$ category for at least 15 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 4 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 5 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 6 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for 9 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold regions defined for the EWK analysis in the $2\ell^{\mathrm{sc}}$ category for at least 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 15 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 40 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 12 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 60 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 11 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 80 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 100 GeV jet $p_{\mathrm{T}}$ threshold region in the $1\ell$ category for at least 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 20 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 10 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 40 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 8 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 60 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 80 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 7 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

The observed data event yields and the corresponding estimates for the backgrounds in the different $b$-jet multiplicity bins for the 100 GeV jet $p_{\mathrm{T}}$ threshold region in the $2\ell^{\mathrm{sc}}$ category for at least 6 jets. The background is estimated by including all bins in the fit. All uncertainties, which may be correlated across the bins, are included in the total background uncertainty.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $1\ell$ category, as well as the observed and expected 95% CL model-independent upper limits on product of cross-section, acceptance and efficiency (in fb). The parameters of the model are determined in a fit to a reduced set of bins, corresponding to the model-independent fit discussed in the text.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $2\ell^{\mathrm{sc}}$ category, as well as the observed and expected 95% CL model-independent upper limits on product of cross-section, acceptance and efficiency (in fb). The parameters of the model are determined in a fit to a reduced set of bins, corresponding to the model-independent fit discussed in the text.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Expected 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Observed 95% CL exclusion contour for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Expected 95% CL exclusion contour for the stop pair production model with bino LSP.

Observed 95% CL exclusion contour for the stop pair production model with bino LSP.

Expected 95% CL exclusion contour for the stop pair production model with higgsino LSP.

Observed 95% CL exclusion contour for the stop pair production model with higgsino LSP.

Expected 95% CL exclusion contour for the stop pair production model with wino LSP.

Observed 95% CL exclusion contour for the stop pair production model with wino LSP.

Expected 95% CL excluded cross section for the RPV model with electroweakino prodiction with higgsino LSP hypothesis.

Observed 95% CL excluded cross section for the RPV model with electroweakino prodiction with higgsino LSP hypothesis.

Expected 95% CL excluded cross section for the RPV model with electroweakino prodiction with wino LSP hypothesis.

Observed 95% CL excluded cross section for the RPV model with electroweakino prodiction with wino LSP hypothesis.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $1\ell$ category with 20 and 40 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $1\ell$ category with 60, 80 and 100 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $2\ell^{\mathrm{sc}}$ category with 20 and 40 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Data event yields compared with the expected contributions from relevant background sources, in the discovery signal regions defined for the $2\ell^{\mathrm{sc}}$ category with 60, 80 and 100 GeV jet $p_{\mathrm{T}}$ thresholds. The uncertainties across backgrounds can exhibit strong anticorrelations.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with bino LSP.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with wino LSP.

Expected 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Observed 95% CL excluded cross section for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP.

Expected 95% CL excluded cross section for the stop pair production model with bino LSP.

Observed 95% CL excluded cross section for the stop pair production model with bino LSP.

Expected 95% CL excluded cross section for the stop pair production model with wino LSP.

Observed 95% CL excluded cross section for the stop pair production model with wino LSP.

Expected 95% CL excluded cross section for the stop pair production model with higgsino LSP.

Observed 95% CL excluded cross section for the stop pair production model with higgsino LSP.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model in the $1\ell$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model in the $1\ell$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP in the $1\ell$ category.

Acceptance and efficiency for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with higgsino LSP in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the stop model with higgsino LSP in the $1\ell$ category.

Acceptance and efficiency for the stop model with higgsino LSP in the $2\ell^{\mathrm{sc}}$ category.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 4 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 5 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 6 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 7 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 8 jets, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $2\ell^{\mathrm{sc}}$ category, considering $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 4 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 5 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 6 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 7 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $1\ell$ category with 8 jets, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Acceptance and efficiency for the or the electroweakino production model in the EWK analysis discovery SR for the $2\ell^{\mathrm{sc}}$ category, considering $\tilde{\chi}^0_1 \tilde{\chi}^0_2$ production.

Cut flow for the gluino $\tilde{g}\rightarrow q\bar{q}\tilde{\chi}^0_1\rightarrow q\bar{q}q\bar{q} \ell\nu$ model with $m_{\tilde{g}} = 2$ TeV and $m_{\tilde{\chi}^0_1} = 1$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the gluino $\tilde{g}\rightarrow \bar{t}\tilde{t} \rightarrow \bar{t}\bar{b}\bar{s}$ model with with $m_{\tilde{g}} = 1.6$ TeV and $m_{\tilde{t}} = 1$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the gluino $\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1 \rightarrow t\bar{t} tbs$ model with $m_{\tilde{g}} = 2.2$ TeV and $m_{\tilde{\chi}^0_1} = 1.05$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the stop model with $m_{\tilde{t}} = 1.175$ TeV and $m_{\tilde{\chi}^0_1} = 0.7$ TeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the electroweakino production model, considering only the production of $\tilde{\chi}^{\pm}_1 \tilde{\chi}^0_1$, with $m(\tilde{\chi}^{\pm}_1,\tilde{\chi}^0_1)= 250$ GeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. In the $2\ell^{\mathrm{sc}}$ category no events are expected, as only one lepton is expected to be produced in the decay.Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

Cut flow for the electroweakino production model, considering only the production of $\tilde{\chi}^0_1 \tilde{\chi}^0_2$, with $m(\tilde{\chi}^0_1,\tilde{\chi}^0_2)= 250$ GeV. The column labelled $\mathrm{N}_{\mathrm{raw}}$ shows the number of generated events, while $\mathrm{N}_{\mathrm{events}}$ shows the expected number of events with a luminosity of 139fb$^{−1}$. The last column shows the cut flow efficiency with respect to all weighted events. The events are skimmed by requiring at least one electron or muon that satisfies very loose identification criteria, where the lepton satisfies $p_{\mathrm{T}} > 25$ GeV. The efficiency of this cut is considered in the quoted efficiency of the lepton trigger requirement. Selections with negligible inefficiencies on the given sample, such as data quality requirements, are not displayed. Selections that have not been evaluated in the analysis or are not applicable are denoted with a dash (--).

The centrality dependencies of the $\Delta\gamma\{\psi_\mathrm{ZDC}\}$ for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the A for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the a for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the A/a using full-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the A/a using sub-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $f_\mathrm{CME}$ using full-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $f_\mathrm{CME}$ using sub-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $\Delta\gamma_\mathrm{CME}$ using full-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The centrality dependencies of the $\Delta\gamma_\mathrm{CME}$ using sub-event method for Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $f_\mathrm{CME}$ from different methods for 20-50% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $f_\mathrm{CME}$ from different methods for 50-80% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $\Delta\gamma_\mathrm{CME}$ from different methods for 20-50% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The $\Delta\gamma_\mathrm{CME}$ from different methods for 50-80% centrality Au+Au collision at $\sqrt{s_{\rm NN}}$=200 GeV.

The contribution from different systematic uncertainties to the measured $t\bar{t}t\bar{t}$ production cross section, grouped in categories.

The results of the fitted signal strength $\mu$ in the 1L/2LOS and 2LSS/3L combined channel.

The results of fitted inclusive ${t\bar{t}t\bar{t}}$ cross-section in the 1L/2LOS and 2LSS/3L combined channel.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 1L,$\geq$9j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 1L,$\geq$9j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 2LOS,$\geq$7j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 2LOS,$\geq$7j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of number of jets in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of number of jets in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of number of jets in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of number of jets in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,$\geq$5b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,$\geq$5b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bL signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bL signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bH signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bH signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,$\geq$5b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,$\geq$5b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,$\geq$4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,$\geq$4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bL signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bL signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bH signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bH signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,$\geq$4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,$\geq$4b signal region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,4b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,4b control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,$\geq$5b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,$\geq$5b control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,$\geq$4b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,$\geq$4b control region after the fit.

95% CL exclusion limits on $\mathcal{B}(H\rightarrow aa \rightarrow b\bar{b}b\bar{b})$ for $m_a = 45$ GeV.

95% CL exclusion limits on $\mathcal{B}(H\rightarrow aa \rightarrow b\bar{b}b\bar{b})$ for $m_a = 55$ GeV.

The fraction of $a$ boson decays matched to reconstructed displaced vertices passing all vertex selections in signal MC.

The extrapolated signal selection efficiency as a function of $c\tau_{a}$.

$v_3$ vs $p_T$, p+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination

$v_3$ vs $p_T$, d+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination

$v_3$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, BBCS-FVTXS-CNT detector combination

$v_2$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_2$ vs $p_T$, d+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_2$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_2$ ratio vs $p_T$, p/d/3He+Au at 200 GeV, 0-5% central

$v_3$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_3$ vs $p_T$, d+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$v_3$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXS-CNT-FVTXN detector combination

$c_2$ vs $\eta$, p+Au at 200 GeV, 0-5% central

$c_2$ vs $\eta$, d+Au at 200 GeV, 0-5% central

$c_2$ vs $\eta$, 3He+Au at 200 GeV, 0-5% central

$c_3$ vs $\eta$, p+Au at 200 GeV, 0-5% central

$c_3$ vs $\eta$, d+Au at 200 GeV, 0-5% central

$c_3$ vs $\eta$, 3He+Au at 200 GeV, 0-5% central

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, BBCS-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, BBCS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, FVTXS-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, CNT-BBCS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, CNT-FVTXS correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+p at 200 GeV, minimum bias, CNT-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, p+Au at 200 GeV, 0-5% central, CNT-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, d+Au at 200 GeV, 0-5% central, CNT-FVTXN correlation

$C(\Delta\phi)$ vs $\Delta\phi$, 3He+Au at 200 GeV, 0-5% central, CNT-FVTXN correlation

$c_n$ vs detector combination, p+p at 200 GeV, Minimum Bias

$c_n$ vs $p_T$, p+p at 200 GeV, Minimum Bias, BBCS-CNT detector combination

$c_n$ vs $p_T$, p+p at 200 GeV, Minimum Bias, FVTXS-CNT detector combination

$c_n$ vs $p_T$, p+p at 200 GeV, Minimum Bias, FVTXN-CNT detector combination

$c_n$ vs detector combination, p+Au at 200 GeV, 0-5% central

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, BBCS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXN-CNT detector combination

$c_n$ vs detector combination, p+Au at 200 GeV, 0-5% central

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, BBCS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXS-CNT detector combination

$c_n$ vs $p_T$, p+Au at 200 GeV, 0-5% central, FVTXN-CNT detector combination

$c_n$ vs detector combination, p+Au at 200 GeV, 0-5% central

$c_n$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, BBCS-CNT detector combination

$c_n$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXS-CNT detector combination

$c_n$ vs $p_T$, 3He+Au at 200 GeV, 0-5% central, FVTXN-CNT detector combination

Differential cross section as a function of $m_{\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $p_{T,\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $a_{T,\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $\phi_{\eta}*$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $\pi-\Delta\phi_{\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $|cos\theta*|^{(CS)}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

The centrality dependence of Global $\Lambda$-hyperon Polarization at $\sqrt{s_{\rm NN}}=3$ GeV.

The $p_{\rm T}$ dependence of Global $\Lambda$-hyperon Polarization at $\sqrt{s_{\rm NN}}=3$ GeV.

The rapidity dependence of Global $\Lambda$-hyperon Polarization at $\sqrt{s_{\rm NN}}=3$ GeV.

$K^-$ (a), invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\phi$ meson (b) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\Xi^-$ (c) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 10--40\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$K^-$ (a), invariant yields as a function of $m_T-m_0$ for various rapidity regions in 40--60\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

$\phi$ meson (b) invariant yields as a function of $m_T-m_0$ for various rapidity regions in 40--60\% central Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. Statistics and systematic uncertainties are added quadratic here for plotting. Solid and dashed black lines depict $m_T$ exponential function fits to the measured data points with arbitrate scaling factors in each rapidity windows.

Rapidity density distributions of $K^-$ (squares), $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

Rapidity density distributions of $\phi$ meson (circles) $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

Rapidity density distributions of $\Xi^-$ (diamonds) $p_T$-integrated yields $dN/dy$ in 0--10\% (a), 10--40\% (b) and 40--60\% central (c) Au+Au collisions at ${\sqrt{s_{\mathrm{NN}}} = \mathrm{3\,GeV}}$. %The full symbols show the measured data, while the open ones are data reflected with respect to $y=0$ in the center-of-mass frame. Solid lines depict Gaussian function fits to the data points.

$\phi/K^-$ (a) and $\phi/\Xi^-$ (b) ratio as a function of collision energy, $\sqrt{s_{\mathrm{NN}}}$. The solid black circles show the measurements presented here in 0-10\% centrality bin, while empty markers in black or grey are used for data from various other energies and/or collision systems. The grey solid line represents a THERMUS calculation based on the Grand Canonical Ensemble (GCE) while the dotted lines depict calculations based on the Canonical Ensemble (CE) with different parameters of strangeness correlation radius ($r_c$). The green dashed line, green shaded band and the solid red line show transport model calculations from the public versions $\mathrm{UrQMD}^{1}$, modified $\mathrm{UrQMD}^{2}$ and SMASH, respectively.

Rapidity($y$) dependence of $v_1$ (top panels) and $v_2$ (bottom panels) of proton and $\Lambda$ baryons (left panels), pions (middle panels) and kaons (right panels) in 10-40% centrality for the $\sqrt{s_{NN}}$ = 3GeV Au+Au collisions. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The UrQMD and JAM results are shown as bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. The value of the incompressibility $\kappa$ = 380 MeV is used in the mean-field option. More detailed model descriptions and data comparisons can be found in Supplemental Material.

Rapidity($y$) dependence of $v_1$ (top panels) and $v_2$ (bottom panels) of proton and $\Lambda$ baryons (left panels), pions (middle panels) and kaons (right panels) in 10-40% centrality for the $\sqrt{s_{NN}}$ = 3GeV Au+Au collisions. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The UrQMD and JAM results are shown as bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. The value of the incompressibility $\kappa$ = 380 MeV is used in the mean-field option. More detailed model descriptions and data comparisons can be found in Supplemental Material.

Rapidity($y$) dependence of $v_1$ (top panels) and $v_2$ (bottom panels) of proton and $\Lambda$ baryons (left panels), pions (middle panels) and kaons (right panels) in 10-40% centrality for the $\sqrt{s_{NN}}$ = 3GeV Au+Au collisions. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The UrQMD and JAM results are shown as bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. The value of the incompressibility $\kappa$ = 380 MeV is used in the mean-field option. More detailed model descriptions and data comparisons can be found in Supplemental Material.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

$v_2$ scaled by the number of constituent quarks, $v_2/n_q$ , as a function of scaled transverse kinetic energy ($(m_T − m_0)/n_q$) for pions, kaons and protons from Au+Au collisions in 10-40% centrality at $\sqrt{s_{NN}}$ = 3, 27, and 54.4 GeV for positive charged particles (left panel) and negative charged particles (right panel). Colored dashed lines represent the scaling fit to data in 7.7, 14.5, 27, 54.4, and 200 GeV Au+Au collisions from STAR experiment at RHIC [43–45]. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Collision energy dependence of (top panel) directed flow slope $dv_1/dy|_{y=0}$ for p, $\Lambda, \pi$(combined from $\pi^{\pm}$), K(combined from $K^{\pm}$ and $K_S^0$), $\phi$ and $\Xi^−$, and (bottom panel) elliptic flow $v_2$ for p, $\pi$(combined from $\pi^{\pm}$) in heavy- ion collisions. Collision centrality for all data from RHIC is 10-40%, except for 4.5 GeV where 10-30% is for $dv_1/dy|_{y=0}$ and 0-30% is for $v_2$. Statistical and systematic uncertainties are shown as bars and gray bands, respectively. Some uncertainties are smaller than the data points. The JAM and UrQMD results are shown as colored bands:golden, red and blue bands stand for JAM mean-field, UrQMD mean-field and UrQMD cascade mode, respectively. For clarity the x-axis value of the data points have been shifted.

Distributions of $s_{\text{T}}$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $m_{\text{T}}(\tau)$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $\Sigma m_{\text{T}}(b_{1,2})$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $p_{\text{T}}(\tau)$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Observed event yields in data ('Observed') and expected event yields for SM background processes obtained from the background-only fit ('Total bkg.' and rows below) in the signal regions of the di-tau and single-tau channels. The quoted uncertainties include both the statistical and systematic uncertainties and are truncated at zero yield. By construction, no $t\bar{t}$ (2 real $\tau$) events can pass the selections in the single-tau channel. As the individual uncertainties are correlated, they do not add in quadrature to equal the total background uncertainty.

From left to right: upper limits at the 95% confidence level (CL) on the visible cross section ($\sigma_\text{vis}$) and on the number of signal events ($S_{\text{obs}}^{95}$). The third column ($S_{\text{exp}}^{95}$) shows the upper limit at the 95% CL on the number of signal events, given the expected number (and $\pm 1\,\sigma$ excursions on the expectation) of background events. The last two columns indicate the confidence level observed for the background-only hypothesis ($\text{CL}_{b}$), the discovery $p$-value ($p(s=0)$) and the significance ($Z$). In the di-tau SR, where fewer events are observed than predicted by the fitted background estimate, the $p$-value is capped at 0.5.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for the minimal-coupling scenario. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for the minimal-coupling scenario. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for vector leptoquarks with additional gauge couplings. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for vector leptoquarks with additional gauge couplings. The limits are derived from the binned single-tau signal region.

Exclusion contours at the 95% confidence level for the stop-stau signal model as a function of the masses of the top squark $m(\tilde{t}_{1})$ and of the tau slepton $m(\tilde{\tau}_{1})$. Expected and observed limits are shown for the present search in comparison to observed limits from previous ATLAS analyses based on data from Run-1 of the LHC at $\sqrt{s} = 8$ TeV [Eur. Phys. J. C 76 (2016)] and on a partial dataset from Run 2 at $\sqrt{s} = 13$ TeV [Phys. Rev. D 98 (2018) 032008]. The green band indicates the limit on the mass of the tau slepton (for a massless LSP) from the LEP experiments.

Exclusion contours at the 95% confidence level for the stop-stau signal model as a function of the masses of the top squark $m(\tilde{t}_{1})$ and of the tau slepton $m(\tilde{\tau}_{1})$. Expected and observed limits are shown for the present search in comparison to observed limits from previous ATLAS analyses based on data from Run-1 of the LHC at $\sqrt{s} = 8$ TeV [Eur. Phys. J. C 76 (2016)] and on a partial dataset from Run 2 at $\sqrt{s} = 13$ TeV [Phys. Rev. D 98 (2018) 032008]. The green band indicates the limit on the mass of the tau slepton (for a massless LSP) from the LEP experiments.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{u}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for up-type leptoquarks $\text{LQ}_{3}^{\text{u}})$ with charge $+2/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{u}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for up-type leptoquarks $\text{LQ}_{3}^{\text{u}})$ with charge $+2/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{d}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for down-type leptoquarks $\text{LQ}_{3}^{\text{d}})$ with charge $-1/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{d}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for down-type leptoquarks $\text{LQ}_{3}^{\text{d}})$ with charge $-1/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Upper limits on the signal cross section at the 95 % confidence level for the stop-stau signal model.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with up-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with down-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with minimal coupling (MC).

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with additional gauge couplings (YM).

Acceptance of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the di-tau SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau one-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau multi-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.