The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, 0 forward jet, $t\leq 0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, 0 forward jet, $t>0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, $\geq 1$ forward jet, $t\leq 0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, $\geq 1$ forward jet, $t>0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the $\geq 2$ b jet, $t\leq 0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the $\geq 2$ b jet, $t>0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the $\geq 2$ b jet, different flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the $\geq 2$ b jet, same flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the 1 b jet, different flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the 1 b jet, same flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The model-independent 95\% CL limits on production cross section for new physics processes for the scalar interaction. The expected limit is shown by the black dashed line with the 68 and 95\% CL uncertainty bands in green and yellow, respectively, while the observed limit is shown by the solid black line. Theoretical LO cross section values for the DM model and their associated uncertainties are also presented (grey line).

The model-independent 95\% CL limits on production cross section for new physics processes for the pseudoscalar interaction. The expected limit is shown by the black dashed line with the 68 and 95\% CL uncertainty bands in green and yellow, respectively, while the observed limit is shown by the solid black line. Theoretical LO cross section values for the DM model and their associated uncertainties are also presented (grey line).

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the QCD control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the tt(1l) control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the W+jets control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the Z(2l) control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the W+jets control region in the single lepton channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the tt(2l) control region in the single lepton channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $(N_{\text{jet}},N_{\text{bjet}})$ distribution of the ttZ control region in the dileptonic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit NN distribution of the Drell-Yan control region in the dileptonic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

Predicted signal yields for the 2016, 2017, and 2018 data-taking periods, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Cross sections calculated to leading order (LO) accuracy for the tt+DM scalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the tt+DM pseudoscalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the tt+DM scalar mediator signal process considering dileptonic decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the tt+DM pseudoscalar mediator signal process considering dileptonic decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM t-channel scalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel scalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM t-channel pseudoscalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel pseudoscalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel scalar mediator signal process considering dileptonic decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel pseudoscalar mediator signal process considering dileptonic decays from the visible particles in the final state.

Observed best fit differential cross section for the $|y_{H}|$ observable

Observed best fit differential cross section for the $m_{jj}$ (GeV) observable

Observed best fit differential cross section for the $\Delta\eta_{jj}$ observable.

Observed best fit differential cross section for the $\tau_{C}^{j}$ (GeV) observable

Bin-to-bin correlation matrix for the $p_{T}^{H}$ spectrum

Bin-to-bin correlation matrix for the $N_{jets}$ spectrum

Bin-to-bin correlation matrix for the $p_{T}^{j0}$ spectrum

Bin-to-bin correlation matrix for the $|y_{H}|$ spectrum

Bin-to-bin correlation matrix for the $m_{jj}$ spectrum

Bin-to-bin correlation matrix for the $\Delta\eta_{jj}$ spectrum

Bin-to-bin correlation matrix for the $\tau_{C}^{j}$ spectrum

Rotation matrix for the PCA of the SMEFT Wilson coefficients

Correlation matrix of the linear combinations of Wilson coefficients obtained from the PCA, obtained by fitting the observed data to the $p_{T}^{H}$ spectra of all decay channels.

Summary of observed and expected confidence intervals for the first eigenvectors obtained from the PCA of the SMEFT Wilson coefficients.

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

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}$

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator distributions constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions and for pp collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 1$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=1$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $120 < p_{\mathrm{T,jet}} < 140$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $140 < p_{\mathrm{T,jet}} < 160$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $160 < p_{\mathrm{T,jet}} < 180$ GeV. The results are shown for different centrality bins in PbPb collisions.

The energy-energy correlator PbPb/pp ratios constructed with charged particles with $p_{\mathrm{T}} > 2$ GeV for energy weight $n=2$ and jet $p_{\mathrm{T}}$ selection $180 < p_{\mathrm{T,jet}} < 200$ GeV. The results are shown for different centrality bins in PbPb collisions.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

The double ratios of PbPb to pp single ratios with $p_{\mathrm{T}}^{\mathrm{ch}}$ > 2 GeV and $p_{\mathrm{T}}^{\mathrm{ch}}$ > 1 GeV.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{e u}|$, versus $S_{e u}$ mass for scalar LQs coupled to electrons.

Observed and Expected UL exclusions on the $BR(H\to S)$ of generic signals with $m_{A'} = 1.0\;GeV$ and $BR(A' \rightarrow \pi\pi) = 1.0$.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $S_{e d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Observed and Expected UL exclusions on the number of 'signal' yields in sideband A from the model agnostic fit, assuming no 'SUEP-like' signal contamination in the CRs.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{e d}|$, versus $S_{e d}$ mass for scalar LQs coupled to electrons.

Postfit values for SR in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $S_{\mu u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Postfit values for PR in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{\mu u}|$, versus $S_{\mu u}$ mass for scalar LQs coupled to muons.

Postfit values for CRDY in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $S_{\mu d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Postfit values for CRTT in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{\mu d}|$, versus $S_{\mu d}$ mass for scalar LQs coupled to muons.

Per-selection efficiencies on Generic - $m_{A'} = 1.0\;GeV$ $BR(A' \rightarrow \pi\pi) = 1.0$ Hadronic - $m_{A'} = 0.7\;GeV$ and $BR(A' \rightarrow ee) = BR(A' \rightarrow \mu\mu) = 0.15$ $BR(A' \rightarrow \pi\pi) = 0.7$ Leptonic - $m_{A'} = 0.5\;GeV$ and $BR(A' \rightarrow ee) = BR(A' \rightarrow \mu\mu) = 0.2$ and $BR(A' \rightarrow \pi\pi) = 0.6$ samples. Selections: 2 OSSF leptons, one SUEP cluster, $60\leq Z_{m}\leq 120$ GeV, $p_{T}^Z\geq 25$ GeV, veto events with b-tagged jets, $p_{T}^{SUEP} > 60\;GeV$, and $dR^{Ak4}_{Ak15} < 1.5$. Efficiencies are the ratio of the histogram size at the selection step to the previous step (except twoLepton vs. initial). The final column is the total cumulative efficiency (dR vs. initial). Note that these selection numbers were all calculated using samples generated with 600,000 events for each signal point considered.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $V_{e u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Interpolated observed and expected UL exclusions on $BR(H\to S) = 1$ for generic signals with $m_{A'} = 1.0\;GeV$ and $BR(A' \rightarrow \pi\pi) = 1.0$. The discreate heatmap points can be found in the UL Generic (Observed & Expected) figure in this sandbox

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{e u}|$, versus $V_{e u}$ mass for vector LQs coupled to electrons.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $V_{e d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{e d}|$, versus $V_{e d}$ mass for vector LQs coupled to electrons.

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $V_{\mu u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{\mu u}|$, versus $V_{\mu u}$ mass for vector LQs coupled to muons.

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $V_{\mu d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{\mu d}|$, versus $V_{\mu d}$ mass for vector LQs coupled to muons.

Distributions of the number of hits in the cluster (Nhits) for the CSC category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the VLL mass for the pseudoscalar mean proper decay length c$\tau$ = 0.025 m. The pseudoscalar mass is 2 GeV.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the pseudoscalar mean proper decay length c$\tau$ for VLL Mass of 700 GeV. The pseudoscalar mass is 2 GeV.

The 95% CL observed upper limits on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

The 95% CL observed exclusion contour on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the CSC category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the DT category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.