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.

The groomed momentum balance distribution of b jets.

The jet fragmentation function distribution of b jets.

The b jet to inclusive jet ratio of the groomed jet radius distributions.

The b jet to inclusive jet ratio of the groomed momentum balance distributions.

v_{3}{4}/v_{3}{2} as a function of centrality in XeXe and PbPb collisions and ratios of XeXe/PbPb.

Normalized symmetric cumulant NSC(2,3) as a function of centrality in XeXe and PbPb collisions and ratios of XeXe/PbPb.

Normalized symmetric cumulant NSC(2,4) as a function of centrality in XeXe and PbPb collisions and ratios of XeXe/PbPb.

Three-harmonic symmetric cumulant SC(2,3,4) and normalized symmetric cumulant NSC(2,3,4) as a function of centrality in XeXe and PbPb collisions.

Three-harmonic symmetric cumulant SC(2,3,5) as a function of centrality in XeXe and PbPb collisions.

Three-harmonic normalized symmetric cumulant NSC(2,3,5) as a function of centrality in XeXe collisions.

Three-harmonic normalized symmetric cumulant NSC(2,3,5) as a function of centrality in PbPb collisions.

Three-harmonic symmetric cumulant SC(3,4,5) as a function of centrality in XeXe and PbPb collisions.

Three-harmonic symmetric cumulant SC(2,4,6) as a function of centrality in XeXe and PbPb collisions.

Six-particle normalized mixed harmonic cumulant nMHC(v_{2}^{4}, v_{3}^{2} as a function of centrality in XeXe and PbPb collisions.

Six-particle normalized mixed harmonic cumulant nMHC(v_{2}^{2}, v_{3}^{4} as a function of centrality in XeXe and PbPb collisions.

Six-particle normalized mixed harmonic cumulant nMHC(v_{2}^{4}, v_{4}^{2} as a function of centrality in XeXe and PbPb collisions.

Six-particle normalized mixed harmonic cumulant nMHC(v_{2}^{2}, v_{4}^{4} as a function of centrality in XeXe and PbPb collisions.

Eight-particle normalized mixed harmonic cumulant nMHC(v_{2}^{6}, v_{3}^{2} as a function of centrality in XeXe and PbPb collisions.

Eight-particle normalized mixed harmonic cumulant nMHC(v_{2}^{4}, v_{3}^{4} as a function of centrality in XeXe and PbPb collisions.

Eight-particle normalized mixed harmonic cumulant nMHC(v_{2}^{2}, v_{3}^{6} as a function of centrality in XeXe and PbPb collisions.

Eight-particle normalized mixed harmonic cumulant nMHC(v_{2}^{6}, v_{4}^{2} as a function of centrality in XeXe and PbPb collisions.

Eight-particle normalized mixed harmonic cumulant nMHC(v_{2}^{4}, v_{4}^{4} as a function of centrality in XeXe and PbPb collisions.

Eight-particle normalized mixed harmonic cumulant nMHC(v_{2}^{2}, v_{4}^{6} as a function of centrality in XeXe and PbPb collisions.

Two-particle correlations v_{4}{2} as a function of centrality in XeXe and PbPb collisions and ratios of XeXe/PbPb.

Four-particle cumulants v_{2}{4} as a function of centrality in XeXe and PbPb collisions.

Four-particle cumulants v_{3}{4} as a function of centrality in XeXe and PbPb collisions.

Observed differential dijet spectrum using full Run-2 data.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for quark-quark type dijet resonances.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for quark-gluon type dijet resonances.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for gluon-gluon type dijet resonances.

The observed and expected 95% CL upper limits on the universal quark coupling $g_{q}'$ as a function of resonance mass for a leptophobic Z' resonance that only couples to quarks.

Distributions of DNN score for the data and post-fit backgrounds (stacked histograms), in the SRs of the ZV channel for the b tag (left) and the b veto (right) channels, for the resolved (merged) category in the first (second) row. The post-fit VBS EW ZV signal is shown overlaid as a red solid line. The overflow is included in the last bin. The lower panels show the ratios of the data to the pre-fit background prediction and post-fit background yield as red open squares and blue points, respectively. The gray band in the lower panels indicates the systematic component of the post-fit background uncertainty. The vertical bars on the data points represent statistical uncertainties. The last bin includes overflow.

Distributions of M(ZV) for the data and post-fit backgrounds (stacked histograms), in the SRs of the merged b tag (left) and the merged b veto (right) channels. The template for one signal hypothesis is shown overlaid as a red solid line. The overflow is included in the last bin. The lower panels show the ratios of the data to the pre-fit background prediction and post-fit background yield as red open squares and blue points, respectively. The gray band in the lower panels indicates the systematic component of the post-fit background uncertainty. The vertical bars on the data points represent statistical uncertainties. The last bin includes overflow.

Distributions of M(ZV) for the data and post-fit backgrounds (stacked histograms), in the SRs of the merged b tag (left) and the merged b veto (right) channels. The template for one signal hypothesis is shown overlaid as a red solid line. The overflow is included in the last bin. The lower panels show the ratios of the data to the pre-fit background prediction and post-fit background yield as red open squares and blue points, respectively. The gray band in the lower panels indicates the systematic component of the post-fit background uncertainty. The vertical bars on the data points represent statistical uncertainties. The last bin includes overflow.

Distributions of M(WV) for the data and post-fit backgrounds (stacked histograms), in the SRs of the electron (left) and muon (right) channels in the merged regime. The template for one signal hypothesis is shown overlaid as a red solid line. The overflow is included in the last bin. The lower panels show the ratios of the data to the pre-fit background prediction and post-fit background yield as red open squares and blue points, respectively. The gray band in the lower panels indicates the systematic component of the post-fit background uncertainty. The vertical bars on the data points represent statistical uncertainties. The last bin includes overflow.

Distributions of M(WV) for the data and post-fit backgrounds (stacked histograms), in the SRs of the electron (left) and muon (right) channels in the merged regime. The template for one signal hypothesis is shown overlaid as a red solid line. The overflow is included in the last bin. The lower panels show the ratios of the data to the pre-fit background prediction and post-fit background yield as red open squares and blue points, respectively. The gray band in the lower panels indicates the systematic component of the post-fit background uncertainty. The vertical bars on the data points represent statistical uncertainties. The last bin includes overflow.

Observed and expected 95% CL intervals on the parameters of quartic operators in the WV and ZV channels, and their combination.

Measured and expected signal strength μ for EW ZV production.

Observed (black markers) and expected distributions of number of jets after event selection. The hatched (grey) areas denote the systematic (total) uncertainties in the expected yields. Events from all data-taking periods and all channels are combined. The lower panel of each plot shows the ratio between observed and expected distributions. The last bin includes all events above the plotted range. The entry Background corresponds to the sum of all the SM predictions.

Observed (black markers) and expected distributions of number of b tagged jets after event selection. The hatched (grey) areas denote the systematic (total) uncertainties in the expected yields. Events from all data-taking periods and all channels are combined. The lower panel of each plot shows the ratio between observed and expected distributions. The last bin includes all events above the plotted range. The entry Background corresponds to the sum of all the SM predictions.

Difference between $p_{\mathrm{T,gen.}}^\text{miss}$ and $p_{\mathrm{T,rec.}}^\text{miss}$ as a function of the $p_{\mathrm{T,gen.}}^\text{miss}$ for signal events. The mean difference between the generated and the $p_{\mathrm{T,rec.}}^\text{miss}$ per bin is shown as solid line. The results for $p_{\mathrm{T}}^\text{miss}$ corrected by the DNN regression (light blue), derived with the PUPPI algorithm (red), and the PF algorithm (orange) are shown. A simulated sample for the 2018 data-taking period is used.

Difference between $p_{\mathrm{T,gen.}}^\text{miss}$ and $p_{\mathrm{T,rec.}}^\text{miss}$ as a function of the $p_{\mathrm{T,gen.}}^\text{miss}$ for signal events. The dashed line shows the standard deviation $\sigma$, which corresponds to the resolution of the $p_{\mathrm{T,gen.}}^\text{miss}$ and $p_{\mathrm{T,rec.}}^\text{miss}$ difference. The results for $p_{\mathrm{T}}^\text{miss}$ corrected by the DNN regression (light blue), derived with the PUPPI algorithm (red), and the PF algorithm (orange) are shown. A simulated sample for the 2018 data-taking period is used.

Difference between $p_{\mathrm{T,gen.}}^\text{miss}$ and $p_{\mathrm{T,rec.}}^\text{miss}$ as a function of the number of primary vertices for signal events. The mean difference between the generated and the $p_{\mathrm{T,rec.}}^\text{miss}$ per bin is shown as solid line. The results for $p_{\mathrm{T}}^\text{miss}$ corrected by the DNN regression (light blue), derived with the PUPPI algorithm (red), and the PF algorithm (orange) are shown. A simulated sample for the 2018 data-taking period is used.

Difference between $p_{\mathrm{T,gen.}}^\text{miss}$ and $p_{\mathrm{T,rec.}}^\text{miss}$ as a function of the $p_{\mathrm{T,gen.}}^\text{miss}$ for signal events. The dashed line shows the standard deviation $\sigma$, which corresponds to the resolution of the $p_{\mathrm{T,gen.}}^\text{miss}$ and $p_{\mathrm{T,rec.}}^\text{miss}$ difference. The results for $p_{\mathrm{T}}^\text{miss}$ corrected by the DNN regression (light blue), derived with the PUPPI algorithm (red), and the PF algorithm (orange) are shown. A simulated sample for the 2018 data-taking period is used.

Results of the $\chi^{2}$ tests for the absolute and normalized differential cross section measurements for each of the predictions. The $\chi^{2}$ values including uncertainties on the predictions are given in parentheses in the table in the paper and in a separate column here.

The p-value of the $\chi^{2}$ tests for the absolute and normalized cross section measurements for each of the predictions. The resulting p-values including uncertainties on the predictions are given in parentheses in the table in the paper and in a separate column here.

Breakdown of the relative contribution from experimental uncertainties in the differential cross section measurement as a function of $p_{\text{T}}^{\nu\nu}$. The statistical uncertainty (dark grey) takes the size of the available simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

Breakdown of the relative contribution from theory uncertainties in the differential cross section measurement as a function of $p_{\text{T}}^{\nu\nu}$. The statistical uncertainty (dark grey) takes the size of the available simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

Breakdown of the relative contribution from experimental uncertainties in the differential cross section measurement as a function of $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$. The statistical uncertainty (dark grey) takes the size of the available simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

Breakdown of the relative contribution from theory uncertainties in the differential cross section measurement as a function of $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$. The statistical uncertainty (dark grey) takes the size of the available simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

Breakdown of the relative contribution from experimental uncertainties in the differential cross section measurement as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0<$\phi$<0.56. Each group of four bins corresponds to $p_{\text{T}}^{\nu\nu}$ bins, and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of the available simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

Breakdown of the relative contribution from experimental uncertainties in the differential cross section measurement as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0.56<$\phi$<1.08. Each group of four bins corresponds to $p_{\text{T}}^{\nu\nu}$ bins, and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of the available simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

Breakdown of the relative contribution from experimental uncertainties in the differential cross section measurement as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 1.08<$\phi$<3.14. Each group of four bins corresponds to $p_{\text{T}}^{\nu\nu}$ bins, and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of the available simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

Breakdown of the relative contribution from theory uncertainties in the differential cross section measurement as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0<$\phi$<0.56. Each group of four bins corresponds to $p_{\text{T}}^{\nu\nu}$ bins, and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of the available simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

Breakdown of the relative contribution from theory uncertainties in the differential cross section measurement as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0.56<$\phi$<1.08. Each group of four bins corresponds to $p_{\text{T}}^{\nu\nu}$ bins, and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of the available simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

Breakdown of the relative contribution from theory uncertainties in the differential cross section measurement as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 1.08<$\phi$<3.14. Each group of four bins corresponds to $p_{\text{T}}^{\nu\nu}$ bins, and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of the available simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

The observed (black markers) and simulated distributions of $p_{\mathrm{T,DNN}}^\text{miss}$ is shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty in the simulation. The lower panel of each plot shows the ratio of observed data to the expectation from MC simulation in each bin. The last bin includes all events above the plotted range.

The observed (black markers) and simulated distributions of $\text{min}[\Delta\phi(\vec{p}_{\text{T, DNN}}^{\text{miss}},\vec{p}_{\text{T}}^{\ell})]$ is shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty in the simulation. The lower panel of each plot shows the ratio of observed data to the expectation from MC simulation in each bin. The last bin includes all events above the plotted range.

The observed (black markers) and simulated two-dimensional distribution of $p_{\mathrm{T,DNN}}^\text{miss}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T, DNN}}^{\text{miss}},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0<$\phi$<0.28 is shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty in the simulation. The lower panel of each plot shows the ratio of observed data to the expectation from MC simulation in each bin. The last bin includes all events above the plotted range.

The observed (black markers) and simulated two-dimensional distribution of $p_{\mathrm{T,DNN}}^\text{miss}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T, DNN}}^{\text{miss}},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0.28<$\phi$<0.56 is shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty in the simulation. The lower panel of each plot shows the ratio of observed data to the expectation from MC simulation in each bin. The last bin includes all events above the plotted range.

The observed (black markers) and simulated two-dimensional distribution of $p_{\mathrm{T,DNN}}^\text{miss}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T, DNN}}^{\text{miss}},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0.56<$\phi$<0.82 is shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty in the simulation. The lower panel of each plot shows the ratio of observed data to the expectation from MC simulation in each bin. The last bin includes all events above the plotted range.

The observed (black markers) and simulated two-dimensional distribution of $p_{\mathrm{T,DNN}}^\text{miss}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T, DNN}}^{\text{miss}},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0.82<$\phi$<1.08 is shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty in the simulation. The lower panel of each plot shows the ratio of observed data to the expectation from MC simulation in each bin. The last bin includes all events above the plotted range.

The observed (black markers) and simulated two-dimensional distribution of $p_{\mathrm{T,DNN}}^\text{miss}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T, DNN}}^{\text{miss}},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 1.08<$\phi$<2.14 is shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty in the simulation. The lower panel of each plot shows the ratio of observed data to the expectation from MC simulation in each bin. The last bin includes all events above the plotted range.

The observed (black markers) and simulated two-dimensional distribution of $p_{\mathrm{T,DNN}}^\text{miss}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T, DNN}}^{\text{miss}},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 2.14<$\phi$<3.2 is shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty in the simulation. The lower panel of each plot shows the ratio of observed data to the expectation from MC simulation in each bin. The last bin includes all events above the plotted range.

Result of the closure test based on simulation accounting for potential BSM contributions based on a top squark pair production scenario with a stop mass of 525 GeV and a neutralino mass of 350 GeV. The potential BSM contribution is scaled by a factor of ten. The test is performed for the $p_{\text{T}}^{\nu\nu}$ using the nominal (black), the regularized (orange), and the bin-by-bin unfolding (purple), based on all data-taking periods combined. The unfolded distributions are compared to the expected distribution (red) based on the sum of the nominal dileptonic $t\bar{t}$ signal sample and BSM signal with the corresponding $\chi^2/\text{ndf}$ values given in the legend in square brackets. The distribution used for the response matrix, is shown in blue. The error bars correspond to the statistical uncertainty. The lower part shows the ratios between the unfolded and the expected distributions.

Result of the closure test based on simulation. The potential BSM contribution is scaled by a factor of ten. The test is performed for the 2D measurement of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ in the azimuthal angle range 0<$\phi$<0.56 using the nominal (black), the regularized (orange), and the bin-by-bin unfolding (purple), based on all data-taking periods combined. The unfolded distributions are compared to the expected distribution (red) based on the sum of the nominal dileptonic $t\bar{t}$ signal sample and BSM signal with the corresponding $\chi^2/\text{ndf}$ values given in the legend. The nominal distribution, used for the response matrix, is shown in blue. The error bars correspond to the statistical uncertainty. The lower part shows the ratios between the unfolded and the expected distributions.

Result of the closure test based on simulation. The potential BSM contribution is scaled by a factor of ten. The test is performed for the 2D measurement of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ in the azimuthal angle range 0.56<$\phi$<1.08 using the nominal (black), the regularized (orange), and the bin-by-bin unfolding (purple), based on all data-taking periods combined. The unfolded distributions are compared to the expected distribution (red) based on the sum of the nominal dileptonic $t\bar{t}$ signal sample and BSM signal with the corresponding $\chi^2/\text{ndf}$ values given in the legend. The nominal distribution, used for the response matrix, is shown in blue. The error bars correspond to the statistical uncertainty. The lower part shows the ratios between the unfolded and the expected distributions.

Result of the closure test based on simulation. The potential BSM contribution is scaled by a factor of ten. The test is performed for the 2D measurement of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ in the azimuthal angle range 1.08<$\phi$<3.14 using the nominal (black), the regularized (orange), and the bin-by-bin unfolding (purple), based on all data-taking periods combined. The unfolded distributions are compared to the expected distribution (red) based on the sum of the nominal dileptonic $t\bar{t}$ signal sample and BSM signal with the corresponding $\chi^2/\text{ndf}$ values given in the legend. The nominal distribution, used for the response matrix, is shown in blue. The error bars correspond to the statistical uncertainty. The lower part shows the ratios between the unfolded and the expected distributions.

The measured differential signal cross section (black markers) as a function of $p_{\text{T}}^{\nu\nu}$ is shown. The theoretical predictions from POWHEG+PYTHIA (dark red), POWHEG+HERWIG (orange), MC@NLO+PYTHIA (purple), and the fixed-order NLO (light blue) and NNLO (brown) calculations are compared to the measurement. The total (statistical) uncertainty in the measurement is shown as an orange (dark grey) band. The lower panel of each plot shows the ratio between theoretical predictions and the measurement.

The measured differential signal cross section (black markers) as a function of $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ is shown. The theoretical predictions from POWHEG+PYTHIA (dark red), POWHEG+HERWIG (orange), MC@NLO+PYTHIA (purple), and the fixed-order NLO (light blue) and NNLO (brown) calculations are compared to the measurement. The total (statistical) uncertainty in the measurement is shown as an orange (dark grey) band. The lower panel of each plot shows the ratio between theoretical predictions and the measurement.

The measured two-dimensional differential signal cross section (black markers) as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0<$\phi$<0.56 is shown. The theoretical predictions from POWHEG+PYTHIA (dark red), POWHEG+HERWIG (orange), MC@NLO+PYTHIA (purple), and the fixed-order NLO (light blue) and NNLO (brown) calculations are compared to the measurement. The total (statistical) uncertainty in the measurement is shown as an orange (dark grey) band. The lower panel of each plot shows the ratio between theoretical predictions and the measurement.

The measured two-dimensional differential signal cross section (black markers) as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0.56<$\phi$<1.08 is shown. The theoretical predictions from POWHEG+PYTHIA (dark red), POWHEG+HERWIG (orange), MC@NLO+PYTHIA (purple), and the fixed-order NLO (light blue) and NNLO (brown) calculations are compared to the measurement. The total (statistical) uncertainty in the measurement is shown as an orange (dark grey) band. The lower panel of each plot shows the ratio between theoretical predictions and the measurement.

The measured two-dimensional differential signal cross section (black markers) as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 1.08<$\phi$<3.14 is shown. The theoretical predictions from POWHEG+PYTHIA (dark red), POWHEG+HERWIG (orange), MC@NLO+PYTHIA (purple), and the fixed-order NLO (light blue) and NNLO (brown) calculations are compared to the measurement. The total (statistical) uncertainty in the measurement is shown as an orange (dark grey) band. The lower panel of each plot shows the ratio between theoretical predictions and the measurement.

The normalized measured differential signal cross section (black markers) as a function of $p_{\text{T}}^{\nu\nu}$ is shown. The theoretical predictions from POWHEG+PYTHIA (dark red), POWHEG+HERWIG (orange), MC@NLO+PYTHIA (purple), and the fixed-order NLO (light blue) and NNLO (brown) calculations are compared to the measurement. The total (statistical) uncertainty in the measurement is shown as an orange (dark grey) band. The lower panel of each plot shows the ratio between theoretical predictions and the measurement.

The normalized measured differential signal cross section (black markers) as a function of $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ is shown. The theoretical predictions from POWHEG+PYTHIA (dark red), POWHEG+HERWIG (orange), MC@NLO+PYTHIA (purple), and the fixed-order NLO (light blue) and NNLO (brown) calculations are compared to the measurement. The total (statistical) uncertainty in the measurement is shown as an orange (dark grey) band. The lower panel of each plot shows the ratio between theoretical predictions and the measurement.

The normalized measured two-dimensional differential signal cross section (black markers) as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0<$\phi$<0.56 is shown. The theoretical predictions from POWHEG+PYTHIA (dark red), POWHEG+HERWIG (orange), MC@NLO+PYTHIA (purple), and the fixed-order NLO (light blue) and NNLO (brown) calculations are compared to the measurement. The total (statistical) uncertainty in the measurement is shown as an orange (dark grey) band. The lower panel of each plot shows the ratio between theoretical predictions and the measurement.

The normalized measured two-dimensional differential signal cross section (black markers) as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 0.56<$\phi$<1.08 is shown. The theoretical predictions from POWHEG+PYTHIA (dark red), POWHEG+HERWIG (orange), MC@NLO+PYTHIA (purple), and the fixed-order NLO (light blue) and NNLO (brown) calculations are compared to the measurement. The total (statistical) uncertainty in the measurement is shown as an orange (dark grey) band. The lower panel of each plot shows the ratio between theoretical predictions and the measurement.

The normalized measured two-dimensional differential signal cross section (black markers) as a function of $p_{\text{T}}^{\nu\nu}$ and $\text{min}[\Delta\phi(\vec{p}_{\text{T}}^{\nu\nu},\vec{p}_{\text{T}}^{\ell})]$ for the azimuthal angle range 1.08<$\phi$<3.14 is shown. The theoretical predictions from POWHEG+PYTHIA (dark red), POWHEG+HERWIG (orange), MC@NLO+PYTHIA (purple), and the fixed-order NLO (light blue) and NNLO (brown) calculations are compared to the measurement. The total (statistical) uncertainty in the measurement is shown as an orange (dark grey) band. The lower panel of each plot shows the ratio between theoretical predictions and the measurement.

The d$^{2}\sigma$/dydp$_{\mathrm{T}}$ production cross section of D$^{0}$ for $8 < p_{\mathrm{T}} < 12$ GeV in ultraperipheral PbPb collisions.

Observed upper limits at 95% CL on the product of the cross section and branching fraction squared, $\sigma (\tilde{\mathrm{t}} \, \bar{\tilde{\mathrm{t}}}) \, \mathcal{B}^{2} ( \tilde{\mathrm{t}} \rightarrow \mathrm{c} \tilde{\chi}_1^0 )$ (left panel), are shown using the color scale where the $\tilde{\mathrm{t}}$ mass is on the $x$-axis and the mass difference between the $\tilde{\mathrm{t}}$ and the LSP is on the $y$-axis. The expected lower mass limits (magenta line) together with their $\pm 1\sigma$ uncertainties (magenta dashed lines) and the observed lower mass limits (black line) are indicated for 100% branching fractions.

Chargino-neutralino production. Observed upper limits at 95% CL on the product of the cross section and the two branching fractions, $\sigma (\tilde{\chi}_1^\pm \tilde{\chi}_2^0) \, \mathcal{B} ( \tilde{\chi}_1^\pm \rightarrow \mathrm{W}^{\pm} \tilde{\chi}_1^0 ) \, \mathcal{B} ( \tilde{\chi}_2^0 \rightarrow \mathrm{Z} \tilde{\chi}_1^0 )$, are shown using the color scale where the $\tilde{\chi}_1^\pm/\tilde{\chi}_2^0$ mass is on the $x$-axis and the mass difference between the $\tilde{\chi}_1^\pm/\tilde{\chi}_2^0$ and the LSP is on the $y$-axis. For these results, based on the TChiWZ simplified model, the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ masses are set equal. The expected lower mass limits (magenta line) together with their $\pm 1\sigma$ uncertainties (magenta dashed lines) and the observed lower mass limits (black line) are indicated for 100% branching fractions for Wino-like cross-sections (left panel) and for Higgsino-like cross-sections (right panel).

Chargino-neutralino production. Observed upper limits at 95% CL on the product of the cross section and the two branching fractions, $\sigma (\tilde{\chi}_1^\pm \tilde{\chi}_2^0) \, \mathcal{B} ( \tilde{\chi}_1^\pm \rightarrow \mathrm{W}^{\pm} \tilde{\chi}_1^0 ) \, \mathcal{B} ( \tilde{\chi}_2^0 \rightarrow \mathrm{Z} \tilde{\chi}_1^0 )$, are shown using the color scale where the $\tilde{\chi}_1^\pm/\tilde{\chi}_2^0$ mass is on the $x$-axis and the mass difference between the $\tilde{\chi}_1^\pm/\tilde{\chi}_2^0$ and the LSP is on the $y$-axis. For these results, based on the TChiWZ simplified model, the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ masses are set equal. The expected lower mass limits (magenta line) together with their $\pm 1\sigma$ uncertainties (magenta dashed lines) and the observed lower mass limits (black line) are indicated for 100% branching fractions for Wino-like cross-sections (left panel) and for Higgsino-like cross-sections (right panel).

Chargino pair production. The left panel shows the observed upper limits at 95% CL on the product of the cross section and the branching fraction squared, $\sigma (\tilde{\chi}_1^+ \tilde{\chi}_1^- ) \, \mathcal{B}^{2} ( \tilde{\chi}_1^\pm \rightarrow \mathrm{W}^{\pm} \tilde{\chi}_1^0 ) $ are shown using the color scale where the $\tilde{\chi}_1^\pm$ mass is on the $x$-axis and the mass difference between the $\tilde{\chi}_1^\pm $ and the LSP is on the $y$-axis. The expected lower mass limits (magenta line) together with their $\pm 1\sigma$ uncertainties (magenta dashed lines) and the observed lower mass limits (black line) are indicated for 100% branching fractions for Wino-like cross-sections. The right panel shows the results for chargino pair production with decays as in the TChiSlepSnu model with democratic decay via an intermediate sneutrino or charged slepton ($\widetilde{\ell}^{\pm}_{\mathrm{L}}$) with mass halfway between the chargino and the lightest neutralino. These model predictions also assume Wino-like cross sections.

Chargino pair production. The left panel shows the observed upper limits at 95% CL on the product of the cross section and the branching fraction squared, $\sigma (\tilde{\chi}_1^+ \tilde{\chi}_1^- ) \, \mathcal{B}^{2} ( \tilde{\chi}_1^\pm \rightarrow \mathrm{W}^{\pm} \tilde{\chi}_1^0 ) $ are shown using the color scale where the $\tilde{\chi}_1^\pm$ mass is on the $x$-axis and the mass difference between the $\tilde{\chi}_1^\pm $ and the LSP is on the $y$-axis. The expected lower mass limits (magenta line) together with their $\pm 1\sigma$ uncertainties (magenta dashed lines) and the observed lower mass limits (black line) are indicated for 100% branching fractions for Wino-like cross-sections. The right panel shows the results for chargino pair production with decays as in the TChiSlepSnu model with democratic decay via an intermediate sneutrino or charged slepton ($\widetilde{\ell}^{\pm}_{\mathrm{L}}$) with mass halfway between the chargino and the lightest neutralino. These model predictions also assume Wino-like cross sections.

Slepton pair production. Observed 95% CL upper limits on the product of the cross section and branching fraction squared for direct slepton pair production followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). Slepton $\tilde{l}_{\mathrm{L}/\mathrm{R}}$ indicates the scalar supersymmetric partner of left- and right-handed electrons and muons. The limit is shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. The median expected exclusion regions for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the product of the cross section and branching fraction squared for direct slepton pair production followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). Slepton $\tilde{l}_{\mathrm{L}/\mathrm{R}}$ indicates the scalar supersymmetric partner of left- and right-handed electrons and muons. The limit is shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. The median expected exclusion regions for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the product of the cross section and branching fraction squared for direct slepton pair production followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). Slepton $\tilde{l}_{\mathrm{L}/\mathrm{R}}$ indicates the scalar supersymmetric partner of left- and right-handed electrons and muons. The limit is shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. The median expected exclusion regions for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

Parametric models of the signal process for mX=1000GeV, mY=800GeV in their most sensitive SR. The histograms are normalized to unity. The acronym 'dof' stands for the numbers of degrees of freedom of the parametric model. The signal is modeled using a double-sided Crystal Ball (DCB) function defined as: DCB$(x)$ = \[ \begin{cases} N \cdot A_1 \cdot (B_1 - x_s)^{-m_1}, & x_s \leq -\beta_1 \\ N \cdot e^{-\frac{1}{2} x_s^2}, & -\beta_1 < x_s < \beta_2 \\ N \cdot A_2 \cdot (B_2 + x_s)^{-m_2}, & x_s \geq \beta_2 \end{cases} \] with \(x_s = \frac{x - \mu}{\sigma}\), and: \[ A_1 = \left( \frac{m_1}{\beta_1} \right)^{m_1} e^{-\frac{1}{2} \beta_1^2}, \quad B_1 = \frac{m_1}{\beta_1} - \beta_1 \] \[ A_2 = \left( \frac{m_2}{\beta_2} \right)^{m_2} e^{-\frac{1}{2} \beta_2^2}, \quad B_2 = \frac{m_2}{\beta_2} - \beta_2 \] The DCB parameters for this signal model are: \[ \begin{aligned} N &= 0.35127, & \mu &= 798.50175, & \sigma &= 7.32360 \\ \beta_1 &= 1.50996, & m_1 &= 5.30414, & \beta_2 &= 1.56253, & m_2 &= 18.3172 \end{aligned} \]

Background-only fit (red line) and signal+background fit (blue line) for the mass point hypothesis of mX=280GeV, mY=90GeV. The red line in the lower panel shows the background-only fit, which is by definition zero, and it is added as a visual aid. The most sensitive SR is shown. The choice of the background functional form is determined by the maximum likelihood fit. Category 0 is the most sensitive signal region. Drell-Yan background component is fitted as \[f_{\text{DY}}(m) = \frac{0.530965}{\sqrt{2\pi} \cdot 2.81859} \cdot \exp\left( -\frac{(m - 90.7199)^2}{2 \cdot 2.81859^2} \right)\] In the both the background-only fit and signal plus background fit, the non-resonant component of the pdf is choosing the following function. \[f(m) = \frac{1}{2.6086 \times 10^{-6}} \cdot m^{-3.82346}\] The fitted normalization factor for the function is: 23.5522 Signal fitted r value is 1.21391 in all the category simultaneous fit. The signal is fitted with a double-sided Crystal Ball (DCB) function defined as: DCB$(x)$ = \[ \begin{cases} N \cdot A_1 \cdot (B_1 - x_s)^{-m_1}, & x_s \leq -\beta_1 \\ N \cdot e^{-\frac{1}{2} x_s^2}, & -\beta_1 < x_s < \beta_2 \\ N \cdot A_2 \cdot (B_2 + x_s)^{-m_2}, & x_s \geq \beta_2 \end{cases} \] with \(x_s = \frac{x - \mu}{\sigma}\), and: \[ A_1 = \left( \frac{m_1}{\beta_1} \right)^{m_1} e^{-\frac{1}{2} \beta_1^2}, \quad B_1 = \frac{m_1}{\beta_1} - \beta_1 \] \[ A_2 = \left( \frac{m_2}{\beta_2} \right)^{m_2} e^{-\frac{1}{2} \beta_2^2}, \quad B_2 = \frac{m_2}{\beta_2} - \beta_2 \] The DCB parameters are: 2016: \[ \begin{aligned} N &= 0.6422, & \mu &= 89.9423, & \sigma &= 0.8881 \\ \beta_1 &= 1.3471, & m_1 &= 2.7908, & \beta_2 &= 1.4211, & m_2 &= 19.9999 \end{aligned} \] 2017: \[ \begin{aligned} N &= 0.8047, & \mu &= 89.9948, & \sigma &= 1.03024 \\ \beta_1 &= 1.6449, & m_1 &= 1.9425, & \beta_2 &= 1.7434, & m_2 &= 19.9982 \end{aligned} \] 2018: \[ \begin{aligned} N &= 1.1200, & \mu &= 90.0211, & \sigma &= 0.9284 \\ \beta_1 &= 1.1778, & m_1 &= 4.4273, & \beta_2 &= 1.5851, & m_2 &= 7.1198 \end{aligned} \] The plot is filled with a bin size of 5.84GeV.

Background-only fit (red line) and signal+background fit (blue line) for the mass point hypothesis of mX=280GeV, mY=90GeV. The red line in the lower panel shows the background-only fit, which is by definition zero, and it is added as a visual aid. The second most sensitive SR is shown. The choice of the background functional form is determined by the maximum likelihood fit. Category 1 is the second most sensitive signal region. Drell-Yan background component is fitted as \[f_{\text{DY}}(m) = \frac{0.532195}{\sqrt{2\pi} \cdot 3.10187} \cdot \exp\left( -\frac{(m - 89.7677)^2}{2 \cdot 3.10187^2} \right)\] In the both the background-only fit and signal plus background fit, the non-resonant component of the pdf is choosing the following function. \[f(m) = \frac{1}{14.0606} \cdot \exp\left(-\left(0.1 \cdot \frac{m}{100} + 1.73 \cdot \left(\frac{m}{100}\right)^2\right)\right)\] The fitted normalization factor for the function is: 21.6386 Signal fitted r value is 1.21391 in all the category simultaneous fit. The signal is fitted with a double-sided Crystal Ball (DCB) function defined as: DCB$(x)$ = \[ \begin{cases} N \cdot A_1 \cdot (B_1 - x_s)^{-m_1}, & x_s \leq -\beta_1 \\ N \cdot e^{-\frac{1}{2} x_s^2}, & -\beta_1 < x_s < \beta_2 \\ N \cdot A_2 \cdot (B_2 + x_s)^{-m_2}, & x_s \geq \beta_2 \end{cases} \] with \(x_s = \frac{x - \mu}{\sigma}\), and: \[ A_1 = \left( \frac{m_1}{\beta_1} \right)^{m_1} e^{-\frac{1}{2} \beta_1^2}, \quad B_1 = \frac{m_1}{\beta_1} - \beta_1 \] \[ A_2 = \left( \frac{m_2}{\beta_2} \right)^{m_2} e^{-\frac{1}{2} \beta_2^2}, \quad B_2 = \frac{m_2}{\beta_2} - \beta_2 \] The DCB parameters are: 2016: \[ \begin{aligned} N &= 0.3784, & \mu &= 89.8969, & \sigma &= 0.9189 \\ \beta_1 &= 1.1920, & m_1 &= 3.7178, & \beta_2 &= 1.3318, & m_2 &= 10.9787 \end{aligned} \] 2017: \[ \begin{aligned} N &= 0.4914, & \mu &= 89.8859, & \sigma &= 1.0728 \\ \beta_1 &= 1.3515, & m_1 &= 3.2846, & \beta_2 &= 1.4371, & m_2 &= 20.0000 \end{aligned} \] 2018: \[ \begin{aligned} N &= 0.6766, & \mu &= 89.9327, & \sigma &= 1.1213 \\ \beta_1 &= 1.9188, & m_1 &= 1.5353, & \beta_2 &= 1.9907, & m_2 &= 12.4664 \end{aligned} \] The plot is filled with a bin size of 5.84GeV.

Observed(expected) upper limits on the $\sigma\mathcal{B}$ for the X to $Y(\gamma\gamma)H(bb)$ signal with the different mass hypotheses are shown with solid(dashed) lines. The green and the yellow bands define the $\pm 68\%$ and $\pm 95\%$ uncertainty bands, respectively. For visualization purposes, the upper limit for mass points with different $m_X$ has been multiplied with a constant factor quoted on the right of each band. Mass points with $m_X$ ranging from 600 to 1000GeV are shown.

Observed(expected) upper limits on the $\sigma\mathcal{B}$ for the X 5p $Y(\gamma\gamma)H(bb)$ signal with the different mass hypotheses are shown with solid(dashed) lines. The green and the yellow bands define the $\pm 68\%$ and $\pm 95\%$ uncertainty bands, respectively. For visualization purposes, the upper limit for mass points with different $m_X$ has been multiplied with a constant factor quoted on the right of each band. Mass points with $m_X$ ranging from 240 to 550GeV are shown.

Observed(expected) upper limits on the $\sigma\mathcal{B}$ for the X to $Y(\gamma\gamma)H(bb)$ signal with the different $m_Y$ hypotheses are shown with solid(dashed) lines. The inner (green) band and the outer (yellow) band indicate the regions containing $68\%$ and $95\%$, respectively, of the distribution of limits expected under the background-only hypothesis. This plot shows the set of mass points with the lowest $m_X = 240GeV$

Observed(expected) upper limits on the $\sigma\mathcal{B}$ for the X to $Y(\gamma\gamma)H(bb)$ signal with the different $m_Y$ hypotheses are shown with solid(dashed) lines. The inner (green) band and the outer (yellow) band indicate the regions containing $68\%$ and $95\%$, respectively, of the distribution of limits expected under the background-only hypothesis. This plot shows the set of mass points with the highest $m_{X} = 1000GeV$