Signal-selection efficiency as a function of the di-tau invariant mass squared $q^2$ for simulated events in the SR. The error bars indicate the statistical uncertainty.

The unfolded dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The unfolded dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at backward and midrapidity regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with both jets at the midrapidity regions.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and forward regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and backward regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at backward and midrapidity regions, respectively.

Mean xj ratio, $\langle x_{j} \rangle$, in the range $10-60$ as unfolded for different multiplicities, for different leading and subleading jet combinations for the midrapidity, forward, and backward regions

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with both jets at the midrapidity regions.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at midrapidity and forward regions, respectively.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at midrapidity and backward regions, respectively.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at backward and midrapidity regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with both jets at the midrapidity regions.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and forward regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and backward regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at backward and midrapidity regions, respectively.

Mean xj ratio, $\langle x_{j} \rangle$, in the range $10-60$ as reconstructed (raw) for different multiplicities, for different leading and subleading jet combinations for the midrapidity, forward, and backward regions

$p_{T}$ dependence of $v_{2}$ for $\pi^{+}$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{+}$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{-}$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{-}$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{-}$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{-}$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $\pi^{-}$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{+}$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{+}$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{+}$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{+}$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{+}$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{-}$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{-}$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{-}$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{-}$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{-}$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{0}_{S}$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{0}_{S}$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{0}_{S}$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{0}_{S}$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $K^{0}_{S}$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $p$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $p$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $p$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $p$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $p$ in Au+Au collisions at 4.5 GeV

$p_{T}$ dependence of $v_{2}$ for $\Lambda$ in Au+Au collisions at 3 GeV

$p_{T}$ dependence of $v_{2}$ for $\Lambda$ in Au+Au collisions at 3.2 GeV

$p_{T}$ dependence of $v_{2}$ for $\Lambda$ in Au+Au collisions at 3.5 GeV

$p_{T}$ dependence of $v_{2}$ for $\Lambda$ in Au+Au collisions at 3.9 GeV

$p_{T}$ dependence of $v_{2}$ for $\Lambda$ in Au+Au collisions at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{+}$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{+}$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{+}$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{+}$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{+}$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{-}$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{-}$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{-}$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{-}$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\pi^{-}$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{+}$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{+}$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{+}$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{+}$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{+}$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{-}$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{-}$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{-}$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{-}$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{-}$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{0}_{S}$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{0}_{S}$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{0}_{S}$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{0}_{S}$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $K^{0}_{S}$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $p$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $p$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $p$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $p$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $p$ at 4.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\Lambda$ at 3 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\Lambda$ at 3.2 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\Lambda$ at 3.5 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\Lambda$ at 3.9 GeV

The number of constituent quarks $n_q$ scaled $v_2$ as a function of $n_q$ scaled transverse kinetic energy for $\Lambda$ at 4.5 GeV

The energy dependence of $p_T$ integrated $v_2$ for $\pi^{\pm}$, $K^{\pm}$, $K^{0}_{S}$, $p$, and $\Lambda$ in 10-40\% centrality from Au+Au collisions at $\sqrt{s_{\text{NN}}}$ = 3.0, 3.2, 3.5, 3.9, and 4.5 GeV.

The energy dependence of $n_{q}$ scaled $v_{2}$ ratios of $v_{2}^{q}(\pi^{+}) / v_{2}^{q}(K^{+})$ and $v_{2}^{q}(p) / v_{2}^{q}(K^{+})$ at ($m_T - m_0$)/$n_q$ = 0.4 GeV/$c^2$

UrQMD model calculations for net-proton cumulant ratios and proton factorial cumulant ratios in Au+Au collisions from $\sqrt{s_{NN}}$ = 7.7 - 200 GeV for 0-5$\%$ centrality. Same choices of centrality definition are used as used for experimental data.

Significance of deviations of central 0-5$\%$ data from various choices of baselines or references.

Cumulant ratios C2/C1 and C4/C2 of net-proton distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV as a function of multiplicity (RefMult3X).

Cumulant ratios C2/C1 and C4/C2 of net-proton distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV as a function of average multiplicity (RefMult3X) for 8 centrality bins from 0-10$\%$ to 70-80$\%$. Data with and without CBWC are given.

Cumulant ratios C2/C1 and C4/C2 of net-proton distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV as a function of average multiplicity (RefMult3X) for 16 centrality bins from 0-5$\%$ to 75-80$\%$. Data with and without CBWC are given.

Reference multiplicity distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 11.5 GeV from BES-II. Three choices of multiplicity distributions (RefMult3 scaled by efficiency factor of $\epsilon$ = 0.83 or 83$\%$ to mimic RefMult3 of BES-I, RefMult3, and RefMult3X) are presented.

Cumulant ratios C2/C1 and C4/C2 of net-proton distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 11.5 GeV from BES-II as a function of centrality obtained with various centrality definitions.

Reference multiplicity distributions (RefMult3, RefMult3X, and RefMult3A) in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV from UrQMD model. The RefMult3 and RefMult3X distributions have been scaled by an efficiency factor of $\epsilon$ = 77.9$\%$ to mimic the distribution in experimenal data.

Cumulant ratios C2/C1 and C4/C2 of net-proton distributions in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV from UrQMD model as a function of centrality obtained with various centrality definitions (RefMult3, RefMult3X, and RefMult3A). The RefMult3 and RefMult3X distributions have been scaled by an efficiency factor of $\epsilon$ = 77.9$\%$ to mimic the distribution in experimenal data.

Collision energy dependence of net-proton C4/C2 in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7 - 27 GeV from BES-II obtained with RefMult3 centrality definition.

Difference between net-proton C4/C2 (0-5$\%$) data from various choices of baselines or references as a function of collision energy from $\sqrt{s_{NN}}$ = 7.7 - 200 GeV.

The distributions of the variables used in the simultaneous fit for the SR. The black points with error bars represent the data and their statistical uncertainties, whereas the shaded band represents the predicted uncertainties. The bottom panel in each figure shows the ratio of the number of events observed in data to that of the total SM prediction. The last bin of each plot has been extended to include the overflow contribution.

Expected and observed 95% upper limits on the product of the cross section and branching fraction$\sigma(\mathrm{pp} \ ightarrow W\mathrm{a})\mathcal{B}(W ightarrow \ell^{+} v_{\ell})\mathcal{B}(\mathrm{a} \rightarrow Z\gamma)\mathcal{B}(Z \rightarrow \ell^{+} \ell^{-})$ s a function of the ALP mass from 110 to 400 GeV

Expected and observed 95% upper limits on the photophobic ALP model parameter $1/f_{\mathrm{a}}$ as a function of ALP mass reinterpreted from $1/f_{\mathrm{a}}=2~\mathrm{TeV}^{-1}$

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the SR (t-pass) and SR (t-fail) after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. A representative mono-top signal (vector coupling scenario) with a mediator mass of 1 TeV, a DM candidate mass of 150 GeV, and a cross section of 1 pb is overlaid as an orange line. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the SR (t-pass) and SR (t-fail) after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. A representative mono-top signal (vector coupling scenario) with a mediator mass of 1 TeV, a DM candidate mass of 150 GeV, and a cross section of 1 pb is overlaid as an orange line. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W(e$\nu$) (t-pass) and W(e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W(e$\nu$) (t-pass) and W(e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W($\mu\nu$) (t-pass) and W($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W($\mu\nu$) (t-pass) and W($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z(ee) (t-pass) and Z(ee) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z(ee) (t-pass) and Z(ee) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z($\mu\mu$) (t-pass) and Z($\mu\mu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z($\mu\mu$) (t-pass) and Z($\mu\mu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the t$\bar{t}$ (e$\nu$) (t-pass) and t$\bar{t}$ (e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the t$\bar{t}$ (e$\nu$) (t-pass) and t$\bar{t}$ (e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the tt ($\mu\nu$) (t-pass) and tt ($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the tt ($\mu\nu$) (t-pass) and tt ($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the $\gamma$ (t-pass) and $\gamma$ (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the $\gamma$ (t-pass) and $\gamma$ (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Upper limits at 95% CL on $\sigma\mathcal{B}$ of mono-top production presented in the two- dimensional plane spanned by the mediator and DM candidate masses for a mediator mass between 200 and 2250 GeV and a DM candidate mass between 1 and 1125 GeV only consider- ing on-shell decays of the mediator to the DM candidates. The mediator has vector couplings to quarks and DM candidates in the upper plot and axial-vector couplings in the lower plot. The median expected exclusion range is indicated by a black solid line, demonstrating the search sensitivity of the analysis. The 68% probability interval of the expected exclusion is shown in black dashed lines. Contours of theory predictions for constant values of $\sigma\mathcal{B}$ are shown in gray dashed lines. The observed exclusion contour of mediator and DM candidate masses is rep- resented by the red solid line. The exclusion contour obtained from measurements of the DM relic density $\Omega_{\mathrm{nbm}}h^2$ by the Planck Collaboration is shown in the gray solid line.

Upper limits at 95% CL on $\sigma\mathcal{B}$ of mono-top production presented in the two- dimensional plane spanned by the mediator and DM candidate masses for a mediator mass between 200 and 2250 GeV and a DM candidate mass between 1 and 1125 GeV only consider- ing on-shell decays of the mediator to the DM candidates. The mediator has vector couplings to quarks and DM candidates in the upper plot and axial-vector couplings in the lower plot. The median expected exclusion range is indicated by a black solid line, demonstrating the search sensitivity of the analysis. The 68% probability interval of the expected exclusion is shown in black dashed lines. Contours of theory predictions for constant values of $\sigma\mathcal{B}$ are shown in gray dashed lines. The observed exclusion contour of mediator and DM candidate masses is rep- resented by the red solid line. The exclusion contour obtained from measurements of the DM relic density $\Omega_{\mathrm{nbm}}h^2$ by the Planck Collaboration is shown in the gray solid line.