Improvements on dipole shower colour

The dipole formalism provides a powerful framework from which parton showers can be constructed. In a recent paper, we proposed a dipole shower with improved colour accuracy and in this paper we show how it can be further improved. After an explicit check at $\mathcal{O}(\alpha_{\mathrm{s}}^{2})$ we confirm that our original shower performs as it was designed to, i.e. inheriting its handling of angular-ordered radiation from a coherent branching algorithm. We also show how other dipole shower algorithms fail to achieve this. Nevertheless, there is an $\mathcal{O}(\alpha_{\mathrm{s}}^{2})$ topology where it differs at sub-leading $N_{\mathrm{c}}$ from a coherent branching algorithm. This erroneous topology can contribute a leading logarithm to some observables and corresponds to emissions that are ordered in $k_t$ but not angle. We propose a simple, computationally efficient way to correct this and assign colour factors in accordance with the coherence properties of QCD to all orders in $\alpha_{\mathrm{s}}$.


Introduction
Parton showers typically are constructed using one of two basic approaches: angularordered showers (based on the coherent branching formalism) and dipole showers. Angular ordering is a very powerful approach, providing next-to-leading logarithmic (NLL) accuracy in some observables, 1 but it fails to capture physics salient to the description of multi-jet final states in hadron colliders and non-global observables. By comparison, dipole showers are typically restricted to leading-colour accuracy but they can be applied across the board. In recent literature, much attention has been focused on improving the framework upon which dipole showers are constructed [4][5][6][7][8][9][10][11][12][13][14]. Substantial progress has been made demonstrating their capacity for NLL resummation [1,15] and methods for partially addressing sub-leading colour have also been proposed, by extending dipole showers beyond leading-N c colour flows [16][17][18][19]. In a recent paper [1], we constructed a dipole shower that has the virtue of inheriting some of the colour dynamics of an angular-ordered shower, which improved sub-leading colour accuracy. In this paper we perform a fixed-order cross-check of that approach. We do so by comparing the improved shower's assignment of colour factors to the corresponding exact e + e − matrix elements, computed with second-order QCD corrections. Motivated by these calculations, we are able to further improve our dipole shower's description of colour, in a way that is applicable to evolution with an arbitrary number of emissions.
In [1] we derived an improved dipole shower in the context of e + e − → qq collisions 2 , starting from an algorithm for the evolution of QCD amplitudes first presented in [20]. The shower can be understood by considering a few key features of angular-ordered and dipole showers. When a shower emits a parton, three new degrees of freedom (DoF) are introduced, describing the new parton's energy and direction. Angular-ordered showers average over one of the DoF (a contextually defined azimuth) which allows the effects of QCD coherence to become manifest. In turn, this reduces the shower to a Markovian sequence of parton decays (1 → 2 transitions). Thus the final-state partons produced by the shower have a unique branching history with colour factors assigned in accordance with QCD coherence. The angular-ordered approach is very powerful; by harnessing QCD coherence, NLL resummation can be achieved for a broad class of global observables [2]. However, averaging a DoF limits the approach.
In contrast to angular ordering, the dipole approach retains full dependence on the DoF of each emitted parton. Instead it approximates the colour structures in the shower by emitting partons from colour-connected dipoles. This restricts a basic dipole shower to leading-colour accuracy. Thus a dipole shower is built from a Markovian sequence of 2 → 3 transitions and, as a result, dipole showers lack a unique branching history of parent partons and their decay products. However, a branching history can be constructed by introducing a dipole partitioning, which probabilistically assigns the emitted parton to one of the two parent partons in the dipole. Modern dipole showers use this partitioning to assign colour factors and momentum conservation, and to facilitate hard-process matching. In effect, our approach in [1] was to define a partitioning so that, after averaging over azimuths, each branching history and its relative weight matches with a corresponding branching history generated by an angular-ordered shower. Through this link, we could assign colour factors beyond leading colour in the dipole shower. As we show in this paper, when applied naively (as was done in [1]) this procedure does not completely eliminate sub-leading colour errors in the dipole shower for some observables, even at LL accuracy. The problem arises since the k t -ordered dipole shower necessarily involves branching histories disordered in angle: soft, large-angle emissions can appear anywhere in the branching history. These particular branching histories complicate any attempt to assign colour factors in a dipole shower (a point previously noted in [21]) and were not completely accounted for in our original approach. In this paper, we solve this problem by introducing dynamical colour factors, i.e. we fix the LL, sub-leading colour errors in event shape observables and increase the shower's sensitivity to full-colour NLLs (falling short of full-colour NLL resummation).
The rest of the paper is structured as follows. After a review of the double emission matrix element in Section 2, we repeat the calculation for our original dipole shower in Section 3 and compare the two. We find that the shower works as intended, i.e. the colour factors assigned to partons whose emissions are ordered in angle agree with those of the fixed-order result. However, for emissions unordered in angle, the shower has only leading-colour accuracy. The understanding brought about by the fixed-order analysis allows us to construct a new method for the correct assignment of dynamical colour factors for emissions unordered in angle. The specific partitioning we introduced in [1] plays a crucial role in the construction of the new colour factors. The approach we take involves altering shower kernels by introducing a dynamic colour factor that is a function of the branching history. This method involves a computational complexity that asymptotically grows logarithmically with the parton multiplicity. Finally, to illustrate the importance of the dipole partitioning, we compute the O(α 2 s ) difference between exact squared matrix elements and those calculated using a dipole shower employing a different (Catani-Seymour [22]) partitioning. We find that, in the limit the emissions are strongly ordered in energy and angle, the O(α 2 s ) difference does not vanish, with the possibility of a LL, sub-leading colour error, as was noted in [21]. For specific observables (e.g. thrust) this error may be removed at order O(α 2 s ) from a dipole shower employing a Catani-Seymour-type dipole partitioning by using our dynamic colour factors. However, the error will likely re-emerge at higher-orders.

A recap of the O(α s ) QCD squared matrix element
First we recap the calculation of the O(α 2 s ) e + e − → qqgg squared matrix element when the gluons are either soft or collinear. Figure 1 illustrates our labelling of the partons and the angles between them. This calculation is essentially a recap of Section 5.5 in Ellis, Stirling and Webber [23] and Chapter 4 in Dokshitzer, Khoze, Mueller and Troyan [24]. To start, we will only take the limit that lab-frame energies satisfy E q E j , E i , E k (i.e. the pure q θ jq i k j θ ik θ ij Figure 1. One of the Feynman diagrams contributing to the O(α 2 s ) e + e − → qqgg matrix element used to compute Eq. (2.2). In the present work we calculate these in the soft approximation for which the second gluon is assumed to have energy much less than the first emission. soft limit for q). Thus our starting point is 3 and where |M 1 ( p i , p j , p k )| 2 is the O(α s ) squared matrix element. At leading colour (LC) we have This can be interpreted as a sum of emissions from two independent dipoles, (ij) and (jk), and is the basic result on which dipole showers and the Banfi-Marchesini-Smye (BMS) equation [3] are built, see also the discussion in [25] for a more detailed analysis in the case of more general processes. Without approximating colour, we can simplify the matrix element by only keeping terms which are logarithmically enhanced in the two-jet limit (i.e. terms that diverge as θ ij /θ ik → 0). To do this, we write each w ab as q is the azimuth as measured around the direction of p a . We define the following shorthand for averaging the azimuths: ab only depends on parton b via the theta function constraining the angle of emission. We now consider the limit θ ij θ ik = π, whence we can assume where (ij) refers to the momentum p (ij) = p i + p j , which is approximately on-shell in the collinear limit we consider. By employing this and similar relations we can simplify the matrix element: The four contributions are illustrated in Figure 2 and Table 1 tabulates each term in Eq. (2.9) across the entire emission phase-space for q. Table 1. The contributions to Eq. (2.9), the azimuthally-averaged squared matrix element in the limit θ ij θ ik . Terms where wide-angle, soft physics has been lost as a result of the collinear approximation are underlined.
It is important to note that when deriving Eq. (2.9), terms such as were set to zero by approximating the direction of i and j with a combined direction (ij). For finite θ ij , these terms are only subject to energy divergences and therefore are negligible so long as we insist that a collinear logarithm is picked up. However, because we azimuthally averaged and neglected these pieces, some wide-angle physics is lost which is otherwise captured at LC in Eq. (2.3) (and consequently in dipole showers and BMS evolution). These soft poles are crucial to a complete description of non-global logarithms. Regions of phase-space for which some wide-angle physics has been set to zero are underlined in Table 1. Crucially, all wide-angle soft physics is included in the limit θ ij → 0. Eq. (2.9) can be generalised to include the situation where the parton energies are no longer strongly ordered by introducing hard-collinear physics: For instance, in the limits defining the rows 1 through 3 of Table 1, Eq. (2.9) with hardcollinear physics is ji dz j P gg + P where P ab is an unregularised splitting function and 1 − z m ≈ E q /E m with E m the energy of parton m before the emission 4 . Likewise, Eq. (2.9) in the limit defined in the last row of Table 1 becomes 3 Computing the squared matrix element with the dipole shower Now we want to compute the same squared matrix element using our k t -ordered dipole shower. We wish to test that it correctly reproduces the terms in Table 1 and the hardcollinear physics in Eqs. (2.12) and (2.13). The relevant contributions are pictured in Figure 3. Consider a generic k t -ordered dipole shower, for which emission from a dipole (a, b) at a given slice in k t is generated by an emission probability of dProb ln k ab where P d q→qg (z) = P d qq (z) and P d g→gg (z) = P d gg (z). We neglect g → qq transitions, which are sub-leading in colour and only contribute a NLL for doubly-logarithmic observables 5 . P d qq (z) and P d gg (z) are the usual dipole splitting functions, stripped of their colour factors: and where C a = C F or C A /2 if parton a is a quark or gluon. The g ab are dipole partitioning functions, they define how colour factors and momentum conservation should be distributed across the two members of a dipole and are functions of the momenta of all partons emitted so far. Functions g ab can be smooth or discontinuous functions of the parton momenta.
Since we are neglecting momentum conservation in this section, g ab + g ba = 1. The relevant kinematic variables are and φ is an azimuth so that k ab ⊥ = k ab ⊥ (sin φ n 1 + cos φ n 2 ) where n 1,2 are two mutually orthogonal and normalised transverse vectors in the (a, b) dipole zero-momentum frame.
After two emissions, the shower gives

O(α 2 s ) with emissions ordered in angle
We will first consider whether our dipole shower can recreate the physics in rows 1 through 3 of Table 1. The diagrams contributing to this limit will be produced in our shower when the parton transverse momenta and angles are concurrently ordered. For now we will neglect recoil and hard-collinear pieces. Keeping only the soft parts, we have Our dipole shower is built using the partitioning where the sum over i is a sum over all partons in the event. T plays the role of projecting the lab frame energy when it is contracted with a momentum vector. Roughly speaking, this way of partitioning a dipole corresponds to splitting the dipole in half in the laboratory frame; it is defined specifically to ensure and therefore After azimuthal averaging: As we are working in the limits defined in rows 1 through 3 of Table 1, and in the soft limit for q, the k t ordering theta functions all saturate and can be removed. Thus, we find This is the same as the fixed-order result given in the previous section.
We will now relax the soft approximation and check whether our dipole shower correctly includes the hard-collinear physics too, i.e. that it reconstructs Eq. (2.12). We can once again start from Eq. (3.3): where the 'hard pieces' part depends on the splitting function: Since the collinear limit requires θ iq or θ jq θ ij θ ik , we can let z ab ≈ z a where 1 − z a ≈ E q /E a and where E a is the energy of parton a before the emission is generated. Thus 'hard pieces' does not have any azimuthal dependence 6 and so azimuthal averaging proceeds as before. We find that Once again, the collinear limit results in P jk . The small-angle approximation for q saturates the ordering theta functions and so This is equivalent to the fixed-order result of Eq. (2.12). An important part of this section was to assume we can neglect recoil and that further emissions do not modify momenta in such a way that these correctly computed matrix elements are destroyed. As we showed explicitly in Section 3.1 of [1], our global recoil does not mess the computation of NLLs at this order. This is further discussed in Section 4.1.

O(α 2 s ) with emissions unordered in angle
In the previous section we validated our dipole shower's ability to reproduce rows 1 through 3 of Table 1. In this section we wish to test the shower's ability to reproduce the last row of Table 1 and the LC limit in Eq. (2.3), which is applicable across all the limits considered in the table and also when θ ij θ ik . To start we will test our dipole shower in the limit that q is soft whilst θ ij θ ik but θ iq ≈ θ jq > θ ij (i.e. we will compare against row 4 in Table 1). We describe these emissions as unordered in angle since they are produced in the shower with angles out of order; the k t and angle of these are emissions are not concurrently ordered. However, these emissions can still have a strong angular hierarchy allowing them to produce a LL, i.e. θ ij θ iq ≈ θ jq θ ik . The region of phase-space which has this hierarchy is highly restricted, due to the opposing k t ordering, but is nevertheless present and its mistreatment can induce a (small) LL error in some observables, for instance thrust [21]. At the end of this section we will check the crucial soft, wide-angle limit, where parton q is soft but all angles are unconstrained.
We will begin from Eq. (3.9), which was derived from our shower by only assuming parton q is soft. Employing θ ij θ ik allows us to replace P Now we take the limit that θ iq , θ jq θ ij , thus P (ij)k which gives rise to We should add to this the contribution where parton q is emitted first. Doing so gives (3.17) Comparing with row 4 of Table 1 we see an N −2 c suppressed error in the colour factor of the P (ij)k term. This error is due to the parton angles being disordered and is not present in an angular-ordered shower. However, in our dipole shower the disordered configuration is present and important, since it is required to get the correct wide-angle, soft physics beyond the two-jet limit.
Following the same logic as before, it is simple to show that our dipole shower includes the hard-collinear physics in Eq. (2.13) with the same N −2 c suppressed error as in Eq. (3.17). Finally, a good dipole shower should encode Eq. (2.3) and therein BMS evolution in the limit that the emission is soft. Starting from Eq. (3.8) and taking the leading colour limit: This is equal to which is equal to Eq. (2.3), modulo the use of k t instead of energy as the ordering variable, which does not hinder the logarithmic accuracy [28,29]. Hence the dipole shower correctly handles wide-angle soft radiation in the LC approximation. To go beyond the LC approximation generally requires amplitude-level evolution [20,25,[30][31][32][33].

Summary
In this section we have evaluated the accuracy at which our original dipole shower recreates the squared matrix elements summarised in Section 2. In summary, when parton q is emitted from parton i or k the matrix element is reproduced without error. When parton q is emitted from parton j there is an N −2 c suppressed error. We can look at different limits of the phase-space for partons j and q and evaluate the colour accuracy of our shower in each limit as follows: represents the re-arrangements of (a) that can be made in the limit θ ij 1. These diagrams correspond to those of an angular-ordered shower. The red factor is the N −2 c suppressed error produced by our original dipole shower.
1. θ ij 1: (a) θ jq θ ij : in this region an angular-ordered shower has full colour accuracy and our shower agrees with an angular-ordered shower (see rows 1 and 3 of Table 1).
(b) θ jq ∼ θ ij : in this region an angular-ordered shower cannot recreate the complete matrix element and our shower only guarantees LC accuracy in the soft limit. This region does not contain a strong angular hierarchy so at most can contribute a NLL and is suppressed further in event shape observables only sensitive to perturbations from the two-jet limit, for instance thrust.
(c) θ jq ≈ θ iq θ ij : in this region an angular-ordered shower has full colour accuracy and our shower currently lacks complete agreement with an angular-ordered shower beyond LC (row 4 of Table 1). This is the region we will address in Section 4.
2. θ ij ∼ 1: angular ordering cannot describe this region and our shower only guarantees LC accuracy. 7 7 In this limit (which is potentially subject to all manner of soft and non-global logarithms), it is difficult to make statements on the logarithmic accuracy of the shower beyond the leading accuracy of soft logarithms achievable through the BMS equation [34], which is embedded in the dipole shower approach. Though, with this in mind, Dasgupta et al. [15] have demonstrated LC NLL accuracy in non-global observables for dipole showers with carefully constructed global recoils and lab-frame based dipole partitionings. Our shower has both these properties and our fixed-order tests of the shower [1] are consistent with their results. Note that Dasgupta et al.'s definition of NLL accuracy encompasses NLL in the exponent but is also applicable to logs that do not resum into an exponential form such as non-global logs.
In Figure 4 we illustrate the origin of the N −2 c suppressed error: the erroneous factor is shown in red. The diagram in this figure is sufficient to enable us to read off the correct colour factor, and we make heavy use of this perspective in what follows.

Colour factors for emissions unordered in angle
In the previous section we computed the double-emission matrix elements squared corresponding to e + e − → qqgg, comparing the result from our dipole shower formalism with the relevant limits of the exact matrix element. We showed that, when the two emissions are strongly ordered in angle (with one emission collinear in the direction of one of the hemispheres), the matrix elements calculated from our dipole shower were correct except when a gluon is emitted with an angle larger than the opening angle of its parent dipole. In such a configuration, the coherent branching calculation would correctly assign a colour factor C F , whilst the dipole shower gives C A /2 (see Eq. (3.17)). At this order, we can correct the colour factor by replacing C i in Eq. (3.1) with a dynamic colour factor of where δ iq (δ ig ) is one when the parton i is a quark (gluon), and zero otherwise. We stress that this correction leads to the correct result only because our way of partitioning is able to encode angular ordering via (d cos θ aq ) dφ which ensures that the error is localized in the colour factor of Eq. (3.17). Our partitioning satisfies this requirement exactly. 8 It is not too difficult to generalize to higher orders, and the solution is particularly straightforward in the absence of g → qq branchings, which will be discussed at the end of this section (see also [35]). Figures 5 and 6 illustrate errors that occur in the case of three emissions. They highlight a key feature: the colour factor of the last emission is incompatible with coherence only when it is emitted at an angle larger than the angular extent of the colour charge distribution of the chain of partons leading to the emission. Figure 7 shows the generalisation to an arbitrary fixed order 9 . As a consequence of using a partitioning which defines a unique branching history of 1 → 2 transitions, the collection of partons in an event can be divided into m branches for an m parton hard process. Each branch contains one of the hard-process partons and the radiation emitted from it. Each parton in the branch can also be assigned a unique sub-branch consisting of Figure 5. Diagram (a) is generated by our dipole shower and is chosen because it contains an incorrect colour factor. Diagram (b) represents the re-arrangements of (a) corresponding to an angular-ordered shower. The red colour factor is the N −2 c suppressed error produced by our original dipole shower. Figure 6. A second possible ordering of angles that also leads to an incorrect colour factor. the parton and its "parental chain", see Figure 8. We only need to modify colour factors for gluons which cannot probe the largest angle in their sub-branch. We do this by extending the definition of C ij to where J is the hard parton in the sub-branch and L is the parton in the sub-branch emitted at the largest angle. θ LJ is the angle between L and J. 10 One should use these newly defined dynamic colour factors in both emissions and in the Sudakov form factors, i.e. so that the two are related by unitarity. 11 The computation of the dynamic colour factor grows at most linearly as the shower progresses and on average logarithmically. 12 In summary we have constructed a dipole shower which encodes the physics of QCD coherence just as in an angular-ordered shower. The resulting dipole shower, at LC, reproduces BMS evolution and, after using the CMW running coupling [37], will match the NLL-accurate dipole showers with global recoils discussed in [15]. In all, we expect our dipole shower to achieve complete LL accuracy in any observable for which an angularordered shower can also be used to resum LLs. In the case of e + e − → qq, the NLLs of some observables (i.e. thrust) do not directly depend on g → qq transitions 13 in which case 10 θLJ = π for an emission with a sub-branch of length 2. 11 All current dipole shower implementations [4,6,36] could directly employ our algorithm using existing methods such as the Sudakov veto algorithm. 12 The average sub-branch length for a multiplicity, n, of partons in the branch is n i=1 1 i ≤ ln n + 1. 13 Such transitions are restricted to secondary branchings, the remnants of which can be resummed into the CMW coupling and are otherwise rendered trivial by the angular-ordering constraint [2].  Figure 7. The generalisation of Figures 5 and 6 to an arbitrary fixed order. Cones i and J represent a unspecified number of parton branchings, each at angles smaller than θ LJ , which is the largest angle in q's sub-branch. As before, diagram (a) is generated by our dipole shower and contains an incorrect colour factor associated with the emission of q. Diagram (b) represents the re-arrangements of (a) corresponding to an angular-ordered shower. The red colour factor is the N −2 c suppressed error produced by our original dipole shower. our shower is accurate to NLL. Our methodology of assigning colour factors by mapping branching histories onto those of an angular-ordered shower could be generalised to assign the correct colour factors after including g → qq transitions. 14 However, because these transitions introduce more quark lines into the parton cascade, there would be the need to correct incorrect factors of 2C F . This would worsen the computational efficiency.

The effects of momentum conservation
In the paper so far we only briefly mentioned momentum conservation, which is vital for any implementation in an event generator, and needs to be treated very carefully. Bad implementations of momentum conservation have the potential to modify the phase-space boundaries of partons in the cascade or the matrix elements, leading to NLL errors [21]. In a dipole shower, emissions are on-shell and their momentum is typically expressed using three components: momentum longitudinal to the emitter, momentum transverse to the emitting dipole and momentum in the 'backwards' direction (collinear to the other parton in the dipole). A momentum map is used after an emission to ensure energy-momentum conservation in the shower by distributing 'recoil' across the partons in an event whilst keeping the partons are on-shell. In [1] we presented a momentum map with the idea of being as simple as possible whilst preserving the matrix elements computed by the shower. In the map, longitudinal recoil is trivially handled correctly (it is conserved between the emission and the parent parton as dictated by the dipole partitioning) and does not spoil anything. The other components are handled by a Lorentz boost and a global re-scaling of every momentum in the event after the emission. The emission kernels are invariant under both of these (as both z ab and dk ab ⊥ /k ab ⊥ are invariant under boosts and re-scalings). Thus only the phase-space is modified by the momentum map, not matrix elements. In Section 3.1 of [20] we showed that the changes to the phase-space due to recoil will generally not produce a log-enhanced term at O(α 2 s ) and that, for global two-jet observables such as thrust, artifacts in the phase-space from the recoil after iterated emissions produce terms beyond NLL. Alternative global momentum maps with similar constructions have also been studied in [15] where the NLL accuracy of the maps was demonstrated for a wide range of observables. The momentum maps in [15] were designed so that their action preserved key features of the Lund plane [39,40] (for instance preserving the separation between emissions on the plane). They have the added benefit of conserving 'backwards' components of momentum locally in a dipole, minimising the affect of the map on the phase-space available to partons in the shower. Any of these global prescriptions could be implemented into our shower without effecting the results in this paper.

Errors from in other dipole shower
In this section we want to emphasize the role of the dipole partitioning to our findings. To eliminate sub-leading colour errors, the partitioning function g ab must satisfy ab + negligible. (5.1) In Appendix B we discuss the term labelled 'negligible'; the remainder after azimuthal averaging when compared to the strict angular ordering result. Our dipole algorithm was carefully constructed to not produce such a contribution at all. Note that the demand of P [a] ab being proportional to a theta function cannot be satisfied with a zero remainder if g ab is positive definite and only zero at a finite number of points in the phase-space. On top of this, since P

[a]
ab has no dependence on the energies of the partons in the dipole, any partitioning that retains such a dependence after azimuthal averaging will result in a non-zero contribution remainder.
An interesting example to illustrate how wrong results can be obtained is that of Catani-Seymour (CS) dipole factorisation. The errors due to using a CS factorisation to construct the dipole partitioning have been previously noted in [21]. Here we give a complementary discussion. The CS partitioning contains both the issues described in the previous paragraph; the partitioning function is positive definite and has strong dependence on parton energies after azimuthal averaging. The partitioning that generates Catani-Seymour dipole factorisation is where η is the dipole-frame rapidity of parton q (η → ∞ as p q /E q → p a /E a and η → −∞ as p q /E q → p b /E b ). We must compute Using the basis p a = E a (1, 0, 0, 1), where Note D > sin θ ab sin θ aq for all momentum configurations. It is therefore easily shown that For all momentum configurations other than θ ab = π and E b = E a this results in ab . (5.6) Using this azimuthal averaging of the Catani-Seymour partitioning we can compute the azimuthally-averaged squared matrix element in the limit that emissions are strongly ordered in angle and energy: As k ab ⊥ ≈ k ca ⊥ ≈ E q θ aq in this limit and since energies are strongly ordered, the k t ordering theta functions are saturated and we find: We can subtract this from the correct result (for rows 1 through 3 of Table 1) to find the error: kj ) Of course the error vanishes if C F = C A /2. The error becomes large when E j E i ≈ E k . In this limit (1 − cos θ jb ) (1 − cos θ jq ) | cos θ jb − cos θ jq | , (5.10) when θ ij ≈ θ jq and θ jk ≈ θ jq 15 , and where b = i, k. Also in this limit 11) once again this is only valid when θ ij ≈ θ jq and θ jk ≈ θ jq (equivalently θ iq ≈ 0 and θ kq ≈ 0). Thus Note that the C A /C F piece contains a non-cancelling collinear pole when θ jq → 0 and so is capable of generating logarithms in observables that probe secondary emissions even in the limit of a strong angular hierarchy, where θ jq θ ij θ (ij)k , since the numerator of ji goes as O(θ 2 ij ) whilst the numerator of P [j] ji goes as O (1). Also note that this error cannot be fixed by using the dynamic colour factors C iJ (θ iq , θ LJ ) since in the limit we are considering the dynamic colour factors reduce exactly to the usual colour factors already present in Eq. (5.12).
We can also compare the error made using a CS-partitioned dipole shower with row 4 from Table 1. Here we find This error is potentially LL, since with a strong hierarchy in emission energies and angles the functions W are singular and so capable of generating double logarithms. Of course this too vanishes if C F = C A /2. If the dipole shower instead used colour factors C iJ (θ iq , θ LJ ) this limit would be improved since the error would instead become (5.14) 15 In the region where θij ≈ θjq terms in Ej/Ei are not negligible as they screen the divergence θij = θjq.
However, this improvement may not extend to higher orders since θ LJ as computed with the CS dipole shower branching history will not necessarily equal θ LJ as computed from a branching history matched to the angular-ordered description. This problem, combined with Eq. (5.12), is sufficient for us to assert that CS dipole showers employing the dynamic colour factors C iJ (θ iq , θ LJ ) will still be subject to LL errors in some observables that angular-ordering can completely describe at LL.

Conclusions
We have performed a fixed-order cross-check of the dipole shower presented in [1] and shown that the shower performs as it was designed to: the shower inherits its handling of collinear radiation from an angular-ordered shower whilst improving over angular ordering in the case of the leading colour, wide-angle soft radiation. We also highlight a limitation of our original approach, showing how the dipole shower will not assign correct colour factors to emissions disordered in angle, though they will be correct at leading colour. We then introduced a new method for correcting these colour factors. The new method is efficient: the computation time on average grows logarithmically with parton multiplicity. Using this method, our shower will match the LL accuracy of an angular-ordered shower in cases where an angular-ordered shower has LL accuracy. When enhanced with the CMW running coupling [37], our shower will include all leading logarithms and leading-colour, next-toleading logarithms in the two-jet limit for continuously-global observables. As it stands, the shower will not be capable of the full-colour NLL resummation of global observables, due to the absence of full colour g → qq transitions. These transitions could be included as described in [1] but would generate sub-leading colour NLL errors: however, they could be included at full colour by extending the methods outlined in Section 4.

A Drawing planar diagrams at arbitrary order
In this appendix we demonstrate that the planar diagrams representing re-arrangements of our dipole shower into an angular ordered shower, in Figure 4, are not just a feature of our dipole shower at O(α 2 s ) but rather can continue to be used at higher orders if we continue to assume that the branching history produced by our shower has a strong hierarchy in angles. We do not assume a hierarchy in angles that is concurrently ordered with their k t . At a scale k ⊥ , a given n-parton state produced by our dipole shower has a weight at a point in the n-parton phase-space dS n (k ⊥ ). We consider dressing this state with one further gluon, q, produced by the shower. This gives an (n + 1)-parton state: where c.c. means a and b are colour connected in the n-parton state and J is the hard parton that initiated a's branch. All other symbols have the same definitions as in the previous sections. We can azimuthally average exactly as in Section 3, and find ab P d a→aq (z ab )dz ab Just as we have already demonstrated at O(α 2 s ), the weight assigned to the (n+1)-state after azimuthal averaging uses the same LC emission kernels as an angular-ordered approach. We can make this very explicit by exchanging the sum over colour lines with a sum over parton indices. To illustrate this, at LC we find ab P d a→aq (z ab )dz ab δ(k ab ⊥ − k ⊥ ) We can exchange the non-singular dependence on b and c in Eqs. (A.2) and (A.3) with that of J (or the other hard parton J if either b or c are in the opposing hemisphere). Similarly, for the non-singular region θ aq θ aJ we can exchange the dependence on a with J so that θ aq ≈ θ Jq . Thus, just as in an angular-ordered framework, emissions are generated with a weight C P [a] aJ P d a→aq (z aJ ) when they can probe the jet and C P

[J]
JJ P d J→Jq (z JJ ) when they cannot. At LC, C = C A /2 when a is a quark and, when a is a gluon, C = C A if q can probe the jet (determined by the angular ordering constraint embedded in P [a] aJ 16 ) and C = C A /2 when q cannot. It is these properties that our planar diagrams are defined to encapsulate, validating their usage at arbitrary higher orders. The planar diagrams led us to define C aJ (θ aq , θ LJ ) so that the sub-leading N c terms are included in accordance with Figure 7.

B Current limitations of our dipole shower
An important part of our dipole shower is its partitioning. However, the form of our partitioning, g ab (defined through Eqs. (3.5) and (3.7)), might cause complications in a computational implementation of our shower. In this appendix we will discuss the issues and possible solutions.
A dipole shower is fully differential in the parton phase-space and so emits by sampling from the distribution g ab w ab to populate a 3-dimensional phase-space for each parton. 17 However, g ab w ab has two undesirable properties (illustrated in Figure 9): firstly g ab w ab is negative in some portions of the emission phase-space, introducing negative weights into the shower; secondly g ab w ab contains an integrable singularity when θ aq = θ ab < π and φ (a) q = 0 (i.e. q is in the plane of partons a and b). Both of these features can be handled in a modern dipole shower: the Herwig dipole shower already contains all the necessary machinery [6], as do others [42]. However, both features will hinder numerical convergence. Fortunately the two features counter balance each other: g ab w ab is most negative when θ aq = (1 + ε)θ ab , for ε 1 whilst strictly positive, and φ (a) q = 0. The negative weights and integrable singularity are linked such that, when θ aq = θ ab , g ab w ab azimuthally averages to a well behaved quantity, A simple solution to the two issues would be, in regions bounded by θ aq = θ ab ± δθ (for δθ/θ ab 1), to sample emissions according to the azimuthally averaged distribution, Eq. (B.1). This would entail sampling emissions from a discontinuous distribution but would alleviate the undesirable features whilst only introducing a power correction in δθ to azimuthal correlations in the shower.
Alternatively, one might use an alternative partitioning,g ab , free from negative weights and integrable singularities, that satisfies (d cos θ aq ) dφ Thisg ab would be suitable for use with our proposed dynamic colour factors and retain our shower's accuracy concerning LC NLL physics. It is possible that a pre-existing partitioning employed by another parton shower might already achieve this. We have demonstrated that the Catani-Seymour partitioning [22] does not satisfy this requirement but there are others on the market that we have not tested [15,[42][43][44]]. An acceptable partitioning should at least satisfy (d cos θ aq ) dφ = 0, i.e. q is in the plane of the dipole (a, b). For θaq < θ ab , q is 'inside' the dipole. The negative weights and integrable divergence are present when θaq = θ ab . Figure 9. Graphs of g ab w ab as a function of θ ab , θ aq , φ (a) q as measured in the lab frame.