# Higgs production with large transverse momentum in hadronic collisions at next-to-leading order *footnote **footnote *Work partly supported by the EU Fourth Framework Programme ‘Training and Mobility of Researchers’, Network ‘Quantum Chromodynamics and the Deep Structure of Elementary Particles’, contract FMRX-CT98-0194 (DG 12 - MIHT) and the Swiss National Foundation.

## Abstract

Inclusive associated production of a light Higgs boson ($`m_\mathrm{H}m_t`$) with one jet in $`pp`$ collisions is studied in next-to-leading order QCD. Transverse momentum ($`p_\mathrm{T}30\mathrm{GeV}`$) and rapidity distributions of the Higgs boson are calculated for the LHC in the large top-quark mass limit. It is pointed out that, as much as in the case of inclusive Higgs production, the $`K`$-factor of this process is large ($`1.6`$) and depends weakly on the kinematics in a wide range of transverse momentum and rapidity intervals. Our result confirms previous suggestions that the production channel $`p+pH+\mathrm{jet}`$ $`\gamma +\gamma +`$ jet gives a measurable signal for Higgs production at the LHC in the mass range $`100140\mathrm{GeV}`$, crucial also for the ultimate test of the Minimal Supersymmetric Standard Model.

preprint:   ETH-TH/99-06 February 1999   

Recent results from LEP and the SLC indicate that the Higgs boson of the Standard Model might be light. A fit to the precision data has given the values $`m_\mathrm{H}=76_{47}^{+85}\mathrm{GeV}`$, corresponding to $`m_\mathrm{H}262\mathrm{GeV}`$ at the $`95\%`$ confidence level, whereas a direct search at LEP200 gives the lower limit as $`90\mathrm{GeV}m_\mathrm{H}`$ . In addition, a crucial theoretical upper limit exists on the mass of the light neutral scalar Higgs boson of the Minimal Supersymmetric Standard Model $`m_h130\mathrm{GeV}`$. It is, therefore, significant that one attempts to get the best possible signals in the light mass range of $`100\mathrm{GeV}m_\mathrm{H}140\mathrm{GeV}`$ at the LHC. Simulation studies carried out by ATLAS and CMS have shown that assuming a low integrated luminosity of $`310^4\mathrm{pb}^1`$, even in the case of the “gold-plated” decay channel into two photons, the signal significance ($`S/\sqrt{B}`$) is only around 5 . This conclusion depends on the value of the $`K`$-factor (a conservative value $`K`$=1.5 was used) and some plausible assumptions on the size of the background. The calculation of the next-to-leading order (NLO) corrections to the background is not yet complete  and the contribution of the NNLO subprocess $`gg\gamma \gamma `$ is large . The complete NNLO analysis is extremely laborious, but appears to be feasible. Actually, for the full numerical control the background has to be calculated in NNNLO, which is completely beyond the scope of presently available techniques. Fortunately, the ambiguity in the value of the background to the signal is suppressed by the square root appearing in the definition of the signal significance. This situation, which is not completely satisfactory, can be improved by studying the $`\gamma +\gamma +\mathrm{jet}(\mathrm{s})`$ final statesThe study of Higgs production in association with a jet was first suggested in the context of improving $`\tau `$ reconstruction in the $`\tau ^+\tau ^{}`$ decay channel . ; this offers several advantages. The photons are more energetic than in the case of the inclusive channel and the reconstruction of the jet in the calorimeter allows a more precise determination of the interaction vertex, improving the efficiency and mass resolution. Furthermore the existence of a jet in the final state allows for a new type of event selection and a more efficient background suppression. In addition, the necessary control of the background contributions can probably be already achieved by the inclusion of the NLO corrections (the matrix elements are known ). In a recent phenomenological study  it has been found that these advantages appear to be able to compensate the loss in production rates, provided one gets a large $`K`$-factor also for this process. The presentation of the NLO QCD corrections for this process is the main purpose of this letter.

The production process $`ggH`$ is given by loop diagrams in the Born approximation, since the gluons interact with the Higgs boson via virtual quark loops . The exact calculation of the NLO corrections is rather complex . Fortunately, the effective field theory approach  obtained in the large top mass limit with effective gluon–gluon–Higgs coupling gives an accurate approximation (with or without QCD corrections) with an error less than 5%, provided $`m_\mathrm{H}2m_\mathrm{t}`$ . It has been checked in LO, by an explicit calculation, that the approximation remains valid also for the production of Higgs bosons with large transverse momentum, provided both $`m_\mathrm{H}`$ and $`p_\mathrm{T}`$ are smaller than $`m_\mathrm{t}`$ . It is therefore plausible to assume that the approximation remains valid also if we include NLO QCD corrections. Recently, in this approximation and using the helicity method, the transition amplitudes relevant to the NLO corrections have been analytically calculated for all the contributing subprocesses (loop corrections  and bremsstrahlung ).

The available NLO matrix elements contain soft and collinear singularities and therefore do not allow for a direct numerical evaluation of the physical cross section. In the past few years, exploiting the universal structure of the soft and collinear contributions, several efficient algorithms have been suggested to obtain finite cross section expressions from the singular NLO matrix elements. We have used the method of ref.  and implemented it into a numerical Monte Carlo style program which allows to calculate any infrared-safe physical quantity for the inclusive production of a Higgs boson with one jet in NLO accuracy.

In this paper we report some of our results obtained for proton–proton collisions with $`\sqrt{S}=14\mathrm{TeV}`$. For the strong coupling constant at NLO (LO) we use the standard two-loop (one-loop) form with $`\mathrm{\Lambda }_{QCD}`$ set to the value used in the analysis of the parton distribution function under consideration. Our default choice for the factorization and renormalization scales is $`Q_0^2=(m_\mathrm{H}^2+p_\mathrm{T}^2)`$, where $`p_\mathrm{T}`$ is the transverse momentum of the Higgs boson. Here, unless stated, we consider the case of $`m_\mathrm{H}=120`$ GeV and $`p_\mathrm{T}>30`$ GeV in the kinematical region where the perturbative result can be applied without having to consider low-$`p_\mathrm{T}`$ resummation effects.

Most of our curves have been obtained with MRST (ft08a) parton distribution functions, but we will also show some results using CTEQ(4M) and GRV98 . To compare the leading with the NLO results, for consistency, we use the corresponding LO parton distributions from each set. We shall discuss only results for the inclusive production cross section of a Higgs boson with large transverse momentum, although as we mentioned above the Monte Carlo program allows to study any infrared-safe quantity, including the implementation of different jet algorithms and experimental cuts.

In Fig. 1(a) we show the $`p_\mathrm{T}`$ distribution of the NLO and LO cross sections using MRST parton densities at three different scales $`Q=\mu Q_0`$, with $`\mu =0.5,1,2`$. In this figure one can see three important points. First, the radiative corrections are large; second, there is a reduction in the scale dependence when going from LO to NLO; third, the improvement in the scale dependence is still not completely satisfactory. The same features can be observed in more detail in Fig. 1(b), where the LO and NLO cross sections integrated for $`p_\mathrm{T}`$ larger than 30 and 70 GeV are shown as a function of the renormalization/factorization scale. Both the LO and NLO cross sections increase monotonically with decreasing $`\mu `$ scale, down to the limiting value where perturbative QCD can still be applied. In NLO the scale dependence has a maximum and its position characterizes the stability of the NLO perturbative results. In our case, as a result of the very large positive radiative corrections, the position of the maximum is shifted down to small $`\mu `$ values where the perturbative treatment is not valid, indicating that the stability of the NLO result is not completely satisfactory. In the usual range of variation of $`\mu `$ from 0.5 to 2, however, the LO scale uncertainty amounts up to $`\pm 35\%`$, whereas at NLO this is reduced to $`\pm 20\%`$, indicating the relevance of the QCD corrections. This feature of the scale dependence in NLO is very much the same as the one in the case of inclusive Higgs production .

In Fig. 2 we show the ratio

$$K=\frac{\mathrm{\Delta }\sigma _{NLO}}{\mathrm{\Delta }\sigma _{LO}}$$
(1)   
of the next-to leading and the LO cross sections ($`K`$-factor) for the three different sets of parton distributions: MRST, CTEQ and GRV98, as a function of the transverse momentum and the rapidity of the Higgs boson. We can see that the $`K`$-factor is in the range 1.5–1.6 and it is almost constant (within 15% accuracy) for a large range of $`p_\mathrm{T}`$ and $`y`$. In the $`p_\mathrm{T}`$ distribution the variation never exceeds 10%, whereas it is a bit larger for the $`y`$ distribution at large $`|y|`$.

The ratios of the NLO cross sections

$$R=\frac{\mathrm{\Delta }\sigma _{\mathrm{CTEQ},\mathrm{GRV}}}{\mathrm{\Delta }\sigma _{\mathrm{MRST}}}$$
(2)   
computed by using CTEQ and GRV98 over the one obtained by using MRST parton densities are also shown. From there, it is possible to see that the differences in the $`K`$-factors basically come from variations in the LO cross sections, mostly because of the value of $`\mathrm{\Lambda }_{QCD}`$ used in each set.

We conclude that the properties of the $`K`$-factors found for large transverse momentum Higgs production are very similar to the ones obtained for the total inclusive Higgs production. They are about the same size, they show the same scale dependence, and the $`K`$-factor changes mildly with changing $`m_\mathrm{H},y`$ and $`p_\mathrm{T}`$. Its large value and its surprising independence from the kinematics might be interpreted as an evidence for some universal origin of the large radiative corrections . This requires further theoretical understanding.

We have not included in our analysis the contributions from electroweak reactions, which can increase the cross section about 10% when suitable cuts are applied . Nevertheless it is worth noticing that since QCD corrections to electroweak boson fusion are substantially smaller than the ones corresponding to gluon fusion, the significance of the electroweak contributions is reduced at NLO.

In Fig. 3 we show NLO cross section values for the physics signal $`p+pH+\mathrm{jet}\gamma +\gamma +\mathrm{jet}`$ as a function of the Higgs mass (with a reference value for the branching ratio given by Br$`(H\gamma \gamma )=2.1810^3`$ for $`m_\mathrm{H}=120`$ GeV ). For comparison, the cross section values of the physics signal $`p+pH\gamma +\gamma `$ are also shown. From there it is possible to see that the loss in production rate due to the transverse momentum cut of $`p_\mathrm{T}>30`$ GeV is less than a factor of 2 for the range of masses considered.

In conclusion, we have pointed out that, much as in the case of inclusive Higgs production, the cross section values of the associated production of a Higgs boson with a jet are increased by a $`K`$-factor of 1.5–1.6 given by NLO QCD radiative corrections. Our result confirms previous suggestions that the production channel $`p+pH+\mathrm{jet}`$ $`\gamma +\gamma +`$ jet gives a measurable signal for Higgs production at the LHC in the mass range $`100`$$`140\mathrm{GeV}`$, crucial also for the ultimate test of the Minimal Supersymmetric Standard Model.   

We are grateful to M. Spira for discussions. One of us (DdeF) would like to thank   
S. Frixione for helpful comments.