CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-TOP-24-001 ; CERN-EP-2025-201
Measurement of the dineutrino system kinematic variables in dileptonic top quark pair production in proton-proton collisions at $ \sqrt{s}= $ 13 TeV
Submitted to J. High Energy Phys.
Abstract: Differential top quark pair production cross sections are measured in the dilepton final states $ \mathrm{e}^+\mathrm{e}^- $, $ \mu^{+}\mu^{-} $, and $ \mathrm{e}^{\pm}\mu^{\mp} $, as a function of kinematic variables of the two-neutrino system: the transverse momentum $ p_{\mathrm{T}}^{\nu\nu} $ of the dineutrino system, the minimum distance in azimuthal angle between $ {\vec p}_{\mathrm{T}}^{\,\nu\nu} $ and leptons, and in two dimensions in bins of both observables. The measurements are performed using CERN LHC proton-proton collisions at $ \sqrt{s}= $ 13 TeV, recorded by the CMS detector between 2016 and 2018, corresponding to an integrated luminosity of 138 fb$ ^{-1} $. The measured cross sections are unfolded to the particle level using an unregularized least squares method. Results are compared with predictions by the standard model of particle physics, and found to be in agreement with theoretical calculations as well as Monte Carlo simulations.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
In the left diagram, the SM $ \mathrm{t} \overline{\mathrm{t}} $ production process is sketched, whereas the right diagram shows the production of a hypothetical top squark pair, $ \tilde{\mathrm{t}}_{1}\overline{\tilde{\mathrm{t}}}_{1} $, with both top squarks decaying to a top quark and a neutralino, $ \tilde{\chi}_{1}^{0} $. This analysis focuses on signatures, where both of the W bosons decay leptonically.

png pdf
Figure 1-a:
In the left diagram, the SM $ \mathrm{t} \overline{\mathrm{t}} $ production process is sketched, whereas the right diagram shows the production of a hypothetical top squark pair, $ \tilde{\mathrm{t}}_{1}\overline{\tilde{\mathrm{t}}}_{1} $, with both top squarks decaying to a top quark and a neutralino, $ \tilde{\chi}_{1}^{0} $. This analysis focuses on signatures, where both of the W bosons decay leptonically.

png pdf
Figure 1-b:
In the left diagram, the SM $ \mathrm{t} \overline{\mathrm{t}} $ production process is sketched, whereas the right diagram shows the production of a hypothetical top squark pair, $ \tilde{\mathrm{t}}_{1}\overline{\tilde{\mathrm{t}}}_{1} $, with both top squarks decaying to a top quark and a neutralino, $ \tilde{\chi}_{1}^{0} $. This analysis focuses on signatures, where both of the W bosons decay leptonically.

png pdf
Figure 2:
Observed (black markers) and expected (filled histograms) distributions of leading lepton $ p_{\mathrm{T}} $ (upper left), leading jet $ p_{\mathrm{T}} $ (upper right), the number of jets (lower left), and the number of b tagged jets (lower right), after event selection. The hatched (grey) areas denote the systematic (total) uncertainties on the expected yields. Vertical and horizontal error bars on the data points represent the statistical uncertainty and the bin width, respectively. 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.

png pdf
Figure 2-a:
Observed (black markers) and expected (filled histograms) distributions of leading lepton $ p_{\mathrm{T}} $ (upper left), leading jet $ p_{\mathrm{T}} $ (upper right), the number of jets (lower left), and the number of b tagged jets (lower right), after event selection. The hatched (grey) areas denote the systematic (total) uncertainties on the expected yields. Vertical and horizontal error bars on the data points represent the statistical uncertainty and the bin width, respectively. 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.

png pdf
Figure 2-b:
Observed (black markers) and expected (filled histograms) distributions of leading lepton $ p_{\mathrm{T}} $ (upper left), leading jet $ p_{\mathrm{T}} $ (upper right), the number of jets (lower left), and the number of b tagged jets (lower right), after event selection. The hatched (grey) areas denote the systematic (total) uncertainties on the expected yields. Vertical and horizontal error bars on the data points represent the statistical uncertainty and the bin width, respectively. 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.

png pdf
Figure 2-c:
Observed (black markers) and expected (filled histograms) distributions of leading lepton $ p_{\mathrm{T}} $ (upper left), leading jet $ p_{\mathrm{T}} $ (upper right), the number of jets (lower left), and the number of b tagged jets (lower right), after event selection. The hatched (grey) areas denote the systematic (total) uncertainties on the expected yields. Vertical and horizontal error bars on the data points represent the statistical uncertainty and the bin width, respectively. 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.

png pdf
Figure 2-d:
Observed (black markers) and expected (filled histograms) distributions of leading lepton $ p_{\mathrm{T}} $ (upper left), leading jet $ p_{\mathrm{T}} $ (upper right), the number of jets (lower left), and the number of b tagged jets (lower right), after event selection. The hatched (grey) areas denote the systematic (total) uncertainties on the expected yields. Vertical and horizontal error bars on the data points represent the statistical uncertainty and the bin width, respectively. 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.

png pdf
Figure 3:
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} $ (upper) and the number of primary vertices (lower) for signal events. The mean difference between the generated and the $ p_{\mathrm{T,rec.}}^\text{miss} $ per bin is shown as solid line, while the dashed line shows the standard deviation $ \sigma $, which corresponds to the resolution. The results for $ p_{\mathrm{T}}^\text{miss} $ corrected by the DNN regression (thick light blue), derived with the PUPPI algorithm (thin red), and the PF algorithm (thinner orange) are shown. A simulated sample for the 2018 data-taking period is used.

png pdf
Figure 3-a:
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} $ (upper) and the number of primary vertices (lower) for signal events. The mean difference between the generated and the $ p_{\mathrm{T,rec.}}^\text{miss} $ per bin is shown as solid line, while the dashed line shows the standard deviation $ \sigma $, which corresponds to the resolution. The results for $ p_{\mathrm{T}}^\text{miss} $ corrected by the DNN regression (thick light blue), derived with the PUPPI algorithm (thin red), and the PF algorithm (thinner orange) are shown. A simulated sample for the 2018 data-taking period is used.

png pdf
Figure 3-b:
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} $ (upper) and the number of primary vertices (lower) for signal events. The mean difference between the generated and the $ p_{\mathrm{T,rec.}}^\text{miss} $ per bin is shown as solid line, while the dashed line shows the standard deviation $ \sigma $, which corresponds to the resolution. The results for $ p_{\mathrm{T}}^\text{miss} $ corrected by the DNN regression (thick light blue), derived with the PUPPI algorithm (thin red), and the PF algorithm (thinner orange) are shown. A simulated sample for the 2018 data-taking period is used.

png pdf
Figure 4:
Breakdown of the relative uncertainties from experimental (left) and theoretical (right) uncertainties on the differential cross section measurement as a function of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and both variables (lower). In the last case, each group of four bins corresponds to $ p_{\mathrm{T}}^{\nu\nu} $ bins, and $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of available event counts in both the simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

png pdf
Figure 4-a:
Breakdown of the relative uncertainties from experimental (left) and theoretical (right) uncertainties on the differential cross section measurement as a function of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and both variables (lower). In the last case, each group of four bins corresponds to $ p_{\mathrm{T}}^{\nu\nu} $ bins, and $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of available event counts in both the simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

png pdf
Figure 4-b:
Breakdown of the relative uncertainties from experimental (left) and theoretical (right) uncertainties on the differential cross section measurement as a function of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and both variables (lower). In the last case, each group of four bins corresponds to $ p_{\mathrm{T}}^{\nu\nu} $ bins, and $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of available event counts in both the simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

png pdf
Figure 4-c:
Breakdown of the relative uncertainties from experimental (left) and theoretical (right) uncertainties on the differential cross section measurement as a function of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and both variables (lower). In the last case, each group of four bins corresponds to $ p_{\mathrm{T}}^{\nu\nu} $ bins, and $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of available event counts in both the simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

png pdf
Figure 4-d:
Breakdown of the relative uncertainties from experimental (left) and theoretical (right) uncertainties on the differential cross section measurement as a function of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and both variables (lower). In the last case, each group of four bins corresponds to $ p_{\mathrm{T}}^{\nu\nu} $ bins, and $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of available event counts in both the simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

png pdf
Figure 4-e:
Breakdown of the relative uncertainties from experimental (left) and theoretical (right) uncertainties on the differential cross section measurement as a function of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and both variables (lower). In the last case, each group of four bins corresponds to $ p_{\mathrm{T}}^{\nu\nu} $ bins, and $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of available event counts in both the simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

png pdf
Figure 4-f:
Breakdown of the relative uncertainties from experimental (left) and theoretical (right) uncertainties on the differential cross section measurement as a function of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and both variables (lower). In the last case, each group of four bins corresponds to $ p_{\mathrm{T}}^{\nu\nu} $ bins, and $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ bin edges are indicated by vertical solid lines. The statistical uncertainty (dark grey) takes the size of available event counts in both the simulated samples and the recorded data into account. The last bin includes all events above the plotted range.

png pdf
Figure 5:
The observed (black markers) and simulated distributions of $ p_{\mathrm{T,DNN}}^\text{miss} $ (upper left), $ \min[\Delta\phi(\vec{p}_{\mathrm{T,DNN}}^{\,\text{miss}},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (upper right), and the 2D distribution (lower) with both observables are shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty on the simulation. Vertical and horizontal error bars on the data points represent the statistical uncertainty and the bin width, respectively. 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.

png pdf
Figure 5-a:
The observed (black markers) and simulated distributions of $ p_{\mathrm{T,DNN}}^\text{miss} $ (upper left), $ \min[\Delta\phi(\vec{p}_{\mathrm{T,DNN}}^{\,\text{miss}},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (upper right), and the 2D distribution (lower) with both observables are shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty on the simulation. Vertical and horizontal error bars on the data points represent the statistical uncertainty and the bin width, respectively. 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.

png pdf
Figure 5-b:
The observed (black markers) and simulated distributions of $ p_{\mathrm{T,DNN}}^\text{miss} $ (upper left), $ \min[\Delta\phi(\vec{p}_{\mathrm{T,DNN}}^{\,\text{miss}},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (upper right), and the 2D distribution (lower) with both observables are shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty on the simulation. Vertical and horizontal error bars on the data points represent the statistical uncertainty and the bin width, respectively. 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.

png pdf
Figure 5-c:
The observed (black markers) and simulated distributions of $ p_{\mathrm{T,DNN}}^\text{miss} $ (upper left), $ \min[\Delta\phi(\vec{p}_{\mathrm{T,DNN}}^{\,\text{miss}},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (upper right), and the 2D distribution (lower) with both observables are shown. Events from all data-taking periods and all channels are combined. The hatched (grey) areas indicate the systematic (total) uncertainty on the simulation. Vertical and horizontal error bars on the data points represent the statistical uncertainty and the bin width, respectively. 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.

png pdf
Figure 6:
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_{\mathrm{T}}^{\nu\nu} $ (upper) and 2D measurement (lower) using the nominal (black circles), the regularized (orange triangles), and the bin-by-bin unfolding (purple squares), based on all data-taking periods combined. The unfolded distributions are compared to the expected distribution (red solid line) based on the sum of the nominal dileptonic $ \mathrm{t} \overline{\mathrm{t}} $ signal sample and BSM signal with the corresponding $ \chi^2 $/(number of degrees of freedom) values given in the legend in square brackets. The distribution used for the response matrix is shown in blue (dotted line). The error bars correspond to the statistical uncertainty. The lower part shows the ratios between the unfolded and the expected distributions.

png pdf
Figure 6-a:
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_{\mathrm{T}}^{\nu\nu} $ (upper) and 2D measurement (lower) using the nominal (black circles), the regularized (orange triangles), and the bin-by-bin unfolding (purple squares), based on all data-taking periods combined. The unfolded distributions are compared to the expected distribution (red solid line) based on the sum of the nominal dileptonic $ \mathrm{t} \overline{\mathrm{t}} $ signal sample and BSM signal with the corresponding $ \chi^2 $/(number of degrees of freedom) values given in the legend in square brackets. The distribution used for the response matrix is shown in blue (dotted line). The error bars correspond to the statistical uncertainty. The lower part shows the ratios between the unfolded and the expected distributions.

png pdf
Figure 6-b:
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_{\mathrm{T}}^{\nu\nu} $ (upper) and 2D measurement (lower) using the nominal (black circles), the regularized (orange triangles), and the bin-by-bin unfolding (purple squares), based on all data-taking periods combined. The unfolded distributions are compared to the expected distribution (red solid line) based on the sum of the nominal dileptonic $ \mathrm{t} \overline{\mathrm{t}} $ signal sample and BSM signal with the corresponding $ \chi^2 $/(number of degrees of freedom) values given in the legend in square brackets. The distribution used for the response matrix is shown in blue (dotted line). The error bars correspond to the statistical uncertainty. The lower part shows the ratios between the unfolded and the expected distributions.

png pdf
Figure 7:
The measured differential signal cross sections (black markers) as functions of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and double-differential in both observables (lower) are 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 cross sections in the last $ p_{\mathrm{T}}^{\nu\nu} $ bin of the upper (lower) plot are determined for $ p_{\mathrm{T}}^{\nu\nu} > $ 410 (200) GeV and normalised to a bin width of 90 (200) GeV for visualisation purposes. The total (statistical) uncertainty on 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.

png pdf
Figure 7-a:
The measured differential signal cross sections (black markers) as functions of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and double-differential in both observables (lower) are 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 cross sections in the last $ p_{\mathrm{T}}^{\nu\nu} $ bin of the upper (lower) plot are determined for $ p_{\mathrm{T}}^{\nu\nu} > $ 410 (200) GeV and normalised to a bin width of 90 (200) GeV for visualisation purposes. The total (statistical) uncertainty on 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.

png pdf
Figure 7-b:
The measured differential signal cross sections (black markers) as functions of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and double-differential in both observables (lower) are 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 cross sections in the last $ p_{\mathrm{T}}^{\nu\nu} $ bin of the upper (lower) plot are determined for $ p_{\mathrm{T}}^{\nu\nu} > $ 410 (200) GeV and normalised to a bin width of 90 (200) GeV for visualisation purposes. The total (statistical) uncertainty on 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.

png pdf
Figure 7-c:
The measured differential signal cross sections (black markers) as functions of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and double-differential in both observables (lower) are 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 cross sections in the last $ p_{\mathrm{T}}^{\nu\nu} $ bin of the upper (lower) plot are determined for $ p_{\mathrm{T}}^{\nu\nu} > $ 410 (200) GeV and normalised to a bin width of 90 (200) GeV for visualisation purposes. The total (statistical) uncertainty on 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.

png pdf
Figure 8:
The measured normalized differential signal cross sections (black markers) as functions of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and double-differential in both observables (lower) are 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 cross sections in the last $ p_{\mathrm{T}}^{\nu\nu} $ bin of the upper (lower) plot are determined for $ p_{\mathrm{T}}^{\nu\nu} > $ 410 (200) GeV and normalised to a bin width of 90 (200) GeV for visualisation purposes. The total (statistical) uncertainty on 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.

png pdf
Figure 8-a:
The measured normalized differential signal cross sections (black markers) as functions of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and double-differential in both observables (lower) are 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 cross sections in the last $ p_{\mathrm{T}}^{\nu\nu} $ bin of the upper (lower) plot are determined for $ p_{\mathrm{T}}^{\nu\nu} > $ 410 (200) GeV and normalised to a bin width of 90 (200) GeV for visualisation purposes. The total (statistical) uncertainty on 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.

png pdf
Figure 8-b:
The measured normalized differential signal cross sections (black markers) as functions of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and double-differential in both observables (lower) are 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 cross sections in the last $ p_{\mathrm{T}}^{\nu\nu} $ bin of the upper (lower) plot are determined for $ p_{\mathrm{T}}^{\nu\nu} > $ 410 (200) GeV and normalised to a bin width of 90 (200) GeV for visualisation purposes. The total (statistical) uncertainty on 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.

png pdf
Figure 8-c:
The measured normalized differential signal cross sections (black markers) as functions of $ p_{\mathrm{T}}^{\nu\nu} $ (upper), $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $ (middle), and double-differential in both observables (lower) are 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 cross sections in the last $ p_{\mathrm{T}}^{\nu\nu} $ bin of the upper (lower) plot are determined for $ p_{\mathrm{T}}^{\nu\nu} > $ 410 (200) GeV and normalised to a bin width of 90 (200) GeV for visualisation purposes. The total (statistical) uncertainty on 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.
Tables

png pdf
Table 1:
Data and MC simulation yields after the event selection, combined for all data-taking periods and split by channels. The uncertainties on the expected yields include systematic and statistical uncertainties (discussed in Section 6). The relative contribution in percent of each process to the total expected yield of a channel is given in parentheses.

png pdf
Table 2:
Definition of the fiducial phase space for the same-flavor and the different-flavor channels.

png pdf
Table 3:
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 the uncertainties on the predictions are given in parentheses.

png pdf
Table 4:
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.
Summary
Differential cross section measurements have been presented for the top quark pair production in the dileptonic channel in proton-proton (pp) collisions at a center-of-mass energy of 13 TeV using observables based on the dineutrino kinematic properties. The measurements are performed based on the pp collision data recorded by the CMS detector at the CERN LHC between 2016 and 2018, corresponding to an integrated luminosity of 138 fb$ ^{-1} $. The differential cross sections are measured as functions of the transverse momentum of the dineutrino system $ p_{\mathrm{T}}^{\nu\nu} $, and the minimal azimuthal angle between the dineutrino system and a charged lepton $ \min[\Delta\phi({\vec p}_{\mathrm{T}}^{\,\nu\nu},{\vec p}_{\mathrm{T}}^{\,\ell})] $, as well as both observables in two dimensions. To improve the resolution of the missing transverse momentum, which serves as a measure of $ p_{\mathrm{T}}^{\nu\nu} $ in signal events, a dedicated deep neural network regression has been developed. The method significantly improves the resolution of both the magnitude and the azimuthal angle of the missing transverse momentum. The absolute and normalized differential cross section results are obtained based on an unregularized least squares unfolding method. The differential cross sections are compared to predictions based on Monte Carlo simulation as well as two fixed-order theoretical calculations, corresponding to next-to-leading order and next-to-next-to-leading order (NNLO) accuracy in quantum chromodynamics. A remarkable agreement between the different theoretical predictions and the measured differential cross sections is observed. For both one-dimensional measurements, the best overall description is provided by the two POWHEG predictions, while for the two-dimensional measurement the best agreement is observed for the NNLO fixed-order calculation. However, the differences between the five predictions are mostly small, such that none of the predictions is significantly disfavored by the measured differential cross sections. These results constitute the first differential cross section measurements based on the dineutrino kinematic properties in top quark pair production.
References
1 ATLAS Collaboration Measurements of top quark pair relative differential cross-sections with ATLAS in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 7 TeV EPJC 73 (2013) 2261 1207.5644
2 CMS Collaboration Measurement of differential top-quark pair production cross sections in $ {\mathrm{p}\mathrm{p}} $ colisions at $ \sqrt{s}= $ 7 TeV EPJC 73 (2013) 2339 CMS-TOP-11-013
1211.2220
3 ATLAS Collaboration Measurements of normalized differential cross sections for $ \mathrm{t} \overline{\mathrm{t}} $ production in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 7 TeV using the ATLAS detector PRD 90 (2014) 072004 1407.0371
4 ATLAS Collaboration Differential top-antitop cross-section measurements as a function of observables constructed from final-state particles using $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 7 TeV in the ATLAS detector JHEP 06 (2015) 100 1502.05923
5 ATLAS Collaboration Measurement of top quark pair differential cross-sections in the dilepton channel in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 7 and 8 TeV with ATLAS PRD 94 (2016) 092003 1607.07281
6 CMS Collaboration Measurement of the differential cross section for top quark pair production in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s} = $ 8 TeV EPJC 75 (2015) 542 CMS-TOP-12-028
1505.04480
7 ATLAS Collaboration Measurements of top-quark pair differential cross-sections in the lepton+jets channel in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 8 TeV using the ATLAS detector EPJC 76 (2016) 538 1511.04716
8 ATLAS Collaboration Measurement of the differential cross-section of highly boosted top quarks as a function of their transverse momentum in $ \sqrt{s}= $ 8 TeV proton-proton collisions using the ATLAS detector PRD 93 (2016) 032009 1510.03818
9 CMS Collaboration Measurement of the $ \mathrm{t} \overline{\mathrm{t}} $ production cross section in the all-jets final state in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 8 TeV EPJC 76 (2016) 128 CMS-TOP-14-018
1509.06076
10 CMS Collaboration Measurement of the integrated and differential $ \mathrm{t} \overline{\mathrm{t}} $ production cross sections for high-$ p_{\mathrm{T}} $ top quarks in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 8 TeV PRD 94 (2016) 072002 CMS-TOP-14-012
1605.00116
11 CMS Collaboration Measurement of double-differential cross sections for top quark pair production in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 8 TeV and impact on parton distribution functions EPJC 77 (2017) 459 CMS-TOP-14-013
1703.01630
12 ATLAS Collaboration Measurement of the top-quark mass in $ \mathrm{t} \overline{\mathrm{t}} $ +1-jet events collected with the ATLAS detector in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 8 TeV JHEP 11 (2019) 150 1905.02302
13 CMS Collaboration Measurement of differential cross sections for top quark pair production using the lepton+jets final state in proton-proton collisions at 13 TeV PRD 95 (2017) 092001 CMS-TOP-16-008
1610.04191
14 ATLAS Collaboration Measurement of jet activity produced in top-quark events with an electron, a muon and two b-tagged jets in the final state in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector EPJC 77 (2017) 220 1610.09978
15 ATLAS Collaboration Measurements of top-quark pair differential cross-sections in the $ {\mathrm{e}\mu} $ channel in $ pp $ collisions at $ \sqrt{s}= $ 13 TeV using the ATLAS detector EPJC 77 (2017) 292 1612.05220
16 CMS Collaboration Measurement of normalized differential $ \mathrm{t} \overline{\mathrm{t}} $ cross sections in the dilepton channel from $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV JHEP 04 (2018) 060 CMS-TOP-16-007
1708.07638
17 CMS Collaboration Measurement of differential cross sections for the production of top quark pairs and of additional jets in lepton+jets events from $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV PRD 97 (2018) 112003 CMS-TOP-17-002
1803.08856
18 CMS Collaboration Measurements of $ \mathrm{t} \overline{\mathrm{t}} $ differential cross sections in proton-proton collisions at $ \sqrt{s}= $ 13 TeV using events containing two leptons JHEP 02 (2019) 149 CMS-TOP-17-014
1811.06625
19 CMS Collaboration Measurement of $ \mathrm{t} \overline{\mathrm{t}} $ normalised multi-differential cross sections in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV, and simultaneous determination of the strong coupling strength, top quark pole mass, and parton distribution functions EPJC 80 (2020) 658 CMS-TOP-18-004
1904.05237
20 ATLAS Collaboration Measurements of top-quark pair differential and double-differential cross-sections in the \ell+jets channel with $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV using the ATLAS detector EPJC 79 (2019) 1028 1908.07305
21 ATLAS Collaboration Measurement of the $ \mathrm{t} \overline{\mathrm{t}} $ production cross-section and lepton differential distributions in $ {\mathrm{e}\mu} $ dilepton events from $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector EPJC 80 (2020) 528 1910.08819
22 ATLAS Collaboration Measurements of top-quark pair single- and double-differential cross-sections in the all-hadronic channel in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV using the ATLAS detector JHEP 01 (2021) 033 2006.09274
23 CMS Collaboration Measurement of differential $ \mathrm{t} \overline{\mathrm{t}} $ production cross sections in the full kinematic range using lepton+jets events from proton-proton collisions at $ \sqrt{s}= $ 13 TeV PRD 104 (2021) 092013 CMS-TOP-20-001
2108.02803
24 S. Weinberg Supersymmetry at ordinary energies. Masses and conservation laws PRD 26 (1982) 287
25 W. Beenakker et al. NNLL resummation for stop pair-production at the LHC JHEP 05 (2016) 153 1601.02954
26 CMS Collaboration Search for top squark pair production using dilepton final states in $ {\mathrm{p}\mathrm{p}} $ collision data collected at $ \sqrt{s}= $ 13 TeV EPJC 81 (2021) 3 CMS-SUS-19-011
2008.05936
27 CMS Collaboration Combined searches for the production of supersymmetric top quark partners in proton-proton collisions at $ \sqrt{s}= $ 13 TeV EPJC 81 (2021) 970 CMS-SUS-20-002
2107.10892
28 M. Czakon, A. Mitov, and R. Poncelet NNLO QCD corrections to leptonic observables in top-quark pair production and decay JHEP 05 (2021) 212 2008.11133
29 CMS Collaboration HEPData record for this analysis link
30 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
31 CMS Collaboration Development of the CMS detector for the CERN LHC \mboxRun 3 JINST 19 (2024) P05064 CMS-PRF-21-001
2309.05466
32 CMS Collaboration Performance of the CMS Level-1 trigger in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JINST 15 (2020) P10017 CMS-TRG-17-001
2006.10165
33 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
34 CMS Collaboration Performance of the CMS high-level trigger during LHC \mboxRun 2 JINST 19 (2024) P11021 CMS-TRG-19-001
2410.17038
35 CMS Collaboration Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC JINST 16 (2021) P05014 CMS-EGM-17-001
2012.06888
36 CMS Collaboration Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at $ \sqrt{s}= $ 13 TeV JINST 13 (2018) P06015 CMS-MUO-16-001
1804.04528
37 CMS Collaboration Description and performance of track and primary-vertex reconstruction with the CMS tracker JINST 9 (2014) P10009 CMS-TRK-11-001
1405.6569
38 CMS Collaboration Precision luminosity measurement in proton-proton collisions at $ \sqrt{s}= $ 13 TeV in 2015 and 2016 at CMS EPJC 81 (2021) 800 CMS-LUM-17-003
2104.01927
39 CMS Collaboration CMS luminosity measurement for the 2017 data-taking period at $ \sqrt{s}= $ 13 TeV CMS Physics Analysis Summary, 2018
CMS-PAS-LUM-17-004
CMS-PAS-LUM-17-004
40 CMS Collaboration CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s}= $ 13 TeV CMS Physics Analysis Summary, 2019
CMS-PAS-LUM-18-002
CMS-PAS-LUM-18-002
41 NNPDF Collaboration Unbiased global determination of parton distributions and their uncertainties at NNLO and at LO NPB 855 (2012) 153 1107.2652
42 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
43 T. Sjostrand et al. An introduction to PYTHIA8.2 Comput. Phys. Commun. 191 (2015) 159 1410.3012
44 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA8 tunes from underlying-event measurements EPJC 80 (2020) 4 CMS-GEN-17-001
1903.12179
45 GEANT4 Collaboration GEANT 4---a simulation toolkit NIM A 506 (2003) 250
46 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
47 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
48 S. Alioli, P. Nason, C. Oleari, and E. Re A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG box JHEP 06 (2010) 043 1002.2581
49 J. Alwall et al. The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations JHEP 07 (2014) 079 1405.0301
50 R. Frederix and S. Frixione Merging meets matching in MC@NLO JHEP 12 (2012) 061 1209.6215
51 J. Bellm et al. HERWIG 7.0/ HERWIG++ 3.0 release note EPJC 76 (2016) 196 1512.01178
52 CMS Collaboration Development and validation of HERWIG 7 tunes from CMS underlying-event measurements EPJC 81 (2021) 312 CMS-GEN-19-001
2011.03422
53 J. Alwall et al. Comparative study of various algorithms for the merging of parton showers and matrix elements in hadronic collisions EPJC 53 (2008) 473 0706.2569
54 S. Alioli, P. Nason, C. Oleari, and E. Re NLO single-top production matched with shower in POWHEG: $ s $- and $ t $-channel contributions JHEP 09 (2009) 111 0907.4076
55 E. Re Single-top $ {\mathrm{W}\mathrm{t}} $-channel production matched with parton showers using the POWHEG method EPJC 71 (2011) 1547 1009.2450
56 M. Czakon and A. Mitov top++: a program for the calculation of the top-pair cross-section at hadron colliders Comput. Phys. Commun. 185 (2014) 2930 1112.5675
57 M. Czakon, P. Fiedler, and A. Mitov Total top-quark pair-production cross section at hadron colliders through $ \mathcal{O}({\alpha_\mathrm{S}}^4) $ PRL 110 (2013) 252004 1303.6254
58 M. Cacciari et al. Top-pair production at hadron colliders with next-to-next-to-leading logarithmic soft-gluon resummation PLB 710 (2012) 612 1111.5869
59 P. Bärnreuther, M. Czakon, and A. Mitov Percent level precision physics at the Tevatron: Next-to-next-to-leading order QCD corrections to $ \mathrm{q}\overline{\mathrm{q}}\to{\mathrm{t}\overline{\mathrm{t}}} +{\mathrm{X}} $ PRL 109 (2012) 132001 1204.5201
60 M. Czakon and A. Mitov NNLO corrections to top-pair production at hadron colliders: the all-fermionic scattering channels JHEP 12 (2012) 054 1207.0236
61 M. Czakon et al. Top-pair production at the LHC through NNLO QCD and NLO EW JHEP 10 (2017) 186 1705.04105
62 R. Gavin, Y. Li, F. Petriello, and S. Quackenbush FEWZ 2.0: A code for hadronic Z production at next-to-next-to-leading order Comput. Phys. Commun. 182 (2011) 2388 1011.3540
63 Y. Li and F. Petriello Combining QCD and electroweak corrections to dilepton production in FEWZ PRD 86 (2012) 094034 1208.5967
64 N. Kidonakis Two-loop soft anomalous dimensions for single top quark associated production with a $ \mathrm{W^-} $ or $ \mathrm{H}^{-} $ PRD 82 (2010) 054018 1005.4451
65 M. Aliev et al. hathor: Hadronic top and heavy quarks cross section calculator Comput. Phys. Commun. 182 (2011) 1034 1007.1327
66 P. Kant et al. hathor for single top-quark production: Updated predictions and uncertainty estimates for single top-quark production in hadronic collisions Comput. Phys. Commun. 191 (2015) 74 1406.4403
67 T. Gehrmann et al. $ {\mathrm{W^+}\mathrm{W^-}} $ production at hadron colliders in next to next to leading order QCD PRL 113 (2014) 212001 1408.5243
68 J. M. Campbell, R. K. Ellis, and C. Williams Vector boson pair production at the LHC JHEP 07 (2011) 018 1105.0020
69 CMS Collaboration Technical proposal for the Phase-II upgrade of the Compact Muon Solenoid CMS Technical Proposal CERN-LHCC-2015-010, CMS-TDR-15-02, 2015
link
70 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
71 CMS Collaboration Performance of the CMS electromagnetic calorimeter in $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV JINST 19 (2024) P09004 CMS-EGM-18-002
2403.15518
72 CMS Collaboration Performance of CMS muon reconstruction from proton-proton to heavy ion collisions JINST 19 (2024) P09012 CMS-MUO-21-001
2404.17377
73 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ k_{\mathrm{T}} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
74 M. Cacciari, G. P. Salam, and G. Soyez FASTJET user manual EPJC 72 (2012) 1896 1111.6097
75 CMS Collaboration Pileup mitigation at CMS in 13 TeV data JINST 15 (2020) P09018 CMS-JME-18-001
2003.00503
76 E. Bols et al. Jet flavour classification using DeepJet JINST 15 (2020) P12012 2008.10519
77 CMS Collaboration Performance of missing transverse momentum reconstruction in proton-proton collisions at $ \sqrt{s}= $ 13 TeV using the CMS detector JINST 14 (2019) P07004 CMS-JME-17-001
1903.06078
78 D. Bertolini, P. Harris, M. Low, and N. Tran Pileup per particle identification JHEP 10 (2014) 059 1407.6013
79 Particle Data Group , S. Navas et al. Review of particle physics PRD 110 (2024) 030001
80 CMS Collaboration Identification of heavy-flavour jets with the CMS detector in $ {\mathrm{p}\mathrm{p}} $ collisions at 13 TeV JINST 13 (2018) P05011 CMS-BTV-16-002
1712.07158
81 M. Abadi et al. TensorFlow: Large-scale machine learning on heterogeneous systems link
82 F. Chollet et al. keras link
83 R. Cousins Generalization of chisquare goodness-of-fit test for binned data using saturated models, with application to histograms link
84 CMS Collaboration The CMS statistical analysis and combination tool: combine Comput. Softw. Big Sci. 8 (2024) 19 CMS-CAT-23-001
2404.06614
85 CMS Collaboration Jet energy scale and resolution in the CMS experiment in $ {\mathrm{p}\mathrm{p}} $ collisions at 8 TeV JINST 12 (2017) P02014 CMS-JME-13-004
1607.03663
86 CMS Collaboration Measurement of the inelastic proton-proton cross section at $ \sqrt{s}= $ 13 TeV JHEP 07 (2018) 161 CMS-FSQ-15-005
1802.02613
87 CMS Collaboration Differential cross section measurements for the production of top quark pairs and of additional jets using dilepton events from $ {\mathrm{p}\mathrm{p}} $ collisions at $ \sqrt{s}= $ 13 TeV JHEP 02 (2025) 064 CMS-TOP-20-006
2402.08486
88 A. Andreassen and B. Nachman Neural networks for full phase-space reweighting and parameter tuning PRD 101 (2020) 091901 1907.08209
89 CMS Collaboration Reweighting simulated events using machine-learning techniques in the CMS experiment EPJC 85 (2025) 495 CMS-MLG-24-001
2411.03023
90 CMS Collaboration CMS PYTHIA8 colour reconnection tunes based on underlying-event data EPJC 83 (2023) 587 CMS-GEN-17-002
2205.02905
91 S. Schmitt TUnfold, an algorithm for correcting migration effects in high energy physics JINST 7 (2012) T10003 1205.6201
92 ATLAS Collaboration Inclusive and differential cross-sections for dilepton $ \mathrm{t} \overline{\mathrm{t}} $ production measured in $ \sqrt{s}= $ 13 TeV $ {\mathrm{p}\mathrm{p}} $ collisions with the ATLAS detector JHEP 07 (2023) 141 2303.15340
Compact Muon Solenoid
LHC, CERN