An automated subtraction of NLO EW infrared divergences

The European Physical Journal C, Feb 2018

In this paper a generalisation of the Catani–Seymour dipole subtraction method to next-to-leading order electroweak calculations is presented. All singularities due to photon and gluon radiation off both massless and massive partons in the presence of both massless and massive spectators are accounted for. Particular attention is paid to the simultaneous subtraction of singularities of both QCD and electroweak origin which are present in the next-to-leading order corrections to processes with more than one perturbative order contributing at Born level. Similarly, embedding non-dipole-like photon splittings in the dipole subtraction scheme discussed. The implementation of the formulated subtraction scheme in the framework of the Sherpa Monte-Carlo event generator, including the restriction of the dipole phase space through the \(\alpha \)-parameters and expanding its existing subtraction for NLO QCD calculations, is detailed and numerous internal consistency checks validating the obtained results are presented.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-018-5600-z.pdf

An automated subtraction of NLO EW infrared divergences

Eur. Phys. J. C An automated subtraction of NLO EW infrared divergences Marek Schönherr 0 0 Theoretical Physics Department , CERN, 1211 Geneva 23 , Switzerland In this paper a generalisation of the CataniSeymour dipole subtraction method to next-to-leading order electroweak calculations is presented. All singularities due to photon and gluon radiation off both massless and massive partons in the presence of both massless and massive spectators are accounted for. Particular attention is paid to the simultaneous subtraction of singularities of both QCD and electroweak origin which are present in the next-to-leading order corrections to processes with more than one perturbative order contributing at Born level. Similarly, embedding non-dipole-like photon splittings in the dipole subtraction scheme discussed. The implementation of the formulated subtraction scheme in the framework of the SHERPA Monte-Carlo event generator, including the restriction of the dipole phase space through the α-parameters and expanding its existing subtraction for NLO QCD calculations, is detailed and numerous internal consistency checks validating the obtained results are presented. 1 Introduction . . . . . . . . . . . . . . . . . . . . . 1 2 Catani-Seymour subtraction at NLO EW . . . . . . 2 2.1 Differential subtraction terms . . . . . . . . . . 3 2.2 Integrated subtraction terms . . . . . . . . . . . 4 3 Implementation . . . . . . . . . . . . . . . . . . . . 5 3.1 Identification of dipoles . . . . . . . . . . . . . 5 3.2 External photons . . . . . . . . . . . . . . . . 7 3.3 Process mappings . . . . . . . . . . . . . . . . 7 3.4 Fiducial phase space definition . . . . . . . . . 8 3.5 Flavour scheme conversion . . . . . . . . . . . 9 4 Checks of the implementation . . . . . . . . . . . . 9 5 Conclusions . . . . . . . . . . . . . . . . . . . . . 14 A Differential splitting functions . . . . . . . . . . . . 14 A.1 Final-final dipoles . . . . . . . . . . . . . . . 14 A.2 Final-initial dipoles . . . . . . . . . . . . . . . 15 A.3 Initial-final dipoles . . . . . . . . . . . . . . . 16 A.4 Initial-initial dipoles . . . . . . . . . . . . . . 17 B Integrated splitting functions . . . . . . . . . . . . . 17 B.1 Terms of the I perator . . . . . . . . . . . . . . 17 B.2 Terms of the K and P operators . . . . . . . . . 18 C Precise forms of AiIk and AaKb . . . . . . . . . . . . 20 C.1 The AiIk term . . . . . . . . . . . . . . . . . . 20 C.2 The AaKb term . . . . . . . . . . . . . . . . . . 21 References . . . . . . . . . . . . . . . . . . . . . . . . 22 - Contents 1 Introduction As Run-I of the LHC has been successfully completed, culminating in the celebrated experimental confirmation of the existence of the Higgs boson, Run-II proceeds its data-taking at the unprecedented centre-of-mass energy of 13 TeV. As the much anticipated discovery of signals of beyond-theStandard-Model physics is still lacking, precision tests scrutinising the Standard Model are of prime importance, now and in the foreseeable future. At the same time, new physics searches are looking for increasingly small signals demanding more precise estimates of the Standard Model backgrounds. This expansion of sensitivity of both precision measurements and new physics searches in the multi-TeV region demand an immense improvement of theoretical predictions. This precision can be achieved by the inclusion of next-to and next-to-next-to-leading order (NLO and NNLO) corrections in the strong coupling and next-to-leading order electroweak (EW) corrections. Here it should be noted that both NNLO QCD and NLO EW corrections are expected to be of a similar magnitude for inclusive observables as numerically αs2 ≈ α. On selected differential distributions, however, electroweak corrections can grow much larger. They are dominated by photon emissions in the distributions of final state leptons, for example. In invariant mass spectra of lepton pairs below a resonance, for example, O(1) corrections can be present, in which case a proper resummation should be included [ 1 ]. Similarly, looking at the (multi-)TeV regime, the NLO EW corrections quickly grow considerably, reducing cross sections by a few tens of percent, due to the emergence of large electroweak Sudakov corrections arising as the scattering energies Q2 m2W [ 2–17 ]. In this regime they are larger than even the NLO QCD corrections in many cases and their omission becomes the dominant uncertainty in experimental studies and searches. To this end, benefiting from the well-established techniques developed for the automation of NLO QCD corrections many NLO EW corrections have been calculated recently [ 18–38 ]. To fully automate these computations at NLO EW accuracy in a Monte-Carlo framework all infrared divergences need to be regulated, where various incarnations of subtraction methods have proven to be the methods of choice for practical implementations [ 39–43 ]. Similar subtractions have also been published for NLO EW calculations using a dipole picture [ 1,44–46 ]. Contrary to the QCD case, only the implementation of [46], restricted to photon emissions of fermions, is publicly available though. Besides the generalisation to all divergent splittings at O(α), including photon splittings and photon emissions off massive scalars and vector bosons, this publication addresses the issue of automatically detecting simultaneously occurring QCD and QED singularities and subtracting them consistently. These occur as soon as the Born process is defined at multiple orders O(αsnα N −n). In this case the NLO EW correction to the O(αsnα N −n) process, being of O(αsnα N −n+1), is at the same time the NLO QCD correction to the O(αsn−1α N −n+1) process and will in general exhibit the corresponding infrared singularities. Further, matters of the organisation of the contributing partonic processes and their mapping to reduce the computational complexity along with the provision of infrared safe phase space cuts are discussed. The algorithm is implemented in the AMEGIC [ 47 ] matrix element generator which is part of the SHERPA [ 48 ] Monte-Carlo event generator framework. It bases on the automated subtraction of massless NLO QCD divergences therein [ 49,50 ]. The implementation presented in this publication has, in various preliminary forms, already been used to calculate electroweak corrections to a multitude of important signal and background processes [ 19,21,25,31,32,34,37,51–53 ], highlighting its versatility. The present paper is structured as follows: first, in Sect. 2 the Catani–Seymour dipole subtraction method is reviewed and the general modifications to its differential and integrated subtraction terms are discussed. Section 3 then details its automation in SHERPA’s matrix element generator AMEGIC, highlighting the necessary changes and improvements with respect to [ 49,50 ]. This section also discusses various options implemented for the incorporation of photon splittings, general infrared safe fiducial phase space definitions and flavour scheme conversions. Essential cross checks validating the presented implementation are then presented in Sect. 4 before concluding in Sect. 5. Explicit formulae for all differential and integrated dipoles are given in Appendix A–C. 2 Catani–Seymour subtraction at NLO EW In order to be applicable to NLO EW calculations the wellknown Catani–Seymour dipole subtraction [ 40,41 ] needs to be recast in a suitable form. To highlight the changes from the original formulation for NLO QCD calculations the complete structure of the formalism is reviewed. This subtraction formalism starts from the the expectation value of any infrared safe observable O described at NLO accuracy through O NLO = d (m4) B( (m4)) O( (m4)) + d (md) V( (md)) + C( (md)) O( (md)) + d (d) m+1 R( (md+)1) O( (md+)1) =0 . (2.1) Therein, the Born term B consist of the squared matrix element |Mm |2 = m s1, . . . , sm |s1, . . . , sm m with helicity states sn and further includes all parton densities, parton fluxes, symmetry and averaging factors. The virtual and real corrections, V and R, as well as the collinear counter term, C, are defined analogously. When regulating their respective divergences through dimensional regularisation, they have to be evaluated consistently in d = 4 − 2 dimensions for all singularities to cancel. Only after their summation can the limit → 0 be taken. d (n4) and d (nd) are the four and d dimensional phase space element. As V and C on the one hand side and R on the other are defined on phase spaces of different parton multiplicity, Eq. (2.1) cannot be used for numerical evaluation straight forwardly and it is rewritten as O NLO = d (m4) B( (m4)) O( (m4)) + + d (m4) V( (md)) + C( (md)) + d (1d) D( (md) · (1d)) d (4) m+1 R( (m4+)1) O( (m4+)1) −D( (m4) · (14)) O( (m4)) , =0 (1) FF (3) IF a pa aj pj ıj pk j pi pj k i j pk k (2) FI (4) II a pa a b ıj pi pj i j aj pa pb pj j collinear counterterm on the other side are cancelled separately by the integral of the subtraction term over the oneparticle phase space, rendering the integrals of the second line finite in four dimensions as well. Section 2.1 now describes the construction of the differential subtraction term, D, used to subtract all divergences from the real emission correction, while Sect. 2.2 presents the integrated subtraction terms, ID = d 1 D, that subtracts all divergences from the virtual corrections and the collinear counterterm in its explicit Laurent expansion after being analytically integrated over the factorised one-particle phase space. All infrared divergences that occur at NLO EW are of QED origin. No subtractions of potentially large, but finite, corrections involving the emissions of real and virtual massive electroweak gauge bosons will be considered. The practical implementation described in Sect. 3 follows the general lines of [ 49,50 ]. 2.1 Differential subtraction terms To describe all singular limits of a given real emission matrix element R, related to the real emission term of Eq. (2.1) in a similar fashion as the Born term, it is decomposed into a sum over dipoles D [ 40,41 ] as 2 |Mn+1| = R → D = i, j k=i, j a D j,k + Di j,k + Daj ,b. i, j a a Di j (2.3) Therein, i is the emitter in the final state, j is the emittee, k is the spectator in the final state, a and b are the initial state partons. Each dipole thus encodes the singularity structure caused by the emission of j in the presence of the charge of the spectator. While the divergence associated with the soft emission of j off the dipole ıj − k˜ is partially fractioned into a piece associated with a splitting ıj → i + j in the presence of spectator k and a piece where i and k swap their roles, the divergence associated with the collinear emission of j off ıj is recovered through charge conservation once all dipoles having ıj as emitter are summed over. All four dipole types are depicted in Fig. 1. The individual dipoles take the form [ 40,41,54,55 ] 1 Qı2j k˜ m . . . , ıj , . . . , k˜, Di j,k = − ( pi + p j )2 − mı2j Therein, the charge-correlator is defined as [ 1,44,45 ] 2 Qıj k˜ = ⎧ Qıj Qk˜θıj θk˜ ıj = γ ⎨ Qı2j ⎩ κıj k˜ ıj = γ and κıj k˜ = −1 ∀ıj = γ . k˜=ıj (2.5) The Qıj and Qk˜ are the charges of the emitter and the spectator and their θıj/k˜ are 1(−1), if they are in the final (initial) state. Of course, Qk˜ = Qk and θk˜ = θk . In the case of photon splittings no soft divergence is present. Thus, these splittings have no dipole character. To include them in the dipole-formalism nonetheless and to distribute the recoil in splittings away from the collinear limit, spectator partons need to be assigned. As their only role is to absorb transverse momentum of the splitting process, any other particle may be considered as spectator. Each thus assigned recoil partner may be assigned a weight κıj k˜ , only their sum is constraint by Eq. (2.5) in order to add up to the correct collinear limit. Various options to assign recoil partners are implemented, they are detailed in Sect. 3.2. Since the QED charges are real numbers, the chargecorrelator simply multiplies the matrix element and only leaves the spin-correlators Vi j,k , Viaj , Vaj,k and Va,b as inserj tions in the spin-correlated underlying Born matrix elements. The spin-correlators directly correspond to their QCD counterparts and are detailed in Appendix A. It also defines the initial state momentum rescaling parameters xi j,a , xa j,k and xa j,b, as well as the splitting variables yi j,k , u j and v j . The {αdip} = {α FF, α FI, α IF, α II} parameters serve to restrict the phase space where the individual dipole terms are non-zero and therefore need to be evaluated [ 54,55 ]. They are constructed such that for every αıj k˜ > 0 the singularity is fully subtracted. The introduction of a parameter κ in dipoles where a final state photon splits into a massive fermion pair in the presence of a final state spectator similarly allows a redistribution of finite terms, cf. Appendix A. 2.2 Integrated subtraction terms By the above construction, the integral of the subtraction terms D over the one-particle phase space possesses all the necessary poles to render the second line in Eq. (2.2) finite as → 0. This section now summarises the formulae and findings of [ 40,41 ], translated to the QED case, and discusses their important features. The integrated subtraction terms, together with the collinear counterterm, are commonly reorganised into I, K and P operators through the following identification a,b = a,b dηa dηb + + a b 2 + Paa (xa ; μF ) dxa Ba b( (m4)) · Kaa (xa ; {αdip}) dxb Bab ( (m4)) · Kbb (xb; {αdip}) 2 + Pbb (xb; μF ) O( (m4)). (2.6) The I operator contains the necessary infrared poles to cancel all divergences of the virtual correction such that, after summing both, the limit → 0 can be taken and the integral can be evaluated in four dimensions. The K and P operators are infrared finite by construction and, thus, evaluation in four dimension is unproblematic as well. Please note, that due to the fact that, contrary to the colour correlator in QCD, the charge correlator in QED is a simple real number and thus the convolution of the Born matrix element with the respective insertion operators in QCD becomes a trivial product of the Born matrix element and the operators in QED. The spincorrelation that was still present in the real subtraction terms has been integrated out in full analogy to the QCD case. In the following, the structure of all three operators in the QED case is discussed. The I operator. The I operator contains all flavour-diagonal endpoint contributions and cancels all divergences present in the one-loop matrix elements. It takes the general form α I( , μ2; κ, {αdip}) = − 2π (4π ) (1 − ) × i k=i Iik ( , μ2; κ, {αdip}) (2.7) and contains single and, in case of massless emitters, double poles in . Due to the presence of such poles a dependence on the regularisation scale μ2 enters. It is commonly identified with the renormalisation scale μ2R. The I operator is further dependent on the choice of the {αdip} and κ parameters that effect non-singular terms only. It is further decomposed into dipoles [ 41 ] Iik ( , μ2; κ, {αdip}) = Qik Vik ( , μ2; κ) + i ( , μ2) + γi 1 + ln 2 + Ki + AiIk ({αdip}) + O( ) , μ2 sik (2.8) wherein i labels the emitter and k the spectator. The charge insertion operator, which is a trivial real number in the QED case, is defined in Eq. (2.5). The full crossing invariance of the QCD I operator may be broken in the QED case in the presence of photon splittings as their recoil partner assignment is arbitrary and may involve information on initial or final state particles. Different possible choices are discussed in Sect. 3.2, some of which may break this crossing invariance. The divergences of the I operator are encoded in the functions Vik and i . While the former contains all soft-(quasi)collinear divergences the latter contains the pure (quasi-) collinear ones. They do not only differentiate whether the emitter is a photon or not, but also between different spins of the emitter. Their precise form as well as the the flavour constants γi and Ki are given in Appendix B. AiIk encodes the dependence on the phase space restriction of the individual dipoles {αdip}. Finite terms originating in dipoles involving initial state legs, however, can be pushed into the K operator. Thus, AiIk by convention only depends on α FF. Its precise form is given in Appendix C. The K and P operators. The K and P operators collect all pieces of the integrated dipole terms that are not collected in the I operator and combines them with the collinear counterterms C to give a finite result as → 0. By construction they contain only remainders of splittings where either the emitter or the spectator is in the initial state. Thus, they are comprised of terms arising due to the change of the flavour or the partonic momentum fraction x of an initial state due to a splitting. The K-operator is given by [ 41 ] Kaa (x ; {αdip}) α = 2π − k K aa (x ) − KFS(x ) − It depends on the partonic x , and the flavour change from the Born process initial state flavour a to a of the convolution Eq. (2.6). Therein, the K collect universal terms present in all splitting involving an initial state as either emitter or spectator. Then, while K contains solely remaining terms from final state splittings in the presence of initial state spectators, the K t are their counterparts for initial state splittings in the presence of a final state spectator, i and k running over all final state partons in each case. K˜ contains solely related correlations between both initial states, arising from dipoles where both the emitter and the spectator are in the initial state. The AK terms collect all finite terms arising when any of α FI, α IF or α II is different from unity, thus restricting the phase space of the respective dipoles. Again, the Qi2k are the charge correlators of Eq. (2.5). Finally, KFS contains the factorisation scheme dependence. Currently, both only the MS schemes is supported, setting these terms identically zero. i 3.1 Identification of dipoles The P-operator now collects the remaining initial state collinear singularity from all dipoles involving initial states either as emitters or as spectators and cancels them against the collinear counterterm. Through this counterterm a dependence on the factorisation scale enters. The P operator is given by [ 40, 41 ] 2 Paa (x , μF ) α = 2π P aa (x ) k 2 Qa k log Only initial state splittings are present, either in the presence of a spectator in the final state, which is encoded in the sum of k, or with the opposite initial state b acting as the spectator. It otherwise only depends on the Alterelli–Parisi splitting function detailed in Appendix B. 3 Implementation The implementation of the QED generalisation of the Catani– Seymour dipole subtraction scheme in SHERPA’s matrix element generator AMEGIC proceeds along the lines of [ 49 ]. As in general real and virtual corrections of O (αsn αm ) contain divergences of both QCD and QED origin, both cases are included in this section. In the following, the general structure of the implementation is reviewed. The starting point to construct the involved subtraction terms in the Catani–Seymour subtraction formalism is a given flavour configuration in the Born or the real emission phase space and the perturbative order O(αsn αm ) in accordance with the respective virtual or real correction to be computed. For all parts, on-the-fly variations of both the factorisation scale μF and the renormalisation scale μR are available through an extension of the algorithm detailed in [ 56 ]. Differential subtraction terms. A given real emission configuration {ab} → {1, .., m + 1} at order O (αsn αm ) can in general exhibit both QCD and QED divergences simultaneously. The following therefore describes the identification of both types of dipoles. Thus, all triplets {i, j, k} that can be built from the external particles of the process are tested for the presence of an infrared divergence, QCD or QED, by checking for the existence of a dipole subtraction term. In these triplets i and k may be in the initial or final state while j may be in the final state only. Likewise, i = j , i = k, j = k and triplets that only differ in a permutation of i and j are considered identical. Then, the following steps are executed. 1. Based on the quantum numbers and flavours of the triplet it is decided whether a QCD, a QED or both splitting function can exist. A QCD splitting function can exist only if i , j and k are colour charged, while a QED splitting function can exist only if the charge-correlator Qı2j k˜ does not vanish. The dipole type is determined based on whether i and k are in the initial or final state. A given triplet {i, j, k} may exhibit both QCD and QED divergences, and thus may form both a QCD and a QED dipole. 2. The flavours ıj and k˜ are determined for each possible splitting function. 3. For each possible splitting function {ıj , k˜} → {i, j, k} the underlying Born configuration and its order, O(αsn−1αm ) in case of a QCD splitting function and O(αsnαm−1) in case of a QED splitting function, are determined. If, including the insertion of the appropriate colour-, chargeand spin-correlations, such a process at this order exists, a dipole subtraction term is built. If the above steps do not lead to a valid dipole subtraction term, no divergence can be present. The real emission, a conventional tree-level process, is grouped together with all its subtraction terms into one computational unit and their respective cross sections and observable values, O( m+1) and all O( mi ), are treated as correlated. Integrated subtraction terms. Similar to the above discussed real emission corrections, the virtual correction configurations {ab} → {1, . . . , m} at order O(αsnαm ) in general exhibits poles due to both QCD and QED origins. To subtract them, both QCD and QED integrated subtraction need to be included. In fact, they naturally arise as counterparts to differential subtraction terms constructed for the corresponding real emission correction, as guaranteed by Bloch and Nordsieck [ 57 ] or Kinoshita–Lee–Nauenberg [ 58,59 ] theorem. Consequently, QCD and QED I, K and P operators are constructed. While their QCD version are described in detail in [ 49,50 ], their QED version of Eq. (2.6) are discussed below. The I operator, on the one hand side, has the same initial state flavours and momentum fractions as the virtual correction. The K and P operators on the other hand, resulting from the combination of the integrated subtraction terms and the collinear counterterms, involve a summation over possible initial state flavours and comprise the following general structure in their dependence on the additional x integration variable [g(x )]+ + δ(1 − x )h(x ) + k(x ). Therein, both h(x ) and k(x ) are regular functions in x , while the plus distribution of g(x ) is defined as (3.1) dx f (x ) [g(x )]+ = dx [ f (x ) − f (1)] g(x ). (3.2) 0 1 Hence, the potentially computationally intensive matrix elements in B have to be calculated twice for every phase space point in addition to the flavour summation. To remedy this, the original expression of Eq. (2.6) is reformulated. Dropping the dependence on the {αdip} parameters and explicitly stating the dependence of the underlying Born term on the initial state momentum fractions and parton densities, Bab = fa fb Bab, the K and P operators it can be recast to a,b = a,b with fa (η) being the parton density of flavour a and momentum fraction η in the proton, otherwise absorbed in the Born term B. All PDFs are evaluated at the same scale μF. Gab(η) = 0η dx gab(x ) are the analytically computed divergence free parts of the integral of the functions under the plus distribution. Its divergence at x = 1 is cancelled numerically on the second last line. Effectively, this reformulation results in a redefinition of the PDF for incoming parton a of Bab. As the integrand of the remaining integral over xa is a simple one-dimensional function without a pronounced structure, its numerical evaluation is very stable and can be accomplished by a single point for each summand per phase space point d (m4). The K and P operators for the second incoming parton b are recast similarly. The thus transformed form of the K and P operators require only a single evaluation of the potentially costly matrix elements in B while retaining the number of computations of PDFs needed, speeding up the computation considerably for involved processes. This allows to generate and evaluate the underlying Born matrix element for both the I operators and the K and P operators at the same time. In fact, due to these three operators being simple multiplicative scalars, the common underlying Born matrix element is identical to the standard Born matrix element, allowing for their simultaneous calculation at no extra cost. The operators themselves are then built from dipoles constructed from all doublets {i, k}, i = k, of external partons available in the partonic process for which the charge-correlator Qi2k does not vanish. 3.2 External photons External photons can play different roles in a calculation: they can either be resolved or unresolved. According to this distinction they should be treated differently at NLO EW [ 31,60 ]. Initial state photons are always unresolved at a hadron collider. They thus should be treated in a short-distance scheme, allowing them to split into massless fermions, necessitating a proper subtraction of infrared divergences. Final state photons on the other hand can play both roles. If they are considered resolved, they should be treated in an on-shell scheme and no explicit photon splitting is allowed. Concerning the dipole subtraction discussed in this paper they are thus neutral particles and do not form part of a dipole, except as possible recoil partner of another unresolved photon. A final state unresolved photon, on the other hand, again must be treated in a short-distance scheme, necessitating their splittings to be subtracted. As discussed in Sect. 2.1, the dipole picture is not necessary to capture the divergences of photons splitting in pairs of massless fermions, leptons and quarks, due to the absence of soft divergences that necessitate the correlation with the emissions off other partons of the event. Nonetheless, it offers a possibility to assign one or more recoil partners to absorb the transverse momentum of the splitting, thus fitting these purely collinear splittings into the dipole picture. The choice of spectator is essentially arbitrary, and all other partons of the event offer being viable spectators. Only the following condition has to hold k˜=ıj κıj k˜ = −1 for every splitting photon ıj , cf. 2.5. Therein, the κıj k˜ are arbitrary and possibly dynamic weights assigned to every dipole with spectator k˜. In practise, for initial state as well as final state splittings five choices ckγ have been implemented: ˜ 0. only allow initial state partons as spectators, 1. only allow final state partons as spectators, (3.4) 2. only allow QED charged particles as spectators, 3. only allow QED neutral particles as spectators, 4. allow all particles as spectators. The κıj k˜ are set to the phase space point independent value −ns−p1ec, with nspec the number of assigned spectators, thus trivially fulfilling Eq. (3.4). 3.3 Process mappings In general, physical cross sections include multiple different partonic channels. However, many of these partonic channels share identical squared matrix elements, potentially differing by constant factors. They, thus, do not need to be recomputed for every flavour channel but can be reused. As a typical example, consider the production of a lepton pair in association with two jets. The process g d → e+e− g d shares a common squared matrix element with g s → e+e− g s and g b → e+e− g b on Born level at O(αs2α2). Hence, the squared matrix elements of the latter two partonic processes are mapped on the first, reusing its computed value. In this way, of the 95 partonic channels of this process, only 30 have to be computed. Further mappings of individual graphs and subgraphs are implemented but are not discussed in the following, see [ 47,48 ]. In SHERPA’s matrix element generators various forms of process mappings are implemented to reduce the computational complexity and memory footprint, both for the virtual and real emission corrections, cf. [ 47–50 ]. However, while for NLO QCD corrections to the leading order Born process it is true that if the Born process is mappable onto another existing process, then also both the virtual correction and the insertion-operator-augmented colour-correlated underlying Born process of the integrated subtraction term are mappable to the same process. This is no longer true when considering the NLO EW or NLO QCD corrections to subleading Born processes. In the present implementation of a full NLO QCD and NLO EW subtraction in the matrix element generator AMEGIC the following process mapping strategy is followed for the real subtracted contributions: • The real emission process and its associated dipole subtraction terms are grouped in one computational unit. A given partonic channel of the real emission process can be mapped onto another already existing one if both processes consist of the same diagrams and all involved (internal and external) particles have the same masses and widths and the same underlying interaction (coupling factors may differ by a constant). Is this the case the whole computational unit can be mapped and the result of the mapped-to process can simply be reused. • Individual dipoles can be mapped if the emitter, emittee and spectator indices are identical, and the underlying Born process can be mapped according to the above rules. In this case, the result of the mapped-to process can simply be reused. • Underlying Born matrix elements can be mapped if the Born-level emitter ıj carries the same indices and the underlying Born process itself can be mapped. The spin correlation insertion operator, needed if parton ıj is either a gluon or a photon and described in [ 49 ], is encoded in the calculational routines and necessitates the above restriction. Can the underlying Born process be mapped, only the calculational routines can be shared, reducing the memory footprint, but due to the potentially differing underlying Born momenta its result has to be recomputed. The strategy for the virtual subtracted contributions reads as follows: • The virtual correction process, interfaced from an external virtual correction provider, as well as its associated integrated dipole subtraction terms and the collinear counterterm, taken together and reformulated according to Sect. 3.1, are grouped into one computational unit. If the underlying Born processes of the integrated subtraction contribution (in all orders required) can be mapped according to the above rules and the virtual correction provider confirms that the virtual correction can be mapped onto the same process that SHERPA’s tree-level matrix element generator maps the underlying Born processes onto, the whole computational unit is mapped. In this case, the result of the mapped-to process can simply be reused. The K and P operators, whose internal PDF factors depend on the initial state flavour, still have to be recomputed, but their matrix element coefficients are cached. • If the virtual correction provider cannot confirm the mapping performed for the underlying Born process, only the underlying Born process is mapped and virtual correction is recomputed. Again, the K and P have to be recomputed, but their matrix element coefficients are cached. Here, efficiency is lost if the virtual correction provider uses less efficient process mappings. 3.4 Fiducial phase space definition Phase space restrictions are an essential part of the implementation of a framework for automated NLO calculcations. These cuts, however, need to be applied in an infrared safe way. At NLO, they must not discriminate between a massless parton before and after its collinear splitting or before and after a soft gluon or photon emission. Thus, if QCD singularities are present, massless QCD charged particles must be clustered into jets before any further cuts are applied. Similarly, in case of the presence of QED divergences, massless charged particles must either be dressed with the surrounding photons or be included in the jet algorithm. Subtleties arise in the presence of both QCD and QED singularities simultaneously. Here, usually, only a fully democratic jet finding can consistently treat all singularities, although specialised solutions exist for simplified situations. Massive QCD and QED charged particles may be treated as bare as their mass shields the collinear singularity, but can also be included into jet finding and dressing algorithms. An intermediate scheme which includes only the logarithms of the parton mass needed to regulate the collinear divergences [ 45 ], but are otherwise treated massless in the calculation, is not implemented. To this end, the implementation in SHERPA is equipped with a range of algorithms to define infrared safe quantities on which further restrictions can be applied. Multiple such selectors can be nested. DressedParticleSelector This selector takes a choice of dressing algorithm (cone or sequential recombination) and (flavour-dependent) dressing parameters (cone radius or radial parameter and exponent). All charged particles of the process are then dressed with all photons using the specified algorithm with the given (flavour-dependent) parameters. Their four momenta are added such that four momentum is conserved. The dressed charged particles may no longer be on-shell and the photons used to dress the charged particles are removed from the list of particles. The resulting list of particles and their momenta are then passed to all subselectors. Jet_Selector This selector uses FASTJET [ 61 ] to build jets from a given list of input particles. It takes a list of flavours that are considered as jet finding input particles, the jet finding algorithm and its parameters including phase space boundaries in p⊥, η or y, as well as a minimal and maximal number of jets to be found as arguments. Additionally, clustered jets can be tagged or anti-tagged based on their flavour content, including relative and absolute constituent momentum requirements (e.g. b tagging a jet if one of its constituents is a b-quark or anti-γ -tagging a jet if its photon constitents carry in excess of zthr fraction of the total jet momentum). The clustered jets as well as all particles not used as jet finding input are passed on to all subselectors. Isolation_Selector This selector uses the smooth cone isolation of [ 62 ] to isolate particles of a given flavour against particles of another flavour. It takes both the isolation flavour and rejection flavour list as well as the algorithm parameters as input. It further can be specified how many isolated particles of the given flavour should minimally and maximally be found in a given phase space volume (bounded by p⊥ and η or y ranges), e.g. exactly two isolated photons with p⊥ > 30 GeV and |η| < 2.35. The list of found isolated flavours and all other particles except for those that should be isolated, but are not, are passed to all subselectors. Once an infrared safe definition of particles is found, the standard single- or multi-particle selectors such as PT (implementing transverse momentum requirements on a given flavour), PTmis (missing transverse momentum build from all neutrinos), or DR (angular distance between the two given particles) can be used. Further, SHERPA is equipped with a user hook system, providing the possibility for users to implement arbitrary routines for phase space cuts and dynamically load them at run time, without the need to modify their SHERPA installation. Identified particles should in principle be defined using fragmentation functions, denoted Dij (z) for finding parton i in parton j at momentum fraction z. As, however, all Dii (z) have a δ(1 − z) as leading term in on-shell renormalisation schemes (and only differ by ratios of couplings in other schemes) simplified schemes exist that are applicable to many practical situations. Hence, no fragmentation function is implemented yet. 3.5 Flavour scheme conversion All publicly available PDF sets containing QED effects in their evolution are fitted with five (MRSTqed [ 63 ], CT14qed [ 64 ], LUXqed [ 65,66 ], HKR16 [67], NNPDF3.0qed [ 68 ], NNPDF3.1luxqed [ 69 ]) or six (NNPDF2.3qed [ 70 ]) light flavours. Thus, for next-to-leading order calculations with four or less light flavours the following scheme-conversion terms need to be added for consistency [ 71 ] O (NnLlfO) QCD(μ2R, μ2F) = O (NnLf)O QCD(μ2R, μ2F) + d m × p log m2 i μ2R μR − mi2 − 2 nf i=nlf {ab} m2 agbg log μ2Fi αs TR 3π μF − m2 2 i Ba(nbf)( m ; μ2R, μ2F) O( m ) (3.5) and O (NnLlfO) EW(μ2R, μ2F) = O (NnLf)O EW(μ2R, μ2F) − d m nf i=nlf {ab} α 3π NC Q2 i aγbγ log m2 i μ2F Therein, the O (NnL)O QCD/EW is the expectation value of an arbitrary observable O computed at NLO QCD or NLO EW, respectively, with n-flavour scheme parton densities and the corresponding strong coupling summed over all initial state contributions. Ba(nb) is the Born term, as defined in Sect. 2, in the ab channel. The sums run over all (nf − nlf) flavours of mass mi that are part of the nf-flavour PDF parametrisation but not the nlf-flavour scheme of the calculation and all combinations of partonic channels ab occurring in the nlf-flavour scheme, respectively. p is the power of the strong coupling in the Born process and agbg = δga + δgb, i.e. it takes the values 2 in the gg channel, 1 in all qg and q¯ g channels, and zero otherwise. aγbγ is defined analogously. Each logarithm of course only contributes if the scale is larger than the respective mass. As the electroweak coupling is not taken running in the common renormalisation schemes it is independent of the number of massless flavours in the calculation. For MSlike renormalisation schemes a similar term proportional to its power in the Born process is to be added. 4 Checks of the implementation The implementation of the formalism described in the previous sections in SHERPA needs to be validated. To this end, both its independence of its internal free parameters, {αdip} = {α FF, α FI, α IF, α II}, κ and the choice of spectator in photon splittings, and its agreement with independent implementations for fixed values of these parameters need to be tested. While the latter were carried out in [ 19,21,25,31,32 ] against the private implementations in Munich [ 74 ], in [ 51 ] against MadGraph5 [ 75 ], and in [ 51,53 ] against Recola [76] and found full agreement, the former represent a powerful check of internal consistency. Further, as the {αdip} parameters regulate the phase space coverage of the differential subtraction terms, they can be used to lower the average number of contributing dipoles for a given real emission contribution and, thus, reduce the computational costs of the real-subtracted contribution. There are now two independent aspects of the calculation that need to be checked: (a) The I operator of Sect. 2.2 containing the explicit Laurent expansion in as → 0 needs to reproduce the correct −2 and −1 coefficients in order to cancel all corresponding poles of the virtual matrix elements, leading to a finite integrand in Eq. (2.2) in d = 4. These checks were performed for all possible dipoles in [ 19,21,25,31,32,34,37,51–53 ] and are not repeated here. (b) The expectation value of any infrared observable is independent of the choice of the technical {αdip} and κ parameters as well as the choice ckγ of spectator for photon ˜ splittings. It thus needs to be verified that the sum of the contributions of the 0 coefficient of the integrated subtraction terms and the differential subtraction terms is independent of these parameters and choices. In order to arrive at a finite result in d = 4, cf. Eq. (2.2), the real emission correction and the collinear counterterms are added. The corresponding quantity is defined as O IRD = d (4) m m ; {αdip}, κ, cγ ) ID( (4) k ˜ + C( (m4)) and will be evaluated for processes containing all available dipole configurations in the following. In the following, the inclusive or fiducial cross section contribution σIRD is evaluated for a range of different processes, testing all possible dipole configurations. Throughout the input parameters of Table 1 are used in the Gμ scheme, although several other EW input parameter schemes are available in general. Similarly, all unstable particles are treated in the complex-mass scheme [ 77 ]. Further, the CT14nlo [ 78 ] and CT14qed [ 64 ]1 PDF sets with five active and massless flavours and their corresponding αs parametrisations with αs (m Z ) = 0.118 are used. While the use of the CT14nlo PDF set for NLO EW calculations is, strictly speaking, inconsistent due its missing QED evolution of the initial state quarks, it offers the quantification of the importance of photon initiated processes in these technical comparisons. To the same end, the CT14qed PDF set is not employed using its the best fit value for the intrinsic inelastic photon momentum fraction at the reference scale of Q = 1.295 GeV p0γ = 0.05%, but rather evaluate it for the extremes of its 1σ uncertainty, p0γ = 0 and p0γ = 0.14%, where applicable. In processes where jets need to be constructed to define a fiducial phase space volume, anti-kt jets with R = 0.4 and p⊥ > 30 GeV are used [ 80 ] and both partons and photons are considered as constituents. Similarly, if leptons need to be defined for the same purpose, they are dressed with all photons in a cone of R = 0.1. The number of phase space points in the computation of σIRD({αdip}, κ, cγ ) is kept k ˜ constant for each process, such that the indicated statistical 1 To be precise the CT14nlo and CT14qed_inc_proton PDF sets interfaced through LHAPDF [ 79 ] 6.2.1 are used. Table 1 Numerical values of all input parameters. The gauge boson masses are taken from [ 72 ], while their widths are obtained from stateof the art calculations. The Higgs mass and width are taken from [ 73 ]. The top quark mass is taken from [ 72 ] while its width has been calculated at NLO QCD. In calculations where a massive particle is present as an external state, its width is set to zero uncertainty can be interpreted as a measure of the change of convergence of the subtraction with respect to the variation of the technical parameter choices. Massless dipoles. Contrary to the QCD case, in the Standard Model almost all particles carry QED charges and therefore participate in the construction of dipoles. One notable exception are neutrinos. Therefore, in order to investigate the behaviour of the massless II, IF, FI and FF dipoles, and, thus, the independence of α II, α IF, α FI and α FF separately in this technical validation σIRD is considered for all three different rotations of the interaction of a quark-anti-quark pair with a neutrino-anti-neutrino pair. Hence, besides the σIRD contribution to the O(α) correction to pp → νeν¯e at the LHC at an invariant mass of 13 TeV, σIRD is computed for both the production of at least two jets at a hypothetical νe − ν¯e collider at a centre-of-mass energy of 1 TeV and inclusive single jet production in equally hypothetical νe p deep inelastic scattering (DIS) with the same centre-of-mass energy are calculated. Figure 2 now details σIRD for all three setups. νeν¯e → j j production, detailed in the left hand side plot, comprises only FF dipoles and, thus, only depends on α FF. Varying its value over four orders of magnitude leads leaves the value of σIRD unchanged within the statistical uncertainties. The black line is placed at the central value of the computation with the smallest statistical uncertainty to guide the eye. As is evident, lowering the α parameter too much, i.e. restricting the subtraction to act only on configurations very close to the divergence, results in large cancellations between the real-subtracted and the integrated dipole contributions of Eq. (4.1), degrading the statistical prowess of the calculation. Similarly, pp → νen¯ue production, detailed on the right hand side, comprises only II dipoles and, thus, depends on α II only. Showing the results for all three choices of photon content, contributing directly through the respective γ q/γ q¯ channels as well as indirectly through the impact on the quark PDFs and, through momentum conservation, on the gluon PDF, a similar picture as for νeν¯e → j j emerges. The three coloured lines now indicate the central value of the calculation with the smallest statistical uncertainty for each PDF choice. In each case, σIRD is found to be stable under the variation of α II. Finally, the centre plot shows the correction contribution for the hypothetical DIS process νe p → νe j requiring at least one jet in the lab frame. As FI and IF dipoles always occur in pairs the α FI and α IF dependence is evaluated together here. The resulting σIRD(α IF = α FI) are also found to be stable when varying both parameters simultaneously over four orders of magnitude. Moving away from the simplest configurations, Fig. 3 displays the results for electron-positron pair production. The left hand side plot again displays their production at the hypothetical νe-ν¯e collider used before, finding very similar results and their independence of α FF. The centre plot now, however, displays the production of an electron-positron pair at the LHC. At leading order, this process proceeds through qq¯ → e−e+ and γ γ → e−e+ at O(α2). Consequently, the qq¯ channel exhibits six dipoles of all four types. In the γ γ channel, the number and types of dipoles present depends on the choice of possible photon splitting spectators cγ . To reguk late all LO singularities the fiducial phase space is ˜defined by requiring the dressed electrons to have a transverse momentum of at least 20 GeV and the pair to have an invariant mass of at least 60 GeV. As σIRD now potentially depends on the full set {αdip} no continuous variation over four orders of magnitude is performed. Instead, each of the four parameters is varied independently to 0.001 keeping all others at their default value of 1. These four variations are completed by setting α FF = α FI = α IF = α II = 1 and 0.001. The resulting correction contributions are found to be independent of {αdip}. Finally, the right hand side plot of Fig. 3 displays the dependence of σIRD on the choice of spectators in photon splittings. Only the γ q/γ q¯ channel is considered as both the γ γ and qq¯ channels are independent of this choice in this process. The γ q/γ q¯ channel, however, still receives contributions from photon radiation off quarks in addition to the sought after photon splittings into quark-antiquark pairs. Hence, the minimum invariant mass is raised to 2 TeV to increase the relative importance of the photon PDF, enhancing the photon splitting contribution. The resulting correction contribution at O(α3) is found to be independent of all five choices of photon splitting spectators available. Please note that for this process scheme 0 and 3 as well as scheme 1 and 2 lead to the same set of allowed spectators, respectively, and therefore to identical results. Massive dipoles. Massive particles are so far only allowed in the final state.2 Dipoles involving massive partons, either as emitter or spectator, comprise only three types: FF, FI and IF. Emittees are always considered massless, otherwise no singularity would be present. In Fig. 4 again top-anti-top pair production at a hypothetical νe-ν¯e collider and single top production at a hypothetical νe- p collider is considered in order to study the individual dipoles separately. In the left plot, the α FF (in)dependence of the O(α3) correction contribution σIRD, containing only massive FF dipoles, is shown. It exhibits the familiar picture of decreasing statistical prowess of the calculation with too small α FF, but otherwise consistent results. The right plot details the α IF and α FI dependent massive dipoles in the hypothetical DIS scenario. As before, α IF and α FI are varied simultaneously for this purpose and the independence of the corrections contribution on both parameters is observed. The left plot of Fig. 5 now investigates the dependence of σIRD on the κ parameter. In only arises in FF dipoles of gluons splitting into gluons or massless quarks or photons splitting into massless fermions in the presence of a massive spectator. Consequently, to restrict the number of additional contributions, top-anti-top pair production in association with one jet at the hypothetical νe–ν¯e collider is considered. At LO, this process contributes is defined both at O(αs α2) and O(α3) where the final state photon forms the jet. At NLO, there are contributions at O(αs2α2), O(αs α3) and O(α4). The O(αs α3) contribution, however, contains 2 For advances for NLO calculations with initial state massive particles see [ 44,81,82 ]. pp → e+e− in the γ q/γ q¯ channel on the choice of recoil partner for the initial state photon splittings cγ . The invariant mass of the electron pair k ˜ is raised to increase the contribution of the photon induced channels. The qq¯ and γ γ channels comprise no dipoles with splitting photons Fig. 5 Left: κ-dependence of the sum of integrated subtraction term and differentially subtracted real emission for νeν¯e → t t¯ j . Right: dependence of the sum of the integrated subtraction term and differentially subtracted real emission for νeν¯e → t t¯ j at O(α4) on the choice of recoil partner for the final state photon splittings cγ k ˜ neither gluon nor photon splittings. Due to the relative size of g → gg and g → qq¯ splittings in relation to gluon radiation off the top quarks the κ dependence of the O(αs2α2) contribution is much more pronounced than at O(α4), where photon radiation off the top quarks overwhelms the photon splitting contribution. Nonetheless, at both orders an independence of the σIRD of κ is found. In addition, the right plot investigates the influence on the choice of spectators for the final state photon splitting ocurring at O(α4). No dependence on this choice is observed. Please note that for this process scheme 0 and 3 as well as scheme 1 and 2 lead to the same set of allowed spectators, respectively, and therefore to identical results. Figure 6 now considers top-anti-top pair production at the LHC. This process occurs at LO at O(αs2), O(αs α) and O(α2). Thus, at NLO there exist four contributions, at O(αs3), O(αs2α), O(αs α2) and O(α3). While the O(αs3) and O(α3) terms can be clearly identified as NLO QCD and NLO EW corrections to the LO O(αs2) and O(α2) expressions, respectively, the O(αs2α) and O(αs α2) terms do not possess such a unique characterisation: they act as both NLO QCD corrections and NLO EW corrections to different LO processes. Consequently, their divergence structure contains singularities of both QCD and QED origin. Hence, both QCD and QED dipoles with underlying Born processes of different orders are needed for a full subtraction of all divergences. As explained in Sect. 3.1, this holds true both for the differential and integrated subtraction terms. It therefore serves as an additional check to verify the independence of the result of the {αdip} for each O(αs3−m αm ), m = 0..3, individually. And indeed, the correction contribution σIRD is found independent of {αdip} for each such order. External W bosons. Lastly, it may become necessary to also consider external W bosons (or other massive charged particle with spin > 21 in BSM theories) as stable final state particles, e.g. to reduce the computational complexity for high final state multiplicity processes where off-shell effects and effects in the decays can be ignored or recovered through other means [ 83 ]. In this case the literature does not provide expressions for the respective massive dipole functions. As their mass, however, can be assumed to be large enough to suppress collinear radiation well enough, this only leaves the spin-independent soft photon emission limit. 3 Here, both the expressions for the radiation of a photon off a massive scalar or a massive fermion can be used. 3 This assumption also proves to be correct for on-shell W production at a 100 TeV collider. Fig. 7 {αdip}-dependence of the sum of integrated subtraction term and differentially subtracted real emission for νeν¯e → W +W − and pp → W +W −. Both choices of subtraction terms for external massive Figure 7 details the production of a W +W − pair in the hypothetical νe − ν¯e collider, separating the FF dipoles and their α FF dependence, on the left hand side. As before, σIRD is found to be independent within the statistical accuracy and also independent of the whether the massive scalar or massive fermion subtraction terms are used. The right hand side focusses on their production at the LHC. Again, all dipoles contribute at O(α3), leading again to the afore described sixpoint variation. Also in this case, the result is independent of {αdip} and the choice of scalar or fermionic subtraction term. The FI and IF dipoles cannot be investigated separately, as in the νe p → νet case, due to charge conservation. 5 Conclusions This paper detailed the construction and implementation of the adaptation of the Catani–Seymour subtraction formalism for NLO EW calculations. Besides the translation of the QCD dipole functions to the QED case, several other issues have been addressed. They include the special role photon splittings play in the formalism, embedding extermal massive emitters of spin > 21 into the formalism and the interplay of QCD and QED subtractions for processes exhibiting both kinds of divergences. The resulting general subtraction for NLO EW calculations has been implemented in the SHERPA Monte-Carlo event generator framework. Interfaces to OPENLOOPS, GOSAM and RECOLA to access the needed virtual corrections exist and are fully functional. In addition to the checks against independent implementations on the level of partial and total cross sections performed in previous publications, numerous internal cross checks for independence of technical parameter choices, {αdip} = {α FF, α FI, α IF, α II}, κ and the choice of spectator in photon splittings ckγ , have been presented here. This imple˜ charged vector bosons, using the ones of an external massive scalar (filled symbols) and massive fermions (empty symbols), are compared mentation will become publically available in the near future with the next major SHERPA release and an extension to the COMIX matrix element generator [ 84 ] is foreseen. Acknowledgements MS would like to thank N. Greiner, S. Höche, S. Kallweit, J. Lindert, and S. Schumann for numereous checks of the implementation and S. Pozzorini for many useful discussion. MS acknowledges support by the Swiss National Foundation (SNF) under contract PP00P212855 as well as the European Union’s Horizon 2020 research and innovation programme as part of the Marie SkodowskaCurie Innovative Training Network MCnetITN3 (Grant agreement no. 722104) and the ERC Advanced Grant MC@NNLO (340983). MS would also like to thank the University DuisburgEssen for the kind hospitality during essential parts of the project. All plots are based on RIVET’s [ 85 ] plotting tools. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3. A Differential splitting functions This appendix details the complete functional form of the dipoles introduced in Eq. (2.3), Di j,k , Diaj , D j,k and a Daj ,b. The functional forms of the spin-dependent insertions Vi j,k , Viaj , Vaj,k and Va,b and the kinematic maps of the difj ferential splitting functions are taken from the original QCD case [ 40, 41 ]. They are summarised in the following. A.1 Final–final dipoles Dipoles with both emitter and spectator in the final state take the form 1 2 Di j,k = − ( pi + p j )2 − mı2j Qıj k˜ m . . . , ıj , . . . , k˜, . . . |Vi j,k | . . . , ıj , . . . , k˜, . . . m . Therein, the charge-correlator is defined in Eq. (2.5). All momenta of the dipole are on-shell pi2 = mi2, p2j = m2j , pı2j = mı2j , pk2 = pk2˜ = mk , and the total four momentum flowing through it is given by q = pi + p j + pk = pıj + pk˜ . It is thus invariant under the emission. The momenta of the parton before and after the splitting are connected through the map pk˜ = λ q2, mı2j , m2 k ˜ λ q2, ( pi + p j )2, mk2 1 + 2 q2 + m2 k˜ − mı2j q q2 pıj = q − pk˜ q pk pk − q2 q wherein mk˜ = mk and the Kallen function is defined as λ (a, b, c) = a2 + b2 + c2 − 2ab − 2ac − 2bc. The splitting variables read pi p j , yi j,k = pi p j + pi pk + p j pk pi pk zi = 1 − z j = ( pi + p j ) pk zi(,mj) = zi, j − 1 − vi j,k , 2 and which are defined on the intervals yi j,k ∈ [y−, y+] and zi ∈ [z−, z+]. The boundaries are given by (A.1) (A.2) (A.3) (A.4) massless. The relative velocities between ıj and k˜, i + j and k, and i + j and i are given by vıj,k˜ = vi j,k = vi j,i = ! ! and They are introduced to facilitate the analytic integration and only take effect away from the singular limit. Thus, denoting the Spins of parton ıj in . . . , ıj , . . . |Vi j,k | . . . , ıj , . . . by s, s (if ıj is a fermion), μ, ν (if ıj is a photon), or omitting them (if ıj is a scalar) the dipole functions are defined as s|Vγ f,k |s It is thus clear that only in the case of photons splitting into massless fermions does the dipole splitting function have a non-diagonal structure. In all other cases, the insertion . . . , ıj , . . . |Vi j,k | . . . , ıj , . . . reverts to a simple multiplication of a real-valued splitting function and a standard Born matrix element. The parameter κ controls finite terms and only takes effects in the case where the respective fermionic products of the splittings are massive. Setting κ = 0 somewhat simplifies the spin-dependence of the differential dipoles. y− = 2mi m j 2 qi j,k z 2mi2 − qi j,k yi j,k 2 − = 2(mi2 + m2j + qi2j,k yi j,k ) 2mk q2 − m2 k y+ = 1 − z 2 qi j,k 2mi2 − qi j,k yi j,k 2 + = 2(mi2 + m2j + qi2j,k yi j,k ) 1 − vi j,i vi j,k 1 + vi j,i vi j,k . (A.6) A.2 Final–initial dipoles with qi2j,k = q2 −mi2−m2j −mk2. As can be seen, a divergence, residing at yi j,k = 0, is only present if either i or j are Dipoles with the emitter in the final state and the spectator in the initial state take the form pi2 = mi2, p2j = m2j , pı2j = mı2j , pa2 = pa2˜ = 0, and the total four momentum flowing through it is given by q = pi + p j − pa = pıj − pa˜ . It is thus invariant under the emission. The momenta of the parton before and after the splitting are connected through the map pa˜ = xi j,a pa pıj = q + pa˜ . The splitting variables read xi j,a = 1 − ( pi + p j ) pa pi pa zi = 1 − z j = ( pi + p j ) pa . 1 pi p j − 2 mı2j − m2 j i − m2 and The singularity of the splitting resides at xi j,a = 1 and is only present if either mıj = mi or mıj = m j . Adopting the above convention for labeling the emitter’s spins the dipole functions are defined as s|Vγa f |s δss − ( pi + p j )2 × zi piμ − z j p μj zi piν − z j pνj . |Vγ s | = 8π μ2 α a 2 ms2 2 − xi j,a − z j − 2 − pi p j Again, only photon splittings exhibit any spin correlations. In all other cases the dipole function and the underlying Born matrix element factorise. (A.9) (A.10) (A.11) (A.12) (A.13) (A.14) A.3 Initial–final dipoles Dipoles with the emitter in the initial state and the spectator in the final state take the form and As the splitting partons are all in the initial state and therefore massless, the singularity at u j = 0 is present in any case. Adopting the above convention for labeling the spin of the emitter aj the dipole functions are defined as f s|Vγ,k |s × p μj pkμ u j − uk γ s|V f,k |s A.4 Initial–initial dipoles Dipoles with both emitter and spectator in the initinal state take the form before and after the emission, respectively. Although q acquires transverse momentum after the emission, q2 = q˜ 2 holds, i.e. the mass of the dipole is invariant. The momenta of the partos before and after the splitting are connected through the map paj = xai,b pa pb˜ = pb. While the initial state cannot absorb recoil transverse to the beam access, it is transferred to the final state through a collective boost kı˜μ = kiν δνμ − 2 (q + q˜ )ν (q + q˜ )μ + 2 qν q˜ μ (q + q˜ )2 q2 of all final state partons. The splitting variables read xa j,b = 1 − p j ( pa + pb) pa pb pa p j and v j = pa pb , vb = 1. As the splitting partons are all in the initial state and therefore massless, the singularity at v j = 0 is present in any case. Adopting the above convention for labeling the spin of the emitter aj the dipole functions are defined as s V f,b s | γ | − 2xa j,b(1 − xa j,b) δss μ|V ff,b|ν = 8π μ2 α −gμν xa j,b − 1 − xa j,b 2v j vb xa j,b p j pb (A.21) (A.22) (A.23) (A.24) (A.25) (A.26) i = γ i = γ . (B.1) (B.2) × p μj pbμ v j − vb pνj pbν v j − vb (A.27) As only massless particles are considered as initial states, no splitting functions involving massive scalars are needed. B Integrated splitting functions This appendix summarises the complete functional form of the integrated dipole terms cast in the form of the I, K and P operators [ 40,41 ]. B.1 Terms of the I perator Each dipole contribution to the I operator reads Iik ( , μ2; κ, {αdip}) = Qik Vik ( , μ2; κ) + i ( , μ2) 2 by Wherein, Vi(kS) contains all infrared singularities (double poles for massless fermion splittings, single pole for massive split(NS) that incorporates all tings), and a non-singular piece Vik finite dependences on the masses of the emitter and the spectator. The latter vanishes in the massless case, see below. In the case of photon splittings, only the non-singular part survives, again being non-zero only in the case of massive spectators. The singular part vanishes in case of photon splittings. The singular term Vi(kS) is symmetric in i and k, it is given 1 " qi2k # Vi(kS)( ) = vik sik (mk )) (B.3) sik + γi 1 + ln + Ki + AiIk ({αdip}) + O( ) , Vi(kS)( ) + Vi(kNS)(κ) − π32 cf. Eq. (2.8). Therein, Vik takes the general form Vik ( , μ2; κ) = ⎧ Q2 μ2 ⎨⎪ i sik ⎪⎩ Vγ(NkS)(κ) with ρi/k = vik = and the velocity factors (1 + vik )sik + 2mi2/k sik parametrise the mass dependence. Therein, qi2k = ( pi + pk )2 = sik + mi2 + mk2 is the dipole invariant mass. In d = 4 − 2 dimensions, > 0, the massless limit is approached smoothly. Subsequently taking the limit → 0, however, results in −2 poles only if at least one of both particles making up the dipole is massless. Completely massive dipoles only result in single infrared poles. (NS) contributions encode the addi The non-singular Vik tional finite mi/k > 0 contributions. They vanish in the limit that both mi → 0 and mk → 0 and their precise form can be found in Eq. (6.21)–(6.26) of [ 41 ]. Of particular interest is the term for a photon emitter in the presence of a massive spectator, (NS) Vγ k (κ ) = γγ log qsii2kk − 2 log |qik||qi−k |mk − |qik2|m+k mk " sik # qi2k 2 + 2γγ κ − 3 m2 k log sik 2mk |qik | + mk , with |qik | = qi2k . Thus, choosing κ = 23 somewhat simplifies the integrated subtraction term. Further, Qi2 1 + 21 ln mμ2i2 − 2 i = massive fermion Qi2 1 + ln mμ2i2 − 2 i = massless fermion i = γ i = massive scalar. (B.4) (B.5) ⎪⎩⎪⎪⎪ Therein, the K collect universal terms present in all splitting involving an initial state as either emitter or spectator. They are defined as and are independent of the particle masses. For the photon, the sums run over all charged massless flavours of the model, NC, f being their respective colour multiplicity. Lastly, the dependence of the I operator on the technical parameters {αdip} = {α II, α IF, α FI, α FF} is encoded in the last term AiIk ({αdip}) which is detailed in App. C. B.2 Terms of the K and P operators The K operator of Eq. (2.9) reads α Kaa (x ; {αdip}) = 2π − k Qa2 k K ata ,k (x ) − Qa2 b K˜ aa (x ) + AaKa ({αdip}) . In case of the photon, the sum runs over all massless charged flavours present in the model. The flavour constants γi and Ki are given by (B.7) ⎧ Praeag (x ) = ⎨⎪⎪ −Qa2 (1 + x ) aa Qa2 [x 2 + (1 − x )2] aa ⎪⎪⎩ 0NC,a Qa2 1+(1x−x)2 aaaa = f f = f γ = γ f = γ γ , i + (B.8) (B.9) (B.11) K aa (x ) = Praeag (x ) log + δaa Qa2 1 − x + Pˆ aa (x ) x 2 log 1 − x 1 − x x 5 −δ(1 − x ) γa + Ka − 6 π 2 Qa2 . (B.10) As can be seen, they are independent of the final state and thus independent of any final state particle masses. Here, the regular part of the Altarelli–Parisi splitting functions are given by with NC,a the multiplicity in the QCD representation of parton a (3 for quarks, 1 for leptons). Similarly, the -dependent part of the splitting functions is defined as Pˆ aa (x ) = ⎧ Qa2 (1 − x ) aa = f f ⎪⎪⎨ Qa2 x aa = f γ NC,a Qa2 x (1 − x ) aa = γ f ⎪⎪⎩ 0 aa = γ γ . Thus, for a = a = γ the K term takes the simple form K γ γ (x ) = − 38 γγ δ(1 − x ), and comprises an end-point term only. The factorisation scheme dependent terms KFS vanish in the MS scheme, and the Ki,aa , containing remnants from final state splittings where the initial state a forms the spectator. As it also depends on the final state flavour i , its forms are given as follows. The contributions, listed separately for scalars, fermions and photons, read (B.12) (B.13) K f,γ f (x) = 0 K f, f f (x) = 2 log(1 − x) 1 − x + − log(2 − x) 1 − x a mi + Jγ Q x, √sik Ks,aa (x ) = K f,aa (x ) − δaa Simultaneously, for i = γ Kγ,aa (x ) = −δaa γγ due to the different normalisation of Qi2k in the photon case. The K ata ,k terms comprise of initial state splittings in the presence of a final state spectator. They take the form K ata ,k (x) = Praeag (x) log (1 −(1x−)sxak)s+ak mk2 + γa δaa δ(1 − x) ⎡ × ⎣ log sak − 2mk sak + mk2 + 2mk2 sak + 2mk sak + mk2 + mk ⎤ ⎦ (B.18) Again, for a = a = γ the end-point contribution cannot arise as no corresponding splitting is integrated over, K γt γ,k = 0. Finally, the K˜ terms, containing correlations between both initial state partons, read K˜ aa (x ) = Praeag (x ) log(1 − x ) + Qa2 δaa log(1 − x ) 2 1 − x + − 3 δ(1 − x ) . (B.19) They thus vanish for a = a = γ . Finally, the P operator is given by Paa (x , μ2F ) α Paa (x ) = 2π k cf. Eq. (2.10). Therein, the regularised Altarelli-Parisi splitting functions are given by Paa (x ) = Praeag (x )+δaa The missing function JγaQ (x , y) is given in Eq. (5.58) of [ 41 ]. In the massless limit this reduces to Ki,aa (x ) = −δaa γi Q2 i + + δ(1 − x ) (B.16) Most importantly, Pγ γ (x ) = γγ δ(1 − x ). (B.21) (B.22) C Precise forms of AiIk and AaKb C.1 The AiIk term TheexplicitdependenceoftheI operatoronthe{αdip}parameters, or more precisely only on αFF, is given by [ 54,55,86– 88 ] if i is a massless fermion AiIk({α}) = −ln2 αFF − γi (lnαFF + 1 − αFF) mk = 0 2 1 − y+ + 2xy+ = ln2 (1 − x + y+)(1 + x − y+) − 2ln2 1 + y+ − x 1 + y+ + 21 (μ1i2−+α2FαFyF+Fy)+μdi2 − 2lnαFF +Li2 xx−− −+ ax − Li2 x− + a x− +Li2 xx++ −− bx − Li2 x+ − b x+ −Li2 cc +− ax + Li2 c − a c +ln(c + x)ln (a − c)(x+ − x) (a + x)(x+ + c) (a − c)x+ −lnc ln a(x+ + c) (a + b)(x− − x) −ln(b − x)ln (a + x)(x− − b) (a + b)x− +lnb ln a(x− − b) −ln((a + x)(b − x+))ln(x+ − x) +ln(a(b − x+))ln x+ (a + x)x+x− +lnd ln a(x+ − x)(x− − x) +ln x−x− x ln c + x− − a + x− +21 ln a +a x ln a(a + x)(a + x+)2 ln a(a + x )(a + x+)2 2 (1 − μk )2 − μi 1 − μi2 − μk2 1 if mk = 0 1−μk if mi = 0 1+μk (C.5) (C.7) (C.8) (C.9) y = 21 1 − μk2 + α FF(1 − μk )2 − (1 − μk ) × 2 (α FF(1 − μk ))2 + (1 + μk )2 − 2α FF(1 + μk ) " 2 (1 − α FF) 1 − α FF y+ − The superscripts indicate the general functional form of the contribution. In case of the α FI and α IF depedence the existing plus distributions are modified in the K t terms as follows (..)+ → (..)1−α FI with 0 1 1 1−α FI 2 μi/k = y+ = mi2/k qi2k As can be seen, the massless case can be recovered from the massive one in the limit that mk → 0. The AiIk term could in principle also be shifted into the K operator to leave the I operator invariant. To be applicable to pure final-state calculations, though, it is left where it is. The various abbreviations above are defined as (C.4) dx f (x )[g(x )]1−α FI = dx g(x )[ f (x ) − f (1)]. As the singular point sits at x = 1 this modification amouts to a finite contribution. The α IF on the other hand leads to finite regular functional contribution in the AaKbt term given in [ 87,89 ]. Lastly, the α II dependence reads [ 55 ] AaK˜b(x , α II) = log α II , 1 − x except when a = b = γ , in which case it is zero as the original K˜ γ γ term. 1. D.R. Yennie , S.C. Frautschi , H. Suura , The infrared divergence phenomena and high-energy processes . Ann. Phys . 13 , 379 - 452 ( 1961 ) 2. V.V. Sudakov , Vertex parts at very high-energies in quantum electrodynamics . Sov. Phys. JETP 3 , 65 - 71 ( 1956 ). [Zh. Eksp. Teor. Fiz . 30 , 87 ( 1956 )] 3. W. Beenakker , A. Denner , S. Dittmaier , R. Mertig , T. Sack, Highenergy approximation for on-shell W pair production . Nucl. Phys. B 410 , 245 - 279 ( 1993 ). ( CERN-TH- 6832-93) 4. W. Beenakker , A. Denner , S. Dittmaier , R. Mertig , On shell W pair production in the TeV range . Phys. Lett. B 317 , 622 - 630 ( 1993 ). ( CERN-TH- 6905-93) 5. W. Beenakker , A. Denner , W. Hollik , R. Mertig , T. Sack , D. Wackeroth , Electroweak one loop contributions to top pair production in hadron colliders . Nucl. Phys. B 411 , 343 - 380 ( 1994 ). (MPI-PH-93-20) 6. A. Denner , S. Dittmaier , T. Hahn, Radiative corrections to Z Z → Z Z in the electroweak standard model . Phys. Rev. D 56 , 117 - 134 ( 1997 ). arXiv:hep-ph/9612390 7. V.S. Fadin , L.N. Lipatov , A.D. Martin , M. Melles , Resummation of double logarithms in electroweak high-energy processes . Phys. Rev. D 61 , 094002 ( 2000 ). arXiv:hep-ph/9910338 8. J.H. Kühn , A.A. Penin , V.A. Smirnov , Summing up subleading Sudakov logarithms . Eur. Phys. J. C 17 , 97 - 105 ( 2000 ). arXiv:hep-ph/9912503 9. M. Ciafaloni , P. Ciafaloni , D. Comelli , Bloch-Nordsieck violating electroweak corrections to inclusive TeV scale hard processes . Phys. Rev. Lett . 84 , 4810 - 4813 ( 2000 ). arXiv:hep-ph/0001142 10. M. Hori , H. Kawamura , J. Kodaira , Electroweak Sudakov at two loop level . Phys. Lett. B 491 , 275 - 279 ( 2000 ). arXiv:hep-ph/0007329 11. A. Denner , S. Pozzorini , One loop leading logarithms in electroweak radiative corrections. 1 . Results . Eur. Phys. J. C 18 , 461 - 480 ( 2001 ). arXiv:hep-ph/0010201 12. A. Denner , S. Pozzorini , One loop leading logarithms in electroweak radiative corrections. 2. Factorization of collinear singularities . Eur. Phys. J. C 21 , 63 - 79 ( 2001 ). arXiv:hep-ph/0104127 13. M. Melles , Resummation of angular dependent corrections in spontaneously broken gauge theories . Eur. Phys. J. C 24 , 193 - 204 ( 2002 ). arXiv:hep-ph/0108221 14. W. Beenakker , A. Werthenbach , Electroweak two loop Sudakov logarithms for on-shell fermions and bosons . Nucl. Phys. B 630 , 3 - 54 ( 2002 ). arXiv:hep-ph/0112030 15. B. Feucht , J.H. Kühn , A.A. Penin , V.A. Smirnov , Two loop Sudakov form-factor in a theory with mass gap . Phys. Rev. Lett . 93 , 101802 ( 2004 ). arXiv:hep-ph/0404082 16. B. Jantzen , J.H. Kühn , A.A. Penin , V.A. Smirnov , Two-loop electroweak logarithms . Phys. Rev. D 72 , 051301 ( 2005 ). arXiv:hep-ph/0504111. [Erratum: Phys. Rev. D 74 , 019901 ( 2006 )] 17. U. Baur, Weak boson emission in hadron collider processes . Phys. Rev. D 75 , 013005 ( 2007 ). arXiv:hep-ph/0611241 18. A. Denner , L. Hofer , A. Scharf , S. Uccirati , Electroweak corrections to lepton pair production in association with two hard jets at the LHC . JHEP 01 , 094 ( 2015 ). arXiv: 1411 .0916 [hep-ph] 19. S. Kallweit , J.M. Lindert , P. Maierhöfer , S. Pozzorini , M. Schönherr, NLO electroweak automation and precise predictions for W+multijet production at the LHC . JHEP 04 , 012 ( 2015 ). arXiv: 1412 .5157 [hep-ph] 20. A. Denner , S. Dittmaier , M. Hecht , C. Pasold, NLO QCD and electroweak corrections to W + γ production with leptonic W-boson decays . JHEP 04 , 018 ( 2015 ). arXiv: 1412 .7421 [hep-ph] 21. S. Kallweit , J.M. Lindert , S. Pozzorini , M. Schönherr , P. Maierhöfer, NLO QCD+ EW automation and precise predictions for V+multijet production . Proceedings , 50th Rencontres de Moriond, QCD and high energy interactions: La Thuile , Italy, March 21 -28, 2015 , pp. 121 - 124 ( 2015 ) 22. S. Frixione , V. Hirschi , D. Pagani , H.S. Shao , M. Zaro , Electroweak and QCD corrections to top-pair hadroproduction in association with heavy bosons . JHEP 06 , 184 ( 2015 ). arXiv: 1504 .03446 [hepph] 23. M. Chiesa , N. Greiner , F. Tramontano , Automation of electroweak corrections for LHC processes . J. Phys. G 43 ( 1 ), 013002 ( 2016 ). arXiv: 1507 .08579 [hep-ph] 24. A. Denner , S. Dittmaier , M. Hecht , C. Pasold, NLO QCD and electroweak corrections to Z + γ production with leptonic Z-boson decays . JHEP 02 , 057 ( 2016 ). arXiv: 1510 .08742 [hep-ph] 25. S. Kallweit , J.M. Lindert , P. Maierhöfer , S. Pozzorini , M. Schönherr , NLO QCD+ EW predictions for V + jets including off-shell vector-boson decays and multijet merging . JHEP 04 , 021 ( 2016 ). arXiv: 1511 .08692 [hep-ph] 26. B. Biedermann , A. Denner , S. Dittmaier , L. Hofer , B. Jäger , Electroweak corrections to pp → μ+μ−e+e− + X at the LHC: a Higgs background study . Phys. Rev. Lett . 116 ( 16 ), 161803 ( 2016 ). arXiv: 1601 .07787 [hep-ph] 27. B. Biedermann , M. Billoni , A. Denner , S. Dittmaier , L. Hofer , B. Jäger , L. Salfelder, Next-to-leading-order electroweak corrections to pp → W +W − → 4 leptons at the LHC . JHEP 06 , 065 ( 2016 ). arXiv: 1605 .03419 [hep-ph] 28. A. Denner , M. Pellen, NLO electroweak corrections to off-shell top-antitop production with leptonic decays at the LHC . JHEP 08 , 155 ( 2016 ). arXiv: 1607 .05571 [hep-ph] 29. R. Frederix , S. Frixione , V. Hirschi , D. Pagani , H.-S. Shao, M. Zaro , The complete NLO corrections to dijet hadroproduction . JHEP 04 , 076 ( 2017 ). arXiv: 1612 .06548 [hep-ph] 30. A. Denner , J.-N. Lang , M. Pellen , S. Uccirati , Higgs production in association with off-shell top-antitop pairs at NLO EW and QCD at the LHC . JHEP 02 , 053 ( 2017 ). arXiv: 1612 .07138 [hep-ph] 31. S. Kallweit , J.M. Lindert , S. Pozzorini , M. Schönherr , NLO QCD+ EW predictions for 2 2ν diboson signatures at the LHC . JHEP 11 , 120 ( 2017 ). arXiv: 1705 .00598 [hep-ph] 32. J .M. Lindert et al., Precise predictions for V + jets dark matter backgrounds . Eur. Phys. J. C 77 ( 12 ), 829 ( 2017 ). arXiv: 1705 .04664 [hep-ph] 33. F. Granata , J.M. Lindert , C. Oleari , S. Pozzorini, NLO QCD+ EW predictions for HV and HV+jet production including partonshower effects . JHEP 09 , 012 ( 2017 ). arXiv: 1706 .03522 [hep-ph] 34. M. Chiesa , N. Greiner , M. Schönherr , F. Tramontano , Electroweak corrections to diphoton plus jets . JHEP 10 , 181 ( 2017 ). arXiv: 1706 .09022 [hep-ph] 35. B. Biedermann , A. Denner , M. Pellen , Complete NLO corrections to W+W+ scattering and its irreducible background at the LHC . JHEP 10 , 124 ( 2017 ). arXiv: 1708 .00268 [hep-ph] 36. B. Biedermann , A. Denner , L. Hofer, Next-to-leading-order electroweak corrections to the production of three charged leptons plus missing energy at the LHC . JHEP 10 , 043 ( 2017 ). arXiv: 1708 .06938 [hep-ph] 37. N. Greiner , M. Schönherr , NLO QCD+ EW corrections to diphoton production in association with a vector boson . arXiv:1710 .11514 [hep-ph] 38. R. Frederix , D. Pagani , M. Zaro , Large NLO corrections in t t¯W ± and t t¯ t t¯ hadroproduction from supposedly subleading EW contributions . arXiv:1711 .02116 [hep-ph] 39. S. Frixione , Z. Kunszt , A. Signer , Three jet cross-sections to next-to-leading order . Nucl. Phys. B 467 , 399 - 442 ( 1996 ). arXiv:hep-ph/9512328 40. S. Catani , M.H. Seymour , A general algorithm for calculating jet cross-sections in NLO QCD . Nucl. Phys. B 485 , 291 - 419 ( 1997 ). arXiv:hep-ph/9605323. [Erratum: Nucl. Phys. B 510 , 503 ( 1998 )] 41. S. Catani , S. Dittmaier , M.H. Seymour , Z. Trocsanyi , The dipole formalism for next-to-leading order QCD calculations with massive partons . Nucl. Phys. B 627 , 189 - 265 ( 2002 ). arXiv:hep-ph/0201036 42. D.A. Kosower , Antenna factorization of gauge theory amplitudes . Phys. Rev. D 57 , 5410 - 5416 ( 1998 ). arXiv:hep-ph/9710213 43. Z. Nagy , D.E. Soper , General subtraction method for numerical calculation of one loop QCD matrix elements . JHEP 09 , 055 ( 2003 ). arXiv:hep-ph/0308127 44. S. Dittmaier , A general approach to photon radiation off fermions . Nucl. Phys. B 565 , 69 - 122 ( 2000 ). arXiv:hep-ph/9904440 45. S. Dittmaier , A. Kabelschacht , T. Kasprzik, Polarized QED splittings of massive fermions and dipole subtraction for noncollinear-safe observables . Nucl. Phys. B 800 , 146 - 189 ( 2008 ). arXiv: 0802 .1405 [hep-ph] 46. T. Gehrmann, N. Greiner , Photon radiation with MadDipole . JHEP 12 , 050 ( 2010 ). arXiv: 1011 .0321 [hep-ph] 47. F. Krauss , R. Kuhn , G. Soff, AMEGIC++ 1.0: a matrix element generator in C++ . JHEP 02 , 044 ( 2002 ). arXiv:hep-ph/0109036 48. T. Gleisberg , S. Höche , F. Krauss , M. Schönherr , S. Schumann , F. Siegert , J. Winter , Event generation with SHERPA 1.1 . JHEP 02, 007 ( 2009 ). arXiv: 0811 .4622 [hep-ph] 49. T. Gleisberg , F. Krauss , Automating dipole subtraction for QCD NLO calculations . Eur. Phys. J. C 53 , 501 - 523 ( 2008 ). arXiv: 0709 .2881 [hep-ph] 50. J. Archibald , An automated NLO framework in SHERPA. Ph.D. thesis , Durham U. ( 2011 ) 51. J.R. Andersen et al., Les Houches 2015: physics at TeV colliders Standard Model Working Group Report. 9th Les Houches Workshop on Physics at TeV Colliders (PhysTeV 2015 ) Les Houches , France, June 1-19, 2015 ( 2016 ) 52. S. Alioli et al., Precision studies of observables in pp → W → lνl and pp → γ , Z → l+l − processes at the LHC . Eur. Phys. J. C 77 ( 5 ), 280 ( 2017 ). arXiv: 1606 .02330 [hep-ph] 53. B. Biedermann , S. Bruer , A. Denner , M. Pellen , S. Schumann , J.M. Thompson , Automation of NLO QCD and EW corrections with Sherpa and Recola . Eur. Phys. J. C 77 , 492 ( 2017 ). arXiv: 1704 .05783 [hep-ph] 54. Z. Nagy , Z. Trocsanyi , Next-to-leading order calculation of four jet observables in electron positron annihilation . Phys. Rev. D 59 , 014020 ( 1999 ). arXiv:hep-ph/9806317. [Erratum: Phys. Rev. D 62 , 099902 ( 2000 )] 55. Z. Nagy , Next-to-leading order calculation of three jet observables in hadron hadron collision . Phys. Rev. D 68 , 094002 ( 2003 ). arXiv:hep-ph/0307268 56. E. Bothmann, M. Schönherr , S. Schumann , Reweighting QCD matrix -element and parton-shower calculations . Eur. Phys. J. C 76 ( 11 ), 590 ( 2016 ). arXiv: 1606 .08753 [hep-ph] 57. F. Bloch , A. Nordsieck, Note on the radiation field of the electron . Phys. Rev . 52 , 54 - 59 ( 1937 ). (RX-1199) 58. T. Kinoshita, Mass singularities of Feynman amplitudes . J. Math. Phys. 3 , 650 - 677 ( 1962 ) 59. T.D. Lee , M. Nauenberg , Degenerate Systems and Mass Singularities. Phys. Rev . 133 , B1549 - B1562 ( 1964 ). [ 25 ( 1964 )] 60. L.A. Harland-Lang , V.A. Khoze , M.G. Ryskin , Sudakov effects in photon-initiated processes . Phys. Lett. B 761 , 20 - 24 ( 2016 ). arXiv: 1605 .04935 [hep-ph] 61. M. Cacciari , G.P. Salam , G. Soyez, FastJet user manual . Eur. Phys. J. C 72 , 1896 ( 2012 ). arXiv: 1111 .6097 [hep-ph] 62. S. Frixione, Isolated photons in perturbative QCD . Phys. Lett. B 429 , 369 - 374 ( 1998 ). arXiv:hep-ph/9801442 63. A.D. Martin , R.G. Roberts , W.J. Stirling , R.S. Thorne , Parton distributions incorporating QED contributions . Eur. Phys. J. C 39 , 155 - 161 ( 2005 ). arXiv:hep-ph/0411040 64. C. Schmidt , J. Pumplin , D. Stump , C.P. Yuan , CT14QED parton distribution functions from isolated photon production in deep inelastic scattering . Phys. Rev. D 93 ( 11 ), 114015 ( 2016 ). arXiv: 1509 .02905 [hep-ph] 65. A. Manohar , P. Nason , G.P. Salam , G. Zanderighi, How bright is the proton? A precise determination of the photon parton distribution function . Phys. Rev. Lett . 117 ( 24 ), 242002 ( 2016 ). arXiv: 1607 .04266 [hep-ph] 66. A.V. Manohar , P. Nason , G.P. Salam , G. Zanderighi, The photon content of the proton . JHEP 12 , 046 ( 2017 ). arXiv: 1708 .01256 [hep-ph] 67. L .A. Harland-Lang , V .A. Khoze , M .G. Ryskin, The photon PDF in events with rapidity gaps . Eur. Phys. J. C 76 ( 5 ), 255 ( 2016 ). arXiv: 1601 .03772 [hep-ph] 68. R.D. Ball et al., NNPDF , Parton distributions for the LHC Run II , JHEP 04 , 040 ( 2015 ). arXiv: 1410 .8849 [hep-ph] 69. V. Bertone , S. Carrazza , N.P. Hartland , J. Rojo , Illuminating the photon content of the proton within a global PDF analysis . arXiv:1712 .07053 [hep-ph] 70. R.D. Ball , V. Bertone , S. Carrazza , L. Del Debbio , S. Forte , A. Guffanti , N.P. Hartland , J. Rojo, NNPDF , Parton distributions with QED corrections . Nucl. Phys. B 877 , 290 - 320 ( 2013 ). arXiv: 1308 .0598 [hep-ph] 71. M. Cacciari , M. Greco , P. Nason , The P( T) spectrum in heavy flavor hadroproduction . JHEP 05 , 007 ( 1998 ). arXiv:hep-ph/9803400 72. K.A. Olive et al., Particle Data Group, Review of particle physics. Chin. Phys. C 38 , 090001 ( 2014 ) 73. J.R. Andersen et al., LHC Higgs Cross Section Working Group, Handbook of LHC Higgs cross sections: 3. Higgs properties . arXiv: 1307 .1347 [hep-ph] 74. Munich is the abbreviation of “MUlti-chaNnel Integrator at Swiss (CH) precision”-an parton level NLO generator by S. Kallweit . In preparation 75. J. Alwall , R. Frederix , S. Frixione , V. Hirschi , F. Maltoni , O. Mattelaer , H.S. Shao , T. Stelzer , P. Torrielli , M. Zaro , The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations . JHEP 07 , 079 ( 2014 ). arXiv: 1405 .0301 [hep-ph] 76. S. Actis , A. Denner , L. Hofer , A. Scharf , S. Uccirati , Recursive generation of one-loop amplitudes in the standard model . JHEP 04 , 037 ( 2013 ). arXiv: 1211 .6316 [hep-ph] 77. A. Denner , S. Dittmaier , M. Roth , L.H. Wieders , Electroweak corrections to charged-current e+e− → 4 fermion processes: technical details and further results . Nucl. Phys. B 724 , 247 - 294 ( 2005 ). arXiv:hep-ph/0505042. [Erratum: Nucl. Phys. B 854 , 504 ( 2012 )] 78. S. Dulat, T.-J. Hou , J. Gao , M. Guzzi , J. Huston , P. Nadolsky , J. Pumplin , C. Schmidt , D. Stump , C.P. Yuan , New parton distribution functions from a global analysis of quantum chromodynamics . Phys. Rev. D 93 ( 3 ), 033006 ( 2016 ). arXiv: 1506 .07443 [hep-ph] 79. A. Buckley , J. Ferrando , S. Lloyd , K. Nordström , B. Page , M. Rüfenacht , M. Schönherr , G. Watt, LHAPDF6: parton density access in the LHC precision era . Eur. Phys. J. C 75 ( 3 ), 132 ( 2015 ). arXiv: 1412 .7420 [hep-ph] 80. M. Cacciari , G.P. Salam , G. Soyez, The anti-k(t ) jet clustering algorithm . JHEP 04 , 063 ( 2008 ). arXiv: 0802 .1189 [hep-ph] 81. L. Basso , S. Dittmaier , A. Huss , L. Oggero , Techniques for the treatment of IR divergences in decay processes at NLO and application to the top-quark decay . Eur. Phys. J. C 76 ( 2 ), 56 ( 2016 ). arXiv: 1507 .04676 [hep-ph] 82. F. Krauss , D. Napoletano , Towards a fully massive five-flavour scheme . arXiv:1712 .06832 [hep-ph] 83. M. Schönherr , F. Krauss , Soft photon radiation in particle decays in SHERPA . JHEP 12 , 018 ( 2008 ). arXiv: 0810 .5071 [hep-ph] 84. T. Gleisberg , S. Höche , Comix, a new matrix element generator . JHEP 12 , 039 ( 2008 ). arXiv: 0808 .3674 [hep-ph] 85. A. Buckley , J. Butterworth , L. Lönnblad , D. Grellscheid , H. Hoeth , J. Monk , H. Schulz , F. Siegert , Rivet user manual . Comput. Phys. Commun . 184 , 2803 - 2819 ( 2013 ). arXiv: 1003 .0694 [hep-ph] 86. R. Frederix , T. Gehrmann , N. Greiner , Integrated dipoles with MadDipole in the MadGraph framework . JHEP 06 , 086 ( 2010 ). arXiv: 1004 .2905 [hep-ph] 87. J.M. Campbell , R.K. Ellis , F. Tramontano , Single top production and decay at next-to-leading order . Phys. Rev. D 70 , 094012 ( 2004 ). arXiv:hep-ph/0408158 88. G. Bevilacqua, M. Czakon , C.G. Papadopoulos , R. Pittau , M. Worek, Assault on the NLO Wishlist: pp → t anti-t b anti-b . JHEP 09 , 109 ( 2009 ). arXiv: 0907 .4723 [hep-ph] 89. J.M. Campbell , F. Tramontano , Next-to-leading order corrections to Wt production and decay . Nucl. Phys. B 726 , 109 - 130 ( 2005 ). arXiv:hep-ph/0506289


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-018-5600-z.pdf

Marek Schönherr. An automated subtraction of NLO EW infrared divergences, The European Physical Journal C, 2018, 119, DOI: 10.1140/epjc/s10052-018-5600-z