CMS-PAS-HIG-17-026

CMS-PAS-HIG-17-026
Search for $\mathrm{t\bar{t}H}$ production in the $\mathrm{H\rightarrow b\bar{b}}$ decay channel with leptonic $\mathrm{t\bar{t}}$ decays in proton-proton collisions at $\sqrt{s}=$ 13 TeV with the CMS detector
CMS Collaboration
March 2018

Abstract: A search for the associated production of a standard model Higgs boson with a top quark-antiquark pair ( $\mathrm{t\bar{t}H}$ ) in proton-proton collisions at a center-of-mass energy $\sqrt{s}=$ 13 TeV is presented. The data correspond to an integrated luminosity of 35.9 fb $^{-1}$ recorded with the CMS detector at the CERN LHC in 2016. Candidate $\mathrm{t\bar{t}H}$ events are selected that contain either one or two electrons or muons from the $\mathrm{t\bar{t}}$ decays, and are categorized according to the number of jets. Multivariate techniques are employed to categorize further the events and eventually discriminate between signal and background events. The results are characterized by an observed $\mathrm{t\bar{t}H}$ signal strength relative to the standard model cross section, $\mu = \sigma/\sigma_{\mathrm{SM}}$ , under the assumption of $m_{\mathrm{H}} =$ 125 GeV. A combined fit of multivariate discriminant distributions in all categories results in an observed (expected) upper limit of $\mu <$ 1.5 (0.9) at the 95% confidence level, and a best fit value of $\mu =$ 0.72 $\pm$ 0.45, corresponding to an observed (expected) significance of 1.6 (2.2) standard deviations.
Links: CDS record (PDF) ; inSPIRE record ; CADI line (restricted) ; These preliminary results are superseded in this paper, JHEP 03 (2019) 026. The superseded preliminary plots can be found here.

Figures & Tables	Summary	References	CMS Publications

Figures
png pdf	Figure 1: Representative leading-order Feynman diagrams for $\mathrm{ t\bar{t} H }$ production, including the subsequent decay of the Higgs boson into a b quark-antiquark pair, and the decay of the top quark-antiquark pair into final states with either one (single-lepton channel, left) or two (dilepton channel, right) electrons or muons.
png pdf	Figure 1-a: Representative leading-order Feynman diagram for $\mathrm{ t\bar{t} H }$ production, including the subsequent decay of the Higgs boson into a b quark-antiquark pair, and the decay of the top quark-antiquark pair into final states with either one (single-lepton channel, left) or two (dilepton channel, right) electrons or muons.
png pdf	Figure 1-b: Representative leading-order Feynman diagram for $\mathrm{ t\bar{t} H }$ production, including the subsequent decay of the Higgs boson into a b quark-antiquark pair, and the decay of the top quark-antiquark pair into final states with either one (single-lepton channel, left) or two (dilepton channel, right) electrons or muons.
png pdf	Figure 2: Jet (left) and b-tagged jet (right) multiplicity after the baseline event selection in the single-lepton channel. The uncertainty bands correspond to the total statistic and systematic uncertainties (excluding uncertainties that affect only the rate of the distribution) added in quadrature.
png pdf	Figure 2-a: Jet multiplicity after the baseline event selection in the single-lepton channel. The uncertainty bands correspond to the total statistic and systematic uncertainties (excluding uncertainties that affect only the rate of the distribution) added in quadrature.
png pdf	Figure 2-b: b-tagged jet multiplicity after the baseline event selection in the single-lepton channel. The uncertainty bands correspond to the total statistic and systematic uncertainties (excluding uncertainties that affect only the rate of the distribution) added in quadrature.
png pdf	Figure 3: Jet (left) and b-tagged jet (right) multiplicity after the baseline event selection in the dilepton channel. The uncertainty bands correspond to the total statistical and systematic uncertainties (excluding uncertainties that affect only the rate of the distribution) added in quadrature.
png pdf	Figure 3-a: Jet multiplicity after the baseline event selection in the dilepton channel. The uncertainty bands correspond to the total statistical and systematic uncertainties (excluding uncertainties that affect only the rate of the distribution) added in quadrature.
png pdf	Figure 3-b: b-tagged jet multiplicity after the baseline event selection in the dilepton channel. The uncertainty bands correspond to the total statistical and systematic uncertainties (excluding uncertainties that affect only the rate of the distribution) added in quadrature.
png pdf	Figure 4: Final discriminant shapes in the single-lepton (SL) channel before the fit to data: DNN discriminant in the jet-process categories with $\geq$ 6 jets- $\mathrm{ t\bar{t} H }$ (top left); 5 jets- $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ (top right); 4 jets- $\mathrm{ t \bar{t} } + \text{lf}$ (bottom left) and $\geq$ 6 jets- $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ (bottom right). The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model. The distributions observed in data (markers) are overlayed.
png pdf	Figure 4-a: Final discriminant shape in the single-lepton (SL) channel before the fit to data: DNN discriminant in the jet-process category with $\geq$ 6 jets- $\mathrm{ t\bar{t} H }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model. The distributions observed in data (markers) are overlayed.
png pdf	Figure 4-b: Final discriminant shape in the single-lepton (SL) channel before the fit to data: DNN discriminant in the jet-process category with 5 jets- $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model. The distributions observed in data (markers) are overlayed.
png pdf	Figure 4-c: Final discriminant shape in the single-lepton (SL) channel before the fit to data: DNN discriminant in the jet-process category with 4 jets- $\mathrm{ t \bar{t} } + \text{lf}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model. The distributions observed in data (markers) are overlayed.
png pdf	Figure 4-d: Final discriminant shape in the single-lepton (SL) channel before the fit to data: DNN discriminant in the jet-process category with $\geq$ 6 jets- $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model. The distributions observed in data (markers) are overlayed.
png pdf	Figure 5: Final discriminant shapes in the dilepton (DL) channel before the fit to data: BDT discriminant in the analysis category with ( $\geq$ 4 jets, 3 b-tags) (top row) and MEM discriminant in the analysis categories with ( $\geq$ 4 jets, $\geq$ 4 b-tags) (bottom row) with low (left) and high (right) BDT output. The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model. The distributions observed in data (markers) are overlayed.
png pdf	Figure 5-a: Final discriminant shape in the dilepton (DL) channel before the fit to data: BDT discriminant in the analysis category with ( $\geq$ 4 jets, 3 b-tags). The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model. The distributions observed in data (markers) are overlayed.
png pdf	Figure 5-b: Final discriminant shape in the dilepton (DL) channel before the fit to data: MEM discriminant in the analysis categories with ( $\geq$ 4 jets, $\geq$ 4 b-tags) with low BDT output. The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model. The distributions observed in data (markers) are overlayed.
png pdf	Figure 5-c: Final discriminant shape in the dilepton (DL) channel before the fit to data: MEM discriminant in the analysis categories with ( $\geq$ 4 jets, $\geq$ 4 b-tags) with high BDT output. The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model. The distributions observed in data (markers) are overlayed.
png pdf	Figure 6: Final discriminant shapes in the single-lepton (SL) channel after the fit to data: DNN discriminant in the jet-process categories with $\geq$ 6 jets- $\mathrm{ t\bar{t} H }$ (top left); 5 jets- $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ (top right); 4 jets- $\mathrm{ t \bar{t} } + \text{lf}$ (bottom left) and $\geq$ 6 jets- $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ (bottom right). The error bands include the total uncertainty after the fit to data. The distributions observed in data (markers) are overlayed.
png pdf	Figure 6-a: Final discriminant shape in the single-lepton (SL) channel after the fit to data: DNN discriminant in the jet-process categories with $\geq$ 6 jets- $\mathrm{ t\bar{t} H }$ . The error bands include the total uncertainty after the fit to data. The distribution observed in data (markers) is overlayed.
png pdf	Figure 6-b: Final discriminant shape in the single-lepton (SL) channel after the fit to data: DNN discriminant in the jet-process categories with 5 jets- $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ . The error bands include the total uncertainty after the fit to data. The distribution observed in data (markers) is overlayed.
png pdf	Figure 6-c: Final discriminant shape in the single-lepton (SL) channel after the fit to data: DNN discriminant in the jet-process categories with 4 jets- $\mathrm{ t \bar{t} } + \text{lf}$ . The error bands include the total uncertainty after the fit to data. The distribution observed in data (markers) is overlayed.
png pdf	Figure 6-d: Final discriminant shape in the single-lepton (SL) channel after the fit to data: DNN discriminant in the jet-process categories with $\geq$ 6 jets- $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ . The error bands include the total uncertainty after the fit to data. The distribution observed in data (markers) is overlayed.
png pdf	Figure 7: Final discriminant shapes in the dilepton (DL) channel after the fit to data: BDT discriminant in the analysis category with ( $\geq$ 4 jets, 3 b-tags) (top row) and MEM discriminant in the analysis categories with ( $\geq$ 4 jets, $\geq$ 4 b-tags) (bottom row) with low (left) and high (right) BDT output. The error bands include the total uncertainty after the fit to data. The distributions observed in data (markers) are overlayed.
png pdf	Figure 7-a: Final discriminant shapes in the dilepton (DL) channel after the fit to data: BDT discriminant in the analysis category with ( $\geq$ 4 jets, 3 b-tags). The error bands include the total uncertainty after the fit to data. The distribution observed in data (markers) is overlayed.
png pdf	Figure 7-b: Final discriminant shapes in the dilepton (DL) channel after the fit to data: MEM discriminant in the analysis category with ( $\geq$ 4 jets, $\geq$ 4 b-tags) with low BDT output. The error bands include the total uncertainty after the fit to data. The distribution observed in data (markers) is overlayed.
png pdf	Figure 7-c: Final discriminant shapes in the dilepton (DL) channel after the fit to data: MEM discriminant in the analysis category with ( $\geq$ 4 jets, $\geq$ 4 b-tags) with high BDT output. The error bands include the total uncertainty after the fit to data. The distribution observed in data (markers) is overlayed.
png pdf	Figure 8: Bins of the final discriminants as used in the fit (left), re-ordered by the pre-fit expected signal-to-background ratio (S/B). Each of the shown bins includes multiple bins of the final discriminants with similar S/B. The background is shown post-fit (S+B), with the fitted signal in cyan compared to the expectation for the SM Higgs boson ( $\mu =$ 1) in red. Best-fit values of the signal strength modifiers $\mu$ (right) with their $\pm$ 1 $\sigma$ confidence intervals (outer error bar), also split into their statistical (inner error bar) and systematic components.
png pdf	Figure 8-a: Bins of the final discriminants as used in the fit, re-ordered by the pre-fit expected signal-to-background ratio (S/B). Each of the shown bins includes multiple bins of the final discriminants with similar S/B. The background is shown post-fit (S+B), with the fitted signal in cyan compared to the expectation for the SM Higgs boson ( $\mu =$ 1) in red.
png pdf	Figure 8-b: Best-fit values of the signal strength modifiers $\mu$ with their $\pm$ 1 $\sigma$ confidence intervals (outer error bar), also split into their statistical (inner error bar) and systematic components.
png pdf	Figure 9: Median expected (dashed line) and observed (markers) 95% CL upper limits on $\mu$ . The expected limits are displayed together with $\pm$ 1 $\sigma$ and $\pm$ 2 $\sigma$ confidence intervals (green and yellow bands). Also shown is the observed limit that is expected in case a SM $\mathrm{ t\bar{t} H }$ signal ( $\mu =1$ ) is present in the data (solid red line).
png pdf	Figure 10: Final discriminant (DNN) shapes in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and (from top left to bottom right) $\mathrm{ t\bar{t} H }$ , $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ , $\mathrm{ t \bar{t} } + \text{2b}$ , $\mathrm{ t \bar{t} } + \text{b}$ , $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ and $\mathrm{ t \bar{t} } + \text{lf}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 10-a: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and $\mathrm{ t\bar{t} H }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 10-b: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 10-c: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{2b}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 10-d: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{b}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 10-e: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 10-f: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{lf}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 11: Final discriminant (DNN) shapes in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and (from top left to bottom right) $\mathrm{ t\bar{t} H }$ , $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ , $\mathrm{ t \bar{t} } + \text{2b}$ , $\mathrm{ t \bar{t} } + \text{b}$ , $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ and $\mathrm{ t \bar{t} } + \text{lf}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 11-a: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t\bar{t} H }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 11-b: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 11-c: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{2b}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 11-d: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{b}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 11-e: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 11-f: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{lf}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 12: Final discriminant (DNN) shapes in the single-lepton channel before the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and (from top left to bottom right) $\mathrm{ t\bar{t} H }$ , $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ , $\mathrm{ t \bar{t} } + \text{2b}$ , $\mathrm{ t \bar{t} } + \text{b}$ , $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ and $\mathrm{ t \bar{t} } + \text{lf}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 12-a: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t\bar{t} H }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 12-b: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 12-c: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{2b}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 12-d: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{b}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 12-e: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 12-f: Final discriminant (DNN) shape in the single-lepton channel before the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{lf}$ . The expected background contributions (filled histograms) are stacked, and the expected signal distribution (line) for a Higgs-boson mass of ${m_{\mathrm{H}}} =$ 125 GeV is superimposed. Each contribution is normalized to an integrated luminosity of 35.9 fb $^{-1}$ , and the signal distribution is additionally scaled by a factor of 15 for better readability. The error bands include the total uncertainty of the fit model.
png pdf	Figure 13: Final discriminant (DNN) shapes in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and (from top left to bottom right) $\mathrm{ t\bar{t} H }$ , $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ , $\mathrm{ t \bar{t} } + \text{2b}$ , $\mathrm{ t \bar{t} } + \text{b}$ , $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ and $\mathrm{ t \bar{t} } + \text{lf}$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 13-a: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and $\mathrm{ t\bar{t} H }$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 13-b: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 13-c: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{2b}$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 13-d: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{b}$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 13-e: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and \mathrm{ t \bar{t} } + \mathrm{ c \bar{c} } $. The error bands include the total uncertainty after the fit to data.
png pdf	Figure 13-f: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (4 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{lf}$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 14: Final discriminant (DNN) shapes in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and (from top left to bottom right) $\mathrm{ t\bar{t} H }$ , $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ , $\mathrm{ t \bar{t} } + \text{2b}$ , $\mathrm{ t \bar{t} } + \text{b}$ , $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ and $\mathrm{ t \bar{t} } + \text{lf}$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 14-a: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t\bar{t} H }$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 14-b: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 14-c: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{2b}$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 14-d: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{b}$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 14-e: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 14-f: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with (5 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{lf}$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 15: Final discriminant (DNN) shapes in the single-lepton channel after the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and (from top left to bottom right) $\mathrm{ t\bar{t} H }$ , $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ , $\mathrm{ t \bar{t} } + \text{2b}$ , $\mathrm{ t \bar{t} } + \text{b}$ , $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ and $\mathrm{ t \bar{t} } + \text{lf}$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 15-a: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t\bar{t} H }$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 15-b: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \mathrm{ b \bar{b} }$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 15-c: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{2b}$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 15-d: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{b}$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 15-e: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \mathrm{ c \bar{c} }$ . The error bands include the total uncertainty after the fit to data.
png pdf	Figure 15-f: Final discriminant (DNN) shape in the single-lepton channel after the fit to data, in the ''jet-process'' categories with ( $\geq$ 6 jets, 3 b-tags) and $\mathrm{ t \bar{t} } + \text{lf}$ . The error bands include the total uncertainty after the fit to data.

Tables
png pdf	Table 1: Event yields observed in data and predicted by the simulation after the baseline selection requirements in the single-lepton (SL) channel (at least four jets, at least two of which are b tagged) and dilepton (DL) channel (at least three jets, at least one of which is b tagged). The quoted uncertainties are statistical only.
png pdf	Table 2: Systematic uncertainties considered in the analysis.
png pdf	Table 3: Observed and expected event yields per jet-process category in the single-lepton channel with 4 jets and at least 3 b-tags, prior to the fit to data (after the fit to data). The quoted uncertainties denote the total statistical and systematic uncertainty.
png pdf	Table 4: Observed and expected event yields per jet-process category in the single-lepton channel with 5 jets and at least 3 b-tags, prior to the fit to data (after the fit to data). The quoted uncertainties denote the total statistical and systematic uncertainty.
png pdf	Table 5: Observed and expected event yields per jet-process category in the single-lepton channel with at least 6 jets and at least 3 b-tags, prior to the fit to data (after the fit to data). The quoted uncertainties denote the total statistical and systematic uncertainty.
png pdf	Table 6: Observed and expected event yields per jet-tag category in the dilepton channel, prior to the fit to data (after the fit to data). The quoted uncertainties denote the total statistical and systematic uncertainty.
png pdf	Table 7: Best-fit value of the signal strength modifier $\mu$ and the observed and median expected 95% CL upper limits in the dilepton and the single-lepton channels as well as the combined results. The one standard deviation ( $\pm$ 1 $\sigma$ ) confidence intervals of the expected limit and the best-fit value are also quoted, split into the statistical and systematic components in the latter case. Expected limits are calculated with the asymptotic method [83].
png pdf	Table 8: Contributions of different sources of uncertainties to the result for the fit to the data (observed) and to the expectation from simulation (expected). The quoted uncertainties $\sigma _{\mu}$ on $\mu$ are obtained by fixing the listed uncertainties in the fit, and subtracting the obtained result in quadrature from the result of the full fit. The quadratic sum of the contributions is different from the total uncertainty due to correlations between the nuisance parameters.

Summary

A search for the associated production of a Higgs boson and a top quark-antiquark pair is performed using 35.9 fb

$^{-1}$ value of pp collision data recorded with the CMS detector at a center-of-mass energy of 13 TeV in 2016. Candidate events are selected in final states compatible with the Higgs boson decay

$\text{H} \rightarrow \text{bb}$ and the single-lepton and dilepton decay channels of the

$\mathrm{t\bar{t}}$ pair. Selected events are split into mutually exclusive categories according to their

$\mathrm{t\bar{t}}$ decay channel and jet content. In each category a powerful discriminant is constructed to separate the

$\mathrm{ t\bar{t} H }$ signal from the

$\mathrm{t\bar{t}}$ -dominated background, based on different multivariate analysis techniques (boosted decision trees, matrix element method, deep neural networks). An observed (expected) upper limit on the

$\mathrm{ t\bar{t} H }$ production cross section relative to the SM expectations of

$\mu =$ 1.5 (0.9) at the 95% confidence level is obtained. The best-fit value of

$\mu$ is 0.72

$\pm$ 0.24 (stat)

$\pm$ 0.38 (syst). These results correspond to a significance of 1.6 standard deviations above the background-only hypothesis.

References
1	ATLAS Collaboration	Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC	PLB 716 (2012), no. 1, 1--29	1207.7214
2	CMS Collaboration	Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC	PLB 716 (2012), no. 1, 30--61	CMS-HIG-12-028 1207.7235
3	CMS Collaboration	Evidence for the direct decay of the 125 GeV Higgs boson to fermions	NP 10 (2014), no. 5, 557--560	CMS-HIG-13-033 1401.6527
4	ATLAS Collaboration	Evidence for the Higgs-boson Yukawa coupling to tau leptons with the ATLAS detector	JHEP 04 (2015) 117	1501.04943
5	CMS Collaboration	Observation of the Higgs boson decay to a pair of tau leptons		CMS-HIG-16-043 1708.00373
6	ATLAS Collaboration	Measurements of Higgs boson production and couplings in diboson final states with the ATLAS detector at the LHC	PLB 726 (2013), no. 1-3, 88--119	1307.1427
7	CMS Collaboration	Precise determination of the mass of the Higgs boson and tests of compatibility of its couplings with the standard model predictions using proton collisions at 7 and 8 TeV	EPJC. 75 (2015), no. 5, 212	CMS-HIG-14-009 1412.8662
8	ATLAS Collaboration	Evidence for the spin-0 nature of the Higgs boson using ATLAS data	PLB 726 (2013), no. 1-3, 120--144	1307.1432
9	CMS Collaboration	Constraints on the spin-parity and anomalous HVV couplings of the Higgs boson in proton collisions at 7 and 8 TeV	PRD 92 (2015), no. 1, 012004	CMS-HIG-14-018 1411.3441
10	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	Submitted to Eur.\ Phys.\ J.\ C. (2017)	CMS-TOP-17-009 1710.10614
11	LHC Higgs Cross Section Working Group Collaboration	Handbook of LHC Higgs Cross Sections: 4. Deciphering the Nature of the Higgs Sector		1610.07922
12	G. Burdman, M. Perelstein, and A. Pierce	Large Hadron Collider tests of a little Higgs model	PRL 90 (2003), no. 24, 241802	hep-ph/0212228
13	T. Han, H. E. Logan, B. McElrath, and L.-T. Wang	Phenomenology of the little Higgs model	PRD 67 (2003), no. 9, 095004	hep-ph/0301040
14	M. Perelstein, M. E. Peskin, and A. Pierce	Top quarks and electroweak symmetry breaking in little Higgs models	PRD 69 (2004), no. 7, 075002	hep-ph/0310039
15	H.-C. Cheng, I. Low, and L.-T. Wang	Top partners in little Higgs theories with T-parity	PRD 74 (2006), no. 5, 055001	hep-ph/0510225
16	H.-C. Cheng, B. A. Dobrescu, and C. T. Hill	Electroweak symmetry breaking and extra dimensions	NPB. 589 (2000), no. 1-3, 249--268	hep-ph/9912343
17	M. Carena, E. Ponton, J. Santiago, and C. E. M. Wagner	Light Kaluza Klein States in Randall-Sundrum Models with Custodial SU(2)	NPB. 759 (2006), no. 1-2, 202--227	hep-ph/0607106
18	R. Contino, L. Da Rold, and A. Pomarol	Light custodians in natural composite Higgs models	PRD 75 (2007), no. 5, 055014	hep-ph/0612048
19	G. Burdman and L. Da Rold	Electroweak Symmetry Breaking from a Holographic Fourth Generation	JHEP 12 (2007) 086	0710.0623
20	C. T. Hill	Topcolor: Top quark condensation in a gauge extension of the standard model	PLB 266 (1991), no. 3, 419--424
21	A. Carmona, M. Chala, and J. Santiago	New Higgs Production Mechanism in Composite Higgs Models	JHEP 07 (2012) 049	1205.2378
22	CMS Collaboration	Search for the associated production of the Higgs boson with a top-quark pair	JHEP 09 (2014) 087	CMS-HIG-13-029 1408.1682
23	ATLAS Collaboration	Search for the associated production of the Higgs boson with a top quark pair in multilepton final states with the ATLAS detector	PLB 749 (2015) 519--541	1506.05988
24	CMS Collaboration	Search for a standard model Higgs boson produced in association with a top-quark pair and decaying to bottom quarks using a matrix element method	EPJC. 75 (2015), no. 6, 251	CMS-HIG-14-010 1502.02485
25	ATLAS Collaboration	Search for the Standard Model Higgs boson produced in association with top quarks and decaying into $b\bar{b}$ in pp collisions at $\sqrt{s} =$ 8 TeV with the ATLAS detector	EPJC. 75 (2015), no. 7, 349	1503.05066
26	LHC Higgs Cross Section Working Group Collaboration	Handbook of LHC Higgs Cross Sections: 1. Inclusive Observables		1101.0593
27	ATLAS Collaboration	Evidence for the associated production of the Higgs boson and a top quark pair with the ATLAS detector	Submitted to: PRD (2017)	1712.08891
28	ATLAS Collaboration	Search for the Standard Model Higgs boson produced in association with top quarks and decaying into a $b\bar{b}$ pair in $pp$ collisions at $\sqrt{s} =$ 13 TeV with the ATLAS detector	Submitted to: PRD (2017)	1712.08895
29	T. J. Hastie, R. J. Tibshirani, and J. H. Friedman	The elements of statistical learning : data mining, inference, and prediction	Springer series in statistics. Springer, New York, NY, 2. ed., corr. at 10. print. edition, 2013 ISBN 978-0-387-84857-0
30	P. C. Bhat	Multivariate analysis methods in particle physics	Annual Review of Nuclear and Particle Science 61 (2011), no. 1, 281
31	A. Hocker et al.	TMVA: Toolkit for Multivariate Data Analysis	PoS ACAT (2007) 040	physics/0703039
32	K. Kondo	Dynamical Likelihood Method for Reconstruction of Events With Missing Momentum. 1: Method and Toy Models	J. Phys. Soc. Jap. 57 (1988) 4126--4140
33	D0 Collaboration	A precision measurement of the mass of the top quark	Nature 429 (2004) 638--642	hep-ex/0406031
34	CMS Collaboration	The CMS experiment at the CERN LHC	JINST 3 (2008), no. 8, S08004	CMS-00-001
35	CMS Collaboration	The CMS trigger system	JINST 12 (2017) P01020	CMS-TRG-12-001 1609.02366
36	GEANT4 Collaboration	GEANT4---a simulation toolkit	NIMA 506 (2003), no. 3, 250
37	S. Frixione, P. Nason, and C. Oleari	Matching NLO QCD computations with parton shower simulations: the POWHEG method	JHEP 11 (2007) 070	0709.2092
38	E. Re	Single-top Wt-channel production matched with parton showers using the POWHEG method	EPJC 71 (2011) 1547	1009.2450
39	NNPDF Collaboration	Parton distributions for the LHC Run II	JHEP 04 (2015) 040	1410.8849
40	T. Sjostrand et al.	An introduction to PYTHIA 8.2	CPC 191 (2015) 159	1410.3012
41	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
42	R. Frederix and S. Frixione	Merging meets matching in MC@NLO	JHEP 12 (2012) 061	1209.6215
43	CMS Collaboration	Underlying event tunes and double parton scattering	CDS
44	P. Skands, S. Carrazza, and J. Rojo	Tuning PYTHIA 8.1: the Monash 2013 Tune	EPJC 74 (2014), no. 8	1404.5630
45	N. Kidonakis	Two-loop soft anomalous dimensions for single top quark associated production with $\mathrm{W^-}$ or $\mathrm{H^-}$	PRD 82 (2010), no. 5, 054018	1005.4451
46	M. Aliev et al.	HATHOR: HAdronic Top and Heavy quarks crOss section calculatoR	CPC 182 (2011) 1034--1046	1007.1327
47	P. Kant et al.	HatHor for single top-quark production: Updated predictions and uncertainty estimates for single top-quark production in hadronic collisions	CPC 191 (2015) 74--89	1406.4403
48	F. Maltoni, D. Pagani, and I. Tsinikos	Associated production of a top-quark pair with vector bosons at NLO in QCD: impact on $t \bar{t} H$ searches at the LHC		1507.05640
49	J. M. Campbell, R. K. Ellis, and C. Williams	Vector boson pair production at the LHC	JHEP 07 (2011) 018	1105.0020
50	W. Beenakker et al.	Higgs radiation off top quarks at the Tevatron and the LHC	PRL 87 (2001), no. 20, 201805	hep-ph/0107081
51	W. Beenakker et al.	NLO QCD corrections to $\mathrm{t\bar{t}}\mathrm{H}$ production in hadron collisions	NPB 653 (2003), no. 1-2, 151	hep-ph/0211352
52	S. Dawson, L. H. Orr, L. Reina, and D. Wackeroth	Associated top quark Higgs boson production at the LHC	PRD 67 (2003), no. 7, 071503	hep-ph/0211438
53	S. Dawson et al.	Associated Higgs production with top quarks at the large hadron collider: NLO QCD corrections	PRD 68 (2003), no. 3, 034022	hep-ph/0305087
54	A. Djouadi, J. Kalinowski, and M. Spira	HDECAY: A program for Higgs boson decays in the standard model and its supersymmetric extension	CPC 108 (1998), no. 1, 56	hep-ph/9704448
55	A. Djouadi, M. M. Muhlleitner, and M. Spira	Decays of supersymmetric particles: The Program SUSY-HIT (SUspect-SdecaY-Hdecay-InTerface)	Acta Phys. Polon. B 38 (2007) 635	hep-ph/0609292
56	A. Bredenstein, A. Denner, S. Dittmaier, and M. M. Weber	Precise predictions for the Higgs-boson decay $\mathrm{H} \to \mathrm{W}\mathrm{W}/\mathrm{Z}\mathrm{Z} \to 4$ leptons	PRD 74 (2006), no. 1, 013004	hep-ph/0604011
57	A. Bredenstein, A. Denner, S. Dittmaier, and M. M. Weber	Radiative corrections to the semileptonic and hadronic Higgs-boson decays $\mathrm{H} \to \mathrm{W}\mathrm{W}/\mathrm{Z}\mathrm{Z} \to 4$ fermions	JHEP 02 (2007) 080	hep-ph/0611234
58	M. Cacciari et al.	Top-pair production at hadron colliders with next-to-next-to-leading logarithmic soft-gluon resummation	PLB 710 (2012), no. 4-5, 612	1111.5869
59	P. Baernreuther et al.	Percent level precision physics at the tevatron: First genuine nnlo qcd corrections to $q\bar{q}\rightarrow t\bar{t} +x$	PRL 109 (2012), no. 13, 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 and A. Mitov	Nnlo corrections to top-pair production at hadron colliders: the quark-gluon reaction	JHEP 01 (2013) 080	1210.6832
62	M. Beneke et al.	Hadronic top-quark pair production with nnll threshold resummation	NPB 855 (2012), no. 3, 695	1109.1536
63	M. Czakon, P. Fiedler, and A. Mitov	Total Top-Quark Pair-Production Cross Section at Hadron Colliders Through $O({\alpha_S}^4)$	PRL 110 (2013), no. 25, 252004	1303.6254
64	M. Czakon and A. Mitov	Top++: A Program for the Calculation of the Top-Pair Cross-Section at Hadron Colliders	CPC 185 (2014), no. 11	1112.5675
65	CMS Collaboration	Identification of heavy-flavour jets with the CMS detector in pp collisions at 13 TeV	Submitted to: JINST (2017)	CMS-BTV-16-002 1712.07158
66	CMS Collaboration	Particle-flow reconstruction and global event description with the CMS detector	JINST 12 (2017) P10003	CMS-PRF-14-001 1706.04965
67	M. Cacciari, G. P. Salam, and G. Soyez	The anti- $k_t$ jet clustering algorithm	JHEP 04 (2008) 063	0802.1189
68	M. Cacciari, G. P. Salam, and G. Soyez	FastJet User Manual	EPJC72 (2012) 1896	1111.6097
69	M. Cacciari, G. P. Salam, and G. Soyez	The catchment area of jets	JHEP 04 (2008) 005	0802.1188
70	CMS Collaboration	Jet energy scale and resolution in the CMS experiment in pp collisions at 8 TeV	JINST 12 (2017), no. 02, P02014	CMS-JME-13-004 1607.03663
71	J. H. Friedman	Stochastic gradient boosting	Computational Statistics \& Data Analysis 38 (2002), no. 4, 367, . Nonlinear Methods and Data Mining
72	J. Kennedy and R. Eberhart	Particle swarm optimization	in Proceedings of the IEEE International Conference on neural networks, volume 4, pp. 1942--1948 Nov
73	K. El Morabit	A study of the multivariate analysis of higgs boson production in association with a top quark-antiquark pair in the boosted regime at the cms experiment	ms, Karlsruher Institut f\"ur Technologie (KIT), 2015 Karlsruher Institut f\"ur Technologie (KIT), Masterarbeit
74	CMS Collaboration	CMS Luminosity Measurements for the 2016 Data Taking Period	CMS-PAS-LUM-17-001	CMS-PAS-LUM-17-001
75	CMS Collaboration	Measurements of $t\bar{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	PLB776 (2018) 355--378	CMS-TOP-16-010 1705.10141
76	G. Bevilacqua, M. V. Garzelli, and A. Kardos	$t\bar{t}b\bar{b}$ hadroproduction with massive bottom quarks with PowHel		1709.06915
77	T. Jezo, J. M. Lindert, N. Moretti, and S. Pozzorini	New NLOPS predictions for $\boldsymbol{t\bar{t}+b}$ -jet production at the LHC		1802.00426
78	T. Gleisberg et al.	Event generation with SHERPA 1.1	JHEP 02 (2009) 007	0811.4622
79	F. Cascioli, P. Maierhofer, and S. Pozzorini	Scattering Amplitudes with Open Loops	PRL 108 (2012) 111601	1111.5206
80	CMS Collaboration	Investigations of the impact of the parton shower tuning in Pythia 8 in the modelling of $\mathrm{t\overline{t}}$ at $\sqrt{s}=$ 8 and 13 TeV
81	R. J. Barlow and C. Beeston	Fitting using finite Monte Carlo samples	CPC 77 (1993), no. 2, 219--228
82	J. S. Conway	Incorporating Nuisance Parameters in Likelihoods for Multisource Spectra	in Proceedings, PHYSTAT 2011 Workshop on Statistical Issues Related to Discovery Claims in Search Experiments and Unfolding, CERN,Geneva, Switzerland 17-20 January 2011 2011	1103.0354
83	G. Cowan, K. Cranmer, E. Gross, and O. Vitells	Asymptotic formulae for likelihood-based tests of new physics	EPJC71 (2011) 1554	1007.1727
84	A. Read	Modified frequentist analysis of search results (the ${CL}_s$ method)	CERN-OPEN-2000-005, CERN
85	T. Junk	Confidence level computation for combining searches with small statistics	NIMA 434 (1999), no. 2-3, 435	hep-ex/9902006
86	B. Efron	The Jackknife, the Bootstrap and Other Resampling Plans	Society for Industrial and Applied Mathematics