5 * Revision 1.1.1.2 1996/10/30 08:32:35 cernlib
8 * Revision 1.2 1996/10/28 15:02:12 pdflib
9 * PDFLIB new version 7.04
11 * Revision 1.1.1.1 1996/04/12 15:29:23 plothow
15 #include "pdf/pilot.h"
16 SUBROUTINE GRVGAH0 (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL)
17 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
19 * 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 *
21 * FOR A DETAILED EXPLANATION SEE : *
22 * M. GLUECK, E.REYA, A.VOGT: DO-TH 91/31 *
24 * THE OUTPUT IS ALWAYS 1./ ALPHA(EM) * X * PARTON DENSITY *
25 * output modified by HPB to be always X * PARTON DENSITY *
27 * THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS *
28 * FOR Q ** 2 BETWEEN MU ** 2 (= 0.25 / 0.30 GEV ** 2 IN LO *
29 * / HO) AND 1.E6 GEV ** 2 AND FOR X BETWEEN 1.E-5 AND 1. *
31 * HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
32 * M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
34 * CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS : *
35 * LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, *
36 * LAMBDA(5) = 0.153, LAMBDA(6) = 0.082 GEV *
37 * HO : LAMBDA(3) = 0.248, LAMBDA(4) = 0.200, *
38 * LAMBDA(5) = 0.131, LAMBDA(6) = 0.053 GEV *
40 * HO DISTRIBUTIONS REFER TO THE DIS(GAMMA) SCHEME, SEE : *
41 * M. GLUECK, E.REYA, A.VOGT: DO-TH 91/26 *
43 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
46 #include "pdf/expdp.inc"
47 + ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL
49 DATA ALPHEM/7.29927D-3/
55 S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
58 C...X * U = X * UBAR :
61 AK = 0.527 + 0.200 * S - 0.107 * S2
62 BK = 7.106 - 0.310 * SS - 0.786 * S2
63 AG = 0.197 + 0.533 * S
64 BG = 0.062 - 0.398 * S + 0.109 * S2
65 C = 0.755 * S - 0.112 * S2
68 ES = 1.752 + 0.866 * S
69 U0 = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
72 C...X * D = X * DBAR :
75 AK = 0.500 + 0.067 * SS - 0.055 * S2
76 BK = 0.376 - 0.453 * SS + 0.405 * S2
77 AG = 0.156 + 0.184 * S
78 BG = 0.0 - 0.528 * S + 0.146 * S2
80 D = 0.379 - 0.301 * S + 0.081 * S2
82 ES = 1.645 + 1.016 * S
83 D0 = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
89 AK = 0.537 - 0.600 * SS
90 BK = 6.389 - 0.953 * S2
91 AG = 0.558 - 0.383 * SS + 0.261 * S2
93 C = -0.222 + 0.078 * S2
94 D = 0.153 + 0.978 * S - 0.209 * S2
96 ES = 3.331 + 0.806 * S
97 G0 = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
99 C...X * S = X * SBAR :
103 AK = 0.622 + 0.332 * S - 0.300 * S2
105 AG = 0.211 - 0.064 * SS - 0.018 * S2
106 BG = -0.215 + 0.122 * S
108 D = 0.0 + 0.253 * S - 0.081 * S2
109 E = 3.990 + 2.014 * S
110 ES = 1.720 + 0.986 * S
111 S0 = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
113 C...X * C = X * CBAR :
117 AK = 1.228 - 0.231 * S
118 BK = 3.806 - 0.337 * S2
119 AG = 0.932 + 0.150 * S
122 D = 0.0 + 0.138 * S - 0.028 * S2
123 E = 5.588 + 0.628 * S
124 ES = 2.665 + 1.054 * S
125 C0 = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
127 C...X * B = X * BBAR :
131 AK = 1.719 - 0.292 * S
132 BK = 0.928 + 0.096 * S
133 AG = 0.845 + 0.178 * S
136 D = -0.191 + 0.151 * S
137 E = 6.089 + 0.282 * S
138 ES = 3.379 + 1.062 * S
139 B0 = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)