# NLO calculations matched with the Nagy-Soper parton shower

###### Abstract

An MCNLO-like matching of NLO QCD calculations with the Nagy-Soper parton shower is briefly summarised. Uncertainties and ambiguities of the matching scheme are shortly discussed. A few results for the production process at the LHC with TeV are also shown. All results have been obtained using the Nagy-Soper parton shower implementation in the Deductor program together with the Helac-NLO framework.

NLO calculations matched with the Nagy-Soper parton shower

Michal Czakon, Heribertus B. Hartanto, Manfred Kraus, Malgorzata Worek^{†}^{†}thanks: Speaker. ^{†}^{†}thanks: Preprint number: TTK-16-04

Institute for Theoretical Particle Physics and Cosmology

RWTH Aachen University

D-52056 Aachen, Germany

E-mails:

\abstract@cs

## 1 Introduction

After an almost two year shutdown, the Large Hadron Collider (LHC) started delivering physics data in 2015. This marked the begin of Run 2 at the LHC and opened the path to an even deeper understanding of the Standard Model (SM) physics and hopefully to new discoveries. Protons are now collided at the LHC at the record breaking energy of TeV, almost double the collision energy of Run 1 ( TeV). On the of December 2015, the ATLAS and CMS experiments both reported a number of results using TeV proton collision data. Even though the amount of data on which analyses are based is still limited, experimentalists have succeeded in producing numerous results for the SM physics, including new results for the Higgs boson and the top quark. Additionally, many searches for physics beyond the SM, yielding many improved limits beyond Run 1, have been presented. In 2016 we shall wait for more (and updated) results. To describe experimental data various theoretical tools are needed, among others, Monte Carlo (MC) generators that are widely used in all experimental analyses. Such programs need to be improved as data become more precise. There are various ways to reach better accuracy of theoretical predictions. In particular, one such possibility is an increase of the order in the fixed order perturbation expansion. Next-to-leading order (NLO) calculations can now be performed in a fully automatic manner. We mention here only the Helac-NLO framework [1], that will be used in our studies. The software has recently been used to calculate NLO QCD corrections to a process with five objects in the final state, i.e. [2]. Adding one more order complicates the picture tremendously. For next-to-next-to-leading order (NNLO) calculations, such automatic frameworks are not yet publicly available. In addition, a list of processes that have been completed at that order comprises processes with two objects in the final state only. The most prominent example being the NNLO calculation for the production process at hadron colliders, where two coloured and massive final state fermions are present [3, 4]. Finally, an important milestone has been recently reached when a first calculation at the next-to-next-to-next-to-leading order () in perturbative QCD for the Higgs production in gluon fusion in the large top-mass limit has been completed [5]. Since for each order one needs to include more Feynman diagrams both with one particle more in the final state and with one loop more in the intermediate state, a formal order-by-order perturbative calculation is limited to a few orders only. Alternatively, to improve the theoretical accuracy one can use parton shower (PS) programs. Parton shower algorithms approximate higher-order corrections by including the leading soft and collinear contributions to all orders. These programs are not able to correctly estimate the radiation of hard jets. They can, however, generate arbitrarily many jets in the final states. In order to improve the simulation of hard jets production in the parton shower, approaches were developed to match parton showers with fixed order NLO calculations [6, 7] and to merge matched calculations for different jet multiplicities both at LO and NLO [8, 9, 14, 10, 13, 11, 12, 15]. Finally, those methods have already started to be applied to the NNLO calculations. Formally, however, in parton showers that are currently being use by experimental collaborations the resummation of soft and collinear contributions is only certified to Leading Logarithmic (LL) accuracy. In practice, many of the NLL contributions are already included, either by the angular ordering of the parton shower or by an optimal scale choice for . Thus, yet another progress on the front of theoretical predictions for the LHC can be accomplished by improving the parton shower itself. This can be done through an inclusion of color and spin correlations in subsequent parton shower emissions. An impact of the colour suppressed terms in parton shower simulations, for the final state parton shower, have already been studied in [16]. First results for the production at LEP energy have shown that standard LEP-observables, e.g. event shapes and jet rates, are only slightly affected. However, for tailored observables, i.e. observables that are sensitive to soft (wide angle) splitting, deviations as large as have been observed. The first step in the direction of the complete algorithm, which consistently takes into account both spin and color correlations in the parton emissions for the collision of hadrons, has been done by Nagy and Soper [17, 18]. A large part of their idea has been implemented in the Deductor program [19]. Even though the exact colour and spin correlations are not yet included in the software, it can already provide results with an extension of the leading colour approximation, which has been implemented instead.

In this proceedings we shall briefly summarise the concept behind the Nagy-Soper parton shower. In the next step we shall outline a matching procedure between the Nagy-Soper parton shower and a fixed order calculation at NLO [20]. To this end, we use the Mc@Nlo approach [6]. When matching NLO calculations with parton shower programs, one needs to avoid the double counting of emissions, which on the one hand can be generated with a jet from the matrix element and on the other hand can appear as a jet from the parton shower. The Mc@Nlo method removes double counting contributions by expanding the parton shower to first order in and compensating for the terms, which are already present at fixed order. Finally, a few results for the production process, which involves non-trivial colour exchange, massive partons and requires cuts already at the Born level, will be shown.

## 2 Nagy-Soper parton shower

In case of a generic process the expectation value for a completely inclusive observable in the Nagy-Soper formalism can be written as follows

(2.0) |

where the sum runs over all final state multiplicities. Here, is the sum of all -particle phase space measures for different flavour sequences . The factor is necessary to account for identical contributions. The parton density functions evaluated at the momentum fraction and factorization scale are denoted by . The factor in the denominator comes from averaging over initial state spins and is the colour factor. The quantity

(2.0) |

is the quantum density and a matrix element can be viewed as a vector in colour spin space and can be resolved into basis vectors and with complex expansion coefficients

(2.0) |

The propagation of the quantum density matrix from some initial shower time, , to some final shower time, , is described by the evolution equation. The initial shower time corresponds to the hard interaction, whereas the final shower time characterises the physical scale at which parton emissions cannot be described perturbatively. The perturbative evolution is described by a unitary operator that obeys

(2.0) |

Here, describes resolved emission and the unresolved/virtual one. The latter can be further decomposed into a color diagonal, , and a color off-diagonal part, , as follows

(2.0) |

Eq. (2) can be solved as

(2.0) |

where is the Sudakov form factor (a number) defined as

(2.0) |

In case of a non-trivial colour evolution, the exponentiation of a non-diagonal matrix is very difficult. Thus, only the colour diagonal part, , is exponentiated, whereas the off-diagonal part, , is treated perturbatively in the same way as . The expectation value of the observable , including shower effects, is provided via

(2.0) |

The evolution in the shower time is always ordered in some chosen kinematic variable to correctly resum leading logarithms of infrared sensitive quantities. In the Nagy-Soper parton shower the following variable has been proposed

(2.0) |

where is the emitter momentum after emission, the emitted parton momentum, the emitter momentum before emission () and is the total final state momentum. The minus/plus sign between and in Eq. (2) applies to an initial/final state splitting.

## 3 Matching Inclusive Processes

There exist several schemes for matching NLO calculations with parton showers, the most popular being Powheg and Mc@Nlo. In order to benefit from the recently implemented subtraction scheme based on the Nagy-Soper parton shower splitting kernels [21], which is implemented in the Helac-Dipoles framework [22], we chose the Mc@Nlo formalism. For a generic process at NLO, one can write the quantum density matrix in a perturbative expansion in , according to

(3.0) |

where the leading order contribution is counted as order in the strong coupling . In addition, and are tree level matrix elements, whereas is the one-loop amplitude. Based on this quantum density matrix, the expectation value of the observable including shower effects suffers from double counting. This can be seen when the iterative solution to the evolution equation is expanded to . Indeed, we obtain

(3.0) |

and contains the first emission contributions twice, once from the real emission quantum density , and once from the parton shower approximation . The problem of double counting can be overcome if the density matrix is modified as follows

(3.0) |

Considering and expanding the evolution equation up to shows that the undesired parton shower contributions are cancelled. Thus, for an infrared safe observable , we can write

(3.0) |

where and can be either or . The parton shower splitting kernels are used to provide subtraction terms for the infrared singularities

(3.0) |

(3.0) |

where the sum runs over all external legs and is the total splitting kernel for a given external leg . The parameter is the shower time and represents the ordering of the emissions. The phase space integration of the additional parton is denoted by . The decomposition of the integrated into two operators and is arbitrary. However, we choose to match the divergencies of the virtual amplitude. In addition, the limit has been taken. This is a source of a mismatch between the fixed order and the shower calculation, however, it is numerically small due to the exponential damping by the Sudakov form factor. The total cross section, which includes parton shower evolution amounts to

(3.0) |

where the following shorthands has been defined

(3.0) | |||

(3.0) |

The whole procedure is divided into two steps. First the samples according to Eq. (3) and (3) are generated and afterwards is applied.

## 4 Matching Exclusive Processes

For processes with massless partons already at leading order, the matching prescription as described in the previous section must be slightly modified by the inclusion of generation cuts. The following replacements

(4.0) | ||||

(4.0) |

can be made, where is a jet function applied during the generation of events, on the momenta and flavours of . However, when the parton shower is applied to these ensembles we can easily see that double counting is still there. Indeed, after expanding the evolution operator we obtain

(4.0) |

where the contribution of has been cancelled by the linear expansion of the Sudakov form factor. This mismatch is cured by enforcing the subtraction terms to fulfill , i.e. by modifying the real subtracted cross section according to

(4.0) |

where we make use of the inverse momentum mapping and obtain

(4.0) |

This modification allows us to introduce restrictions on the functional form of . Expanding the shower evolution yields

(4.0) |

The double counting is removed from Eq. (4), when the following condition is true

(4.0) |

We have used here the fact that . This can be achieved when

(4.0) |

The interpretation of the last equation is straightforward. Generation cuts () need to be more inclusive than the cuts on the final observable (). Let us note here, that the Mc@Nlo matching formalism in combination with the Nagy-Soper parton shower introduces the intrinsic uncertainties. Let us name just the most important ones here. The Nagy-Soper parton shower treats bottom and charm quarks as massive. On the other hand in the NLO calculation we treat them as massless. Thus during the matching procedure masses for the relevant quarks are introduced by the on-shell projection. In addition, PDFs are evolved differently in the NLO calculation and in the shower. In the former case, NLO PDFs are used, whereas in the latter PDFs are evolved using Nagy-Soper splitting kernels. The evolution is, however, of higher order and the NLO accuracy is maintained if the evolutions share a common point e.g. at the low scale. Finally, the choice of the initial parton shower time in the parton evolution is arbitrary. We only require that NLO prediction is recovered for hard emissions. Thus, different choices of that fulfil this condition are possible. In consequence, the quantity should be varied to study the uncertainty of NLO+PS matching.

## 5 Implementation

The matching scheme has been implemented within the Helac-Nlo framework [1]. We have used it in conjunction with the Nagy-Soper parton shower present in Deductor version 1.0.0 [19]. Since Deductor uses spin averaged splitting functions, we only provide unpolarized event samples for showering. We supply only leading colour events even though Deductor works with the approximation, which comprises a full color description for collinear and soft-collinear limits and the one for pure soft limits. The Helac-Nlo Monte Carlo program generates an event sample ready for showering with Deductor. We produce events subprocess by subprocess. First, we use Helac-1Loop [23] to obtain a set of unweighted leading order events with the virtual contributions. This set of weighted events is subsequently reweighted using Helac-Dipoles in order to include the parton shower virtual operator, which corresponds to taking into account the integrated subtraction terms. The real radiation events are generated separately with Helac-Dipoles. The program has been extended to provide unweighted events with positive and negative weights, which is possible, because both the real radiation and the respective subtraction weights correspond to the same phase space point. For each accepted event, we choose the most probable diagonal colour flow configuration. The generated events are stored in a Les Houches file and transferred to Deductor, which requires an on-shell projection for charm and bottom quarks and a translation of the color flow language from the LHE file to internal representation in Deductor in terms of color strings.

## 6 Top quark pair production in association with a jet at NLO+PS

We present now the results for the production at the LHC with TeV. We use stable top quarks and put GeV. The charm and bottom quarks are considered to be massless at fixed order. Results are obtained using the Mstw2008nlo PDF set [24] with five active flavours and the corresponding two-loop running of the strong coupling. The renormalization and factorisation scales are set to , and the starting shower time to

(6.0) |

where and are external momenta, is the total final state momentum and for the central prediction. The anti- jet algorithm with the separation parameter is used to cluster partons with pseudorapidity . The resulting jets are ordered according to their . We also require the tagged jets to have transverse momentum of GeV and rapidity in the range of . Our analysis is restricted to the perturbative parton shower evolution. Decays of unstable particles, hadronization and multiple interactions are not taken into account. The parton shower treats the charm and bottom as massive particles, thus, we use GeV and GeV. Additionally, the Mstw2008nlo PDF set at GeV is provided as the starting point for the evolution in Deductor. We also use the corresponding two-loop running of , and restrict the parton shower to the leading colour approximation. Results presented in the following are accurate up to .

In Fig. 1 we present the transverse momentum of the system together with inclusive jet cross sections.

The variation bands for and have been obtained using the following sets of three parameter values: and , respectively. The lower panels display corresponding relative deviations from the central value, separately for and . We start the discussion with the transverse momentum of the system, that is given in the left panel of Fig. 1. At LO the transverse momentum of this system is zero due to momentum conservation. When NLO contributions are included, this observable diverges as the transverse momentum of the entire system goes to zero. Thus, it can only be reliably described by the fixed order calculation in the high region. The situation is changed when the parton shower is included. In that case, the low behaviour is altered strongly by the Sudakov form factor as can be seen in Fig. 1. We observe a moderate dependence on in the low region, up to a factor of at the lower end of the spectrum. This dependence decreases down to just a few percent around GeV, whereas for moderate values of it is at the level of . The behaviour is reversed for the renormalization and factorisation scale dependence. Visible deviations from the central value occur once the matrix element dominates and grow substantially up to almost at the end of the spectrum.

The inclusive jet cross sections are given in the right panel of Fig. 1. The NLO cross section with exactly one jet, which is given in the first bin, is rather insensitive to , i.e. at the level. The dependence is slightly larger in the second bin, where the two jet cross section, correct only at the LO level, is stored. Starting from the third bin, cross sections are described via the shower evolution alone, thus, we observe fairly large variations for both parameters, and . For example, the scale dependence for the cross section with five jets is found to be around .

## 7 Summary

To summarise, we have outlined a NLO matching scheme for the Nagy-Soper parton shower. We based our construction on the original Mc@Nlo approach. We have presented the general formulation and a few results for top-quark pair production in association with a jet at the LHC. All results given here have been obtained using an implementation within the framework of the public codes Helac-Nlo and Deductor.

## Acknowledgement

This research was supported in part by the German Research Foundation (DFG) under grant WO 1900/1-1 Signals and Backgrounds Beyond Leading Order. Phenomenological studies for the LHC.

## References

- [1] G. Bevilacqua, M. Czakon, M. V. Garzelli, A. van Hameren, A. Kardos, C. G. Papadopoulos, R. Pittau and M. Worek, Comput. Phys. Commun. 184 (2013) 986 [arXiv:1110.1499 [hep-ph]].
- [2] G. Bevilacqua, H. B. Hartanto, M. Kraus and M. Worek, arXiv:1509.09242 [hep-ph].
- [3] P. Baernreuther, M. Czakon and A. Mitov, Phys. Rev. Lett. 109 (2012) 132001 [arXiv:1204.5201 [hep-ph]].
- [4] M. Czakon, P. Fiedler and A. Mitov, Phys. Rev. Lett. 110 (2013) 252004 [arXiv:1303.6254 [hep-ph]].
- [5] C. Anastasiou, C. Duhr, F. Dulat, F. Herzog and B. Mistlberger, Phys. Rev. Lett. 114 (2015) 212001 [arXiv:1503.06056 [hep-ph]].
- [6] S. Frixione and B. R. Webber, JHEP 0206 (2002) 029 [hep-ph/0204244].
- [7] P. Nason, JHEP 0411 (2004) 040 [hep-ph/0409146].
- [8] M. L. Mangano, M. Moretti, F. Piccinini and M. Treccani, JHEP 0701 (2007) 013 [hep-ph/0611129].
- [9] J. Alwall et al., Eur. Phys. J. C 53 (2008) 473 [arXiv:0706.2569 [hep-ph]].
- [10] K. Hamilton, P. Nason and G. Zanderighi, JHEP 1210 (2012) 155 [arXiv:1206.3572 [hep-ph]].
- [11] S. Hoeche, F. Krauss, M. Schonherr and F. Siegert, JHEP 1304, 027 (2013) doi:10.1007/JHEP04(2013)027 [arXiv:1207.5030 [hep-ph]].
- [12] T. Gehrmann, S. Hoche, F. Krauss, M. Schonherr and F. Siegert, JHEP 1301, 144 (2013) doi:10.1007/JHEP01(2013)144 [arXiv:1207.5031 [hep-ph]].
- [13] R. Frederix and S. Frixione, JHEP 1212 (2012) 061 [arXiv:1209.6215 [hep-ph]].
- [14] L. Lonnblad and S. Prestel, JHEP 1303, 166 (2013) doi:10.1007/JHEP03(2013)166 [arXiv:1211.7278 [hep-ph]].
- [15] S. Hoeche, F. Krauss, P. Maierhoefer, S. Pozzorini, M. Schonherr and F. Siegert, Phys. Lett. B 748 (2015) 74 [arXiv:1402.6293 [hep-ph]].
- [16] S. Platzer and M. Sjodahl, JHEP 1207 (2012) 042 [arXiv:1201.0260 [hep-ph]].
- [17] Z. Nagy and D. E. Soper, JHEP 0709 (2007) 114 [arXiv:0706.0017 [hep-ph]].
- [18] Z. Nagy and D. E. Soper, JHEP 1206 (2012) 044 [arXiv:1202.4496 [hep-ph]].
- [19] Z. Nagy and D. E. Soper, JHEP 1406 (2014) 097 [arXiv:1401.6364 [hep-ph]].
- [20] M. Czakon, H. B. Hartanto, M. Kraus and M. Worek, JHEP 1506 (2015) 033 [arXiv:1502.00925 [hep-ph]].
- [21] G. Bevilacqua, M. Czakon, M. Kubocz and M. Worek, JHEP 1310 (2013) 204 [arXiv:1308.5605 [hep-ph]].
- [22] M. Czakon, C. G. Papadopoulos and M. Worek, JHEP 0908 (2009) 085 [arXiv:0905.0883 [hep-ph]].
- [23] A. van Hameren, C. G. Papadopoulos and R. Pittau, JHEP 0909 (2009) 106 [arXiv:0903.4665 [hep-ph]].
- [24] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63 (2009) 189 [arXiv:0901.0002 [hep-ph]].