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.

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

Distributions of the number of hits in the cluster (Nhits) for the DT category in the signal region (SR). The last histogram bin contains all overflow events.

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

Distributions of the number of hits in the cluster (Nhits) for the CSC category in the signal region (SR). The last histogram bin contains all overflow events.

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.

Distributions of the number of hits in the cluster (Nhits) for the DT 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.

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 pseudoscalar mean proper decay length c$\tau$ for VLL Mass of 700 GeV. 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 VLL mass for the pseudoscalar mean proper decay length c$\tau$ = 0.025 m. 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 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 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.

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.

Selection efficiencies in the CSC category signal region (SR) for different signal hypothesis. 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 DT category signal region (SR) for different signal hypothesis. 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.

The cluster reconstruction efficiency, including both DT and CSC clusters, as a function of the simulated r and |z| decay positions of the pseudoscalar into photons in events with MET > 200 GeV, for a VLL mass of 300 GeV and a pseudoscalar mass of 2 GeV, and a range of ctau values uniformly distributed between 0.01 and 0.1 m.

The cluster reconstruction efficiency, including both DT and CSC clusters, as a function of the simulated r and |z| decay positions of the pseudoscalar into photons in events with MET > 200 GeV, for a VLL mass of 500 GeV and a pseudoscalar mass of 2 GeV, and a range of ctau values uniformly distributed between 0.01 and 0.1 m.

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $6.5<p_T<30.0\,GeV$ and $ 1 < y_{CM} < 1.935$ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $6.5<p_T<30.0\,GeV$ and $ -2.865 < y_{CM} < 1.935$ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $3.0<p_T<6.5\,GeV$ and $ -2.865 < y_{CM} < -2 $ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $3.0<p_T<6.5\,GeV$ and $ 1 < y_{CM} < 1.935 $ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

The unfolded jet axis decorrelation distribution,$\frac{1}{N} \frac{dN}{d\Delta j}$, as a function of $\Delta j$ for the $0-10\%$, $10-30\%$, $30-50\%$, and $50-80\%$ centrality bins in the $230 < p_{\mathrm{T}} < 300$ GeV interval.

The elliptic anisotropy $v_2\{4\}$ as a function of $\langle N_{trk}^{offline} \rangle$ with 4-subevent method compared between pPb collisions at 8.16 TeV and PbPb collisions at 5.02 TeV with $p_{T}^{POI}>6$ GeV.

The elliptic anisotropy $v_2\{4\}$ as a function of $\langle N_{trk}^{corrected} \rangle$ with 4-subevent method compared between pPb collisions at 8.16 TeV and PbPb collisions at 5.02 TeV with $p_{T}^{POI}>6$ GeV.

The differential cumulant $d_2\{4\}$ as a function of $p_T$ in generator level HIJING for $60 \le N_{trk}^{gen} <120$ are compared to the experimental pPb results with $60 \le N_{trk}^{offline} <120$ and $185 \le N_{trk}^{offline} <250$.

Expected and observed upper limits for the b-quark-associated Higgs boson production cross section times branching fraction of the decay into a b quark pair at 95% CL as functions of $m_\phi$ for the 2018 FH category. The vertical dashed lines indicate the boundaries of usage of the different fit ranges, as reflected in the rightmost column of Table 2.

Expected and observed upper limits for the b-quark-associated Higgs boson production cross section times branching fraction of the decay into a b quark pair at 95% CL as functions of $m_\phi$, corresponding to the Run 2 combination. The vertical dashed line separates the mass range where only the 2017 SL category contributes on its left, from the region where also the 2017 FH and 2018 FH categories contribute on its right. The expected limits from the 2017 SL, 2017 FH, and 2018 FH datasets as well as from the previously published result based on the 2016 dataset are also shown as colored lines.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = +1$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = -1$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = -2$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = -3$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $m_h^{mod+}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the hMSSM scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as a function of $m_{A}$ for cos(beta-alpha) = 0.1, for the 2HDM $\it Type-II$ scenario.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for masses of $m_{A} = m_{H} = 300$ GeV, for the 2HDM $\it Type-II$ scenario.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as a function of $m_{A}$ for cos(beta-alpha) = 0.1, for the 2HDM $\it Flipped$ scenario.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for masses of $m_{A} = m_{H} = 300$ GeV, for the 2HDM $\it Flipped$ scenario.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 140$ GeV.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 600$ GeV.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 900$ GeV.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 1200$ GeV.

Simulated signal shapes for three representative values of the Higgs boson mass $m_\phi$ in the 2017 SL channel.

Simulated signal shapes for three representative values of the Higgs boson mass $m_\phi$ in the 2017 FH channel.

Simulated signal shapes for three representative values of the Higgs boson mass $m_\phi$ in the 2018 FH channel.

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the SL category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the SL category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the SL category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

Measured inclusive fiducial H->ZZ->4l cross section as a function of the center-of-mass energy.

Postfit reconstructed distributions of the Higgs boson transverse momentum.

Postfit reconstructed distributions of the Higgs boson rapidity.

Higgs fiducial cross section in bins of pT of the Higgs boson.

Higgs fiducial cross section in bins of absolute value of rapidity of the Higgs boson.

Correlation matrix of fit in different final states

Correlation matrix of fit in pT(4l) bins

Correlation matrix of fit in abs(y(4l)) bins

Cutflow for VBF Z' to WW in ditau 2016 channel for different signal scenarios

Cutflow for VBF Z' to WW in ditau 2017 channel for different signal scenarios

Cutflow for VBF Z' to WW in ditau 2018 channel for different signal scenarios

Cutflow for VBF Z' to tautau in mutau 2016 channel for different signal scenarios

Cutflow for VBF Z' to tautau in mutau 2017 channel for different signal scenarios

Cutflow for VBF Z' to tautau in mutau 2018 channel for different signal scenarios

Cutflow for VBF Z' to WW in mutau 2016 channel for different signal scenarios

Cutflow for VBF Z' to WW in mutau 2017 channel for different signal scenarios

Cutflow for VBF Z' to WW in mutau 2018 channel for different signal scenarios

Cutflow for VBF Z' to tautau in etau 2016 channel for different signal scenarios

Cutflow for VBF Z' to tautau in etau 2017 channel for different signal scenarios

Cutflow for VBF Z' to tautau in etau 2018 channel for different signal scenarios

Cutflow for VBF Z' to WW in etau 2016 channel for different signal scenarios

Cutflow for VBF Z' to WW in etau 2017 channel for different signal scenarios

Cutflow for VBF Z' to WW in etau 2018 channel for different signal scenarios

Cutflow for VBF Z' to tautau in emu 2016 channel for different signal scenarios

Cutflow for VBF Z' to tautau in emu 2017 channel for different signal scenarios

Cutflow for VBF Z' to tautau in emu 2018 channel for different signal scenarios

Cutflow for VBF Z' to WW in emu 2016 channel for different signal scenarios

Cutflow for VBF Z' to WW in emu 2017 channel for different signal scenarios

Cutflow for VBF Z' to WW in emu 2018 channel for different signal scenarios

Observed and expected background yields for different mass ranges in the mutau channel. The quoted uncertainty includes both the statistical and the systematic components. The last bin contains the overflow bin.

Observed and expected background yields for different mass ranges in the ditau channel. The quoted uncertainty includes both the statistical and the systematic components. The last bin contains the overflow bin.

Observed and expected background yields for different mass ranges in the etau channel. The quoted uncertainty includes both the statistical and the systematic components. The last bin contains the overflow bin.

Observed and expected background yields for different mass ranges in the emu channel. The quoted uncertainty includes both the statistical and the systematic components. The last bin contains the overflow bin.

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to tautau, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 0 and $k_V = 0.1.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to tautau, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 0 and $k_V = 0.25.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to tautau, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 0 and $k_V = 0.50.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to tautau, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 0 and $k_V = 0.75.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to tautau, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 0 and $k_V = 1.0.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to tautau, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 1 and $k_V = 0.1.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to tautau, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 1 and $k_V = 0.25.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to tautau, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 1 and $k_V = 0.50.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to tautau, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 1 and $k_V = 0.75.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to tautau, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 1 and $k_V = 1.0.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to WW, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 0 and $k_V = 0.1.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to WW, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 0 and $k_V = 0.25.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to WW, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 0 and $k_V = 0.50.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to WW, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 0 and $k_V = 0.75.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to WW, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 0 and $k_V = 1.0.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to WW, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 1 and $k_V = 0.1.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to WW, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 1 and $k_V = 0.25.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to WW, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 1 and $k_V = 0.50.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to WW, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 1 and $k_V = 0.75.$

Combined $95\%$ confidence level expected and observed upper limits on the product of the $\mathrm{Z'}$ cross section and the branching fraction for a $\mathrm{Z'}$ boson decaying to WW, as a function of $\mathrm{Z'}$ mass and assuming g_(1,2) = 1 and $k_V = 1.0.$

Cross sections for different signal models for VBF Z'. To calculate cross section for branching fraction (br) 0.20 and kV = 0.1, subsitute the values in (br {x} $\kappa_{V}^{2}$ {x} xsec)

Combined 95% CL lower limits on m(Z') as a function of the Z' branching fraction to tautau for g_(1,2)=0 scenario. The observed limits corresponding to $\kappa_V$ equal to 0.1, 0.5, and 1.0 are shown here.

Combined 95% CL lower limits on m(Z') as a function of the Z' branching fraction to tautau for g_(1,2)=1 scenario. The observed limits corresponding to $\kappa_V$ equal to 0.1, 0.5, and 1.0 are shown here.

Combined 95% CL lower limits on m(Z') as a function of the Z' branching fraction to WW for g_(1,2)=0 scenario. The observed limits corresponding to $\kappa_V$ equal to 0.1, 0.5, and 1.0 are shown here.

Combined 95% CL lower limits on m(Z') as a function of the Z' branching fraction to WW for g_(1,2)=1 scenario. The observed limits corresponding to $\kappa_V$ equal to 0.1, 0.5, and 1.0 are shown here.