5 * Revision 1.1.1.2 1996/10/30 08:32:04 cernlib
8 * Revision 1.1.1.1 1996/04/12 15:30:16 plothow
12 #include "pdf/pilot.h"
14 ccccccccccccccccccccccccccccccccccccccccccccccccccccccc
15 DOUBLE PRECISION function WHIT5G(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 WHIT5 gluon (U100)
44 A0g = 1.000000d+01+s*(-3.400000d+01)+s2*( 6.900000d+01)
45 $ +s3*(-7.530000d+01)+s4*( 3.230000d+01)
46 B0g = s*(-1.126000d+00)+s2*( 9.260000d-01)
48 C0g = 9.000000d+00+s*( 4.810000d-01)+s2*( 3.200000d+00)
50 A1g = s*( 6.020000d-01)+s2*( 5.090000d-01)
51 $ +s3*(-2.054000d+00)+s4*( 1.392000d+00)
52 AA1g= s*(-9.220000d-02)+s2*(-1.899000d+00)
53 $ +s3*( 4.180000d+00)+s4*(-2.494000d+00)
54 B1g =-2.895000d-01+s*( 3.760000d-01)+s2*(-1.719000d+00)
56 C1g = 1.439000d+00+s*(-5.570000d-01)+s2*( 3.660000d-01)
57 $ +s3*( 7.330000d-01)+s4*(-7.620000d-01)
61 cccccc WHIT5 gluon (O100)
63 A0g = 1.995000d+00+s*(-3.260000d+00)+s2*( 1.818000d+00)
64 $ +s3*( 1.711000d+00)+s4*(-4.990000d+00)
65 B0g =-4.660000d-01+s*(-6.100000d-01)+s2*( 1.691000d+00)
66 $ +s3*(-6.680000d+00)+s4*( 1.019000d+01)
67 C0g = 1.075000d+01+s*( 5.420000d+00)+s2*( 6.550000d+00)
68 $ +s3*(-2.297000d+01)+s4*( 1.867000d+01)
69 A1g = 2.992000d-01+s*( 1.179000d+00)+s2*(-1.915000d+00)
70 $ +s3*( 7.260000d+00)+s4*(-1.839000d+01)
71 AA1g=-1.600000d-01+s*(-1.114000d+00)+s2*( 2.939000d+00)
72 $ +s3*(-6.660000d+00)+s4*( 1.923000d+01)
73 B1g =-4.830000d-01+s*( 7.550000d-01)+s2*(-3.800000d+00)
74 $ +s3*( 1.075000d+01)+s4*(-1.993000d+01)
75 C1g = 1.297000d+00+s*(-1.669000d-01)+s2*( 1.906000d+00)
79 WHIT5G = prsccf/alinv/x*
80 $ ( A0g * x**B0g * x1**C0g
81 $ +(A1g+AA1g*x) * x**B1g * x1**C1g )