]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:28:34 cernlib | |
6 | * Version 7.04 | |
7 | * | |
8 | * Revision 1.1.1.1 1996/04/12 15:29:22 plothow | |
9 | * Version 7.01 | |
10 | * | |
11 | * | |
12 | #include "pdf/pilot.h" | |
13 | * | |
14 | *...NLO PARAMETRIZATION (DIS) : | |
15 | * | |
16 | SUBROUTINE GRV94DI (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZTB,ZGL) | |
17 | IMPLICIT REAL (A - Y) | |
18 | C+SEQ, IMPDP. | |
19 | #include "pdf/expdp.inc" | |
20 | + ZX,ZQ,ZUV,ZDV,ZDEL,ZUDB,ZUB,ZDB,ZSB,ZCB,ZBB,ZTB,ZGL | |
21 | REAL X, Q | |
22 | X = ZX | |
23 | Q = ZQ | |
24 | MU2 = 0.34 | |
25 | LAM2 = 0.248 * 0.248 | |
26 | Q2 = Q*Q | |
27 | ARG1 = MAX(MU2,Q2)/LAM2 | |
28 | ARG2 = MU2/LAM2 | |
29 | S = LOG ( LOG( ARG1 ) / LOG( ARG2 )) | |
30 | DS = SQRT (S) | |
31 | S2 = S * S | |
32 | S3 = S2 * S | |
33 | *...UV : | |
34 | NU = 2.484 + 0.116 * S + 0.093 * S2 | |
35 | AKU = 0.563 - 0.025 * S | |
36 | BKU = 0.054 + 0.154 * S | |
37 | AU = -0.326 - 0.058 * S - 0.135 * S2 | |
38 | BU = -3.322 + 8.259 * S - 3.119 * S2 + 0.291 * S3 | |
39 | CU = 11.52 - 12.99 * S + 3.161 * S2 | |
40 | DU = 2.808 + 1.400 * S - 0.557 * S2 + 0.119 * S3 | |
41 | UV = GR94FV (X, NU, AKU, BKU, AU, BU, CU, DU) | |
42 | ZUV = UV | |
43 | *...DV : | |
44 | ND = 0.156 - 0.017 * S | |
45 | AKD = 0.299 - 0.022 * S | |
46 | BKD = 0.259 - 0.015 * S | |
47 | AD = 3.445 + 1.278 * S + 0.326 * S2 | |
48 | BD = -6.934 + 37.45 * S - 18.95 * S2 + 1.463 * S3 | |
49 | CD = 55.45 - 69.92 * S + 20.78 * S2 | |
50 | DD = 3.577 + 1.441 * S - 0.683 * S2 + 0.179 * S3 | |
51 | DV = GR94FV (X, ND, AKD, BKD, AD, BD, CD, DD) | |
52 | ZDV = DV | |
53 | *...DEL : | |
54 | NE = 0.099 + 0.019 * S + 0.002 * S2 | |
55 | AKE = 0.419 - 0.013 * S | |
56 | BKE = 1.064 - 0.038 * S | |
57 | AE = -44.00 + 98.70 * S - 14.79 * S2 | |
58 | BE = 28.59 - 40.94 * S - 13.66 * S2 + 2.523 * S3 | |
59 | CE = 84.57 - 108.8 * S + 31.52 * S2 | |
60 | DE = 7.469 + 2.480 * S - 0.866 * S2 | |
61 | DEL = GR94FV (X, NE, AKE, BKE, AE, BE, CE, DE) | |
62 | ZDEL= DEL | |
63 | *...UDB : | |
64 | ALX = 1.215 | |
65 | BEX = 0.466 | |
66 | AKX = 0.326 + 0.150 * S | |
67 | BKX = 0.956 + 0.405 * S | |
68 | AGX = 0.272 | |
69 | BGX = 3.794 - 2.359 * DS | |
70 | CX = 2.014 | |
71 | DX = 7.941 + 0.534 * DS - 0.940 * S + 0.410 * S2 | |
72 | EX = 3.049 + 1.597 * S | |
73 | ESX = 4.396 - 4.594 * DS + 3.268 * S | |
74 | UDB = GR94FW(X, S, ALX, BEX, AKX, BKX, AGX, BGX, CX, DX, EX, ESX) | |
75 | ZUDB= UDB | |
76 | *...UB | |
77 | ZUB = (ZUDB - ZDEL)/2.D0 | |
78 | *...DB | |
79 | ZDB = (ZUDB + ZDEL)/2.D0 | |
80 | *...SB : | |
81 | ALS = 0.175 | |
82 | BES = 0.344 | |
83 | AKS = 1.415 - 0.641 * DS | |
84 | AS = 0.580 - 9.763 * DS + 6.795 * S - 0.558 * S2 | |
85 | BS = 5.617 + 5.709 * DS - 3.972 * S | |
86 | DST = 13.78 - 9.581 * S + 5.370 * S2 - 0.996 * S3 | |
87 | EST = 4.546 + 0.372 * S2 | |
88 | ESS = 5.053 - 1.070 * S + 0.805 * S2 | |
89 | SB = GR94FWS (X, S, ALS, BES, AKS, AS, BS, DST, EST, ESS) | |
90 | ZSB = SB | |
91 | *...GL : | |
92 | ALG = 1.258 | |
93 | BEG = 1.846 | |
94 | AKG = 2.423 | |
95 | BKG = 2.427 + 1.311 * S - 0.153 * S2 | |
96 | AG = 25.09 - 7.935 * S | |
97 | BG = -14.84 - 124.3 * DS + 72.18 * S | |
98 | CG = 590.3 - 173.8 * S | |
99 | DG = 5.196 + 1.857 * S | |
100 | EG = -1.648 + 3.988 * S - 0.432 * S2 | |
101 | ESG = 3.232 - 0.542 * S | |
102 | GL = GR94FW (X, S, ALG, BEG, AKG, BKG, AG, BG, CG, DG, EG, ESG) | |
103 | ZGL = GL | |
104 | *... | |
105 | ZCB = 0.D0 | |
106 | ZBB = 0.D0 | |
107 | ZTB = 0.D0 | |
108 | *... | |
109 | RETURN | |
110 | END |