
we adopt the Minimal Prescription, but we stress that from a practical point of view, dif-
ferences between these prescriptions become only relevant at extremely large values x [120],
a region where no experimental data is available.
2.2 Numerical implementation
In this section we discuss the numerical implementation of the N-soft threshold resum-
mation described above. For the PDF fits performed in this work, the processes that
we are interested in are DIS (both neutral and charged currents), lepton-pair invariant
mass and rapidity distribution for Drell-Yan production, and inclusive top pair production
cross section. For DIS and Drell-Yan, we use a new version of the public code ResHiggs,
written originally [36] to perform threshold resummation of Higgs inclusive cross section,
including several improvements with respect to standard N-soft resummation, and later
extended [37] to also perform (improved) resummation in the Soft-Collinear Effective The-
ory formalism. Because of the inclusion of additional processes, the new version of this
code changes name from ResHiggs to TROLL, standing for TROLL Resums Only Large-x
Logarithms, publicly available at the webpage [121]. To give continuity with the original
code ResHiggs, the first version of TROLL is v3.0. For top pair production we use the
public code Top++ [7].
The code TROLL is designed to provide only the contribution of the resummation,
while the fixed-order calculation is obtained from a separate code (in our case, the same
FKgenerator code used for the NNPDF3.0 fits). More specifically, the output of TROLL is
∆
j
K
N
k
LL
, defined as the difference between a resummed K-factor at N
j
LO+N
k
LL and a
fixed-order K-factor at N
j
LO, such that
σ
N
j
LO+N
k
LL
= σ
N
j
LO
+ σ
LO
× ∆
j
K
N
k
LL
, (2.11)
where all the cross sections appearing in the above equation are evaluated with a common
N
j
LO+N
k
LL PDF set. Internally, ∆
j
K
N
k
LL
is computed by subtracting off the expansion
of the resummed coefficient up to O(α
j
s
) from the coefficient itself, multiplying this by
the parton luminosity, computing the inverse Mellin transform and finally dividing by the
parton luminosity.
There are several advantages in using the ∆K-factors defined in eq. (2.11): the fixed-
order normalisation is irrelevant, the K-factor is much less sensitive to the input PDFs than
the cross section itself, and finally, since the resummed contribution has the same kinematic
structure as the Born cross section (soft radiation does not change the kinematics), the
effect of phase space constraints like kinematic cuts are correctly taken into account if they
are applied to the LO cross section in eq. (2.11).
We note that in DIS, Target Mass Corrections (TMCs) at next-to-leading twist are
included in the resummation according to the same prescription used in the NNPDF fitting
code [88], which amounts to multiplying the Mellin transforms of the partonic coefficient
functions by an N-dependent factor. No TMCs are included for the fixed-target DY data.
In figure 1 we show the ∆K-factors for the neutral current DIS structure function
F
2
(x, Q), as a function of x, for Q = 2 GeV and Q = 30 GeV. The plot on the left
– 7 –