]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:29:05 cernlib | |
6 | * Version 7.04 | |
7 | * | |
8 | * Revision 1.2 1996/10/28 15:04:20 pdflib | |
9 | * PDFLIB new version 7.04 | |
10 | * | |
11 | * Revision 1.1.1.1 1996/04/12 15:29:26 plothow | |
12 | * Version 7.01 | |
13 | * | |
14 | * | |
15 | #include "pdf/pilot.h" | |
16 | FUNCTION PDZXMT (IST, LP, XX, QQ, IRT) | |
17 | C | |
18 | #include "pdf/w50511.inc" | |
19 | #include "pdf/expdp.inc" | |
20 | + AIEX, FXIEX, PDFD, ONED | |
21 | REAL PDF | |
22 | C | |
23 | PARAMETER (D0=0.0, D1=1.0, D2=2.0, D3=3.0, D4=4.0, D10=10.0) | |
24 | PARAMETER (NEX = 3, MXFL = 6, NPN = 2, NST = 11) | |
25 | C | |
26 | DIMENSION | |
27 | 1 AC(0:NEX, 0:NPN, -MXFL:2, NST), A(0:NEX), T(0:NPN), FX(0:NEX), | |
28 | 1 ALM(NST), Q02(NST), MEX(NST), MPN(NST), MQRK(NST), IORD(NST) | |
29 | C | |
30 | DATA ONED/1.D0/ | |
31 | C | |
32 | DATA MEX, MPN, MQRK / NST*3, NST*2, NST*6 / | |
33 | C Set S1: PDF parameters fro | |
34 | DATA IORD(1)/2/, ALM(1)/0.212/, Q02(1)/4.00/ | |
35 | DATA (((AC(IEX,IPN,IFL,1), IEX=0,3), IPN=0,2), IFL=2,-6,-1) | |
36 | > / 1.34, 0.15, | |
37 | > 5.30, -1.96, -0.57, 0.16, 0.43, 1.08, -0.08, -0.02, | |
38 | > 0.06, -0.03, 1.62, 0.11, 3.68, -1.94, -0.33, 0.14, | |
39 | > 0.53, 0.87, -0.10, -0.01, 0.03, 0.02, 1.88, -0.33, | |
40 | > 7.52, -1.34, -2.78, 0.10, -1.13, 2.92, 0.13, -0.04, | |
41 | > 0.04, -0.49, -0.99, -0.33, 8.53, -1.55, -1.54, 0.03, | |
42 | > -1.08, 2.02, 0.10, -0.03, 0.39, -0.39, -0.99, -0.33, | |
43 | > 8.53, -1.55, -1.54, 0.03, -1.08, 2.02, 0.10, -0.03, | |
44 | > 0.39, -0.39, -0.99, -0.33, 8.53, -1.55, -1.54, 0.03, | |
45 | > -1.08, 2.02, 0.10, -0.03, 0.39, -0.39, -3.98, -0.15, | |
46 | > 7.46, 0.35, 0.72, -0.06, 0.96, 0.89, -0.63, 0.00, | |
47 | > -0.30, -0.04, -6.28, -0.18, 6.56, 0.65, 2.62, 0.02, | |
48 | > 1.40, 1.13, -1.18, -0.03, -0.38, -0.16,-13.08, -0.40, | |
49 | > 15.35, -0.43, 8.54, 0.31,-11.83, 3.18, -2.70, -0.12, | |
50 | > 4.16, -0.82 / | |
51 | C Set B1: PDF parameters fr | |
52 | DATA IORD(2)/2/, ALM(2)/0.194/, Q02(2)/4.00/ | |
53 | DATA (((AC(IEX,IPN,IFL,2), IEX=0,3), IPN=0,2), IFL=2,-6,-1) | |
54 | > / 1.30, 0.19, | |
55 | > 5.24, -1.81, -0.57, 0.15, 0.44, 1.06, -0.09, -0.02, | |
56 | > 0.05, -0.02, 1.59, 0.14, 3.65, -1.81, -0.34, 0.13, | |
57 | > 0.53, 0.86, -0.10, -0.01, 0.03, 0.02, 1.48, -0.14, | |
58 | > 6.75, -0.50, -2.49, -0.11, -0.54, 2.13, 0.04, 0.03, | |
59 | > -0.15, -0.24, -1.08, -0.13, 8.40, -0.88, -1.33, -0.21, | |
60 | > -0.51, 1.18, -0.03, 0.06, 0.07, -0.05, -1.08, -0.13, | |
61 | > 8.39, -0.88, -1.33, -0.21, -0.50, 1.18, -0.03, 0.06, | |
62 | > 0.07, -0.05, -1.08, -0.13, 8.39, -0.88, -1.33, -0.21, | |
63 | > -0.50, 1.18, -0.03, 0.06, 0.07, -0.05, -4.22, -0.02, | |
64 | > 7.29, 0.90, 0.88, -0.17, 1.08, 0.50, -0.69, 0.03, | |
65 | > -0.39, 0.08, -6.42, -0.09, 6.47, 1.03, 2.67, -0.03, | |
66 | > 1.39, 1.00, -1.21, -0.02, -0.42, -0.14,-12.92, -0.36, | |
67 | > 15.74, -0.30, 8.33, 0.32,-12.73, 3.35, -2.68, -0.13, | |
68 | > 4.51, -0.91 / | |
69 | C Set B2: PDF parameters fr | |
70 | DATA IORD(3)/2/, ALM(3)/0.191/, Q02(3)/4.00/ | |
71 | DATA (((AC(IEX,IPN,IFL,3), IEX=0,3), IPN=0,2), IFL=2,-6,-1) | |
72 | > / 1.38, 0.18, | |
73 | > 5.40, -1.91, -0.59, 0.16, 0.42, 1.11, -0.08, -0.02, | |
74 | > 0.06, -0.03, 1.64, 0.09, 3.74, -2.02, -0.33, 0.14, | |
75 | > 0.54, 0.88, -0.10, -0.01, 0.03, 0.02, 1.52, -0.72, | |
76 | > 7.75, -2.18, -2.71, 0.45, -1.56, 3.75, 0.15, -0.15, | |
77 | > 0.16, -0.76, -0.85, -0.82, 9.19, -2.76, -1.43, 0.35, | |
78 | > -0.92, 2.56, -0.03, -0.09, 0.12, -0.40, -0.85, -0.82, | |
79 | > 9.19, -2.76, -1.43, 0.35, -0.92, 2.56, -0.03, -0.10, | |
80 | > 0.12, -0.40, -0.85, -0.82, 9.19, -2.76, -1.43, 0.35, | |
81 | > -0.92, 2.56, -0.03, -0.10, 0.12, -0.40, -3.74, -0.58, | |
82 | > 9.63, -1.09, 0.21, 0.24, -1.13, 2.10, -0.50, -0.07, | |
83 | > 0.25, -0.33, -6.07, -0.52, 8.33, -0.52, 2.33, 0.22, | |
84 | > 0.28, 1.91, -1.15, -0.07, -0.28, -0.31,-12.08, -0.73, | |
85 | > 21.14, -1.92, 7.31, 0.54,-19.17, 4.59, -2.35, -0.18, | |
86 | > 6.64, -1.25 / | |
87 | C Set E1: PDF parameters fr | |
88 | DATA IORD(4)/2/, ALM(4)/0.155/, Q02(4)/4.00/ | |
89 | DATA (((AC(IEX,IPN,IFL,4), IEX=0,3), IPN=0,2), IFL=2,-6,-1) | |
90 | > / 1.43, 0.16, | |
91 | > 6.17, -1.94, -0.65, 0.16, 0.43, 1.12, -0.08, -0.02, | |
92 | > 0.06, -0.02, 1.69, 0.11, 3.69, -1.99, -0.33, 0.14, | |
93 | > 0.54, 0.90, -0.11, -0.01, 0.03, 0.02, 2.11, -0.33, | |
94 | > 7.93, -1.51, -3.01, 0.10, -1.40, 3.14, 0.18, -0.04, | |
95 | > 0.09, -0.55, -0.84, -0.32, 8.96, -1.70, -1.65, 0.02, | |
96 | > -1.24, 2.15, 0.12, -0.03, 0.45, -0.43, -0.84, -0.32, | |
97 | > 8.96, -1.70, -1.65, 0.02, -1.24, 2.15, 0.12, -0.03, | |
98 | > 0.45, -0.43, -0.84, -0.32, 8.96, -1.70, -1.65, 0.02, | |
99 | > -1.24, 2.15, 0.12, -0.03, 0.45, -0.43, -3.87, -0.15, | |
100 | > 7.83, 0.21, 0.85, -0.07, 1.00, 0.93, -0.73, 0.00, | |
101 | > -0.36, -0.03, -6.09, -0.17, 6.75, 0.54, 2.81, 0.01, | |
102 | > 1.74, 1.15, -1.34, -0.03, -0.56, -0.16,-12.56, -0.38, | |
103 | > 14.62, -0.41, 8.69, 0.30,-11.27, 3.19, -2.93, -0.12, | |
104 | > 4.29, -0.87 / | |
105 | C 1/2 strange sea; PDF parameters fro | |
106 | DATA IORD(5)/2/, ALM(5)/0.237/, Q02(5)/4.00/ | |
107 | DATA (((AC(IEX,IPN,IFL,5), IEX=0,3), IPN=0,2), IFL=2,-6,-1) | |
108 | > / 1.42, 0.16, | |
109 | > 5.40, -1.99, -0.59, 0.17, 0.41, 1.12, -0.08, -0.02, | |
110 | > 0.06, -0.03, 1.68, 0.08, 3.75, -2.09, -0.33, 0.15, | |
111 | > 0.53, 0.89, -0.10, -0.01, 0.03, 0.02, 0.90, -0.17, | |
112 | > 5.27, -0.20, -1.86, -0.10, 0.43, 1.67, -0.09, 0.02, | |
113 | > -0.26, -0.14, -1.48, -0.13, 7.83, -0.38, -0.89, -0.19, | |
114 | > -0.06, 0.68, -0.12, 0.04, 0.01, 0.05, -1.48, -0.13, | |
115 | > 7.83, -0.38, -0.89, -0.19, -0.05, 0.68, -0.13, 0.04, | |
116 | > 0.00, 0.05, -2.26, -0.15, 7.47, -0.23, -0.90, -0.10, | |
117 | > -0.61, 1.22, -0.06, 0.01, 0.28, -0.16, -4.68, -0.06, | |
118 | > 5.55, 1.13, 0.92, -0.12, 1.16, 0.50, -0.62, 0.01, | |
119 | > -0.26, 0.03, -6.83, -0.12, 5.24, 1.19, 2.68, -0.01, | |
120 | > 1.14, 0.93, -1.13, -0.03, -0.24, -0.13,-14.41, -0.28, | |
121 | > 11.48, 0.65, 9.65, 0.15, -7.50, 1.99, -2.98, -0.06, | |
122 | > 2.54, -0.43 / | |
123 | C Set S1-MSbar: PDF parameters fro | |
124 | DATA IORD(6)/2/, ALM(6)/0.212/, Q02(6)/4.0/ | |
125 | DATA (((AC(IEX,IPN,IFL,6), IEX=0,3), IPN=0,2), IFL=2,-6,-1) | |
126 | > / 1.75, 0.11, | |
127 | > 6.20, -2.35, -1.02, 0.26, -0.41, 1.68, 0.05, -0.06, | |
128 | > 0.29, -0.24, 2.03, 0.06, 4.43, -2.35, -0.78, 0.24, | |
129 | > -0.18, 1.52, 0.03, -0.04, 0.22, -0.19, 1.09, -0.24, | |
130 | > 5.97, -0.64, -2.41, 0.08, -0.90, 2.71, -0.12, 0.02, | |
131 | > -0.35, -0.20, -0.14, -0.49, 10.24, -2.57, -1.98, 0.02, | |
132 | > -1.43, 2.32, 0.23, -0.02, 0.44, -0.47, -0.14, -0.49, | |
133 | > 10.24, -2.57, -1.98, 0.02, -1.44, 2.32, 0.23, -0.02, | |
134 | > 0.45, -0.47, -0.15, -0.49, 10.23, -2.57, -1.98, 0.02, | |
135 | > -1.44, 2.32, 0.23, -0.02, 0.45, -0.47, -2.36, -0.49, | |
136 | > 9.00, -1.74, -1.42, 0.44, -0.46, 3.93, 0.21, -0.22, | |
137 | > 0.29, -1.34, -2.19, -1.07, 11.30, -4.85, -3.86, 1.56, | |
138 | > -7.20, 10.51, 1.57, -0.73, 3.85, -4.36,-24.77, 7.52, | |
139 | > -99.51, 36.02,-23.00, 0.48,-16.45, 16.51, 34.44, -6.26, | |
140 | > 97.19,-40.40 / | |
141 | C Set B1-MSbar: PDF parameters fro | |
142 | DATA IORD(7)/2/, ALM(7)/0.194/, Q02(7)/4.0/ | |
143 | DATA (((AC(IEX,IPN,IFL,7), IEX=0,3), IPN=0,2), IFL=2,-6,-1) | |
144 | > / 1.66, 0.18, | |
145 | > 6.04, -2.09, -0.94, 0.18, -0.25, 1.42, 0.03, -0.03, | |
146 | > 0.23, -0.14, 2.00, 0.09, 4.40, -2.24, -0.81, 0.24, | |
147 | > -0.20, 1.53, 0.05, -0.05, 0.25, -0.23, 0.92, -0.07, | |
148 | > 5.79, -0.01, -2.28, -0.16, -0.68, 1.93, -0.07, 0.06, | |
149 | > -0.23, -0.11, -0.60, -0.13, 9.31, -1.18, -1.76, -0.27, | |
150 | > -0.94, 1.31, 0.13, 0.09, 0.21, -0.10, -0.60, -0.13, | |
151 | > 9.31, -1.18, -1.76, -0.27, -0.94, 1.31, 0.13, 0.09, | |
152 | > 0.21, -0.10, -0.60, -0.13, 9.31, -1.18, -1.76, -0.27, | |
153 | > -0.94, 1.31, 0.14, 0.09, 0.21, -0.10, -2.94, -0.19, | |
154 | > 7.94, -0.46, -1.12, 0.16, -0.05, 2.93, 0.15, -0.13, | |
155 | > 0.27, -1.05, -2.95, -0.62, 9.97, -3.00, -3.21, 0.99, | |
156 | > -6.33, 8.42, 1.38, -0.51, 3.71, -3.61, -3.88, -0.78, | |
157 | > 3.80, -2.37, -1.59, -0.07, 2.13, 0.48, -0.05, 0.40, | |
158 | > 0.96, 2.30 / | |
159 | C Set B2-MSbar: PDF parameters fro | |
160 | DATA IORD(8)/2/, ALM(8)/0.191/, Q02(8)/4.0/ | |
161 | DATA (((AC(IEX,IPN,IFL,8), IEX=0,3), IPN=0,2), IFL=2,-6,-1) | |
162 | > / 1.77, 0.13, | |
163 | > 6.28, -2.30, -0.98, 0.23, -0.34, 1.60, 0.03, -0.04, | |
164 | > 0.26, -0.18, 2.04, 0.03, 4.48, -2.47, -0.75, 0.26, | |
165 | > -0.15, 1.52, 0.02, -0.05, 0.21, -0.19, 0.74, -0.59, | |
166 | > 6.31, -1.37, -2.44, 0.42, -1.62, 3.56, 0.07, -0.15, | |
167 | > 0.18, -0.77, -0.43, -0.86, 10.16, -3.14, -1.96, 0.43, | |
168 | > -1.91, 3.14, 0.20, -0.14, 0.53, -0.68, -0.43, -0.86, | |
169 | > 10.16, -3.14, -1.96, 0.43, -1.91, 3.14, 0.20, -0.14, | |
170 | > 0.53, -0.68, -0.43, -0.86, 10.16, -3.14, -1.96, 0.43, | |
171 | > -1.91, 3.14, 0.20, -0.14, 0.53, -0.68, -3.07, -0.66, | |
172 | > 8.57, -1.68, -1.03, 0.45, -0.32, 3.48, 0.04, -0.17, | |
173 | > 0.17, -0.98, -4.44, -0.68, 6.90, -1.82, -0.13, 0.50, | |
174 | > 1.46, 4.11, -0.23, -0.18, -0.53, -1.16, -7.03, -1.13, | |
175 | > 8.56, -5.47, 1.10, 1.07, -2.33, 9.08, -0.41, -0.35, | |
176 | > 0.87, -2.66 / | |
177 | C Set E1-MSbar: PDF parameters fro | |
178 | DATA IORD(9)/2/, ALM(9)/0.155/, Q02(9)/4.0/ | |
179 | DATA (((AC(IEX,IPN,IFL,9), IEX=0,3), IPN=0,2), IFL=2,-6,-1) | |
180 | > / 1.79, 0.12, | |
181 | > 7.03, -2.29, -1.05, 0.24, -0.38, 1.63, 0.03, -0.04, | |
182 | > 0.27, -0.18, 2.12, 0.02, 4.46, -2.57, -0.85, 0.32, | |
183 | > -0.28, 1.82, 0.07, -0.08, 0.29, -0.33, 1.58, -0.28, | |
184 | > 6.84, -1.08, -2.68, 0.05, -0.93, 2.76, 0.01, 0.00, | |
185 | > -0.26, -0.32, -0.10, -0.43, 10.43, -2.49, -2.29, 0.09, | |
186 | > -2.14, 2.80, 0.35, -0.06, 0.73, -0.67, -0.10, -0.43, | |
187 | > 10.43, -2.49, -2.29, 0.09, -2.14, 2.80, 0.35, -0.06, | |
188 | > 0.73, -0.67, -0.11, -0.43, 10.43, -2.48, -2.29, 0.09, | |
189 | > -2.14, 2.80, 0.35, -0.06, 0.73, -0.67, -2.53, -0.35, | |
190 | > 8.67, -1.24, -1.16, 0.26, -0.10, 3.26, 0.12, -0.15, | |
191 | > 0.27, -1.06, -3.91, -0.44, 6.85, -1.56, -0.19, 0.38, | |
192 | > 2.15, 4.07, -0.24, -0.17, -0.74, -1.24, -6.57, -0.90, | |
193 | > 7.27, -5.07, 1.15, 0.95, -0.28, 9.02, -0.48, -0.33, | |
194 | > 0.28, -2.75 / | |
195 | C (MSbar) 1/2 strange sea; PDF parameters fro | |
196 | DATA IORD(10)/2/, ALM(10)/0.237/, Q02(10)/4.00/ | |
197 | DATA (((AC(IEX,IPN,IFL,10), IEX=0,3), IPN=0,2), IFL=2,-6,-1) | |
198 | > / 1.84, 0.12, | |
199 | > 6.34, -2.40, -0.97, 0.22, -0.34, 1.53, 0.03, -0.04, | |
200 | > 0.25, -0.16, 2.08, 0.02, 4.53, -2.51, -0.66, 0.19, | |
201 | > -0.04, 1.24, -0.02, -0.01, 0.15, -0.05, 0.31, -0.10, | |
202 | > 4.18, 0.34, -1.84, -0.10, 0.05, 1.64, -0.06, 0.01, | |
203 | > -0.12, -0.16, -1.13, -0.15, 8.43, -0.64, -1.26, -0.16, | |
204 | > -0.39, 1.01, -0.01, 0.03, 0.05, -0.06, -1.13, -0.15, | |
205 | > 8.43, -0.64, -1.26, -0.16, -0.39, 1.01, -0.01, 0.03, | |
206 | > 0.05, -0.06, -1.82, -0.18, 7.94, -0.56, -1.40, -0.06, | |
207 | > -0.82, 1.65, 0.09, -0.01, 0.30, -0.31, -3.69, -0.15, | |
208 | > 5.72, 0.26, -0.47, 0.04, 0.93, 1.85, -0.10, -0.05, | |
209 | > -0.11, -0.50, -5.06, -0.25, 4.42, -0.14, 0.39, 0.16, | |
210 | > 2.38, 2.72, -0.35, -0.08, -0.63, -0.75, -9.92, -0.38, | |
211 | > -1.27, -1.60, 4.60, 0.24, 9.17, 4.40, -1.53, -0.08, | |
212 | > -2.88, -1.08 / | |
213 | C Set B0: PDF parameters fr | |
214 | DATA IORD(11)/1/, ALM(11)/0.144/, Q02(11)/4.00/ | |
215 | DATA (((AC(IEX,IPN,IFL,11), IEX=0,3), IPN=0,2), IFL=2,-6,-1) | |
216 | > / 1.38, 0.16, | |
217 | > 5.40, -1.97, -0.62, 0.19, 0.59, 1.24, -0.10, -0.02, | |
218 | > 0.03, -0.05, 1.67, 0.08, 3.75, -2.09, -0.33, 0.17, | |
219 | > 0.70, 0.98, -0.13, -0.01, 0.00, 0.02, 1.52, -0.25, | |
220 | > 7.01, -0.79, -3.17, -0.01, -0.90, 2.90, 0.25, 0.00, | |
221 | > -0.08, -0.54, -0.81, -0.07, 9.19, -0.89, -1.13, -0.46, | |
222 | > 0.35, 0.33, -0.26, 0.16, -0.49, 0.40, -0.81, -0.07, | |
223 | > 9.19, -0.89, -1.13, -0.46, 0.35, 0.33, -0.26, 0.16, | |
224 | > -0.49, 0.40, -0.81, -0.07, 9.19, -0.89, -1.13, -0.46, | |
225 | > 0.35, 0.33, -0.26, 0.16, -0.49, 0.40, -3.62, -0.06, | |
226 | > 8.30, 0.16, 0.03, -0.21, -0.60, 1.26, -0.48, 0.05, | |
227 | > 0.25, -0.15, -6.16, -0.11, 6.49, 0.71, 2.37, -0.05, | |
228 | > 1.28, 1.37, -1.24, -0.02, -0.41, -0.26,-12.68, -0.35, | |
229 | > 14.87, -0.17, 8.36, 0.28,-12.56, 3.39, -2.89, -0.12, | |
230 | > 4.75, -0.96 / | |
231 | ||
232 | IRT = 0 | |
233 | IFL = LP | |
234 | C | |
235 | PDZXMT = 0. | |
236 | IBLP = ABS(LP) | |
237 | IF(IBLP .GE. 4 .AND. QQ .LE. 1.5) RETURN | |
238 | IF(IBLP .GE. 5 .AND. QQ .LE. 5.0) RETURN | |
239 | C | |
240 | TM = TMAS | |
241 | X = XX | |
242 | Q0 = SQRT (Q02(IST)) | |
243 | ALAM = ALM(IST) | |
244 | SQ = LOG ( LOG( MAX(Q0,QQ) /ALAM) / LOG(Q0/ALAM) ) | |
245 | IF(IBLP .EQ. 6) THEN | |
246 | SQ = LOG ( LOG( MAX(Q0,QQ) /ALAM * (100./TM)) / LOG(Q0/ALAM)) | |
247 | ENDIF | |
248 | C | |
249 | FX(0) = EXP(D1) | |
250 | FX(1) = X | |
251 | FX(2) = 1.- X | |
252 | FX(3) = LOG (1.+ 1./X) | |
253 | ||
254 | PDFD = ONED | |
255 | DO 20 IEX = 0, MEX(IST) | |
256 | A(IEX) = AC(IEX, 0, IFL, IST) | |
257 | SQN=SQ | |
258 | DO 21 IPN = 1, MPN(IST) | |
259 | A(IEX) = A(IEX) + AC(IEX, IPN, IFL, IST) * SQN | |
260 | SQN = SQN*SQ | |
261 | 21 CONTINUE | |
262 | FXIEX = FX(IEX) | |
263 | AIEX = A(IEX) | |
264 | PDFD = PDFD * FXIEX **(AIEX) | |
265 | 20 CONTINUE | |
266 | C | |
267 | C PDZXMT = PDF / X | |
268 | PDF = PDFD | |
269 | IF(PDF.GT.0.) PDZXMT = PDF | |
270 | C | |
271 | RETURN | |
272 | C | |
273 | ENTRY VLAMBD (ISET, IORDER) | |
274 | C | |
275 | IORDER = IORD (ISET) | |
276 | VLAMBD = ALM (ISET) | |
277 | C | |
278 | RETURN | |
279 | C ************************* | |
280 | END |