Automated evaluation of electroweak Sudakov logarithms in Sherpa

We present an automated implementation for the calculation of one-loop double and single Sudakov logarithms stemming from electroweak radiative corrections within the Sherpa event generation framework, based on the derivation in[1]. At high energies, these logarithms constitute the leading contributions to the full NLO electroweak corrections. As examples, we show applications for relevant processes at both the LHC and future hadron colliders, namely on-shell W boson pair production, EW-induced dijet production and electron-positron production in association with four jets, providing the first estimate of EW corrections at this multiplicity.


Introduction
In perturbative electroweak (EW) theory, higher-order corrections are given by emissions of virtual and real gauge bosons and are known to have a large effect in the hard tail of observables at the LHC and future colliders [2]. In contrast to massless gauge theories, where real emission terms are necessary to regulate infrared divergences, in EW theory the weak gauge bosons are massive and therefore provide a natural lower scale cut-off, such that finite logarithms appear instead. Moreover, since the emission of an additional weak boson leads to an experimentally different signature with respect to the Born process, virtual logarithms are of physical significance without the inclusion of real radiation terms. The structure of such logarithmic contributions, referred to as Sudakov logarithms [3], and their factorisation properties were derived in full generality by Denner and Pozzorini [1,4,5,6] at leading and next-to-leading logarithmic accuracy at both one-and two-loop order in the EW coupling expansion. In particular, they have shown that in the high-energy scattering limit, at least at one-loop order, these logarithmic corrections can be factored as a sum over pairs of external legs in an otherwise process-independent way, hence providing a straightforward algorithm for computing them for any given process. The high-energy limit requires all invariants formed by pairs of particles to be large compared to the scale set by the weak gauge boson masses. In this sense we expect effects due to electroweak corrections in general, and to Sudakov logarithms in particular, to give rise to large effects in the hard tail of observables that have the dimension of energy, such as for example the transverse momentum of a final state particle. In this work we present an implementation for computing EW Sudakov corrections in the form they are presented in [1] within the SHERPA event generator framework [7]. While Sudakov corrections have been SL contributions coming from a diagram of the type shown in Fig. 1. The former correspond to the case of the virtual gauge boson V becoming soft and collinear to one of the external legs i or j of the matrix element M, i.e. when |r ij | ≡ |p i + p j | 2 m 2 W . Similarly, as |p i + p j | 2 ∼ 2 p i · p j , this limit also encodes the case of a collinear, non-soft (and vice-versa) gauge boson which instead gives rise to SL contributions. In turn, these must also be included, since they can numerically be of a similar size as the DL one, see the discussion and the detailed derivation in [1,5]. The reason for this can also be understood intuitively: while the diagram in Fig. 1 is the only source of double logarithms there are several contributions giving rise to single logarithms, which we discuss later in this section. In all cases, following the conventions in [1], we write DL and SL terms in the following way, respectively: (2.1) Note that the use of the W mass as a reference scale is only a convention and the same form is indeed also used for contributions coming from Z bosons and photons. Remainder logarithms containing ratios of m W and m Z or the vanishing photon mass are taken care of by additional logarithmic terms, as explained in more details below. After having evaluated all relevant diagrams, and using Ward identities, one can show that the DL and SL contributions factorise off the Born matrix element M, and can be written in terms of a sum over pairs of external legs only. Hence, we define these corrections as a multiplicative factor K NLL to the squared matrix element |M(Φ)| 2 at a given phase space point Φ: where we have written M h ≡ M h (Φ) and the sums go over helicity configurations h and the different types of logarithmic contributions c. Following the original notation, these contributions are divided into Leading and Subleading Soft and Collinear (LSC/SSC) logarithms, Collinear or Soft (C) logarithms and those arising from Parameter Renormalisation (PR). More details about their respective definitions are given below. Since δ c is in general a tensor in SU(2) × U(1), the δ c M h are usually not proportional to the original matrix element M h . The tensor structure come from the fact that in general the various δ c are proportional to EW vertices, which in turn means that a single leg or pairs of legs can get replaced with particles of a different kind, according to the relevant Feynman rules. As in [1], we denote the leg index i as a shorthand for the external fields φ i . Denoting with {i} the set of all external fields, we therefore have h . In our implementation, the construction and evaluation of these additional amplitudes is taken care of by an interface to the automated tree-level COMIX matrix element generator [13], which is available within the SHERPA framework. Before evaluating such amplitudes, the energy is re-distributed among the particles to put all (replaced) external legs on-shell again. The required auxiliary processes are automatically set up during the initialisation of the framework. Since the construction of additional processes can be computationally non-trivial, we have taken care of reducing their number by re-using processes that are shared between different contributions c, and by using crossing relations when only the order of external fields differs. In our implementation we consistently neglect purely electromagnetic logarithms which arise from the gap between zero (or the fictitious photon mass) and m W . In SHERPA such logarithms can be treated in either of two ways. First, one can compute the purely electromagnetic NLO correction to the desired process, consisting of both virtual and real photon emission, which gives the desired logarithm at fixed order. Alternatively, one can resum soft-photon emission to all orders, which in SHERPA is achieved through the formalism of Yennie, Frautschi and Suura (YFS) [33] or by a QED parton shower, as is discussed e.g. in [35]. In contrast, logarithms originating from the running between the W mass and the Z mass are included, as we explain in the next paragraphs, where we discuss the individual contributions δ c .
Leading Soft-Collinear logarithms: LSC The leading corrections are given by the soft and collinear limit of a virtual gauge boson and are proportional to L(|r kl |). By writing and neglecting the last term on the right-hand side, which is of order O(1), one has split the soft-collinear contribution into a leading soft-collinear (LSC) and a angular-dependent subleading soft-collinear (SSC) one, corresponding to the first and second term on the right-hand side, respectively. The full LSC correction is then given by where the sum runs over the n external legs, and the coefficient matrix on the right-hand side is given by C ew and I Z are the electroweak Casimir operator and the Z gauge coupling, respectively. The second term δ Z appears as an artefact of writing L(s) and l(s) in terms of the W mass even for Z boson loops, and hence the inclusion of this term takes care of the gap between the Z and the W mass at NLL accuracy. We denote the remaining terms with the superscript "LSC". Note that this term is in general non-diagonal, since C ew mixes photons and Z bosons. As is explicit in Eq. (2.2), coefficients need to be computed per helicity configuration; in the case of longitudinally polarised vector boson appearing as external particles, they need to be replaced with the corresponding Goldstone bosons using the Goldstone Boson Equivalence Theorem. This is a consequence of the scaling of longitudinal polarisation vectors with the mass of the particle, which prevents the direct application of the eikonal approximation used to evaluate the high-energy limit, as is detailed in [1]. To calculate such matrix elements we have extended the default Standard Model implementation in SHERPA with a full Standard Model including Goldstone bosons generated through the UFO interface of SHERPA [36,37]. We have tested this implementation thoroughly against RECOLA [38]. ‡ Subleading Soft-Collinear logarithms: SSC The second term in Eq. 2.3 gives rise to the angulardependent subleading part of the corrections from soft-collinear gauge-boson loops. It can be written as a sum over pairs of external particles: where the coefficient matrices for the different gauge bosons V a = A, Z, W ± are Note that while the photon couplings I A are diagonal, the eigenvalues for I Z and I ± ≡ I W ± can be nondiagonal, leading again to replacements of external legs as described in the LSC case. ‡ Note that in particular this adds the possibility for SHERPA users to compute Goldstone bosons contributions to any desired process, independently of whether this is done in the context of calculating EW Sudakov corrections.
Collinear or soft single logarithms: C The two sources that provide either collinear or soft logarithms are the running of the field renormalisation constants, and the collinear limit of the loop diagrams where one external leg splits into two internal lines, one of which being a vector boson V a . The ensuing correction factor can be written as a sum over external legs: For chiral fermions f , the coefficient matrix is given by We label the chirality of the fermion by κ, which can take either the value L (left-handed) or R (righthanded). The label σ specifies the isospin, its values ± refer to up-type quarks/neutrinos or down-type quarks/leptons, respectively. The sine (cosine) of the Weinberg angle is denoted by s w (c w ). Note that we further subdivide the collinear contributions for fermion legs into "Yukawa" terms formed by terms proportional to the ratio of the fermion masses over the W mass, which we denote by the superscript "Yuk", and the remaining collinear contributions denoted by the superscript "C". These Yukawa terms only appear for external fermions, such that for all other external particles ϕ, we have δ C ϕ ϕ = δ C ϕ ϕ . For charged and neutral transverse gauge bosons, must be used, where the b ew are combinations of Dynkin operators that are proportional to the one-loop coefficients of the β-function for the running of the gauge-boson self-energies and mixing energies. Their values are given in terms of s w and c w in [1]. Longitudinally polarised vector bosons are again replaced with Goldstone bosons. When using the matrix element on the right-hand side of Eq. (2.7) in the physical EW phase, the following (diagonal) coefficient matrices must be used for charged and neutral longitudinal gauge bosons: 11) where N C = 3 is the number of colour charges.
Parameter Renormalisation logarithms: PR The last contribution we consider is the one coming from the renormalisation of EW parameters, such as boson and fermion masses and the QED coupling α, and all derived quantities. The way we extract these terms, is by running all EW parameters up to a given scale, µ EW , which in all cases presented here corresponds to the partonic scattering centre of mass energy, and re-evaluate the matrix element value with these evolved parameters. We then take the ratio of this 'High Energy' (HE) matrix element with respect to the original value, such that, calling {p ew } the complete set of EW parameters, The evolution of each EW parameter is obtained through 13) and the exact expressions for δp ew can be found in Eqs. (5.4-5.22) of [1]. The right hand side of Eq. (2.12) corresponds to the original derivation of Denner and Pozzorini, while the left hand side is the actual implementation we have used. The two differ by terms that are formally of higher order, (α log µ 2 EW /m 2 W ) 2 . In fact, although they are logarithmically enhanced, they are suppressed by an additional power of α with respect to the leading terms considered here, α log 2 µ 2 EW /m 2 W .
Generating event samples with EW Sudakov logarithmic contributions After having discussed the individual contributions c, we can return to Eq. (2.2), for which we now have all the ingredients to evaluate it for an event with the phase-space point Φ. Defining the relative logarithmic contributions ∆ c , we can rewrite it as (2.14) In the event sample, the relative contributions ∆ c are given in the form of named HEPMC weights [14] (details on the naming will be given in the manual of the upcoming SHERPA release). This is done to leave the user freedom on how to combine such weights with the main event weight.
In the context of results, Sec. 3, we employ these corrections in either of two ways. One way is to include them at fixed order, where B is the Born contribution. The alternative is to exploit the fact that Sudakov EW logarithms exponentiate (see Refs. [22]- [31]) to construct a resummed fully differential cross section following the approach discussed in [32].
Matching to higher orders and parton showers SHERPA internally provides the possibility to the user to obtain NLO corrections, both of QCD [39] and EW [15] origin, and to further generate fully showered [40,41] and hadronised events [7]. In addition one is able to merge samples with higher multiplicities in QCD through the CKKW [42], or the MEPS@NLO [43] algorithms. In all of the above cases (except the NLO EW), the corrections implemented here can be simply applied using the K factor methods in SHERPA to one or all the desired processes of the calculations, as there is no double counting between EW Sudakov corrections and pure QCD ones. A similar reasoning can be applied for the combination with a pure QCD parton shower. Although we do not report the result of this additional check here, we have tested the combination of the EW Sudakov corrections with the default shower of SHERPA. Technical checks and physics applications for the combination with matching and merging schemes, and for the combination with QED logarithms, are left for a future publication. If one aims at matching Sudakov logarithms to higher-order EW effects, such as for example combining fixed order NLO EW results with a resummed NLL Sudakov correction, it is for now required for the user to manually do the subtraction of double counting that one encounters in these cases. An automation of such a scheme is also outside the scope of this publication, and will be explored in the future.

Results
Before discussing our physics applications, we report a number of exact comparisons to other existing calculations of NLL EW Sudakov logarithms. A subset of the results used for this comparison for pp → V j processes is shown in App. A. They all agree with reference ones on a sub-percent level over the entire probed transverse momentum range from p T,V = 100 GeV to 2000 GeV. In addition to this, the implementation has been guided by passing a number of tests based on a direct comparison of tabulated numerical values for each contribution discussed in Sec. 2 in the high-energy limit, that are given for several electron-positron collision processes in [1] and for the pp → W Z and pp → W γ processes in [44]. Our final implementation passes all these tests.
In the remainder of this section, we present a selection of physics results obtained using our implementation and where possible show comparisons with existing alternative calculations. The aim is twofold, first we wish to show the variety of processes that can be computed with our implementation, and second we want to compare to existing alternative methods to obtain EW corrections to further study the quality of the approximation. In particular, where available, we compare our NLL predictions to full NLO EW corrections, as well as to the EW virt approximation defined in [21]. It is important to note that we only consider EW corrections to purely EW processes here, such that there are no subleading Born corrections that might otherwise complicate the comparison to a full EW NLO calculation or to the EW virt approximation. We consider diboson production, dijet production, and Z production in association with 4 jets. In all cases, we focus on the discussion of (large) transverse momentum distributions, since they are directly sensitive to the growing effect of the logarithmic contributions when approaching the high-energy limit.
The logarithmic contributions are applied to parton-level LO calculations provided by the COMIX matrix element generator implemented in SHERPA. The parton shower and non-perturbative parts of the simulation are disabled, including the simulation of multiple interactions and beam remnants, and the hadronisation. Also higher-order QED corrections to the matrix element and by YFS-type resummation are turned off. The contributions to the NLL corrections as defined in Sec. 2 are calculated individually and combined aposteriori, such that we can study their effects both individually and combined. We consider the fixed-order and the resummed option for the combination, as detailed in Eqs. (2.15) and (2.16), respectively.
For the analysis, we use the Rivet 2 [45] framework, and events are passed to analysis using the HEPMC event record library. Unless otherwise specified, simulations are obtained using the NNPDF3.1, next-tonext-to-leading order PDF set [46], while in processes where we include photon initiated processes we instead use the NNPDF3.1 LUX PDF set [47]. In all cases, PDFs are obtained through the LHAPDF [48] interface implemented in SHERPA. When jets appear in the final state, we cluster them using the anti-k T algorithm [49] with a jet radius parameter of R = 0.4, through an interface to FASTJET [50]. The CKM matrix in our calculation is equal to the unit matrix, i.e. no mixing of the quark generations is allowed. Electroweak parameters are determined using tree-level relations using a QED coupling value of α(m Z ) = 1/128.802 and the following set of masses and decay widths, if not explicitly mentioned otherwise: Note that α is not running in the nominal calculation as running effects are all accounted for by the PR contributions.

W W production in pp collisions at 13 and 100 TeV
Our first application is the calculation of EW Sudakov effects in on-shell W boson pair production at hadron colliders, which has lately been experimentally probed at the LHC [51,52], e.g. to search for anomalous gauge couplings. We compare the Sudakov EW approximation to both the full NLO EW calculation and the EW virt approximation. The latter has been applied to W W production in [35,53]. In addition to that, EW corrections for W W production have also been studied in [54,55,56,57,58,59], and NNLL EW Sudakov corrections have been calculated in [54]. We have performed this study at current LHC energies, √ s = 13 TeV, as well as at a possible future hadron collider with √ s = 100 TeV. In all cases, we include photon induced channels (they can be sizeable at large energies [55]), and we set the renormalisation and factorisation scales to µ 2 R,F = 1 2 (2m 2 W + i p 2 T,i ), where the sum runs over the final-state W bosons, following the choice for gauge-boson pair production in [60]. Lastly, we set the widths of the Z and W boson consistently to 0, as the W is kept on-shell in the matrix elements. In Fig. 2 we report results for the transverse momentum of each W boson in the final state for both centre of mass energies. Plots are divided into four panels, of which the first one collects results for the various predictions. The second and third panels report the ratio to the leading order and to the EW virt  In addition, the virtual approximation EW virt and a variant of the NLO EW calculation with additional jets vetoed are also included. The ratio plots show the ratios to the LO and to the EW virt calculation, and the relative size of each NLL contribution. approximation, respectively. The aim is to show the general behaviour of EW corrections in the tail of distributions in the former, while the latter serves as a direct comparison between the Sudakov and the EW virt approximations. Finally the fourth panel shows the relative impact of the individual contributions ∆ c appearing in Eq. (2.14). Looking at the second panel, we see that all but the full NLO EW calculation show a strong suppression in the p T > m W region, reaching between −70 % and −90 % at p T = 2 TeV. This effect is the main effect we discuss in this work, and it is referred to as Sudakov suppression. To explicitly confirm that this behaviour originates from virtual contributions of EW nature, we compare the Sudakov LO+NLL curve to the NLO EW and EW virt approximations. Indeed, the latter only takes into account virtual corrections and the minimal amount of integrated counter terms to render the cross section finite. The Sudakov approximation is close to that for both centre of mass energies (see also the third panel), with deviations of the order of a few percent. It begins to deviate more at the end of the spectrum, with similar behaviours observed for both collider setups. However, with a Sudakov suppression of 70 % and more, we are already in a regime where the relative corrections become O(1) and a fixed logarithmic order description becomes invalid. In fact, at p T,W 3 TeV, the LO+NLL becomes negative both at √ s = 13 TeV and 100 TeV, the same being true for the EW virt approximation. Note in particular, that this is the main reason for choosing to show the p T distribution only up to p T = 2 TeV in Fig. 2 for both collider setups. It is also clear from the second panel that for both setups, the full NLO EW calculation does not show such large suppressions, and in the context of the question whether to use EW Sudakovs or the EW virt to approximate the full corrections, this may be worrisome. However, in the full calculation we have included the real emission matrix elements, which also show a logarithmic enhancement at high p T , as e.g. discussed in [56]. In this case, the real emission contribution almost entirely cancels the Sudakov suppression at √ s = 13 TeV, and even overcompensates it by about 20 % in the p T,W 1 TeV region for √ s = 100 TeV. To show that this is the case, we have also reported a jet-vetoed NLO EW simulation, which indeed again shows the high p T suppression as expected. Moreover, we have included a prediction labelled "LO+NLL (resum)", where we exponentiate the Sudakov contribution using Eq. (2.16). It is similar to the NLL approximation, but resumming these logarithms leads to a smaller suppression in the large p T tail, which is reduced by about 20 % at p T = 2 TeV, thus increasing the range of validity compared to the non-exponentiated prediction. Moreover, this agrees well qualitatively with the NNLL result reported in [54], and suggests that even higher-order logarithmic effects should be rather small in comparison in the considered observable range. Lastly, we compare the individual LL and NLL contributions. As expected, we find that the double logarithmic LSC term is the largest contribution and drives the Sudakov suppression. Some single logarithmic terms, in particular the C terms, also give a sizeable contribution, reducing the net suppression. This confirms that the inclusion of single logarithmic terms is needed in order to provide accurate predictions in the Sudakov approximation, with deviations of the order of O(10 %) with respect to the EW virt calculation. Comparing the individual contributions for the two collider energies, we see qualitatively similar effects. As can be expected from a larger admixture of bb initial states, the ∆ Yuk is enhanced at the larger collider energy.

EW-induced dijet production in pp collisions at 13 TeV
For the second comparison in this section, we simulate purely EW dijet production in hadronic collisions at √ s = 13 TeV at the Born-level perturbative order O(α 2 α 0 s ). As for the case of diboson production we add photon initiated channels, and we also include photons in the set of final-state partons, such that γγ and γj production is also part of our sample. Partons are thus clustered into jets which are then sorted by their p T . We select events requiring at least two jets, the leading jet (in p T ) to have a p T > 60 GeV and the subleading jet to have a p T > 30 GeV. Note in particular that for all but the real emission case in the full NLO EW simulation, this corresponds to imposing a p T cut on the two generated partons. The renormalisation and factorisation scales are set to µ 2 R,F =Ĥ 2 T = ( i p T,i ) 2 , where the sum runs over final-state partons. We compare our LO+NLL EW results, as in the previous subsection, with the LO, the NLO EW and the EW virt predictions. NLO EW corrections for dijet production have been first discussed in [61,62], while [63] discusses those corrections in the context of the 3-to-2-jet ratio R 32 . In this context, we only consider EW corrections to the Born process described above, i.e. we consider the O(α 0 s α 3 ) contributions, which is a subset of the contributions considered in the above references. In Fig. 3 (left), we present the transverse momentum distribution of the leading jet p T given by the LO+NLL EW calculation, again both at fixed order and resummed, see Eqs. (2.15) and (2.16). As for W pair production, the plot is divided into four panels. However, this time we do not include jet-vetoed NLO results, as in this case we do not observe large real emission contributions. The panels below the main plot give the ratio to the LO calculation, the ratio to the EW virt approximation, and the ratio of each NLL contribution to the LO calculation.
Compared to W pair production, we observe a smaller but still sizeable Sudakov suppression, reaching approximately −40 % at p T = 2 TeV. The NLO contributions not included in the NLO EW virt (i.e. mainly real emission terms) are small and cancel the Sudakov suppression only by a few percent and as such both the Sudakov and the EW virt approximations agree well throughout the p T spectrum with the full NLO calculation. The same is true for the resummed NLL case which gives a Sudakov suppression of −30 % for p T = 2 TeV. Note that a small step can be seen for p T ∼ m W in the Sudakov approximation. The reason for this is that we force, during the simulation, all Sudakov contributions to be zero when at least one of the invariants formed by the scalar product of the external momenta is below the W mass, as Sudakov corrections are technically only valid in the high-energy limit. This threshold behaviour can be disabled, giving a smoother transition between the LO and the LO+NLL, as can be seen for example in Fig. 2. Comparing the individual NLL contributions we find again that the double logarithmic LSC term is the largest contribution, its size is however reduced by a third compared to the diboson case, since the prefactor C ew in Eq. (2.5) is smaller for quarks and photons compared to W bosons. Among the single logarithmic terms the C and PR terms are the most sizeable over the whole range, and are of a similar size compared to the diboson case, while the SSC contribution only give a small contribution. Again, subleading terms must be included to approximate the NLO calculation at the observed 10 % level, although we observe an almost entire accidental cancellation, in this case, of the SL terms.

Off-shell Z production in association with 4 jets in pp collisions at 13 TeV
As a final example, we present for the first time the LO+NLL calculation for e + e − production in association with four additional jets from proton-proton collisions at √ s = 13 TeV. The process is one of the key benchmark processes at the LHC to make precision tests of perturbative QCD and is a prominent background constituent for several Standard Model and New Physics processes. In this case, we neglect photon induced contributions, to better compare to the NLL effects in Z plus one-jet production presented in App. A, which in turn is set up as a direct comparison to [8]. For the same reason, in this case, we only consider QCD partons in the final state. The factorisation and renormalisation scale are set to µ 2 F = µ 2 R =Ĥ 2 T = (M T, + i p T,i ) 2 /4.0, where M T, is the transverse mass of the electron-positron pair and the sum runs over the final-state partons. This choice is inspired by [64], where the full next-to-leading QCD calculation is presented. As we consider here an off-shell Z, the invariant mass formed by its decay products is on average only slightly above the W mass threshold. This may cause an issue as this is one of the invariants considered in the definition of the high-energy limit discussed in Sec. 2, and therefore this limit is strictly speaking not fulfilled. In turn, this can introduce sizeable logarithms in particular in Eq. (2.6). However, in practice we only see a small number of large K factors, that only contribute a negligible fraction of the overall cross section. For set-ups at a larger centre of mass energy, one should monitor this behaviour closely, as the average value of s in Eq. (2.6) would increase too. We therefore foresee that a more careful treatment might be required then, e.g. by vetoing EW Sudakov K factors whenever |r kl | s. In Fig. 3 (right), we present the LO+NLL EW calculation of the transverse momentum distribution of the reconstructed Z boson. To our knowledge, there is no existing NLO EW calculation to compare it against § , hence we do not have a ratio-to-NLO plot in this case. However, we do show the ratio-to-NLL plot and the plot that shows the different contributions of the LO+NLL calculation, as for the other processes. The sizeable error bands give the MC errors of the LO calculation. Note that the errors of the LO and the LO+NLL calculation are fully correlated, since the NLL terms are completely determined by the phase-space point of the LO event, and the same LO event samples are used for both predictions. Hence the reported ratios are in fact very precise. With the aim of additionally making a comparison to the Z + j result reported in Fig. 4 (right panel), we opt for a linear x axis, in contrast with the other results of this section. In both cases we see a similar overall LO+NLL effect, reaching approximately −40 % for p T 2 TeV. This in turn implies that the effects of considering an off-shell Z (as opposed to an on-shell decay) as well as the additional number of jets have very little effect on the size of the Sudakov corrections. Finding similar EW corrections for processes that only differ in their number of additional QCD emissions can be explained by the fact that § Note that while such a calculation has not yet been published, existing tools, including SHERPA in conjunction with OpenLoops [65], are in principle able to produce such a set-up. the sum of EW charges of the external lines are equal in this case. As has recently been noted in [53], this can be deduced from the general expressions for one-loop corrections in [1] and from soft-collinear effective theory [28,30]. Although the overall effect for Zj and Z + 4j is found to be very similar here, the individual contributions partly exhibit a different behaviour between the two, with the SSC terms becoming negative in the four-jet case and thus switching sign, and the C terms becoming a few percent smaller. It is in general noticeable that the SSC terms exhibit the strongest shape differences among all processes considered in this study. Finally, similarly to the previous studied cases, the resummed result gives a slightly reduced Sudakov suppression, reaching approximately −30 % for p T 2 TeV, implying that in this case, higher logarithmic contributions should be small.  Figure 3: The transverse momentum of the leading jet in EW-induced dijet production in proton-proton collisions (including photon channels), and for the reconstructed Z boson in e + e − plus four jets production, For the dijet production, LO and NLO calculations are shown, whereas for the Z plus jets production only the LO is shown. These baseline calculations are compared with the results of the LO+NLL calculation, both at fixed-order and resummed. In the dijet case, the virtual approximation EW virt is shown in addition. The ratio plots show the ratios to the LO and the EW virt calculations, and the relative size of each NLL contribution.

Conclusions
We have presented for the first time a complete, automatic and fully general implementation of the algorithm presented in [1] to compute double and single logarithmic contributions of one-loop EW corrections, dubbed as EW Sudakov logarithms. These corrections can give rise to large shape distortions in the high-energy tail of distributions, and are therefore an important contribution in order to improve the accuracy of the prediction in this region. Sudakov logarithms can provide a good approximation of the full next-to-leading EW corrections. An exponentiation of the corrections can be used to resum the logarithmic effects and extend the region of validity of the approximation. In our implementation, each term contributing to the Sudakov approximation is returned to the user in the form of an additional weight, such that the user can study them individually, add their sum to the nominal event weight, or exponentiate them first. Our implementation will be made available with the next upcoming major release of the SHERPA Monte Carlo event generator. We have tested our implementation against an array of existing results in the same approximation, and for a variety of processes against full NLO EW corrections and the virtual approximation EW virt , which is also available in SHERPA. A selection of such tests is reported in this work, where we see that indeed EW Sudakov logarithms give rise to large contributions and model the full NLO corrections well when real emissions are small. We stress that our implementation is not limited to the examples shown here, but it automatically computes such corrections for any Standard Model process. In terms of final-state multiplicity it is only limited by the available computing resources.
In a future publication we plan to apply the new implementation in the context of state-of-the-art event generation methods, in particular to combine them with the matching and merging of higher-order QCD corrections and the QCD parton shower, while also taking into account logarithmic QED corrections using the YFS or QED parton shower implementation in SHERPA. This will allow for an automated use of the method for the generation of event samples in any LHC or future collider context, in the form of optional additive weights. As discussed towards the end of Sec. 2, we foresee that this is a straightforward exercise, since there is no double counting among the different corrections, and applying differential K factors correctly within these methods is already established within the SHERPA framework. The result should also allow for the inclusion of subleading Born contributions, as the EW virt scheme in SHERPA does. The automated combination of (exponentiated) EW Sudakov logarithms with fixed-order NLO EW corrections is another possible follow-up, given the presence of phenomenologically relevant applications.
A Validation plots for V plus jet production in pp collisions at 13 TeV As an additional validation for our implementation, we show in Fig. 4 the LO+NLL calculation for Z, W + and γ production, each in association with one jet. The setups, the analysis, and the binning are configured to closely match the choices made in [8,11,9], respectively, and the reader is referred to these publications for further details. In particular, the photon is not considered to be part of the parton content of the incoming protons, and the MRST01 LO PDF set has been used [66]. We have confirmed that our NLL predictions agree with the reference NLL ones at a sub-percent level over the whole range of the observables.