5 * Revision 1.1.1.2 1996/10/30 08:28:40 cernlib
8 * Revision 1.1.1.1 1996/04/12 15:29:23 plothow
12 #include "pdf/pilot.h"
13 SUBROUTINE GRVGAHO (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL)
14 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
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 *
18 * FOR A DETAILED EXPLANATION SEE : *
19 * M. GLUECK, E.REYA, A.VOGT: DO-TH 91/31 *
21 * THE OUTPUT IS ALWAYS 1./ ALPHA(EM) * X * PARTON DENSITY *
22 * output modified by HPB to be always X * PARTON DENSITY *
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. *
28 * HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
29 * M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
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 *
37 * HO DISTRIBUTIONS REFER TO THE DIS(GAMMA) SCHEME, SEE : *
38 * M. GLUECK, E.REYA, A.VOGT: DO-TH 91/26 *
40 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
43 #include "pdf/expdp.inc"
44 + ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL
45 DATA ALPHEM/7.29927D-3/
52 S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
55 C...X * U = X * UBAR :
58 AK = 0.449 - 0.025 * S - 0.071 * S2
59 BK = 5.060 - 1.116 * SS
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
65 ES = 1.661 + 0.874 * S
66 UH = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
69 C...X * D = X * DBAR :
72 AK = 0.442 - 0.132 * S - 0.058 * S2
73 BK = 5.437 - 1.916 * SS
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
79 ES = 1.841 + 0.990 * S
80 DH = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
86 AK = 0.530 - 0.742 * SS + 0.025 * S2
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
93 ES = 3.071 - 0.378 * S
94 GH = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
96 C...X * S = X * SBAR :
100 AK = 1.770 - 0.735 * SS - 0.079 * S2
102 AG = 0.084 - 0.023 * S
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)
110 C...X * C = X * CBAR :
114 AK = 1.142 - 0.175 * S
116 AG = 0.504 + 0.317 * S
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)
124 C...X * B = X * BBAR :
128 AK = 1.953 - 0.391 * S
129 BK = 1.657 - 0.161 * S
130 AG = 1.076 + 0.034 * S
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)