For x, the points lie between a minimum and a maximum x, evenly distributed in LX
LX = DLOG(X/(1.0D0 - X))
so
X = DEXP(LX)/(1.0D0+DEXP(LX))
For mu, the points lie between a minimum and a maximum mu, evenly distributed in LMU
LMU = DLOG(MU)
so
MU = DEXP(LMU)
DLOG( pdf + 1.0D-16 )
where pdf is the value of the parton distribution function at this point. The tiny value 1.0D-16 is added so that the table can handle zero values for the parton distributions. (We assume that the parton distribution functions are never negative.)
! CTEQ3M parton distributions generated by DES 14 August 1995
! The comment area follows.
! Each comment line starts with an exclamation point as initial character.
! The first line after the comment line has to be PARTON DATA.
PARTON DATA
VERSION = 110496
LAMBDA = 0.1580000000
NFL_LAMBDA = 5
XMIN = 0.1000000000E-04
XMAX = 0.9500000000
N_XPOINTS = 32
MUMIN = 1.600000000
MUMAX = 1000.000000
N_MUPOINTS = 32
11.901 11.896 14.505 11.899 11.902 11.204 -36.841 -36.841 -36.841
12.050 12.048 14.772 12.050 12.052 11.488 10.176 -36.841 -36.841
12.189 12.188 15.003 12.190 12.191 11.728 10.887 -36.841 -36.841
12.319 12.318 15.207 12.320 12.320 11.936 11.313 -36.841 -36.841
.
.
.
There is a header section, followed by data:
DLOG( pdf + 1.0D-16 ).
Each row corresponds to a different (x,mu) lattice point.
The columns corresponds to different partons:
dbar ubar gluon u d s c b tSome parton distribution sets may include only five flavors. In this case, the numbers in the top quark column correspond to
pdf = 0.
That is, the column contains DLOG( 0 + 1.0D-16 ) = -36.841.
The lattice points are written in the order
(x1,mu1), (x1,mu2), (x1,mu3), ... (x2,mu1), (x2,mu2), (x2,mu3), ... (x3,mu1), (x3,mu2), (x3,mu3), ... ... , ... , ... , ...That is, they are written with the code (assuming six flavors)
DO 30 N=1,NX
DO 30 M=1,NMU
WRITE(99,50)F(-2,N,M),F(-1,N,M),F(0,N,M),F(1,N,M),
> F(2,N,M),F(3,N,M),F(4,N,M),F(5,N,M),F(6,N,M)
50 FORMAT(9G12.5)
C
30 CONTINUE
where
F(NPARTON,N,M) = DLOG( pdf + 1.0D-16 )
PARTONIN(NPARTON,X,MU)
returns the parton distribution function for parton number
NPARTON at momentum fraction X and scale
MU. Here NPARTON is given by
dbar ubar gluon u d s c b t -2 -1 0 1 2 3 4 5 6in the case of six flavors. The number of flavors expected from
PARTONIN can be changed and is specified in the line
NFL = 6