CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-TOP-20-008
A profile likelihood approach to measure the top quark mass in the lepton+jets channel at $\sqrt{s}= $ 13 TeV
Abstract: The mass of the top quark is measured in 36 fb$^{-1}$ of LHC proton-proton collision data collected with the CMS detector at $\sqrt{s}= $ 13 TeV. The measurement uses a sample of $\mathrm{t\bar{t}}$ events containing one isolated muon or electron and at least four jets in the final state. For each event, the mass is reconstructed from a kinematic fit of the decay products to a $\mathrm{t\bar{t}}$ hypothesis. A profile likelihood method is applied to up to five observables to extract the top quark mass. The impact of systematic uncertainty sources is reduced by including their effect as nuisance parameters in the likelihood. The top quark mass is measured to be 171.77 $\pm$ 0.38 GeV.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
The top quark mass distribution before (left) and after (right) the $ {P_\text {gof}} \ge $ 0.2 selection and the kinematic fit. For the simulated $\mathrm{t\bar{t}}$ events, the jet-parton assignments are classified as correct, wrong, and unmatched permutations as described in the text. The uncertainty bands contain statistical uncertainties in the simulation, normalization uncertainties due to luminosity and cross-section, and all weight-based uncertainties. A value of $ {m_{\mathrm{t}}^\text {gen}} = $ 172.5 GeV is used in the simulation.

png pdf
Figure 1-a:
The top quark mass distribution before (left) and after (right) the $ {P_\text {gof}} \ge $ 0.2 selection and the kinematic fit. For the simulated $\mathrm{t\bar{t}}$ events, the jet-parton assignments are classified as correct, wrong, and unmatched permutations as described in the text. The uncertainty bands contain statistical uncertainties in the simulation, normalization uncertainties due to luminosity and cross-section, and all weight-based uncertainties. A value of $ {m_{\mathrm{t}}^\text {gen}} = $ 172.5 GeV is used in the simulation.

png pdf
Figure 1-b:
The top quark mass distribution before (left) and after (right) the $ {P_\text {gof}} \ge $ 0.2 selection and the kinematic fit. For the simulated $\mathrm{t\bar{t}}$ events, the jet-parton assignments are classified as correct, wrong, and unmatched permutations as described in the text. The uncertainty bands contain statistical uncertainties in the simulation, normalization uncertainties due to luminosity and cross-section, and all weight-based uncertainties. A value of $ {m_{\mathrm{t}}^\text {gen}} = $ 172.5 GeV is used in the simulation.

png pdf
Figure 2:
The distributions of the reconstructed W boson mass for the $ {P_\text {gof}} \ge $ 0.2 category (left) and of the invariant mass of the lepton and the jet assigned to the semileptonic decaying top quark for the category $ {P_\text {gof}} < $ 0.2 category (right). The uncertainty bands contain statistical uncertainties in the simulation, normalization uncertainties due to luminosity and cross-section, and all weight-based uncertainties. A value of $ {m_{\mathrm{t}}^\text {gen}} = $ 172.5 GeV is used in the simulation.

png pdf
Figure 2-a:
The distributions of the reconstructed W boson mass for the $ {P_\text {gof}} \ge $ 0.2 category (left) and of the invariant mass of the lepton and the jet assigned to the semileptonic decaying top quark for the category $ {P_\text {gof}} < $ 0.2 category (right). The uncertainty bands contain statistical uncertainties in the simulation, normalization uncertainties due to luminosity and cross-section, and all weight-based uncertainties. A value of $ {m_{\mathrm{t}}^\text {gen}} = $ 172.5 GeV is used in the simulation.

png pdf
Figure 2-b:
The distributions of the reconstructed W boson mass for the $ {P_\text {gof}} \ge $ 0.2 category (left) and of the invariant mass of the lepton and the jet assigned to the semileptonic decaying top quark for the category $ {P_\text {gof}} < $ 0.2 category (right). The uncertainty bands contain statistical uncertainties in the simulation, normalization uncertainties due to luminosity and cross-section, and all weight-based uncertainties. A value of $ {m_{\mathrm{t}}^\text {gen}} = $ 172.5 GeV is used in the simulation.

png pdf
Figure 3:
The distributions of $ {m_{\ell \mathrm{b}}^\text {reco}} / {m_{\mathrm{t}}^\text {fit}} $, the invariant mass of the lepton and the jet assigned to the semileptonic decaying top quark divided by the top quark mass form the kinematic fit, (left) and of ${R_{\mathrm{b} \mathrm{q}}^\text {reco}}$, the ratio of the scalar sum of the transverse momenta of the two b-tagged jets, and the two leading untagged jets (right), both for the $ {P_\text {gof}} \ge $ 0.2 category. The uncertainty bands contain statistical uncertainties in the simulation, normalization uncertainties due to luminosity and cross-section, and all weight-based uncertainties. A value of $ {m_{\mathrm{t}}^\text {gen}} = $ 172.5 GeV is used in the simulation.

png pdf
Figure 3-a:
The distributions of $ {m_{\ell \mathrm{b}}^\text {reco}} / {m_{\mathrm{t}}^\text {fit}} $, the invariant mass of the lepton and the jet assigned to the semileptonic decaying top quark divided by the top quark mass form the kinematic fit, (left) and of ${R_{\mathrm{b} \mathrm{q}}^\text {reco}}$, the ratio of the scalar sum of the transverse momenta of the two b-tagged jets, and the two leading untagged jets (right), both for the $ {P_\text {gof}} \ge $ 0.2 category. The uncertainty bands contain statistical uncertainties in the simulation, normalization uncertainties due to luminosity and cross-section, and all weight-based uncertainties. A value of $ {m_{\mathrm{t}}^\text {gen}} = $ 172.5 GeV is used in the simulation.

png pdf
Figure 3-b:
The distributions of $ {m_{\ell \mathrm{b}}^\text {reco}} / {m_{\mathrm{t}}^\text {fit}} $, the invariant mass of the lepton and the jet assigned to the semileptonic decaying top quark divided by the top quark mass form the kinematic fit, (left) and of ${R_{\mathrm{b} \mathrm{q}}^\text {reco}}$, the ratio of the scalar sum of the transverse momenta of the two b-tagged jets, and the two leading untagged jets (right), both for the $ {P_\text {gof}} \ge $ 0.2 category. The uncertainty bands contain statistical uncertainties in the simulation, normalization uncertainties due to luminosity and cross-section, and all weight-based uncertainties. A value of $ {m_{\mathrm{t}}^\text {gen}} = $ 172.5 GeV is used in the simulation.

png pdf
Figure 4:
Left: Comparison of the expected total uncertainty on ${m_{\mathrm{t}}}$ in the combined lepton+jets channel and for the different observable-category sets defined in Table 1. Right: The difference between the measured and generated ${m_{\mathrm{t}}}$ values divided by the uncertainty reported by the fit from pseudo-experiments without (red) or with (blue) the statistical nuisance parameters ${\hat{\beta}}$ and ${\hat{\omega}}$ in the ML fit. The $\mu $ and $\sigma $ parameters of Gaussian functions (red and blue lines) fit to the histograms are included in the legend.

png pdf
Figure 4-a:
Left: Comparison of the expected total uncertainty on ${m_{\mathrm{t}}}$ in the combined lepton+jets channel and for the different observable-category sets defined in Table 1. Right: The difference between the measured and generated ${m_{\mathrm{t}}}$ values divided by the uncertainty reported by the fit from pseudo-experiments without (red) or with (blue) the statistical nuisance parameters ${\hat{\beta}}$ and ${\hat{\omega}}$ in the ML fit. The $\mu $ and $\sigma $ parameters of Gaussian functions (red and blue lines) fit to the histograms are included in the legend.

png pdf
Figure 4-b:
Left: Comparison of the expected total uncertainty on ${m_{\mathrm{t}}}$ in the combined lepton+jets channel and for the different observable-category sets defined in Table 1. Right: The difference between the measured and generated ${m_{\mathrm{t}}}$ values divided by the uncertainty reported by the fit from pseudo-experiments without (red) or with (blue) the statistical nuisance parameters ${\hat{\beta}}$ and ${\hat{\omega}}$ in the ML fit. The $\mu $ and $\sigma $ parameters of Gaussian functions (red and blue lines) fit to the histograms are included in the legend.

png pdf
Figure 5:
Left: Measurement of ${m_{\mathrm{t}}}$ in the the three different channels for the different sets of observables and categories. Right: Dependence of the 5D result on the assumed correlation between the FSR PS scales per branching in the lepton + jets channel.

png pdf
Figure 5-a:
Left: Measurement of ${m_{\mathrm{t}}}$ in the the three different channels for the different sets of observables and categories. Right: Dependence of the 5D result on the assumed correlation between the FSR PS scales per branching in the lepton + jets channel.

png pdf
Figure 5-b:
Left: Measurement of ${m_{\mathrm{t}}}$ in the the three different channels for the different sets of observables and categories. Right: Dependence of the 5D result on the assumed correlation between the FSR PS scales per branching in the lepton + jets channel.

png pdf
Figure 6:
Distribution of ${m_{\mathrm{t}}^\text {fit}}$ (top) and the additional observables (bottom) that are the input to the ML fit and their post-fit probability density functions. The green and yellow bands represent the 1-sigma and 2-sigma error bands.

png pdf
Figure 6-a:
Distribution of ${m_{\mathrm{t}}^\text {fit}}$ (top) and the additional observables (bottom) that are the input to the ML fit and their post-fit probability density functions. The green and yellow bands represent the 1-sigma and 2-sigma error bands.

png pdf
Figure 6-b:
Distribution of ${m_{\mathrm{t}}^\text {fit}}$ (top) and the additional observables (bottom) that are the input to the ML fit and their post-fit probability density functions. The green and yellow bands represent the 1-sigma and 2-sigma error bands.

png pdf
Figure 7:
Measurement of ${m_{\mathrm{t}}}$ in the combined lepton+jets channel using the 5D set of observables and categories. The left plot shows the post-fit pulls on the most important nuisance parameters and the numbers quote the post-fit uncertainty on the nuisance parameter. The right plot shows their pre-fit and post-fit impacts. The post-fit impacts include the contribution from the nuisance parameters accounting for the limited size of simulation samples (MC stat.). The average of these post-fit impacts is printed on the right. The rows are sorted by the size of the post-fit impact.
Tables

png pdf
Table 1:
The overall list of different input histograms and their inclusion in a certain histogram set. A histogram marked with "x" is included in a set (measurement).
Summary
The mass of the top quark is measured in 36 fb$^{-1}$ of LHC proton-proton collision data collected with the CMS detector at $ \sqrt{s} = $ 13 TeV. The measurement uses a sample of $\mathrm{t\bar{t}}$ events containing one isolated muon or electron and at least four jets in the final state. For each event the mass is reconstructed from a kinematic fit of the decay products to a $\mathrm{t\bar{t}}$ hypothesis. A likelihood method is applied to up to five observables to extract the top quark mass and constrain the influence of systematic effects which are included as nuisance parameters in the likelihood. Consistent results are obtained for measurements with different sets of observables and the top quark mass is measured to be 171.77 $\pm$ 0.38 GeV, including 0.04 GeV statistical uncertainty. This result denotes a considerable improvement compared to all previously published top quark mass measurements and supersedes the previously published measurement in this channel on the same data set.
References
1 CDF Collaboration Observation of top quark production in $ \bar{p}p $ collisions PRL 74 (1995) 2626 hep-ex/9503002
2 D0 Collaboration Observation of the top quark PRL 74 (1995) 2632 hep-ex/9503003
3 G. Degrassi et al. Higgs mass and vacuum stability in the standard model at NNLO JHEP 08 (2012) 1 1205.6497
4 F. Bezrukov, M. Y. Kalmykov, B. A. Kniehl, and M. Shaposhnikov Higgs boson mass and new physics JHEP 10 (2012) 140 1205.2893
5 The ALEPH, CDF, D0, DELPHI, L3, OPAL, SLD Collaborations, the LEP Electroweak Working Group, the Tevatron Electroweak Working Group, and the SLD electroweak and heavy flavour groups Precision Electroweak Measurements and Constraints on the Standard Model technical report, ALEPH 1012.2367
6 M. Baak et al. The electroweak fit of the standard model after the discovery of a new boson at the LHC EPJC 72 (2012) 2205 1209.2716
7 M. Baak et al. The global electroweak fit at NNLO and prospects for the LHC and ILC EPJC 74 (2014) 3046 1407.3792
8 Particle Data Group, P. A. Zyla et al. Review of particle physics Prog. Theor. Exp. Phys. 2020 (2020) 083C01
9 CDF and D0 Collaborations Combination of CDF and D0 results on the mass of the top quark using up 9.7 fb$^{-1}$ at the Tevatron FERMILAB-CONF-16-298-E
10 ATLAS Collaboration Measurement of the top quark mass in the $ t\bar{t}\rightarrow $ lepton+jets channel from $ \sqrt{s}= $ 8 TeV ATLAS data and combination with previous results EPJC 79 (2019), no. 4, 290 1810.01772
11 CMS Collaboration Measurement of the top quark mass using proton-proton data at $ {\sqrt{s}} = $ 7 and 8 TeV PRD 93 (2016) 072004 CMS-TOP-14-022
1509.04044
12 A. H. Hoang What Is the Top Quark Mass? Ann. Rev. Nucl. Part. Sci. 70 (2020) 225 2004.12915
13 CMS Collaboration Measurement of the top quark mass with lepton+jets final states using pp collisions at $ \sqrt{s}= $ 13 TeV EPJC 78 (2018) 891 CMS-TOP-17-007
1805.01428
14 DELPHI Collaboration Measurement of the mass and width of the W boson in $ {\rm e}^{+}{\rm e}^{-} $ collisions at $ \sqrt{s} = $ 161 -- 209 GeV EPJC 55 (2008) 1 0803.2534
15 CMS Collaboration Measurement of the top-quark mass in $ \mathrm{t\bar{t}} $ events with lepton+jets final states in pp collisions at $ \sqrt{s}= $ 7 TeV JHEP 12 (2012) 105 CMS-TOP-11-015
1209.2319
16 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA-8 tunes from underlying-event measurements EPJC 80 (2020), no. 1, 4 CMS-GEN-17-001
1903.12179
17 E. Bols et al. Jet flavour classification using DeepJet JINST 15 (2020), no. 12, P12012 2008.10519
18 CMS Collaboration Precision luminosity measurement in proton-proton collisions at $ \sqrt{s} = $ 13 TeV in 2015 and 2016 at CMS EPJC 81 (2021), no. 9, 800 CMS-LUM-17-003
2104.01927
19 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
20 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
21 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
22 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
23 T. Sjostrand, S. Mrenna, and P. Z. Skands PYTHIA 6.4 Physics and Manual JHEP 05 (2006) 026 hep-ph/0603175
24 T. Sjostrand et al. An introduction to PYTHIA 8.2 CPC 191 (2015) 159 1410.3012
25 M. Botje et al. The PDF4LHC Working Group Interim Recommendations 1101.0538
26 S. Alekhin et al. The PDF4LHC Working Group Interim Report 1101.0536
27 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017), no. 10, 663 1706.00428
28 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
29 E. Re Single-top Wt-channel production matched with parton showers using the POWHEG method EPJC 71 (2011) 1547 1009.2450
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 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
32 R. Frederix and S. Frixione Merging meets matching in MC@NLO JHEP 12 (2012) 061 1209.6215
33 P. Skands, S. Carrazza, and J. Rojo Tuning PYTHIA 8.1: the Monash 2013 Tune EPJC 74 (2014) 3024 1404.5630
34 GEANT4 Collaboration GEANT4--a simulation toolkit NIMA 506 (2003) 250
35 M. Czakon and A. Mitov Top++: A program for the calculation of the top-pair cross-section at hadron colliders CPC 185 (2014) 2930 1112.5675
36 Y. Li and F. Petriello Combining QCD and electroweak corrections to dilepton production in FEWZ PRD 86 (2012) 094034 1208.5967
37 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 1406.4403
38 N. Kidonakis NNLL threshold resummation for top-pair and single-top production Phys. Part. Nucl. 45 (2014) 714 1210.7813
39 T. Sjostrand, S. Mrenna, and P. Z. Skands A Brief Introduction to PYTHIA 8.1 CPC 178 (2008) 852 0710.3820
40 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
41 CMS Collaboration Performance of the DeepJet b tagging algorithm using 41.9/fb of data from proton-proton collisions at 13 TeV with Phase 1 CMS detector CDS
42 D0 Collaboration Direct measurement of the top quark mass at D0 PRD 58 (1998) 052001 hep-ex/9801025
43 ATLAS Collaboration Measurement of the top quark mass in the $ \mathrm{t\bar{t}}\to $ dilepton channel from $ \sqrt{s}= $ 8 TeV ATLAS data PLB 761 (2016) 350 1606.02179
44 CMS Collaboration Measurement of the $ \mathrm{t}\overline{\mathrm{t}} $ production cross section, the top quark mass, and the strong coupling constant using dilepton events in pp collisions at $ \sqrt{s} = $ 13 TeV EPJC 79 (2019), no. 5, 368 CMS-TOP-17-001
1812.10505
45 ATLAS Collaboration Measurement of the top quark mass in the $ t\bar{t}\rightarrow \text{lepton+jets} $ and $ t\bar{t}\rightarrow $ dilepton channels using $ \sqrt{s}= $ 7 TeV ATLAS data EPJC 75 (2015), no. 7, 330 1503.05427
46 CMS Collaboration Measurement of the top quark mass in the all-jets final state at $ \sqrt{s} = $ 13 TeV and combination with the lepton+jets channel EPJC 79 (2019), no. 4, 313 CMS-TOP-17-008
1812.10534
47 CMS Collaboration Measurement of the inelastic proton-proton cross section at $ \sqrt{s}= $ 13 TeV JHEP 07 (2018) 161 CMS-FSQ-15-005
1802.02613
48 CMS Collaboration Measurement of the t-channel single-top-quark production cross section and of the V$_{tb}$ CKM matrix element in pp collisions at $ \sqrt{s} = $ 8 TeV JHEP 06 (2014) 090 CMS-TOP-12-038
1403.7366
49 CMS Collaboration Cross section measurement of $ t $-channel single top quark production in pp collisions at $ \sqrt s = $ 13 TeV PLB 772 (2017) 752 CMS-TOP-16-003
1610.00678
50 CMS Collaboration Measurement of the production cross section of a W boson in association with two b jets in pp collisions at $ \sqrt{s} = $ 8 TeV Eur. Phys. J C. 77 (2017) 92 CMS-SMP-14-020
1608.07561
51 CMS Collaboration Measurements of the associated production of a Z boson and b jets in pp collisions at $ {\sqrt{s}} = $ 8 TeV EPJC 77 (2017) 751 CMS-SMP-14-010
1611.06507
52 CMS Collaboration Measurement of the WZ production cross section in pp collisions at $ \sqrt(s) = $ 13 TeV PLB 766 (2017) 268 CMS-SMP-16-002
1607.06943
53 CMS Collaboration Measurements of the pp $ \rightarrow $ ZZ production cross section and the Z $ \rightarrow $ 4$\ell $ branching fraction, and constraints on anomalous triple gauge couplings at $ \sqrt{s} = $ 13 TeV EPJC 78 (2018) 165 CMS-SMP-16-017
1709.08601
54 M. Bahr et al. Herwig++ physics and manual EPJC 58 (2008) 639 0803.0883
55 DELPHI Collaboration A study of the b-quark fragmentation function with the DELPHI detector at LEP I and an averaged distribution obtained at the Z Pole EPJC 71 (2011) 1557 1102.4748
56 ALEPH Collaboration Study of the fragmentation of b quarks into B mesons at the Z peak PLB 512 (2001) 30 hep-ex/0106051
57 S. Dulat et al. New parton distribution functions from a global analysis of quantum chromodynamics PRD 93 (2016), no. 3, 033006 1506.07443
58 L. A. Harland-Lang, A. D. Martin, P. Motylinski, and R. S. Thorne Parton distributions in the LHC era: MMHT 2014 PDFs EPJC 75 (2015), no. 5, 204 1412.3989
59 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 CMS-PAS-TOP-16-021 CMS-PAS-TOP-16-021
60 S. Mrenna and P. Skands Automated Parton-Shower Variations in Pythia 8 PRD 94 (2016), no. 7, 074005 1605.08352
61 M. Czakon, D. Heymes, and A. Mitov High-precision differential predictions for top-quark pairs at the LHC PRL 116 (2016) 082003 1511.00549
62 M. Czakon et al. Top-pair production at the LHC through NNLO QCD and NLO EW JHEP 10 (2017) 186 1705.04105
63 S. Catani et al. Top-quark pair production at the LHC: Fully differential QCD predictions at NNLO JHEP 07 (2019) 100 1906.06535
64 CMS Collaboration Measurement of differential cross sections for top quark pair production using the lepton+jets final state in proton-proton collisions at 13 TeV PRD 95 (2017) 092001 CMS-TOP-16-008
1610.04191
65 CMS Collaboration Measurement of normalized differential $ \mathrm{t}\overline{\mathrm{t}} $ cross sections in the dilepton channel from pp collisions at $ \sqrt{s}= $ 13 TeV JHEP 04 (2018) 060 CMS-TOP-16-007
1708.07638
66 J. R. Christiansen and P. Z. Skands String formation beyond leading colour JHEP 08 (2015) 003 1505.01681
67 S. Argyropoulos and T. Sjostrand Effects of color reconnection on $ \mathrm{t\bar{t}} $ final states at the LHC JHEP 11 (2014) 043 1407.6653
68 W. Verkerke and D. P. Kirkby The RooFit toolkit for data modeling eConf C0303241 (2003) MOLT007 physics/0306116
69 F. James and M. Roos Minuit: A System for Function Minimization and Analysis of the Parameter Errors and Correlations CPC 10 (1975) 343
70 P. K. Suetin Chebyshev polynomials Springer-Verlag
Compact Muon Solenoid
LHC, CERN