Constraining anomalous Higgs boson couplingin H+\gamma production

Constraining anomalous Higgs boson coupling
in + production

Liaoshan Shi School of Physics, Sun Yat-sen University, Guangzhou 510275, China
   Zhijun Liang    Bo Liu Institute of High Energy Physics, Chinese Academy of Science, Beijing 100049, China
   Zhenhui He School of Physics, Sun Yat-sen University, Guangzhou 510275, China

Higgs boson production in association with a photon (+) offers a promising channel to test the Higgs boson to photon coupling at various energy scales. Its potential sensitivity to anomalous couplings of the Higgs boson has not been explored with proton-proton collision data. In this paper, we reinterpret the latest ATLAS + resonance search results within the Standard Model effective field theory (EFT) framework, using 36.1 fb of proton-proton collision data recorded with the ATLAS detector at TeV. Constraints on Wilson coefficients of dimension-six EFT operators related to the Higgs boson to photon coupling are provided for the first time in + final state at the LHC.

preprint: Submitted to Chinese Physics Cthanks: Supported by the National Natural Science Foundation of China (11875278), Beijing Municipal Science & Technology Commission (Z181100004218003), the National Key Research and Development Program of China (2018YFA0404001), and the Hundred Talent Program of the Chinese Academy of Sciences (Y6291150K2).

I Introduction

After the discovery of the Higgs boson HIGG-2012-27 (); CMS-HIG-12-028 (), measurements of the Higgs boson couplings to the other fundamental particles become crucial tests of the nature of the Higgs boson. In the Standard Model (SM), coupling of the Higgs boson to photon is forbidden at the tree level, and is induced by heavy particle loops in, e.g., and processes. The Higgs-photon coupling has been extensively studied in various Higgs boson decay channels including and with LHC data recorded by the ATLAS and CMS experiments HIGG-2015-02 (); Aaboud:2018xdt (); Aaboud:2017vzb (); ATL-PHYS-PUB-2017-018 (); CMS-HIG-14-018 (); CMS-HIG-15-002 (); CMS-HIG-17-011 (); Englert:2015hrx (); Ellis:2018gqa ().

Apart from the Higgs boson decay channels involving photons, Higgs boson production in association with a photon can also be used to measure the Higgs-photon coupling. The production cross section is predicted to be very small in the SM, but anomalous couplings introduced by models beyond the SM (BSM) can have significant effects on it. The process was considered as a promising and clean channel at LEP Djouadi:1996ws (); PhysRevD.52.3919 (), and has been used by the DELPHI collaboration to search for anomalous couplings of the Higgs boson to vector bosons Abreu:1999vt (). At the LHC, potential sensitivity of the process to anomalous Higgs-photon couplings has been discussed in Ref. Khanpour:2017inb (). It is predicted that some Wilson coefficients of dimension-six operators related to Higgs-photon couplings can be probed down to with 300 fb collision data at 14 TeV. There is no particular analysis measuring the anomalous Higgs-photon couplings via this channel using the LHC data. ATLAS and CMS collaborations have reported the results on heavy resonance searches in 13 TeV collision data Aaboud:2018fgi (); CMS-PAS-EXO-17-019 (). Besides the resonance models, their results are also sensitive to non-resonant production and to the anomalous coupling between Higgs boson and photon. However, these results have not been interpreted as limits on the anomalous Higgs-photon coupling.

In this paper, the latest resonance search results from the ATLAS collaboration Aaboud:2018fgi () are reinterpreted within the SM effective field theory (EFT) framework, and are presented as constraints on Wilson coefficients of dimension-six EFT operators. The study is based on a collision dataset of 36.1 fb at TeV.

This paper is organized as follows. Section II gives a short overview of the EFT framework, and a brief description of the signal Monte Carlo generation for the reinterpretation. Section III describes the analysis strategy. Section IV presents the constraints on Wilson coefficients of dimension-six EFT operators that are obtained in the channel, and compares them with existing results from other measurements. Our conclusions are summarized in Section V.

Ii Effective field theory

In the SM effective field theory approach, effects of BSM interactions are parametrized using higher-dimension operators in addition to the SM Lagrangian. Leading contributions at collider energies are expected to originate from dimension-six operators. A general effective Lagrangian with dimension-six operators takes the form


with Wilson coefficients describing the strengths of the BSM interactions.

We focus on a set of dimension-six operators known as the strongly-interacting light Higgs (SILH) Lagrangian Giudice:2007fh (). It is written as

where , and are the gauge field strength tensors, and is the Higgs doublet. Among all the Wilson coefficients in the SILH Lagrangian, , and are related to the anomalous Higgs-photon coupling through a direct or vertex. With the presence of these BSM vertices, additional tree level diagrams, in particular an -channel diagram via a virtual photon or Z boson as the mediator, can contribute to the process and lead to a large relative change in its production cross section. Therefore, the process is a sensitive probe to explore the anomalous Higgs-photon coupling Khanpour:2017inb ().

A public implementation of the SILH Lagrangian is available in a general Higgs Effective Lagrangian (HEL) Alloul:2013naa (); Contino:2013kra (). The HEL model is implemented in FeynRules Alloul:2013bka (), comprising 39 dimension-six operators and their corresponding Wilson coefficients. Its Universal FeynRules Output Degrande:2011ua () has been interfaced to the MadGraph5_aMC@NLO Alwall:2014hca () event generator. In this work, the HEL model is used with all the other Wilson coefficients fixed at 0 except , and . The production cross section is computed at different values of , and , using MadGraph5_aMC@NLO v2.6.2 with NNPDF2.3 Ball:2012cx () parton distribution functions. We then parametrize the signal cross section as functions of the Wilson coefficients according to the computation. Figure 1 presents a two-dimensional parametrization of the signal cross section parametrized as functions of every two of the three Wilson coefficients with the third coefficient fixed at 0. Monte Carlo event samples are also generated with the same configurations.

Figure 1: Two-dimensional parametrization of the signal cross section of the process with different values of , and . Besides the two parameters indicated in each plot, the third parameter is fixed at 0.

Iii Analysis strategy

The ATLAS resonance search Aaboud:2018fgi () is carried out to search for heavy resonances decaying to a SM Higgs boson and a photon, using the decay of the Higgs boson. In its signal region, both the selected photon and the Higgs boson are highly boosted (with large momenta). The search is performed by looking for a bump on a smooth background of the invariant mass spectrum . As reported by the ATLAS paper, the mass spectrum observed is consistent with the background-only hypothesis and no evidence of new resonances is found.

This highly boosted signature is of particular interest for probing the anomalous Higgs-photon coupling, as the BSM signal contribution may show longer tails extending up to TeV scale in the and photon distributions, while the SM expectation drops more steeply Khanpour:2017inb (). Instead of performing a bump hunt on the spectrum as in the original ATLAS paper, we perform a counting experiment with the published spectrum to constrain anomalous coupling of the Higgs boson. According to the ATLAS paper Aaboud:2018fgi (), 138 events have been observed in signal region , consistent with the expected number of background events . We reinterpret the ATLAS data as follows.

The expected number of events in signal region can be expressed as , with and being the expected number of signal and background events, respectively. To constrain the Wilson coefficients parameters, we construct a likelihood function assuming that the number of observed events follows a Poisson distribution with an expectation value :


Here is treated as a nuisance parameter. It is constrained by a Gaussian term with a mean value and a standard deviation . Both and are obtained from the background fit results in the ATLAS paper Aaboud:2018fgi (). The expected number of signal events depends on the Wilson coefficients . It can be further expressed as:


The integrated luminosity of the ATLAS data sample is 36.1 fb. A Standard Model branching ratio for a 125 GeV Higgs boson deFlorian:2016spz () is used. The signal efficiency accounts for event loss due to detector effects, reconstruction and selection efficiencies in the ATLAS analysis. It is determined by applying the efficiency table published in the ATLAS paper Aaboud:2018fgi () to the simulated spectra in the signal Monte Carlo samples and evaluate the overall efficiency. The production cross section is computed in terms of Wilson coefficients , and , as described in Section II. Wilson coefficients , and are treated as the parameters of interest (POIs).

Constraints on the Wilson coefficients are obtained by evaluating the profiled likelihood ratio assuming the asymptotic approximation Cowan:2010js ():


Here the numerator is the conditional maximum-likelihood function, with being the value of nuisance parameter that maximize the likelihood function for a given set of values of the Wilson coefficients . The denominator is the unconditional maximum-likelihood function with and being the maximum-likelihood estimates of and , respectively.

Iv Results and discussions

A one-dimensional likelihood scan is performed to obtain constraints on each of the three Wilson coefficients in the EFT framework with the other two fixed at 0. Constraints on , and are shown in Figure 2. The 68% and 95% confidence intervals are shown in Table 1.

Figure 2: One-dimensional likelihood scan of the Wilson coefficients (left), (middle) and (right) in the EFT framework with all the other coefficients fixed at 0. The 95% (68%) confidence interval is indicated by the red (green) line. Constraints are obtained from data in mass range of 800  3.2 .
Parameter 68% C.L. 95% C.L.
[-0.061, 0.064] [-0.087, 0.090]
[-0.167, 0.161] [-0.236, 0.231]
[-0.162, 0.167] [-0.230, 0.236]
Table 1: 68% and 95% confidence intervals on the Wilson coefficients , and in the EFT framework.

Two-dimensional likelihood scans are also performed and the confidence regions are shown in Figure 3. Besides the two Wilson coefficients indicated in the plot, the remaining one is fixed at 0 during the scan.

Figure 3: Two-dimensional likelihood scan of the Wilson coefficients in the EFT framework. Besides the two parameters indicated in each plot, the third parameter is fixed at 0. The 95% (68%) confidence region is indicated by the red (green) contour. Constraints are obtained from data in mass range of 800  3.2 . The SM expectation at (0, 0) is also shown.

We compare the channel results with those obtained in the combined and channels based on the same collision dataset collected by the ATLAS experiment. The 68% C.L. intervals from the combined channels ATL-PHYS-PUB-2017-018 () are:

The limits on is much stringent, but the limits on and are in the same order of magnitude compared to the channel results. These results demonstrate excellent sensitivity of the production process to some of the Wilson coefficients in the EFT framework. Combination of the channel with other channels will further improve the sensitivity on the Higgs boson anomalous couplings.

Limits can be further improved by considering shape information from differential distributions instead of doing a simple counting experiment. In the channel, the 95% C.L. observed limit of has been improved to after including differential distributions Aaboud:2018xdt (), which is four times better than the limit achieved from the channel in our study. In the channel, improvements to the sensitivity are also anticipated by including additional information from the and photon distributions, but we consider a shape analysis beyond the scope of this paper.

V Conclusions

We present an interpretation of the recent ATLAS resonance search results with 36.1 fb collision data at TeV on the search for Higgs boson anomalous coupling in the final state. We provide constraints on Wilson coefficients of dimension-six EFT operators for the first time in final state with collision data. These results demonstrate excellent physics potential of the production process. With differential cross sections measured for the process in the future, the constraints can be further improved.


  • [1] ATLAS Collaboration. Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC. Phys. Lett. B, 716:1, 2012.
  • [2] CMS Collaboration. Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC. Phys. Lett. B, 716:30, 2012.
  • [3] ATLAS Collaboration. Constraints on non-Standard Model Higgs boson interactions in an effective Lagrangian using differential cross sections measured in the decay channel at with the ATLAS detector. Phys. Lett. B, 753:69, 2016.
  • [4] ATLAS Collaboration. Measurements of Higgs boson properties in the diphoton decay channel with 36 fb of collision data at TeV with the ATLAS detector. Phys. Rev. D, 98:052005, 2018.
  • [5] ATLAS Collaboration. Measurement of the Higgs boson coupling properties in the decay channel at = 13 TeV with the ATLAS detector. JHEP, 03:095, 2018.
  • [6] ATLAS Collaboration. Constraints on an effective Lagrangian from the combined and channels using 36.1 fb of 13 TeV collision data collected with the ATLAS detector. (ATL-PHYS-PUB-2017-018), Nov 2017.
  • [7] CMS Collaboration. Constraints on the spin-parity and anomalous HVV couplings of the Higgs boson in proton collisions at and . Phys. Rev. D, 92:012004, 2015.
  • [8] 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 collision data at and . JHEP, 08:045, 2016.
  • [9] CMS Collaboration. Constraints on anomalous Higgs boson couplings using production and decay information in the four-lepton final state. Phys. Lett. B, 775:1, 2017.
  • [10] Christoph Englert, Roman Kogler, Holger Schulz, and Michael Spannowsky. Higgs coupling measurements at the LHC. Eur. Phys. J. C, 76(7):393, 2016.
  • [11] John Ellis, Christopher W. Murphy, Verónica Sanz, and Tevong You. Updated Global SMEFT Fit to Higgs, Diboson and Electroweak Data. JHEP, 06:146, 2018.
  • [12] A. Djouadi, V. Driesen, W. Hollik, and J. Rosiek. Associated production of Higgs bosons and a photon in high-energy e+ e- collisions. Nucl. Phys. B, 491:68–102, 1997.
  • [13] Ali Abbasabadi, David Bowser-Chao, Duane A. Dicus, and Wayne W. Repko. Higgs-boson–photon associated production at eē colliders. Phys. Rev. D, 52:3919–3928, Oct 1995.
  • [14] DELPHI Collaboration. Search for the Higgs boson in events with isolated photons at LEP-2. Phys. Lett. B, 458:431–446, 1999.
  • [15] Hamzeh Khanpour, Sara Khatibi, and Mojtaba Mohammadi Najafabadi. Probing Higgs boson couplings in H+ production at the LHC. Phys. Lett. B, 773:462–469, 2017.
  • [16] ATLAS Collaboration. Search for heavy resonances decaying to a photon and a hadronically decaying boson in collisions at with the ATLAS detector. Phys. Rev. D, 98(3):032015, 2018.
  • [17] CMS Collaboration. Search for narrow H-gamma resonances in proton-proton collisions at sqrt(s) = 13 TeV. (CMS-PAS-EXO-17-019), 2018.
  • [18] G. F. Giudice, C. Grojean, A. Pomarol, and R. Rattazzi. The Strongly-Interacting Light Higgs. JHEP, 06:045, 2007.
  • [19] Adam Alloul, Benjamin Fuks, and Verónica Sanz. Phenomenology of the Higgs Effective Lagrangian via FEYNRULES. JHEP, 04:110, 2014.
  • [20] Roberto Contino, Margherita Ghezzi, Christophe Grojean, Margarete Muhlleitner, and Michael Spira. Effective Lagrangian for a light Higgs-like scalar. JHEP, 07:035, 2013.
  • [21] Adam Alloul, Neil D. Christensen, Céline Degrande, Claude Duhr, and Benjamin Fuks. FeynRules 2.0 - A complete toolbox for tree-level phenomenology. Comput. Phys. Commun., 185:2250–2300, 2014.
  • [22] Celine Degrande, Claude Duhr, Benjamin Fuks, David Grellscheid, Olivier Mattelaer, and Thomas Reiter. UFO - The Universal FeynRules Output. Comput. Phys. Commun., 183:1201–1214, 2012.
  • [23] J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer, H. S. Shao, T. Stelzer, P. Torrielli, and M. Zaro. The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations. JHEP, 07:079, 2014.
  • [24] Richard D. Ball et al. Parton distributions with LHC data. Nucl. Phys. B, 867:244–289, 2013.
  • [25] D. de Florian et al. Handbook of LHC Higgs Cross Sections: 4. Deciphering the Nature of the Higgs Sector. 2016.
  • [26] Glen Cowan, Kyle Cranmer, Eilam Gross, and Ofer Vitells. Asymptotic formulae for likelihood-based tests of new physics. Eur. Phys. J. C, 71:1554, 2011. [Erratum: Eur. Phys. J. C,73: 2501(2013)].
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description