]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:32:01 cernlib | |
6 | * Version 7.04 | |
7 | * | |
8 | * Revision 1.1.1.1 1996/04/12 15:30:15 plothow | |
9 | * Version 7.01 | |
10 | * | |
11 | * | |
12 | #include "pdf/pilot.h" | |
13 | c | |
14 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccc | |
15 | DOUBLE PRECISION function WHIT3G(x,Q2) | |
16 | c input: x,Q2 | |
17 | c output: clg | |
18 | c (gluon dist.) | |
19 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccc | |
20 | implicit none | |
21 | c arg | |
22 | DOUBLE PRECISION Q2,x | |
23 | c const | |
24 | DOUBLE PRECISION q42it,q52it,lam42,lam52 | |
25 | DOUBLE PRECISION alinv | |
26 | c local | |
27 | DOUBLE PRECISION A0g,B0g,C0g,A1g,AA1g,B1g,C1g | |
28 | DOUBLE PRECISION s,s2,s3,s4,prsccf | |
29 | DOUBLE PRECISION x1 | |
30 | c parameters | |
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 | |
35 | c | |
36 | c begin | |
37 | x1=1.0d0-x | |
38 | c | |
39 | if(Q2.lt.100.0d0) then | |
40 | c under 100 GeV^2 | |
41 | c | |
42 | cccccc WHIT3 gluon (U100) | |
43 | c | |
44 | A0g = 8.000000d+00+s*(-2.864000d+01)+s2*( 5.590000d+01) | |
45 | $ +s3*(-5.760000d+01)+s4*( 2.366000d+01) | |
46 | B0g = s*(-9.870000d-01)+s2*( 5.100000d-01) | |
47 | $ +s3*(-6.670000d-02) | |
48 | C0g = 1.500000d+01+s*( 3.310000d-01)+s2*( 3.500000d+00) | |
49 | $ +s3*( 8.920000d-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) | |
55 | $ +s3*(-5.070000d-01) | |
56 | C1g = 1.244000d+00+s*( 5.880000d-01)+s2*(-1.228000d+00) | |
57 | $ +s3*( 8.090000d-01) | |
58 | else | |
59 | c over 100 GeV^2 | |
60 | c | |
61 | cccccc WHIT3 gluon (O100) | |
62 | c | |
63 | A0g = 1.270000d+00+s*(-2.817000d+00)+s2*( 5.740000d+00) | |
64 | $ +s3*(-1.327000d+01)+s4*( 1.268000d+01) | |
65 | B0g =-4.610000d-01+s*(-8.170000d-01)+s2*( 3.320000d+00) | |
66 | $ +s3*(-1.296000d+01)+s4*( 1.893000d+01) | |
67 | C0g = 1.721000d+01+s*( 1.257000d+00)+s2*( 5.050000d+01) | |
68 | $ +s3*(-2.761000d+02)+s4*( 4.900000d+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) | |
76 | $ +s3*(-2.072000d+00) | |
77 | endif | |
78 | c | |
79 | WHIT3G = prsccf/alinv/x* | |
80 | $ ( A0g * x**B0g * x1**C0g | |
81 | $ +(A1g+AA1g*x) * x**B1g * x1**C1g ) | |
82 | c | |
83 | return | |
84 | end |