]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:28:40 cernlib | |
6 | * Version 7.04 | |
7 | * | |
8 | * Revision 1.1.1.1 1996/04/12 15:29:23 plothow | |
9 | * Version 7.01 | |
10 | * | |
11 | * | |
12 | #include "pdf/pilot.h" | |
13 | SUBROUTINE GRVGAHO (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL) | |
14 | * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * | |
15 | * * | |
16 | * G R V - P H O T O N - P A R A M E T R I Z A T I O N S * | |
17 | * * | |
18 | * FOR A DETAILED EXPLANATION SEE : * | |
19 | * M. GLUECK, E.REYA, A.VOGT: DO-TH 91/31 * | |
20 | * * | |
21 | * THE OUTPUT IS ALWAYS 1./ ALPHA(EM) * X * PARTON DENSITY * | |
22 | * output modified by HPB to be always X * PARTON DENSITY * | |
23 | * * | |
24 | * THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS * | |
25 | * FOR Q ** 2 BETWEEN MU ** 2 (= 0.25 / 0.30 GEV ** 2 IN LO * | |
26 | * / HO) AND 1.E6 GEV ** 2 AND FOR X BETWEEN 1.E-5 AND 1. * | |
27 | * * | |
28 | * HEAVY QUARK THRESHOLDS Q(H) = M(H) : * | |
29 | * M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV * | |
30 | * * | |
31 | * CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS : * | |
32 | * LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, * | |
33 | * LAMBDA(5) = 0.153, LAMBDA(6) = 0.082 GEV * | |
34 | * HO : LAMBDA(3) = 0.248, LAMBDA(4) = 0.200, * | |
35 | * LAMBDA(5) = 0.131, LAMBDA(6) = 0.053 GEV * | |
36 | * * | |
37 | * HO DISTRIBUTIONS REFER TO THE DIS(GAMMA) SCHEME, SEE : * | |
38 | * M. GLUECK, E.REYA, A.VOGT: DO-TH 91/26 * | |
39 | * * | |
40 | * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * | |
41 | C | |
42 | IMPLICIT REAL (A - Y) | |
43 | #include "pdf/expdp.inc" | |
44 | + ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL | |
45 | DATA ALPHEM/7.29927D-3/ | |
46 | REAL X, Q | |
47 | X = ZX | |
48 | Q = ZQ | |
49 | MU2 = 0.3 | |
50 | LAM2 = 0.248 * 0.248 | |
51 | Q2 = Q*Q | |
52 | S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2)) | |
53 | SS = SQRT (S) | |
54 | S2 = S * S | |
55 | C...X * U = X * UBAR : | |
56 | AL = 0.583 | |
57 | BE = 0.688 | |
58 | AK = 0.449 - 0.025 * S - 0.071 * S2 | |
59 | BK = 5.060 - 1.116 * SS | |
60 | AG = 0.103 | |
61 | BG = 0.319 + 0.422 * S | |
62 | C = 1.508 + 4.792 * S - 1.963 * S2 | |
63 | D = 1.075 + 0.222 * SS - 0.193 * S2 | |
64 | E = 4.147 + 1.131 * S | |
65 | ES = 1.661 + 0.874 * S | |
66 | UH = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES) | |
67 | ZUV = UH * ALPHEM | |
68 | ZUB = ZUV | |
69 | C...X * D = X * DBAR : | |
70 | AL = 0.591 | |
71 | BE = 0.698 | |
72 | AK = 0.442 - 0.132 * S - 0.058 * S2 | |
73 | BK = 5.437 - 1.916 * SS | |
74 | AG = 0.099 | |
75 | BG = 0.311 - 0.059 * S | |
76 | C = 0.800 + 0.078 * S - 0.100 * S2 | |
77 | D = 0.862 + 0.294 * SS - 0.184 * S2 | |
78 | E = 4.202 + 1.352 * S | |
79 | ES = 1.841 + 0.990 * S | |
80 | DH = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES) | |
81 | ZDV = DH * ALPHEM | |
82 | ZDB = ZDV | |
83 | C...X * G : | |
84 | AL = 1.161 | |
85 | BE = 1.591 | |
86 | AK = 0.530 - 0.742 * SS + 0.025 * S2 | |
87 | BK = 5.662 | |
88 | AG = 0.533 - 0.281 * SS + 0.218 * S2 | |
89 | BG = 0.025 - 0.518 * S + 0.156 * S2 | |
90 | C = -0.282 + 0.209 * S2 | |
91 | D = 0.107 + 1.058 * S - 0.218 * S2 | |
92 | E = 0.0 + 2.704 * S | |
93 | ES = 3.071 - 0.378 * S | |
94 | GH = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES) | |
95 | ZGL = GH * ALPHEM | |
96 | C...X * S = X * SBAR : | |
97 | SF = 0.0 | |
98 | AL = 0.635 | |
99 | BE = 0.456 | |
100 | AK = 1.770 - 0.735 * SS - 0.079 * S2 | |
101 | BK = 3.832 | |
102 | AG = 0.084 - 0.023 * S | |
103 | BG = 0.136 | |
104 | C = 2.119 - 0.942 * S + 0.063 * S2 | |
105 | D = 1.271 + 0.076 * S - 0.190 * S2 | |
106 | E = 4.604 + 0.737 * S | |
107 | ES = 1.641 + 0.976 * S | |
108 | SH = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES) | |
109 | ZSB = SH * ALPHEM | |
110 | C...X * C = X * CBAR : | |
111 | SF = 0.820 | |
112 | AL = 0.926 | |
113 | BE = 0.152 | |
114 | AK = 1.142 - 0.175 * S | |
115 | BK = 3.276 | |
116 | AG = 0.504 + 0.317 * S | |
117 | BG = -0.433 | |
118 | C = 3.334 | |
119 | D = 0.398 + 0.326 * S - 0.107 * S2 | |
120 | E = 5.493 + 0.408 * S | |
121 | ES = 2.426 + 1.277 * S | |
122 | CH = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES) | |
123 | ZCB = CH * ALPHEM | |
124 | C...X * B = X * BBAR : | |
125 | SF = 1.297 | |
126 | AL = 0.969 | |
127 | BE = 0.266 | |
128 | AK = 1.953 - 0.391 * S | |
129 | BK = 1.657 - 0.161 * S | |
130 | AG = 1.076 + 0.034 * S | |
131 | BG = -2.015 | |
132 | C = 1.662 | |
133 | D = 0.353 + 0.016 * S | |
134 | E = 5.713 + 0.249 * S | |
135 | ES = 3.456 + 0.673 * S | |
136 | BH = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES) | |
137 | ZBB = BH * ALPHEM | |
138 | c | |
139 | RETURN | |
140 | END |