Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hq}^{(3)}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hu}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $c_{Hd}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(1)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hq}^{(3)}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hu}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(2)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $c_{Hd}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while allowing the other coefficients to float freely at each point of the sca.

Observed two-dimensional likelihood scans for $g^{(2)}_{ZZ}$ vs. $g^{(4)}_{ZZ}$ while other coefficients are fixed at their SM values.

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the muon channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{T}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the invisible channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{T}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the invisible channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the electron channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Distributions in $m_{Z'}^{rec}$ for data in the SRs, together with fits of the background functions under the background-only hypothesis for the electron channel. The number of observed events in each bin is divided by the bin width. The signal predictions are shown for different Z' boson masses, normalized to an arbitrary cross section of 1 fb. In the panels below the distributions, the ratios of data to the background function are displayed. The shaded green areas represent the statistical uncertainty from the fit. The $\chi^2$ values per number of degrees of freedom ($\chi^2$/n.d.f.) and the corresponding $p$-values are provided for each fit.

Expected upper limits at 95% CL on the product of the production cross section and the branching fraction as a function of the Z' boson mass.Each line corresponds to the three different final states and their combined result.

Expected upper limits at 95% CL on the product of the production cross section and the branching fraction as a function of the Z' boson mass.Each line corresponds to the three different final states and their combined result.

Expected and observed upper limits at 95% CL on the product of the production cross section and the branching fraction as functions of the Z' boson mass from the combination of all final states. The limits are compared with the predictions of the HVT model and the expected limits (shown by red curves) from a previous analysis.

Expected and observed upper limits at 95% CL on the product of the production cross section and the branching fraction as functions of the Z' boson mass from the combination of all final states. The limits are compared with the predictions of the HVT model and the expected limits (shown by red curves) from a previous analysis.

Prediction from the HVT model.

Prediction from the HVT model.

Observed upper limits at 95% CL on gF for different Z' boson masses as functions of the product of $g_{H}$ with the sign of $g_{F}$. The two benchmark scenarios of the HVT model are shown by the black markers.

Observed upper limits at 95% CL on gF for different Z' boson masses as functions of the product of $g_{H}$ with the sign of $g_{F}$. The two benchmark scenarios of the HVT model are shown by the black markers.

Distributions for data and expected signal and background contributions of the MVA score for the µ + jets channel in the 3j1b category, before the maximum likelihood fit. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty in the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-prediction ratio. The first and last bins in each distribution include underflow and overflow events, respectively.

Observed and predicted number of events in each of the eight categories of the signal region, before the maximum likelihood fit. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty in the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-prediction ratio.

Distributions for the e + jets final state before the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization before the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the e + jets final state after the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization after the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the µ + jets final state before the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization before the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Distributions for the µ + jets final state after the maximum likelihood fit, MVA score bins for the 3j1b category and ∆R med ( j, j 0 ) bins for the other categories. The vertical error bars represent the statistical uncertainty in the data, and the shaded band the uncertainty of the prediction. All uncertainties considered in the analysis are included in the uncertainty band. The lower panels show the data-to-simulation ratio, where the simulations are adjusted according to their normalization after the fit. The first and last bins in each distribution include underflow and overflow events, respectively.

Covariance matrix for the lepton+jets ML fit.

Covariance matrix for the combination ML fit.

Distributions in $S_T$ in the SR for the electron channel, after a background-only fit to the data. The signal distributions are scaled to the cross section predicted by the theory. The hatched bands show the post-fit uncertainty band, combining all sources of uncertainty. The ratio of data to the background predictions is shown in the panels below the distributions.

The $G_{\text{i}}^{\text{Strips}}$ distribution in the PASS region for events passing the event selection and with $55 < p_{\mathrm{T}} < 200$ GeV.

The $G_{\text{i}}^{\text{Strips}}$ distribution in the FAIL region for events passing the event selection and with $p_{\mathrm{T}} > 200$ GeV.

The $G_{\text{i}}^{\text{Strips}}$ distribution in the PASS region for events passing the event selection and with $p_{\mathrm{T}} > 200$ GeV.

Mass spectrum predicted in the signal region defined by $G_{\text{i}}^{\text{Strips}} > 0.22$ and $p_{\mathrm{T}} > 70$ GeV.

Cross section limits for gluino and supersymmetric top R-hadrons for both background prediction methods.

Cross section limits for supersymmetric tau models for both background prediction methods.

Cross section limits for DY-produced tau prime models (single and multicharged) for both background prediction methods.

Cross section limits for Z prime to multicharged tau prime model for both background prediction methods. All Z prime models assume a narrow width, and a 100% branching fraction to 600 GeV tau primes.

2D exclusion showing the observed cross section limit as a function of the multicharged tau prime mass and Z prime mass for the ionization method.

2D exclusion showing the observed cross section limit as a function of the multicharged tau prime mass and Z prime mass for the mass method.

Cumulative selection efficiency for the data and for two signal hypotheses.

Expected and observed mass limits obtained using 2017-2018 data for various HSCP candidate models,for the two background estimate methods.

Mass windows used in the mass method as a function of the signal target mass for the signal samples assuming a charge of $1e$ and $2e$, as appropriate to the signal model.

Mass spectrum predicted in the validation region defined by $0.018 < G_{\text{i}}^{\text{Strips}} < 0.057$ and $p_{\mathrm{T}} > 70$ GeV. The thresholds used in the $G_{\text{i}}^{\text{Strips}}$ requirement represent the 50% and 90% quantile of the distribution.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $\text{abs(}\eta\text{)}<0.3$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.3<\text{abs(}\eta\text{)}<0.6$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.6<\text{abs(}\eta\text{)}<0.9$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $0.9<\text{abs(}\eta\text{)}<1.2$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $1.2<\text{abs(}\eta\text{)}<2.1$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{g}$ signals as a function of $\beta$ for $2.1<\text{abs(}\eta\text{)}<2.4$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $\text{abs(}\eta\text{)}<0.3$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.3<\text{abs(}\eta\text{)}<0.6$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.6<\text{abs(}\eta\text{)}<0.9$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $0.9<\text{abs(}\eta\text{)}<1.2$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $1.2<\text{abs(}\eta\text{)}<2.1$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency for the $\tilde{\tau}$ signals as a function of $\beta$ for $2.1<\text{abs(}\eta\text{)}<2.4$. The up and down variations are conservatively estimated assuming a delay of 1.5 ns in the muon chambers (equivalent to the time resolution of the chambers) and are used to evaluate the signal systematic uncertainties.

Trigger efficiency as a function of $\beta$, for the $\tilde{g}$ signals.

Trigger efficiency as a function of $\beta$, for the $\tilde{\tau}$ signals.

Cross section limits for $\tilde{g}$ R-hadrons obtained with the ionization method.

Cross section limits for $\tilde{g}$ R-hadrons obtained with the mass method.

Cross section limits for $\tilde{t}$ R-hadrons obtained with the ionization method.

Cross section limits for $\tilde{t}$ R-hadrons obtained with the mass method.

Cross section limits for $\tilde{\tau}$ obtained with the ionization method.

Cross section limits for $\tilde{\tau}$ obtained with the mass method.

Cross section limits for $\tilde{\tau}$ production within the GMSB SPS7 modelobtained with the ionization method.

Cross section limits for $\tilde{\tau}$ production within the GMSB SPS7 model obtained with the mass method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 1e$ with the ionization method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 1e$ with the mass method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 2e$ with the ionization method.

Cross section limits for DY-produced $\tau'$ with $abs(Q) = 2e$ with the mass method.

Cross section limits for the production of $Z'$ boson decaying into a pair of $\tau'$ fermions of charge $2e$ (with a branching fraction equal to 1 and a fixed $\tau'$ mass of 600 GeV), obtained with the ionization method.

Cross section limits for the production of $Z'$ boson decaying into a pair of $\tau'$ fermions of charge $2e$ (with a branching fraction equal to 1 and a fixed $\tau'$ mass of 600 GeV), obtained with the mass method.

Two-dimensional exclusion showing the observed cross section limit as a function of the masses of the $\tau'$ (on the $x$ axis) and of the $Z'$ boson (on the $y$ axis), for the ionization method.

Two-dimensional exclusion showing the observed cross section limit as a function of the masses of the $\tau'$ (on the $x$ axis) and of the $Z'$ boson (on the $y$ axis), for the mass method.

Expected background yield, expected signal yield, and observed data for the mass method for 2017.

Expected background yield, expected signal yield, and observed data for the mass method for 2018.

Selection efficieny for GMSB supersymmetric $\tau$.

Selection efficiency for gluino R-hardon samples.

Selection efficiency for supersymmetric top R-hadrons.

Selection efficiency for pair-produced supersymmetric $\tau$.

Selection efficiency for $\tau'$ with $|Q| = 1e$.

Selection efficiecny for $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ (mass = 3 TeV) going to $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ (mass 5 TeV) goes to $\tau'$ with $|Q| = 2e$.

Selection efficiency for $Z'$ going to $\tau'$ (mass 600 GeV) with $|Q| = 2e$.

Selection efficiency for $Z'$ going to $\tau'$ (mass 1 TeV) with charge $|Q| = 2e$.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH3 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 and ISR starting scale SpaceShower:pT0Ref = 1 GeV at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH3 and ISR starting scale SudakovCommon:pTmin=0.7 GeV at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-proton center-of-mass energy 38.8 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-proton center-of-mass energy 8000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-nucleon center-of-mass energy 8160 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-nucleon center-of-mass energy 13000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 38.8 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 8000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-nucleon center-of-mass energy 8160 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 13000 GeV versus the dilepton mass range of the process.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Covariance for full matrix measurement all bins from $m(t\bar{t})$ fit

Results for full matrix measurement inclusive from $p_{T}(t)$

Covariance for full matrix measurement inclusive from $p_{T}(t)$

Results for full matrix measurement all bins from $p_{T}(t)$ fit

Covariance for full matrix measurement all bins from $p_{T}(t)$ fit

Results for full matrix measurement differential in $m(t\bar{t})$

Covariance for full matrix measurement differential in $m(t\bar{t})$

Results for full matrix measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Covariance for full matrix measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Results for full matrix measurement differential in $p_{T}(t)$

Covariance for full matrix measurement differential in $p_{T}(t)$

Results for $D$ measurement inclusive from $m(t\bar{t})$

Covariance for $D$ measurement inclusive from $m(t\bar{t})$

Results for $D$ measurement inclusive from $p_{T}(t)$

Covariance for $D$ measurement inclusive from $p_{T}(t)$

Results for $\tilde{D}$ measurement inclusive in $m(t\bar{t})$

Covariance for $\tilde{D}$ measurement inclusive in $m(t\bar{t})$

Results for $\tilde{D}$ measurement inclusive in $p_{T}(t)$

Covariance for $\tilde{D}$ measurement inclusive in $p_{T}(t)$

Results for $D$ measurement all bins from from $m(t\bar{t})$ fit

Covariance for $D$ measurement all bins from from $m(t\bar{t})$ fit

Results for $D$ measurement all bins from from $p_{T}(t)$ fit

Covariance for $D$ measurement all bins from from $p_{T}(t)$ fit

Results for $\tilde{D}$ measurement all bins from from $m(t\bar{t})$ fit

Covariance for $\tilde{D}$ measurement all bins from from $m(t\bar{t})$ fit

Results for $\tilde{D}$ measurement all bins from from $p_{T}(t)$ fit

Covariance for $\tilde{D}$ measurement all bins from from $p_{T}(t)$ fit

Results for $D$ measurement differential in $m(t\bar{t})$

Covariance for $D$ measurement differential in $m(t\bar{t})$

Results for $\tilde{D}$ measurement differential in $m(t\bar{t})$

Covariance for $\tilde{D}$ measurement differential in $m(t\bar{t})$

Results for $\tilde{D}$ measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Covariance for $\tilde{D}$ measurement differential in $m(t\bar{t})$ for $|cos(\theta)| < 0.4$

Results for $D$ measurement differential in $p_{T}(t)$

Covariance for $D$ measurement differential in $p_{T}(t)$

Results for $\tilde{D}$ measurement differential in $p_{T}(t)$

Covariance for $\tilde{D}$ measurement differential in $p_{T}(t)$

The Bs meson $p_{\rm{T}}$-dependent RAA in PpPp. The measurment was carried out inside a fiducial region respecting ($p_{\rm{T}}$<10 & 1.5<|y|<2.4) and ($p_{\rm{T}}$>10 & |y|<2.4).

Distribution of the RF output for events in the 2j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

The distribution of the Subleading jet $p_{T}$ for events in the 2j2b region. The number of observed events (points) and estimated signal and background events (filled histograms) from the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events after the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) after the fit, with the bands giving the corresponding uncertainties.

Distribution of the Subleading jet $p_{T}$ for events in the 2j2b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

Observed and theoretical cross sections. In the observed, the first uncertainty is statistical, the second is the systematic, and the third the luminosity. In the theoretical, the first uncertainty is due to scale variations, the second due to the choice of PDF. The theoretical cross section has been computed at approximate third-order in QCD, with the addition of third-order corrections of soft gluon emission terms, assuming a top quark mass of 172.5 GeV and using the PDF4LHC21 PDF set.

The distribution of the RF output for events in the 1j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) from the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events after the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) after the fit, with the bands giving the corresponding uncertainties.

The number of estimated signal and background events after the fit in the 1j1b, 2j1b, and 2j2b regions compared to the observed number of events. The total uncertainties in the estimated events after the fit are given.

The distribution of the RF output for events in the 2j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) from the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events after the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) after the fit, with the bands giving the corresponding uncertainties.

Normalised fiducial differential tW production cross section as a function of the $Leading$ $lepton$ $p_{T}$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

The distribution of the Subleading jet $p_{T}$ for events in the 2j2b region. The number of observed events (points) and estimated signal and background events (filled histograms) from the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events after the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) after the fit, with the bands giving the corresponding uncertainties.

Normalised fiducial differential tW production cross section as a function of the $p_{Z}(e^{\pm}, \mu^{\pm}, j)$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Observed and theoretical cross sections. In the observed, the first uncertainty is statistical, the second is the systematic, and the third the luminosity. In the theoretical, the first uncertainty is due to scale variations, the second due to the choice of PDF. The theoretical cross section has been computed at approximate third-order in QCD, with the addition of third-order corrections of soft gluon emission terms, assuming a top quark mass of 172.5 GeV and using the PDF4LHC21 PDF set.

Normalised fiducial differential tW production cross section as a function of the $Jet$ $p_{T}$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

The number of estimated signal and background events after the fit in the 1j1b, 2j1b, and 2j2b regions compared to the observed number of events. The total uncertainties in the estimated events after the fit are given.

Normalised fiducial differential tW production cross section as a function of the $m(e^{\pm}, \mu^{\pm}, j)$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $Leading$ $lepton$ $p_{T}$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $p_{Z}(e^{\pm}, \mu^{\pm}, j)$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Normalised fiducial differential tW production cross section as a function of the $Jet$ $p_{T}$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

The p-values from the goodness-of-fit tests comparing the six differential cross section measurements with the predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5 aMC@NLO (aMC) DR, DR2, DS, and DS with a dynamic factor + PYTHIA 8. The complete covariance matrix from the results and the statistical uncertainties in the predictions are taken into account.

Normalised fiducial differential tW production cross section as a function of the $m(e^{\pm}, \mu^{\pm}, j)$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Distribution of the RF output for events in the 1j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

Normalised fiducial differential tW production cross section as a function of the $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Distribution of the RF output for events in the 2j1b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

Normalised fiducial differential tW production cross section as a function of the $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$. The horizontal bars on the points show the bin width. Predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5_aMC@NLO (aMC) DR, DR2, DS and DS with a dynamic factor + PYTHIA 8 are also shown. The grey band represents the statistical uncertainty and the yellow band the total uncertainty. In the lower panels, the ratio of the predictions to the data is shown.

Distribution of the Subleading jet $p_{T}$ for events in the 2j2b region. The number of observed events (points) and estimated signal and background events (filled histograms) before the maximum likelihood fit are shown. The vertical bars on the points represent the statistical uncertainty in the data, and the hatched band the total uncertainty in the estimated events before the fit. The lower panels display the ratio of the data to the sum of the estimated events (points) before the fit, with the bands giving the corresponding uncertainties.

The p-values from the goodness-of-fit tests comparing the six differential cross section measurements with the predictions from POWHEG (PH) DR and DS + PYTHIA 8 (P8), POWHEG DR + HERWIG 7 (H7), MADGRAPH5 aMC@NLO (aMC) DR, DR2, DS, and DS with a dynamic factor + PYTHIA 8. The complete covariance matrix from the results and the statistical uncertainties in the predictions are taken into account.

Response matrix between detector and particle level for the $Leading$ $lepton$ $p_{T}$.

Response matrix between detector and particle level for the $Leading$ $lepton$ $p_{T}$.

Response matrix between detector and particle level for the $p_{Z}(e^{\pm}, \mu^{\pm}, j)$.

Response matrix between detector and particle level for the $p_{Z}(e^{\pm}, \mu^{\pm}, j)$.

Response matrix between detector and particle level for the $Jet$ $p_{T}$.

Response matrix between detector and particle level for the $Jet$ $p_{T}$.

Response matrix between detector and particle level for the $m(e^{\pm}, \mu^{\pm}, j)$.

Response matrix between detector and particle level for the $m(e^{\pm}, \mu^{\pm}, j)$.

Response matrix between detector and particle level for the $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$.

Response matrix between detector and particle level for the $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$.

Response matrix between detector and particle level for the $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$.

Response matrix between detector and particle level for the $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $Leading$ $lepton$ $p_{T}$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $Leading$ $lepton$ $p_{T}$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $p_{Z}(e^{\pm}, \mu^{\pm}, j)$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $p_{Z}(e^{\pm}, \mu^{\pm}, j)$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $Jet$ $p_{T}$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $Jet$ $p_{T}$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $m(e^{\pm}, \mu^{\pm}, j)$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $m(e^{\pm}, \mu^{\pm}, j)$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$ in units of 1.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $\Delta\phi(e^{\pm}, \mu^{\pm})/\pi$ in units of 1.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$ in units of $1/GeV^{2}$.

Covariance matrix including all uncertainties of the normalised differential cross section in bins of $m_{T}(e^{\pm}, \mu^{\pm}, j, p_{T}^{miss})$ in units of $1/GeV^{2}$.

Observed and predicted $p_T^{miss}$ ($U$) distributions in the 1b category, in the Signal Region. The bottom plot shows the ratio of observed and predicted (both prefit and postfit) distributions along with an uncertainty band that includes both the systematic and statistical components.

Observed and predicted $p_T^{miss}$ ($U$) distributions in the 1b category, in the $Z(\ell\ell)$ Control Region. The bottom plot shows the ratio of observed and predicted (both prefit and postfit) distributions along with an uncertainty band that includes both the systematic and statistical components.

Observed and predicted $p_T^{miss}$ ($U$) distributions in the 1b category, in the $W(\ell\nu)$ Control Region. The bottom plot shows the ratio of observed and predicted (both prefit and postfit) distributions along with an uncertainty band that includes both the systematic and statistical components.

Observed and predicted $\cos\Theta$* distributions in the 2b category, in the Signal Region. The bottom plot shows the ratio of observed and predicted (both prefit and postfit) distributions along with an uncertainty band that includes both the systematic and statistical components.

Observed and predicted $\cos\Theta$* distributions in the 2b category, in the $Z(\ell\ell)$ Control Region. The bottom plot shows the ratio of observed and predicted (both prefit and postfit) distributions along with an uncertainty band that includes both the systematic and statistical components.

Observed and predicted $\cos\Theta$* distributions in the 2b category, in the $t\bar{t}$ Control Region. The bottom plot shows the ratio of observed and predicted (both prefit and postfit) distributions along with an uncertainty band that includes both the systematic and statistical components.

The $95\%$ CL upper limit on the signal strength modifier of DM produced in association with a pair of bottom quarks for $m_A=600$ GeV, $\sin\theta=0.7$, $m_\chi=1$ GeV, and $\tan\beta=35$, for the combination of SR1 and SR2. The green and yellow bands show the $\pm 1$ and $\pm 2$ standard deviations from expected limits. The mass points below the red line are excluded.

$95\%$ CL upper limit on the signal strength in the $m_a-\tan\beta$ plane. The shaded area above the curve is excluded.

$95\%$ CL upper limit on the signal strength in the $m_a-\sin\theta$ plane. The shaded area above the curve is excluded.

$95\%$ CL upper limit on the signal strength in the $m_a-m_\chi$ plane. The shaded area below the curve is excluded.