Observed and expected 95% CL upper limits on the cross section for the resonant production of a new spin-2 particle $\mathrm{X}^{(2)}$ which decays to Higgs boson pairs, $\sigma(\mathrm{pp} \to \mathrm{X}^{(2)} \to \mathrm{HH})$, given for different values of $m_\mathrm{X}$ in the range 260-1000 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$ which subsequently decay to a pair of photons and a pair of tau leptons, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH} \to \gamma\gamma\tau\tau)$. The limits are shown for different values of $m_\mathrm{Y}$ in the range 50-800 GeV and at particular values of $m_\mathrm{X}$ in the range 300-1000 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$ which subsequently decay to a pair of photons and a pair of tau leptons, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH} \to \gamma\gamma\tau\tau)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-1000 GeV and at particular values of $m_\mathrm{Y}$ in the range 50-800 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$ which subsequently decay to a pair of photons and a pair of tau leptons, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH} \to \gamma\gamma\tau\tau)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-1000 GeV and $m_\mathrm{Y}$ in the range 50-800 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{Y}$ in the range 70-125 GeV and at particular values of $m_\mathrm{X}$ in the range 300-1000 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-800 GeV and at particular values of $m_\mathrm{Y}$ in the range 70-125 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-1000 GeV and $m_\mathrm{Y}$ in the range 70-125 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{Y}$ in the range 125-800 GeV and at particular values of $m_\mathrm{X}$ in the range 300-1000 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-800 GeV and at particular values of $m_\mathrm{Y}$ in the range 125-800 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-1000 GeV and $m_\mathrm{Y}$ in the range 125-800 GeV.

Comparison of the shapes extrapolated with a linear or quadratic extrapolation of the QCD contribution in the W+ signal region.

Comparison of the shapes extrapolated with a linear or quadratic extrapolation of the QCD contribution in the W- signal region.

Post-fit distributions of the transverse mass in the W+ signal region.

Post-fit distributions of the transverse mass in the W- signal region.

Post-fit distributions of the dimuon mass in the Z boson signal region.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W+ signal region in the transverse mass bin 12 < mT < 18 GeV. The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W- signal region in the transverse mass bin 12 < mT < 18 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W+ signal region in the transverse mass bin 102 < mT < 108 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (larger value) and quadratic (smaller value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W- signal region in the transverse mass bin 102 < mT < 108 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

The table compares theoretical and measured values of the product of the fiducial cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value. The relative uncertainties on the measured values are given separately with and without the contribution from the luminosity uncertainty.

The table compares theoretical and measured values of the product of the total cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value. The relative uncertainties on the measured values are given separately with and without the contribution from the luminosity uncertainty.

The table compares theoretical and measured ratios of the product of the fiducial cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value.

The table compares theoretical and measured ratios of the product of the total cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value.

Comparison of the measured and theoretical cross sections for Z, W+, W-, and W production at different center-of-mass energies. The theoretical predictions are generated with DYTurbo in combination with the NNPDF 3.1 PDF set.

Data vs. the ABCD method background prediction for 2017 in $\mathrm{SR_{b}^{mm}}$. Events are divided by the bin width, hence fractional data counts. Error bars show statistical uncertainties of data. The blue band shows the propagated uncertainty of all individual fit variations in a given bin, which we consider to be uncorrelated. The lower panels show the ratio of the observed data to the background estimation.

Data vs. the ABCD method background prediction for 2018 in $\mathrm{SR_{b}^{mm}}$. Events are divided by the bin width, hence fractional data counts. Error bars show statistical uncertainties of data. The blue band shows the propagated uncertainty of all individual fit variations in a given bin, which we consider to be uncorrelated. The lower panels show the ratio of the observed data to the background estimation.

Data vs. the ABCD method background prediction for 2016 in $\mathrm{SR_{b+\textrm{j}/b}^{mm}}$. Events are divided by the bin width, hence fractional data counts. Error bars show statistical uncertainties of data. The blue band shows the propagated uncertainty of all individual fit variations in a given bin, which we consider to be uncorrelated. The lower panels show the ratio of the observed data to the background estimation.

Data vs. the ABCD method background prediction for 2017 in $\mathrm{SR_{b+\textrm{j}/b}^{mm}}$. Events are divided by the bin width, hence fractional data counts. Error bars show statistical uncertainties of data. The blue band shows the propagated uncertainty of all individual fit variations in a given bin, which we consider to be uncorrelated. The lower panels show the ratio of the observed data to the background estimation.

Data vs. the ABCD method background prediction for 2018 in $\mathrm{SR_{b+\textrm{j}/b}^{mm}}$. Events are divided by the bin width, hence fractional data counts. Error bars show statistical uncertainties of data. The blue band shows the propagated uncertainty of all individual fit variations in a given bin, which we consider to be uncorrelated. The lower panels show the ratio of the observed data to the background estimation.

The asymptotic 95% CL limits on the product of cross section ($\sigma$), branching fraction into a pair of muons ($\mathcal{B}$) to dimuons, acceptance (A), and efficiency ($\epsilon$) in $\mathrm{SR_{b}^{mm}}$. A combination of both limits depends on the relative acceptance contributions in each region and is omitted here.

The asymptotic 95% CL limits on the product of cross section ($\sigma$), branching fraction into a pair of muons ($\mathcal{B}$) to dimuons, acceptance (A), and efficiency ($\epsilon$) in $\mathrm{SR_{b+\textrm{j}/b}^{mm}}$. A combination of both limits depends on the relative acceptance contributions in each region and is omitted here.

Cutflow for signal samples in 2016. Samples are divided by generated mass, signal region, and the b to s quark flavor changing coupling of the $Z^\prime$, dbs. For some mass values, multiple values of this coupling strength have been generated. It is described in this minimal lagrangian: $\mathcal{L} \supset Z^{\prime}_{\eta} \Big[ g_{\mu}\, \bar{\mu}\, \gamma^{\eta}\, \mu + \big( g_b\, \delta_{bs}\, \bar{s}\, \gamma^{\eta}\, P_{L}\, b + \text{h.c.} \big) \Big]$. The cutflow shows the selection into signal regions, and then the effects of specialized HTLT and RelMET selections, which grow with greater mass.

Cutflow for signal samples in 2017. Samples are divided by generated mass, signal region, and the b to s quark flavor changing coupling of the $Z^\prime$, dbs. For some mass values, multiple values of this coupling strength have been generated. It is described in this minimal lagrangian: $\mathcal{L} \supset Z^{\prime}_{\eta} \Big[ g_{\mu}\, \bar{\mu}\, \gamma^{\eta}\, \mu + \big( g_b\, \delta_{bs}\, \bar{s}\, \gamma^{\eta}\, P_{L}\, b + \text{h.c.} \big) \Big]$. The cutflow shows the selection into signal regions, and then the effects of specialized HTLT and RelMET selections, which grow with greater mass.

Cutflow for signal samples in 2018. Samples are divided by generated mass, signal region, and the b to s quark flavor changing coupling of the $Z^\prime$, dbs. For some mass values, multiple values of this coupling strength have been generated. It is described in this minimal lagrangian: $\mathcal{L} \supset Z^{\prime}_{\eta} \Big[ g_{\mu}\, \bar{\mu}\, \gamma^{\eta}\, \mu + \big( g_b\, \delta_{bs}\, \bar{s}\, \gamma^{\eta}\, P_{L}\, b + \text{h.c.} \big) \Big]$. The cutflow shows the selection into signal regions, and then the effects of specialized HTLT and RelMET selections, which grow with greater mass.