$k_T$-factorization approach to the Higgs boson production in $ZZ^*\to 4\ell$ channel at the LHC

TL;DR

This paper employs the $k_T$-factorization approach to calculate the differential cross section of Higgs boson production in the $ZZ^* o 4 ext{l}$ channel at the LHC, comparing results with experimental data and collinear factorization predictions.

Contribution

It introduces a $k_T$-factorization calculation for Higgs production in the $ZZ^*$ decay channel using CCFM evolution and off-shell matrix elements, providing a novel comparison with experimental and collinear results.

Findings

Results are close to NNLO+NNLL collinear predictions.

Calculations agree with experimental measurements at 8 TeV and 13 TeV.

The approach offers an alternative to collinear factorization for Higgs production.

Abstract

We calculated a differential cross section of the Higgs boson production in the $h\to ZZ^*\to 4\ell$ decay channel within the framework of $k_T$-factorization. Results are obtained using an off-shell matrix element for the $g^*g^*\to h\to ZZ^*\to 4\ell$ process together with Ciafaloni-Catani-Fiorani-Marchesini (CCFM) evolution equations for an unintegrated gluon distribution function. We have presented a comparison of our results with the latest experimental measurements at $\sqrt{S}$ = 8 TeV and $\sqrt{S}$ = 13 TeV from the ATLAS and CMS collaborations at the LHC. In addition to this, we have compared our results with the results from the collinear factorization formalism calculated up to next-to-next-to-leading order plus next-to-next-to-leading logarithm (NNLO + NNLL) accuracy obtained using the HRes code for the Higgs boson production in the gluon-gluon fusion process. Our estimates are consistently close to NNLO + NNLL results obtained using a collinear factorization formalism and are also in agreement with experimental measurements.

Tables

Table 1. Table 1: Total inclusive cross section ( σ t ​ o ​ t superscript 𝜎 𝑡 𝑜 𝑡 \sigma^{tot} ) for the Higgs boson production in gluon-gluon fusion channel.

$\sqrt{S}$ (TeV)	${\sigma }^{tot}$ (pb)
	($k_T$-factorization)	(Collinear factorization)
8	27.12	6.11
13	66.40	15.10

Equations

$${\mathcal{L}}_{ggh}=\frac{{\alpha }_s}{12\pi }(\surd 2G_F{)}^{1/2}{G}_{\mu \nu }^a{G}^{a\mu \nu }h,$$

$$T_{ggh}^{\mu \nu ,ab}(k_1,k_2)={\delta }^{ab}\frac{{\alpha }_S}{3\pi }(\surd 2G_F{)}^{1/2}[k_2^{\mu }{k}_1^{\nu }-(k_1\cdot k_2)g^{\mu \nu }].$$

$$\overline{\sum }{ϵ}_{\mu }^{\ast k}{ϵ}_{\nu }^k\simeq \frac{k_{\perp \mu }{k}_{\perp \nu }}{k_{\perp }^2}.$$

$$\overline{|\mathcal{M}{|}^2}=\frac{2}{9}\frac{{\alpha }_S^2}{{\pi }^2}\frac{m_Z^4}{v^4}\frac{[\hat{s}+({\sum }_{i=1}^2{\mathbf{k}}_{i\perp }{)}^2{]}^2}{[(\hat{s}-m_h^2{)}^2-m_h^2{\Gamma }_h^2]}{\cos }^2\phi \times \frac{[g_L^4+g_R^4](p_1\cdot p_3)(p_2\cdot p_4)+2g_L^2{g}_R^2(p_1\cdot p_4)(p_2\cdot p_3)}{[(2p_1\cdot p_2-m_Z^2{)}^2+m_Z^2{\Gamma }_Z^2][(2p_3\cdot p_4-m_Z^2{)}^2+m_Z^2{\Gamma }_Z^2]},$$

$$g_L=$$

$$g_R=$$

$$\sigma =\int \prod _{i=1}^2\frac{f_g(x_i,{\mathbf{k}}_{i\perp }^2,{\mu }_F^2)}{x_i^2{S}^2}d{\mathbf{k}}_{i\perp }^2\frac{d{\phi }_i}{2\pi }\prod _{f=1}^3{d}^2{\mathbf{p}}_{f\perp }d{y}_{f}d{y}_4\frac{\overline{|\mathcal{M}{|}^2}}{2^{12}{\pi }^5},$$

$$x_1=\sum _{f=1}^4\frac{|{\mathbf{p}}_{f\perp }|}{\sqrt{S}}{e}^{y_f},x_2=\sum _{f=1}^4\frac{|{\mathbf{p}}_{f\perp }|}{\sqrt{S}}{e}^{-y_f},$$

$$\sum _{i=1}^2{\mathbf{k}}_{i\perp }=\sum _{f=1}^4{\mathbf{p}}_{f\perp }.$$

$${\mathcal{L}}_{ggh}=\frac{{\alpha }_s}{12\pi v}{G}_{\mu \nu }^a{G}^{a\mu \nu }h,$$

$$\mathrm{i}\mathcal{M}={ϵ}_{\mu }^{\ast k_1}{ϵ}_{\nu }^{\ast k_2}\frac{[\mathrm{i}{T}_{ggh}^{\mu \nu ,ab}(k_1,k_2)]\mathrm{i}[\mathrm{i}{T}_{h\to Z{Z}^{\ast }\to 4\ell }(p_1,p_2,p_3,p_4)]}{(\hat{s}-m_h^2)+i{m}_h{\Gamma }_h},$$

$$T_{ggh}^{\mu \nu ,ab}(k_1,k_2)={\delta }^{ab}\frac{{\alpha }_s}{3\pi v}[k_2^{\mu }{k}_1^{\nu }-(k_1\cdot k_2)g^{\mu \nu }],$$

$$T_{h\to Z{Z}^{\ast }\to 4\ell }(p_1,p_2,p_3,p_4)=-2\frac{m_Z^2}{v}\frac{\bar{u}(p_1){\gamma }^{\rho }[g_L^{(1)}{P}_L+g_R^{(1)}{P}_R]v(p_2)}{[(p_1+p_2{)}^2-m_Z^2+i{m}_Z{\Gamma }_Z]}\times \frac{\bar{u}(p_3){\gamma }_{\rho }[g_L^{(2)}{P}_L+g_R^{(2)}{P}_R]v(p_4)}{[(p_3+p_4{)}^2-m_Z^2+i{m}_Z{\Gamma }_Z]},$$

$$\displaystyle g^{(i)}_{L}=\frac{g_{W}}{c_{W}}\Big{(}-\frac{1}{2}+s^{2}_{W}\Big{)},\quad g^{(i)}_{R}=\frac{g_{W}}{c_{W}}s^{2}_{W},\quad i=1,2.$$

$$\overline{|\mathcal{M}{|}^2}=$$

$$|T_{h\to Z{Z}^{\ast }\to 4\ell }(p_1,p_2,p_3,p_4){|}^2=64\frac{m_Z^4}{v^2}\times \frac{[g_L^4+g_R^4](p_1\cdot p_3)(p_2\cdot p_4)+2g_L^2{g}_R^2(p_1\cdot p_4)(p_2\cdot p_3)}{[(2p_1\cdot p_2-m_Z^2{)}^2+m_Z^2{\Gamma }_Z^2][(2p_3\cdot p_4-m_Z^2{)}^2+m_Z^2{\Gamma }_Z^2]},$$

$$|T_{ggh}^{\mu \nu ,ab}(k_1,k_2){|}^2=\overline{{\delta }^{ab}{\delta }^{ab}}\frac{1}{9}\frac{{\alpha }_s^2}{{\pi }^2}\frac{1}{v^2}\times \overline{\sum }{ϵ}_{\mu }^{\ast k_1}{ϵ}_{{\mu }^{\prime }}^{k_1}\overline{\sum }{ϵ}_{\nu }^{\ast k_2}{ϵ}_{{\nu }^{\prime }}^{k_2}[k_2^{\mu }{k}_1^{\nu }-(k_1\cdot k_2)g^{\mu \nu }]\times [k_2^{{\mu }^{\prime }}{k}_1^{{\nu }^{\prime }}-(k_1\cdot k_2)g^{{\mu }^{\prime }{\nu }^{\prime }}].$$

$$\begin{array}{rl}{k}_1=& z_1{P}_1+{\bar{z}}_1{P}_2+k_{1\perp },\\ k_2=& z_2{P}_1+{\bar{z}}_2{P}_2+k_{2\perp },\end{array}\text{where}\begin{array}{r}{k}_{1\perp }=(0,{\mathbf{k}}_{1\perp },0),\\ k_{2\perp }=(0,{\mathbf{k}}_{2\perp },0).\end{array}$$

$$P_{1,2}=\frac{\sqrt{S}}{2}(1,0,\pm 1),$$

$$\begin{array}{rl}{k}_1=& z_1{P}_1+k_{1\perp },\\ k_2=& {\bar{z}}_2{P}_2+k_{2\perp }.\end{array}$$

$$k_1\cdot k_2=z_1{\bar{z}}_2{P}_1\cdot P_2+k_{1\perp }\cdot k_{2\perp },$$

$$\hat{s}=(k_1+k_2{)}^2=2z_1{\bar{z}}_2{P}_1\cdot P_2-({\mathbf{k}}_{1\perp }+{\mathbf{k}}_{2\perp }{)}^2.$$

$$\overline{\sum }{ϵ}_{\mu }^{\ast k}{ϵ}_{\nu }^k\simeq \frac{k_{\perp \mu }{k}_{\perp \nu }}{k_{\perp }^2},$$

$$\overline{\sum }{ϵ}_{\mu }^{\ast k_1}{ϵ}_{{\mu }^{\prime }}^{k_1}$$

$$=$$

$$\overline{|\mathcal{M}{|}^2}=\frac{2}{9}\frac{{\alpha }_S^2}{{\pi }^2}\frac{m_Z^4}{v^4}\frac{[\hat{s}+({\mathbf{k}}_{1\perp }+{\mathbf{k}}_{2\perp }{)}^2{]}^2}{(\hat{s}-m_h^2{)}^2+m_h^2{\Gamma }_h^2}{\cos }^2\phi \times \frac{[g_L^4+g_R^4](p_1\cdot p_3)(p_2\cdot p_4)+2g_L^2{g}_R^2(p_1\cdot p_4)(p_2\cdot p_3)}{[(2p_1\cdot p_2-m_Z^2{)}^2+m_Z^2{\Gamma }_Z^2][(2p_3\cdot p_4-m_Z^2{)}^2+m_Z^2{\Gamma }_Z^2]}.$$

$$p=(E_T\cosh y,{\mathbf{p}}_T,E_T\sinh y),$$

$$\frac{d^3\mathbf{p}}{p_{{}_0}}=d^2{\mathbf{p}}_{T}dy.$$

$$p=(E_T\cosh y,p_T\cos \phi ,p_T\sin \phi ,E_T\sinh y).$$

$$\frac{d^3\mathbf{p}}{p_{{}_0}}=p_{T}d{p}_{T}d\phi dy.$$

$$\displaystyle\sigma=\!\!\int\!\!\frac{1}{x_{1}x_{2}S}\prod_{i=1}^{2}\frac{dx_{i}}{x_{i}}f_{g}(x_{i},{\bf k}^{2}_{i\perp},\mu^{2}_{F})\frac{d^{2}{\bf k}_{i\perp}}{\pi}\times\prod_{f=1}^{4}\frac{d^{3}{\bf p}_{f}}{(2\pi)^{3}2p_{f0}}\overline{|{\cal M}|^{2}}\times(2\pi)^{4}\delta^{(4)}\bigg{(}\sum_{i=1}^{2}k_{i}-\sum_{f=1}^{4}p_{i}\bigg{)}.$$

$$\int \frac{d^2{\mathbf{k}}_{i\perp }}{\pi }=\int \frac{k_{i\perp }d{k}_{i\perp }d{\phi }_i}{\pi }=\int d{\mathbf{k}}_{i\perp }^2\frac{d{\phi }_i}{2\pi }.$$

Full text

I Introduction

The Higgs boson discovery at the LHC, by ATLAS and CMS collaborations Aad et al. (2012); Chatrchyan et al. (2012), has enabled experimental measurements to be taken to investigate its properties. The ATLAS and CMS collaborations have performed an improved measurement of the Higgs boson mass, considering an invariant mass spectra of the $h\to\gamma\gamma$ and $h\to ZZ^{*}\to 4\ell$ decay channels Aad et al. (2014a); Khachatryan et al. (2015a). Further studies of spin and the parity quantum number of the Higgs boson have established that it is a neutral scalar boson, with a mass equal to 125.09 GeV, rather than a pseudoscalar boson Chatrchyan et al. (2013); Aad et al. (2013); Patrignani et al. (2016). Its coupling strength to vector bosons and to fermions is studied by analyzing various decay modes of the Higgs boson Khachatryan et al. (2015a, b); Aad et al. (2015, 2016). Establishing various aspects of the Higgs boson’s properties and coupling strength allows us to study other aspects of it.

A dominant channel for the inclusive Higgs boson production at the LHC is gluon-gluon fusion Georgi et al. (1978); Graudenz et al. (1993); Spira et al. (1995). Hence, the Higgs boson production at the LHC can be effectively used to understand the gluon dynamics inside a proton. The gluon density $xf_{g}(x,\mu^{2}_{F})$ in a proton is a function of the Bjorken variable $x$ and the hard scale $\mu^{2}_{F}$. The scale evolution of parton densities, in general, is described using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation Gribov and Lipatov (1972); Lipatov (1975); Altarelli and Parisi (1977); Dokshitzer (1977), where large logarithmic terms proportional to $\ln\mu^{2}_{F}$ are resummed up to all orders.

The factorization theorem in perturbative quantum chromodynamics (pQCD) allows us to write a convolution of the matrix element of the short distance process and the universal parton distribution functions to obtain an inclusive cross section for a given scattering process Collins et al. (1989). The QCD collinear factorization theorem is based on the collinear approach where the parton distribution function depends on the longitudinal momentum fraction $x$ and the hard scale $\mu^{2}_{F}$. The Higgs boson production cross section at leading order and higher order QCD corrections to it, up to next-to-next-to-next-to-leading order (N3LO), have been computed within the collinear factorization framework Dawson (1991); Djouadi et al. (1991); Harlander and Kilgore (2002); Anastasiou and Melnikov (2002); Ravindran et al. (2003); Anastasiou et al. (2015). However, it should be noted that the NNLO and N3LO results that were obtained so far are using an effective theory and in the heavy top quark mass limit. Study of the Higgs boson’s transverse momentum spectrum resummed at NNLL accuracy is shown in de Florian et al. (2012); Bozzi et al. (2003, 2006). Recently, the state-of-the-art predictions for the Higgs boson’s transverse momentum at the LHC, at next-to-next-to-next-to-leading-logarithmic accuracy (N3LL) matched, at NNLO is presented in Ref. Bizon et al. (2017).

For the inclusive Higgs boson production at the LHC, the longitudinal momentum fraction of the incident gluons is small ($x_{1}x_{2}\sim$ 0.0089-0.0175). This domain of small longitudinal momentum fraction ($x$) is still in the perturbative regime where is it expected that collinear factorization should break down because the large logarithmic term proportional to $1/x$ becomes dominant Andersson et al. (2002); Andersen et al. (2004, 2006). The contribution from the terms proportional to $1/x$ is taken into account in the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation Kuraev et al. (1976, 1977); Balitsky and Lipatov (1978). An unintegrated parton densities (uPDFs) obeying BFKL evolution, convoluted with an off-shell matrix element within a generalized factorization is called $k_{T}$-factorization Gribov et al. (1983); Levin et al. (1991); Catani et al. (1990a, 1991); Collins and Ellis (1991). The evolution equation is valid for both small $x$ and large $x$ is given by the Ciafaloni-Catani-Fiorani-Marchesini (CCFM) evolution equation Ciafaloni (1988); Catani et al. (1990b, c); Marchesini (1995). CCFM evolution is equivalent to BFKL evolution in the limit of very small $x$ and is equivalent to the DGLAP evolution for a large $x$ region.

In this work, we have not implemented the reggeized parton approach Lipatov and Vyazovsky (2001); Bogdan and Fadin (2006); Hentschinski (2012); Chachamis et al. (2012); Hentschinski and Sabio Vera (2012) based on the Lipatov effective action formalism, which ensures a gauge invariance of the off-shell amplitude Lipatov (1995, 1997). However, our investigation is based on the assumption that the off-shell partonic amplitudes being gauge-invariant in a small-$x$ limit. The approach, which is based on Lipatov effective action formalism that has been employed recently in the calculation using $k_{T}$-factorization approach for the inclusive prompt photon production at LHC Lipatov and Malyshev (2016).

In this paper, the inclusive Higgs boson production within the $k_{T}$-factorization approach, together with CCFM evolution equations have been studied and demonstrated importance of higher order corrections included within the $k_{T}$-factorization Lipatov and Zotov (2005). The authors of Ref. Lipatov et al. (2014) have shown that $k_{T}$-factorization gives a description of an experimental data from ATLAS experiment for the differential cross section of the Higgs boson production in the diphoton decay channel. They have calculated a leading order (LO) matrix element for the partonic subprocess $gg\to h\to\gamma\gamma$ considering gluons to be off-shell. Inclusive Higgs boson production analysis based on off-shell gluon-gluon fusion, and considering $H\to\gamma\gamma$, $H\to ZZ^{*}\to 4\ell$ (where $\ell=e,\mu$) and $H\to W^{+}W^{-}\to e^{\pm}\mu^{\mp}\nu\bar{\nu}$ decay channel is given in Ref. Abdulov et al. (2018).

The ATLAS and CMS collaboration at the LHC presented a measurement of a fiducial differential cross section of the Higgs boson decay into four-leptons at $\sqrt{S}$ = 8 TeV Aad et al. (2014b); Khachatryan et al. (2016) and $\sqrt{S}$ = 13 TeV Aaboud et al. (2017). We compare the results obtained using the $k_{T}$-factorization approach with the recent ATLAS and CMS data. We have evaluated the off-shell matrix element for the partonic subprocess $g^{*}g^{*}\to h\to ZZ^{*}\to 4\ell,\ell=e,\mu$. Convolution of the off-shell matrix element of partonic subprocess with CCFM uPDFs Jung (2004) is used to obtain a differential cross section.

This article is organized as follows: We discuss in detail the formalism behind our study and the necessary expressions for further numerical analysis in Section II. In Section III, we give the results of our numerical simulation. Here we also discuss the details of the analyses that have gone into the study. Finally, add we conclude and draw inferences from the analysis in Section IV.

II Formalism

In the present section, we briefly discuss the formalism we have used in our study. The details are in the Appendix A and Appendix B. In Section I, we have mentioned that to explore the effects of $k_{T}$-factorization, we need to take the initial state partons to be off-shell. In calculating the off-shell matrix element for the process $g^{*}\!g^{*}\!\to\!h\!\to\!ZZ^{*}\!\to\!4\ell$

(see Fig. 1), we have used the effective field theory approach. The effective Lagrangian in the large top quark mass limit, $m_{t}\to\infty$, for the Higgs boson coupling to gluon is Ellis et al. (1976); Shifman et al. (1979)

[TABLE]

where $\alpha_{s}$ and $G_{F}$ are the strong and Fermi coupling constants, respectively. $G^{a}_{\mu\nu}$ is the gluon field strength tensor and $h$ is the Higgs scalar field. The effective $ggh$ triangle vertex (see Eq. 10) thus becomes

[TABLE]

The non-zero transverse momentum for an initial gluon leads to the corresponding polarization sum Gribov et al. (1983); Levin et al. (1991)

[TABLE]

Using Eqs. 2 and 3, we derived the off-shell matrix element for the hard scattering process. The matrix element thus obtained is given by Eq. 23 as

[TABLE]

where

[TABLE]

Here $\Gamma_{h}$ and $\Gamma_{Z}$ are the total decay widths of the Higgs boson and $Z$ boson, respectively. ${\bf k}_{i\perp}\!\!$ are the intrinsic transverse momenta of the initial gluons. $\varphi$ is the azimuthal angle between ${\bf k}_{1\perp}$ and ${\bf k}_{2\perp}$. $m_{h}$ and $m_{Z}$ are the Higgs boson and $Z$ boson masses, respectively. The partonic centre of mass energy is denoted by $\hat{s}$. $\theta_{W}$ and $g_{W}$ are the weak mixing angle and the coupling of weak interaction, respectively.

Finally, we arrived at the hadronic cross section for the off-shell hard scattering amplitude of Eq. 4 within the framework of $k_{T}$-factorization as 111The processes contributing to $4$ lepton final states are $e^{+}e^{-}e^{+}e^{-}$, $\mu^{+}\mu^{-}\mu^{+}\mu^{-}$, $e^{+}e^{-}\mu^{+}\mu^{-}$ and $\mu^{+}\mu^{-}e^{+}e^{-}$. We have added these contributions with proper weight. The first two processes have two pairs of identical particles. Hence their phase space has to be multiplied by a factor $(1/2)\times(1/2)=1/4$ to get rid of overcounting. (see Eq. 37)

[TABLE]

with the longitudinal momentum fractions $x_{1}$ and $x_{2}$ of initial gluons to be

[TABLE]

and the transverse momenta

[TABLE]

In Eq. 5, $\varphi_{1,2}$ are the azimuthal angle of ${\bf k}_{1\perp,2\perp}$. $y$ and ${\bf p}_{\perp}\!\!$ are rapidities and transverse momenta of the final state leptons, respectively. The hadronic centre of mass energy is denoted by $S$.

III Results and Discussion

With all the calculational tools at our disposal, we proceed to perform a numerical calculation using Eq. 5 together with the off-shell hard scattering amplitude given in Eq. 4. We estimate the cross section of the Higgs boson production as a function of transverse momentum ($p_{T}$) and rapidity ($y$) of the Higgs boson in the four-lepton decay channel. Results are obtained using CCFM A0 the set of uPDFs Jung (2004) which is commonly used for such phenomenological studies. Recently a fit to a high precision data from deeply inelastic scattering at the HERA is performed using a $k_{T}$-factorization and CCFM evolution Hautmann and Jung (2014). A transverse momentum dependent gluon density function including experimental and theoretical uncertainties were obtained. The application of these unintegrated gluon densities to vector boson + jet production processes at LHC is given in Ref. Dooling et al. (2014). Unintegrated gluon densities including experimental and theoretical uncertainties are given in the CCFM JH2013-set in TMDlib library Hautmann et al. (2014); Connor et al. (2016). We have used the TMDlib library to calculate our results using CCFM JH2013-set. For our phenomenological study, we have used the CCFM JH2013-set2 which is determined from the fit to both structure functions $F^{(charm)}_{2}$ and $F_{2}$ data whereas CCFM JH2013-set1 is determined from the fit to inclusive $F_{2}$ data only.

Total decay width and mass of the Higgs boson is set to be equal to 4.0 MeV and 125.09 GeV, respectively Patrignani et al. (2016). We have implemented kinematical cuts on the rapidity and transverse momentum of leptons used by ATLAS and CMS experiments in their measurements. For the ATLAS experiment, the absolute value of rapidity is $|\eta|<2.5$ and the leading transverse momentum of the lepton is $p_{T}<20$ GeV. The transverse momenta of sub-leading leptons are $p_{T}<15,10,7$ GeV. Similarly for the CMS experiment, the absolute value of rapidity is $|\eta|<2.5$ and the ordered transverse momenta of the leptons are $p_{T}<20,10,7,7$ GeV. The cross section in Eq. 5 depends on the renormalization and factorization scales $\mu_{R}$ and $\mu_{F}$, respectively. The scale uncertainty in the cross section is estimated by varying the scale between $\mu_{R}=\mu_{F}=m_{h}/2$ and $\mu_{R}=\mu_{F}=2m_{h}$.

We have also calculated a total inclusive cross section for the Higgs boson production with $k_{T}\to 0$ and averaged over an azimuthal angle of the Higgs boson. This result of the inclusive cross section is equivalent to the cross section obtained using the collinear factorization approach at LO while using collinear parton densities. We have used the Martin-Stirling-Thorne-Watt (MSTW) set Martin et al. (2009) for collinear parton densities. We have also obtained results of the total inclusive cross section for the Higgs production with $k_{T}$-factorization formalism using the CCFM JH2013-set2 of uPDFs.

In Table 1, we have given our results for the total inclusive cross section for the Higgs production, using both collinear and $k_{T}$-factorization framework. Our results, for total inclusive cross section, obtained using collinear approach, are consistent with the results obtained in Ref. Anastasiou et al. (2012) at $\sqrt{S}$ = 8 TeV. The results obtained with $k_{T}$-factorization is close to next-to-leading order (NNLO) results given in Ref. Anastasiou et al. (2012) at $\sqrt{S}$ = 8 TeV. The cross section estimates given here are for a gluon-gluon fusion process only. The inclusive cross section for the Higgs boson production can be obtained using a hadron level Monte Carlo event generator called CASCADE Jung et al. (2010). CASCADE uses the CCFM evolution equation in the initial state with the off-shell parton level matrix element.

We have presented our results in Figs. 2, 3 and 4. Figs. 2, 3 and 4 shows the result of the differential cross section for the Higgs boson production in the four-lepton decay channel at $\sqrt{S}$ = 8 TeV and $\sqrt{S}$ = 13 TeV, respectively. We have compared our results of the differential cross section obtained using the $k_{T}$-factorization approach with experimental measurements from the ATLAS at $\sqrt{S}$ = 8 TeV, 13 TeV and CMS collaboration at $\sqrt{S}$ = 8 TeV Aad et al. (2014b); Khachatryan et al. (2016). The solid (red) and dashed (purple) histogram corresponds to our results obtained using the CCFM JH2013-set2 set of uPDFs and CCFM A0 set of uPDFs respectively. We also see that the results obtained using CCFM JH2013-set2 has a better agreement with experimetal measurements than A0 set.

Our results are plotted against state-of-the-art results for the cross section calculated up to next-to-leading order plus next-to-leading logarithm (NLO + NLL) and next-to-next-to-leading order plus next-to-next-leading logarithm (NNLO + NNLL) obtained using the HRes tool de Florian et al. (2012); Grazzini and Sargsyan (2013) within the collinear factorization framework. Our results of both the differential cross section in $p_{T}$ and $y$ using the $k_{T}$-factorization framework with CCFM unintegrated PDFs are consistently close to NNLO + NNLL results at $\sqrt{S}$ = 8 TeV and $\sqrt{S}$ = 13 TeV. This can be explained considering the fact that the main part of higher order corrections are included in the $k_{T}$- factorization approach Lipatov and Zotov (2005); Hautmann (2002); Ryskin et al. (2001). For the $p_{T}$ distribution, we are using the convention that NLO + NLL and NNLO + NNLL results are labelled as LO + NLL and NLO + NNLL, respectively considering that the $p_{T}$ distribution is non-zero at NLO.

IV Conclusions

In this paper, add we present a phenomenological study of the Higgs boson production in the four-lepton decay channel within the $k_{T}$-factorization framework. Here CCFM unintegrated parton densities were convoluted with the hard matrix element, considering initial gluons to be off- shell. We present a comparison of our results with experimental measurements. Our results are evaluated using same experimental conditions (i.e., the same $p_{T}$ and $y$ cuts were used for our estimates as given by the experimental results) for both the ATLAS and CMS at $\sqrt{S}$ = 8 TeV and at $\sqrt{S}$ = 13 TeV, respectively.

Further comparison of our estimates with the state-of-the-art results of a differential cross section within collinear factorization up to NLO + NLL and NNLO + NNLL obtained using HRes code is presented. We have also estimated a total inclusive cross section for the Higgs boson production within both collinear factorization and $k_{T}$-factorization framework. Within $k_{T}$-factorization approach, we have compared the results obtained using CCFM JH2013-set2 and CCFM A0 uPDF set. Our results for the differential cross section with CCFM JH2012-set2 are close to the NNLO + NNLL results obtained using the HRes tool.

Our results show that the observed $p_{T}$ distribution of the final state can be generated at leading order subprocesses, using unintegrated gluon distributions. Moreover, gluons in the initial state have finite transverse momenta, which results in the transverse momenta of the final state. The total inclusive cross section estimated using $k_{T}$-factorization is close to the NNLO results obtained using collinear factorization. The main reason for this behavior is that the main part of higher order correction in collinear pQCD is already included in the $k_{T}$-factorization Watt et al. (2004) framework.

The higher order corrections within $k_{T}$-factorization at the parton level would be an interesting study, to see any additional effect. The cross section for the Higgs boson production has been calculated using a mixture of LO and NLO partonic diagrams and unintegrated PDFs from the $k_{T}$- factorization. Considering the effect of the transverse momentum of the initial gluon on the transverse momentum distribution of the final state, our study, add as well as further studies in this direction could impose constraints on uPDFs of gluons.

Appendix A Amplitude of $g^{*}\!g^{*}\!\to\!h\!\to\!ZZ^{*}\!\to\!4\ell$

In this appendix, add we give details of our calculations that went into our analysis. Fig. 1 shows the assignment of the momenta. In the limit of large top quark mass $m_{t}\to\infty$, the effective Lagrangian for the Higgs boson coupling to gluons given in Eq. 1 can be written as

[TABLE]

where $v$ is the vacuum expectation value of the scalar field. The amplitude of the process $g^{*}\!g^{*}\!\to\!h\!\to\!ZZ^{*}\!\to\!\ell\bar{\ell}\ell^{\prime}\bar{\ell}^{\prime}$ is

[TABLE]

where

[TABLE]

where

[TABLE]

Hence we shall denote them only by $g_{L}$ and $g_{R}$ hereafter.

Now using the above expressions into Eq. 9 we can get $\overline{|{\cal M}|^{2}}$ for the process as

[TABLE]

where

[TABLE]

whereas

[TABLE]

To calculate the gluon part of the above amplitude we use the Sudakov decomposition of momenta as followed in Ref. Catani et al. (1991); Collins and Ellis (1991). Accordingly, we can write the gluon momenta $k_{1,2}$ as

[TABLE]

In the above equation, $k_{1\perp,2\perp}$ are vectors transverse to the momenta $P_{1,2}$ of the incoming hadrons. It is convenient to take both $P_{1,2}$ to be light-like (i.e. $P^{2}_{1,2}=0$). Because of this, $P_{1,2}$ are slightly different from the actual momenta of the incoming hadrons. Let us take $P_{1,2}$ to be as follows

[TABLE]

where $\sqrt{S}$ is the centre of mass energy of the hadrons. In the high energy limit, the introduction of strong ordering in longitudinal momenta gives

[TABLE]

Therefore using Eqs. 17 and 18 we get

[TABLE]

In the last step, add we have used the definition of $k_{1\perp,2\perp}$ introduced in Eq. 16.

With the help of the above expressions, we are now ready to calculate the gluon part of Eq. 15. Using the polarization sum for the off-shell gluons given by

[TABLE]

we can calculate

[TABLE]

Therefore putting the expression of Eq. 22 into Eq. 15 we get the total spin averaged squared amplitude of the process $g^{*}\!g^{*}\!\to\!h\!\to\!ZZ^{*}\!\to\!\ell\bar{\ell}\ell^{\prime}\bar{\ell}^{\prime}$ as

[TABLE]

Appendix B Phase Space Calculation

In this appendix, we give details of the phase space calculations related to our analysis. Let us express the 4-momentum $p$ in a 3-component vector in terms of the transverse momentum ${\bf p}_{T}$ and the rapidity $y$ as follows

[TABLE]

where the transverse energy $E_{T}=\sqrt{{\bf p}^{2}_{T}+m^{2}}$ and $m$ is the mass. We can write the measure of the phase space integration as

[TABLE]

Let us further express the 4-momentum $p$ in a 4-component vector in terms of $p_{T}(\equiv|{\bf p}_{T}|)$, $y$ and the azimuthal angle $\varphi$ as

[TABLE]

Now the measure of integration over the phase space becomes

[TABLE]

The above can be easily understood if we express the 2-component vector ${\bf p}_{T}$ in polar coordinates $(p_{T}\cos\varphi,p_{T}\sin\varphi)$. Then we can write $d^{2}{\bf p}_{T}=p_{T}dp_{T}d\varphi$.

Now let us turn our attention to the phase space integration of a 4-body final state process in a hadron collider. The hadronic cross section in the $k_{T}$-factorization approach for the off-shell hard scattering amplitude of Appendix A is (see Eq. 2.1 of Ref. Szczurek et al. (2014))

[TABLE]

Introducing the azimuthal angles $\varphi_{1,2}$ of the off-shell gluons into Eq. 28 we can write the gluon phase space as

[TABLE]

Let us write the 4-momenta of the partons, using Eq. 18, as

[TABLE]

Here $z_{1}\to x_{1}$ and $\bar{z}_{2}\to x_{2}$. Also the 4-momenta of the leptons are

[TABLE]

Using the above designation of the gluon and lepton momenta we can simplify the delta function of Eq. 28 as

[TABLE]

Putting the results from Eqs. 29 and 33 into Eq. 28 and integrating over $x_{1}$ and $x_{2}$ we get

[TABLE]

where

[TABLE]

Now integration over ${\bf p}_{4\perp}$ in Eq. 34 gives us

[TABLE]

In the last step above we have used $d^{2}{\bf p}_{f\perp}=2\pi p_{f\perp}dp_{f\perp}=\pi d{\bf p}^{2}_{f\perp},(f=1,2,3)$.

Acknowledgements

The authors would like to extend their sincere gratitude to V. Ravindran for continuous support and fruitful discussions during the course of this work and A. Lipatov for providing their codes of a similar study for diphoton distributions. VR would like to thank Giancarlo Ferrera for useful discussions on HRes tool and related topics. RI would also like to acknowledge the hospitality provided by the Institute of Mathematical Sciences, Chennai, India where a part of the work has been done. RI thank the SERB-DST, India for the research grant EMR/2015/000333.