Eur. Phys. J. C (2018) 78:119 Page 5 of 24 119
presence of photon splittings as their recoil partner assign-
ment is arbitrary and may involve information on initial or
final state particles. Different possible choices are discussed
in Sect. 3.2, some of which may break this crossing invari-
ance.
The divergences of the I operator are encoded in the func-
tions V
ik
and
i
. While the former contains all soft-(quasi-
)collinear divergences the latter contains the pure (quasi-)
collinear ones. They do not only differentiate whether the
emitter is a photon or not, but also between different spins
of the emitter. Their precise form as well as the the flavour
constants γ
i
and K
i
are given in Appendix B. A
I
ik
encodes the
dependence on the phase space restriction of the individual
dipoles {α
dip
}. Finite terms originating in dipoles involving
initial state legs, however, can be pushed into the K opera-
tor. Thus, A
I
ik
by convention only depends on α
FF
. Its precise
form is given in Appendix C.
The K and P operators. The K and P operators collect all
pieces of the integrated dipole terms that are not collected in
the I operator and combines them with the collinear coun-
terterms C to give a finite result as → 0. By construction
they contain only remainders of splittings where either the
emitter or the spectator is in the initial state. Thus, they are
comprised of terms arising due to the change of the flavour
or the partonic momentum fraction x of an initial state due
to a splitting.
The K-operator is given by [41]
K
aa
(x ;{α
dip
})
=
α
2π
K
aa
(x ) − K
FS
(x ) −
i
Q
2
ia
K
i,aa
(x )
−
k
Q
2
a
k
K
t
aa
,k
(x ) − Q
2
a
b
˜
K
aa
(x ) + A
K
aa
({α
dip
})
.
(2.9)
It depends on the partonic x, and the flavour change from the
Born process initial state flavour a to a
of the convolution
Eq. (2.6). Therein, the
K collect universal terms present in
all splitting involving an initial state as either emitter or spec-
tator. Then, while K contains solely remaining terms from
final state splittings in the presence of initial state specta-
tors, the K
t
are their counterparts for initial state splittings
in the presence of a final state spectator, i and k running over
all final state partons in each case.
˜
K contains solely related
correlations between both initial states, arising from dipoles
where both the emitter and the spectator are in the initial
state. The A
K
terms collect all finite terms arising when any
of α
FI
, α
IF
or α
II
is different from unity, thus restricting the
phase space of the respective dipoles. Again, the Q
2
ik
are the
charge correlators of Eq. (2.5). Finally, K
FS
contains the fac-
torisation scheme dependence. Currently, both only the
MS
schemes is supported, setting these terms identically zero.
The P-operator now collects the remaining initial state
collinear singularity from all dipoles involving initial states
either as emitters or as spectators and cancels them against
the collinear counterterm. Through this counterterm a depen-
dence on the factorisation scale enters. The P operator is
given by [40,41]
P
aa
(x,μ
2
F
)
=
α
2π
P
aa
(x)
k
Q
2
a
k
log
μ
2
F
xs
ak
+ Q
2
a
b
log
μ
2
F
xs
ab
.
(2.10)
Only initial state splittings are present, either in the presence
of a spectator in the final state, which is encoded in the sum of
k, or with the opposite initial state b acting as the spectator.
It otherwise only depends on the Alterelli–Parisi splitting
function detailed in Appendix B.
3 Implementation
The implementation of the QED generalisation of the Catani–
Seymour dipole subtraction scheme in S
HERPA’s matrix ele-
ment generator A
MEGIC proceeds along the lines of [49]. As
in general real and virtual corrections of O(α
n
s
α
m
) contain
divergences of both QCD and QED origin, both cases are
included in this section. In the following, the general struc-
ture of the implementation is reviewed.
3.1 Identification of dipoles
The starting point to construct the involved subtraction terms
in the Catani–Seymour subtraction formalism is a given
flavour configuration in the Born or the real emission phase
space and the perturbative order O(α
n
s
α
m
) in accordance with
the respective virtual or real correction to be computed. For
all parts, on-the-fly variations of both the factorisation scale
μ
F
and the renormalisation scale μ
R
are available through
an extension of the algorithm detailed in [56].
Differential subtraction terms. A given real emission con-
figuration {ab}→{1, .., m + 1} at order O(α
n
s
α
m
) can in
general exhibit both QCD and QED divergences simultane-
ously. The following therefore describes the identification of
both types of dipoles. Thus, all triplets {i, j, k} that can be
built from the external particles of the process are tested for
the presence of an infrared divergence, QCD or QED, by
checking for the existence of a dipole subtraction term. In
these triplets i and k may be in the initial or final state while
j may be in the final state only. Likewise, i = j, i = k, j = k
and triplets that only differ in a permutation of i and j are
considered identical. Then, the following steps are executed.
123