CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-TOP-22-009
Inclusive and differential cross section measurements of $ \mathrm{t}\mathrm{\bar{t}}\mathrm{b}\mathrm{\bar{b}} $ production in the lepton+jets channel at $ \sqrt{s}= $ 13 TeV with the CMS detector
Abstract: Measurements of inclusive and differential cross sections of the associated production of top quark and b quark pairs, $ \mathrm{t}\mathrm{\bar{t}}\mathrm{b}\mathrm{\bar{b}} $, are presented. The results are based on LHC data from proton-proton collisions at a centre-of-mass energy of $ \sqrt{s} = $ 13 TeV, corresponding to an integrated luminosity of 138 fb$ ^{-1} $. The cross sections are measured in the lepton+jets decay channel of the top quark pair, using events containing exactly one isolated electron or muon and at least five jets. These cross sections are measured in four fiducial phase space regions, targeting different aspects of the $ \mathrm{t}\mathrm{\bar{t}}\mathrm{b}\mathrm{\bar{b}} $ process. The jets arising from the b quark pair are probed using two complementary approaches, one relying on kinematic variables at the stable particle level and another based on a dedicated multivariate identification algorithm exploiting the simulated event history. Distributions are unfolded to the particle level through maximum likelihood fits, and compared with predictions from several event generators. The inclusive cross section measurements of the fiducial phase space regions are the most precise measurements of $ \mathrm{t}\mathrm{\bar{t}}\mathrm{b}\mathrm{\bar{b}} $ production so far. With one exception, the inclusive cross sections predicted by the considered event generators are 10-50% lower than the measured values in the different phase space regions. The differential cross sections show varying degrees of compatibility with the theoretical predictions, whereby no generator simultaneously describes all the measured distributions.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
Jet (left) and b-tagged jet (right) multiplicity in the 5j3b selection prior to any fit, shown for both lepton channels and all data periods combined. The contributions from simulation have been scaled by a common factor to match the yield in data. The $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{X} $ contribution includes the $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{H} $, $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{W} $, and $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{Z} $ processes. The shaded bands include all a-priori uncertainties described in Sec. 7, including the $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{B} $ cross section uncertainty estimated from the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ simulation. Only effects on the shape of the distributions are considered. The last bins also contain the overflow.

png pdf
Figure 1-a:
Jet (left) and b-tagged jet (right) multiplicity in the 5j3b selection prior to any fit, shown for both lepton channels and all data periods combined. The contributions from simulation have been scaled by a common factor to match the yield in data. The $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{X} $ contribution includes the $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{H} $, $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{W} $, and $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{Z} $ processes. The shaded bands include all a-priori uncertainties described in Sec. 7, including the $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{B} $ cross section uncertainty estimated from the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ simulation. Only effects on the shape of the distributions are considered. The last bins also contain the overflow.

png pdf
Figure 1-b:
Jet (left) and b-tagged jet (right) multiplicity in the 5j3b selection prior to any fit, shown for both lepton channels and all data periods combined. The contributions from simulation have been scaled by a common factor to match the yield in data. The $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{X} $ contribution includes the $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{H} $, $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{W} $, and $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{Z} $ processes. The shaded bands include all a-priori uncertainties described in Sec. 7, including the $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{B} $ cross section uncertainty estimated from the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ simulation. Only effects on the shape of the distributions are considered. The last bins also contain the overflow.

png pdf
Figure 2:
Number of jets b-tagged at the tight working point in the 5j3b (left) and 6j4b selections (right) prior to any fit, shown for all lepton channels and years combined. The contributions from simulation have been scaled by a common factor to match the yield in data. The shaded bands include all uncertainties described in Sec. 7, including the $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{B} $ cross section uncertainty estimated from the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ simulation. Only effects on the shape of the distributions are considered. The last bins also contain the overflow. The vertical dashed lines indicate the ancillary regions.

png pdf
Figure 2-a:
Number of jets b-tagged at the tight working point in the 5j3b (left) and 6j4b selections (right) prior to any fit, shown for all lepton channels and years combined. The contributions from simulation have been scaled by a common factor to match the yield in data. The shaded bands include all uncertainties described in Sec. 7, including the $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{B} $ cross section uncertainty estimated from the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ simulation. Only effects on the shape of the distributions are considered. The last bins also contain the overflow. The vertical dashed lines indicate the ancillary regions.

png pdf
Figure 2-b:
Number of jets b-tagged at the tight working point in the 5j3b (left) and 6j4b selections (right) prior to any fit, shown for all lepton channels and years combined. The contributions from simulation have been scaled by a common factor to match the yield in data. The shaded bands include all uncertainties described in Sec. 7, including the $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{B} $ cross section uncertainty estimated from the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ simulation. Only effects on the shape of the distributions are considered. The last bins also contain the overflow. The vertical dashed lines indicate the ancillary regions.

png pdf
Figure 3:
Response matrix for $ \Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\text{extra}}) $ in the 6j4b phase space. The $ x $ ($ y $) axes show the generator- (detector-)level observables. The upper figure includes the ancillary variable, unrolled on the same axis as the detector-level observable, so that the binning of the detector-level observable, stacked vertically, is repeated twice. For the lower figure, the ancillary variables are projected out to more easily visualize the correspondence between true and reconstructed values. The coloured bins show the finer binning used at reconstructed level (bins split in two), while the numbers show the values one would obtain when using the same binning at generator and detector level.

png pdf
Figure 3-a:
Response matrix for $ \Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\text{extra}}) $ in the 6j4b phase space. The $ x $ ($ y $) axes show the generator- (detector-)level observables. The upper figure includes the ancillary variable, unrolled on the same axis as the detector-level observable, so that the binning of the detector-level observable, stacked vertically, is repeated twice. For the lower figure, the ancillary variables are projected out to more easily visualize the correspondence between true and reconstructed values. The coloured bins show the finer binning used at reconstructed level (bins split in two), while the numbers show the values one would obtain when using the same binning at generator and detector level.

png pdf
Figure 3-b:
Response matrix for $ \Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\text{extra}}) $ in the 6j4b phase space. The $ x $ ($ y $) axes show the generator- (detector-)level observables. The upper figure includes the ancillary variable, unrolled on the same axis as the detector-level observable, so that the binning of the detector-level observable, stacked vertically, is repeated twice. For the lower figure, the ancillary variables are projected out to more easily visualize the correspondence between true and reconstructed values. The coloured bins show the finer binning used at reconstructed level (bins split in two), while the numbers show the values one would obtain when using the same binning at generator and detector level.

png pdf
Figure 4:
Effect of the considered sources of uncertainties on the measurement of the normalized differential cross section of the $ H_{\mathrm{T}} $ of b jets in the 5j3b phase space, obtained by combining the impacts of associated nuisance parameters according to Eq. (2). The ordering of the various sources is similar for other observables and in the other phase space regions. The last bin of the distribution is not shown, since it has no associated parameter of interest but is constrained by the other bins as described in Section 6.3. The category "other theory" includes b fragmentation, top $ p_{\mathrm{T}} $ modelling, PDF, $ \mathrm{h}_{\mathrm{damp}} $, colour reconnection, and underlying event uncertainties. The category "other experimental" includes pileup and luminosity uncertainties.

png pdf
Figure 5:
Measured inclusive cross sections for each considered phase space, compared to predictions from different $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ simulation approaches shown as coloured symbols. The blue colour is reserved for models using massive b quarks and NLO QCD $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ matrix elements, while red is used for the inclusive $ \mathrm{t} \bar{\mathrm{t}} $ generators at NLO in QCD with massless b quarks. The right panel shows the ratios between the predicted and measured cross sections.

png pdf
Figure 6:
The $ |\eta| $ of the third hardest b jet in $ p_{\mathrm{T}} $ ($ | \eta(\mathrm{b}_{3}) | $) in the 5j3b phase space (top) and the $ |\eta| $ of the sub-leading additional b jet ($ | \eta(\mathrm{b}_{2}^\text{extra}) | $) in the 6j4b phase space (bottom) after the fit to data, shown for both lepton channels and all data periods combined. The distributions are shown separately for each ancillary region, as defined in Sec. 6.1. The shaded bands include all uncertainties described in Sec. 7 after profiling the nuisance parameters in the fit, estimated by sampling the predicted yields from the fit covariance matrix. The blue line shows the sum of the predicted yields for all processes before the fit to data, using the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ samples and its corresponding cross section for the signal. The last bins contain the overflow.

png pdf
Figure 6-a:
The $ |\eta| $ of the third hardest b jet in $ p_{\mathrm{T}} $ ($ | \eta(\mathrm{b}_{3}) | $) in the 5j3b phase space (top) and the $ |\eta| $ of the sub-leading additional b jet ($ | \eta(\mathrm{b}_{2}^\text{extra}) | $) in the 6j4b phase space (bottom) after the fit to data, shown for both lepton channels and all data periods combined. The distributions are shown separately for each ancillary region, as defined in Sec. 6.1. The shaded bands include all uncertainties described in Sec. 7 after profiling the nuisance parameters in the fit, estimated by sampling the predicted yields from the fit covariance matrix. The blue line shows the sum of the predicted yields for all processes before the fit to data, using the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ samples and its corresponding cross section for the signal. The last bins contain the overflow.

png pdf
Figure 6-b:
The $ |\eta| $ of the third hardest b jet in $ p_{\mathrm{T}} $ ($ | \eta(\mathrm{b}_{3}) | $) in the 5j3b phase space (top) and the $ |\eta| $ of the sub-leading additional b jet ($ | \eta(\mathrm{b}_{2}^\text{extra}) | $) in the 6j4b phase space (bottom) after the fit to data, shown for both lepton channels and all data periods combined. The distributions are shown separately for each ancillary region, as defined in Sec. 6.1. The shaded bands include all uncertainties described in Sec. 7 after profiling the nuisance parameters in the fit, estimated by sampling the predicted yields from the fit covariance matrix. The blue line shows the sum of the predicted yields for all processes before the fit to data, using the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ samples and its corresponding cross section for the signal. The last bins contain the overflow.

png pdf
Figure 7:
The $ H_{\mathrm{T}} $ of all light jets in the 6j3b3l phase space (top) and the azimuthal angle between the hardest remaining light jet and the softest b jet ($ |\Delta\phi({\mathrm{l_{j1}}^{\text{extra}},\mathrm{b}_{\text{soft}}})| $) in the 7j4b3l phase space (bottom) after the fit to data, shown for both lepton channels and all data periods combined. The distributions are shown separately for each ancillary region, as defined in Sec. 6.1. The shaded bands include all uncertainties described in Sec. 7 after profiling the nuisance parameters in the fit, estimated by sampling the predicted yields from the fit covariance matrix. The blue line shows the sum of the predicted yields for all processes before the fit to data, using the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ samples and its corresponding cross section for the signal. The last bins contain the overflow.

png pdf
Figure 7-a:
The $ H_{\mathrm{T}} $ of all light jets in the 6j3b3l phase space (top) and the azimuthal angle between the hardest remaining light jet and the softest b jet ($ |\Delta\phi({\mathrm{l_{j1}}^{\text{extra}},\mathrm{b}_{\text{soft}}})| $) in the 7j4b3l phase space (bottom) after the fit to data, shown for both lepton channels and all data periods combined. The distributions are shown separately for each ancillary region, as defined in Sec. 6.1. The shaded bands include all uncertainties described in Sec. 7 after profiling the nuisance parameters in the fit, estimated by sampling the predicted yields from the fit covariance matrix. The blue line shows the sum of the predicted yields for all processes before the fit to data, using the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ samples and its corresponding cross section for the signal. The last bins contain the overflow.

png pdf
Figure 7-b:
The $ H_{\mathrm{T}} $ of all light jets in the 6j3b3l phase space (top) and the azimuthal angle between the hardest remaining light jet and the softest b jet ($ |\Delta\phi({\mathrm{l_{j1}}^{\text{extra}},\mathrm{b}_{\text{soft}}})| $) in the 7j4b3l phase space (bottom) after the fit to data, shown for both lepton channels and all data periods combined. The distributions are shown separately for each ancillary region, as defined in Sec. 6.1. The shaded bands include all uncertainties described in Sec. 7 after profiling the nuisance parameters in the fit, estimated by sampling the predicted yields from the fit covariance matrix. The blue line shows the sum of the predicted yields for all processes before the fit to data, using the nominal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ samples and its corresponding cross section for the signal. The last bins contain the overflow.

png pdf
Figure 8:
Correlations between the parameters of interest $ \vec{\mu} $ in the fit for $ | \eta(\mathrm{b}_{3}) | $ in the 5j3b phase space.

png pdf
Figure 9:
Predicted and observed normalized differential cross sections in the 5j3b fiducial phase space, for the inclusive jet multiplicity (upper left), the b-jet multiplicity (upper right), the inclusive jet $ H_{\mathrm{T}} $ (middle left, $ H^{j}_{\mathrm{T}} $), the $ H_{\mathrm{T}} $ of b jets (middle right, $ H^{\mathrm{b}}_{\mathrm{T}} $), the $ |\eta| $ of the third b jet (lower left), and the $ p_{\mathrm{T}} $ of the third b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ N_{\mathrm{b}} $, $ N_{\text{jets}} $, $ H_{\mathrm{T}} $, and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 9-a:
Predicted and observed normalized differential cross sections in the 5j3b fiducial phase space, for the inclusive jet multiplicity (upper left), the b-jet multiplicity (upper right), the inclusive jet $ H_{\mathrm{T}} $ (middle left, $ H^{j}_{\mathrm{T}} $), the $ H_{\mathrm{T}} $ of b jets (middle right, $ H^{\mathrm{b}}_{\mathrm{T}} $), the $ |\eta| $ of the third b jet (lower left), and the $ p_{\mathrm{T}} $ of the third b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ N_{\mathrm{b}} $, $ N_{\text{jets}} $, $ H_{\mathrm{T}} $, and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 9-b:
Predicted and observed normalized differential cross sections in the 5j3b fiducial phase space, for the inclusive jet multiplicity (upper left), the b-jet multiplicity (upper right), the inclusive jet $ H_{\mathrm{T}} $ (middle left, $ H^{j}_{\mathrm{T}} $), the $ H_{\mathrm{T}} $ of b jets (middle right, $ H^{\mathrm{b}}_{\mathrm{T}} $), the $ |\eta| $ of the third b jet (lower left), and the $ p_{\mathrm{T}} $ of the third b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ N_{\mathrm{b}} $, $ N_{\text{jets}} $, $ H_{\mathrm{T}} $, and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 9-c:
Predicted and observed normalized differential cross sections in the 5j3b fiducial phase space, for the inclusive jet multiplicity (upper left), the b-jet multiplicity (upper right), the inclusive jet $ H_{\mathrm{T}} $ (middle left, $ H^{j}_{\mathrm{T}} $), the $ H_{\mathrm{T}} $ of b jets (middle right, $ H^{\mathrm{b}}_{\mathrm{T}} $), the $ |\eta| $ of the third b jet (lower left), and the $ p_{\mathrm{T}} $ of the third b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ N_{\mathrm{b}} $, $ N_{\text{jets}} $, $ H_{\mathrm{T}} $, and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 9-d:
Predicted and observed normalized differential cross sections in the 5j3b fiducial phase space, for the inclusive jet multiplicity (upper left), the b-jet multiplicity (upper right), the inclusive jet $ H_{\mathrm{T}} $ (middle left, $ H^{j}_{\mathrm{T}} $), the $ H_{\mathrm{T}} $ of b jets (middle right, $ H^{\mathrm{b}}_{\mathrm{T}} $), the $ |\eta| $ of the third b jet (lower left), and the $ p_{\mathrm{T}} $ of the third b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ N_{\mathrm{b}} $, $ N_{\text{jets}} $, $ H_{\mathrm{T}} $, and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 9-e:
Predicted and observed normalized differential cross sections in the 5j3b fiducial phase space, for the inclusive jet multiplicity (upper left), the b-jet multiplicity (upper right), the inclusive jet $ H_{\mathrm{T}} $ (middle left, $ H^{j}_{\mathrm{T}} $), the $ H_{\mathrm{T}} $ of b jets (middle right, $ H^{\mathrm{b}}_{\mathrm{T}} $), the $ |\eta| $ of the third b jet (lower left), and the $ p_{\mathrm{T}} $ of the third b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ N_{\mathrm{b}} $, $ N_{\text{jets}} $, $ H_{\mathrm{T}} $, and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 9-f:
Predicted and observed normalized differential cross sections in the 5j3b fiducial phase space, for the inclusive jet multiplicity (upper left), the b-jet multiplicity (upper right), the inclusive jet $ H_{\mathrm{T}} $ (middle left, $ H^{j}_{\mathrm{T}} $), the $ H_{\mathrm{T}} $ of b jets (middle right, $ H^{\mathrm{b}}_{\mathrm{T}} $), the $ |\eta| $ of the third b jet (lower left), and the $ p_{\mathrm{T}} $ of the third b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ N_{\mathrm{b}} $, $ N_{\text{jets}} $, $ H_{\mathrm{T}} $, and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 10:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the inclusive jet $ H_{\mathrm{T}} $ (upper left), the $ H_{\mathrm{T}} $ of b jets (upper right), the $ |\eta| $ of the third b jet (middle left), the $ p_{\mathrm{T}} $ of the third b jet (middle right), the $ |\eta| $ of the fourth b jet (lower left), and the $ p_{\mathrm{T}} $ of the fourth b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 10-a:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the inclusive jet $ H_{\mathrm{T}} $ (upper left), the $ H_{\mathrm{T}} $ of b jets (upper right), the $ |\eta| $ of the third b jet (middle left), the $ p_{\mathrm{T}} $ of the third b jet (middle right), the $ |\eta| $ of the fourth b jet (lower left), and the $ p_{\mathrm{T}} $ of the fourth b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 10-b:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the inclusive jet $ H_{\mathrm{T}} $ (upper left), the $ H_{\mathrm{T}} $ of b jets (upper right), the $ |\eta| $ of the third b jet (middle left), the $ p_{\mathrm{T}} $ of the third b jet (middle right), the $ |\eta| $ of the fourth b jet (lower left), and the $ p_{\mathrm{T}} $ of the fourth b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 10-c:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the inclusive jet $ H_{\mathrm{T}} $ (upper left), the $ H_{\mathrm{T}} $ of b jets (upper right), the $ |\eta| $ of the third b jet (middle left), the $ p_{\mathrm{T}} $ of the third b jet (middle right), the $ |\eta| $ of the fourth b jet (lower left), and the $ p_{\mathrm{T}} $ of the fourth b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 10-d:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the inclusive jet $ H_{\mathrm{T}} $ (upper left), the $ H_{\mathrm{T}} $ of b jets (upper right), the $ |\eta| $ of the third b jet (middle left), the $ p_{\mathrm{T}} $ of the third b jet (middle right), the $ |\eta| $ of the fourth b jet (lower left), and the $ p_{\mathrm{T}} $ of the fourth b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 10-e:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the inclusive jet $ H_{\mathrm{T}} $ (upper left), the $ H_{\mathrm{T}} $ of b jets (upper right), the $ |\eta| $ of the third b jet (middle left), the $ p_{\mathrm{T}} $ of the third b jet (middle right), the $ |\eta| $ of the fourth b jet (lower left), and the $ p_{\mathrm{T}} $ of the fourth b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 10-f:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the inclusive jet $ H_{\mathrm{T}} $ (upper left), the $ H_{\mathrm{T}} $ of b jets (upper right), the $ |\eta| $ of the third b jet (middle left), the $ p_{\mathrm{T}} $ of the third b jet (middle right), the $ |\eta| $ of the fourth b jet (lower left), and the $ p_{\mathrm{T}} $ of the fourth b jet (lower right). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 11:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the average $ \Delta\mathrm{R} $ of all possible $ \mathrm{b} \bar{\mathrm{b}} $\ pairs (upper left), the largest invariant mass of any $ \mathrm{b} \bar{\mathrm{b}} $\ pair (upper right), the invariant mass (middle left), $ \Delta\mathrm{R} $ (middle right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the extra b-jet pair. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 11-a:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the average $ \Delta\mathrm{R} $ of all possible $ \mathrm{b} \bar{\mathrm{b}} $\ pairs (upper left), the largest invariant mass of any $ \mathrm{b} \bar{\mathrm{b}} $\ pair (upper right), the invariant mass (middle left), $ \Delta\mathrm{R} $ (middle right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the extra b-jet pair. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 11-b:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the average $ \Delta\mathrm{R} $ of all possible $ \mathrm{b} \bar{\mathrm{b}} $\ pairs (upper left), the largest invariant mass of any $ \mathrm{b} \bar{\mathrm{b}} $\ pair (upper right), the invariant mass (middle left), $ \Delta\mathrm{R} $ (middle right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the extra b-jet pair. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 11-c:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the average $ \Delta\mathrm{R} $ of all possible $ \mathrm{b} \bar{\mathrm{b}} $\ pairs (upper left), the largest invariant mass of any $ \mathrm{b} \bar{\mathrm{b}} $\ pair (upper right), the invariant mass (middle left), $ \Delta\mathrm{R} $ (middle right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the extra b-jet pair. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 11-d:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the average $ \Delta\mathrm{R} $ of all possible $ \mathrm{b} \bar{\mathrm{b}} $\ pairs (upper left), the largest invariant mass of any $ \mathrm{b} \bar{\mathrm{b}} $\ pair (upper right), the invariant mass (middle left), $ \Delta\mathrm{R} $ (middle right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the extra b-jet pair. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 11-e:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the average $ \Delta\mathrm{R} $ of all possible $ \mathrm{b} \bar{\mathrm{b}} $\ pairs (upper left), the largest invariant mass of any $ \mathrm{b} \bar{\mathrm{b}} $\ pair (upper right), the invariant mass (middle left), $ \Delta\mathrm{R} $ (middle right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the extra b-jet pair. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 11-f:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the average $ \Delta\mathrm{R} $ of all possible $ \mathrm{b} \bar{\mathrm{b}} $\ pairs (upper left), the largest invariant mass of any $ \mathrm{b} \bar{\mathrm{b}} $\ pair (upper right), the invariant mass (middle left), $ \Delta\mathrm{R} $ (middle right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the extra b-jet pair. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 12:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the $ |\eta| $ (upper left) and $ p_{\mathrm{T}} $ (upper right) of the first extra b jet, the $ |\eta| $ (middle left) and $ p_{\mathrm{T}} $ (middle right) of the second extra b jet, and the inclusive jet multiplicity (lower left). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ N_{\text{jets}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 12-a:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the $ |\eta| $ (upper left) and $ p_{\mathrm{T}} $ (upper right) of the first extra b jet, the $ |\eta| $ (middle left) and $ p_{\mathrm{T}} $ (middle right) of the second extra b jet, and the inclusive jet multiplicity (lower left). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ N_{\text{jets}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 12-b:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the $ |\eta| $ (upper left) and $ p_{\mathrm{T}} $ (upper right) of the first extra b jet, the $ |\eta| $ (middle left) and $ p_{\mathrm{T}} $ (middle right) of the second extra b jet, and the inclusive jet multiplicity (lower left). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ N_{\text{jets}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 12-c:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the $ |\eta| $ (upper left) and $ p_{\mathrm{T}} $ (upper right) of the first extra b jet, the $ |\eta| $ (middle left) and $ p_{\mathrm{T}} $ (middle right) of the second extra b jet, and the inclusive jet multiplicity (lower left). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ N_{\text{jets}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 12-d:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the $ |\eta| $ (upper left) and $ p_{\mathrm{T}} $ (upper right) of the first extra b jet, the $ |\eta| $ (middle left) and $ p_{\mathrm{T}} $ (middle right) of the second extra b jet, and the inclusive jet multiplicity (lower left). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ N_{\text{jets}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 12-e:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the $ |\eta| $ (upper left) and $ p_{\mathrm{T}} $ (upper right) of the first extra b jet, the $ |\eta| $ (middle left) and $ p_{\mathrm{T}} $ (middle right) of the second extra b jet, and the inclusive jet multiplicity (lower left). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ N_{\text{jets}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 13:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the invariant mass (upper left), $ \Delta\mathrm{R} $ (upper right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the additional b-jet pair not originating from decaying top quarks. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 13-a:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the invariant mass (upper left), $ \Delta\mathrm{R} $ (upper right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the additional b-jet pair not originating from decaying top quarks. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 13-b:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the invariant mass (upper left), $ \Delta\mathrm{R} $ (upper right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the additional b-jet pair not originating from decaying top quarks. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 13-c:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the invariant mass (upper left), $ \Delta\mathrm{R} $ (upper right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the additional b-jet pair not originating from decaying top quarks. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 13-d:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the invariant mass (upper left), $ \Delta\mathrm{R} $ (upper right), $ p_{\mathrm{T}} $ (lower left), and $ |\eta| $ (lower right) of the additional b-jet pair not originating from decaying top quarks. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ \mathrm{m}_{\mathrm{b}\mathrm{b}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 14:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the $ |\eta| $ (left) and $ p_{\mathrm{T}} $ (right) of the first (upper row) and second (lower row) additional b of the b-jet pair not originating from decaying top quarks. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 14-a:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the $ |\eta| $ (left) and $ p_{\mathrm{T}} $ (right) of the first (upper row) and second (lower row) additional b of the b-jet pair not originating from decaying top quarks. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 14-b:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the $ |\eta| $ (left) and $ p_{\mathrm{T}} $ (right) of the first (upper row) and second (lower row) additional b of the b-jet pair not originating from decaying top quarks. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 14-c:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the $ |\eta| $ (left) and $ p_{\mathrm{T}} $ (right) of the first (upper row) and second (lower row) additional b of the b-jet pair not originating from decaying top quarks. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 14-d:
Predicted and observed normalized differential cross sections in the 6j4b fiducial phase space, for the $ |\eta| $ (left) and $ p_{\mathrm{T}} $ (right) of the first (upper row) and second (lower row) additional b of the b-jet pair not originating from decaying top quarks. The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols. For $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 15:
Predicted and observed normalized differential cross sections in the 6j3b3l (left) and 7j4b3l (right) fiducial phase space regions, for the $ H_{\mathrm{T}} $ of light jets (upper row), the $ p_{\mathrm{T}} $ of the extra light jet (middle row), and the $ \Delta\phi $ between the extra light jet and the softest b jet (lower row). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 15-a:
Predicted and observed normalized differential cross sections in the 6j3b3l (left) and 7j4b3l (right) fiducial phase space regions, for the $ H_{\mathrm{T}} $ of light jets (upper row), the $ p_{\mathrm{T}} $ of the extra light jet (middle row), and the $ \Delta\phi $ between the extra light jet and the softest b jet (lower row). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 15-b:
Predicted and observed normalized differential cross sections in the 6j3b3l (left) and 7j4b3l (right) fiducial phase space regions, for the $ H_{\mathrm{T}} $ of light jets (upper row), the $ p_{\mathrm{T}} $ of the extra light jet (middle row), and the $ \Delta\phi $ between the extra light jet and the softest b jet (lower row). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 15-c:
Predicted and observed normalized differential cross sections in the 6j3b3l (left) and 7j4b3l (right) fiducial phase space regions, for the $ H_{\mathrm{T}} $ of light jets (upper row), the $ p_{\mathrm{T}} $ of the extra light jet (middle row), and the $ \Delta\phi $ between the extra light jet and the softest b jet (lower row). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 15-d:
Predicted and observed normalized differential cross sections in the 6j3b3l (left) and 7j4b3l (right) fiducial phase space regions, for the $ H_{\mathrm{T}} $ of light jets (upper row), the $ p_{\mathrm{T}} $ of the extra light jet (middle row), and the $ \Delta\phi $ between the extra light jet and the softest b jet (lower row). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 15-e:
Predicted and observed normalized differential cross sections in the 6j3b3l (left) and 7j4b3l (right) fiducial phase space regions, for the $ H_{\mathrm{T}} $ of light jets (upper row), the $ p_{\mathrm{T}} $ of the extra light jet (middle row), and the $ \Delta\phi $ between the extra light jet and the softest b jet (lower row). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 15-f:
Predicted and observed normalized differential cross sections in the 6j3b3l (left) and 7j4b3l (right) fiducial phase space regions, for the $ H_{\mathrm{T}} $ of light jets (upper row), the $ p_{\mathrm{T}} $ of the extra light jet (middle row), and the $ \Delta\phi $ between the extra light jet and the softest b jet (lower row). The data are represented by points, with inner (outer) vertical bars indicating the systematic (total) uncertainties, also represented as grey (yellow) bands. Cross section predictions obtained at particle level from different simulation approaches are shown as coloured symbols with different shapes. For $ H_{\mathrm{T}} $ and $ p_{\mathrm{T}} $, the last bins contain the overflow.

png pdf
Figure 16:
Observed $ z $ score for each of the theoretical predictions, given the unfolded normalized differential cross sections and their covariances. A lower value indicates a better agreement between prediction and measurement. The dashed line at $ z= $ 2 indicates a $ p $-value of 5%. Predictions for which the $ z $ score exceeds the visible range of the figure are marked with arrows ($ \rightarrow $).

png pdf
Figure A1:
$ |\eta| $ of the b jet with third highest $ p_{\mathrm{T}} $ in the 5j3b phase space.

png pdf
Figure A2:
$ |\eta| $ of the subleading extra b jet in the 6j4b phase space.

png pdf
Figure A3:
$ H_{\mathrm{T}} $ of light jets in the 6j3b3l phase space.

png pdf
Figure A4:
$ \Delta\phi $ between leading light jet and softest b jet in the 7j4b3l phase space.

png pdf
Figure B1:
$ |\eta| $ of the b jet with third highest $ p_{\mathrm{T}} $ in the 5j3b phase space.

png pdf
Figure B2:
$ |\eta| $ of the subleading extra b jet in the 6j4b phase space.

png pdf
Figure B3:
$ H_{\mathrm{T}} $ of light jets in the 6j3b3l phase space.

png pdf
Figure B4:
$ \Delta\phi $ between leading light jet and softest b jet in the 7j4b3l phase space.

png pdf
Figure C1:
$ |\eta| $ of the subleading extra b jet in the 6j4b phase space.

png pdf
Figure C2:
$ H_{\mathrm{T}} $ of light jets in the 6j3b3l phase space.

png pdf
Figure C3:
$ \Delta\phi $ between leading light jet and softest b jet in the 7j4b3l phase space.
Tables

png pdf
Table 1:
Description of all measured observables for each of the four fiducial phase space regions. Observables marked as checkmark rely on the definition of additional b jets, and do not fully correspond to the 6j4b fiducial phase space defined at particle level, but also require the presence of b jets without top (anti)quarks in their simulated history.

png pdf
Table 2:
Summary of the systematic uncertainty sources in the inclusive and differential $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ cross section measurements. The first column lists the source of the uncertainty. The second (third) column indicates the treatment of correlations between the uncertainties between the data-taking periods (processes), where a checkmark means fully correlated, $ \sim $ means partially correlated (i.e.,, contains sub-sources that are either fully correlated or uncorrelated), $ \times $ means uncorrelated, and --- means not applicable.

png pdf
Table 3:
Contributions of the considered sources of uncertainty to the total uncertainty in the inclusive cross sections. For each source, the impacts of the corresponding nuisance parameters on the total cross section are combined, taking into account their correlation in the fit. The numbers show relative uncertainties (in %). The statistical uncertainty is obtained as the difference, in quadrature, between the total uncertainty and the sum of all systematic uncertainties.

png pdf
Table 4:
Measured and predicted inclusive cross sections in the four considered phase space regions (in fb).

png pdf
Table D1:
Observed $ z $ score for each of the theoretical predictions in the 5j3b phase space, given the unfolded data and covariance matrix. For the determination of the $ z $ score only the measurement uncertainties are considered.

png pdf
Table D2:
Observed z-score for each of the theoretical predictions in the 6j4b phase space, given the unfolded data and covariance matrix. For the determination of the $ z $ score only the measurement uncertainties are considered.

png pdf
Table D3:
Observed $ z $ score for each of the theoretical predictions in the 6j4b phase space of the observables related to the $ \mathrm{b}\mathrm{b}^{\text{extra}} $ pair, given the unfolded data and covariance matrix. For the determination of the $ z $ score only the measurement uncertainties are considered.

png pdf
Table D4:
Observed $ z $ score for each of the theoretical predictions in the 6j4b phase space of the observables related to the $ \mathrm{b}\mathrm{b}^{\text{add.}} $ pair, given the unfolded data and covariance matrix. For the determination of the $ z $ score only the measurement uncertainties are considered.

png pdf
Table D5:
Observed $ z $ score for each of the theoretical predictions in the 6j3b3l phase space, given the unfolded data and covariance matrix. For the determination of the $ z $ score only the measurement uncertainties are considered.

png pdf
Table D6:
Observed $ z $ score for each of the theoretical predictions in the 7j4b3l phase space, given the unfolded data and covariance matrix. For the determination of the $ z $ score only the measurement uncertainties are considered.
Summary
Measurements of inclusive and normalized differential cross sections of $ \mathrm{t} \bar{\mathrm{t}} $ production in association with b jets, for events containing an electron or a muon, have been presented. These measurements use pp collision data recorded by the CMS detector at $ \sqrt{s} = $ 13 TeV and correspond to an integrated luminosity of 138 fb$^{-1}$. The inclusive cross sections are measured in four phase space regions requiring different jet, b-jet, and light-jet multiplicities. With total uncertainties of 6--18%, depending on the phase space, these are the most precise measurements of the $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ cross section to date. The uncertainties are dominated by systematic sources, with the leading uncertainties originating from the calibration of the b tagging and of the jet energy scale, and from the choice of renormalization scale in the signal $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ and background $ \mathrm{t} \bar{\mathrm{t}} $ processes. Differential cross section measurements are performed as a function of several observables in the aforementioned phase space regions. These observables mainly target b jets as well as additional light jets produced in association with the top quark pairs. In the 6j4b phase space, the additional b-jet radiation is probed with two different approaches. The first approach uses observables defined purely at the particle level, without any reference to the top quark decay chains, by selecting the two b jets closest in $ \Delta\mathrm{R} $. The second approach uses explicitly at generator level the b jets which do not originate from top quark decays, and identifies those jets at the detector level with a deep neural network discriminant. The differential measurements have relative uncertainties in the range of 2--50%, depending on the phase space and the observable. The results are compared to the predictions of several event generators, and it is found that none of them simultaneously describe all measured distributions in the various phase space regions. In the more inclusive phase space with five jets and three b jets, the agreement between data and predictions is generally poor, while in the phase space with six jets and four b jets, corresponding to the case in which the two additional b jets in $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ production are resolved, most predictions are compatible with the data within the larger experimental uncertainties. These measurements of $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ production will help to further tune and refine the theoretical predictions and better assess the validity of the theoretical uncertainties estimated from the various $ {\mathrm{t}\bar{\mathrm{t}}} \mathrm{b}\bar{\mathrm{b}} $ event generators.
References
1 F. Buccioni, S. Kallweit, S. Pozzorini, and M. F. Zoller NLO QCD predictions for $ \mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}} $ production in association with a light jet at the LHC JHEP 12 (2019) 015 1907.13624
2 T. Ježo, J. M. Lindert, N. Moretti, and S. Pozzorini New NLOPS predictions for $ \mathrm{t}\bar{\mathrm{t}} $ + b-jet production at the LHC EPJC 78 (2018) 502 1802.00426
3 ATLAS Collaboration Search for the standard model Higgs boson decaying into $ \mathrm{b}\bar{\mathrm{b}} $ produced in association with top quarks decaying hadronically in pp collisions at $ \sqrt{s}= $ 8 TeV with the ATLAS detector JHEP 05 (2016) 160 1604.03812
4 ATLAS Collaboration Search for the standard model Higgs boson produced in association with top quarks and decaying into a $ \mathrm{b}\bar{\mathrm{b}} $ pair in pp collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector PRD 97 (2018) 072016 1712.08895
5 ATLAS Collaboration Observation of Higgs boson production in association with a top quark pair at the LHC with the ATLAS detector PLB 784 (2018) 173 1806.00425
6 CMS Collaboration Search for $ \mathrm{t}\bar{\mathrm{t}}\mathrm{H} $ production in the all-jet final state in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 06 (2018) 101 CMS-HIG-17-022
1803.06986
7 CMS Collaboration Observation of $ \mathrm{t}\bar{\mathrm{t}}\mathrm{H} $ production PRL 120 (2018) 231801 CMS-HIG-17-035
1804.02610
8 CMS Collaboration Search for $ \mathrm{t}\bar{\mathrm{t}}\mathrm{H} $ production in the $ \mathrm{H} \to \mathrm{b}\bar{\mathrm{b}} $ decay channel with leptonic $ \mathrm{t}\bar{\mathrm{t}} $ decays in proton-proton collisions at $ \sqrt{s}= $ 13 TeV JHEP 03 (2019) 026 CMS-HIG-17-026
1804.03682
9 ATLAS Collaboration Search for four-top-quark production in the single-lepton and opposite-sign dilepton final states in pp collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector PRD 99 (2019) 052009 1811.02305
10 ATLAS Collaboration Evidence for $ \mathrm{t}\bar{\mathrm{t}}\mathrm{t}\bar{\mathrm{t}} $ production in the multilepton final state in proton-proton collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector EPJC 80 (2020) 1085 2007.14858
11 ATLAS Collaboration Measurement of the $ \mathrm{t}\bar{\mathrm{t}}\mathrm{t}\bar{\mathrm{t}} $ production cross section in $ pp $ collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector JHEP 11 (2021) 118 2106.11683
12 CMS Collaboration Search for standard model production of four top quarks in the lepton + jets channel in pp collisions at $ \sqrt{s}= $ 8 TeV JHEP 11 (2014) 154 CMS-TOP-13-012
1409.7339
13 CMS Collaboration Search for standard model production of four top quarks with same-sign and multilepton final states in proton-proton collisions at $ \sqrt{s}= $ 13 TeV EPJC 78 (2018) 140 CMS-TOP-17-009
1710.10614
14 CMS Collaboration Search for standard model production of four top quarks in proton-proton collisions at $ \sqrt{s}= $ 13 TeV PLB 772 (2017) 336 CMS-TOP-16-016
1702.06164
15 CMS Collaboration Search for the production of four top quarks in the single-lepton and opposite-sign dilepton final states in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JHEP 11 (2019) 082 CMS-TOP-17-019
1906.02805
16 CMS Collaboration Search for production of four top quarks in final states with same-sign or multiple leptons in proton-proton collisions at $ \sqrt{s}= $ 13 TeV EPJC 80 (2020) 75 CMS-TOP-18-003
1908.06463
17 CMS Collaboration Evidence for four-top quark production in proton-proton collisions at $ \sqrt{s} $ = 13 TeV Submitted to Phys. Lett. B, 2023 CMS-TOP-21-005
2303.03864
18 Q.-H. Cao, S.-L. Chen, and Y. Liu Probing Higgs width and top quark Yukawa coupling from $ \mathrm{t}\bar{\mathrm{t}}\mathrm{H} $ and $ \mathrm{t}\bar{\mathrm{t}}\mathrm{t}\bar{\mathrm{t}} $ productions PRD 95 (2017) 053004 1602.01934
19 Q.-H. Cao et al. Limiting top quark-Higgs boson interaction and Higgs-boson width from multitop productions PRD 99 (2019) 113003 1901.04567
20 M. Worek Next-to-leading order QCD corrections to $ \mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}} $ production at the LHC Acta Phys. Polon. B 40 (2009) 2937 0910.4080
21 G. Bevilacqua et al. Assault on the NLO wishlist: $ \mathrm{p}\mathrm{p} \rightarrow \mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}} $ JHEP 09 (2009) 109 0907.4723
22 A. Bredenstein, A. Denner, S. Dittmaier, and S. Pozzorini NLO QCD corrections to $ \mathrm{p}\mathrm{p} \rightarrow \mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}}$+X at the LHC PRL 103 (2009) 012002 0905.0110
23 A. Bredenstein, A. Denner, S. Dittmaier, and S. Pozzorini NLO QCD corrections to $ \mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}} $ production at the LHC: 2. Full hadronic results JHEP 03 (2010) 021 1001.4006
24 M. Worek On the next-to-leading order QCD $ \mathcal{K} $-factor for $ \mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}} $ production at the TeVatron JHEP 02 (2012) 043 1112.4325
25 A. Denner, J.-N. Lang, and M. Pellen Full NLO QCD corrections to off-shell $ \mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}} $ production PRD 104 (2021) 056018 2008.00918
26 G. Bevilacqua et al. $ \mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}} $ at the LHC: on the size of corrections and b-jet definitions JHEP 08 (2021) 008 2105.08404
27 S. Höche et al. Next-to-leading order QCD predictions for top-quark pair production with up to two jets merged with a parton shower PLB 748 (2015) 74 1402.6293
28 M. V. Garzelli, A. Kardos, and Z. Trócsányi Hadroproduction of $ \mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}} $ final states at LHC: predictions at NLO accuracy matched with parton shower JHEP 03 (2015) 083 1408.0266
29 S. Höche et al. Next-to-leading order QCD predictions for top-quark pair production with up to three jets EPJC 77 (2017) 145 1607.06934
30 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
31 R. Frederix and S. Frixione Merging meets matching in MC@NLO JHEP 12 (2012) 061 1209.6215
32 S. Frixione, P. Nason, and G. Ridolfi A positive-weight next-to-leading-order Monte Carlo for heavy flavour hadroproduction JHEP 09 (2007) 126 0707.3088
33 F. Cascioli et al. NLO matching for $ \mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}} $ production with massive b-quarks PLB 734 (2014) 210 1309.5912
34 G. Bevilacqua, M. V. Garzelli, and A. Kardos $ \mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}} $ hadroproduction with massive bottom quarks with PowHel 1709.06915
35 ATLAS Collaboration Study of heavy-flavor quarks produced in association with top-quark pairs at $ \sqrt{s}= $ 7 TeV using the ATLAS detector PRD 89 (2014) 072012 1304.6386
36 ATLAS Collaboration Measurements of fiducial cross-sections for $ \mathrm{t}\bar{\mathrm{t}} $ production with one or two additional b-jets in pp collisions at $ \sqrt{s}= $ 8 TeV using the ATLAS detector EPJC 76 (2016) 11 1508.06868
37 ATLAS Collaboration Measurements of inclusive and differential fiducial cross-sections of $ \mathrm{t}\bar{\mathrm{t}} $ production with additional heavy-flavour jets in proton-proton collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector JHEP 04 (2019) 046 1811.12113
38 CMS Collaboration Measurement of $ \mathrm{t}\bar{\mathrm{t}} $ production with additional jet activity, including b quark jets, in the dilepton decay channel using pp collisions at $ \sqrt{s}= $ 8 TeV EPJC 76 (2016) 379 CMS-TOP-12-041
1510.03072
39 CMS Collaboration Measurement of the cross section ratio $ \sigma_{\mathrm{t}\bar{\mathrm{t}}\mathrm{b}\bar{\mathrm{b}}} / \sigma_\mathrm{\mathrm{t}\bar{\mathrm{t}} \mathrm{jj}} $ in pp collisions at $ \sqrt{s}= $ 8 TeV PLB 746 (2015) 132 CMS-TOP-13-010
1411.5621
40 CMS Collaboration Measurements of $ \mathrm{t}\bar{\mathrm{t}} $ cross sections in association with b jets and inclusive jets and their ratio using dilepton final states in pp collisions at $ \sqrt{s}= $ 13 TeV PLB 776 (2018) 355 CMS-TOP-16-010
1705.10141
41 CMS Collaboration Measurement of the $ \mathrm{t\bar{t}}\mathrm{b\bar{b}} $ production cross section in the all-jet final state in pp collisions at $ \sqrt{s}= $ 13 TeV PLB 803 (2020) 135285 CMS-TOP-18-011
1909.05306
42 CMS Collaboration Measurement of the cross section for $ \mathrm{t\bar{t}}\mathrm{b\bar{b}} $ production with additional jets and b jets in pp collisions at $ \sqrt{s}= $ 13 TeV JHEP 07 (2020) 125 CMS-TOP-18-002
2003.06467
43 CMS Collaboration First measurement of the cross section for top quark pair production with additional charm jets using dileptonic final states in pp collisions at $\sqrt{s}=$ 13 TeV PLB 820 (2021) 136565 CMS-TOP-20-003
2012.09225
44 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
45 CMS Tracker Group Collaboration The CMS Phase-1 Pixel Detector Upgrade JINST 16 (2021) P02027 2012.14304
46 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
47 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
48 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
49 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
50 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
51 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
52 E. Re Single-top $ \mathrm{W}\mathrm{t} $-channel production matched with parton showers using the POWHEG method EPJC 71 (2011) 1547 1009.2450
53 H. B. Hartanto, B. Jager, L. Reina, and D. Wackeroth Higgs boson production in association with top quarks in the POWHEG BOX PRD 91 (2015) 094003 1501.04498
54 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
55 P. Artoisenet, R. Frederix, O. Mattelaer, and R. Rietkerk Automatic spin-entangled decays of heavy resonances in Monte Carlo simulations JHEP 03 (2013) 015 1212.3460
56 NNPDF Collaboration Parton distributions for the LHC Run II JHEP 04 (2015) 040 1410.8849
57 T. Sjöstrand et al. An introduction to PYTHIA 8.2 Comput. Phys. Commun. 191 (2015) 159 1410.3012
58 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
59 GEANT4 Collaboration GEANT 4---a simulation toolkit NIM A 506 (2003) 250
60 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
61 N. Kidonakis Top quark production 1311.0283
62 Y. Li and F. Petriello Combining QCD and electroweak corrections to dilepton production in the framework of the FEWZ simulation code PRD 86 (2012) 094034 1208.5967
63 T. Ježo and P. Nason On the treatment of resonances in next-to-leading order calculations matched to a parton shower JHEP 12 (2015) 065 1509.09071
64 F. Buccioni et al. OpenLoops 2 EPJC 79 (2019) 866 1907.13071
65 Bahr, M. and others Herwig++ Physics and Manual EPJC 58 (2008) 639 0803.0883
66 Bellm, Johannes and others Herwig 7.0/Herwig++ 3.0 release note EPJC 76 (2016) 196 1512.01178
67 CMS Collaboration Development and validation of HERWIG 7 tunes from CMS underlying-event measurements EPJC 81 (2021) 312 CMS-GEN-19-001
2011.03422
68 T. Gleisberg et al. Event generation with SHERPA 1.1 JHEP 02 (2009) 007 0811.4622
69 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
70 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
CDS
71 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
72 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
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 Jet energy scale and resolution in the CMS experiment in pp collisions at 8 TeV JINST 12 (2017) P02014 CMS-JME-13-004
1607.03663
76 CMS Collaboration Pileup mitigation at CMS in 13 TeV data JINST 15 (2020) P09018 CMS-JME-18-001
2003.00503
77 CMS Collaboration Identification of heavy-flavour jets with the CMS detector in pp collisions at 13 TeV JINST 13 (2018) P05011 CMS-BTV-16-002
1712.07158
78 E. Bols et al. Jet flavour classification using DeepJet JINST 1 (2020) 5 2008.10519
79 M. Cacciari and G. P. Salam Pileup subtraction using jet areas PLB 659 (2008) 119 0707.1378
80 CMS Collaboration Simulation of the Silicon Strip Tracker pre-amplifier in early 2016 data CMS Detector Performance Note CMS-DP-2020-045
CDS
81 Y. Le Cun et al. Handwritten digit recognition with a back-propagation network in Proceedings of the 2nd International Conference on Neural Information Processing Systems, NIPS'89
82 S. Hochreiter and J. Schmidhuber Long short-term memory Neural Comput. 9 (1997) 1735
83 F. Chollet et al. Keras link
84 A. Martìn et al. TensorFlow: Large-scale machine learning on heterogeneous systems link
85 D. P. Kingma and J. Ba Adam: A method for stochastic optimization link
86 N. Srivastava et al. Dropout: A simple way to prevent neural networks from overfitting J. Mach. Learn. Res. 15 (2014) 1929
87 Smoothfit team Smoofit: a package for smooth binned likelihood fits link
88 J. Bradbury et al. JAX: composable transformations of Python+NumPy programs (2018)
link
89 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
90 CMS Collaboration CMS luminosity measurement for the 2017 data-taking period at $ \sqrt{s} $ = 13 TeV CMS Physics Analysis Summary, 2018
link
CMS-PAS-LUM-17-004
91 CMS Collaboration CMS luminosity measurement for the 2018 data-taking period at $ \sqrt{s} $ = 13 TeV CMS Physics Analysis Summary, 2019
link
CMS-PAS-LUM-18-002
92 M. Czakon et al. Top-pair production at the LHC through NNLO QCD and NLO EW JHEP 10 (2017) 186 1705.04105
93 S. Argyropoulos and T. Sjöstrand Effects of color reconnection on $ \mathrm{t}\bar{\mathrm{t}} $ final states at the LHC JHEP 11 (2014) 1407.6653
94 J. R. Christiansen and P. Z. Skands String formation beyond leading colour JHEP 08 (2015) 1505.01681
95 ATLAS Collaboration Measurement of $ b $-quark fragmentation properties in jets using the decay $ B^{\pm} \to J/\psi K^{\pm} $ in $ pp $ collisions at $ \sqrt{s} $ = 13 TeV with the ATLAS detector JHEP 12 (2021) 131 2108.11650
96 R. J. Barlow and C. Beeston Fitting using finite Monte Carlo samples Comput. Phys. Commun. 77 (1993) 219
97 J. S. Conway Incorporating Nuisance Parameters in Likelihoods for Multisource Spectra PHYSTAT 201 (2011) 115 1103.0354
Compact Muon Solenoid
LHC, CERN