5 * Revision 1.1.1.2 1996/10/30 08:28:43 cernlib
8 * Revision 1.1.1.1 1996/04/12 15:29:24 plothow
12 #include "pdf/pilot.h"
13 C. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
15 C. * G R V - P A R A M E T R I Z A T I O N S *
17 C. * adapted by H. Plothow-Besch August 1991 *
20 C. * FOR A DETAILED EXPLANATION SEE : *
21 C. * M. GLUECK, E.REYA, A.VOGT: DO-TH 91/07 *
23 C. * THE PARAMETRIZATIONS ARE FITTED TO OUR PARTON DISTRIBUTIONS *
24 C. * FOR Q ** 2 BETWEEN MU ** 2 (= 0.25 / 0.30 GEV ** 2 IN LO *
25 C. * / HO) AND 1.E8 GEV ** 2 AND FOR X BETWEEN 1.E-5 AND 1. *
26 C. * REGIONS, WHERE THE DISTRIBUTION UNDER CONSIDERATION IS NEG- *
27 C. * LIGIBLE, I.E. BELOW ABOUT 1.E-4, WERE EXCLUDED FROM THE FIT. *
29 C. * HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
30 C. * M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
32 C. * CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS : *
33 C. * LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, *
34 C. * LAMBDA(5) = 0.153, LAMBDA(6) = 0.082 GEV *
35 C. * HO : LAMBDA(3) = 0.248, LAMBDA(4) = 0.200, *
36 C. * LAMBDA(5) = 0.131, LAMBDA(6) = 0.053 GEV *
38 C. * HO DISTRIBUTION REFER TO THE MS-BAR SCHEME OF BARDEEN ET AL. *
40 C. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
42 C...LO PARAMETRIZATIONS :
43 SUBROUTINE GRVLO (ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL)
45 #include "pdf/expdp.inc"
46 + ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL
53 S = ALOG (ALOG( MAX( MU2,Q2)/LAM2) / ALOG(MU2/LAM2))
57 NUD = 0.663 + 0.191 * S - 0.041 * S2 + 0.031 * S3
59 AGUD = -1.97 + 6.74 * S - 1.96 * S2
60 BUD = 24.4 - 20.7 * S + 4.08 * S2
61 DUD = 2.86 + 0.70 * S - 0.02 * S2
62 UDV = GRVFV (X, NUD, AKUD, AGUD, BUD, DUD)
64 ND = 0.579 + 0.283 * S + 0.047 * S2
65 AKD = 0.523 - 0.015 * S
66 AGD = 2.22 - 0.59 * S - 0.27 * S2
67 BD = 5.95 - 6.19 * S + 1.55 * S2
68 DD = 3.57 + 0.94 * S - 0.16 * S2
69 DV = GRVFV (X, ND, AKD, AGD, BD, DD)
79 AGG = 0.0 + 4.879 * S - 1.383 * S2
80 BGG = 25.92 - 28.97 * S + 5.596 * S2
81 CG = -25.69 + 23.68 * S - 1.975 * S2
82 DG = 2.537 + 1.718 * S + 0.353 * S2
83 EG = 0.595 + 2.138 * S
85 GL = GRVFW (X, S, ALG, BEG, AKG, BKG, AGG, BGG, CG, DG, EG, ESG)
87 C...X * UBAR = X * DBAR :
90 AKU = 0.412 - 0.171 * S
91 BKU = 0.566 - 0.496 * S
94 CU = 1.029 + 1.785 * S - 0.459 * S2
95 DU = 4.696 + 2.109 * S
96 EU = 3.838 + 1.944 * S
98 UDB = GRVFW (X, S, ALU, BEU, AKU, BKU, AGU, BGU, CU, DU, EU, ESU)
100 C...X * SBAR = X * S :
104 AKS = 2.082 - 0.577 * S
105 AGS = -3.055 + 1.024 * S ** 0.67
106 BS = 27.4 - 20.0 * S ** 0.154
108 EST = 4.33 + 1.408 * S
109 ESS = 8.27 - 0.437 * S
110 SB = GRVFWS (X, S, SS, ALS, BES, AKS, AGS, BS, DS, EST, ESS)
112 C...X * CBAR = X * C :
118 BC = 4.24 - 0.804 * S
119 DC = 3.46 + 1.076 * S
120 EC = 4.61 + 1.490 * S
121 ESC = 2.555 + 1.961 * S
122 CB = GRVFWS (X, S, SC, ALC, BEC, AKC, AGC, BC, DC, EC, ESC)
124 C...X * BBAR = X * B :
131 DB = 2.929 + 1.396 * S
132 EB = 4.71 + 1.514 * S
133 ESB = 4.02 + 1.239 * S
134 BB = GRVFWS (X, S, SBO, ALB, BEB, AKB, AGB, BBO, DB, EB, ESB)
136 C...X * TBAR = X * T :