11 IMPLICIT DOUBLE PRECISION(a-
h, o-
z)
15 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
16 common/
pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
21 dimension xppr(-6:6),q2min(16)
22 DATA q2min/ 2.56d0, 2.56d0, 2.56d0, 0.4d0, 0.4d0, 0.4d0,
23 &1.0d0, 1.0d0, 2*0d0, 0.25d0, 5d0, 5d0, 4d0, 4d0, 0d0/
32 IF(nset.EQ.9.OR.nset.EQ.10) nset=6
34 IF(
mstp(57).EQ.0)
THEN
37 q2l=
max(q2min(nset),q2)
40 IF(nset.GE.1.AND.nset.LE.3)
THEN
42 qrt=sqrt(
max(1d0,q2l))
55 ELSEIF(nset.GE.4.AND.nset.LE.6)
THEN
58 CALL
pygrvl(
x, q2l, uv, dv, del, udb, sb, chm, bot, gl)
59 ELSEIF(nset.EQ.5)
THEN
60 CALL
pygrvm(
x, q2l, uv, dv, del, udb, sb, chm, bot, gl)
62 CALL
pygrvd(
x, q2l, uv, dv, del, udb, sb, chm, bot, gl)
67 xppr(-1)=0.5d0*(udb+del)
68 xppr(-2)=0.5d0*(udb-del)
78 ELSEIF(nset.EQ.7)
THEN
82 qrt=sqrt(
max(1d0,min(1d8,q2l)))
83 xin=
max(1d-6,min(1d0,
x))
94 xppr(
i)=xin*sumudb*ratudb/(1d0+ratudb)
96 xppr(
i)=xin*sumudb/(1d0+ratudb)
99 IF(
i.LT.0) xppr(-
i)=xppr(
i)
103 ELSEIF(nset.EQ.8)
THEN
105 qrt=sqrt(
max(1d0,min(1d8,q2l)))
106 xin=
max(1d-6,min(1d0,
x))
117 xppr(
i)=xin*sumudb*ratudb/(1d0+ratudb)
119 xppr(
i)=xin*sumudb/(1d0+ratudb)
122 IF(
i.LT.0) xppr(-
i)=xppr(
i)
126 ELSEIF(nset.GE.11.AND.nset.LE.15)
THEN
132 ELSEIF(nset.EQ.16)
THEN