The midpoint between dipole and parton showers

We present a new parton-shower algorithm. Borrowing from the basic ideas of dipole cascades, the evolution variable is judiciously chosen as the transverse momentum in the soft limit. This leads to a very simple analytic structure of the evolution. A weighting algorithm is implemented, that allows to consistently treat potentially negative values of the splitting functions and the parton distributions. We provide two independent, publicly available implementations for the two event generators Pythia and Sherpa.


Introduction
Parton showers and fragmentation models have been used for more than three decades to predict the dynamics of multi-particle final states in collider experiments [1]. More recently, the traditional approaches implemented in HERWIG [2], PYTHIA [3] and SHERPA [4] were supplemented by methods based on dipole and antenna factorization [5,6,7,8,9]. A characteristic feature of these new shower programs is the description of QCD coherence in the color dipole picture [10], which has first been implemented in the ARIADNE Monte Carlo [11]. In this article we present a dipole-like parton shower similar to existing ones, but we focus on the simplest implementation and enforce sum rules and DGLAP collinear anomalous dimensions. We choose ordering variables based on transverse momenta in the soft approximation, while existing dipole-like shower models employ collinear transverse momenta. As such, the model is a hybrid of dipole and parton shower. These choices will eventually allow to compare with analytic approaches, such as CSS [12] and SCET [13]. In the past decade, the matching of parton showers to NLO calculations [14,15] and the merging of LO [16] and NLO matched results for different jet multiplicity [17,18] was in the focus of interest of the majority of Monte-Carlo developers [19]. Comparably few efforts were made to provide publicly available implementations of parton-showers [5,6,8,9] or to improve their formal accuracy [20], and even fewer of the new parton showers have made their way into complete event generators used by experiments. When comparing results of matched and merged calculations, it is therefore often unclear whether a particular difference stems from mismodeling in the parton shower, from differences in the matching or merging algorithm, or simply from technical problems. Similarly, when comparing the results of different event generators at the hadron level it is often unclear whether differences should be ascribed to the hadronization model, to the simulation of multiple scattering / rescattering effects, or to the parton shower. We intend to remedy this situation to some extent, by providing two implementations of one and the same algorithm, to be used with the two different event generation frameworks PYTHIA [21] and SHERPA [22]. We subject our codes to rigorous scrutiny by comparing their predictions at the sub-permille level. This paper is organized as follows: Section 2 reviews the basic parton-shower formalism. Section 3 explains the construction principles of our new parton shower, which we call DIRE (acronym for DIpole REsummation). Section 4 contains the validation of the numerical implementation, and Sec. 5 presents a comparison of the predictions from DIRE with experimental measurements. Section 6 contains some concluding remarks. arXiv:1506.05057v2 [hep-ph] 19 Sep 2015 2 Parton-shower formalism The evolution of parton densities and fragmentation functions in the collinear limit is governed by the DGLAP equations [23]: where P ab are the regularized evolution kernels. Assume that we define P ab in terms of unregularized kernels, P ab , restricted to all but an ε-environment around the soft-collinear pole, plus an endpoint contribution.
For finite ε, the endpoint subtraction can be interpreted as the approximate virtual plus unresolved real corrections, which are included in the parton shower by enforcing unitarity. The precise value of ε is defined in terms of an infrared cutoff on the evolution variable, using four-momentum conservation. When ignoring momentum conservation, this cutoff can be taken to zero, which allows us to identify [P ba (z)] + as the ε → 0 limit of P ba (z, ε). For 0 < ε 1, Eq. (2.1) changes to Using the Sudakov form factor one can define the generating functional for splittings of parton a as In this context, Π a (x, t, µ 2 ) is the probability that the parton does not undergo a branching process between the two scales µ 2 and t [24]. Equation (2.3) can now be written in the simple form The generalization to an n-parton state can involve multiple PDFs and fragmentation functions: (2.8) In this context, we have extended the argument of the generating functional toΦ n , which denotes the nparton phase-space configuration, including all light-cone momentum fractions, x i and z j , for initial-state (IS) and final-state (FS) partons. F also depends on all parton flavors, denoted by a. If we do not fix the momenta of final-state hadrons, the fragmentation functions can be integrated over z j , leading to the simplified formula The change fromΦ n to Φ n signals that F has become independent of z j . An observable-dependent generating functional for the parton shower can now be defined recursively as (2.10) The first term includes all virtual corrections and unresolved real emissions, resummed into a no-branching probability. The second term describes a single branching, followed by further parton evolution. Both terms can be generated simultaneously by implementing the veto algorithm [25] for Eq. (2.6). We have introduced an observable, O, that measures the kinematics of the final state. In general, this observable will act differently on the no-emission term and on the emission term. In the trivial case that O = 1, Eq. (2.10) returns the unitarity constraint, F a (Φ n , t, t ; 1) = 1. Generating a branching in the parton shower involves selecting a new color topology for the n + 1-particle state. For non-trivial color configurations, F will therefore depend on the color assignment in the large-N c limit. While it is in principle necessary to keep track of this dependence, we omit any notation relating to color in order to simplify our final formulae. The selection of color topologies proceeds as in existing dipolelike parton showers, which is described in great detail in [6]. It is straightforward to extend our notation in this regard. The choice of evolution variable is crucial. At leading color the soft radiation pattern emerges from the coherent gluon radiation off "color dipoles" that are spanned by the two partons at opposite ends of a color string [26]. This mandates the choice of an evolution variable which treats these two partons democratically. In other words, the evolution variable should be identical no matter whether one or the other of two colorconnected partons is considered the radiator. This will be the guiding principle for its selection in Sec. 3. The splitting functions, P (z), are formally defined in the collinear limit, and they do not reflect the soft radiation pattern outside the collinear region. Traditionally, this problem is dealt with by imposing angular ordering constraints on the final-state phase space [2]. Alternatively one can use the approach of Catani and Seymour [27], and introduce a t-dependence in the splitting functions that restores the correct soft anomalous dimension at one-loop order [6,9]. We will use this procedure in the next section. It is important that the modified splitting functions satisfy the sum rules, which are enforced by Eq. (2.2) and by the corresponding flavor sum rule [24]. The new splitting functions may also be negative in the non-singular phase-space region. This requires a modification of the Sudakov veto algorithm [28,29,15], and it entails an analytic event weight to ensure that both emission-and no-emission probabilities are accounted for. We find that, in our parton-shower approach, the variance of this weight is small. In fact, for both final-final and initial-initial dipoles momentum conservation guarantees that no negative weights can arise from the splitting functions. Negative weights may however appear in initial-initial configurations due to negative values of the PDFs.

Construction of the DIRE shower
A basic branching process is sketched in Fig. 1. In this case we consider initial-state evolution. We employ the kinematics from Ref. [27,30], which we summarize in Appendix A. For initial-state splitters with initialstate spectator, all particles typically have zero on-shell mass, which greatly simplifies the calculation. Their momenta can be parametrized in terms of the light-cone momenta p a and p b , using the standard Sudakov decomposition [31] p µ aj = z p µ a + In the Catani-Seymour approach [27], the correct soft anomalous dimension is obtained after replacing the soft enhanced term in the splitting functions by a partial fraction of the soft eikonal for the color dipole defined by the splitting parton and its spectator. Schematically this can be done as follows: If we define the evolution variable of our parton shower to be a scaled transverse momentum, t = (z − v) k 2 ⊥ , then the soft-enhanced term in Eq. (3.2) is conveniently expressed as Note that the evolution variable has the desired symmetry property, i.e. it is symmetric in emitter and spectator momentum. More precisely, our evolution variable is the exact inverse of the soft eikonal. As such, it is different from the hardness parameter, k 2 ⊥ . Consequently, the parton shower will fill the entire final-state phase space, even for factorization scales much smaller than the hadronic center of mass energy. We define the evolution kernels using the replacement of the soft enhanced term in Eq. (3.2). Additionally, we require the collinear anomalous dimension to be unchanged. Imposing flavor and momentum sum rules, these two requirements determine the complete set of leading-order spin-averaged splitting functions as: The corresponding anomalous dimensions are listed in Appendix B. Using the phase-space factorization derived in [27], we obtain the following differential branching probability: the light-cone momentum fraction entering the PDFs is well defined. This variable has changed compared to Eq. (2.1) to account for four-momentum conservation, while Eq. (2.1) remains valid in the collinear limit, v → 0. In addition, the strong coupling is evaluated at the evolution variable, hence the Landau pole is avoided by the infrared cutoff of the parton shower, t 0 (which is of order 1 GeV). Finally, the splitting kinematics are constructed as described in App. A.5.
The technical implementation of Eqs. (3.4) in terms of dipole terms proceeds as described in [6]. That is, we divide the splitting function according to the number of spectator partons in the large-N c limit, and sum over color-adjacent splitter-spectator pairs. The corresponding evolution equations are straightforward extensions of Eq. (3.5), and therefore we do not list them here. The same reasoning applies to all dipole types discussed in the following. Initial-state splitters with final-state spectator are treated along the same lines. The construction of finalstate momenta is described in App. A.3. The kinematics are now defined in terms of the following variables (p k is the final-state spectator momentum) We define the evolution variable, t, and its dimensionless variant, κ 2 , as The evolution variable is symmetric in emitter and spectator momentum, up to a rescaling by 1/z 2 . The differential branching probability is: Final-state splitters with initial-state spectator must have the same evolution variable due to symmetry arguments. Therefore the asymmetric rescaling by 1/z 2 in the IF case must also be applied in the FI case.
The only choice to be made concerns the splitting variable, which is taken to be 1 − u. The differential branching probability is: The construction of final-state momenta is described in App. A.2. For final-state splittings, ij → i, j, with final-state spectator, k, we use the variables The symmetric evolution variable, its scaled variant κ 2 and the splitting variable are defined as The differential branching probability is: where 2 z ± = 1 ± 1 − 4 t 0 /Q 2 . The splitting kinematics are described in App. A.1. Note that the scaled transverse momentum defined in Eq. (3.11) is substantially different from the ones defined in [6,9], which can be written as This variable is symmetric in i and j, but not in i and k, which would be required in order to interpret it as the inverse soft eikonal for gluon radiation off the dipole spanned by i and k. Kinematically, Eq. (3.13) represents the transverse momentum of partons i / j with respect to the anti-collinear direction defined by k. This is what we call a "collinear" transverse momentum. In contrast, Eq. (3.11) can be interpreted as the transverse momentum of the two daughter dipoles (ij) / (kj) in the center-of-mass frame of the decaying dipole [11]. In this case, i defines the collinear, and k defines the anti-collinear direction, making the symmetry explicit. We refer to such a definition as a "soft" transverse momentum. The change in the definition of transverse momenta compared to existing p T -ordered dipole-like parton showers [6,9] also involves changing the splitting variable, in order to reduce the related Jacobians to unity while maintaining Eqs. (3.4), simultaneously for all dipole types. In contrast, the kinematics mapping is identical to the previously published methods [6,32]. If massive quarks are involved in the branching process, we would like to map the evolution variable to the soft enhanced term of the full matrix element, just like in the massless limit. The singularity structure in the soft limit is given in [30]. For the most involved case of two massive radiators, i and k, it leads to an eikonal of the form (3.14) The inverse of this function is difficult to interpret. Its scaling property in the soft limit, however, is completely determined by the first term in Eq. (3.14), whose inverse can therefore be used to define an ordering variable for the evolution of massive partons with γ defined in Appendix A, and with s ik = (p i + p k ) 2 . Here we have defined the massless equivalent of the evolution variable, k (0) 2 ⊥ , and the light-cone momentum fractions, ζ i and ζ j in a Sudakov decomposition of the gluon momentum, p j , along the directions of p i and p k : Equation (3.15) is valid in the soft limit. For practical purposes the denominator p i p k in the evolution variable should be the hard scale of the radiating dipole, which is given by (2 p a p b ) 2 /Q 2 , (2 p a (p i + p k )) 2 /Q 2 and Q 2 for II, IF/FI and FF dipoles, respectively. The splitting functions for massive partons can be taken from Eq. (3.4) and be modified according to [30]. We use the following unregularized massive kernels for final-state splitter with final-or initial-state spectator: (3.17) The relative velocity between two momenta, p and q, is defined as and v ı,k and v ij,k stand for the relative velocities between the emitter parton and the spectator before and after the branching, respectively. The branching probabilities are modified asP ab (z, κ 2 ) → J(z, κ 2 )P ab (z, κ 2 ), where J(z, t) is a spectator-dependent Jacobian factor [30,6]. It is unity, except for the case of final-state splitter with final-state spectator, where using Q 2 = 2 p i p k + 2 (p i + p k )p j , cf. Eq. (3.10). The phase-space boundaries are given by the roots of the Gram determinant While the massless case leads to simple constraints on z, the general massive case generates a rather involved functional form of the z-boundary as a function of t. Algorithmically, it is preferable to use the veto algorithm [25] to implement this constraint, or to use Eqs. (3.11) and (3.7) and evaluate the constraint in collinear variables, where [30] In final-state splittings with initial-state spectator the PDF is evaluated at For initial-state splitter with final-state spectator the mass-dependent splitting functions arê

Validation
In this section we validate the numerical implementation of the DIRE parton shower. The two event generation frameworks PYTHIA [3] and SHERPA [4] are used to construct two entirely independent Monte Carlo programs. Aside from a thorough cross-check of the implementation, this allows, for the first time, to simulate Deep Inelastic Scattering in PYTHIA 8. We employ the CT10nlo PDF set [33], and the corresponding value of the strong coupling. Following standard practice to improve the logarithmic accuracy of the parton shower, the soft enhanced term of the splitting functions is rescaled by 1 + α s (t)/(2π)K, where K = (67/18 − π 2 /6) C A − 10/9 T R n f [34].  Fig. 2 left), e + p → e + jet (Fig. 3 left), and pp → e + e − (Fig. 4 left). These three cases cover all possible dipole types, i.e. final-state splitter with final-state spectator, final-state splitter with initial-state spectator and vice versa, and initial-state splitter with initial-state spectator. Gluon splitting functions are validated in the reactions τ + τ − → hadrons (Fig. 2 right), τ + p → τ + jet (Fig. 3 right), and pp → τ + τ − (Fig. 4 left), all mediated by Higgs-boson exchange using HEFT [35].

Results
In this section we compare DIRE predictions from SHERPA [22] with experimental data. When applicable, we use the CT10nlo PDF set [33] and the corresponding strong coupling. We employ the kinematics scheme from App. A.3. Our results include the simulation of QED radiation in the case of Drell-Yan lepton pair production [36], and hadronization in the case of e + e − → hadrons [37]. Otherwise they are given at the parton level in order to exhibit the features of the DIRE shower only. Analyses are performed with RIVET [38]. Figure 5 shows predictions from the DIRE parton shower for differential jet rates in the Durham scheme compared to experimental results from the JADE and OPAL collaborations [39]. The perturbative region is to the right of the plot, and y ∼ 2.8·10 −3 corresponds to the b-quark mass. The simulation of nonperturbative effects dominates the predictions below ∼ 10 −4 . We observe that, in the perturbative region, the results are in excellent agreement with the experimental measurements. Figure 6 shows a comparison for event shapes measured by the ALEPH collaboration [40]. The perturbative region is to the right of the plot, except for the thrust distribution, where it is to the left. We notice some deviation in the predictions for jet broadening and for the C-parameter. However, these deviations are mostly within the 2σ uncertainty of the experimental measurements, and they occur close to the nonperturbative region. It can also be expected that the simulations improve upon including matrix-element corrections or when merging the DIRE shower with higher-multiplicity calculations. This has been demonstrated, for example, in [17]. Figure 7 shows angular correlations in comparison to ATLAS data from [41], and the transverse momentum spectrum of the Drell-Yan lepton pair in comparison to ATLAS data from [42]. It is well known that pure parton shower predictions are insufficient to describe these measurements. Therefore, we merge our parton shower with 1-jet matrix elements using the CKKW-L procedure [16]. In order to assess the related uncertainties, we vary the merging cut by a factor 2 around the central value of Q cut = 10 GeV. The associated uncertainty band is shown in light red. The size of the variation is comparable to the statistical uncertainties, which are displayed as error bars on the Monte-Carlo prediction. Figure 8 shows di-jet azimuthal decorrelations in different regions of jet transverse momentum. We compare DIRE predictions with experimental results from the ATLAS collaboration [43]. This observable tests for higher-order effects in some detail [44].

Conclusions
We presented a new dipole-like parton-shower algorithm, that is constructed along very simple arguments: Firstly, the ordering variable should exhibit a symmetry in emitter and spectator momenta, such that the dipole-like picture can be re-interpreted as a dipole picture in the soft limit. At the same time, the splitting functions are regularized in the soft anti-collinear region using partial fractioning of the soft eikonal in the Catani-Seymour approach. They are then modified to satisfy the ordinary sum rules in the collinear limit. This leads to an invariant formulation of the parton-shower algorithm, which is in complete analogy to the standard DGLAP case. We computed the anomalous dimensions, which match previous results for dipolelike parton showers. We presented first phenomenologically relevant predictions using the new algorithm, and we observe very good agreement with experimental data from LEP and LHC experiments.

A Parton-shower kinematics
The precise algorithm for constructing the splitting kinematics depends on the type of splitter and spectator parton. There are four separate cases. Note that initial-state partons are assumed to be massless for collinear PDF evolution to be valid. In practically implemented parton-shower algorithms they are often taken massive instead, in order to obtain a better description of experimental data. Therefore we give the kinematics formulae with full mass-dependence, including initial-state parton masses. 2 A.1 Final-state splitter with final-state spectator 1. Determine the new momentum of the spectator parton as 2. Construct the new momentum of the emitter parton, p i , as The parametersz i and k 2 ⊥ = −k 2 ⊥ of this decomposition are given bȳ In the massless case, this algorithm reduces to the following [6] p wherez i = (z − y ij,k )/(1 − y ij,k ) and y ij,k = κ 2 /(1 − z), cf. Sec. 3.
A.2 Final-state splitter with initial-state spectator 1. Determine the new momentum of the spectator parton as 2. Proceed as in Sec. A.1, except that m k → m a and p k → −p a .
In the massless case, this algorithm reduces to the following [6] p µ i =z ip where q =p aj −p k , q = q +p k⊥ and s jk = (1 − 1/x jk,a ) (q 2 − m 2 a ) + (m 2 j + m 2 k )/x jk,a . 2. Proceed as in Sec. A.1, except thatz j → u j , m k → m a , m j → m k , p k → −p a and p j → p k .

A.4 Initial-state splitter with final-state spectator (global recoil)
This scheme can be chosen as an alternative to the one in Sec. A.3. The algorithm is equivalent to the method outlined in [28,32]. It alleviates a formal problem with transverse momentum resummation in Drell-Yan type processes [45], but it is numerically less stable, and therefore employed only if |x jka − u j | > ε, where ε ∼ 10 −4 .
1. Determine the new momentum of the spectator parton as where q =p k −p aj and s aj = u j /x jk,a (q 2 − m 2 k ) + (m 2 a + m 2 j ) (1 − u j )/x jk,a . 2. Construct the momentum of the emitted particle, p j , as The parametersz j and k 2 ⊥ = −k 2 ⊥ of this decomposition are given bȳ 3. Boost p a and all final state particles into the frame where p a is aligned along the beam direction, with p b , the opposite-side beam particle, unchanged.
In the massless case, this algorithm reduces to the following [9] 3 (A.12) Note that, following the arguments in [9], the light-cone momentum fraction entering the PDF is still given by A.5 Initial-state splitter with initial-state spectator 1. Determine the new momentum of the splitting parton as where q =p a + p b and s ab = (q 2 − m 2 j )/x j,ab + (1 − 1/x j,ab ) (m 2 a + m 2 b ). 2. Construct the momentum of the emitted parton, p j , as The parametersz aj and k 2 ⊥ = −k 2 ⊥ of this decomposition are given bȳ 3. Boost all remaining final-state particles into the frame defined by p a + p b − p j , using the algorithm defined in Sec. (5.5) of [27].
In the massless case, this algorithm reduces to the following [6] p µ a = 1 x j,abp µ aj , where x j,ab = z −ṽ j andṽ j = κ 2 /(1 − z), cf. Sec. 3, and where k i stands for any final-state momentum, including EW particles. The Lorentz transformation, Λ, is given by The inverted kinematics needed for matching the DIRE shower to NLO computations and for merging with higher-order tree-level calculations are given in [30]. The combined momenta are obtained from Eqs. (A.1), (A.5), (A.7) and (A.13) by swapping momenta and masses before and after emission. In the case of final-state emitter with final-state spectator, for example, this amounts to the replacement p k ↔p k and s ij ↔ m 2 ij .

C Momentum mapping and q T spectra in Drell-Yan type processes
This section presents a brief phenomenological analysis of the different recoil strategies described in Apps. A.3 and A.4. We investigate the impact on the transverse momentum (q T ) spectrum of the Drell-Yan lepton pair at the LHC. It was pointed out [45] that the non-singlet initial-state parton evolution is generated in dipole-like parton showers by initial-state splitters with final-state spectator, except for the first branching, which stems from an initial-state splitter with initial-state spectator. The kinematics mapping described in App. A.3 then results in the Drell-Yan lepton pair acquiring its entire transverse momentum in the first branching 4 . Comparing to analytical resummation in impact parameter space [12], one is led to conclude that only leading logarithms can be resummed in this scheme [46]. However, comparing to resummation in transverse momentum space [47], the modified leading logarithmic structure characteristic for parton showers emerges from Eq. (2.10). A kinematics mapping more appropriate for comparison with analytical approaches is given by App. A.4. Here we restrict ourselves to analyzing its impact on the q T spectrum at the LHC. Similar analyses have been performed elsewhere for a variety of other observables using standard dipole-like parton showers based on Catani-Seymour subtraction [28,32]. Figure 9 shows DIRE predictions from SHERPA for the two different kinematics schemes described in Apps. A.3 and A.4. To highlight the differences in the resummation, we present pure parton shower results. Correspondingly, experimental data are omitted. The global recoil scheme from A.4 shows a small difference at low q T . Similar conclusions were reached in [28,32].