Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Signal yields for two charge scenarios considered in the analysis, as well as their associated uncertainties.

Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Signal yields for two charge scenarios considered in the analysis, as well as their associated uncertainties.

Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the early 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 6%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text. Comparing with respect to the binomial fit starting at $N_{\text{hits}}^{\text{low dE/dx}} = 2$, and not $N_{\text{hits}}^{\text{low dE/dx}} = 1$, is needed to account for the fact that early 2016 data is more strongly affected low dE/dx by instrumental effects that widen the N hits distribution.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the early 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 6%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text. Comparing with respect to the binomial fit starting at $N_{\text{hits}}^{\text{low dE/dx}} = 2$, and not $N_{\text{hits}}^{\text{low dE/dx}} = 1$, is needed to account for the fact that early 2016 data is more strongly affected low dE/dx by instrumental effects that widen the N hits distribution.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the late 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 78%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the late 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 78%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2017 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 65%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2017 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 65%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2018 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 9%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2018 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 9%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Exclusion region (hatched) at 95% CL in the FCP charge-mass plane for the considered signal. The expected exclusion is shown with the associated 1 and 2 standard deviations $\sigma$ bands. Signal points at charges 0.9, 0.8, 2/3, 0.5, and 1/3 e are connected by straight lines to guide the eye. This is a conservative interpolation. Previous exclusions from CMS [Phys. Rev. D 87 (2013) 092008), JHEP 07 (2013) 122] as well as OPAL [Phys. Lett. B 572 (2003) 8] are given for comparison.

Exclusion region (hatched) at 95% CL in the FCP charge-mass plane for the considered signal. The expected exclusion is shown with the associated 1 and 2 standard deviations $\sigma$ bands. Signal points at charges 0.9, 0.8, 2/3, 0.5, and 1/3 e are connected by straight lines to guide the eye. This is a conservative interpolation. Previous exclusions from CMS [Phys. Rev. D 87 (2013) 092008), JHEP 07 (2013) 122] as well as OPAL [Phys. Lett. B 572 (2003) 8] are given for comparison.

Particle-level distributions of ungroomed AK4 pTD2 in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 thrust in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 thrust in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 width in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 width in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 LHA in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 LHA in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 LHA in 408 < PT < 1500 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 LHA in 1000 < PT < 4000 GeV in the central dijet region.

Particle-level distributions of ungroomed AK8 LHA in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK8 LHA in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 LHA (charged-only) in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 LHA (charged-only) in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of groomed AK4 LHA in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of groomed AK4 LHA in 120 < PT < 150 GeV in the central dijet region.

Mean of ungroomed LHA for AK4 jets as a function of PT in the Z+jet region.

Mean of ungroomed LHA for AK4 jets as a function of PT in the central dijet region.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of dimuon invariant mass. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Forward-backward ratios for $15<m_{\mu\mu}<60$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Forward-backward ratios for $60<m_{\mu\mu}<120$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of dimuon invariant mass.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV.

The observed DNN output distribution for data collected in 2016 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution for data collected in 2017 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution for data collected in 2018 in the VBF-SB region compared to the post-fit background estimate from SM processes. The predicted backgrounds are obtained from a S+B fit performed across analysis regions and years. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution in the VBF-SB regions for the combination of 2016, 2017, and 2018 data, compared to the post-fit prediction from SM processes. The total post-fit and pre-fit uncertainties on the background prediction are also reported.

The observed DNN output distribution in the VBF-SR region for the combination of 2016, 2017, and 2018 data, compared to the post-fit prediction from SM processes. The post-fit distributions for the Higgs boson signal produced via ggH and VBF modes with mass of 125.38 GeV, along with the total post-fit and pre-fit uncertainties on the background prediction, are also reported.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ggH category, with dimuon mass between 110–150 GeV, are considered. The expected distributions for ggH, VBF, and other signal processes are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat3 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat4 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ggH-cat5 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ttH hadronic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ttH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHhad-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHhad-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHhad-cat3 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHlep-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ttHlep-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the WH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

The observed BDT output distribution compared to the prediction from the simulation of various SM background processes. Dimuon events passing the event selection requirements of the ZH leptonic category, with dimuon mass between 110–150 GeV, are considered. The expected distributions from different production modes of the Higgs boson are also reported. The background uncertainty is obtained by summing in quadrature the uncertainty due to the limited size of the simulated samples with those arising from experimental and theoretical sources of systematic uncertainty.

Comparison between the data and the total background extracted from a S+B fit performed in the WH-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the WH-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the WH-cat3 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ZH-cat1 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Comparison between the data and the total background extracted from a S+B fit performed in the ZH-cat2 subcategory. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Observed local p-values as a function of the Higgs boson mass hypothesis extracted from the combined fit across event categories as well as from each individual production category.

Expected local p-values as a function of the Higgs boson mass hypothesis extracted from the combined fit across event categories as well as from each individual production category. They are calculated using the background-only expectation obtained from the S+B fit and injecting a signal with mass of 125.38 GeV and strenght of one.

Signal strength modifiers measured for a Higgs boson with mass of 125.38 GeV in each production category as well as from the combined fit.

Scan of the profiled likelihood ratio as a function of $\mu_{ggH,t\overline{t}H}$ and $\mu_{VBF,VH}$.

The 68% CL contour for the profiled likelihood scan as a function of $\mu_{ggH,t\overline{t}H}$ and $\mu_{VBF,VH}$.

The 95% CL contour for the profiled likelihood scan as a function of $\mu_{ggH,t\overline{t}H}$ and $\mu_{VBF,VH}$.

The mass distribution for the weighted combination of VBF-SB and VBF-SR events. Each event is weighted proportionally to the S/(S+B) ratio, calculated as a function of the mass-decorrelated DNN output. The total post-fit and pre-fit uncertainties on the background prediction, along with the best fit signal contribution, are also reported.

Comparison between the data and the total background extracted for the S/(S+B) weighted combination of all event categories. The one and two standard deviation bands include the uncertainties in the background component of the fit. The best fit signal contribution with mass of 125.38 GeV is also reported.

Observed and expected local p-values as a function of the Higgs boson mass hypothesis extracted from the combined fit performed on data recorded at centre-of-mass energies of 7, 8, and 13 TeV. The expected p-values are calculated using the background-only expectation obtained from the S+B fit and injecting a signal with mass of 125.38 GeV and strenght of one.

The observed profile likelihood ratio as a function of $\kappa_{\mu}$ for a Higgs boson mass of 125.38 GeV, obtained from a combined fit with Ref.[10] in the resolved $\kappa$-framework. The best fit value for $\kappa_{\mu}$ is 1.07 and the corresponding observed 68% CL interval is $0.85 < \kappa_{\mu} < 1.29$.

The best fit estimates for the reduced coupling modifiers extracted for fermions and weak bosons from the resolved $\kappa$-framework compared to their corresponding prediction from the SM. The error bars represent 68% CL intervals for the measured parameters.

The distribution of Top jet $ au_{3}/ au_{2}$ for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of W jet $ au_{2}/ au_{1}$ for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of m_{W} for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of Top jet DeepCSV score for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of m_{t} for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

Distributions of $m_{tW}$ in the signal region for an interval of $m_{t}$ (65-105 GeV).

Distributions of $m_{tW}$ in the signal region for an interval of $m_{t}$ (105-225 GeV).

Distributions of $m_{tW}$ in the signal region for an interval of $m_{t}$ (225-285 GeV).

Distributions of $m_{tW}$ in themultijet control region of the signal region for an interval of $m_{t}$ (65-105 GeV).

Distributions of $m_{tW}$ in themultijet control region of the signal region for an interval of $m_{t}$ (105-225 GeV).

Distributions of $m_{tW}$ in themultijet control region of the signal region for an interval of $m_{t}$ (225-285 GeV).

Distributions of $m_{t}$ in the ttbar measurement region for an interval of $m_{tar{t}}$ (1200-1300 GeV).

Distributions of $m_{t}$ in the ttbar measurement region for an interval of $m_{tar{t}}$ (1300-1800 GeV).

Distributions of $m_{t}$ in the ttbar measurement region for an interval of $m_{tar{t}}$ (1800-3000 GeV).

Distributions of $m_{t}$ in themultijet control region of the ttbar measurement region for an interval of $m_{tar{t}}$ (1200-1300 GeV).

Distributions of $m_{t}$ in themultijet control region of the ttbar measurement region for an interval of $m_{tar{t}}$ (1300-1800 GeV).

Distributions of $m_{t}$ in themultijet control region of the ttbar measurement region for an interval of $m_{tar{t}}$ (1800-3000 GeV).

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the HP VBF category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the HP ggF/DY category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the LP VBF category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

The distributions of the transfer factor ($\alpha$) versus $m_T$ in the LP ggF/DY category is shown. The last bin corresponds to the value obtained by integrating events above the penultimate bin.

Distributions of mT for high-purity ggF/DY CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for high-purity VBF CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity ggF/DY CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity VBF CR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for high-purity ggF/DY SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for high-purity VBF SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity ggF/DY SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Distributions of mT for low-purity VBF SR events after performing background-only fits. The last bin corresponds to the yields integrated above the penultimate bin. Both the ggF and VBF-produced 1 TeV graviton signals (normalized to 10 fb) are shown in the plot. Nonresonant background corresponds to the Z+jets and W+jets events. Resonant background corresponds to the ttbar, single top, and di/tri-boson events.

Expected and observed 95% CL upper limits on the product of the radion (R) production (gluon-gluon fusion) cross section and the R -> ZZ branching fraction versus the radion mass.

Expected and observed 95% CL upper limits on the product of the radion (R) production (vector boson fusion) cross section and the R -> ZZ branching fraction versus the radion mass.

Expected and observed 95% CL upper limits on the product of the W' (W') production (Drell-Yan) cross section and the W' -> WZ branching fraction versus the W' mass.

Expected and observed 95% CL upper limits on the product of the W' (W') production (vector boson fusion) cross section and the W' -> WZ branching fraction versus the W' mass.

Expected and observed 95% CL upper limits on the product of the graviton (G) production (gluon-gluon fusion) cross section and the G -> ZZ branching fraction versus the graviton mass.

Expected and observed 95% CL upper limits on the product of the graviton (G) production (vector boson fusion) cross section and the G -> ZZ branching fraction versus the graviton mass.

Muon $p_T$ distribution using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$p_T^{miss}$ distribution in the muon channel using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$M_T$ distribution for the E+INVISIBLE system using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$M_T$ distribution for the MU+INVISIBLE system using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the electron decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the muon decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the combined electron and muon decay channels. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the electron decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the muon decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the combined electron and muon decay channels. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the coupling strength ratio, gWprime/gW as a function of the mass of the WPRIME boson, for the electron channel. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. The dotted line represents the case of SSM coupling, gWprime, being equal to gW, the SM coupling.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the coupling strength ratio, gWprime/gW, as a function of the mass of the WPRIME boson, for the muon channel. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. The dotted line represents the case of SSM coupling, gWprime , being equal to gW, the SM coupling.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the electron channel for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the muon channel for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the combined electron and muon channels for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the ${\lambda}_{231}$ coupling in the RPV SUSY model with a STAU mediator, as a function of the mass of STAU, for the electron channel and for the production coupling, ${\lambda}^{'}_{3ij}$=0.5 . The couplings ${\lambda}^{'}_{3ij}$, ${\lambda}_{231}$, and ${\lambda}_{132}$ are defined in Fig. 1. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. Limits for other values of ${\lambda}^{'}_{3ij}$ are obtained multiplying the ones shown by the ratio ${\lambda}^{'}_{3ij}$/0.5

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the ${\lambda}_{132}$ coupling in the RPV SUSY model with a stau mediator, as a function of the mass of the stau, for the muon channel and for the production coupling, ${\lambda}^{'}_{3ij}$ = 0.5 . The couplings ${\lambda}^{'}_{3ij}$, ${\lambda}_{231}$, and ${\lambda}_{132}$ are defined in Fig. 1. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. Limits for other values of ${\lambda}^{'}_{3ij}$ are obtained multiplying these ones by the ratio ${\lambda}^{'}_{3ij}$/0.5

Region in the oblique (W,Y) parameter phase space allowed by the current analysis at 95% CL. The combination of the electron and muon channel distributions from the 2017 and 2018 data sets at sqrt(s) = 13 TeV is used, having 101 fb-1 of luminosity.