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.

Distributions of $m_T$ in the $W^{+}$ signal selection for mu final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for e final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for mu final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for e final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_T$ in the $W^{-}$ signal selection for mu final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to\tau\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for ee final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for mumu final states for the pp collisions at $\sqrt{s}=$ 5TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for ee final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Distributions of $m_{ll}$ in the Z signal selection for mumu final states for the pp collisions at $\sqrt{s}=$ 13TeV after the maximum likelihood fit. The EW backgrounds include the contributions from DY, $W\to l\nu$, and diboson processes.

Figure 4 Fiducial cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 13TeV

Figure 5 Fiducial cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 5.02TeV

Figure 6 Fiducial cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes between 13TeV and 5.02TeV

Figure 7 Total cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 13TeV

Figure 8 Total cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes at 5.02TeV

Figure 9 Total cross sections and ratios for $W\to\ell\nu$ and $Z\to\ell\ell$ processes between 13TeV and 5.02TeV

Summary of theoretical predictions of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at ppbar collisions.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at UA1 and UA2.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at D0.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at CDF.

Summary of theoretical predictions of the total inclusive cross sections for W boson production times leptonic branching fractions at pp collisions.

Summary of experimental measurements of the total inclusive cross sections for Z boson production times leptonic branching fractions at pp collisions.

Summary of experimental measurements of the total inclusive cross sections for W and Z boson production times leptonic branching fractions at CMS.

The model-independent 95\% \CL expected and observed upper limits set on ${\sigma(\PP\to 2\Pa+\PX)\mathcal{B}^2(\Pa\to 2\PGm)\alphaGen}$ over the range $0.21 < \MPa < 9\GeV$ for the combined 2016, 2017, and 2018 analyses. Mass ranges that overlap with \JPsi and \PgU resonances are excluded from the search

The 95\% \CL observed upper limits on the effective coupling $\CahLambda$ of the ALP to the SM Higgs assuming the branching fraction of the ALP to muons is 1 (blue) and 0.1 (orange), for both the expected (dashed) and observed (solid) limits

The 95\% \CL observed upper limits on the effective coupling $\cll/\Lambda$ of the ALP to the SM leptons. The orange shaded region represents the parameter space excluded by this search under three choices of $\CahLambda$ $1\TeV^{-2}$ (solid), $0.1\TeV^{-2}$ (dashed), and $0.01\TeV^{-2}$ (dotted)

The 95\% \CL observed upper limits on $\varepsilon^{2} \mathcal{B}(\ZDark \to \sDark \overline{\sDark}) \mathcal{B}^{2}(\sDark \rightarrow 2\PGm)$ for the vector portal model as a function of the dark scalar mass $\MsDark$ and dark vector boson mass $\MZDark$. The yellow and green bands indicate one and two standard deviation values of the limit, respectively. Because the model-independent limits are calculated only up to a dimuon mass of 60\GeV, the \sDark mass considered for these limits is below 60\GeV

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 90~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 125~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits for $\sigma(\PP\to\PhOneTwo\to 2\PaOne) \mathcal{B}^2(\PaOne\to 2\PGm)$ for the NMSSM as a function of $\MPaOne$ for $\MPhOneTwo$ = 150~GEV. The data presented here reflect the results of the 2017 and 2018 datasets combined with the previously published results of the 2016 dataset in \cite{CMS:2018jid}

The 95\% \CL observed upper limits (black solid curves) from this search as interpreted in the dark SUSY scenario for the process $\PP \to \Ph \to 2\nOne \to 2\gammaDark + 2\nDark \to 4\PGm + \PX$, with \mbox{$\MnOne = 60\GeV$} and $\MnDark = 1\GeV$. The limits are presented in the plane of $\varepsilon$ and $\MgammaDark$. The color gradient represents different branching fraction assumptions for \mbox{$\mathcal{B}(\Ph \to 2\nOne \to 2\gammaDark + 2\nDark)$}, ranging from dark orange (0.05\%) to light orange (10\%). The degradation of the limit around 1\GeV is attributed to the drop in the dimuon branching fraction $\mathcal{B}(\gammaDark \to 2\mu)$ due to the dimuon resonance of the $\phi$ meson~\cite{Batell:2009yf}

Transverse momentum spectrum of charged hadrons measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes selected with the V0M estimator at forward rapidity (top figure, upper panel)

Transverse momentum spectrum of $\pi^{+} + \pi^{-}$ measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

Transverse momentum spectrum of $K^{+} + K^{-}$ measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

Transverse momentum spectrum of $p + \overline{p}$ measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

Transverse momentum spectrum of charged hadrons measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

The $p_{\rm T}$-differential $(K^{+} + K^{-})/(\pi^{+}+\pi^{-})$ particle ratio for two extremes of flattenicity event classes measured in multiplicity-integrated events in pp collisions at $\sqrt{s}$ = 13 TeV

The $p_{\rm T}$-differential $(p + \overline{p})$/($\pi^{+}+\pi^{-}$) particle ratio for two extremes of flattenicity event classes measured in multiplicity-integrated events in pp collisions at $\sqrt{s}$ = 13 TeV

The $p_{\rm T}$-differential $(K^{+} + K^{-})/(\pi^{+}+\pi^{-})$ particle ratio for two extremes of flattenicity event classes measured in the 0-1% V0M event class in pp collisions at $\sqrt{s}$ = 13 TeV

The $p_{\rm T}$-differential $(p + \overline{p})$/($\pi^{+}+\pi^{-}$) particle ratio for two extremes of flattenicity event classes measured in the 0-1% V0M event class in pp collisions at $\sqrt{s}$ = 13 TeV

Transverse momentum-integrated $(K^{+} + K^{-})/(\pi^{+}+\pi^{-})$ particle ratio as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Transverse momentum-integrated $(p + \overline{p})$/($\pi^{+}+\pi^{-}$) particle ratio as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Average transverse momentum of $\pi^{+}+\pi^{-}$ as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Average transverse momentum of $K^{+} + K^{-}$ as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Average transverse momentum of $p + \overline{p}$ as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

'200 GeV, Centrality: 80-100%'

'54.4 GeV, Centrality: 40-60%'

'54.4 GeV, Centrality: 60-80%'

'54.4 GeV, Centrality: 80-100%'

'200 GeV, Centrality: 80-100%'

'$0.4 < M_{ee} < 0.76 GeV/c^{2}$, Centrality: 40-60%'

'$0.4 < M_{ee} < 0.76 GeV/c^{2}$, Centrality: 60-80%'

'$0.4 < M_{ee} < 0.76 GeV/c^{2}$, Centrality: 80-100%'

'$0.76 < M_{ee} < 1.2 GeV/c^{2}$, Centrality: 40-60%'

'$0.76 < M_{ee} < 1.2 GeV/c^{2}$, Centrality: 60-80%'

'$0.76 < M_{ee} < 1.2 GeV/c^{2}$, Centrality: 80-100%'

'$1.2 < M_{ee} < 2.6 GeV/c^{2}$, Centrality: 40-60%'

'$1.2 < M_{ee} < 2.6 GeV/c^{2}$, Centrality: 60-80%'

'$1.2 < M_{ee} < 2.6 GeV/c^{2}$, Centrality: 80-100%'

'54.4 GeV, Centrality: 40-60%'

'54.4 GeV, Centrality: 60-80%'

'54.4 GeV, Centrality: 80-100%'

'200 GeV, Centrality: 80-100%'

'Centrality: 40-60%'

'Centrality: 60-80%'

'Centrality: 80-100%'

'54.4 GeV, Centrality: 40-60%'

'54.4 GeV, Centrality: 60-80%'

'54.4 GeV, Centrality: 80-100%'

'200 GeV, Centrality: 80-100%'

'$-2<cos(4\Delta\phi)>$ as a function of $p_{T}$ @ 200 GeV, Centrality: 80-100%'

'$-2<cos(4\Delta\phi)>$ as a function of $p_{T}$ @ 200 GeV, UPC'

$\chi_{c0}(4500) \to J/\psi(\to \mu^+ \mu^-)\phi(\to K^+ K^-)$ diffractive production cross-section, multiplied by $\mathcal{B}(J/\psi \to \mu^+ \mu^-)$ and $\mathcal{B}(\phi \to K^+ K^-)$ branching ratios.

$\chi_{c1}(4685) \ \mathrm{or} \ \chi_{c0}(4700) \to J/\psi(\to \mu^+ \mu^-)\phi(\to K^+ K^-)$ diffractive production cross-section, multiplied by $\mathcal{B}(J/\psi \to \mu^+ \mu^-)$ and $\mathcal{B}(\phi \to K^+ K^-)$ branching ratios.

Nonresonant $J/\psi(\to \mu^+ \mu^-)\phi(\to K^+ K^-)$ diffractive production cross-section, multiplied by $\mathcal{B}(J/\psi \to \mu^+ \mu^-)$ and $\mathcal{B}(\phi \to K^+ K^-)$ branching ratios.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the narrow spin-0 radion model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the narrow width spin-2 bulk graviton model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the 10%-width spin-0 radion model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the 10%-width spin-2 bulk graviton model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.