5 * Revision 1.1.1.2 1996/10/30 08:31:59 cernlib
8 * Revision 1.1.1.1 1996/04/12 15:30:10 plothow
12 #include "pdf/pilot.h"
14 ccccccccccccccccccccccccccccccccccccccccccccccccccccccc
15 DOUBLE PRECISION function WHIT1G(x,Q2)
19 ccccccccccccccccccccccccccccccccccccccccccccccccccccccc
24 DOUBLE PRECISION q42it,q52it,lam42,lam52
25 DOUBLE PRECISION alinv
27 DOUBLE PRECISION A0g,B0g,C0g,A1g,AA1g,B1g,C1g
28 DOUBLE PRECISION s,s2,s3,s4,prsccf
31 parameter(lam42=0.16d0, lam52=0.091411319d0)
32 parameter(Q42IT=4.0d0, Q52IT=100.0d0)
33 parameter(alinv=137.036d0)
34 common /scale/ s,s2,s3,s4,prsccf
39 if(Q2.lt.100.0d0) then
42 cccccc WHIT1 gluon (U100)
44 A0g = 2.000000d+00+s*(-3.280000d+00)+s2*( 2.894000d+00)
45 $ +s3*(-1.561000d+00)+s4*( 8.180000d-01)
46 B0g = s*(-7.610000d-01)+s2*(-4.900000d-02)
48 C0g = 3.000000d+00+s*( 1.586000d+00)+s2*(-9.490000d-01)
50 A1g = s*( 4.610000d-01)+s2*( 1.041000d-01)
51 $ +s3*(-1.753000d-02)+s4*(-2.717000d-01)
52 AA1g= s*( 9.680000d-03)+s2*(-4.170000d-01)
53 $ +s3*(-3.950000d-01)+s4*( 8.430000d-01)
54 B1g =-4.140000d-01+s*(-6.060000d-02)+s2*( 2.847000d-01)
56 C1g = 1.244000d+00+s*( 5.880000d-01)+s2*(-1.228000d+00)
61 cccccc WHIT1 gluon (O100)
63 A0g = 7.840000d-01+s*(-2.238000d+00)+s2*( 1.617000d+01)
64 $ +s3*(-6.250000d+01)+s4*( 8.390000d+01)
65 B0g =-4.030000d-01+s*(-1.307000d+00)+s2*( 8.780000d+00)
66 $ +s3*(-3.580000d+01)+s4*( 5.350000d+01)
67 C0g = 4.450000d+00+s*( 1.027000d+00)+s2*( 4.460000d+01)
68 $ +s3*(-1.600000d+02)+s4*( 1.816000d+02)
69 A1g = 3.010000d-01+s*( 1.275000d+00)+s2*(-1.563000d+00)
70 $ +s3*( 4.100000d+00)+s4*(-1.337000d+01)
71 AA1g=-1.305000d-01+s*(-1.245000d+00)+s2*( 2.438000d+00)
72 $ +s3*(-2.539000d+00)+s4*( 1.273000d+01)
73 B1g =-4.890000d-01+s*( 9.550000d-01)+s2*(-4.400000d+00)
74 $ +s3*( 1.022000d+01)+s4*(-1.713000d+01)
75 C1g = 1.331000d+00+s*(-2.481000d-01)+s2*( 1.950000d+00)
79 WHIT1G = prsccf/alinv/x*
80 $ ( A0g * x**B0g * x1**C0g
81 $ +(A1g+AA1g*x) * x**B1g * x1**C1g )