text:
# 1 Introduction

## 1 Introduction

The standard model(SM) of electromagnetic and strong interactions has been well established through precision measurements of a large variety of observables in high-energy experiments. However, there is still one missing part of the SM to be confirmed by experiments, i.e. the Higgs boson. In order to search for this missing part and also search for direct signals beyond the SM, the LHC experiments are under construction at CERN. Though LHC has employed a proton-proton colliding machine due to their high discovery potential for new (heavy) particles, a capability for precision measurements must also be seriously considered. For extracting physically meaningful information from a large amount of experimental data, suffered from huge QCD backgrounds, a deep understanding of the behavior of QCD backgrounds is essentially important.

The prediction power of a lowest calculation of QCD processes is very limited due to a large coupling constant of the strong interaction and an ambiguity of the renormalization energy scale. Predictions based on the NLO (or more) matrix elements combined with some all-order summation are desirable. Much effort to calculate NNLO corrections is being made for LHC experiments. Moreover, in order to suit for precision measurements, a simple correction factor (so-called $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{K}$}`$-factor) for the lowest calculation is clearly not sufficient. One must construct event generators including NLO matrix elements with higher order re-summation even for processes with multi-particle final states. When one tries to do this task, one may encounter the following difficulties:

* the number of diagram contributing to processes with many particles in a final state is very large,
* many parton-level processes must be combined to construct a proton-(anti-)proton collision,
* numerical instability in the treatment of a collinear singularity may appear,
* a careful treatment of matching between matrix elements and the re-summation part must be required to avoid double counting, and
* negative weight of events in some part of phase space is not avoidable.

In order to overcome these difficulties, lots of works have already been done. In this report, we propose a new method to construct event-generators with next-to-leading order QCD matrix-elements and leading-logarithmic parton showers. Our solution for the above problems will be given in this report.

## 2 Matrix elements

A calculation of the matrix elements of hard-scattering processes is not a simple task when multiple partons are produced. At the tree level, the GRACE system, an automatic system for generating Feynman diagram and a FORTRAN source-code to evaluate the amplitude, is established. The system can in principle treat any number of external particles, and has been used for up to six fermions in the final state within a practical CPU time. In the GRACE system, matrix elements are calculated numerically using a CHANEL library based on a helicity amplitude. CHANEL contains routines to evaluate such things as: wave-functions/spinors at external states, interaction vertices, and particle propagator.

Based on this method, we have already published an event generator of four $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}`$-quark generation processes at proton-(anti-)proton collisions at tree-level named GR@PPA\_4b, which include all possible tree-level diagram to create four $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}`$-quarks in final state with an interface to PYTHIA.

For the loop diagram, effective vertices of two- and three-point functions have been prepared and implemented in the system. At first amplitudes of the loop diagram are evaluated in $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}`$-dimensional space-time. The dimensionless coupling constant is introduced as $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{/}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}`$. After the loop-momentum integration, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{U}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{V}$}`$ divergences are eliminated by subtracting $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{/}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{U}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{V}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\gamma }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{E}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{ln}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}`$ where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{U}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{V}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{/}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\gamma }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{E}$}`$ is the Euler’s constant. This is so-called $`\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{M}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{S}$}}`$ renormalization scheme. Then we are left with the infra-red divergence which is also regularized using the dimensional regularization by setting $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}`$, where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{>}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$.

Expressions of effective vertices after the $`\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{M}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{S}$}}`$ renormalization are summarized in Appendix A. For more than three-point diagram, an extended CHANEL is prepared. It can evaluate a fermion current including any number of gamma matrices. An amplitude of loop diagram will be calculated based on the fermion current connected by boson propagators. When a fermion loop is included in the diagram, numerical calculation will be performed after taking its trace. A bosonic loop can also be treated by the CHANEL.

## 3 Loop integral

As two- and three-point functions are embedded in the system as effective vertices, loop integrals which must be done in each process start from four-point functions. The loop-momentum integration of the four-point diagram can be expressed as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{U}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{V}$}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}}\frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{A}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{A}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{A}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{A}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (1)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{A}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{A}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{A}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{A}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}`$ is a loop momentum and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}`$’s are external momentum taken to be incoming. After a momentum integration using an usual procedure of Feynman-parameter method, the integral can be represented as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Gamma }}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{J}$}}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{m}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{m}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{m}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{m}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (2)
where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{m}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$, and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$. In general we need tensor integrals

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{J}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{;}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ (3)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
here we set all external particles on-shell (massless). A method to perform this integration is well established. The scaler integrations (when $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$) without an infrared divergence are given in  and those with IR divergence are given in . Tensor integrals can be obtained from a linear combination of scaler integrals . Even though the method is well established, we decided to develop our own method for the tensor integration of an IR divergent case, in order to avoid a possible numerical instability due to some cancellation among IR divergent scaler integrals.

The tensor integration can be done analytically and be represented by a finite number of terms with Beta and Hypergeometric functions as;

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{J}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{;}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{B}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}}`$ (4)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\times }$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{[}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{B}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\times }$}`$ $`{}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{}^{}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \frac{\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}}}{\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}}{\overset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{B}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\times }$}`$ $`{}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{}^{}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \frac{\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}}}{\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{]}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{.}$}`$
When the numerator is unity ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{y}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{z}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$), the scalar integration can be obtained as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{J}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{;}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{B}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}}`$ (5)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\times }$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left[}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}{}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{}^{}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \frac{\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}}}{\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}{}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{}^{}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \frac{\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}}}{\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right]}$}`$
A explicit structure of the infra-red singularity can be obtained with Laurent expansion of above formula with respect to $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}`$. A full set of formulae and a computer program for the numerical calculation can be found in

For more than four-point functions, a reduction formula is used in the system.

## 4 Real emission

### 4.1 Soft/Collinear treatment

At first, the conventional method of phase-space slicing is used to treat soft and collinear singularities. It is explained for the case of an initial-state radiation in this report. The final-state radiation can be treated in a similar way. Let’s consider two colored partons, whose momenta are $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}`$ and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}`$, scattered into the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}`$-body final state of colorless particles as a Born process. For the radiative correction of this process, we must treat processes with one additional colored parton, whose momentum is $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}`$, emitted in addition to the Born process. The matrix elements of the real emission processes must be integrated in ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}`$+1)-body phase-space, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}`$, in a $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}`$-dimensional space-time. The space-time dimension is set to be $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}`$ in this report. A collinear region in the ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}`$)-body phase-space of the final particles in the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}`$-dimensional space-time is defined as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{o}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{<}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (6)
as shown in figure 1, where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ is some cut-off value. It must be sufficiently small so as not to be observed experimentally. Final results must be independent of this value. In the collinear region, matrix elements can be approximated as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\overset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{16}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝙿}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (7)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝙿}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ is a splitting function for a given parton splitting at a leading-logarithmic order, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ is the transverse momentum of the radiated parton, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}`$ is a color factor of a given branching and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}`$ the energy scale of the splitting. The CM-energy of the Born process is $`\widehat{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$, where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$. Here, we can set $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}`$ to be four in the matrix element of the Born process, since it has no IR singularity. In the same approximation, the phase space is expressed as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\overset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}`$ (8)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\times }$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{16}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Gamma }}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{.}$}`$
The total cross section of real-radiation processes in the collinear region can be obtained as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{o}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{o}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (9)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left[}$}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right]}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (10)
where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}`$ means a visible region of the phase space, such as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{o}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{.}$}`$ (11)
Then, the collinear cross section can be obtained after the integration with respect to $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{o}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (12)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (13)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{B}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Gamma }}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝙿}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (14)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (15)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝙿}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{ln}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Theta }}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (16)
where

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\overset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (17)
is the Born cross-section at a CM energy of $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}`$ in four-dimensional space-time, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}`$ is a factorization energy scale, and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Theta }}$}`$ is a step function. For example in case of a splitting for a quark to quark+gluon, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{o}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}`$ after the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}`$ expansion becomes

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{o}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left[}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right]}$}`$ (18)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\varphi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left[}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{ln}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{ln}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right]}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{ln}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Theta }}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
where

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\varphi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝙿}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (19)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{ln}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{/}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{.}$}`$ (20)
The IR-divergent terms in the collinear cross section can be canceled out after summing up terms from virtual corrections. After combining with virtual corrections, there still exists the collinear-divergent term, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\varphi }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$. This term is thrown away by hand, because it is counted in the PDF or PS on the initial partons. Since the $`\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{M}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{S}$}}`$ renormalization scheme is used in this calculation, the same scheme must be used in the PDF or PS.

### 4.2 Visible jet cross-section

The last part of the NLO cross-section calculation is a real emission of the additional parton into the visible region with exact matrix elements. Those matrix elements of ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}`$+1)-parton production at the final state in four-dimensional space-time can be automatically generated using the GRACE system. Cross sections are obtained by integrating those matrix elements under the four-dimensional phase-space as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{.}$}`$ (21)
The space-time dimension is set to be four ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$) hereafter, since there is no IR-divergence in $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}`$. When the threshold energy ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$) is set to be sufficiently small, this cross section can be larger than the Born cross section due to the large magnitude of the coupling constant ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}`$). A parton-level calculation has no problem, except for the large higher-order correction. However, if one tries to calculate the cross sections of a proton (anti-)proton collision, one has to combine the matrix elements with the PDF or PS to construct a proton from partons. The PDF and PS include leading-log(LL) terms of the initial-state parton emission. If the matrix element is combined with the PDF or PS very naively, a double-counting of these LL-terms must be unavoidable. Our proposal to solve this problem is to subtract the LL-terms from the exact matrix elements as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left[}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right]}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (22)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝙿}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{16}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{.}$}`$ (23)
The second term of the integrand is the LL-approximation of the real-emission matrix-elements under the collinear condition. There is nothing but ’$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}`$’ in the last subsection. Then the real-emission cross section ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}`$) can be expressed as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{o}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (24)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (25)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{.}$}`$ (26)
One can expect a numerically stable calculation using Eq.(26) rather than Eq.(24), because cancellation between $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}`$ and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}`$ occurred before the phase-space integration. When the threshold value of $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ is set to be sufficiently small, the integrand in $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}}`$ is very close to zero, because the LL-approximation is very precise around the collinear region. Then, the result of integration in $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}}`$ is independent of the threshold value of $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$.

This subtraction may distort the experimental observables of additional jet distributions. However, a subtracted LL-part will be recovered after adding the PS applied in the Born process. The exact matrix element of real-radiation processes gives only non-logarithmic terms to the visible distributions in the LL-subtraction method. The relation between the PS and the LL-subtraction method is as follows. The LL-approximation cross section in $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}}`$ is

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ (27)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒟}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (28)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒟}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{/}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{m}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝙿}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (29)
where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{m}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ is the maximum value of the transverse momentum. $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒟}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ is just the first term of the all-order re-summation of the leading logarithmic terms in the PS. Actually, some special case of the Sudakov form factor can be obtained from $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒟}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$as:

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Pi }}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{m}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{/}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒟}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{.}$}`$ (30)
However, there is one essential difference between $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒟}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ and the Sudakov form factor in the PS: the upper bound of $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ integration. In the PS, this integration is performed up to the energy scale determined by the Born process instead of their maximum value in kinetically allowed region. Then, if $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}}`$ is subtracted in the full phase space of the ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}`$+1)-particle final state, some part of the LL cross section cannot be recovered from the Born cross-section with the PS. Then, an appropriate restriction on the phase-space integration must be applied, such as

$`\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Omega }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{v}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Phi }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left[}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left|}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right|}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Theta }}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{B}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right]}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (31)
where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{B}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ is the energy scale of the Born process, which must be the same as the factorization energy-scale of the PS.

## 5 Event generation with a parton shower

### 5.1 Conventional method

The calculations of cross sections for proton-(anti) proton collisions are usually performed as follows:

1. matrix elements of a given process at a parton level are prepared,
2. a probability to observe a parton in a proton at given energy scale and momentum fraction is obtained from some parton distribution function(PDF),
3. a numerical integration is performed in phase-space with some initial energy distribution in the PDF,
4. for initial-state partons, a (backward) parton shower is applied from the given (high) energy scale to a low energy-scale of partons,
5. for final-state partons, a parton shower is also applied from a given energy scale to close to the hadron energy-scale,
6. some non-perturbative effects according to string fragmentation or color clustering are taken into account, and
7. physical hadrons are formed based on some hadronization model.

In this method, the parton shower(PS) is used in a supplementary way, just for generating multi-partons with a finite transverse momentum.

### 5.2 An $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$-deterministic PS

The parton shower is a method to solve a DGLAP evolution equation using a Monte Carlo method. The PDF’s are also obtained by solving the DGLAP equation with experimental data-fitting. When an initial distribution of partons is given at some energy scale, the PS can reproduce a consistent result as the PDF. For singlet partons, we employ an evolution scheme based on momentum distributions rather than those on particle-number distributions, which is used for non-singlet processes. The momentum-distribution evolution-scheme has been proposed by Tanaka and Munehisa and shows numerically stable results with high efficiency.

The PDF can give the weight of partons when the momentum fraction ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$) and the energy scale are given by the users. On the other hand, in the forward evolution scheme of the PS, the momentum fraction of a parton ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$) is determined only after evolution takes place. This method is very inefficient, for instance, for a narrow-resonance production. In order to cure this inefficiency, we have developed an $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$-deterministic PS. In a Monte-Carlo procedure in the PS, the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ evolution is controlled by the Sudakov form-factor, which gives a non-branching probability of partons, and the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$ determination is done independently according to a splitting function. It is not necessary to determine the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$ value at each branching. After preforming all branching procedures using the Sudakov form-factor, the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$ value at each branching is determined to give a total $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$, being a given value from outside of the PS. In this PS, each event has a different weight according to the splitting functions and initial distribution of the PDF. When the number of branchings is $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}`$, the weight of this event ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{W}$}`$) can be obtained as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{W}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{W}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}}}{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\overset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{P}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{P}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{D}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (32)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{W}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\overset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{P}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (33)
where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{P}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ is a splitting function, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{P}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{D}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ is parton distribution function, and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$ is a momentum fraction before the parton-shower evolution. The momentum fraction after the evolution is $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}`$. In order to take this momentum fraction ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$) as an independent variable given by users, unity

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}{\displaystyle \colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\delta }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\overset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝑑}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}`$ (34)
is multiplied to Eq.(32) and integration by $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$ is done first as

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{W}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{W}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}}}{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\overset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{P}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{P}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{D}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\delta }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\overset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (35)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{W}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}}}{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\overset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{P}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{P}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{D}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (36)
$`\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{.}$}`$ (37)
If we employ an appropriate transformation of the input random numbers, we can expect a sufficient efficiency.

A numerical test of the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$-deterministic PS is done by comparing results with the PDF of Cteq5L. The initial distributions of partons are taken from Cteq5L at the energy scale of the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}`$-quark mass. Then, after evolution to an energy scale of 100 GeV, the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$ distributions from the PS are compared to those in Cteq5L. One can see in Fig.2 the PS reproduced $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$ distributions in Cteq5L for both of singlet and non-singlet partons.

Our procedure to generate QCD events involves the following three steps, :

1. numerical integration is performed in the phase-space including $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$ integration,
2. the probability for the existence of each parton in the proton at some low-energy scale is obtained from the PDF, and
3. the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$-deterministic PS performs evolution from a low energy scale to the energy scale of a hard scattering while generating multi-partons with a finite transverse momentum,

instead of steps 2, 3, and 4 in the conventional procedure.

### 5.3 Negative weight treatment

Though the LL-subtraction method makes the negative-positive cancellation in the collinear region very mild, the NLO amplitude may still have negative values in some region in the phase-space. This negative weight is unavoidable due to the perturbative calculation truncated at a finite order of expansion. Especially this problem is very serious in QCD due to its large expansion coefficient. Even if there are events with negative weight, experimentally measurable distributions can be positive with a realistic detector resolution (or a bin-width of histograms). Then we decided to employ a ’quasi-unweighted’ event generation as follows.

The BASES system can treat integrand with negative values and prepare a table of a probability density according to the absolute value of the integrand. The SPRING will generate events based on that table with unit-weight with sign $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}`$ or $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}`$ according to a sign of integrand at given phase point. This sign will be kept during detector simulations and user analysis programs. Finally some bins of histograms will be incremented or decremented by unit amount due to their sign.

## 6 A test of the LL-subtraction method

The LL-subtraction method is numerically tested for a process of $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}`$ in proton-anti proton collision at the CM energy of 2 TeV. Here only non-singlet $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}`$-quark is used in this test to avoid additional complex problem. Matrix elements of the Born process and the real radiation process ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{l}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{o}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{n}$}`$) are generated by GRACE. Numerical integration is done using a BASES system. A non-singlet $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}`$-quark distribution at the energy scale of 4.6 GeV is taken from CTEQ5L. The parton evolution from this energy scale to that of hard-scattering, i.e. an invariant mass of the muon-pair ($`\widehat{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}}`$), is done using the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$-deterministic PS for both of the Born and radiative processes. An additional cut of $`\sqrt{\widehat{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{>}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{40}$}`$ GeV is also applied.

The distributions of transverse momenta of gluons are shown in Figure 2. For the Born process, the largest $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ from the PS is filled to histograms if it has $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ greater than 1 GeV. For the radiative process, that from matrix elements is filled when it has passed the same cut as mentioned above. If one does not take care of double-counting of the LL-terms in the matrix elements and the PS, the distributions of $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{t}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{/}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ show larger values compared with the PS in the Born process, as shown in the left histogram of Figure 2. In order to avoid double-counting, we required that $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ of the matrix elements be greater than those of any gluons from the PS in the radiative process. This double-counting rejection leads to a good agreement of the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ distribution around the low-$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ region, as shown in the middle histogram. However, the PS still emits an insufficient number of gluons around a very high-$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ region. If the upper bound of the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ integration is set to be $`\widehat{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}}`$ in the calculation with the exact matrix elements, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ distributions shows a good agreement, as shown in right histogram. This behavior tells us that the parton-shower can reproduce gluon distributions only for the energy scale less than $`\widehat{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}}`$. This result confirmed the necessity of the step function ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Theta }}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{B}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$) in Eq.(31).

The distributions from $`\stackrel{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{~}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\sigma }$}}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{u}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}}`$ combined with those of the PS on the Born process are compared with those from exact matrix-elements with double-counting rejection. In both of the gluon energy and the transverse-momentum distributions, the LL-subtraction method can give consistent results with those of the exact matrix elements, as shown in Figure 3.

## 7 Example:Drell-Yan processes

In order to demonstrate how our proposed method can really work for physical processes in hadron colliders, an example of an event generator for the Drell-Yan process in a proton anti-proton ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}}`$) collider of TEVATRON at the next-to-leading order QCD is constructed. All possible Feynman diagram contributing to the process $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}`$(+additional parton) as shown in figure 5 are taken into account in the generator. All quarks and gluons except top-quark in an initial (anti-)proton are summed up using the PDF of Cteq5L at $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}`$-quark mass of 4.5 GeV and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ evolution to the energy scale of a hard scattering has been done using the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$-deterministic parton shower at leading logarithmic order.

FORTRAN programs to calculate matrix-elements are automatically generated using GRACE including loop- and real-radiation diagram as well as tree ones. All quark masses are set to be massless in the matrix-element calculations. A QCD coupling constant is calculated using a leading-order formula with $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Lambda }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{C}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{D}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{146}$}`$ MeV which is obtained by parameter fittings in Cteq5L. The factorization scale is set to be equal to the energy of a $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}`$-pair system. As it is known well for the Drell-Yan process, the renormalization energy scale does not explicitly appear in the correction terms. Main correction terms are absorbed in the running coupling constant and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ evolution of the parton distribution. Remaining terms with the large logarithm ($`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{ln}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{/}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Z}$}`$) is canceled out after summing up virtual and collinear correction terms. This renormalization scale independence is numerically checked in our program. Infrared divergent part proportional to $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{/}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{/}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\epsilon }$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}}`$ are left in the program and numerically confirmed that these divergent coefficients are canceled out after summing up virtual and collinear parts of the matrix elements at the order of $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{10}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{20}$}}`$. No experimental cut except the energy of a $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}`$-pair system to be greater than 10 GeV is applied. The total cross section of the Drell-Yan process at 2 TeV of the $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}}`$ CM energy with above cut is obtained to be $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1.026}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\times }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{10}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}`$ pb at the tree level and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1.288}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\times }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{10}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}`$ pb at the NLO, which gives a $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{K}$}`$-factor of 1.256. In the cross section of $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1.288}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\times }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{10}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}`$ pb at the NLO, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1.807}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\times }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{10}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ pb comes from the visible cross section after the LL-subtraction. It is rather good behavior as perturbative correction of QCD. Numerical results of the virtual and collinear correction of $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}\overline{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}}`$ initial state are compared with those based on formulae given by Altarelli et al. and are confirmed consistent. It is confirmed that the cross sections are independent from the value of $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{Q}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$ defined in section 4.2 and agree with those with simple phase-space slicing method as expected.

Next let us look at transverse energy distributions of the largest $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{E}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ jet and a $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}`$-pair system. Origin of the transverse energy has two sources, one is finite transverse momentum of the PS and the other is that from real-radiation matrix-elements. In the LL-subtraction method, these two kinds of distributions can be smoothly combined without any double-counting. As shown in figure 6, contribution from the parton shower applied on the NLO matrix elements (solid histograms) is dominated around small $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{E}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ region. On the other hand a high $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{E}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ events are dominated by the real-radiation correction (dashed histograms). These two regions are smoothly connected as shown by stars (\*) which show the sum of solid and dashed histograms. If one applied PDF simply on the the real-radiation matrix-elements without any care of double counting, small $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{E}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ events are clearly over estimated as shown by dot-dashed histograms. On the other hand, high $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{E}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}`$ events are consistent with those from LL-subtraction method, because they are free from both double-counting and smearing due to the PS.

## 8 Conclusions

A new method to construct event-generators based on next-to-leading order QCD matrix-elements and an $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{x}$}`$-deterministic parton shower is proposed. Matrix elements of loop diagram as well as those of a tree level can be generated by the GRACE system. A soft/collinear singularity is treated using the leading-log subtraction method. It has been demonstrated that the LL-subtraction method can give good agreement with the exact matrix elements without any double-counting problem. The NLO event generator for the Drell-Yan processes has been constructed based on our method and showed smooth distributions of transverse momenta combining the PS and exact matrix-elements of real radiation processes without any double-counting problem.

Authors would like to thank Dr. S. Idalia for continuous discussions on this subject with him and his useful suggestions.

This work was supported in part by the Ministry of Education, Science and Culture under the Grant-in-Aid No. 11206203 and 14340081.

Appendix

## Appendix A Effective vertices

### A.1 $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{S}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{U}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ Color factor

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝐓}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{:}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{the}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{fundamental}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{representation}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{:}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{the}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{structure}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{constants}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{[}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝐓}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝐓}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{]}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝐓}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝐓}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝐓}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{C}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{C}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{r}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝐓}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝐓}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{T}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{R}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\delta }$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$
$`{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{C}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{G}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\delta }$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\delta }$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{.}$}`$
### A.2 tree vertex

* quark-gluon vertex: $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝐓}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}`$
* three-gluon vertex: $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}}`$
* four-gluon vertex: $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{b}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{d}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{e}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{cyclic}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{permutation}}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$

### A.3 quark self-energy

### A.4 gluon self-energy

### A.5 quark-gluon vertex

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Lambda }}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝐓}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{C}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{F}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{C}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{G}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\gamma }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathit{}}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝐓}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{C}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{G}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\gamma }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{I}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathit{}}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$
#### A.5.1 case I

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$, and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{ln}}$}\frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}`$, where $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ or $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$.

#### A.5.2 case II

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{p}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}`$, and $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{ln}}$}\frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}`$.

### A.6 three gluon vertex

$$\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{0}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{L}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{ln}}$}\frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}$$

$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mathrm{\Lambda }}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{a}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\alpha }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{s}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\pi }$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left[}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{C}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{G}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\left(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒢}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒢}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒢}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right)}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{+}$}{\displaystyle \frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{N}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{f}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒢}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\right]}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{,}$}`$ (38)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒢}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}`$ $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}`$ $`{\displaystyle \underset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}{\overset{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}}}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{c}$}_{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{i}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒫}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{j}$}`$ (39)
$`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒢}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}`$:gluon loop, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒢}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}`$:ghost loop, $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒢}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}`$:gluon loop (fish type), $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒢}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{4}$}`$:quark loop

Independent momenta

* $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒫}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}`$
* $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒫}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{g}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}`$
* $`\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{𝒫}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{=}$}\frac{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{(}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{1}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{k}$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{)}$}^{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{\mu }$}_\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{3}$}}}{\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{q}$}^\colorbox[rgb]{1,1,1}{$\textcolor[rgb]{0,0,0}{2}$}}`$