Parton distribution functions in Monte Carlo factorisation scheme

A next step in development of the KrkNLO method of including complete NLO QCD corrections to hard processes in a LO parton-shower Monte Carlo is presented. It consists of a generalisation of the method, previously used for the Drell–Yan process, to Higgs-boson production. This extension is accompanied with the complete description of parton distribution functions in a dedicated, Monte Carlo factorisation scheme, applicable to any process of production of one or more colour-neutral particles in hadron–hadron collisions.

The main advantage of the KrkNLO method with respect to other, older methods of matching the fixed-order NLO calculations with PSMCs (MC@NLO and POWHEG) is its simplicity. This simplicity stems from the fact that the entire NLO corrections are implemented using a simple positive multiplicative MC weight. However, in order to profit from it, one has to use in the KrkNLO method parton distribution functions (PDFs) in a special, so-called Monte Carlo (MC) factorisation scheme and PSMC has to fulfil some minimum quality criteria. Most of modern PSMCs [9][10][11][12][13] are good enough for the KrkNLO method.
Construction of PDFs in the MC factorisation scheme (FS) has evolved step by step: in Ref. [1] it was defined for gluonstrahlung only (albeit for two different processes, DY and DIS). In Ref. [2], the KrkNLO PDFs in the MC FS were defined and numerically constructed including also gluon to quark transitions/splittings, relevant for the complete NLO corrections in the DY process, which at the LO level has only quarks and antiquarks in the initial state. PDFs in the MC scheme in Ref. [2] were defined in terms of the standard MS PDFs, and constructed numerically by transforming the MS PDFs into MC-scheme PDFs, before they were plugged into PSMC used in the KrkNLO method.
However, in Ref. [2] certain elements in the transition matrix K , transforming the MS PDFs into the MC-scheme PDFs could be omitted, because they were not relevant (i.e. of a NNLO class) for the DY process. These elements of the transition matrix have to be added for any process with initialstate gluons, such as the Higgs-boson production elaborated in the present work. They will be defined and applied in the following, such that the complete transition matrix K transforming the MS PDFs into the MC-scheme PDFs will be specified for the first time. It will be argued that PDFs in such a MC-scheme can serve in the KrkNLO method for any process at a hadron-hadron collider in which a colour-neutral single or multiple system of heavy particles is produced. For other processes, with one or more coloured partons in the final state at LO level, the KrkNLO method with PDFs in the MC scheme may also work, but this subject is reserved for the forthcoming publications.
The MC factorisation scheme is a complete scheme, such that NLO coefficient functions for any hard process under consideration are known, hence PDFs in the MC FS can be fitted directly to experimental DIS and DY data. However, at present, we obtain them from PDFs in the MS scheme and leave out direct fitting to data for the future developments.
On the methodological side, as seen in Refs. [1,2], the essence of the KrkNLO method is that certain NLO correction terms in an unintegrated/exclusive form present in the MS scheme, which are proportional to unphysical Dirac-delta terms in transverse momentum of emitted real partons, are removed in the KrkNLO methodology by means of redefinition of PDFs from the MS to MC scheme. These 'pathological' terms are preventing the use of a simple multiplicative MC weight for implementing NLO corrections in the MS scheme in real-emission phase space, and they complicate implementation of the MC@NLO and POWHEG methods. These peculiar terms can be determined and calculated either by means of studying the NLO corrections to hard process (coefficient functions), or, alternatively, by means of integrating soft-collinear counter-terms (similar to these in the Catani-Seymour method [14]), which define the MC-scheme PDFs in d = 4 + 2ε dimensions. 1 We are going to calculate them using both methods, obtaining the same results.
Last but not least, the NLO calculations for the DY process of Ref. [2] were also compared with the NNLO calculations of MCFM [15], leading to the conclusion that they are closer to the latter than the results of the MC@NLO and POWHEG methods.
The outline of the paper is the following: in Sect. 2 the KrkNLO method is characterised briefly. In Sect. 3 all distributions needed for implementation of the KrkNLO method for Higgs-boson production in gluon-gluon fusion are elaborated, including also many analytical crosschecks and a necessary update of the virtual corrections in soft-collinear counter-terms used in Ref. [2] for the Z /γ * (DY) process. Section 4 presents numerical results for PDFs in the MC scheme. Then the first numerical results for the total cross section from the KrkNLO method for the Higgs production at the LHC are shown in Sect. 5. Finally, in Sect. 6 we summarise the paper and discuss future prospects of our work. In Appendix A the formulae for the NLO coefficient functions of the DIS process in the MC scheme are provided.

The method
The KrkNLO method was formulated in a few variants. For instance, in the version of Ref. [1], the MC weight implementing the NLO corrections sums the contributions from all relevant partons generated in PSMC next to the hard process "democratically", such that it works equally well for PSMCs based on angular ordering or virtuality ordering, contrary to POWHEG which requires adding extra gluons to a PSMC event. In the present work, we are going to follow the variant of KrkNLO discussed in Ref. [2], in which the NLOcorrecting MC weight uses only one parton, the one closest to the hard process in the transverse momentum, that is, the first parton generated in the backward evolution (BEV) in the PSMC algorithm with k T -ordering.
In any case, in the KrkNLO method, the entire event of PSMC is preserved and reweighted, contrary to POWHEG and MC@NLO where the parton attributed to the hard process is generated outside PSMC and, only later on, the remaining partons are provided by PSMC. Obviously, this puts certain minimum quality requirements on the PSMC: (i) the first parton in the BEV algorithm has to be generated with the distribution which has a correct soft and collinear limit and (ii) its phase space in momentum and flavour space has to be covered completely, without empty regions. Luckily, the above requirement is fulfilled by all modern PSMCs for initial-state emissions discussed in this work.
It is worth to comment in advance on the apparent use in the following of the soft-collinear counter-terms (dipoles) of the Catani-Seymour (CS) subtraction scheme [14]. Their role is twofold: (1) the CS dipoles serve us as a useful benchmark, as they provide a reference model for QCD distributions of real emissions featuring the exact soft and collinear limits and (2) the CS scheme helps us in a proper inclusion of the NLO virtual corrections. However, let us point out immediately an important difference between the MC and CS scheme: the CS dipoles do not include virtual corrections, while soft-collinear counter-terms (SCCTs) of the KrkNLO do include them, albeit not calculated from Feynman diagrams, but deduced from PDF momentum sum rules. The role of the SCCTs in the KrkNLO methodology is also much richer than that of the dipoles in the CS scheme-our SCCTs not only provide subtractions of soft-collinear singularities in real-emission phase space, but they are also used to define PDFs in the MC factorisation scheme. Moreover, their sums are required to coincide with the corresponding sums of real-parton distributions in PSMC. 2 Fig. 1 The LO Feynman diagram for the process of Higgs boson production in gluon-gluon fusion. The effective vertex (black dot) corresponds to a quark loop with summation over all quarks, in which the top-quark mass is set to infinity while the masses of the other quarks are set to zero

Higgs production in gluon-gluon fusion
In the following we are going to collect all distributions needed for implementation of the KrkNLO method for the gluon-fusion Higgs production in hadron-hadron collisions. Elements of the matrix transforming PDFs from the MS to MC scheme will also be obtained as a byproduct.
We start necessarily from the leading order (LO) process see Fig. 1, where Q = p 1 + p 2 . The LO matrix element squared, in the limit m t → ∞ and neglecting all other quarks contributions, reads where v 2 = ( √ 2G F ) −1 is the Higgs vacuum expectation value (VEV) squared. Hence, the LO total cross section takes the form For all NLO subprocesses (channels) where a and b are incoming partons (gluons and/or quarks), while c is an outgoing parton (quark or gluon) we shall use the same parametrisation of the kinematics in terms of the following Sudakov variables: For the gg-channel NLO subprocess shown in Fig. 2, the matrix element squared reads (3.7) For the qg-channel NLO subprocess shown in Fig. 3, one obtains Finally, for the qq channel This last process, unlike the previous ones, is not generated by the backward-evolution PSMC starting from the gg → H hard process, hence in the KrkNLO method, its contribution cannot be treated by NLO-reweighting of events generated by the main branch of the LO PSMC algorithm. It has to be added as an extra tree-level LO process to PSMC. Moreover, it is free of collinear and soft singularities. This poses no problem as most of present-day PSMCs implement such a process.

CS dipoles and MC matrix elements
In the following we shall elaborate mainly on the hadronhadron collision producing the Higgs boson or Z /γ * (Drell-Yan process). However, components of the KrkNLO method defined here will also be applicable to any LO process a +ā → X and the corresponding a + b → X + c where Fig. 2 The NLO Feynman diagrams for real-parton radiation in the process of Higgs-boson production in gluon-gluon fusion: the gg channel Fig. 3 The NLO Feynman diagram for real-parton radiation in the process of Higgs-boson production in gluon-gluon fusion: the gq channel X = H, Z /γ , W ± , Z Z, W + W − or any other colour-neutral heavy object; a, b = q,q, g are initial coloured partons and c is an additional parton emitted at the NLO level. In the following formulation of the KrkNLO-method components, the CS dipoles will serve us as useful auxiliary objects. They are formed by an initial-state (on-shell) emitter a from one hadron and a spectator parton b from another hadron, 3 see Fig. 5. Following closely the notation of the CS work [14], the emitter a splits into an off-shell ac =b entering into the hard process and an emitted parton c. The CS dipoles D (ac,b) relevant for processes of our interest are 3 The role of the spectator is to provide for momentum and colour conservation.

Fig. 4
The NLO Feynman diagram for real-parton radiation in the process of Higgs-boson production in gluon-gluon fusion: the qq channel proportional toP ac,a , the DGLAP kernel for the a → ac splitting. 4 For the processes of the annihilation aā → X at the LO level, such as the Higgs production and the DY process, in each NLO channel ab → cX we must have ac =b in the NLO splitting. In other words, the NLO splitting in the annihilation processes is fully determined by a and b. 5 The above rules are illustrated in Fig. 5 and possible indices are listed in Table 1 for the emission from the incoming line a. 6 Let us first define explicitly the MC distributions (matrix elements) and the CS dipoles representing the initial-state  Table 1 List of indices labelling the CS or MC soft-collinear counterterms for all the NLO channels (except the qq channel in Higgs production) of annihilation processes. Indices (a, b) denote initial partons (channel), while (ac, b) are labelling the CS/MC counter-terms, with a being an emitter and b a spectator real-parton emissions for the Higgs production process in d = 4 + 2ε dimensions: (A) For the g + g → H + g channel a typical/representative distribution of PSMC, summing the emissions from both incoming gluons, is where the g → g splitting function is given bŷ with soft partition functions α α+β and β α+β separating the soft singularity evenly between two incoming emit-ters. Indices (1) and (2) are used to distinguish the above two dipoles. (B) For the g + q → H + q channel we have (with a single soft-collinear pole the soft partition functions are not needed): where the q → g splitting function readŝ (C) Finally, for the g +q → H +q channel, the CS dipole and MC distribution is the same as the previous one for quarks.
The above distributions agree with these used in the POWHEG-method construction of Ref. [16].
For the sake of completeness, let us collect the CS dipoles and MC distributions already known from Refs. [1,2], with the q → q and g → q splittings. They will be needed in the following to define the transition matrix K from the MS to MC factorisation scheme for all PDFs.
(A) For the q +q → Z + g channel, the MC distribution reads wherê 18) and the soft partition function is used again: (3.19) (B) For the q + g → Z + q channel we have (the soft partition function in the MC distribution is not necessary): It should be stressed that all the above MC distributions and CS dipoles are basically in the exclusive (unintegrated) form.
All the above relations between the MC distributions and the exclusive MC/CS counter-terms for any annihilation processes can be summarised in a compact formula as follows: where translation from the indices (ab) to (abc) is unique for a given annihilation process and for a given initial parton splitting, as demonstrated explicitly in Table 1 for the splitting of the initial parton a, see also Fig. 5. Moreover, on the RHS of the above relation only one of D's is nonzero, except for the c = g case (gluonstrahlung), but in this case both D's are equal. Hence, there is in practice one-to-one correspondence (ab) ↔ (abc) for all annihilation processes, to be often exploited in the following section.

Integrated CS dipoles and counter-terms of MC scheme
For the purpose of installing virtual parts (using PDF momentum sum rules) in the MC distributions (soft-collinear counter-terms) and defining the K -matrix for transforming PDFs from the MS to MC scheme, we need to integrate partly all distributions defined in the previous subsection, keeping A z-dependent differential cross section corresponding to the real-emission MC matrix elements can be expressed in the following way: whereˆ MC ab,R (z, ) is the MC real-emission function corresponding to the partly integrated MC distribution of the previous subsection for a given process: a+b → X +c. The inte-gration element dˆ can be expressed in terms of the Sudakov variables as follows: The above expressions are defined in d = 4 + 2 dimensions in order to regularise, in the usual way, the soft and collinear singularities of the real-parton radiation.
Using the exact NLO matrix element, one can similarly write, for each channel ab, a regularised partly integrated NLO cross section for real-parton emission: (3.25) Following Eq. (3.22), one may also define the relation of the integrated MC distribution to the individual integrated soft-collinear counter-terms: whereˆ R are the corresponding integrals dˆ D as in Eq. (3.23). However, contrary to the CS counter-terms, the counter-termsˆ MC of the MC scheme (and theˆ MC radiation functions as well) will also include virtual corrections, calculated using the momentum sum rules, see next subsections for details. Let us calculate all the above objects in more detail for the gg → Hg channel and then, skipping details of analytical integration, for other channels.

gg → Hg channel
A real-emission part of the MC radiation function results from the following integration: 7 7 We employ here and in the following a shorthand notation δ x=y ≡ δ(x − y). (3.27) (3.28) The NLO real correction according to Ref. [17] readŝ (3.29) From the above equations we readily obtain the NLO real coefficient function in the MC scheme: The same expression is obtained in 4 dimensions by means of performing first the MC-dipole subtraction and then integrating the finite result over the phase space: (3.31) A virtual correction to the above MC radiation function gg is calculated from the momentum sum rules: where n f is the number of fermions. The first part in the above virtual correction resulting from integration over the first term in brackets on RHS reads as follows: In order to calculate the second part to the virtual correction in RHS of Eq. (3.32) we need to know first the following MC radiation functions for the g → q transition, e.g. from the process q + g → Z + q: where |M MC qg→Zq | 2 is shown in Eq. (3.20). Using the above result we can cross-check the formula for the gluon-channel MC radiation function of the DY process calculated previously in 4 dimensions in Ref. [2]. For the exact NLO contribution Ref. [18] provideŝ Then the resulting coefficient function for the DY process in the MC scheme reads which agrees with our previous result, given in Ref. [2].
For the sake of completeness, the corresponding coefficient function in the MS factorisation scheme reads and the transition-matrix element transforming part of gluon MS PDF into the quark PDF in the MC scheme is given by (3.38) The above was obtained by comparing the NLO coefficient functions for the DY process in the MS and MC schemes. However, exactly the same result can be obtained alternatively from the difference of the soft-collinear counterterms in these two schemes: where the universal MC-scheme counter-term corresponding to the g → q transition is given bŷ After this brief detour to the DY process, we can now complete the calculation of the virtual correction to the MC radiation function for the gg → Hg channel. Using Eq. (3.34), the second term in RHS of Eq. (3.32) is calculated: The complete result for virtual correction to the gg → Hg MC radiation function, obtained from the momentum sum rule of Eq. (3.32), reads as follows: where T f = n f T R . The complete MC radiation function for gg → Hg process is obtained finally in the following explicit form: Let us also calculate the coefficient function in the MC scheme for the gg → Hg channel. Using the exact NLO virtual correction of Ref. [17]: we obtain the following virtual part of the coefficient function in the MC scheme (with the usual μ 2 = Q 2 assignment): Combining the real and virtual contributions of Eqs. (3.30) and (3.46), the NLO coefficient function for the gg → Hg process in the MC factorisation scheme reads The analogous coefficient function in the MS factorisation scheme is obtained from Eqs. (3.29) and (3.45), after the standard subtraction of the MS soft-collinear counter-terms (with μ 2 = Q 2 ) reads as follows: (3.48) With all the above results at hand we are also ready to determine the element g → g of the transition matrix for transforming PDFs from the MS to MC scheme: The same K MC gg (z) it can also be obtained from the difference of the collinear counter-terms: where the universal MC counter-termˆ MC gg (z, ) corresponding to the g → g transition can be expressed in terms of the MC radiation function of Eq.
is the corresponding MS counter-term.

gq → Hq channel
The channel g + q → H + q is easier because only real correction contributes at NLO. The corresponding MC radiation function can be readily obtained from the integral MC gq (z, ) = The exact NLO correction taken from Ref. [17] readŝ (3.54) Combining the above two functions, the coefficient function for gq → Hq process in the MC factorisation scheme reads (3.55) Exactly the same result can be obtained also from the following integral in 4 dimensions: On the other hand, in the MS factorisation scheme (keeping μ 2 = Q 2 ), from Eq. (3.54) we can obtain (after the standard subtraction) the following coefficient function: (3.57) At this point we are able to define another element of the matrix transforming the MS gluon PDF into the gluon PDF of the MC-scheme (3.58) Alternatively, the same K MC gq (z) can also be obtained as a difference of the soft-collinear counter-terms in the MC and MS schemes: is the corresponding counter-term in the MS scheme.

Revisiting qq → Zg channel
In Ref. [2] the virtual correction to the MC counter-term in the qq channel was calculated from the quark-number conservation sum rule (as minus the integral over z of the real correction). This was justified for the DY process, for which the gluon PDF did not get corrected at NLO from the MS to MC factorisation scheme. Now, since we deal with the complete set of parton-parton transitions, including the transformation/correction of the gluon PDF, we have to rely on the momentum sum rule. For the pertinent channel this amounts tô Using the formula for the MC real-radiation function from Appendix B of Ref. [2]: we can calculate the first part of the above virtual correction as follows: For the second part, using Eq. (3.53), we obtain (3.65) Thus the full virtual correction readŝ (3.66) After combining it with the real correction of Eq. (3.63) we obtain a complete MC radiation function: (3.67) The corresponding NLO correction reads [18] Then, for the coefficient function in the MC factorisation scheme, we obtain (3.69) The above expression differs from the one given in Ref. [2], by a constant term: For completeness, let us also write the corresponding coefficient function in the MS factorisation scheme: and the qq transformation-matrix element to the quark PDF in the MC scheme:

(3.73)
This can also be expressed in a form similar to the corresponding formula for the gg channel, cf. Eq. (3.49):

(3.74)
This is the q → q PDF transition-matrix element from the MS to MC scheme. Similarly as in the previous cases, it can also be obtained from the respective counter-terms: where the universal MC counter-term corresponding to the q → q transition can be related to the MC radiation function is the corresponding MS counter-term.

PDFs in MC scheme
In Ref. [2], were the KrkNLO method was applied to the Drell-Yan process, it was sufficient to transform the MS PDF of quarks and antiquarks. The difference between the MS and MC PDFs for the gluon was an NNLO effect, and hence is beyond the claimed accuracy.
Here, for the Higgs production process, the gluon PDF also has to be transformed to the MC scheme. Having calculated all the necessary ingredients in the previous section, we define this transformation as follows: where K MC gg (z) is given in Eq. (3.49) and K MC gq (z) in Eq. (3.58). However, virtual parts of the transformation matrix in the quark sector now has also changed due to the necessary use of the momentum sum rules. Hence, the entire transformation rule now takes the form where The above formulae can be used for numerical computation of the MC-scheme quark and gluon PDFs from the available parametrisation of the MS PDFs. Alternatively, PDFs in the MC scheme can be fitted directly to DIS and other data, provided the NLO coefficient functions in the MC scheme are known. For DIS they are listed in Appendix A.
We assume that PDFs in the MC scheme satisfy the same momentum sum rule as PDFs in the MS scheme: Inserting in the above formula the expressions for g MC and q MC from Eqs. (4.2), we obtain the following momentum sum rules for the factorisation-scheme transformation matrix elements: The above, of course, results from the momentum sum rules imposed on the MC and MS soft-collinear counterterms, however, it constitutes a useful cross-check of the consistency of the MC scheme.
Looking at the elements of the transition matrix K in Eq. (4.3) one can see that the terms ∼ ln(1 − z) and ∼ ln z are absorbed in the MC-scheme PDFs. As a result the NLO coefficient functions for the DY process and the Higgs-boson production are much simpler than the corresponding ones in the MS scheme, cf. Eqs. • What is the purpose of MC factorisation scheme? It is defined such that the (z)δ(k T ) terms due to emission from initial partons disappear completely from the real NLO corrections in the exclusive/unintegrated form, even before PSMC gets involved. • Why is the above vital in the KrkNLO scheme? Without eliminating such terms it is not possible to include the NLO corrections using simple multiplicative MC weights on top of distributions generated by PSMC. • How to determine elements of the transition matrix K MC ab ? They can be deduced from the difference of soft-collinear counter-terms of the MC and MS scheme or from inspection of the NLO corrections in a few simple processes with initial quarks and gluons in the LO hard process. We have done it both ways.
• Will the same PDFs in the MC scheme eliminate ∼ δ(k T ) terms for all processes? This is a question about the universality of the MC factorisation scheme. For all processes similar to the DY or Higgs-production process, with produced colour-neutral final-state objects, the answer is positive.
In Fig. 6, we present examples of numerical results for the PDFs of quarks and gluon in the MC scheme obtained from PDFs in MS scheme using transformation of Eqs.  Ref. [2] unchanged with respect to MS. As discussed earlier, these types of MC PDFs is sufficient for the Drell-Yan process and it was used in our previous work [2]. Hence, we show them here for comparison.
One can see that the differences between the MC and MS PDFs are noticeable. In particular, the MC quarks are up to 20% smaller at low and moderate x, while they get above the MS distributions at large x. For DY and Higgs production, the latter has consequences only at large rapidities of the bosons. At the same time, we notice that the gluon is larger in the MC scheme at low and moderate x. Hence, the changes in quarks and the gluon have a chance to compensate each other and, indeed, as we checked explicitly, the momentum sum rules (4.4) are numerically satisfied for our MC PDFs.
Other quark flavours, when transformed to the MC scheme, exhibit similar changes to those shown in Fig. 6 for the u and d quarks.
Finally, let us comment briefly on the process-independence (universality) of the MC factorisation scheme and the KrkNLO method. If we treat Eq. (4.2) as a definition of PDFs in the MC scheme, then their universality is just inherited from the MS scheme. The universality of the KrkNLO method is more involved and it would imply that by means of adoption of these PDFs and a careful choice of the exclusive/unintegrated MC distributions for the initial-state splittings, we are able to eliminate from the NLO real corrections all terms proportional to δ(β) f (z) or δ(α) f (z), which means that we can impose the NLO real corrections with the multiplicative MC weights in d = 4 dimensions on top of the PSMC distributions. We are able to state that the above is true for all annihilation process into colour-neutral objects. This can be deduced from analysing the CS counter-terms (which are compatible with the modern PSMCs), where both the emitter and the spectator are in the initial state. They are universal within the class of the above annihilation processes and therefore the KrkNLO method features the same property. The answer to the question whether extending this argument to other processes, with one or more coloured partons in the final state at the LO level, is not trivial and the relevant study is reserved to next dedicated publication. 8 .

NLO cross sections for Higgs production in KrkNLO method
MC weights of the KrkNLO method for the Higgs-boson production in gluon-gluon fusion are very simple, even simpler than those for Drell-Yan process, where they depend on the 8 The analysis in Ref. [1] for the DIS process, albeit limited to the gluonstrahlung NLO subprocess, gives hope for a possible positive answer.
angles of the Z /γ * decay products. For the g + g → H + g subprocess we have whereas for the g + q → H + q channel, the real weight reads For the process with exchanged initial-state partons we have W qg R (α, β) = W gq R (β, α). Virtual+soft-real corrections can be read off from the formulae of the coefficient functions given in Sect. 3. They are just constant terms multiplied by the δ(1 − z) function. In the KrkNLO method they should be included multiplicatively in a parton shower generator for the corresponding process, i.e. the Born-level cross section should be multiplied by the weight The above weights are implemented on top of the CSdipole-based PSMC algorithm of Herwig 7 [19,20] in the socalled "power shower" mode [21] which allows for complete coverage of the phase space in momentum and flavour space, without empty regions. 9 In a way this is analogous to the one described in Ref. [2] for the Drell-Yan process. In Ref. [2] we provide a detailed discussion of the PSMC algorithm for the case of gluonstrahlung where two CS dipoles contribute (cf. Sects. 3.1 and 3.3). Then we prove that applying to such a PSMC an appropriate MC weight according to the KrkNLO method indeed reproduces the NLO differential cross section (cf. Sect. 3.4). All this can be adapted to the current case of the Higgs-boson production, replacing only incoming quarks by gluons. We therefore do not repeat such a discussion herethe interested reader is recommended to check the above paper.
For the numerical evaluation of the cross sections at the LHC for the proton-proton collision energy of and the G μ -scheme [23] for the electroweak sector. To compute the hadronic cross section we also use the MSTW2008 LO set of parton distribution functions [24], and take the renormalisation and the factorisation scales to be μ 2 R = μ 2 F = M 2 H , where M H is the Higgs-boson mass. We also set the Higgs boson to be stable for simplicity.
In Table 2 we show the results for the total cross sections for the Higgs-boson production in gluon-gluon fusion obtained with KrkNLO and MC@NLO. The two methods are matched to the dipole parton shower implemented in Herwig 7 [19,20].
We see that the two methods give slightly different (∼ 3.5%) total cross sections, which come from formally higherorder terms, i.e. beyond the NLO approximation. The relevant distributions and detailed comparisons with MC@NLO, POWHEG and the NNLO calculations from the HNNLO program [25,26] are presented in another publication [27].

Summary and outlook
In this work, we have presented all the ingredients of the KrkNLO method needed for its implementation for the Higgs-boson production process in gluon-gluon fusion. In particular, the complete definitions of PDFs in the MC scheme, together with their numerical distributions, have been provided. Hence, PDFs in the MC FS can be fitted directly to experimental DIS and DY data. We have also presented the first result for the total cross section for the Higgs production. More distributions, comparisons with MC@NLO, POWHEG and the NNLO calculations are presented in a separate paper [27]. A dedicated study of the process-independence (universality) of the KrkNLO method and the MC factorisation scheme is also reserved for the future work.
The current state of NNLO+PS [28][29][30][31][32][33][34] represents a clear progress in matching fixed-order QCD calculations with PSMCs, however they are still limited to certain classes of observables. The other natural extension for KrkNLO is NNLO+NLOPS, where NLOPS is a PSMC that implements the NLO evolution kernels in the fully exclusive form and thus provides the full set of the soft-collinear counter-terms for the hard process. These coefficient functions can be used in fitting the MC PDFs to experimental DIS data.