Groomed jet mass as a direct probe of collinear parton dynamics

We study the link between parton dynamics in the collinear limit and the logarithmically enhanced terms of the groomed jet mass distribution, for jets groomed with the modified mass-drop tagger (mMDT). While the leading-logarithmic (LL) result is linked to collinear evolution with leading-order splitting kernels, here we derive the NLL structure directly from triple-collinear splitting kernels. The calculation we present is a fixed-order calculation in the triple-collinear limit, independent of resummation ingredients and methods. It therefore constitutes a powerful cross-check of the NLL results previously derived using the SCET formalism and provides much of the insight needed for resummation within the traditional QCD approach.


Introduction
The study of the substructure of jets is now an established and highly active component of LHC phenomenology. Following early work on the subject [1,2], the true power of substructure methods for new particle discoveries in the boosted regime was first revealed over a decade ago [3]. This in turn led to an explosion of interest in the subject and rapid development of tools for tagging and grooming jets as well as their direct exploitation in experimental searches and studies at the LHC (see [4][5][6][7] for reviews and further references). a e-mail: dpa@m.scnu.edu.cn b e-mail: mrinal.dasgupta@manchester.ac.uk c e-mail: basem.el-menoufi@manchester.ac.uk d e-mail: mguzzi@kennesaw.edu e e-mail: jack.helliwell@manchester.ac.uk Jet substructure has also proved to be a fertile field for the development of calculations and concepts in QCD, aimed at improving the theoretical understanding of substructure. For a review of such work we refer the reader to Refs. [8,9] and references therein. An especially important development in this context has been the invention of methods and observables that can be perturbatively calculated to high precision in a hadron collider environment and receive only modest non-perturbative corrections. This is a significant development for jet based studies which are generally subject to large non-perturbative uncertainties from hadronisation and the underlying event (in addition to pile-up), as well as perturbative issues such as the presence of non-global logarithms where the resummation accuracy is still generally limited to leading-logarithms and the use of the leading N c limit [10,11].
A well-known set of observables, that eliminate some and ameliorate other issues obstructing precision for hadron collider jets, are those involving groomed jets where grooming is implemented through the modified mass-drop tagger (mMDT) [12,13] or its subsequent generalisation SoftDrop [14]. The use of a z cut parameter as a threshold below which soft particles can be groomed away, eliminates non-global logarithms (NGLs) that otherwise appear in the ungroomed jet mass. Although NGLs still appear as logarithms in z cut , the phenomenologically relevant values of z cut ∼ 0.1 imply that resummation of logarithms of z cut is not strictly necessary. Moreover grooming via SoftDrop also considerably reduces the impact of non-perturbative effects relative to the ungroomed jet mass [12], virtually eliminating the underlying event at high p T and leaving a modest hadronisation contribution amenable to analytic studies. These developments mean that theoretical calculations for jet observables with mMDT/SoftDrop grooming can be carried out to high precision and accurate comparisons can be made to LHC data, which is a program that has already been successfully estab-lished with several recent theoretical studies, experimental measurements and phenomenology [15][16][17][18][19][20][21][22][23][24][25][26][27][28].
In this paper we focus on the mMDT (equivalently Soft-Drop with β = 0) jet mass distribution in the region of small ρ, ρ z cut , where ρ is the normalised jet mass which for hadron collider jets is generally defined as ρ = m 2 p 2 T R 2 , with m 2 the squared jet-mass, p T the jet's transverse momentum and R the jet radius parameter. In addition to the removal of non-global logarithms and a considerable reduction in non-perturbative contributions mentioned above, the mMDT is unique in that the leading-logarithmic resummed result, for the jet mass ρ, is single-logarithmic i.e. contains only terms α n s ln n ρ where the logarithms in ρ are of pure collinear origin. Such a structure is of course well-known for collinear evolution, both space-like and time-like, of partons via the DGLAP evolution equations for structure functions and fragmentation functions and for the small R limit of jet cross-sections [29] but is exceptional for observables like jet masses which in the ungroomed case receive double logarithmic terms at each perturbative order.
The leading-logarithmic (LL) resummed result for the mMDT jet mass, amounting to a resummation of the singlelogarithmic terms in ρ, was first derived in Ref. [12]. The main result of that paper was written in the limit of small z cut i.e. neglecting terms power suppressed in z cut . A result including finite z cut corrections was also presented in the appendix of Ref. [12] where complications arising from the leading parton flavour changing were accounted for through a matrix structure for the resulting resummation. An NLL result, in the small z cut limit, was computed for the first time in Ref. [15] and most recently a resummation to NNLL accuracy in the same limit has been carried out [23]. 1 Phenomenological studies have also been carried out for mMDT/SoftDrop jet masses where resummed calculations, supplemented by non-perturbative corrections, have been compared directly to experimental data. In the mMDT context the resummed calculations refer to LL calculations including finite z cut effects and matching to full next-to-leading order (NLO) calculations [16,17] as well as calculations with NLL resummation but without finite z cut corrections and matched to LO calculations [15]. Generally good agreement with data has been observed in each case and we refer the reader to Refs. [26][27][28] for further details.
In the present article we take a fresh look at the NLL structure of the mMDT jet mass where previous work on the subject has been entirely within the context of Soft-Collinear Effective Theory (SCET) calculations [15,19,21,23,24]. 1 Our description of the resummation accuracy for the mMDT is consistent with the original reference Ref. [12] while more recent work uses a convention inspired by the general structure of SoftDrop jet mass, where single logarithms are referred to as NLL, in spite of the absence of double logarithmic terms. In this convention the work of Ref. [15] is NNLL while that of Ref. [23] achieves NNNLL accuracy.
The original work of Ref. [15] recycled SCET calculations for soft functions with an energy veto [30,31] and calculated a correction for C/A jet clustering [32,33] to derive the result for the mMDT. In the current paper we shall approach this question from the viewpoint of QCD calculations by directly addressing the collinear nature of the ln ρ terms in the jet mass distribution. Thus while a leadinglogarithmic resummed result can be derived in a strongly ordered picture with collinear emissions widely separated in angle, the NLL structure with terms α n s ln n−1 ρ stems from lifting the strong ordering on a pair of emissions. At order α 2 s this leads us to consider two emissions (i.e. three partons) within a jet, which are at small angles θ 2 ρ z cut 1 but with no relative strong ordering between them. Such configurations are related to an NLL DGLAP evolution picture and the QCD matrix element can be approximated by triplecollinear splitting functions [34][35][36]. In the strongly-ordered limit the triple-collinear splitting functions reduce to a product of LO splitting functions thereby restoring the LL picture.
There has been much recent interest and work towards the incorporation of higher-order splitting functions, including the triple-collinear and double-soft limits, in the context of parton shower algorithms [37][38][39]. With the emergence of new research clarifying and pushing forward the logarithmic accuracy of parton showers [40,41], it is also of interest to understand in more detail the connection between higherorder splitting functions and the logarithmically enhanced terms for QCD observables.
With the above points in mind, in this paper we compute the O α 2 s ln ρ NLL term, for quark initiated jets, using the triple-collinear splitting functions and phase-space. Our calculation is a pure fixed-order computation, with the main approximation being the use of the collinear limit. In particular we do not make use of known ingredients from any resummation approach. We show that the expected leadinglogarithmic terms emerge from our calculations along with recovering the proper argument of the running coupling in the soft limit and the constant K related to the emergence of the CMW or "physical" coupling [42][43][44]. Apart from these standard ingredients we also derive the collinear NLL terms finding agreement with the SCET results [15,19,24]. Given our derivation directly from the QCD matrix-elements and independent of any previous SCET work, this constitutes a strong completely independent check on the main results involved in the NLL resummation for mMDT jet mass. We shall also shed further light on the collinear structure we find here in terms of standard QCD resummation elements which should, we hope, enable the resummation of the mMDT jet mass also within a QCD resummation formalism e.g. through a suitable modification of the NNLL ARES method [43,45] or directly by embedding the triple-collinear splitting within a strongly-ordered parton cascade. We note that calculations based on triple-collinear splitting functions in the context of jet physics have also been carried out in Refs. [46,47]. For early work on integrals of triple-collinear splittings in the context of initial state space-like splittings we refer the reader to Ref. [48].
Although we carry out our present calculation in the context of e + e − annihilation precisely as in Ref. [15] to control process dependent pieces of the calculation, the final NLL results we extract derive purely from the collinear limit and hence are universal. We confine our attention here to quark initiated jets though similar calculations can be performed for gluon jets. Also, while for the purposes of the present paper, we shall work in the small z cut approximation purely for simplicity, we stress that our approach can be extended to obtain the finite z cut corrections at the NLL level. Although finite z cut NLL corrections are expected to be numerically small they should be of increased relevance in the context of the recent NNLL calculation performed in the small z cut limit [23].
The current paper is organised as follows. We start in section 2 by defining the normalised jet mass observable ρ, for which we use the same e + e − definition as in Ref. [15] i.e. study the groomed heavy hemisphere mass, and write down the leading-logarithmic resummed formula with running CMW coupling. In section 3.1 we compute the leadingorder result in the small ρ limit and identify the ρ independent constant term C 1 α s that multiplies the Sudakov factor in the LL formula. We then expand the LL result to order α 2 s and note that while this expansion also produces terms that are NLL in ρ, such terms originate within an LL context, being still related to the strongly-ordered picture. In section 4 we carry out our α 2 s calculations for the jet mass differential distribution, ρ dΣ dρ , considering various real-emission configurations including emissions in opposite hemispheres and those in the same hemisphere. We compute separately the contributions from gluon emission with a C 2 F colour factor and from gluon decay with C F C A , C F T R n f and C F (C F − C A /2) terms and combine the divergent real-emission results with those from virtual corrections. At this stage we compare our results to the expansion of the LL formula and identify the process-independent NLL terms. In section 5 we compare our results for each colour channel to the SCET results finding agreement in each case. We finally comment on the nature of our results within the QCD resummation context and briefly mention some prospects for further work. Various formulae relevant to the derivation of our results are listed in the appendices.

Observable definition and leading-log resummation
We are interested in the jet mass distribution of a QCD jet after the application of grooming via the modified mass-drop tagger (mMDT) or equivalently SoftDrop (β = 0). In the current article, for reasons of convenience, we shall work in the context of e + e − collisions, though the NLL pure collinear terms we will eventually extract are process independent and hence apply also to jets at hadron colliders.
We will compute the standard heavy hemisphere jet mass observable extensively studied in e + e − collisions but with the modification that we compute the heavy jet mass after running the mMDT procedure on the particles in each hemisphere. This was also the observable studied in the first mMDT NLL resummed calculation performed in Ref. [15]. In general one may separate the event into two hemispheres in different ways e.g. by clustering to two jets as in Ref. [15] or as is traditional by using the thrust axis. At the level of our calculations and for extracting the terms we seek, we are insensitive to the precise details since in the soft and/or collinear limit we can take the hemispheres to be defined by the directions of the initial quark-anti-quark pair.
The mMDT algorithm, as applied to e + e − collisions, involves declustering particles in a hemisphere of an e + e − event using the Cambridge/Aachen (C/A) algorithm, as one does for a hadron collider jet. We produce, at each stage, two branches i and j and we require If a branching fails this requirement we reject the softer branch and proceed with declustering the harder branch until the condition passes or we end up with a massless hard parton. We then define where we select the larger of the left or right hemisphere squared invariant masses (M 2 L and M 2 R respectively) and normalise to (Q/2) 2 which corresponds to the squared energy of a hemisphere in the Born limit.
We shall work in the formal limit ρ z cut 1, which means that we will examine the structure of ln ρ enhanced terms, but shall neglect power corrections in z cut therin. A leading log (LL) resummation formula for the mMDT jet mass distribution, based on an independent emission picture with emissions strongly ordered in angles, was first provided in Ref. [12] (see Eq. (7.2) of Ref. [12] for the result). While that result applies directly to a jet produced in hadron collisions, it can be easily modified to the case of the heavy groomed hemisphere mass in e + e − annihilation. We express our leading-logarithmic resummed result in terms of the integrated distribution, Σ (ρ) = ρ 0 where we defined the splitting function p qq (z) = (1 + z 2 )/(1 − z) and σ 0 is the Born cross section.
The above result modifies the result of Ref. [12] by replacing p 2 T by Q 2 /4 and is written in terms of the splitting function p qq (z) rather than p gq (z) so that the variable 1 − z in the above result corresponds to the variable z in Eq. (7.2) of Ref. [12]. We have also inserted an additional factor of two in the Sudakov exponent which takes into account the fact that we are considering the result due to two hard partons in the left and right hemispheres respectively, rather than just a single parton initiating a jet. This is accounted for by our definition of the splitting function which has a factor of two relative to that defined in Ref. [12]. The above result is labelled LL since it correctly resums terms α n s L n , with L = ln ρ, i.e. the leading logarithms in ρ. Additionally a fixed-coupling calculation of the exponent in Eq. (3) reveals a term of the form α s ln 2 z cut , which although subleading in ρ, correctly resums terms α s ln 2 z cut n [12].
In Eq. (3) we have defined the coupling in the CMW or "physical" scheme [42][43][44] and have set the argument of the running coupling to be the transverse momentum squared of a soft and collinear emission relative to the direction of its hard parent parton, expressed in terms of the invariant mass ρ and the energy fraction of the emitted gluon Q 2 (1 − z), since we have k 2 t = (1 − z)ρ Q 2 /4. 2 Both the specification of the CMW scheme and the precise details of the argument of the running coupling (beyond the fact that it scales with ρ) in fact produce terms that are formally NLL in ρ but with logarithmic enhancements in z cut . These terms are an intrinsic part of the leading-logarithmic resummation framework and hence naturally belong in our LL formula. In a similar spirit we have also included in Eq. (3), a ρ independent constant order α s coefficient C 1 , which is process dependent and on physical grounds factorises from the leading-logarithmic Sudakov resummation. On expanding the exponent, multiplication by the C 1 term results in an α 2 s L term. While this term should be reproduced by our order α 2 s calculations, it is process dependent and unrelated to the collinear NLL structure we ultimately seek to extract. Hence the explicit identification of C 1 is needed, to account for the role of this piece in our final result.
We then expect that our eventual order α 2 s result should contain all terms generated by the expansion of (3) and additionally all terms of genuine NLL origin (i.e. unrelated to the strongly-ordered in angle LL dynamics). These terms should arise from collinear physics, and be independent of z cut and of the specific hard process. 2 Recall that in the soft-collinear limit ρ = m 2 /(

Leading order calculation and expansion of LL result
Throughout this article we shall work in the resummation region where ρ z cut and hence ignore the presence of a transition point at ρ ≈ z cut , beyond which the groomed jet mass result becomes coincident with the plain ungroomed mass. We shall consider for simplicity that z cut 1, so that we can neglect powers of z cut , although this is not a requirement on the validity of our triple-collinear calculations.
Our first step will be to determine the ρ independent coefficient C 1 that appears in Eq. (3). We do this in the following subsection by performing a leading-order calculation in the soft and collinear limit.

Leading-order calculation
Here we derive, in the small ρ and z cut limit as defined before, the order α s contribution to Σ (ρ). At this order we have to consider a single real emission and the one-loop virtual correction to qq production. To handle divergences in the real emission calculation, we perform the calculation in conventional dimensional regularisation (CDR) in d = 4 − 2ε dimensions, and combine with the virtual correction before taking ε → 0.
Since we are interested in the region of small jet mass and small z cut , we can work in the soft and/or collinear limit, perform the calculation for a single hemisphere and then double the result. Considering the gluon to be in the same hemisphere as the quark (for example) and applying the mMDT, we have a situation where the z cut condition passes or fails. If it passes then one obtains a non-zero hemisphere mass, while if it fails the hemisphere mass vanishes.
When the z cut passes, we are considering a pair of relatively energetic partons which produce a small jet mass, implying that the angle between partons is small, θ 2 ∼ ρ z cut 1, which allows the use of the collinear approximation. In this region the emission probability is just the q → qg splitting function in 4 − 2ε dimensions p qq (z, ε) = p qq (z) − ε(1 − z), and the result follows by applying the standard d-dimensional collinear phase space [49]: where the label r 1 indicates the first of three distinct realemission terms that we encounter, and we have provided an explicit factor of 2 to take account of both hemispheres, which shall apply to all the real-emission terms we compute in this section. The constraint arising from restricting the jet mass is written in terms of a step function involving ρ, while the z cut condition removes any soft singularity in the integral and we obtain just a collinear pole alongside finite corrections. We have not explicitly written, for brevity, the dependence on z cut in the argument of Σ 1,r 1 which is left implicit here and throughout the article.
Next we turn to the region where the z cut condition fails, the softer particle is groomed away, and we obtain a massless hemisphere. This can happen when the quark goes soft for z < z cut or when the gluon goes soft for z > 1 − z cut , and in either case the jet mass vanishes so that we can replace the step function in Eq. (4) by unity. The region corresponding to a soft quark also contributes only a collinear divergence and finite terms which are suppressed by powers of z cut which we neglect. Hence retaining just the collinear divergence, corresponding to a pole in ε produced by the θ 2 integral, we can write this contribution as where we performed the θ 2 integral and discarded all finite terms owing to their suppression with z cut . Accordingly we have dropped the ε dependence of the constant prefactor multiplying the integral, since this only leads to finite power suppressed in z cut corrections. Lastly there is the region z > 1 − z cut corresponding to a soft gluon emission, which is removed by grooming. As there is no longer a constraint on the angle, one has to replace the soft-enhanced part of the splitting function with the full eikonal function to account for the wide-angle region. For soft-regular pieces, however, one merely needs the collinear pole as any finite contributions are power suppressed in z cut . Hence, using p qq (z, ε) = 2/(1−z)−(1+z)− ε(1 − z), to separate the soft enhanced and regular terms of the splitting function, we write Σ 1,r 3 = Σ 1,r 3 soft + Σ 1,r 3 coll. , where explicitly In the expression for Σ 1,r 3 soft the upper limit of the cos θ integral corresponds to θ = π/2 for the boundary of the hemisphere under consideration, while in the expression for Σ 1,r 3 coll. we have retained only the singular contribution as finite corrections vanish with z cut .
Our final ingredient is the well-known virtual correction to qq production (see e.g. [49]): We express our results in terms of the renormalised coupling in the MS scheme, α s (µ 2 R ), using the relation where we have the standard MS factor and we choose µ R = Q/2. 3 Carrying out the necessary integrals we obtain the result for the real emission term Σ 1,r = Σ 1,r 1 + Σ 1,r 2 + Σ 1,r 3 soft + Σ 1,r 3 coll. , Combining with the virtual correction we obtain the leadingorder result, Σ 1 (ρ) = Σ 1,r + V (ε): where the argument of the running coupling reflects our choice of µ R . We have written the above result separating the contributions that arise from expanding the LL formula (i.e the ln 1/ρ and ln 2 z cut contributions that are associated to the fixed-coupling Sudakov exponent) from the contributions we shall associate to C 1 . This allows us to identify We note that defining the observable as v = ρ/4, corresponding to a normalisation to Q 2 , would result in the elimination of the term proportional to ln 2 and hence the ln z cut dependence from C 1 , but the resummation of ln z cut terms is beyond the scope of our results.

LL result at order α 2 s
On performing the full order α 2 s calculation, to our NLL in ρ accuracy, we should recover all terms produced by the expansion of Eq. (3) in addition to terms that are unrelated to the leading logarithmic structure, which will then act as one of the checks on our results. To this end we report below the expansion of Eq. (3) for the differential distribution ρ dΣ dρ (to correspond to the calculations of the following sections). The leading-order result is given by (for our choice of µ R = Q/2), To obtain the order α 2 s result one needs to compute the Sudakov exponent with running coupling and switch from the CMW scheme to the MS scheme, which gives The above result contains a term which is LL in ρ originating from the exponentiation of the leading-order result. It also contains NLL in ρ terms, corresponding to ρ independent terms in ρ dΣ 2 /dρ generated by an interplay of the LL exponent with C 1 and by fixing the scale and scheme of the running coupling. In particular the ln 2 z cut terms in the C F C A and C F T R n f channels derive from the 1 − z factor in the argument of the running coupling while the ln z cut term in the same channels is generated by changing from the CMW scheme to the MS scheme as can be seen through their coefficient, proportional to K = 67 18 − π 2 6 C A − 10 9 T R n f . In addition to these terms, the 1 − z factor in the argument of the running coupling and the CMW coefficient K applied to the full splitting function, rather than just its soft enhanced piece, are also responsible for producing z cut independent NLL terms, which go beyond the strict jurisdiction of the LL formula. The full set of such terms will be identified through the calculation we perform here and can be properly accommodated within a consistent NLL resummation formula.
In the next sections we shall derive the full result at order α 2 s through to NLL accuracy, compare it to the expectations from Eq. (16) and derive the z cut independent NLL corrections.

NLL at O α 2
s : the triple-collinear limit At LL accuracy, for the mMDT jet mass, we have a picture of successive collinear parton branchings which are strongly ordered in angle with each branching being described by a leading-order (LO) 1 → 2 splitting function. Thus at order α 2 s , the real emission matrix-element simply involves a product of two LO splitting functions. To obtain NLL accuracy, at order α 2 s , one has to consider three partons that are comparably collinear i.e. the opening angle between any two partons is small θ 2 i j 1 but there is no strong ordering so that θ 2 12 ∼ θ 2 13 ∼ θ 2 23 . Such configurations are described by the 1 → 3 collinear splitting of an initial parton, and the matrix-element involves triple-collinear splitting functions. In the strongly-ordered limit, the triple-collinear splitting functions reduce to a product of LO splitting functions (in general after azimuthal averaging) thus restoring the LL picture.
For our current calculations, the relevant functions are the unpolarised triple-collinear splitting functions for a quark (or anti-quark) initiated 1 → 3 splitting, denoted P abq for a splitting q → qab with a and b representing parton flavours, which were first computed in Refs. [34][35][36], and are listed in the appendix. For the q → ggq splittings there is both a gluon emission contribution with a C 2 F colour factor, and a gluon decay contribution with a C F C A colour factor. The q → qqq splitting arises from gluon decay and has a C F T R n f colour factor as well as a contribution from an identical particle interference contribution, involving quarks of the same flavour in the final state, which has a colour factor C F (C F −C A /2) i.e. vanishes in the leading N c limit. Identical considerations apply for the decay of an initial anti-quark. In the following subsections, we consider the gluon emission and decay contributions in turn.

Gluon emission contribution
Here we study the emission of two gluons from the initial qq system, associated with a C 2 F colour factor. The emitted gluons can either be in the same or in opposite hemispheres, with the latter case being simply related to the leading-order calculations we have already performed. We deal with each contribution in turn below.

Emissions in opposite hemispheres
Consider a gluon emitted in each of the "right" and "left" hemispheres containing the quark and anti-quark respectively. Let us assume that the right hemisphere is heavier after grooming and that its groomed jet mass is ρ. This implies that the branching in the right hemisphere must pass the z cut requirement corresponding to 1 − z cut > z > z cut and that it must set a mass ρ, while in the left hemisphere the mass must be below ρ for it to be lighter, and hence the grooming can either retain both or remove one of the two particles. The fact that in the right hemisphere the grooming passes, coupled with the limit we are working in, with ρ z cut , allows us to use the collinear approximation, so that the branching in the right hemisphere factorises from the dynamics of the left hemisphere. Here we do not require the triple-collinear splitting functions, as the emissions in opposite hemispheres are well separated.
For the left hemisphere the constraint on the mass to be below ρ simply gives us Σ 1,r /2 with Σ 1,r the real emission result already computed in the previous section (see Eq. (12). For the right hemisphere the distribution can be simply calculated in the collinear limit using the LO splitting function and the collinear 1 → 2 phase space. Finally a factor of two accounts for the case when the left hemisphere is heavier after grooming.
For compactness, here and in the sections below, we define the quantity F (ρ), where Σ 2 is the order α 2 s contribution to Σ (ρ). The result for the emissions in opposite hemispheres can then be written as a product of two leading-order factors 4 : Note that although the LO jet mass distribution for the right hemisphere is a finite quantity, we have retained its ε dependence in the result above, since Σ 1,r contains double and single ε poles and finite terms (in the limit ε → 0) are generated by the ε expansion.

Emissions in the same hemisphere
When two emissions, i.e. three partons, are in the same hemisphere one has to consider the action of the mMDT taking into proper account the Cambridge/Aachen clustering sequence. This involves considering different angular regions Fig. 1 An illustration of the parameterisation used for kinematic variables in the gluon emission process, relevant to the triple-collinear limit calculation discussed in the main text.
where the two emissions can be clustered separately to the quark (or anti-quark according to the hemisphere in question) or are first clustered together and then clustered to the quark. It proves convenient to divide the calculation into two pieces: a first piece where in all angular regions we apply the mMDT as if the gluons are clustered separately to the quark and a second piece which restores the correct action of the mMDT in the angular region where the two gluons are clustered. We are then led to consider the following distinct cases: -Larger-angle gluon passes z cut : Neglecting the clustering of the two gluons, the mMDT declustering produces two branches, consisting of the larger-angle gluon and a massive branch with the quark and the smaller-angle gluon. When the first declustering passes the z cut condition all three partons are retained. The angle between the branches is small, being set by θ 2 ∼ ρ/z cut and hence the three partons are constrained to be within a small angular region and the triple-collinear limit generates the full result. The result for the real emission calculation will be divergent due to the smaller-angle emission becoming soft and/or collinear and will contain double and single ε poles. -Larger-angle gluon fails z cut : Another relevant situation is that the larger-angle gluon is soft and hence the first declustering fails the z cut condition. 5 In this case the soft gluon is groomed away while the tagger then declusters the second gluon and we require the second declustering to pass the z cut condition to obtain a massive hemisphere.
In this case only the smaller-angle gluon is constrained by the jet mass to be collinear to the quark, while the first emission can be at a large angle. Hence we need to modify the triple-collinear splitting functions to match the correct soft large-angle emission pattern for the first gluon, precisely as we did in section 3.1. The result is divergent due to the soft divergence produced by the larger-angle gluon, though it contains only a single pole in ε. -Correction for gluon clustering: In the region where the angle between the two gluons, θ 12 , is the smallest an-gle the gluons are clustered in the C/A algorithm. Hence the first declustering produces a massive branch with the two gluons and a massless branch i.e. the quark. If the quark is soft the z cut condition can fail and the tagger recurses down the massive two-gluon branch. However such configurations with a soft quark are suppressed by powers of z cut and can be ignored consistent with our intended accuracy. The two-gluon branch must also pass the z cut as grooming it away would lead to a massless jet. Hence we can always take the first declustering to pass the z cut condition so that all three partons are retained, implying that the triple-collinear limit is once again the relevant one. In the angular region where the gluons are clustered, we shall subtract the contributions already included in the first two scenarios described above, and shall add the correct constraints just discussed. The difference between the correct and subtraction terms is finite and can be calculated in four dimensions.
For calculations in the triple-collinear limit, we work in terms of the energy fractions z i , defined wrt the initiating parton's energy, and which satisfy ∑ i z i = 1, and the angles θ i j between any two partons i and j, such that θ i j 1. The triple-collinear phase space in 4 − 2ε dimensions may be written as 6 [50] To make contact with the LO splitting functions in the strongly-ordered limit, it is useful to parameterise the variables z i in terms of variables z and z p as depicted in Fig. 1. Our general method for integrating the triple-collinear functions and extracting its divergences is briefly described in Appendix B. It involves systematic subtraction of soft and collinear divergences, via a series expansion around divergent limits, to obtain the pole structure and a pure finite contribution which we integrate numerically in four dimensions. Our results shall thus be partly analytical (stemming from performing an ε expansion of the coefficients of the pole terms) and partly numerical in nature.
We first provide the details for results neglecting the clustering of gluons starting from the contribution where the larger-angle emission passes the z cut , which we label F pass (ρ, ε). We take θ 13 to be the larger angle and hence for the first declustering to pass the z cut condition we have that 1 − z cut > z > z cut . The smaller angle gluon is not examined for the z cut condition and can be arbitrarily soft and/or collinear leading to divergences in both limits i.e. from z p → 1 and θ 23 → 0. The relevant splitting function and prefactor is specified by Eqs. (A.1), (A.3) and we obtain where i j is the squared invariant mass of the three parton system and δ ρ (1, 2, 3) is an abbreviated notation for the condition that the normalised hemisphere jet mass ρ involves all three partons i.e. the condition where the factor of ρ in front of the delta function takes care of the fact that we are studying the logarithmic derivative ρdΣ /dρ. We shall also use the notation Θ z cut (a|b) to denote the condition that two branches a and b, made up of one or two partons, pass the z cut condition. Partons not included in a and b fail the z cut condition and are removed by grooming. Thus Θ z cut (1|23) in Eq. (21) indicates that both branches i.e. gluon 1 and the massive branch with gluon 2 and the quark (with index 3) pass the z cut condition. More explicitly we have Θ z cut (1|23) = Θ (z < 1 − z cut )Θ (z > z cut ), amounting to simply a cut on z.
In Eq. (21) we have also introduced the renormalised MS coupling α s (µ 2 R ) via the use of Eq. (11) leading to the appearance of the S −2 ε factor and chosen µ R = Q/2, though for brevity we have not explicitly written the argument of the running coupling above. We have introduced a factor of two to account for the other hemisphere containing the branching of the anti-quark.
Carrying out the integrals with the method discussed in the appendix, the result can be expressed in the following form: where one notes the presence of a double pole coming from the soft z p → 1 and collinear θ 2 23 → 0 limit and where singlepole contributions are separated into the contributions from soft (z p → 1) and collinear (θ 2 23 → 0) divergences alone given, respectively by the H soft and H coll. functions. We have The function H fin. (z) represents a finite contribution whose precise analytic form we have not extracted. Instead we study this finite contribution by direct numerical integration over the triple-collinear phase-space in 4 dimensions. The result for the integration of H fin. (z) gives a constant as z cut → 0. The result that we obtain using integration with Suave [51], setting z cut = 0 is, Next, we study the situation where the larger-angle gluon, i.e. emission 1, fails the z cut condition and is groomed away, corresponding to 1 > z > 1 − z cut . This leaves the mass to be set by the smaller-angle emission 2, ρ = z 2 z 3 θ 2 23 = z 2 z p (1 − z p )θ 2 23 . This emission must survive grooming i.e. 1 − z cut > z p > z cut and hence θ 2 23 1 for ρ z cut . The softness of emission 1, with energy proportional to 1 − z, implies that any terms regular in the limit z → 1 produce power suppressed corrections in z cut which we can neglect. Hence the only contribution comes from the singularity as z → 1, which produces an ε pole and associated finite corrections.
We start by considering the triple-collinear splitting function P (ab) q→g 1 g 2 q 3 and its integral over the phase-space, as for the previous case where emission 1 passes the z cut condition. However since now emission 1 fails the z cut condition and is groomed away, it is not constrained to be collinear and has a range of angular integration going from θ 13 ∼ θ 23 1 up to the boundary of the hemisphere at θ 13 = π/2. Near the lower limit of integration the triple-collinear approximation is valid, but to account correctly for soft emission at large angles we have to modify the angular dependence of the integral. This is precisely the same modification we made to account for soft large-angle emission for the calculation of C 1 (c.f. Eq. (7)). After neglecting pieces which contribute only an O(z cut ) term on integration, we find the result where δ ρ (2, 3) is the condition that emission 2 and the quark labeled 3 contribute to the hemisphere invariant mass ρ. In fact, one can directly reach the same equation by realizing that the emission probability of a hard-collinear gluon completely factorizes from that of a soft gluon, i.e. the gluons are emitted independently in this region of phase space. Therefore, the total emission probability is a product of an eikonal function and a LO splitting function. This factorized structure is manifest in Eq. (26). The integral over θ 2 23 is trivially performed using the delta function constraint which sets ρ = z 2 z p (1 − z p )θ 2 23 ≈ z p (1 − z p )θ 2 23 , where we have used the fact that z ∼ 1 corresponding to the softness of emission 1. 7 We have modified the angular dependence so that at small θ 13 we obtain the result arising from the triplecollinear splitting functions but for θ 13 ∼ 1 we have the correct angular dependence for a soft emission emitted off the qq dipole. We have also introduced the renormalised MS coupling and choose µ R = Q/2 as before. Evaluating the integrals we obtain: Finally we account for the correct action of the tagger when emissions 1 and 2 i.e. the two gluons are clustered first in the C/A algorithm and then the gluon pair is clustered to the quark. On applying the tagger one first encounters two branches, consisting of the quark and the massive gluon pair respectively. If the quark fails the z cut condition, one would then follow the branch consisting of the gluon pair. However such configurations with a soft quark are suppressed by powers of z cut and hence ignored. On the other hand configurations where the massive gluon branch fails the z cut condition would lead to a massless hemisphere. Hence we only need to study the situation where both branches pass the z cut condition and all three partons are retained. The opening angle between the branches is small, being set by ρ/z cut , which implies that all three partons are collinear and we can use purely triple-collinear kinematics. To correct our earlier results, we simply need to calculate the difference between the correct configuration described here and our simplified treatment included as part of F fail (ρ, ε) and F pass (ρ, ε).
The relevant angular region for the calculation is θ 2 12 < θ 2 23 , which corresponds to the C/A clustering of the two gluons, since we already have the ordering θ 2 13 > θ 2 23 . Our clustering correction takes the form where in the second line the first term in parentheses represents the correct treatment of the tagger while the second and third terms correspond to the removal of the gluon clustering region from F pass (ρ, ε) and F fail (ρ, ε) respectively. In the angular region relevant to their clustering, the smallest angle is that between the two gluons, so there is no collinear divergence in the gluon emission channel. We have potential soft divergences as each of z and z p tend to 1 (i.e. z 1 or z 2 vanish), but in both those limits the correct calculation cancels with the subtraction terms in the second line of Eq. (28), so that the result is purely finite and we can set ε = 0 in the integrals that follow.
In the correct treatment, i.e. the first step function on the second line of Eq. (28), the condition that the first declustering passes the z cut corresponds to 1 − z cut > z 3 > z cut , where z 3 = zz p . There are then two regions for the integration over energy fractions, according to the range of z: For the region z cut < z < 1 − z cut , this corresponds to the region in z where gluon 1 passes the z cut in our simplified calculations, so we need to correct the calculation of F pass (ρ, ε) and the third term in the second line of Eq. (28) vanishes.
In the region 1 > z > 1 − z cut corresponding to the grooming away of gluon 1 in the simplified calculation, we need to correct the treatment of F fail (ρ, ε) and here the second term in Eq. (28) vanishes.
In the former case one has a particularly simple situation since both in the simplified calculation and the correct treatment all three partons are retained and contribute to the jet mass via the δ ρ (1, 2, 3) condition. This results in an identical angular integration in both cases and the difference between the correct and simplified treatment is purely due to the different limits on z p . Denoting the result of the angular integrals by I(z, z p ), we can write the clustering correction as: where the subtracted term above represents the removal of the simplified calculation in the gluon clustering region and where I(z, z p ) has been written with a factor of two to account for both hemispheres and the angular integration extends only over the region relevant to the clustering of emissions 1 and 2. The simplified and correct calculations differ only in how soft the quark is allowed to be and in the limit z cut → 0, the result after integration over z and z p vanishes with z cut . This has been verified directly by numerical integration. Hence the correction term, given by Eq. (30), can be neglected in our approximation. The situation in the region 1 − z cut < z < 1, where we derive a correction to F fail (ρ, ε), is somewhat more subtle. In the simplified version of the calculation, emission 1 is groomed away and the mass is set by emissions 2 and 3 which leads to a different constraint given by δ ρ (2, 3) on the angular integration compared to the correct version where all three partons are retained with δ ρ (1, 2, 3), so one obtains a different resultĨ(z, z p ) given by: Due to the fact that we are restricted to the angular region where the two gluons would be clustered, all angles are constrained to be small and we can use the triple-collinear limit for obtainingĨ(z, z p ), ignoring the wide-angle modification required for the full calculation of F fail (ρ, ε).
One key point here is that in the limit z → 1, or equivalently z 1 → 0, where there is a potential soft divergence, the condition δ ρ (1, 2, 3) reduces to δ ρ (2, 3) so that the difference betweenĨ(z, z p ) and I(z, z p ) vanishes, leading to a finite result as already observed above. Moreover in the z cut → 0 limit we have that z → 1 over the full integration range so that one may simply replaceĨ(z, z p ) with I(z, z p ) up to finite corrections suppressed by z cut . Doing so leads to F clust.
where in writing the third line we have replacedĨ(z, z p ) with I(z, z p ) and have exploited the fact that the difference between the lower limits of z p integration in the correct and subtracted term, corresponding to the region of a soft quark, leads only to terms power suppressed in z cut . We can numerically evaluate the integrals for a given z cut value and on decreasing z cut , to reduce the impact of power suppressed terms, we find the result converges to a constant. For our smallest value z cut = 10 −5 , using the numerical method suave [51], we obtain F clust.
The smallest value of z cut was chosen so that the error on the numerical integration was larger than the difference between the central values for the lowest and second lowest z cut values.
The fact that the result for F clust.
tends to a constant as z cut → 0 is related to the behaviour of the integrand in the soft limit for both emissions z p → 1 and z → 1, and this allows us to also extract the result analytically. A series expansion of I(z, z p ) around z = 1 reveals a leading behaviour ∝ 1 (1−z)(1−z p ) , which derives from the soft limit of the matrix-element, and generates the full result in the limit z cut → 0. It is straightforward to perform the integrals analytically to obtain: where we have expressed the result in terms of the Polygamma function ψ (1) (x). Note that one can further write to obtain a compact result in terms of the Clausen function Cl 2 (x) [52]. This analytic result is consistent with the value obtained numerically for z cut → 10 −5 , keeping in mind that the latter includes power suppressed in z cut terms varying as z cut ln 2 z cut . We also note that ignoring a region of phase space which only contributes a power of z cut to F clust , the limits on the energy fraction integrals of Eq. (33) can be rewritten in terms of z 1 and z 2 : It is now apparent that clustering the two emissions together only leads to differences from our simplified treatment of the tagger where both emissions would separately have failed the z cut condition, but together lead to a cluster which passes the z cut condition. Our result for the F clust.
C 2 F turns out to be precisely the same as the result calculated previously, for the corresponding contribution to the non-cusp global soft anomalous dimension for the mMDT jet mass in SCET [15,19]. While our starting point using the triple-collinear splitting functions goes beyond just the soft limit, the observation made above that the relevant limit for F clust. C 2 F is the limit when emissions 2 and 3 are additionally soft, explains the agreement with the soft limit calculations of Refs. [15,19]. However it is worth stressing that our approach based on triplecollinear splitting functions remains valid beyond the soft limit and hence can also be used to compute the finite z cut corrections that we have neglected in the present article.

Virtual corrections and combined result
Here we combine the results for double-real emission with the one-real one-virtual corrections to generate a finite result. The one-real-one-virtual terms are provided in appendix Appendix D. We define the integral of the one-real-onevirtual term over z as V C 2 F , given by with V We find that after cancellation of all the singular contributions we are left with where the numerically quoted value 1.866 ± 0.002 represents the contribution to F pass (ρ, ε) arising from H fin. (z) (see Eqs. (21), (25)). It is then evident that the terms in Eq. (40) that depend on ln ρ and ln z cut are in precise agreement with those expected from the C 2 F term in the expansion of the leading-log resummed result, i.e. Eq. (16). In addition there is a constant contribution corresponding to an α 2 s ln ρ NLL term in Σ 2 (ρ). We shall analyse the constant contribution in more detail, after including another C 2 F term coming from the gluon decay terms computed in the next sub-section.

Gluon decay contributions
Here we consider the contributions that are associated to the decay of a collinear gluon, emitted off the initiating quark, into a qq pair and a gluon pair associated with C F T R n f and C F C A factors respectively. For the case of the gluon decay to qq with a quark initiated jet i.e. a q → qqq process, there is also an interference contribution from identical fermions in the final state with a colour factor C F (C F −C A )/2, which contributes to the overall results for the C 2 F and C F C A channels. We shall first discuss this piece and then turn to the C F T R n f and C F C A terms.
The identical fermion contribution is simple to compute since it is finite both for the angular and energy integrals. The calculation can therefore be easily carried out numerically in four dimensions. The relevant splitting function is given in Eq. (A.8) and we set ε → 0. Moreover due to the fact that the splitting function is regular in the energy fractions, the contribution from the region of integration where any parton has energy fraction z < z cut is suppressed with z cut . For this reason the clustering and grooming sequence does not matter as the result in the small z cut limit comes from a configuration when all three partons contribute to the jet mass ρ and have energy fraction z > z cut i.e. the ungroomed limit.
To obtain the leading term, which is a constant in the small z cut limit, we set z cut = 0 and numerically perform the integral using our general rescaling method discussed in the appendix. We then have where we performed the integral numerically with Suave and the result includes a factor of 2 to take account of both hemispheres as well an identical particle 1/2! phase-space symmetry factor. We believe that our result here coincides with an older calculation for the identical fermion contribution that enters initial state splittings, by Grazzini and de Florian, who obtained an analytical result which they subtract to construct the relevant non-singlet contribution (see Eq. (71) of Ref. [48]). We shall return to this result, its analytical form and its combination with the C 2 F and C F C A results, when summarising our results. Fig. 2 An illustration of the kinematic parametrisation in the C F T R n f gluon decay channel.

C F T R n f contribution
Next we discuss the C F T R n f term again related to gluon decay to a qq pair. The relevant diagram with our parameterisation of the kinematics is shown in Figure 2. One encounters a collinear singularity in the squared matrix element as θ 12 → 0, since the g → qq splitting is regular in the infrared, which leads to a 1/ε pole for the jet mass distribution. We also expect that an analysis of the gluon decay contributions should lead to the emergence of the correct argument of the running coupling i.e. k t in the soft limit, and the factor associated to the physical CMW scheme.
As we did for the C 2 F piece, it proves convenient to break the calculation into two pieces : a first piece that simplifies the action of the tagger and contains the divergent pole structure and a correction term which only has support in non-singular regions, leading to a finite result which can be computed numerically in four dimensions. We will also account for virtual corrections which cancel the pole in the real emission piece, leaving a finite result. The divergence occurs as θ 12 → 0 which is the region where the q andq are clustered first in the C/A algorithm, and then the qq pair is clustered to the parton 3. We shall therefore first carry out the calculation always taking the q andq to be clustered together which mistreats regions where emissions 1 or 2 could be first clustered with 3. These regions will be subject to our finite clustering corrections.
For the contribution clustering the q andq from the gluon decay, on declustering we produce two branches consisting of the quark (or anti-quark labeled) 3 and the clustered fermion pair or equivalently the massive parent gluon. If the quark labeled 3 is soft and fails the z cut then we obtain only terms power suppressed in z cut . If the massive gluon branch fails the clustering we obtain a massless hemisphere. Thus we have a situation where in order to obtain a finite result in the z cut → 0 limit, both branches pass the z cut condition so that all three partons are retained and contribute to the jet mass. The condition ρ z cut once again implies that all three partons are collinear and we can apply the triplecollinear splitting function and phase-space. In order to use our integration strategy based on rescaled angular variables, we consider two regions θ 2 12 < θ 2 13 and θ 2 12 > θ 2 13 , with the first mentioned region contributing the divergence. We denote the respective contributions to ρ dΣ 2 dρ by F θ 12 <θ 13 (ρ) and F θ 12 >θ 13 (ρ) where explicitly we have for the former, where a sum over flavours leading to a factor n f is left implicit on the RHS of the above equation. The notation Θ z cut (12|3) denotes that each of the two branches passes the z cut condition equivalent to the constraint on the parent gluon energy 1 − z cut > z > z cut , the result is written in terms of the renormalised MS coupling and a factor of two accounts for both hemispheres.
In terms of the rescaled angular variable y = θ 2 12 /θ 2 13 and the parent energy fraction z, we obtain a result of the form (where G(y, z, ε) is regular as y → 0) In writing the above we first isolated the singular contribution by taking only the leading term in the expansion of G(y, z, ε) around y = 0, and integrated over y to obtain the pole in ε in the first term above, while the second term is finite by construction since (G(y, z, 0) − G(0, z, 0)) vanishes as y → 0. The finite term can be computed in four dimensions, so we set ε = 0. However while the term involving (G(y, z, 0) − G(0, z, 0)) is finite, it has a leading behaviour in the z → 1 limit proportional to 1 1−z , which can be extracted through a series expansion about z = 1. The limit z → 1 corresponds to a soft parent gluon and gives rise to ln z cut terms in the result, which build up the constant K which relates the coupling in the MS scheme to that in the CMW scheme. After separating the 1/(1 − z) term which can be handled analytically, we integrate the remainder of the finite contribution numerically.
We then have where Next we need the contribution from the angular region θ 13 < θ 12 which does not contain any poles. However, as discussed above for θ 12 < θ 13 , there is again a soft enhancement as z → 1 giving rise to a ln z cut term related to the CMW constant K. Explicitly we have: which has no collinear pole i.e.
where, as before, the final term is evaluated numerically with z cut = 0 to remove power corrections. A final ingredient for generating the result is the T R n f piece of the one-real-one-virtual contribution. This is reported in Appendix D (see Eq. (D.27)) and a finite result is obtained on combining the virtual contribution with F θ 13 >θ 12 (ρ). Defining and we can combine the results to obtain where It is worth making some remarks on the form of G(z, ρ). Firstly we note a piece corresponding to the n f term in −2 × p qq b 0 ln (ρ(1 − z)), where b 0 = 11 6 C A − 2 3 T R n f is the first perturbative coefficient of the QCD β function. This term produces an LL contribution from the dependence on ln ρ. In the soft limit, i.e. z → 1 and p qq (z) → 2 1−z it can be absorbed into the leading-order result by correcting the scale of the coupling α s , consistent with the LL formula (3). The term − 20 9 p qq (z) corresponds to the appearance of the n f piece of the CMW constant K, as also anticipated in the LL result. The remaining terms produce a constant in the small z cut limit, on integration over z, which can be combined with the constants we obtained numerically. However before doing so we shall evaluate the correction term due to the proper use of the C/A clustering sequence within the mMDT.
Turning to the clustering correction we first note that in the angular region where θ 12 is smallest our treatment of the tagger, working as if partons 1 and 2 are always clustered, needs no correction. The regions where a correction is needed are when θ 13 is smallest and when θ 23 is smallest which is identical due to the symmetry under 1 ↔ 2. These regions contain no divergences and hence the calculation of the correction term, i.e. the difference between the correct and simplified calculations, can be performed numerically in four dimensions. We shall also work in the limit z cut → 0 to eliminate power-suppressed terms in z cut , explicitly take the case that θ 13 is the smallest angle and double the result to account for θ 23 being smallest. We can then write The above equation reflects that when θ 13 is the smallest angle, there are two configurations that yield a massive hemisphere: when parton 2 passes the z cut condition and all three partons are retained (the first term on the second line), and when parton 2 fails the z cut condition but partons 1 and 3 pass (the second term on the second line). The clustering correction may then be expressed as where the subtraction term corresponds to removal of the simplified contribution where 1 and 2 were taken to be clustered first and our usual notation applies where Θ z cut (a|b) denotes the two branches a and b that pass the z cut condition in each case. 8 The integrals I n f (z, z p ) andĨ n f (z, z p ) arise from the angular integration for the case when all three partons contribute to the jet mass and when only two partons contribute respectively. 9 Explicitly we have 8 This implies also the condition that partons not included in these branches fail the z cut e.g. for parton 2 in Θ zcut (1|3). 9 In the case when a parton is soft enough to be groomed away it is not constrained by the jet mass and can in principle be at a large angle. For correlated emission such configurations, where one of the gluon decay offspring is at a large angle to the other, are dynamically suppressed and only contribute at the level of power corrections in ρ. This implies that the integrals converge within the triple-collinear region which does not need to be modified, which we have also verified numerically.
and I n f (z, z p ) = 4 (8πα s ) 2 s 2 123 P q→q 1 q 2 q 3 dΦ 3 δ ρ (1, 3) where a factor of 4 accounts for both hemispheres and the case when θ 23 is the smallest angle. We can evaluate the integrals numerically by choosing a small z cut to suppress power corrections and find that the result tends to a z cut independent constant on decreasing z cut . On evaluation of the integrals numerically with Suave we obtain, for z cut = 10 −5 as for the C 2 F clustering piece, the result Further insight into the nature of the clustering correction reported above can be obtained via similar considerations to those for the C 2 F clustering correction. Firstly one notes that in the limit parton 2 goes soft and fails the z cut condition the angular integralĨ n f (z, z p ) may be replaced by I n f (z, z p ) up to terms that vanish with z cut . This lets us combine the constraints on the z integrals and, again, with neglect of power corrections in z cut one obtains the conditions on z 1 and z 2 : These conditions are the same as for the C 2 F case, however here the clustering correction enters with a negative sign while a positive correction was noted for the C 2 F term. The reason for this is that our simplified treatment in the C 2 F channel amounted to discarding two emissions that individually failed the z cut , thereby excluding the contribution where they pass the z cut when correctly treated as a cluster since z 1 + z 2 > z cut . Here on the other hand, our simplified picture includes configurations where incorrectly treating emissions as a cluster they pass the z cut condition, while in the correct treatment of clustering where the gluon decay products are not clustered the emissions each fail the z cut leading to a massless hemisphere.
Furthermore, we note once again that our clustering corrections originate in the soft region, albeit still also within the triple-collinear regime and that our calculations can be extended to include finite z cut terms. The numerical value we obtain for the clustering piece, F clust.
C F T R n f , is once again in agreement within errors (and potential z cut ln 2 z cut terms) to that previously obtained for the n f part of the clustering term for the non-cusp global soft anomalous dimension for mMDT in the SCET framework [15,19].
Finally we quote our overall result for the C F T R n f channel, combining the different contributions i.e. performing the integral over z in Eq. (52) and adding in the clustering correction: where the −1.754 is the clustering correction computed above, and the other terms are the full result for our simplified treatment of the tagger. We note that the ln ρ and ln z cut dependent terms in Eq. (60) are in exact agreement with expectations from the expansion of the leading logarithmic result (16). We shall shed further light on the constant term 4π 2 9 + 25 3 − 7.721 ± 0.007 following an analytic calculation in the next subsection.

Analytic calculation using web variables
Our aim is to compute the simplified piece of the differential distribution, i.e. F tot. , which neglects the proper C/A clustering in specific angular regions, using the web variables given in Appendix C. In particular, this allows for a complete analytic extraction of the constant in eq. (52). In addition, the analytic computation offers insight into the structure of the answer for the C F C A channel to be presented in section 4.2.4, which shall become evident when we discuss the final results.
First, we recall consequence of the angular ordering property of soft radiation.
which is written as before in terms of the renormalised MS coupling with µ R = Q/2, and a factor of 2 to account for both hemispheres. Further analysis using our general integration method, outlined in the appendix, gives where As in section 4.2.2 the finite term H finite (z) is enhanced in the limit of a soft parent, z → 1, and produces ln z cut terms, which we wish to separate since they relate to the CMW scheme. As before we use a series expansion about z = 1, to make the decomposition where f (z) is finite as z → 1. The constant, c is evaluated numerically (on integrating over the angular variables and z p ), as is the integral over f (z) and, again using NIntegrate with z cut = 0, we obtain: In the region θ 13 < θ 12 there are no poles in ε and we can perform the calculation setting ε → 0, so we have Again separating the integrand into pieces which diverge as z → 1 and those which do not, we obtain, after numerical integration with z cut = 0, Finally we account for the clustering corrections. This is done as for the n f piece in Eq. (54) and the result is finite as soft divergences cancel in the combination of the correct and simplified treatments. The result can be computed numerically in four dimensions. Using Suave with z cut = 10 −5 we obtain: The same comments apply to the origin of the clustering correction here as for the n f piece, namely it originates from incorrectly allowing, in the simplified result, the gluon pair to pass the clustering due to the fact that the parent passes the z cut . The correct tagging procedure would be applied to the individual gluons instead, which fail the clustering leading to a massless jet and a nil contribution. The result again agrees with previous calculations of the clustering piece in the mMDT SCET non-cusp soft anomalous dimension [15,19] to within errors and potential z cut ln 2 z cut corrections.
Finally we combine all pieces and include the one-realone-virtual correction V C F C A (ρ, ε) = 1−z cut z cut dzV C F C A 1,1 (ρ, z, ε) defined in the appendix to obtain the result ρ dΣ 2 dρ q→qgg,nab.
= C F C A α s 2π 2 11 6 (3 + 4 ln z cut ) ln ρ + 11 ln 2 z cut 3 where the labelling q → qgg, nab. indicates the non-abelian contribution to the q → qgg process and where we have separately written the numerically computed clustering contribution and the numerically computed part of the simplified calculation. The result above is in agreement with our expectations from the expansion of the leading-logarithmic resummed result Eq. (16) for the terms involving ln ρ and ln z cut since the numerical value 8.3091 ± 0.0006 is in good agreement with 134 9 − 2π 2 3 , signalling again that the ln z cut term is associated with the CMW scheme, while the ln ρ and ln 2 z cut terms are associated to the argument of the running coupling i.e. k t in the soft limit. This leaves us to comment on the constant term, other than the clustering correction, which we shall do in the next section, where we shall consider the full C F C A result including that from the C F (C F −C A /2) term.

Structure of NLL results
In this section we discuss the structure of our results for each of the C 2 F , C F C A and C F T R n f channels. As has already been noted for every channel, the result at order α 2 s reproduces the terms expected from the expansion of the LL formula (16) in addition to producing genuine NLL corrections. 12 Therefore we may focus only on the additional terms not produced as part of the LL expansion and hence we write: where F clust.
is the clustering contribution Eq. (35). The corresponding result for the C F T R n f term can be obtained partly in numerical form from Eq. (60) or fully analytically from Eq. (69) and after removal of the LL contribution we get where we used above the analytical result of Eq. (69) and the clustering contribution is given in Eq. (58). The result for the C F C A channel is obtained by combining Eq. (78) with the identical particle contribution in 12 Recall that the expansion of the LL formula also contains formally NLL terms but which can be embedded within the LL strongly ordered dynamics.
Eq. (41) and removing the LL contributions so that we obtain where F clust. C F C A is reported in (77). Leaving aside the clustering corrections for the moment, whose soft (and collinear) origin we have already discussed, we focus on the structure of the rest of the result. It is wellknown that the intensity of collinear radiation from a quark at second order in α s is related to a coefficient in the quark form factor generally referred to as B (2) [54][55][56][57]. While there is not a unique definition of B (2) since it depends on the details of how the full resummation formula is organised, i.e. the resummation scheme, it is always related to the endpoint δ (1−z) contribution to the NLO DGLAP splitting functions via the form (for a quark initiated jet) [58][59][60], where b 0 = 1 6 11C A − 4T R n f and where γ (2) q , the DGLAP endpoint contributions for a quark jet, are [61,62] γ (2) We note that our analytic result for the coefficient of C F T R n f α 2 s /(2π) 2 in Eq. (81) is precisely consistent with the form in Eq. (83) with We conclude that without the clustering corrections our NLL results for the mMDT as defined above, are given by the general collinear form Eq. (83) with the value of X specified in (85). Thus excluding the clustering corrections our NLL result for the mMDT has a simple correspondence to the collinear order α 2 s contribution to the quark form factor. Moreover the results for the C F C A and C F T R n f channels, without the clustering correction, agree exactly with the order α 2 s NNLL term in the expansion of the ungroomed heavy jet mass [43,45,63], a consequence of grooming affecting only soft emissions in the small z cut limit. Finally we note that our overall results in every channel are in agreement with those from previous SCET calculations. An explicit expansion of the SCET results to order α 2 s , for the jet mass distribution, has recently been provided in Ref. [24]. 13

Conclusions
In this article we have revisited the NLL structure of the jet mass distribution for mMDT groomed jets from the viewpoint of its direct connection to the QCD matrix elements in the triple-collinear limit. Previous NLL results have entirely been within the framework of Soft-Collinear Effective Theory (SCET) and hence our work represents an approach which provides strong independent confirmation of the main results involved in the NLL resummation [15,19,30,31].
Our results establish a connection between the NLL groomed jet mass result and the standard ingredients used in QCD resummation. In particular we recover the expected scale of the running coupling in the soft limit, i.e. the k t of a soft emission, and the constant K related to the emergence of the CMW coupling. We further obtain a link between the NLL result and the general form of the B (2) coefficient that controls the intensity of collinear radiation from a quark at order α 2 s and hence enters the quark form factor. Our result also involves a clustering correction in all channels, which stems from our simplifying the action of the mMDT to derive the B (2) -like terms. The clustering corrections come from a region of phase space where we have two soft emissions that if examined individually fail the z cut condition but if examined as a cluster pass the z cut condition. We believe that these results should allow for a resummation of the mMDT NLL corrections within a QCD resummation framework. In fact the B (2) -like pieces are already incorporated in an approach such as ARES [43,45] as they also enter into the NNLL structure of the ungroomed heavy jet mass. 14 It therefore remains to consistently include the 13 Note that we have removed ln 2 terms in the C F C A and C F n f channels present in the results of [24] via the choice of Q 2 /4 in the scale of α s for our leading-order result. 14 In ARES [43,45], for the heavy jet mass the factor X in the definition used for B (2) is equal to zero and the b 0 X term is associated to the functions C (1) hc and δ F rec . clustering corrections within the standard QCD resummation formalism.
For future extensions of our work, one development that is possible to make concerns the inclusion of finite z cut corrections beyond the LL level [12,16,17]. These can be derived through our triple-collinear calculations retaining terms that we have omitted in the present article by taking the small z cut limit of various formulae. Since these additional corrections will be purely finite, they can be computed numerically in four dimensions and incorporated into the resummation framework described in Ref. [12]. While we expect the resulting corrections to be numerically small, they should become relevant to examine in the context of recent developments pushing the mMDT jet mass resummation to NNLL level [23]. It would also be of interest to use our approach to study β = 0 values for SoftDrop and hence to develop an NNLL QCD resummation approach for those observables. Our triple-collinear calculations should also give the insight needed to address other similar collinear problems at the NLL level such as that involving the small jet radius limit of QCD jets, for which an LL resummation formalism was constructed in Ref. [29] but a general resummation approach at NLL is still missing. q 1 q 2q 3 ) or q → q 1q2 q 3 (q →q 1 q 2q3 ), s 123 is the squared invariant mass of the three collinear parton system, and P k 1→3 are the process independent spin averaged triple-collinear splitting functions.
Hereinafter, we report the relevant expressions for the triple-collinear splitting functions [35]. Note that due to charge conjugation invariance, the splitting functions for the antiquark initiated channels can be obtained from the corresponding functions for the quark initiated ones,. i.e.Pq 1 q 2q 3 = Pq 1 q 2 q 3 andPq 1 q 2q3 =Pq 1 q 2 q 3 .
The spin-averaged splitting function for the q → g 1 g 2 q 3 process can be written in terms of the different colour factors: can be expressed in terms of our chosen variables as a condition on x: (B.9) It proves to be convenient to map the integral over x to one with simple limits i.e. 0 and 1 by introducing the change of variables x = u(r 2 − r 1 ) + r 1 . Our integration variables are then u and y for the angular integration, both lying in a range 0 to 1, and the energy fractions z 1 and z 2 (recall that z 3 = 1 − z 1 − z 2 ) or equivalently z and z p . We then have to consider the extraction of ε poles, to separate the integral into divergent and finite terms. Our strategy is to isolate the divergences and exploit the simplification of the integrand in divergent regions, to obtain the divergence structure analytically. This also generates finite terms that do not vanish as ε → 0 , which are obtained via an ε expansion of the factors multiplying the poles. Additionally we also obtain a finite integral leftover from the removal of singular terms, which on the other hand is not a compact expression. However, being finite, it can always be integrated numerically.
Since we study the differential distribution rather than its integral, we work at fixed jet-mass which regulates both soft and collinear divergence. In general that leaves us with at most a 1/ε 2 singularity from an emission that does not set the jet mass. For the calculation of F pass in particular, the larger-angle emission passes z cut and cannot produce any divergence, while the smaller angle emission produces divergences from the soft z p → 1 and collinear y → 0 limits. Setting 1 − z p = v we encounter a general integral of the standard form I(ε) = where G(v, y, ε) is finite as v → 0 as well as y → 0 and integration over the other variables is left implicit so as to focus on the divergences. We can re-express this result in the following form where by construction the first term on the LHS of the above has only a soft pole i.e. as v → 0, the second term has only a collinear pole from y → 0, while the final term has a double pole arising from v → 0 and y → 0. We define for convenience f (v, ε) = (G(v, 0, ε) − G(0, 0, ε)) /v 1+2ε where f (v, ε) is finite as v → 0 and also define h(v, y, ε) = (G(v, y, ε) − G(v, 0, ε)) /y 1+ε , which is finite as y → 0. Then one obtains the form h(0, y, ε)dy dy . (B.12) The final integral on the RHS above is purely finite by construction and can be evaluated in the limit ε → 0 i.e in 4 dimensions. The above result shows explicitly the pole structure that emerges from the integral I(ε). The integrals multiplying the 1/ε poles need only to be evaluated up to order ε terms i.e. one can expand the integrand in ε and retain only terms up to order ε. This strategy gives us all divergent and finite contributions in the limit ε → 0. 1 (C.23) Here, the solid angle dΩ 2−2ε is that of q ⊥ in the transverse plane aligning k ⊥ along one axis. A nice feature of eq. (C.23) is the simplicity of the double-soft limit, z → 1. In fact, if we set z = 1 in the measure one recovers the double-soft phase space, the 4-dimensional limit of which is reported in [64].
Appendix D: One-loop corrections to 1 → 2 collinear splittings In addition to the case of two real emissions, for the jet mass distribution at order α 2 s we also have to consider a real emission that sets the mass ρ alongside a one-loop virtual correction which is divergent and where the divergences cancel against those in the double-real emission case, to leave behind finite terms. The relevant real-virtual contribution to ρdΣ 2 /dρ, may be collectively written in the form where we have separated out the various contributions according to the colour factor i.e. C 2 F , C F C A and C F T R n f terms and in our notation V 1,1 is the one-real, one-virtual correction to qq production, in the approximation of a real emission, which is collinear to the q orq, passes grooming and sets a (normalised) jet mass ρ.
(D. 25) Writing our result in terms of the renormalised MS coupling, accounting for both hemispheres with a factor of two, we have