Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-ej and (c, d) SR-2L-ej of the e+light-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_de) = (2.0 TeV, 1.0) and S̃₁ (m, y_de) = (3.0 TeV, 1.0), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-μj and (c, d) SR-2L-μj of the μ+light-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_sμ) = (2.0 TeV, 1.5) and S̃₁ (m, y_sμ) = (3.0 TeV, 1.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-μj and (c, d) SR-2L-μj of the μ+light-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_sμ) = (2.0 TeV, 1.5) and S̃₁ (m, y_sμ) = (3.0 TeV, 1.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-μj and (c, d) SR-2L-μj of the μ+light-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_sμ) = (2.0 TeV, 1.5) and S̃₁ (m, y_sμ) = (3.0 TeV, 1.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-μj and (c, d) SR-2L-μj of the μ+light-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_sμ) = (2.0 TeV, 1.5) and S̃₁ (m, y_sμ) = (3.0 TeV, 1.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-eb and (c, d) SR-2L-eb of the e+b-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_be) = (1.5 TeV, 2.5) and S̃₁ (m, y_be) = (2.0 TeV, 2.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-eb and (c, d) SR-2L-eb of the e+b-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_be) = (1.5 TeV, 2.5) and S̃₁ (m, y_be) = (2.0 TeV, 2.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-eb and (c, d) SR-2L-eb of the e+b-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_be) = (1.5 TeV, 2.5) and S̃₁ (m, y_be) = (2.0 TeV, 2.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-eb and (c, d) SR-2L-eb of the e+b-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_be) = (1.5 TeV, 2.5) and S̃₁ (m, y_be) = (2.0 TeV, 2.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-μb and (c, d) SR-2L-μb of the μ+b-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_bμ) = (1.5 TeV, 1.5) and S̃₁ (m, y_bμ) = (2.0 TeV, 1.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-μb and (c, d) SR-2L-μb of the μ+b-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_bμ) = (1.5 TeV, 1.5) and S̃₁ (m, y_bμ) = (2.0 TeV, 1.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-μb and (c, d) SR-2L-μb of the μ+b-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_bμ) = (1.5 TeV, 1.5) and S̃₁ (m, y_bμ) = (2.0 TeV, 1.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Data (dots) and post-fit SM distribution (histograms) of m_ℓj in (a, b) SR-1L-μb and (c, d) SR-2L-μb of the μ+b-jet channel obtained by a CR+SR background-only fit for Run 2 and Run 3, respectively. The lower panel shows the ratio of observed data to the total post- and pre-fit SM prediction. The last bin includes the overflow. Uncertainties in the background estimates include both the statistical and systematic uncertainties, with correlations between uncertainties taken into account. The dashed lines show the predicted yields for two benchmark signal models corresponding to S̃₁ (m, y_bμ) = (1.5 TeV, 1.5) and S̃₁ (m, y_bμ) = (2.0 TeV, 1.5), respectively. Note: the values in the table are normalized by the width of corresponding bin

Exclusion limits for minimal models of S̃₁ production with only y_de being non-zero obtained from a simultaneous fit to SR-1L and SR-2L of the e+light-jet channel combining Run-2 and Run-3 data. All limits are computed at 95% CL and the observed (red solid lines) and expected (black dashed lines) exclusion limits are shown in the m(S̃₁) – y_de plane. The observed exclusion should be interpreted as the region above the red line. The yellow inner (green outer) shaded band around the expected limits corresponds to the ±1 σ (±2 σ) variations of the expected limit, accounting for all uncertainties. The observed limit obtained from ATLAS searches for LQ pair production is also shown as dark blue line [13]. Constraints from weak charge measurements of protons and nuclei on y_de couplings derived by Ref. [10] are shown as light magenta line.

Exclusion limits for minimal models of S̃₁ production with only y_sμ being non-zero obtained from a simultaneous fit to SR-1L and SR-2L of the μ+light-jet channel combining Run-2 and Run-3 data. All limits are computed at 95% CL and the observed (red solid lines) and expected (black dashed lines) exclusion limits are shown in the m(S̃₁) – y_sμ plane. The observed exclusion should be interpreted as the region above the red line. The yellow inner (green outer) shaded band around the expected limits corresponds to the ±1 σ (±2 σ) variations of the expected limit, accounting for all uncertainties. The observed limit obtained from ATLAS searches for LQ pair production is also shown as dark blue line [13].

Exclusion limits for minimal models of S̃₁ production with only y_be being non-zero obtained from a simultaneous fit to SR-1L and SR-2L of the e+b-jet channel combining Run-2 and Run-3 data. All limits are computed at 95% CL and the observed (red solid lines) and expected (black dashed lines) exclusion limits are shown in the m(S̃₁) – y_be plane. The observed exclusion should be interpreted as the region above the red line. The yellow inner (green outer) shaded band around the expected limits corresponds to the ±1 σ (±2 σ) variations of the expected limit, accounting for all uncertainties. The observed limit obtained from ATLAS searches for LQ pair production is also shown as dark blue line [13].

Exclusion limits for minimal models of S̃₁ production with only y_bμ being non-zero obtained from a simultaneous fit to SR-1L and SR-2L of the μ+b-jet channel combining Run-2 and Run-3 data. All limits are computed at 95% CL and the observed (red solid lines) and expected (black dashed lines) exclusion limits are shown in the m(S̃₁) – y_bμ plane. The observed exclusion should be interpreted as the region above the red line. The yellow inner (green outer) shaded band around the expected limits corresponds to the ±1 σ (±2 σ) variations of the expected limit, accounting for all uncertainties. The observed limit obtained from ATLAS searches for LQ pair production is also shown as dark blue line [13].

Exclusion limits for minimal models of S̃₁ production obtained from a fit to SR-1L, SR-2L and their combination of (a) the e+light-jet, (b) the μ+light-jet, (c) the e+b-jet and (d) the μ+b-jet channel using Run-2 and Run-3 data. All limits are computed at 95% CL and the observed (solid lines) and expected (dashed lines) exclusion limits are shown in the m(S̃₁) – y plane. The observed exclusion regions should be interpreted as the region above the solid lines. Constraints from weak charge measurements of protons and nuclei  and ATLAS searches for LQ pair production  are shown as light magenta and dark blue areas, respectively.

Exclusion limits for minimal models of S̃₁ production obtained from a fit to SR-1L, SR-2L and their combination of (a) the e+light-jet, (b) the μ+light-jet, (c) the e+b-jet and (d) the μ+b-jet channel using Run-2 and Run-3 data. All limits are computed at 95% CL and the observed (solid lines) and expected (dashed lines) exclusion limits are shown in the m(S̃₁) – y plane. The observed exclusion regions should be interpreted as the region above the solid lines. Constraints from weak charge measurements of protons and nuclei  and ATLAS searches for LQ pair production  are shown as light magenta and dark blue areas, respectively.

Exclusion limits for minimal models of S̃₁ production obtained from a fit to SR-1L, SR-2L and their combination of (a) the e+light-jet, (b) the μ+light-jet, (c) the e+b-jet and (d) the μ+b-jet channel using Run-2 and Run-3 data. All limits are computed at 95% CL and the observed (solid lines) and expected (dashed lines) exclusion limits are shown in the m(S̃₁) – y plane. The observed exclusion regions should be interpreted as the region above the solid lines. Constraints from weak charge measurements of protons and nuclei  and ATLAS searches for LQ pair production  are shown as light magenta and dark blue areas, respectively.

Exclusion limits for minimal models of S̃₁ production obtained from a fit to SR-1L, SR-2L and their combination of (a) the e+light-jet, (b) the μ+light-jet, (c) the e+b-jet and (d) the μ+b-jet channel using Run-2 and Run-3 data. All limits are computed at 95% CL and the observed (solid lines) and expected (dashed lines) exclusion limits are shown in the m(S̃₁) – y plane. The observed exclusion regions should be interpreted as the region above the solid lines. Constraints from weak charge measurements of protons and nuclei  and ATLAS searches for LQ pair production  are shown as light magenta and dark blue areas, respectively.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$ej$ and (b) Run 2 SR-2L-$ej$ and (c) Run 3 SR-1-$ej$ and (d) Run 3 SR-2-$ej$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$ej$ and (b) Run 2 SR-2L-$ej$ and (c) Run 3 SR-1-$ej$ and (d) Run 3 SR-2-$ej$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$ej$ and (b) Run 2 SR-2L-$ej$ and (c) Run 3 SR-1-$ej$ and (d) Run 3 SR-2-$ej$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$ej$ and (b) Run 2 SR-2L-$ej$ and (c) Run 3 SR-1-$ej$ and (d) Run 3 SR-2-$ej$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$ej$ and (b) Run 2 SR-2L-$ej$ and (c) Run 3 SR-1-$ej$ and (d) Run 3 SR-2-$ej$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$ej$ and (b) Run 2 SR-2L-$ej$ and (c) Run 3 SR-1-$ej$ and (d) Run 3 SR-2-$ej$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$ej$ and (b) Run 2 SR-2L-$ej$ and (c) Run 3 SR-1-$ej$ and (d) Run 3 SR-2-$ej$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$ej$ and (b) Run 2 SR-2L-$ej$ and (c) Run 3 SR-1-$ej$ and (d) Run 3 SR-2-$ej$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu j$ and (b) Run 2 SR-2L-$\mu j$ and (c) Run 3 SR-1L-$\mu j$ and (d) Run 3 SR-2L-$\mu j$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu j$ and (b) Run 2 SR-2L-$\mu j$ and (c) Run 3 SR-1L-$\mu j$ and (d) Run 3 SR-2L-$\mu j$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu j$ and (b) Run 2 SR-2L-$\mu j$ and (c) Run 3 SR-1L-$\mu j$ and (d) Run 3 SR-2L-$\mu j$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu j$ and (b) Run 2 SR-2L-$\mu j$ and (c) Run 3 SR-1L-$\mu j$ and (d) Run 3 SR-2L-$\mu j$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu j$ and (b) Run 2 SR-2L-$\mu j$ and (c) Run 3 SR-1L-$\mu j$ and (d) Run 3 SR-2L-$\mu j$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu j$ and (b) Run 2 SR-2L-$\mu j$ and (c) Run 3 SR-1L-$\mu j$ and (d) Run 3 SR-2L-$\mu j$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu j$ and (b) Run 2 SR-2L-$\mu j$ and (c) Run 3 SR-1L-$\mu j$ and (d) Run 3 SR-2L-$\mu j$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu j$ and (b) Run 2 SR-2L-$\mu j$ and (c) Run 3 SR-1L-$\mu j$ and (d) Run 3 SR-2L-$\mu j$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$eb$ and (b) Run 2 SR-2L-$eb$ and (c) Run 3 SR-1L-$eb$ and (d) Run 3 SR-2L-$eb$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$eb$ and (b) Run 2 SR-2L-$eb$ and (c) Run 3 SR-1L-$eb$ and (d) Run 3 SR-2L-$eb$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$eb$ and (b) Run 2 SR-2L-$eb$ and (c) Run 3 SR-1L-$eb$ and (d) Run 3 SR-2L-$eb$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$eb$ and (b) Run 2 SR-2L-$eb$ and (c) Run 3 SR-1L-$eb$ and (d) Run 3 SR-2L-$eb$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP of $\tilde{S}$ in (a) Run 2 SR-1L-$eb$ and (b) Run 2 SR-2L-$eb$ and (c) Run 3 SR-1L-$eb$ and (d) Run 3 SR-2L-$eb$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP of $\tilde{S}$ in (a) Run 2 SR-1L-$eb$ and (b) Run 2 SR-2L-$eb$ and (c) Run 3 SR-1L-$eb$ and (d) Run 3 SR-2L-$eb$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP of $\tilde{S}$ in (a) Run 2 SR-1L-$eb$ and (b) Run 2 SR-2L-$eb$ and (c) Run 3 SR-1L-$eb$ and (d) Run 3 SR-2L-$eb$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP of $\tilde{S}$ in (a) Run 2 SR-1L-$eb$ and (b) Run 2 SR-2L-$eb$ and (c) Run 3 SR-1L-$eb$ and (d) Run 3 SR-2L-$eb$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu b$ and (b) Run 2 SR-2L-$\mu b$ and (c) Run 3 SR-1L-$\mu b$ and (d) Run 3 SR-2L-$\mu b$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu b$ and (b) Run 2 SR-2L-$\mu b$ and (c) Run 3 SR-1L-$\mu b$ and (d) Run 3 SR-2L-$\mu b$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu b$ and (b) Run 2 SR-2L-$\mu b$ and (c) Run 3 SR-1L-$\mu b$ and (d) Run 3 SR-2L-$\mu b$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for resonant production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu b$ and (b) Run 2 SR-2L-$\mu b$ and (c) Run 3 SR-1L-$\mu b$ and (d) Run 3 SR-2L-$\mu b$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu b$ and (b) Run 2 SR-2L-$\mu b$ and (c) Run 3 SR-1L-$\mu b$ and (d) Run 3 SR-2L-$\mu b$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu b$ and (b) Run 2 SR-2L-$\mu b$ and (c) Run 3 SR-1L-$\mu b$ and (d) Run 3 SR-2L-$\mu b$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu b$ and (b) Run 2 SR-2L-$\mu b$ and (c) Run 3 SR-1L-$\mu b$ and (d) Run 3 SR-2L-$\mu b$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Values of the signal acceptance multiplied by the selection efficiency ($A\times\epsilon$) for DY+SP production of $\tilde{S}$ in (a) Run 2 SR-1L-$\mu b$ and (b) Run 2 SR-2L-$\mu b$ and (c) Run 3 SR-1L-$\mu b$ and (d) Run 3 SR-2L-$\mu b$. The $A \times \epsilon$ values are calculated as the number of events of reconstructed-level signal simulation divided by the total cross-section of the signal process times the integrated luminosity.

Event selection cutflows of SR-1L-$eb$ for the signal sample with $m(\tilde{S}) = 2$ TeV and $y_{be} = 1.0$ individually for resonant production (RP) and combined Drell-Yan plus single production (DY+SP). The first number in the column for each mass points corresponds to the signal event yield after applying the associated cut while the second number in parentheses represents the efficiency with respect to the previous cut.

Event selection cutflows of SR-2L-$eb$ for the signal sample with $m(\tilde{S}) = 2$ TeV and $y_{be} = 1.0$ individually for resonant production (RP) and combined Drell-Yan plus single production (DY+SP). The first number in the column for each mass points corresponds to the signal event yield after applying the associated cut while the second number in parentheses represents the efficiency with respect to the previous cut.

Event selection cutflows of SR-1L-$ej$ for the signal sample with $m(\tilde{S}) = 2$ TeV and $y_{de} = 1.0$ individually for resonant production (RP) and combined Drell-Yan plus single production (DY+SP). The first number in the column for each mass points corresponds to the signal event yield after applying the associated cut while the second number in parentheses represents the efficiency with respect to the previous cut.

Event selection cutflows of SR-2L-$ej$ for the signal sample with $m(\tilde{S}) = 2$ TeV and $y_{de} = 1.0$ individually for resonant production (RP) and combined Drell-Yan plus single production (DY+SP). The first number in the column for each mass points corresponds to the signal event yield after applying the associated cut while the second number in parentheses represents the efficiency with respect to the previous cut.

Event selection cutflows of SR-1L-$\mu b$ for the signal sample with $m(\tilde{S}) = 2$ TeV and $y_{b\mu} = 1.0$ individually for resonant production (RP) and combined Drell-Yan plus single production (DY+SP). The first number in the column for each mass points corresponds to the signal event yield after applying the associated cut while the second number in parentheses represents the efficiency with respect to the previous cut.

Event selection cutflows of SR-2L-$\mu b$ for the signal sample with $m(\tilde{S}) = 2$ TeV and $y_{b\mu} = 1.0$ individually for resonant production (RP) and combined Drell-Yan plus single production (DY+SP). The first number in the column for each mass points corresponds to the signal event yield after applying the associated cut while the second number in parentheses represents the efficiency with respect to the previous cut.

Event selection cutflows of SR-1L-$\mu j$ for the signal sample with $m(\tilde{S}) = 2$ TeV and $y_{s\mu} = 1.0$ individually for resonant production (RP) and combined Drell-Yan plus single production (DY+SP). The first number in the column for each mass points corresponds to the signal event yield after applying the associated cut while the second number in parentheses represents the efficiency with respect to the previous cut.

Event selection cutflows of SR-2L-$\mu j$ for the signal sample with $m(\tilde{S}) = 2$ TeV and $y_{s\mu} = 1.0$ individually for resonant production (RP) and combined Drell-Yan plus single production (DY+SP). The first number in the column for each mass points corresponds to the signal event yield after applying the associated cut while the second number in parentheses represents the efficiency with respect to the previous cut.

Observed 95% $\text{CL}_\text{S}$ upper limits on the production cross-section times acceptance times branching ratio to jets, $\sigma \cdot A \cdot \text{BR}$, of Gaussian-shaped signals of 5%, 10%, and 15% width relative to their peak mass, $m_G$. Also included are the corresponding expected upper limits predicted for the case the $m_{jj}$ distribution is observed to be identical to the background prediction in each bin and the $1\sigma$ and $2\sigma$ envelopes of outcomes expected for Poisson fluctuations around the background expectation. Limits are derived from the J100 signal region.

Observed 95% $\text{CL}_\text{S}$ upper limits on the the value of the coupling to quarks, $g_q$, of the $Z^\prime$ model as a function of the resonance mass $m_{Z^\prime}$. Also included are the corresponding expected upper limits predicted for the case the $m_{jj}$ distribution is observed to be identical to the background prediction in each bin and the $1\sigma$ and $2\sigma$ envelopes of outcomes expected for Poisson fluctuations around the background expectation. Limits are derived from the J50 signal region.

Observed 95% $\text{CL}_\text{S}$ upper limits on the the value of the coupling to quarks, $g_q$, of the $Z^\prime$ model as a function of the resonance mass $m_{Z^\prime}$. Also included are the corresponding expected upper limits predicted for the case the $m_{jj}$ distribution is observed to be identical to the background prediction in each bin and the $1\sigma$ and $2\sigma$ envelopes of outcomes expected for Poisson fluctuations around the background expectation. Limits are derived from the J100 signal region.

Acceptance of simulated $Z^\prime$ signal events given the analysis event selection as a function of $m_{Z^\prime}$. The solid black line corresponds to the acceptance without a lower $m_{jj}$ requirement applied, the solid blue line includes the $m_{jj} > 344$ GeV requirement of the J50 signal region, and the solid red line includes the $m_{jj} > 481$ GeV requirement of the J100 signal region.

Numbers of events in each of the two considered datasets after consecutively applying the analysis selections.