]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:28:43 cernlib | |
6 | * Version 7.04 | |
7 | * | |
8 | * Revision 1.1.1.1 1996/04/12 15:29:24 plothow | |
9 | * Version 7.01 | |
10 | * | |
11 | * | |
12 | #include "pdf/pilot.h" | |
13 | C. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * | |
14 | C. * * | |
15 | C. * G R V - P A R A M E T R I Z A T I O N S * | |
16 | C. * * | |
17 | C. * adapted by H. Plothow-Besch August 1991 * | |
18 | C. * * | |
19 | C. * * | |
20 | C. * FOR A DETAILED EXPLANATION SEE : * | |
21 | C. * M. GLUECK, E.REYA, A.VOGT: DO-TH 91/07 * | |
22 | C. * * | |
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. * | |
28 | C. * * | |
29 | C. * HEAVY QUARK THRESHOLDS Q(H) = M(H) : * | |
30 | C. * M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV * | |
31 | C. * * | |
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 * | |
37 | C. * * | |
38 | C. * HO DISTRIBUTION REFER TO THE MS-BAR SCHEME OF BARDEEN ET AL. * | |
39 | C. * * | |
40 | C. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * | |
41 | C | |
42 | C...LO PARAMETRIZATIONS : | |
43 | SUBROUTINE GRVLO (ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL) | |
44 | IMPLICIT REAL (A - Y) | |
45 | #include "pdf/expdp.inc" | |
46 | + ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL | |
47 | REAL X, Q | |
48 | X = ZX | |
49 | Q = ZQ | |
50 | MU2 = 0.25 | |
51 | LAM2 = 0.232 * 0.232 | |
52 | Q2 = Q*Q | |
53 | S = ALOG (ALOG( MAX( MU2,Q2)/LAM2) / ALOG(MU2/LAM2)) | |
54 | S2 = S * S | |
55 | S3 = S2 * S | |
56 | C...X * (UV + DV) : | |
57 | NUD = 0.663 + 0.191 * S - 0.041 * S2 + 0.031 * S3 | |
58 | AKUD = 0.326 | |
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) | |
63 | C...X * DV : | |
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) | |
70 | ZDV = DV | |
71 | C...X * UV : | |
72 | UV = UDV - DV | |
73 | ZUV = UV | |
74 | C...X * G : | |
75 | ALG = 0.558 | |
76 | BEG = 1.218 | |
77 | AKG = 1.00 - 0.17 * S | |
78 | BKG = 0.0 | |
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 | |
84 | ESG = 4.066 | |
85 | GL = GRVFW (X, S, ALG, BEG, AKG, BKG, AGG, BGG, CG, DG, EG, ESG) | |
86 | ZGL = GL | |
87 | C...X * UBAR = X * DBAR : | |
88 | ALU = 1.396 | |
89 | BEU = 1.331 | |
90 | AKU = 0.412 - 0.171 * S | |
91 | BKU = 0.566 - 0.496 * S | |
92 | AGU = 0.363 | |
93 | BGU = -1.196 | |
94 | CU = 1.029 + 1.785 * S - 0.459 * S2 | |
95 | DU = 4.696 + 2.109 * S | |
96 | EU = 3.838 + 1.944 * S | |
97 | ESU = 2.845 | |
98 | UDB = GRVFW (X, S, ALU, BEU, AKU, BKU, AGU, BGU, CU, DU, EU, ESU) | |
99 | ZUDB = UDB | |
100 | C...X * SBAR = X * S : | |
101 | SS = 0.0 | |
102 | ALS = 0.803 | |
103 | BES = 0.563 | |
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 | |
107 | DS = 6.22 | |
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) | |
111 | ZSB = SB | |
112 | C...X * CBAR = X * C : | |
113 | SC = 0.888 | |
114 | ALC = 1.01 | |
115 | BEC = 0.37 | |
116 | AKC = 0.0 | |
117 | AGC = 0.0 | |
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) | |
123 | ZCB = CB | |
124 | C...X * BBAR = X * B : | |
125 | SBO = 1.351 | |
126 | ALB = 1.00 | |
127 | BEB = 0.51 | |
128 | AKB = 0.0 | |
129 | AGB = 0.0 | |
130 | BBO = 1.848 | |
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) | |
135 | ZBB = BB | |
136 | C...X * TBAR = X * T : | |
137 | TB=0. | |
138 | ZTB = TB | |
139 | RETURN | |
140 | END |