CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-TOP-17-003
Search for the flavor-changing interactions of the top quark with the Higgs boson in $\mathrm{ H \to b\bar{b} }$ channel at $\sqrt{s}= $ 13 TeV
Abstract: A search for flavor-changing neutral currents (FCNC) in associated production of a top quark with a Higgs boson is presented. Two channels are considered: top quark pair production, with FCNC decay of the top quark or antiquark, and single top production ($\mathrm{ pp \to tH }$). This is the first analysis to probe top-Higgs FCNC couplings in the single top production mode. A final state with one isolated lepton and at least three reconstructed jets, among which at least two are identified as b quark jets, is studied. The data sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$ recorded by the CMS experiment at the LHC in proton-proton collisions at $\sqrt{s}= $ 13 TeV in 2016. No significant deviation is observed from predicted background. Observed (expected) upper limits at 95% confidence level are set on the branching ratios of top quark decays, $B(\mathrm{ t \rightarrow u H }) < $ 0.47% (0.34%) and $B(\mathrm{ t \rightarrow c H }) < $ 0.47% (0.44%), assuming a single non-vanishing FCNC coupling.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
Representative Feynman diagrams for FCNC top-Higgs processes: associated production of the top quark with the Higgs boson (left) and FCNC decay of the top quark in $ {\mathrm{t} {}\mathrm{\bar{t}}} $ (right).

png
Figure 1-a:
Representative Feynman diagrams for FCNC top-Higgs processes: associated production of the top quark with the Higgs boson (left) and FCNC decay of the top quark in $ {\mathrm{t} {}\mathrm{\bar{t}}} $ (right).

png
Figure 1-b:
Representative Feynman diagrams for FCNC top-Higgs processes: associated production of the top quark with the Higgs boson (left) and FCNC decay of the top quark in $ {\mathrm{t} {}\mathrm{\bar{t}}} $ (right).

png pdf
Figure 2:
Comparison between data and MC for the most important BDT input variables: lepton charge (top left), CSVv2 discriminant value for one of the reconstructed b jets assigned to Higgs decay (top right), reconstructed mass of the Higgs boson (bottom left) and the maximum BDT score value from the b jet assignment procedure (bottom right). Each of the signal distributions is normalized to the total number of events in the predicted background to ease the comparison of the shapes. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 2-a:
Comparison between data and MC for the most important BDT input variables: lepton charge (top left), CSVv2 discriminant value for one of the reconstructed b jets assigned to Higgs decay (top right), reconstructed mass of the Higgs boson (bottom left) and the maximum BDT score value from the b jet assignment procedure (bottom right). Each of the signal distributions is normalized to the total number of events in the predicted background to ease the comparison of the shapes. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 2-b:
Comparison between data and MC for the most important BDT input variables: lepton charge (top left), CSVv2 discriminant value for one of the reconstructed b jets assigned to Higgs decay (top right), reconstructed mass of the Higgs boson (bottom left) and the maximum BDT score value from the b jet assignment procedure (bottom right). Each of the signal distributions is normalized to the total number of events in the predicted background to ease the comparison of the shapes. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 2-c:
Comparison between data and MC for the most important BDT input variables: lepton charge (top left), CSVv2 discriminant value for one of the reconstructed b jets assigned to Higgs decay (top right), reconstructed mass of the Higgs boson (bottom left) and the maximum BDT score value from the b jet assignment procedure (bottom right). Each of the signal distributions is normalized to the total number of events in the predicted background to ease the comparison of the shapes. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 2-d:
Comparison between data and MC for the most important BDT input variables: lepton charge (top left), CSVv2 discriminant value for one of the reconstructed b jets assigned to Higgs decay (top right), reconstructed mass of the Higgs boson (bottom left) and the maximum BDT score value from the b jet assignment procedure (bottom right). Each of the signal distributions is normalized to the total number of events in the predicted background to ease the comparison of the shapes. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 3:
BDT output score distributions for different jet categories for Hut training after the fit to data. All background processes are constrained to the SM expectation in the fit. To compare the shapes, each of the signal distributions is normalized to the total number of events in the predicted background. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 3-a:
BDT output score distributions for different jet categories for Hut training after the fit to data. All background processes are constrained to the SM expectation in the fit. To compare the shapes, each of the signal distributions is normalized to the total number of events in the predicted background. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 3-b:
BDT output score distributions for different jet categories for Hut training after the fit to data. All background processes are constrained to the SM expectation in the fit. To compare the shapes, each of the signal distributions is normalized to the total number of events in the predicted background. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 3-c:
BDT output score distributions for different jet categories for Hut training after the fit to data. All background processes are constrained to the SM expectation in the fit. To compare the shapes, each of the signal distributions is normalized to the total number of events in the predicted background. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 3-d:
BDT output score distributions for different jet categories for Hut training after the fit to data. All background processes are constrained to the SM expectation in the fit. To compare the shapes, each of the signal distributions is normalized to the total number of events in the predicted background. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 4:
BDT output score distributions for different jet categories for Hct training after the fit to data. All background processes are constrained to the SM expectation in the fit. Each of the signal distributions is normalized to the total number of events in the predicted background to ease the comparison of the shapes. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 4-a:
BDT output score distributions for different jet categories for Hct training after the fit to data. All background processes are constrained to the SM expectation in the fit. Each of the signal distributions is normalized to the total number of events in the predicted background to ease the comparison of the shapes. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 4-b:
BDT output score distributions for different jet categories for Hct training after the fit to data. All background processes are constrained to the SM expectation in the fit. Each of the signal distributions is normalized to the total number of events in the predicted background to ease the comparison of the shapes. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 4-c:
BDT output score distributions for different jet categories for Hct training after the fit to data. All background processes are constrained to the SM expectation in the fit. Each of the signal distributions is normalized to the total number of events in the predicted background to ease the comparison of the shapes. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 4-d:
BDT output score distributions for different jet categories for Hct training after the fit to data. All background processes are constrained to the SM expectation in the fit. Each of the signal distributions is normalized to the total number of events in the predicted background to ease the comparison of the shapes. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 4-e:
BDT output score distributions for different jet categories for Hct training after the fit to data. All background processes are constrained to the SM expectation in the fit. Each of the signal distributions is normalized to the total number of events in the predicted background to ease the comparison of the shapes. The shaded area corresponds to a total uncertainty on the predicted background.

png pdf
Figure 5:
Excluded signal cross section at 95% CL per event category for Hut (left) and Hct (right).

png pdf
Figure 5-a:
Excluded signal cross section at 95% CL per event category for Hut (left) and Hct (right).

png pdf
Figure 5-b:
Excluded signal cross section at 95% CL per event category for Hut (left) and Hct (right).

png pdf
Figure 6:
Upper limits on $B(\mathrm{ t \rightarrow u H })$ and $B(\mathrm{ t \rightarrow c H })$ at 95% CL.

png pdf
Figure 7:
The best fit signal strength for Hut (left) and Hct (right). The signal strength is restricted to positive values in the fit.

png pdf
Figure 7-a:
The best fit signal strength for Hut (left) and Hct (right). The signal strength is restricted to positive values in the fit.

png pdf
Figure 7-b:
The best fit signal strength for Hut (left) and Hct (right). The signal strength is restricted to positive values in the fit.
Tables

png pdf
Table 1:
Number of events in each category together with its total relative uncertainty as obtained from the fit to data for Hut.

png pdf
Table 2:
Number of events in each category together with its total relative uncertainty as obtained from the fit to data for Hct.
Summary
A search for flavor-changing neutral currents in events with the top quark and the Higgs boson corresponding to a data sample of 35.9 fb$^{-1}$ collected in proton-proton collisions at $\sqrt{s}= $ 13 TeV is presented. This is the first search to probe top-Higgs FCNC couplings in both associated production of the top quark with the Higgs boson and in top quark decays. Observed (expected) upper limits at 95% confidence level are set on the branching ratios of top quark decays, $B(\mathrm{ t \rightarrow u H }) < $ 0.47% (0.34%) and $B(\mathrm{ t \rightarrow c H }) < $ 0.47% (0.44%). Inclusion of the associated production of a single top quark with a Higgs boson in this study provides a $\simeq$20% relative improvement in the upper limit on $B(\mathrm{ t \rightarrow uH })$.
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) 1 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) 30 CMS-HIG-12-028
1207.7235
3 CMS Collaboration Observation of a new boson with mass near 125 gev in pp collisions at $ \sqrt{s} = $ 7 and 8 TeV JHEP 06 (2013) 081 CMS-HIG-12-036
1303.4571
4 ATLAS and CMS Collaborations Measurements of the higgs boson production and decay rates and constraints on its couplings from a combined atlas and cms analysis of the lhc pp collision data at $ \sqrt{s} = $ 7 and 8 TeV JHEP 08 (2016) 045 1606.02266
5 S. L. Glashow, J. Iliopoulos and L. Maiani Weak interactions with lepton-hadron symmetry PRD 2 (1970) 1285
6 G. Eilam, J. L. Hewett, and A. Soni Rare decays of the top quark in the standard and two higgs doublet models PRD 44 (1991) 1473
7 B. Mele, S. Petrarca, and A. Soddu A new evaluation of the $ t \rightarrow ch $ decay width in the standard model PLB 435 (1998) 401 9805498
8 J. A. Aguilar-Saavedra Top flavor-changing neutral interactions: theoretical expectations and experimental detection Acta Phys. Polon. B 35 (2004) 2695 0409342
9 ATLAS Collaboration Search for flavour-changing neutral current top quark decays $ t \rightarrow hq $ in pp collisions at $ \sqrt{s} = $ 8 TeV with the atlas detector JHEP 12 (2015) 061
10 ATLAS Collaboration Search for top quark decays $ t \rightarrow qH $, with $ H \rightarrow \gamma\gamma $, in $ \sqrt{s} = $ 13 TeV pp collisions using the ATLAS detector Submitted to JINST 1707.01404
11 CMS Collaboration Search for top quark decays via higgs-boson-mediated flavor-changing neutral currents in pp collisions at $ \sqrt{s} = $ 8 TeV JHEP 02 (2017) 079 CMS-TOP-13-017
1610.04857
12 CMS Collaboration Cms luminosity measurement for the 2016 data taking period CMS-PAS-LUM-17-001 CMS-PAS-LUM-17-001
13 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
14 M. L. Mangano, M. Moretti, F. Piccinini, and M. Treccani Matching matrix elements and shower evolution for top-quark production in hadronic collisions JHEP 01 (2007) 013 0611129
15 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
16 A. A. et al. FeynRules 2.0 - A complete toolbox for tree-level phenomenology CPC 185 (2014) 2250--2300 1310.1921
17 C. D. et al. UFO - The Universal FeynRules Output CPC 183 (2012) 1201--1214 1108.2040
18 P. Nason A new method for combining nlo qcd with shower monte carlo algorithms JHEP 11 (2004) 040 0409146
19 S. Frixione, P. Nason, and C. Oleari Matching nlo qcd computations with parton shower simulations: the powheg method JHEP 11 (2007) 070 0709.2092
20 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
21 S. Alioli, S.-O. Moch, and P. Uwer Hadronic top-quark pair-production with one jet and parton showering JHEP 01 (2012) 137 1110.5251
22 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
23 M. Botje et al. The pdf4lhc working group interim recommendations 1101.0538
24 A.D Martin et al. Uncertainties on alphas in global pdf analyses and implications for predicted hadronic cross sections EPJC 64 (2009) 653--680 0905.3531
25 J. Gao et al. The ct10 nnlo global analysis of qcd PRD 89 (2014) 033009 1302.6246
26 R.D. Ball et al. Parton distributions with lhc data NPB 867 (2013) 244--289 1207.1303
27 M. Aliev et al. Hathor - hadronic top and heavy quarks cross section calculator CPC 182 (2011) 1034--1046
28 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
29 N. Kidonakis Two-loop soft anomalous dimensions for single top quark associated production with a w- or h- Phys.Rev.D 82 (2010) 054018 1005.4451
30 T. Sjostrand, S. Mrenna, and P. Skands Pythia 6.4 physics and manual JHEP 05 (2006) 026
31 T. Sjostrand, S. Mrenna, and P. Skands A brief introduction to pythia 8.1 Comp.Phys.Commun. 178 (2008) 852
32 CMS Collaboration Event generator tunes obtained from underlying event and multiparton scattering measurements Eur.Phys.J.C 76 (2016) 155 CMS-GEN-14-001
1512.00815
33 CMS Collaboration Investigations of the impact of the parton shower tuning in Pythia 8 in the modelling of $ \mathrm{t\bar{t}} $ at $ \sqrt{s} = $ 8 and 13 TeV CMS-PAS-TOP-16-021 CMS-PAS-TOP-16-021
34 J. Allison et al. Geant4 developments and applications IEEE Trans. Nucl. Sci. 53 (2006) 270
35 A. M. Sirunyan et al. Particle-flow reconstruction and global event description with the CMS detector Submitted to JINST 1706.04965
36 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ k_t $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
37 M. Cacciari, G. P. Salam, and G. Soyez FastJet user manual EPJC 72 (2012) 1896 1111.6097
38 CMS Collaboration Determination of jet energy calibration and transverse momentum resolution in CMS JINST 6 (2011) P11002 1107.4277
39 CMS Collaboration Identification of b quark jets at the CMS Experiment in the LHC Run 2 CMS-PAS-BTV-15-001 CMS-PAS-BTV-15-001
40 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
41 L. Moneta et al. The RooStats Project in 13$^{\text{th}}$ International Workshop on Advanced Computing and Analysis Techniques in Physics Research (ACAT2010) SISSA, 2010 1009.1003
42 A. Read Presentation of search results: the CLs technique J.Phys.G 28 (2002) 2693
43 T. Junk Confidence level computation for combining searches with small statistics NIMA 434 (1999) 435 9902006
44 G. Cowan, K. Cranmer, E. Gross, and O. Vitells Asymptotic formulae for likelihood-based tests of new physics Eur.Phys.J. C 71 (2011) 1554 1007.1727
Compact Muon Solenoid
LHC, CERN