Probing the top-Higgs coupling through the secondary lepton distributions in the associated production of the top-quark pair and Higgs boson at the LHC

We complement the analysis of the anomalous top-Higgs coupling effects on the secondary lepton distributions in the associated production of the top-quark pair and Higgs boson in proton–proton collisions at the LHC of the former work by one of the present authors by taking into account the quark–antiquark production mechanism. We also present simple arguments which explain why the effects of the scalar and pseudoscalar anomalous couplings on the unpolarized cross section of the process are completely insensitive to the sign of either of them.


Introduction
Determination of the coupling of the recently discovered Higgs boson [1,2] to the top quark currently belongs to one of the most challenging tasks of high energy experimental physics. Measurement of the associated production of the top-quark pair and Higgs boson in the clean experimental environment of e + e − collisions was considered in this context already more than two decades ago [3,4], but different projects of the high energy e + e − collider [5][6][7][8][9][10][11][12], despite some of them being more or less intensively discussed for years, are still at a rather early stage of TDR. However, if the LHC performance in next runs is as excellent as it was in run 1 we may expect that the process the search for which, based on run 1 data, were already reported by both the CMS [13] and the ATLAS [14,15] collaborations, will be measured quite precisely. This is why in the past few years the associated production of the top-quark a e-mail: karol.kolodziej@us.edu.pl b e-mail: aleksandra.slapik@gmail.com pair and Higgs boson has invoked quite some interest also from the theoretical side; see, e.g., [16][17][18][19][20][21][22][23][24][25]. It was shown in Ref. [19] that the distributions in rapidity and angles of the secondary lepton that can be produced in the decay of thet-quark of process (1) are quite sensitive to modifications of the top-Higgs coupling. Actually, only the gluon fusion mechanism of tt H production, which is dominant at the LHC energies, and one specific decay channel: t → bW + → bud,t →bW − →bμ −ν μ , and h → bb, were taken into account in Ref. [19], i.e., the hard parton scattering process gg → budbμ −ν μ bb (2) was considered. There are 67300 Feynman diagrams of process (2) already in the leading order (LO) of the standard model (SM) in the unitary gauge, if the Cabibbo-Kobayashi-Maskawa mixing and masses smaller than the b-quark mass m b are neglected. At the same time there are only 32 Feynman diagrams which contribute to the signal cross section of tt H production, two of which are shown in Fig. 1. The remaining 30 signal diagrams are obtained from those depicted by attaching the Higgs boson line of the Hbb vertex to the other t-ort-quark line, or interchanging the b andb quarks in Fig. 1a, b and interchanging the two gluons in Fig. 1b. The diagrams with the Higgs boson line of Hbb vertex attached to either the b-orb-quark line are not counted here, as their contribution to the tt H production signal is suppressed by the mass ratio m b /m t . The effects caused by modifications of the scalar and pseudoscalar couplings of the Higgs boson to top quark were clearly visible in the tt H production signal cross section, but they were to large degree obscured by the irreducible off-resonance background.
In the present work, we complement the analysis of the influence of the anomalous Higgs boson coupling to top quark on the secondary lepton distributions in the process of associated production of the top-quark pair and Higgs Fig. 1 Feynman diagrams of tt H production in process (2). Blobs indicate the top-Higgs coupling Feynman diagrams of tt H production in the process uū → budbμ −ν μ bb. Blobs indicate the top-Higgs coupling boson in proton-proton collisions at the LHC of Ref. [19] by taking into account the quark-antiquark annihilation hardscattering processes with the same final state as that of process (2): with q = u, d. To be more specific, we take into account uū-,ūu-, dd-anddd-scattering processes. Under the same assumptions as those made above for process (2), there are 78068 Feynman diagrams in the LO of SM for each of the qqscattering processes considered. However, only 24 of them contribute to the signal of the tt H production. Examples of the signal diagrams of the process of uū-scattering to the final state of process (3) are shown in Fig. 2. The other signal diagrams can be obtained by attaching the Higgs boson line of the Hbb vertex to the other t-ort-quark line or interchanging the b-andb-quark lines in the diagrams of Fig. 2. Let us note that another 24 diagrams which contain the Feynman propagators of the t-,t-quark and the Higgs boson at a time can be obtained from the signal diagrams just described by the exchange of the u-quark lines between the initial and final state. However, they are not treated as the signal diagrams here, because they contain the gluon, Z 0 or photon propagator in the t-or u-channel and their contribution to the signal cross section is negligible anyway, which has been checked by direct computation.
The rest of the article is organized in the following way. The possible effect of the anomalous top-Higgs coupling on the unpolarized cross section of the process of tt H production at the LHC are analyzed in Sect. 2, our results are presented in Sect. 3 and, finally, some concluding remarks are contained in Sect. 4.

Effects of the anomalous top-Higgs coupling
The most general top-Higgs coupling is given by the following Lagrangian [26]: where GeV, is the top-Higgs Yukawa coupling and the real couplings f and f describe, respectively, the scalar and pseudoscalar departures from the purely scalar top-Higgs Yukawa coupling of SM, which is reproduced for f = 1 and f = 0. The allowed regions of the ( f, f ) plane, according to the analysis of Ref. [20] performed at the 68 and 95 % confidence level, are plotted in Fig. 1 of Ref. [21]. They are derived from the constraints on the Hgg and H γ γ couplings from the Higgs boson production and its decay into γ γ , which among others involve assumptions on the Higgs boson couplings to other fermions and bosons, and hence are model dependent. Therefore, we will not stick to them in the next section, where we will illustrate the effects of f on the process of associated production of the top-quark pair and Higgs boson from which the direct constraints on f and f can be derived.
Let us try to predict the possible effect of the top-Higgs coupling given by Eq. (4) on the unpolarized cross section of the process uū → budbμ −ν μ bb. To this end, let us consider the amplitudes of two dominant diagrams of the tt H production: M (1) a of the diagram depicted in Fig. 2a and M (2) a of the diagram obtained from that of Fig. 2a by attaching the Higgs boson line to thet-quark. They have the following form: where h is a scalar representing a product of the Higgs boson propagator carrying the four-momentum q with the Hbb vertex, u (v) is the Dirac spinor representing the off-shell t-quark (t-quark) of the four-momentum p t ( pt ) that decays into the b-quark (b-quark) and off-shell W + (W − ) boson, ε is a polarization four-vector representing the gluon propagator contracted with the uūg vertex, g s is the strong coupling mass parameter that replaces the mass m t in the top-quark propagator in order to regularize the pole arising if its denominator approaches zero. After some simple algebra Eqs. (5) and (6) can be written in the following form: Now, let us note that, as in the process of tt H production in e + e − collisions that was considered in Refs. [16][17][18], the dominant contribution to the cross section comes from the phase space region, where both the t-quark and thet-quark are close to their mass shells and hence the off-shell spinors u and v should satisfy the following approximate equations: Using Eqs. (9) in (7) and (10) in (8), and neglecting terms ∼ t in the numerators, we get the following approximate expressions for the amplitudes: and for a sum of the two: In order to calculate the sum over polarizations of the squared module of the matrix element pol. |M a | 2 , we take into account the approximate completeness relations for the spinors u and v: pol.
and note that the off-shell polarization four-vectors ε are real, as they are defined in the following way: where the helicity spinors v( p 1 , λ 1 ) and u( p 2 , λ 2 ) of, respectively, theū-and u-quark in initial state, which are calculated according to Eqs. (5) and (6) of Ref. [27], are real if the momenta p 1 and p 2 are antiparallel. Thus pol.
More simplified analytic form of Eq. (16) is irrelevant, as the calculation of the cross section will be performed numerically anyway, but let us note that only the terms on the r.h.s. of Eq. (16) that contain γ 5 may be proportional to the product f f . However, if we use the relation ε / p t / ε / = −ε 2 p t / + 2( p t · ε)ε / in the second and third term and the relation q / p t / q / = −q 2 p t / + 2( p t · q)q / in the second term, and then use the relation in the last term on the r.h.s. of Eq. (16), we see that the dependence on f f , and thus a sensitivity to the sign of either f or f , disappears in the unpolarized cross section of the hardscattering process uū → budbμ −ν μ bb. Let us note that the same arguments can easily be repeated for the amplitudes of the Feynman diagrams of Fig. 1, which dominate the tt H production through the gluon fusion process (2). We would like to stress here that all the above approximations are used for the sake of the argument in this section only and are not used to obtain the full results presented in Sect. 3.

Results
The calculation is performed in the framework of the SM, supplemented with the top-Higgs coupling derived from Lagrangian (4), with the use of carlomat [28], a general purpose program for the MC computation of the lowest order cross sections. The differential cross section of the process is calculated with the use of the following factorization formula: where x 1 and x 2 are the proton momentum fractions carried by partons i and j, respectively, s = x 1 x 2 s is the reduced center of mass energy squared, Q is the factorization scale and we take into account the following pairs of partons (i, j): (g, g), (u,ū), (ū, u), (d,d), (d, d) in the sum. We use MSTW LO parton distribution functions [29] at the factorization scale Q = m 2 t + j p 2 T j , where p T j is the transverse momentum of the final state quark or antiquark of process (18). The calculation is performed separately for the gluon fusion (2) and each of the quark-antiquark hardscattering processes (3). We use the same physical input parameters and cuts (3.2)-(3.7), with m cut bb = 20 GeV in (3.7), as in Ref. [19], and three different combinations of the scalar and pseudoscalar couplings of Lagrangian (4): Let us note that, in order to calculate the total cross section of process (18), a 20-fold phase space integral and a 2-fold integral over parton density functions must be performed, not to mention the additional 9-fold Monte Carlo (MC) integral that replaces the sum over particle helicities, without which the computation would not have been feasible in practice.
The differential cross sections of process (18) at the proton-proton center of mass energy of 14 TeV are plotted in Figs. 3, 4, 5, 6, 7 and 8 as functions of the rapidity and different angular variables of the final state muon, being referred to as the lepton. In Figs. 3, 4, 5, 6, 7 and 8, the left panels show the signal cross sections, which are computed with the signal tt H production diagrams of the hard-scattering processes (2) and (3), as described in the previous sections, and the right panels show the complete LO cross sections, which are computed with the complete set of the LO Feynman diagrams of each of the hard-scattering processes considered. In each of the figures, the SM cross section of process (18)   The individual contributions of the uū-andūu hardscattering processes to the tt H signal differential cross sections of process (18) at √ s = 14 TeV are plotted in Figs. 9, 10 and 11, as functions of the lepton rapidity, cosine of the lepton angle with respect to the beam and cosine of the lepton angle with respect to the Higgs boson in the laboratory (LAB) frame, respectively. The relative effects of the anomalous pseudoscalar coupling f in the plots of Figs. 9, 10 and 11 are approximately the same as in the full signal cross sections plotted in the right panels of Figs. 3, 4, and 8, respectively, and again there is practically no sensitivity to the sign of f . Taking into account the off-resonance background contributions to any of the quark-antiquark hard-scattering processes does not change this conclusion either, i.e., the shapes and relative effect of the anomalous coupling f remain practically the same for all the distributions considered.

Conclusions
We have complemented the analysis of the influence of the anomalous Higgs boson coupling to the top quark on the secondary lepton distributions in the process of associated production of the top-quark pair and Higgs boson in the proton-proton collisions at the LHC of Ref. [19] by taking into account contributions of the quark-antiquark annihilation hard-scattering processes. Although the gluon fusion mechanism dominates the tt H production through process (18) at √ s = 14 TeV, the contribution of quark-antiquark hard-scattering processes (3) is quite substantial and, therefore, should be taken into account in the analyses of data. Moreover, we have explained why the effects of the scalar and pseudoscalar anomalous couplings in the unpolarized cross section of the process are completely insensitive to the sign of either of them.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .