]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:28:44 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 | SUBROUTINE GRVPILO (ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL) | |
14 | * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * | |
15 | * * | |
16 | * G R V - P I 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/16 * | |
20 | * * | |
21 | * THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS * | |
22 | * FOR Q ** 2 BETWEEN MU ** 2 (= 0.25 / 0.30 GEV ** 2 IN LO * | |
23 | * / HO) AND 1.E8 GEV ** 2 AND FOR X BETWEEN 1.E-5 AND 1. * | |
24 | * REGIONS, WHERE THE DISTRIBUTION UNDER CONSIDERATION IS NEG- * | |
25 | * LIGIBLE, I.E. BELOW ABOUT 1.E-4, WERE EXCLUDED FROM THE FIT. * | |
26 | * * | |
27 | * HEAVY QUARK THRESHOLDS Q(H) = M(H) : * | |
28 | * M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV * | |
29 | * * | |
30 | * CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS : * | |
31 | * LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, * | |
32 | * LAMBDA(5) = 0.153, LAMBDA(6) = 0.082 GEV * | |
33 | * HO : LAMBDA(3) = 0.248, LAMBDA(4) = 0.200, * | |
34 | * LAMBDA(5) = 0.131, LAMBDA(6) = 0.053 GEV * | |
35 | * * | |
36 | * HO DISTRIBUTION REFER TO THE MS-BAR SCHEME OF BARDEEN ET AL. * | |
37 | * * | |
38 | * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * | |
39 | C | |
40 | IMPLICIT REAL (A - Y) | |
41 | #include "pdf/expdp.inc" | |
42 | + ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL | |
43 | REAL X, Q | |
44 | X = ZX | |
45 | Q = ZQ | |
46 | MU2 = 0.25 | |
47 | LAM2 = 0.232 * 0.232 | |
48 | Q2 = Q*Q | |
49 | S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2)) | |
50 | DS = SQRT (S) | |
51 | S2 = S * S | |
52 | C...X * VALENCE : | |
53 | NV = 0.519 + 0.180 * S - 0.011 * S2 | |
54 | AKV = 0.499 - 0.027 * S | |
55 | AGV = 0.381 - 0.419 * S | |
56 | DV = 0.367 + 0.563 * S | |
57 | VAP = GRVFVP (X, NV, AKV, AGV, DV) | |
58 | ZUV = VAP | |
59 | ZDV = ZUV | |
60 | C...X * GLUON : | |
61 | ALG = 0.599 | |
62 | BEG = 1.263 | |
63 | AKG = 0.482 + 0.341 * DS | |
64 | BKG = 0.0 | |
65 | AGG = 0.678 + 0.877 * S - 0.175 * S2 | |
66 | BGG = 0.338 - 1.597 * S | |
67 | CG = 0.0 - 0.233 * S + 0.406 * S2 | |
68 | DG = 0.390 + 1.053 * S | |
69 | EG = 0.618 + 2.070 * S | |
70 | ESG = 3.676 | |
71 | GLP = GRVFGP(X,S, ALG, BEG, AKG, BKG, AGG, BGG, CG, DG, EG, ESG) | |
72 | ZGL = GLP | |
73 | C...X * QBAR (SU(3)-SYMMETRIC SEA) : | |
74 | SL = 0.0 | |
75 | ALS = 0.55 | |
76 | BES = 0.56 | |
77 | AKS = 2.538 - 0.763 * S | |
78 | AGS = -0.748 | |
79 | BS = 0.313 + 0.935 * S | |
80 | DS = 3.359 | |
81 | EST = 4.433 + 1.301 * S | |
82 | ESS = 9.30 - 0.887 * S | |
83 | QBP = GRVFQBP (X, S, SL, ALS, BES, AKS, AGS, BS, DS, EST, ESS) | |
84 | ZUDB = QBP | |
85 | ZSB = ZUDB | |
86 | C...X * CBAR = X * C : | |
87 | SC = 0.888 | |
88 | ALC = 1.02 | |
89 | BEC = 0.39 | |
90 | AKC = 0.0 | |
91 | AGC = 0.0 | |
92 | BC = 1.008 | |
93 | DC = 1.208 + 0.771 * S | |
94 | EC = 4.40 + 1.493 * S | |
95 | ESC = 2.032 + 1.901 * S | |
96 | CBP = GRVFQBP (X, S, SC, ALC, BEC, AKC, AGC, BC, DC, EC, ESC) | |
97 | ZCB = CBP | |
98 | C...X * BBAR = X * B : | |
99 | SBO = 1.351 | |
100 | ALB = 1.03 | |
101 | BEB = 0.39 | |
102 | AKB = 0.0 | |
103 | AGB = 0.0 | |
104 | BBO = 0.0 | |
105 | DB = 0.697 + 0.855 * S | |
106 | EB = 4.51 + 1.490 * S | |
107 | ESB = 3.056 + 1.694 * S | |
108 | BBP = GRVFQBP (X, S, SBO, ALB, BEB, AKB, AGB, BBO, DB, EB, ESB) | |
109 | ZBB = BBP | |
110 | C...X * TBAR = X * T : | |
111 | TBP = 0. | |
112 | ZTB = TBP | |
113 | RETURN | |
114 | END |