]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PDF/spdf/pdfset.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Extracting PHOS and EMCAL trackers from the correspondig reconstructors (Yu.Belikov)
[u/mrichter/AliRoot.git] / PDF / spdf / pdfset.F
CommitLineData
21886bb6 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.4 2000/04/19 13:28:34 mclareni
6* Import version 8.04
7*
8* Revision 1.6 2000/04/15 14:10:46 plothow
9* Version 8.03: Nuclear PDFs introduced
10*
11* Revision 1.5 2000/04/09 14:47:29 plothow
12* version 8.02 - add CTEQ4 ref2 PDFs
13*
14* Revision 1.4 2000/04/09 10:06:25 plothow
15* new version 8.00 - MRS 98 and 99, CTEQ5 and GRV98 added
16*
17* Revision 1.3 1997/04/08 15:53:25 pdflib
18* Lambda_QCD for CTEQ4 sets corrected
19*
20* Revision 1.2 1996/10/28 15:03:30 pdflib
21* PDFLIB new version 7.04
22*
23* Revision 1.1.1.1 1996/04/12 15:29:25 plothow
24* Version 7.01
25*
26*
27#include "pdf/pilot.h"
28 SUBROUTINE PDFSET(PARM,VAL)
29C
30C **********************************************************************
31C * *
32C * Routine to set user parameters in the internal COMMON/W50511/, *
33C * COMMON/W50511a/ *
34C * COMMON/W50512/ *
35C * and COMMON/W50513/ *
36C * and provide filling of the internal COMMON/W50518/ *
37C * COMMON/W50518a/ *
38C * COMMON/W505110/ *
39C * COMMON/W505120/ *
40C * COMMON/W505121/ *
41C * and COMMON/W505122/ *
42C * *
43C * Author: H. Plothow-Besch *
44C * CERN-PPE, CH - 1211 Geneva 23, Switzerland *
45C * *
46C **********************************************************************
47C
48#include "pdf/impdp.inc"
49C
50#include "pdf/w5051p0.inc"
51#include "pdf/w5051p1.inc"
52#include "pdf/w5051p2.inc"
53#include "pdf/w5051p3.inc"
54#include "pdf/w5051p7.inc"
55#include "pdf/w50510.inc"
56#include "pdf/w50511.inc"
57#include "pdf/w50511a.inc"
58#include "pdf/w50512.inc"
59#include "pdf/w50513.inc"
60#include "pdf/w50514.inc"
61#include "pdf/w50514w.inc"
62#include "pdf/w50515.inc"
63#include "pdf/w50516.inc"
64#include "pdf/w50517.inc"
65#include "pdf/w50518.inc"
66#include "pdf/w50518a.inc"
67#include "pdf/w50519.inc"
68#include "pdf/w505110.inc"
69#include "pdf/w505120.inc"
70#include "pdf/w505121.inc"
71#include "pdf/w505122.inc"
72#include "pdf/expdp.inc"
73 + PQCDL4(NPTYMX,NGRMAX,NSETMX), PQCDL5(NPTYMX,NGRMAX,NSETMX),
74 + PXMIN(NPTYMX,NGRMAX,NSETMX), PXMAX(NPTYMX,NGRMAX,NSETMX),
75 + PQ2MIN(NPTYMX,NGRMAX,NSETMX), PQ2MAX(NPTYMX,NGRMAX,NSETMX)
76C.
77 PARAMETER (NSDIM = 16)
78 CHARACTER*2 AV(NSDIM)
79 CHARACTER*20 PARM(NCHDIM),STRING
80#include "pdf/expdp.inc"
81 + VAL(NCHDIM)
82C.
83 EQUIVALENCE (VALUE,IVALUE)
84C.
85 DATA (PQCDL4(1,1,K),K=1,ND11) /
86 + 0.0000D0, 0.0450D0, 0.5000D0, 0.4700D0, 0.4000D0,
87 + 0.2000D0, 0.4000D0, 0.2000D0, 0.2900D0, 0.1770D0 /
88 DATA (PQCDL4(1,2,K),K=1,ND12) /
89 + 0.2000D0, 0.2000D0, 0.2000D0, 0.2000D0, 0.2000D0,
90 + 0.3000D0, 0.1600D0, 0.2600D0, 0.3600D0 /
91 DATA (PQCDL4(1,3,K),K=1,ND13) /
92 + 0.1070D0, 0.2500D0, 0.1780D0, 0.0910D0, 0.2280D0,
93 + 0.0910D0, 0.2280D0, 0.1000D0, 0.1900D0, 0.1000D0,
94 + 0.1900D0, 0.1000D0, 0.1000D0, 0.1000D0, 0.1900D0,
95 + 0.1000D0, 0.3000D0, 0.1900D0, 0.1900D0, 0.1900D0,
96 + 0.1900D0, 0.1350D0, 0.1600D0, 0.2000D0, 0.2350D0,
97 + 0.2150D0, 0.2150D0, 0.2150D0, 0.2300D0, 0.2300D0,
98 + 0.2300D0, 0.2300D0, 0.2300D0, 0.2300D0, 0.2300D0,
99 + 0.2300D0, 0.2300D0, 0.2300D0, 0.2310D0, 0.2310D0,
100 + 0.2540D0, 0.2540D0, 0.2300D0, 0.2300D0, 0.1500D0,
101 + 0.2010D0, 0.2660D0, 0.3440D0, 0.4350D0, 0.5420D0,
102 + 0.3440D0, 0.5070D0, 0.2410D0, 0.3440D0, 0.2410D0,
103 + 0.3440D0, 0.2410D0, 0.3440D0, 0.3000D0, 0.3000D0,
104 + 0.3000D0, 0.3000D0, 0.3000D0, 0.3000D0, 0.2290D0,
105 + 0.3830D0, 0.3000D0, 0.3000D0, 0.3000D0, 0.2290D0,
106 + 0.3830D0, 0.1740D0, 0.1740D0, 0.1740D0, 0.1360D0,
107 + 0.2160D0, 0.2925D0, 0.3000D0, 0.3000D0, 0.3000D0,
108 + 0.2290D0, 0.3830D0, 0.3033D0, 0.2903D0, 0.3000D0,
109 + 0.3000D0, 0.3000D0, 0.3000D0, 0.3000D0, 0.3000D0,
110 + 0.3000D0, 0.2290D0, 0.3830D0, 0.3033D0, 0.2903D0,
111 + 0.3000D0, 0.3000D0, 0.3000D0, 0.3000D0, 0.3000D0 /
112 DATA (PQCDL4(1,4,K),K=1,ND14) /
113 + 0.2120D0, 0.1940D0, 0.1910D0, 0.1550D0, 0.2370D0,
114 + 0.2120D0, 0.1940D0, 0.1910D0, 0.1550D0, 0.2370D0,
115 + 0.1440D0, 0.1680D0, 0.2310D0, 0.2310D0, 0.3220D0,
116 + 0.2470D0, 0.1900D0, 0.2130D0, 0.2080D0, 0.2080D0,
117 + 0.3220D0, 0.2350D0, 0.1900D0, 0.2130D0, 0.2080D0,
118 + 0.2080D0, 0.3220D0, 0.2350D0, 0.1770D0, 0.2390D0,
119 + 0.2470D0, 0.2360D0, 0.2980D0, 0.2980D0, 0.2150D0,
120 + 0.2540D0, 0.2980D0, 0.3460D0, 0.4010D0, 0.3030D0,
121 + 0.2610D0, 0.2980D0, 0.2980D0, 0.3850D0, 0.2920D0,
122 + 0.1920D0, 0.3260D0, 0.3260D0, 0.3260D0, 0.3260D0,
123 + 0.3950D0, 0.3090D0, 0.3260D0, 0.3260D0, 0.1920D0,
124 + 0.3260D0 /
125 DATA (PQCDL4(1,5,K),K=1,ND15) /
126 + 0.1600D0, 0.2200D0, 0.2000D0, 0.2000D0, 0.2000D0,
127 + 0.2000D0, 0.2000D0, 0.2000D0, 0.2000D0, 0.2000D0,
128 + 0.2000D0, 0.2000D0, 0.2000D0, 0.2000D0 /
129 DATA (PQCDL4(1,6,K),K=1,ND16) /
130 + 0.2300D0 /
131 DATA (PQCDL4(1,7,K),K=1,ND17) /
132 + 0.2540D0,0.2540D0 /
133C.
134 DATA (PQCDL4(2,1,K),K=1,ND21) /
135 + 0.2000D0, 0.4000D0 /
136 DATA (PQCDL4(2,3,K),K=1,ND23) /
137 + 0.1900D0, 0.1900D0, 0.1900D0 /
138 DATA (PQCDL4(2,5,K),K=1,ND25) /
139 + 0.2000D0, 0.2000D0 /
140 DATA (PQCDL4(2,6,K),K=1,ND26) /
141 + 0.2310D0, 0.1810D0, 0.2810D0 /
142C.
143 DATA (PQCDL4(3,1,K),K=1,ND31) /
144 + 0.3800D0, 0.4400D0 /
145 DATA (PQCDL4(3,2,K),K=1,ND32) /
146 + 0.4000D0, 0.4000D0, 0.4000D0, 0.4000D0 /
147 DATA (PQCDL4(3,3,K),K=1,ND33) /
148 + 0.2000D0, 0.2000D0, 0.2000D0, 0.2000D0 /
149 DATA (PQCDL4(3,4,K),K=1,ND34) /
150 + 0.2000D0, 0.2000D0, 0.2000D0, 0.2000D0, 0.2000D0 /
151 DATA (PQCDL4(3,5,K),K=1,ND35) /
152 + 0.2000D0, 0.2000D0, 0.2000D0, 0.2000D0 /
153 DATA (PQCDL4(3,6,K),K=1,ND36) /
154 + 0.2000D0, 0.2000D0, 0.2000D0 /
155 DATA (PQCDL4(3,8,K),K=1,ND38) /
156 + 0.4000D0, 0.4000D0, 0.4000D0, 0.4000D0, 0.4000D0,
157 + 0.4000D0 /
158 DATA (PQCDL4(3,9,K),K=1,ND39) /
159 + 0.2000D0, 0.2000D0, 0.2000D0, 0.2000D0, 0.2000D0,
160 + 0.2000D0, 0.2000D0, 0.2000D0 /
161C.
162 DATA (PQCDL4(4,1,K),K=1,ND41) /
163 + 0.0000D0 /
164C.
165C.
166 DATA (PQCDL5(1,1,K),K=1,ND11) /
167 + 0.0000D0, 0.0300D0, 0.4111D0, 0.3844D0, 0.3226D0,
168 + 0.1519D0, 0.3226D0, 0.1519D0, 0.2274D0, 0.1330D0 /
169 DATA (PQCDL5(1,2,K),K=1,ND12) /
170 + 0.1519D0, 0.1519D0, 0.1519D0, 0.1519D0, 0.1519D0,
171 + 0.2048D0, 0.1016D0, 0.1745D0, 0.2513D0 /
172 DATA (PQCDL5(1,3,K),K=1,ND13) /
173 + 0.0651D0, 0.1671D0, 0.1144D0, 0.0544D0, 0.1507D0,
174 + 0.0544D0, 0.1507D0, 0.0604D0, 0.1230D0, 0.0604D0,
175 + 0.1230D0, 0.0604D0, 0.0604D0, 0.0604D0, 0.1230D0,
176 + 0.0604D0, 0.2048D0, 0.1230D0, 0.1230D0, 0.1230D0,
177 + 0.1230D0, 0.0842D0, 0.1016D0, 0.1303D0, 0.1559D0,
178 + 0.1412D0, 0.1412D0, 0.1412D0, 0.1522D0, 0.1522D0,
179 + 0.1522D0, 0.1522D0, 0.1522D0, 0.1522D0, 0.1522D0,
180 + 0.1522D0, 0.1522D0, 0.1522D0, 0.1529D0, 0.1529D0,
181 + 0.1700D0, 0.1700D0, 0.1522D0, 0.1522D0, 0.0946D0,
182 + 0.1310D0, 0.1790D0, 0.2388D0, 0.3110D0, 0.3988D0,
183 + 0.2388D0, 0.3698D0, 0.1604D0, 0.2388D0, 0.1604D0,
184 + 0.2388D0, 0.1604D0, 0.2388D0, 0.2048D0, 0.2048D0,
185 + 0.2048D0, 0.2048D0, 0.2048D0, 0.2048D0, 0.1515D0,
186 + 0.2694D0, 0.2048D0, 0.2048D0, 0.2048D0, 0.1515D0,
187 + 0.2694D0, 0.1305D0, 0.1305D0, 0.1305D0, 0.0999D0,
188 + 0.1651D0, 0.1991D0, 0.2048D0, 0.2048D0, 0.2048D0,
189 + 0.1515D0, 0.2694D0, 0.2074D0, 0.1974D0, 0.2048D0,
190 + 0.2048D0, 0.2048D0, 0.2048D0, 0.2048D0, 0.2048D0,
191 + 0.2048D0, 0.1515D0, 0.2694D0, 0.2074D0, 0.1974D0,
192 + 0.2048D0, 0.2048D0, 0.2048D0, 0.2048D0, 0.2048D0 /
193 DATA (PQCDL5(1,4,K),K=1,ND14) /
194 + 0.1390D0, 0.1259D0, 0.1238D0, 0.0981D0, 0.1574D0,
195 + 0.1390D0, 0.1259D0, 0.1238D0, 0.0981D0, 0.1574D0,
196 + 0.1063D0, 0.1250D0, 0.1520D0, 0.1520D0, 0.2200D0,
197 + 0.1640D0, 0.1430D0, 0.1390D0, 0.1350D0, 0.1350D0,
198 + 0.2200D0, 0.1550D0, 0.1430D0, 0.1390D0, 0.1350D0,
199 + 0.1350D0, 0.2200D0, 0.1550D0, 0.1320D0, 0.1580D0,
200 + 0.1640D0, 0.1810D0, 0.2020D0, 0.2020D0, 0.1400D0,
201 + 0.1690D0, 0.2020D0, 0.2390D0, 0.2820D0, 0.2060D0,
202 + 0.1740D0, 0.2020D0, 0.2020D0, 0.0000D0, 0.0000D0,
203 + 0.1460D0, 0.2260D0, 0.2260D0, 0.2260D0, 0.2260D0,
204 + 0.0000D0, 0.0000D0, 0.2260D0, 0.2260D0, 0.1460D0,
205 + 0.2260D0 /
206 DATA (PQCDL5(1,5,K),K=1,ND15) /
207 + 0.1016D0, 0.1684D0, 0.1303D0, 0.1519D0, 0.1530D0,
208 + 0.1310D0, 0.1310D0, 0.1530D0, 0.1530D0, 0.1310D0,
209 + 0.1310D0, 0.1310D0, 0.1530D0, 0.1530D0 /
210 DATA (PQCDL5(1,6,K),K=1,ND16) /
211 + 0.1559D0 /
212 DATA (PQCDL5(1,7,K),K=1,ND17) /
213 + 0.1700D0,0.1700D0 /
214C.
215 DATA (PQCDL5(2,1,K),K=1,ND21) /
216 + 0.1519D0, 0.3226D0 /
217 DATA (PQCDL5(2,3,K),K=1,ND23) /
218 + 0.1230D0, 0.1230D0, 0.1230D0 /
219 DATA (PQCDL5(2,5,K),K=1,ND25) /
220 + 0.1303D0, 0.1519D0 /
221 DATA (PQCDL5(2,6,K),K=1,ND26) /
222 + 0.1529D0, 0.1166D0, 0.1904D0 /
223C.
224 DATA (PQCDL5(3,1,K),K=1,ND31) /
225 + 0.3051D0, 0.3150D0 /
226 DATA (PQCDL5(3,2,K),K=1,ND32) /
227 + 0.3226D0, 0.3226D0, 0.3226D0, 0.3226D0 /
228 DATA (PQCDL5(3,3,K),K=1,ND33) /
229 + 0.1519D0, 0.1519D0, 0.1519D0, 0.1519D0 /
230 DATA (PQCDL5(3,4,K),K=1,ND34) /
231 + 0.1303D0, 0.1519D0, 0.1519D0, 0.1303D0, 0.1519D0 /
232 DATA (PQCDL5(3,5,K),K=1,ND35) /
233 + 0.1303D0, 0.1303D0, 0.1519D0, 0.1530D0 /
234 DATA (PQCDL5(3,6,K),K=1,ND36) /
235 + 0.1303D0, 0.1303D0, 0.1303D0 /
236 DATA (PQCDL5(3,8,K),K=1,ND38) /
237 + 0.3226D0, 0.3226D0, 0.3226D0, 0.3226D0, 0.3226D0,
238 + 0.3226D0 /
239 DATA (PQCDL5(3,9,K),K=1,ND39) /
240 + 0.3226D0, 0.3226D0, 0.3226D0, 0.3226D0, 0.3226D0,
241 + 0.3226D0, 0.3226D0, 0.3226D0 /
242C.
243 DATA (PQCDL5(4,1,K),K=1,ND41) /
244 + 0.0000D0 /
245C.
246C.
247 DATA (PXMIN(1,1,K),K=1,ND11) /
248 + 0.0000D0, 5.000D-3, 5.000D-3, 5.000D-3, 1.000D-3,
249 + 5.000D-3, 5.000D-3, 1.000D-4, 1.000D-4, 5.000D-5 /
250 DATA (PXMIN(1,2,K),K=1,ND12) /
251 + 5.000D-5, 5.000D-5, 5.000D-5, 5.000D-5, 5.000D-5,
252 + 5.000D-5, 5.000D-5, 5.000D-5, 5.000D-5 /
253 DATA (PXMIN(1,3,K),K=1,ND13) /
254 + 1.000D-4, 1.000D-4, 1.000D-4, 1.000D-4, 1.000D-4,
255 + 1.000D-4, 1.000D-4, 1.000D-4, 1.000D-4, 1.000D-4,
256 + 1.000D-4, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
257 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
258 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
259 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
260 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
261 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
262 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
263 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
264 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
265 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
266 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
267 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
268 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
269 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
270 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
271 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
272 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
273 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5 /
274 DATA (PXMIN(1,4,K),K=1,ND14) /
275 + 1.00D-10, 1.00D-10, 1.00D-10, 1.00D-10, 1.00D-10,
276 + 1.00D-10, 1.00D-10, 1.00D-10, 1.00D-10, 1.00D-10,
277 + 1.00D-10, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
278 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
279 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
280 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
281 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
282 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
283 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
284 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
285 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
286 + 1.000D-5 /
287 DATA (PXMIN(1,5,K),K=1,ND15) /
288 + 1.000D-6, 1.000D-6, 1.000D-6, 1.000D-6, 1.000D-6,
289 + 1.000D-6, 1.000D-6, 1.000D-4, 1.000D-4, 1.000D-4,
290 + 1.000D-4, 1.000D-9, 1.000D-9, 1.000D-9 /
291 DATA (PXMIN(1,6,K),K=1,ND16) /
292 + 1.000D-4 /
293 DATA (PXMIN(1,7,K),K=1,ND17) /
294 + 1.000D-6, 1.000D-6 /
295C.
296 DATA (PXMIN(2,1,K),K=1,ND21) /
297 + 5.000D-3, 5.000D-3 /
298 DATA (PXMIN(2,3,K),K=1,ND23) /
299 + 1.000D-5, 1.000D-5, 1.000D-5 /
300 DATA (PXMIN(2,5,K),K=1,ND25) /
301 + 1.000D-5, 1.000D-5 /
302 DATA (PXMIN(2,6,K),K=1,ND26) /
303 + 1.000D-3, 1.000D-3, 1.000D-3 /
304C.
305 DATA (PXMIN(3,1,K),K=1,ND31) /
306 + 1.000D-5, 1.000D-5 /
307 DATA (PXMIN(3,2,K),K=1,ND32) /
308 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5 /
309 DATA (PXMIN(3,3,K),K=1,ND33) /
310 + 1.000D-4, 1.000D-4, 1.000D-4, 1.000D-4 /
311 DATA (PXMIN(3,4,K),K=1,ND34) /
312 + 5.000D-4, 5.000D-4, 5.000D-4, 5.000D-4, 5.000D-4 /
313 DATA (PXMIN(3,5,K),K=1,ND35) /
314 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5 /
315 DATA (PXMIN(3,6,K),K=1,ND36) /
316 + 1.370D-3, 1.370D-3, 1.370D-3 /
317 DATA (PXMIN(3,8,K),K=1,ND38) /
318 + 1.000D-3, 1.000D-3, 1.000D-3, 1.000D-3, 1.000D-3,
319 + 1.000D-3 /
320 DATA (PXMIN(3,9,K),K=1,ND39) /
321 + 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5, 1.000D-5,
322 + 1.000D-5, 1.000D-5, 1.000D-5 /
323C.
324 DATA (PXMIN(4,1,K),K=1,ND41) /
325 + 0.0000D0 /
326C.
327C.
328 DATA (PXMAX(1,1,K),K=1,ND11) /
329 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99998D0, 0.99999D0,
330 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0 /
331 DATA (PXMAX(1,2,K),K=1,ND12) /
332 + 0.95000D0, 0.95000D0, 0.95000D0, 0.95000D0, 0.95000D0,
333 + 0.95000D0, 0.95000D0, 0.95000D0, 0.95000D0 /
334 DATA (PXMAX(1,3,K),K=1,ND13) /
335 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
336 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
337 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
338 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
339 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
340 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
341 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
342 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
343 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
344 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
345 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
346 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
347 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
348 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
349 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
350 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
351 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
352 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
353 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
354 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0 /
355 DATA (PXMAX(1,4,K),K=1,ND14) /
356 + 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0,
357 + 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0,
358 + 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0,
359 + 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0,
360 + 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0,
361 + 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0,
362 + 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0,
363 + 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0,
364 + 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0,
365 + 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0,
366 + 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0, 0.99998D0,
367 + 0.99998D0 /
368 DATA (PXMAX(1,5,K),K=1,ND15) /
369 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
370 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0,
371 + 0.99999D0, 0.99999D0, 0.99999D0, 0.99999D0 /
372 DATA (PXMAX(1,6,K),K=1,ND16) /
373 + 0.99999D0 /
374 DATA (PXMAX(1,7,K),K=1,ND17) /
375 + 0.99999D0, 0.99999D0 /
376C.
377 DATA (PXMAX(2,1,K),K=1,ND21) /
378 + 0.9998D0, 0.9998D0 /
379 DATA (PXMAX(2,3,K),K=1,ND23) /
380 + 0.9998D0, 0.9998D0, 0.9998D0 /
381 DATA (PXMAX(2,5,K),K=1,ND25) /
382 + 0.9998D0, 0.9998D0 /
383 DATA (PXMAX(2,6,K),K=1,ND26) /
384 + 0.9998D0, 0.9998D0, 0.9998D0 /
385C.
386 DATA (PXMAX(3,1,K),K=1,ND31) /
387 + 0.9000D0, 0.9000D0 /
388 DATA (PXMAX(3,2,K),K=1,ND32) /
389 + 0.9998D0, 0.9998D0, 0.9998D0, 0.9998D0 /
390 DATA (PXMAX(3,3,K),K=1,ND33) /
391 + 0.9998D0, 0.9998D0, 0.9998D0, 0.9998D0 /
392 DATA (PXMAX(3,4,K),K=1,ND34) /
393 + 0.9998D0, 0.9998D0, 0.9998D0, 0.9998D0, 0.9998D0 /
394 DATA (PXMAX(3,5,K),K=1,ND35) /
395 + 0.9998D0, 0.9998D0, 0.9998D0, 0.9998D0 /
396 DATA (PXMAX(3,6,K),K=1,ND36) /
397 + 0.9986D0, 0.9986D0, 0.9986D0 /
398 DATA (PXMAX(3,8,K),K=1,ND38) /
399 + 0.9998D0, 0.9998D0, 0.9998D0, 0.9998D0, 0.9998D0,
400 + 0.9998D0 /
401 DATA (PXMAX(3,9,K),K=1,ND39) /
402 + 0.9998D0, 0.9998D0, 0.9998D0, 0.9998D0, 0.9998D0,
403 + 0.9998D0, 0.9998D0, 0.9998D0 /
404C.
405 DATA (PXMAX(4,1,K),K=1,ND41) /
406 + 0.0000D0 /
407C.
408C.
409 DATA (PQ2MIN(1,1,K),K=1,ND11) /
410 + 0.000D0, 0.500D0, 1.800D0, 5.000D0, 4.000D0,
411 + 4.000D0, 4.000D0, 5.000D0, 5.000D0, 4.000D0 /
412 DATA (PQ2MIN(1,2,K),K=1,ND12) /
413 + 10.00D0, 10.00D0, 10.00D0, 10.00D0, 10.00D0,
414 + 10.00D0, 10.00D0, 10.00D0, 10.00D0 /
415 DATA (PQ2MIN(1,3,K),K=1,ND13) /
416 + 5.000D0, 5.000D0, 5.000D0, 5.000D0, 5.000D0,
417 + 5.000D0, 5.000D0, 5.000D0, 5.000D0, 5.000D0,
418 + 5.000D0, 5.000D0, 5.000D0, 5.000D0, 5.000D0,
419 + 5.000D0, 5.000D0, 5.000D0, 5.000D0, 5.000D0,
420 + 5.000D0, 5.000D0, 5.000D0, 5.000D0, 5.000D0,
421 + 5.000D0, 5.000D0, 5.000D0, 5.000D0, 5.000D0,
422 + 5.000D0, 5.000D0, 5.000D0, 5.000D0, 5.000D0,
423 + 5.000D0, 5.000D0, 5.000D0, 5.000D0, 5.000D0,
424 + 5.000D0, 5.000D0, 0.625D0, 0.625D0, 5.000D0,
425 + 5.000D0, 5.000D0, 5.000D0, 5.000D0, 5.000D0,
426 + 5.000D0, 5.000D0, 1.250D0, 1.250D0, 1.250D0,
427 + 1.250D0, 2.000D0, 2.000D0, 1.250D0, 1.250D0,
428 + 1.250D0, 1.250D0, 1.250D0, 1.250D0, 1.250D0,
429 + 1.250D0, 1.250D0, 1.250D0, 1.250D0, 1.250D0,
430 + 1.250D0, 1.250D0, 1.250D0, 1.250D0, 1.250D0,
431 + 1.250D0, 1.250D0, 1.250D0, 1.250D0, 1.250D0,
432 + 1.250D0, 1.250D0, 1.250D0, 1.250D0, 1.250D0,
433 + 1.250D0, 1.250D0, 1.250D0, 1.250D0, 1.250D0,
434 + 1.250D0, 1.250D0, 1.250D0, 1.250D0, 1.250D0,
435 + 1.250D0, 1.250D0, 1.250D0, 1.250D0, 1.250D0 /
436 DATA (PQ2MIN(1,4,K),K=1,ND14) /
437 + 4.000D0, 4.000D0, 4.000D0, 4.000D0, 4.000D0,
438 + 4.000D0, 4.000D0, 4.000D0, 4.000D0, 4.000D0,
439 + 4.000D0, 4.000D0, 4.000D0, 4.000D0, 4.000D0,
440 + 4.000D0, 4.000D0, 4.000D0, 4.000D0, 4.000D0,
441 + 4.000D0, 4.000D0, 4.000D0, 4.000D0, 4.000D0,
442 + 4.000D0, 4.000D0, 4.000D0, 4.000D0, 4.000D0,
443 + 4.000D0, 2.560D0, 2.560D0, 2.560D0, 2.560D0,
444 + 2.560D0, 2.560D0, 2.560D0, 2.560D0, 2.560D0,
445 + 0.490D0, 2.560D0, 1.000D0, 2.560D0, 2.560D0,
446 + 1.000D0, 1.000D0, 1.000D0, 1.000D0, 1.000D0,
447 + 1.000D0, 1.000D0, 1.000D0, 1.000D0, 1.000D0,
448 + 1.000D0 /
449 DATA (PQ2MIN(1,5,K),K=1,ND15) /
450 + 0.200D0, 0.200D0, 0.300D0, 0.250D0, 0.400D0,
451 + 0.400D0, 0.400D0, 0.400D0, 0.400D0, 0.400D0,
452 + 0.400D0, 0.800D0, 0.800D0, 0.800D0 /
453 DATA (PQ2MIN(1,6,K),K=1,ND16) /
454 + 2.0000D0 /
455 DATA (PQ2MIN(1,7,K),K=1,ND17) /
456 + 5.000D0, 5.000D0 /
457C.
458 DATA (PQ2MIN(2,1,K),K=1,ND21) /
459 + 4.0000D0, 4.0000D0 /
460 DATA (PQ2MIN(2,3,K),K=1,ND23) /
461 + 5.0000D0, 5.0000D0, 5.0000D0 /
462 DATA (PQ2MIN(2,5,K),K=1,ND25) /
463 + 0.3000D0, 0.2500D0 /
464 DATA (PQ2MIN(2,6,K),K=1,ND26) /
465 + 2.0000D0, 2.0000D0, 2.0000D0 /
466C.
467 DATA (PQ2MIN(3,1,K),K=1,ND31) /
468 + 10.000D0, 10.000D0 /
469 DATA (PQ2MIN(3,2,K),K=1,ND32) /
470 + 1.0000D0, 1.0000D0, 20.000D0, 200.00D0 /
471 DATA (PQ2MIN(3,3,K),K=1,ND33) /
472 + 4.0000D0, 4.0000D0, 1.0000D0, 4.0000D0 /
473 DATA (PQ2MIN(3,4,K),K=1,ND34) /
474 + 5.3000D0, 5.3000D0, 5.3000D0, 5.3000D0, 5.3000D0 /
475 DATA (PQ2MIN(3,5,K),K=1,ND35) /
476 + 0.3000D0, 0.3000D0, 0.2500D0, 0.6000D0 /
477 DATA (PQ2MIN(3,6,K),K=1,ND36) /
478 + 2.0000D0, 2.0000D0, 2.0000D0 /
479 DATA (PQ2MIN(3,8,K),K=1,ND38) /
480 + 4.0000D0, 4.0000D0, 4.0000D0, 4.0000D0, 4.0000D0,
481 + 4.0000D0 /
482 DATA (PQ2MIN(3,9,K),K=1,ND39) /
483 + 0.3600D0, 0.3600D0, 4.0000D0, 4.0000D0, 0.3600D0,
484 + 0.3600D0, 4.0000D0, 4.0000D0 /
485C.
486 DATA (PQ2MIN(4,1,K),K=1,ND41) /
487 + 0.0000D0 /
488C.
489C.
490 DATA (PQ2MAX(1,1,K),K=1,ND11) /
491 + 0.000D0, 1.000D6, 1.000D6, 1.000D6, 1.000D6,
492 + 1.000D6, 1.000D6, 1.000D8, 1.000D8, 4.000D7 /
493 DATA (PQ2MAX(1,2,K),K=1,ND12) /
494 + 1.000D8, 1.000D8, 1.000D8, 1.000D8, 1.000D8,
495 + 1.000D8, 1.000D8, 1.000D8, 1.000D8 /
496 DATA (PQ2MAX(1,3,K),K=1,ND13) /
497 + 1.310D6, 1.310D6, 1.310D6, 1.310D6, 1.310D6,
498 + 1.310D6, 1.310D6, 1.310D6, 1.310D6, 1.310D6,
499 + 1.310D6, 1.310D6, 1.310D6, 1.310D6, 1.310D6,
500 + 1.310D6, 1.310D6, 1.310D6, 1.310D6, 1.310D6,
501 + 2.680D9, 1.310D6, 1.310D6, 1.310D6, 1.310D6,
502 + 2.680D9, 2.680D9, 2.680D9, 2.680D9, 2.680D9,
503 + 2.680D9, 2.680D9, 2.680D9, 2.680D9, 2.680D9,
504 + 2.680D9, 1.310D6, 1.310D6, 1.310D6, 1.310D6,
505 + 1.310D6, 1.310D6, 1.310D6, 1.310D6, 1.310D6,
506 + 1.310D6, 1.310D6, 1.310D6, 1.310D6, 1.310D6,
507 + 1.310D6, 1.310D6, 1.000D7, 1.000D7, 1.000D7,
508 + 1.000D7, 1.000D7, 1.000D7, 1.000D7, 1.000D7,
509 + 1.000D7, 1.000D7, 1.000D7, 1.000D7, 1.000D7,
510 + 1.000D7, 1.000D7, 1.000D7, 1.000D7, 1.000D7,
511 + 1.000D7, 1.000D7, 1.000D7, 1.000D7, 1.000D7,
512 + 1.000D7, 1.000D7, 1.000D7, 1.000D7, 1.000D7,
513 + 1.000D7, 1.000D7, 1.000D7, 1.000D7, 1.000D7,
514 + 1.000D7, 1.000D7, 1.000D7, 1.000D7, 1.000D7,
515 + 1.000D7, 1.000D7, 1.000D7, 1.000D7, 1.000D7,
516 + 1.000D7, 1.000D7, 1.000D7, 1.000D7, 1.000D7 /
517 DATA (PQ2MAX(1,4,K),K=1,ND14) /
518 + 1.000D8, 1.000D8, 1.000D8, 1.000D8, 1.000D8,
519 + 1.000D8, 1.000D8, 1.000D8, 1.000D8, 1.000D8,
520 + 1.000D8, 1.000D6, 1.000D6, 1.000D6, 1.000D6,
521 + 1.000D6, 1.000D6, 1.000D6, 1.000D6, 1.000D6,
522 + 1.000D6, 1.000D6, 1.000D6, 1.000D6, 1.000D6,
523 + 1.000D6, 1.000D6, 1.000D6, 1.000D6, 1.000D6,
524 + 1.000D6, 1.000D8, 1.000D8, 1.000D8, 1.000D8,
525 + 1.000D6, 1.000D8, 1.000D8, 1.000D8, 1.000D8,
526 + 1.000D8, 1.000D8, 1.000D8, 1.000D8, 1.000D8,
527 + 1.000D8, 1.000D8, 1.000D8, 1.000D8, 1.000D8,
528 + 1.000D8, 1.000D8, 1.000D8, 1.000D8, 1.000D8,
529 + 1.000D8 /
530 DATA (PQ2MAX(1,5,K),K=1,ND15) /
531 + 1.000D8, 1.000D8, 1.000D8, 1.000D8, 1.000D6,
532 + 1.000D6, 1.000D6, 1.000D4, 1.000D4, 1.000D4,
533 + 1.000D4, 1.000D6, 1.000D6, 1.000D6 /
534 DATA (PQ2MAX(1,6,K),K=1,ND16) /
535 + 1.0000D8 /
536 DATA (PQ2MAX(1,7,K),K=1,ND17) /
537 + 1.310D6, 1.310D6 /
538C.
539 DATA (PQ2MAX(2,1,K),K=1,ND21) /
540 + 2.0000D3, 2.0000D3 /
541 DATA (PQ2MAX(2,3,K),K=1,ND23) /
542 + 1.3100D6, 1.3100D6, 1.3100D6 /
543 DATA (PQ2MAX(2,5,K),K=1,ND25) /
544 + 1.0000D6, 1.0000D6 /
545 DATA (PQ2MAX(2,6,K),K=1,ND26) /
546 + 1.0000D8, 1.0000D8, 1.0000D8 /
547C.
548 DATA (PQ2MAX(3,1,K),K=1,ND31) /
549 + 1.0000D4, 1.0000D4 /
550 DATA (PQ2MAX(3,2,K),K=1,ND32) /
551 + 1.0000D4, 50.000D0, 500.00D0, 1.0000D4 /
552 DATA (PQ2MAX(3,3,K),K=1,ND33) /
553 + 1.0000D5, 1.0000D5, 1.0000D5, 1.0000D5 /
554 DATA (PQ2MAX(3,4,K),K=1,ND34) /
555 + 1.0000D8, 1.0000D8, 1.0000D8, 1.0000D8, 1.0000D8 /
556 DATA (PQ2MAX(3,5,K),K=1,ND35) /
557 + 1.0000D6, 1.0000D6, 1.0000D6, 5.0000D4 /
558 DATA (PQ2MAX(3,6,K),K=1,ND36) /
559 + 5.5000D5, 5.5000D5, 5.5000D5/
560 DATA (PQ2MAX(3,8,K),K=1,ND38) /
561 + 2.5000D3, 2.5000D3, 2.5000D3, 2.5000D3, 2.5000D3,
562 + 2.5000D3 /
563 DATA (PQ2MAX(3,9,K),K=1,ND39) /
564 + 5.0000D4, 5.0000D4, 5.0000D4, 5.0000D4, 5.0000D4,
565 + 5.0000D4, 5.0000D4, 5.0000D4 /
566C.
567 DATA (PQ2MAX(4,1,K),K=1,ND41) /
568 + 0.0000D0 /
569C.
570C.
571 DATA AV /'BE','OR','GH','DO','EH','DF','MR','HM','KM',
572 + 'MT','CT','GR','AB','BM','WI','SA'/
573C.
574 DATA ZEROD/0.D0/,PONED/0.1D0/,ONED/1.D0/,ONEDO/1.1D0/,TWOD/2.D0/
575 DATA ZEROMIN/0.0001D0/
576 DATA ISTART/0/
577C.
578 SAVE ISTART, PQCDL4, PQCDL5, PXMIN, PXMAX, PQ2MIN, PQ2MAX,
579 + /W50511/, /W505110/, /W505120/, /W505121/, /W505122/,
580 + /W50514/, /W50514W/, /W50515/, /W50516/, /W50518/
581C.
582 IF (ISTART.EQ.0) IFLSET = 0
583 NAFLAG = 0
584C.
585C Set the PDFLIB parameters to the values of the previous call to PDFSET
586 IF(IFLSET.EQ.1) THEN
587 IFLPRT = IFLPRTP
588 NPTYPE = NPTYPP
589 NGROUP = NGROPP
590 NSET = NSETP
591 NFL = NFLP
592 LO = LOP
593 TMAS = TMASP
594 QCDL4 = QCDL4P
595 QCDL5 = QCDL5P
596 XMIN = XMINP
597 XMAX = XMAXP
598 Q2MIN = Q2MINP
599 Q2MAX = Q2MAXP
600 NATYPE = NATYPP
601 NAGROUP = NAGROPP
602 NASET = NASETP
603 ENDIF
604C.
605C Initialisation of PDFLIB parameters (first call to PDFSET)
606 IF(ISTART.EQ.0) THEN
607C.
608C set maximum dimensions for PDFs per particle type and per author group
609 CALL UFILL(NPGSMX,1,NPTYMX*NGRMAX,0)
610 NPGSMX(1,1) = ND11
611 NPGSMX(1,2) = ND12
612 NPGSMX(1,3) = ND13
613 NPGSMX(1,4) = ND14
614 NPGSMX(1,5) = ND15
615 NPGSMX(1,6) = ND16
616 NPGSMX(1,7) = ND17
617 NPGSMX(1,8) = ND18
618 NPGSMX(1,9) = ND19
619 NPGSMX(2,1) = ND21
620 NPGSMX(2,2) = ND22
621 NPGSMX(2,3) = ND23
622 NPGSMX(2,4) = ND24
623 NPGSMX(2,5) = ND25
624 NPGSMX(2,6) = ND26
625 NPGSMX(2,7) = ND27
626 NPGSMX(2,8) = ND28
627 NPGSMX(2,9) = ND29
628 NPGSMX(3,1) = ND31
629 NPGSMX(3,2) = ND32
630 NPGSMX(3,3) = ND33
631 NPGSMX(3,4) = ND34
632 NPGSMX(3,5) = ND35
633 NPGSMX(3,6) = ND36
634 NPGSMX(3,7) = ND37
635 NPGSMX(3,8) = ND38
636 NPGSMX(3,9) = ND39
637 NPGSMX(4,1) = ND41
638C.
639C set retracted PDF sets invalid
640 CALL UFILL(NSETFL,1,NPTYMX*NGRMAX*NSETMX,0)
641 DO KPTY = 1,NPTYMX
642 DO KGRO = 1,NGRMAX
643 IF(NPGSMX(KPTY,KGRO) .EQ. 0) GOTO 2
644 DO ISET = 1,NPGSMX(KPTY,KGRO)
645 NSETFL(KPTY,KGRO,ISET) = 1
646 IF (KPTY.EQ.1) THEN
647 IF (KGRO.EQ.3) THEN
648 IF (ISET.EQ.8) THEN
649 NSETFL(KPTY,KGRO,ISET) = 0
650 ENDIF
651 ELSEIF (KGRO.EQ.5) THEN
652 IF (ISET.EQ.1 .OR. ISET.EQ.2) THEN
653 NSETFL(KPTY,KGRO,ISET) = 0
654 ENDIF
655 ENDIF
656 ENDIF
657 ENDDO
658 2 CONTINUE
659 ENDDO
660 ENDDO
661C.
662C set character*8 NAME for each PDF
663 DO NPTY=1,NPTYMX
664 DO NGR=1,NGRMAX
665 DO NSET=1,NSETMX
666 SFNAME(NPTY,NGR,NSET) = ' '
667 ENDDO
668 ENDDO
669 ENDDO
670 SFNAME(1,1,1) = 'PRIVATE '
671 SFNAME(1,1,2) = 'BEBC '
672 SFNAME(1,1,3) = 'OR '
673 SFNAME(1,1,4) = 'BEP '
674 SFNAME(1,1,5) = 'GHR '
675 SFNAME(1,1,6) = 'DO-1 '
676 SFNAME(1,1,7) = 'DO-2 '
677 SFNAME(1,1,8) = 'EHLQ-1 '
678 SFNAME(1,1,9) = 'EHLQ-2 '
679 SFNAME(1,1,10) = 'DO-1.1 '
680 SFNAME(1,2,1) = 'DFLM-SV '
681 SFNAME(1,2,2) = 'DFLM-HV '
682 SFNAME(1,2,3) = 'DFLM-SG '
683 SFNAME(1,2,4) = 'DFLM-HG '
684 SFNAME(1,2,5) = 'DFLM-LO '
685 SFNAME(1,2,6) = 'DFLM-HO '
686 SFNAME(1,2,7) = 'DFLM-160'
687 SFNAME(1,2,8) = 'DFLM-260'
688 SFNAME(1,2,9) = 'DFLM-360'
689 SFNAME(1,3,1) = 'MRS-1 '
690 SFNAME(1,3,2) = 'MRS-2 '
691 SFNAME(1,3,3) = 'MRS-3 '
692 SFNAME(1,3,4) = 'MRS-E '
693 SFNAME(1,3,5) = 'MRS-B '
694 SFNAME(1,3,6) = 'MRS-EP '
695 SFNAME(1,3,7) = 'MRS-BP '
696 SFNAME(1,3,8) = 'HMRS-E1.'
697 SFNAME(1,3,9) = 'HMRS-B1.'
698 SFNAME(1,3,10) = 'HMRS-E3.'
699 SFNAME(1,3,11) = 'HMRS-B3.'
700 SFNAME(1,3,12) = 'HMRSE+ '
701 SFNAME(1,3,13) = 'HMRSE- '
702 SFNAME(1,3,14) = 'HMRS-E '
703 SFNAME(1,3,15) = 'HMRS-B '
704 SFNAME(1,3,16) = 'HMRSB100'
705 SFNAME(1,3,17) = 'HMRSB300'
706 SFNAME(1,3,18) = 'KMRS-BM2'
707 SFNAME(1,3,19) = 'KMRS-BM5'
708 SFNAME(1,3,20) = 'KMRS-BM '
709 SFNAME(1,3,21) = 'KMRS-B0 '
710 SFNAME(1,3,22) = 'MRS-B0-1'
711 SFNAME(1,3,23) = 'MRS-B0-2'
712 SFNAME(1,3,24) = 'MRS-B0-3'
713 SFNAME(1,3,25) = 'MRS-B0-4'
714 SFNAME(1,3,26) = 'MRS-S0 '
715 SFNAME(1,3,27) = 'MRS-D0 '
716 SFNAME(1,3,28) = 'MRS-DM '
717 SFNAME(1,3,29) = 'MRS-S0P '
718 SFNAME(1,3,30) = 'MRS-D0P '
719 SFNAME(1,3,31) = 'MRS-DMP '
720 SFNAME(1,3,32) = 'MRS-S0PD'
721 SFNAME(1,3,33) = 'MRS-D0PD'
722 SFNAME(1,3,34) = 'MRS-DMPD'
723 SFNAME(1,3,35) = 'MRS-(H) '
724 SFNAME(1,3,36) = 'MRS-(H)D'
725 SFNAME(1,3,37) = 'MRS-(A) '
726 SFNAME(1,3,38) = 'MRS-(A)F'
727 SFNAME(1,3,39) = 'MRS-Ap '
728 SFNAME(1,3,40) = 'MRS-Ap F'
729 SFNAME(1,3,41) = 'MRS-G '
730 SFNAME(1,3,42) = 'MRS-GFit'
731 SFNAME(1,3,43) = 'MRS-A lQ'
732 SFNAME(1,3,44) = 'MRS-ADlQ'
733 SFNAME(1,3,45) = 'MRS-Ap-1'
734 SFNAME(1,3,46) = 'MRS-Ap-2'
735 SFNAME(1,3,47) = 'MRS-Ap-3'
736 SFNAME(1,3,48) = 'MRS-Ap-4'
737 SFNAME(1,3,49) = 'MRS-Ap-5'
738 SFNAME(1,3,50) = 'MRS-Ap-6'
739 SFNAME(1,3,51) = 'MRS-J '
740 SFNAME(1,3,52) = 'MRS-Jp '
741 SFNAME(1,3,53) = 'MRS-R1 '
742 SFNAME(1,3,54) = 'MRS-R2 '
743 SFNAME(1,3,55) = 'MRS-R3 '
744 SFNAME(1,3,56) = 'MRS-R4 '
745 SFNAME(1,3,57) = 'MRS-R1Ft'
746 SFNAME(1,3,58) = 'MRS-R2Ft'
747 SFNAME(1,3,59) = 'MRRS-ch1'
748 SFNAME(1,3,60) = 'MRRS-ch2'
749 SFNAME(1,3,61) = 'MRRS-ch3'
750 SFNAME(1,3,62) = 'MRST-cgD'
751 SFNAME(1,3,63) = 'MRST-hgD'
752 SFNAME(1,3,64) = 'MRST-lgD'
753 SFNAME(1,3,65) = 'MRSTlasD'
754 SFNAME(1,3,66) = 'MRSThasD'
755 SFNAME(1,3,67) = 'MRST-cg '
756 SFNAME(1,3,68) = 'MRST-hg '
757 SFNAME(1,3,69) = 'MRST-lg '
758 SFNAME(1,3,70) = 'MRST-las'
759 SFNAME(1,3,71) = 'MRST-has'
760 SFNAME(1,3,72) = 'MRST-cgO'
761 SFNAME(1,3,73) = 'MRST-hgO'
762 SFNAME(1,3,74) = 'MRST-lgO'
763 SFNAME(1,3,75) = 'MRSTlasO'
764 SFNAME(1,3,76) = 'MRSThasO'
765 SFNAME(1,3,77) = 'MRST-ht '
766 SFNAME(1,3,78) = 'MRST-cgD'
767 SFNAME(1,3,79) = 'MRST-hgD'
768 SFNAME(1,3,80) = 'MRST-lgD'
769 SFNAME(1,3,81) = 'MRSTlasD'
770 SFNAME(1,3,82) = 'MRSThasD'
771 SFNAME(1,3,83) = 'MRSTqupD'
772 SFNAME(1,3,84) = 'MRSTqdnD'
773 SFNAME(1,3,85) = 'MRSTsupD'
774 SFNAME(1,3,86) = 'MRSTsdnD'
775 SFNAME(1,3,87) = 'MRSTcupD'
776 SFNAME(1,3,88) = 'MRSTcdnD'
777 SFNAME(1,3,89) = 'MRST-cg '
778 SFNAME(1,3,90) = 'MRST-hg '
779 SFNAME(1,3,91) = 'MRST-lg '
780 SFNAME(1,3,92) = 'MRST-las'
781 SFNAME(1,3,93) = 'MRST-has'
782 SFNAME(1,3,94) = 'MRST-qup'
783 SFNAME(1,3,95) = 'MRST-qdn'
784 SFNAME(1,3,96) = 'MRST-sup'
785 SFNAME(1,3,97) = 'MRST-sdn'
786 SFNAME(1,3,98) = 'MRST-cup'
787 SFNAME(1,3,99) = 'MRST-cdn'
788 SFNAME(1,3,100)= 'MRST-d/u'
789 SFNAME(1,4,1) = 'MT-S1DIS'
790 SFNAME(1,4,2) = 'MT-B1DIS'
791 SFNAME(1,4,3) = 'MT-B2DIS'
792 SFNAME(1,4,4) = 'MT-E1DIS'
793 SFNAME(1,4,5) = 'MT-6 DIS'
794 SFNAME(1,4,6) = 'MT-S1MSB'
795 SFNAME(1,4,7) = 'MT-B1MSB'
796 SFNAME(1,4,8) = 'MT-B2MSB'
797 SFNAME(1,4,9) = 'MT-E1MSB'
798 SFNAME(1,4,10) = 'MT-6 MSB'
799 SFNAME(1,4,11) = 'MT-LO '
800 SFNAME(1,4,12) = 'CTEQ-1LO'
801 SFNAME(1,4,13) = 'CTEQ-1M '
802 SFNAME(1,4,14) = 'CTEQ-1MS'
803 SFNAME(1,4,15) = 'CTEQ-1ML'
804 SFNAME(1,4,16) = 'CTEQ-1D '
805 SFNAME(1,4,17) = 'CTEQ-2LO'
806 SFNAME(1,4,18) = 'CTEQ-2M '
807 SFNAME(1,4,19) = 'CTEQ-2MS'
808 SFNAME(1,4,20) = 'CTEQ-2MF'
809 SFNAME(1,4,21) = 'CTEQ-2ML'
810 SFNAME(1,4,22) = 'CTEQ-2D '
811 SFNAME(1,4,23) = 'CTEQ-2pL'
812 SFNAME(1,4,24) = 'CTEQ-2pM'
813 SFNAME(1,4,25) = 'CTEQ2pMS'
814 SFNAME(1,4,26) = 'CTEQ2pMF'
815 SFNAME(1,4,27) = 'CTEQ2pML'
816 SFNAME(1,4,28) = 'CTEQ-2pD'
817 SFNAME(1,4,29) = 'CTEQ-3L '
818 SFNAME(1,4,30) = 'CTEQ-3M '
819 SFNAME(1,4,31) = 'CTEQ-3D '
820 SFNAME(1,4,32) = 'CTEQ-4L '
821 SFNAME(1,4,33) = 'CTEQ-4D '
822 SFNAME(1,4,34) = 'CTEQ-4M '
823 SFNAME(1,4,35) = 'CTEQ-4A1'
824 SFNAME(1,4,36) = 'CTEQ-4A2'
825 SFNAME(1,4,37) = 'CTEQ-4A3'
826 SFNAME(1,4,38) = 'CTEQ-4A4'
827 SFNAME(1,4,39) = 'CTEQ-4A5'
828 SFNAME(1,4,40) = 'CTEQ-4HJ'
829 SFNAME(1,4,41) = 'CTEQ-4lQ'
830 SFNAME(1,4,42) = 'CTEQ-4HQ'
831 SFNAME(1,4,43) = 'CTQ-4HQ1'
832 SFNAME(1,4,44) = 'CTEQ-4F3'
833 SFNAME(1,4,45) = 'CTEQ-4F4'
834 SFNAME(1,4,46) = 'CTEQ-5L '
835 SFNAME(1,4,47) = 'CTEQ-5D '
836 SFNAME(1,4,48) = 'CTEQ-5M '
837 SFNAME(1,4,49) = 'CTEQ-5HJ'
838 SFNAME(1,4,50) = 'CTEQ-5HQ'
839 SFNAME(1,4,51) = 'CTEQ-5F3'
840 SFNAME(1,4,52) = 'CTEQ-5F4'
841 SFNAME(1,4,53) = 'CTEQ-5M1'
842 SFNAME(1,4,54) = 'CTQ-5HQ1'
843 SFNAME(1,4,55) = 'CTEQ-5Lp'
844 SFNAME(1,4,56) = 'CTQ-5M1p'
845 SFNAME(1,5,1) = 'GRV-HO '
846 SFNAME(1,5,2) = 'GRV-LO '
847 SFNAME(1,5,3) = 'GRV-HO '
848 SFNAME(1,5,4) = 'GRV-LO '
849 SFNAME(1,5,5) = 'GRV94-LO'
850 SFNAME(1,5,6) = 'GRV94-MS'
851 SFNAME(1,5,7) = 'GRV94-DI'
852 SFNAME(1,5,8) = 'GRSV-LOS'
853 SFNAME(1,5,9) = 'GRSV-LOV'
854 SFNAME(1,5,10) = 'GRSV-HOS'
855 SFNAME(1,5,11) = 'GRSV-HOV'
856 SFNAME(1,5,12) = 'GRV98-LO'
857 SFNAME(1,5,13) = 'GRV98-MS'
858 SFNAME(1,5,14) = 'GRV98-D '
859 SFNAME(1,6,1) = 'ABFOW '
860 SFNAME(1,7,1) = 'BM-A '
861 SFNAME(1,7,2) = 'BM-B '
862 SFNAME(2,1,1) = 'OW-P-1 '
863 SFNAME(2,1,2) = 'OW-P-2 '
864 SFNAME(2,3,1) = 'SMRS-P-1'
865 SFNAME(2,3,2) = 'SMRS-P-2'
866 SFNAME(2,3,3) = 'SMRS-P-3'
867 SFNAME(2,5,1) = 'GRV-P-HO'
868 SFNAME(2,5,2) = 'GRV-P-LO'
869 SFNAME(2,6,1) = 'ABFKW-P1'
870 SFNAME(2,6,2) = 'ABFKW-P2'
871 SFNAME(2,6,3) = 'ABFKW-P3'
872 SFNAME(3,1,1) = 'DO-G-LO '
873 SFNAME(3,1,2) = 'DO-G-HO '
874 SFNAME(3,2,1) = 'DG-G-1 '
875 SFNAME(3,2,2) = 'DG-G-2 '
876 SFNAME(3,2,3) = 'DG-G-3 '
877 SFNAME(3,2,4) = 'DG-G-4 '
878 SFNAME(3,3,1) = 'LAC-G-1 '
879 SFNAME(3,3,2) = 'LAC-G-2 '
880 SFNAME(3,3,3) = 'LAC-G-3 '
881 SFNAME(3,3,4) = 'GAL-G '
882 SFNAME(3,4,1) = 'GS-G-1 '
883 SFNAME(3,4,2) = 'GS-G-2 '
884 SFNAME(3,4,3) = 'GS-G-3 '
885 SFNAME(3,4,4) = 'GS-96-HO'
886 SFNAME(3,4,5) = 'GS-96-LO'
887 SFNAME(3,5,1) = 'GRV-GLHO'
888 SFNAME(3,5,2) = 'GRV-G-HO'
889 SFNAME(3,5,3) = 'GRV-G-LO'
890 SFNAME(3,5,4) = 'GRS-G-LO'
891 SFNAME(3,6,1) = 'ACFGP-GH'
892 SFNAME(3,6,2) = 'ACFGP-GC'
893 SFNAME(3,6,3) = 'AFG-G-HO'
894 SFNAME(3,8,1) = 'WHIT-G-1'
895 SFNAME(3,8,2) = 'WHIT-G-2'
896 SFNAME(3,8,3) = 'WHIT-G-3'
897 SFNAME(3,8,4) = 'WHIT-G-4'
898 SFNAME(3,8,5) = 'WHIT-G-5'
899 SFNAME(3,8,6) = 'WHIT-G-6'
900 SFNAME(3,9,1) = 'SaS-1G1D'
901 SFNAME(3,9,2) = 'SaS-1G1M'
902 SFNAME(3,9,3) = 'SaS-1G2D'
903 SFNAME(3,9,4) = 'SaS-1G2M'
904 SFNAME(3,9,5) = 'SaS-2G1D'
905 SFNAME(3,9,6) = 'SaS-2G1M'
906 SFNAME(3,9,7) = 'SaS-2G2D'
907 SFNAME(3,9,8) = 'SaS-2G2M'
908 SFNAME(4,1,1) = 'A-EKS98 '
909C.
910C set cross reference tables for NEW PDFLIB format version 4.0
911C (give OLD to get NEW)
912 DO I=1,MODEMX
913 NPTYCR(I) = I/100 + 1
914 IF(I.LE.99) THEN
915 IF(I.LE. 9) THEN
916 NGROCR(I) = 1
917 NSETCR(I) = I
918 ELSEIF(I.GE. 10 .AND. I.LE. 19) THEN
919 NGROCR(I) = 2
920 NSETCR(I) = I - 9
921 ELSEIF(I.GE. 20 .AND. I.LE. 47) THEN
922 NGROCR(I) = 3
923 NSETCR(I) = I - 19
924 ELSEIF(I.GE. 48 .AND. I.LE. 49) THEN
925 NGROCR(I) = -1
926 NSETCR(I) = -1
927 ELSEIF(I.GE. 50 .AND. I.LE. 61) THEN
928 NGROCR(I) = 4
929 NSETCR(I) = I - 49
930 ELSEIF(I.GE. 62 .AND. I.LE. 69) THEN
931 NGROCR(I) = -1
932 NSETCR(I) = -1
933 ELSEIF(I.GE. 70 .AND. I.LE. 73) THEN
934 NGROCR(I) = 5
935 NSETCR(I) = I - 69
936 ELSEIF(I.EQ. 74 ) THEN
937 NGROCR(I) = -1
938 NSETCR(I) = -1
939 ELSEIF(I.EQ. 75 ) THEN
940 NGROCR(I) = 1
941 NSETCR(I) = I - 65
942 ELSEIF(I.GE. 76 .AND. I.LE. 79) THEN
943 NGROCR(I) = -1
944 NSETCR(I) = -1
945 ELSEIF(I.EQ. 80 ) THEN
946 NGROCR(I) = 6
947 NSETCR(I) = I - 79
948 ELSEIF(I.GE. 81 .AND. I.LE. 85) THEN
949 NGROCR(I) = -1
950 NSETCR(I) = -1
951 ELSEIF(I.GE. 86 .AND. I.LE. 87) THEN
952 NGROCR(I) = 7
953 NSETCR(I) = I - 85
954 ELSEIF(I.GE. 88 .AND. I.LE. 99) THEN
955 NGROCR(I) = -1
956 NSETCR(I) = -1
957 ENDIF
958 ELSEIF(I.GE.100 .AND. I.LE.199) THEN
959 IF(I.GE.100 .AND. I.LE.105) THEN
960 NGROCR(I) = -1
961 NSETCR(I) = -1
962 ELSEIF(I.GE.106 .AND. I.LE.107) THEN
963 NGROCR(I) = 1
964 NSETCR(I) = I - 105
965 ELSEIF(I.GE.108 .AND. I.LE.129) THEN
966 NGROCR(I) = -1
967 NSETCR(I) = -1
968 ELSEIF(I.GE.130 .AND. I.LE.132) THEN
969 NGROCR(I) = 3
970 NSETCR(I) = I - 129
971 ELSEIF(I.GE.133 .AND. I.LE.171) THEN
972 NGROCR(I) = -1
973 NSETCR(I) = -1
974 ELSEIF(I.GE.172 .AND. I.LE.173) THEN
975 NGROCR(I) = 5
976 NSETCR(I) = I - 171
977 ELSEIF(I.GE.174 .AND. I.LE.179) THEN
978 NGROCR(I) = -1
979 NSETCR(I) = -1
980 ELSEIF(I.GE.180 .AND. I.LE.182) THEN
981 NGROCR(I) = 6
982 NSETCR(I) = I - 179
983 ELSEIF(I.GE.183 .AND. I.LE.199) THEN
984 NGROCR(I) = -1
985 NSETCR(I) = -1
986 ENDIF
987 ELSEIF(I.GE.200 .AND. I.LE.(MODEMX-1)) THEN
988 IF(I.GE.200 .AND. I.LE.205) THEN
989 NGROCR(I) = -1
990 NSETCR(I) = -1
991 ELSEIF(I.GE.206 .AND. I.LE.207) THEN
992 NGROCR(I) = 1
993 NSETCR(I) = I - 205
994 ELSEIF(I.GE.208 .AND. I.LE.209) THEN
995 NGROCR(I) = -1
996 NSETCR(I) = -1
997 ELSEIF(I.GE.210 .AND. I.LE.213) THEN
998 NGROCR(I) = 2
999 NSETCR(I) = I - 209
1000 ELSEIF(I.GE.214 .AND. I.LE.219) THEN
1001 NGROCR(I) = -1
1002 NSETCR(I) = -1
1003 ELSEIF(I.GE.220 .AND. I.LE.222) THEN
1004 NGROCR(I) = 3
1005 NSETCR(I) = I - 219
1006 ELSEIF(I.GE.223 .AND. I.LE.229) THEN
1007 NGROCR(I) = -1
1008 NSETCR(I) = -1
1009 ELSEIF(I.GE.230 .AND. I.LE.232) THEN
1010 NGROCR(I) = 4
1011 NSETCR(I) = I - 229
1012 ELSEIF(I.GE.233 .AND. I.LE.270) THEN
1013 NGROCR(I) = -1
1014 NSETCR(I) = -1
1015 ELSEIF(I.GE.271 .AND. I.LE.273) THEN
1016 NGROCR(I) = 5
1017 NSETCR(I) = I - 270
1018 ELSEIF(I.GE.274 .AND. I.LE.279) THEN
1019 NGROCR(I) = -1
1020 NSETCR(I) = -1
1021 ELSEIF(I.GE.280 .AND. I.LE.281) THEN
1022 NGROCR(I) = 6
1023 NSETCR(I) = I - 279
1024 ELSEIF(I.GE.282 .AND. I.LE.(MODEMX-1)) THEN
1025 NGROCR(I) = -1
1026 NSETCR(I) = -1
1027 ENDIF
1028 ELSEIF(I.EQ.MODEMX) THEN
1029 NPTYCR(I) = -1
1030 NGROCR(I) = -1
1031 NSETCR(I) = -1
1032 ENDIF
1033 ENDDO
1034C.
1035C Now get opposite cross reference
1036C (give NEW to get OLD)
1037 CALL UFILL(MODECR,1,NPTYMX*NGRMAX*NSETMX,-1)
1038 DO I = 1,NPTYMX
1039 DO J = 1,NGRMAX
1040 IF (I.EQ.1) THEN
1041 IF (J.EQ.1) KT = 0
1042 IF (J.EQ.2) KT = 9
1043 IF (J.EQ.3) KT = 19
1044 IF (J.EQ.4) KT = 49
1045 IF (J.EQ.5) KT = 69
1046 IF (J.EQ.6) KT = 79
1047 IF (J.EQ.7) KT = 85
1048 ELSEIF (I.EQ.2) THEN
1049 IF (J.EQ.1) KT = 105
1050 IF (J.EQ.3) KT = 129
1051 IF (J.EQ.5) KT = 171
1052 IF (J.EQ.6) KT = 179
1053 ELSEIF (I.EQ.3) THEN
1054 IF (J.EQ.1) KT = 205
1055 IF (J.EQ.2) KT = 209
1056 IF (J.EQ.3) KT = 219
1057 IF (J.EQ.4) KT = 229
1058 IF (J.EQ.5) KT = 270
1059 IF (J.EQ.6) KT = 279
1060 ENDIF
1061 IF(I.EQ.4 .OR. NPGSMX(I,J).EQ.0) GOTO 5
1062 DO K = 1,NPGSMX(I,J)
1063 KT = KT + 1
1064 MK = KT
1065 IF (I.EQ.1 .AND. J.EQ.1 .AND. K.EQ.10) MK = 75
1066 IF (I.EQ.1 .AND. J.EQ.3 .AND. K.GE.29) MK = -1
1067 IF (I.EQ.1 .AND. J.EQ.4 .AND. K.GE.12) MK = -1
1068 IF (I.EQ.1 .AND. J.EQ.5 .AND. K.GE. 5) MK = -1
1069 IF (I.EQ.3 .AND. J.EQ.5 .AND. K.GE. 4) MK = -1
1070 IF (I.EQ.3 .AND. J.EQ.6 .AND. K.GE. 3) MK = -1
1071 IF (I.EQ.3 .AND. J.EQ.8 .AND. K.GE. 1) MK = -1
1072 IF (I.EQ.3 .AND. J.EQ.9 .AND. K.GE. 1) MK = -1
1073 MODECR(I,J,K) = MK
1074 ENDDO
1075 5 CONTINUE
1076 ENDDO
1077 ENDDO
1078C Define Printer
1079 IF(N6.LE.0) N6 = L6
1080C Define FIRST flag
1081 FIRST = .TRUE.
1082C Set the default values
1083 IF(IFLPRT.LT.0 .OR. IFLPRT.GT.5) IFLPRT = LFLPRT
1084 NPTYPE = LPTYPE
1085 NGROUP = LGROUP
1086 NSET = LNSET
1087 NFL = IFL
1088 LO = LORD
1089 TMAS = TM
1090 QCDL4 = PQCDL4(NPTYPE,NGROUP,NSET)
1091 QCDL5 = PQCDL5(NPTYPE,NGROUP,NSET)
1092 XMIN = PXMIN(NPTYPE,NGROUP,NSET)
1093 XMAX = PXMAX(NPTYPE,NGROUP,NSET)
1094 Q2MIN = PQ2MIN(NPTYPE,NGROUP,NSET)
1095 Q2MAX = PQ2MAX(NPTYPE,NGROUP,NSET)
1096 NATYPE = LATYPE
1097 NAGROUP = LAGROUP
1098 NASET = LNASET
1099C.
1100C Save the PDFLIB parameters
1101 IFLPRTP = IFLPRT
1102 NPTYPP = NPTYPE
1103 NGROPP = NGROUP
1104 NSETP = NSET
1105 NFLP = NFL
1106 LOP = LO
1107 TMASP = TMAS
1108 QCDL4P = QCDL4
1109 QCDL5P = QCDL5
1110 XMINP = XMIN
1111 XMAXP = XMAX
1112 Q2MINP = Q2MIN
1113 Q2MAXP = Q2MAX
1114 NATYPP = NATYPE
1115 NAGROPP = NAGROUP
1116 NASETP = NASET
1117C. endif from ISTART
1118 ENDIF
1119C.
1120C Validate call to PDFSET
1121 IFLSET = 1
1122C.
1123C Default : NEW version 4.0 of PDFLIB
1124 NEWVER = .TRUE.
1125C.
1126C Now get user defined values
1127 NFLU = 0
1128 LOU = 0
1129 TMASU = ZEROD
1130 QCDL4U = ZEROD
1131 QCDL5U = ZEROD
1132 XMINU = ZEROD
1133 XMAXU = ZEROD
1134 Q2MINU = ZEROD
1135 Q2MAXU = ZEROD
1136 KPTYP = 0
1137 KGROP = 0
1138 KNSET = 0
1139 KMODE = 0
1140 NV = 0
1141 DO 10 I = 1,20
1142 STRING = PARM(I)
1143 CALL CLTOU(STRING)
1144 NSCHECK = 0
1145 DO NS = 1,NSDIM
1146 IF(STRING(1:2).EQ.AV(NS)) NSCHECK = NS
1147 ENDDO
1148 IF( STRING(1:5).NE.'INIT0' .AND. STRING(1:4).NE.'MODE'
1149 + .AND. STRING(1:6).NE.'NPTYPE'.AND. STRING(1:6).NE.'NGROUP'
1150 + .AND. STRING(1:4).NE.'NSET' .AND. STRING(1:3).NE.'NFL'
1151 + .AND. STRING(1:2).NE.'LO' .AND. STRING(1:4).NE.'TMAS'
1152 + .AND. STRING(1:5).NE.'QCDL4' .AND. STRING(1:5).NE.'QCDL5'
1153 + .AND. STRING(1:4).NE.'XMIN' .AND. STRING(1:4).NE.'XMAX'
1154 + .AND. STRING(1:5).NE.'Q2MIN' .AND. STRING(1:5).NE.'Q2MAX'
1155 + .AND. STRING(1:6).NE.'NATYPE'.AND. STRING(1:7).NE.'NAGROUP'
1156 + .AND. STRING(1:5).NE.'NASET'
1157C.
1158 + .AND. STRING(1:2).NE.'NU' .AND. STRING(1:2).NE.'PR'
1159 + .AND. STRING(1:2).NE.'PI' .AND. STRING(1:2).NE.'PH'
1160 + .AND. STRING(1:3).NE.'PRI' .AND. NSCHECK .EQ.0
1161 + .AND. STRING(1:2).NE.'OW' .AND. STRING(1:2).NE.'SM'
1162 + .AND. STRING(1:4).NE.'GRVP' .AND. STRING(1:5).NE.'GRV-P'
1163 + .AND. STRING(1:5).NE.'GRVPI' .AND. STRING(1:6).NE.'GRV-PI'
1164 + .AND. STRING(1:4).NE.'ABFK' .AND. STRING(1:3).NE.'DOG'
1165C + .AND. STRING(1:4).NE.'DO-G' .AND. STRING(1:4).NE.'DOPH'
1166C + .AND. STRING(1:5).NE.'DO-PH' .AND. STRING(1:2).NE.'DG'
1167 + .AND. STRING(1:2).NE.'LA' .AND. STRING(1:2).NE.'GS'
1168 + .AND. STRING(1:3).NE.'GRV' .AND. STRING(1:3).NE.'GRS'
1169 + .AND. STRING(1:4).NE.'WHIT' .AND. STRING(1:3).NE.'SAS'
1170 + .AND. STRING(1:2).NE.'AC' .AND. STRING(1:5).NE.'FIRST')GOTO 20
1171 NV = NV + 1
1172 VALUE = VAL(I)
1173C.
1174 IF(STRING(1:5).EQ.'FIRST' .OR. STRING(1:4).EQ.'MODE' .OR.
1175 + STRING(1:6).EQ.'NPTYPE'.OR. STRING(1:6).EQ.'NGROUP' .OR.
1176 + STRING(1:4).EQ.'NSET' .OR. STRING(1:3).EQ.'NFL' .OR.
1177 + STRING(1:2).EQ.'LO' .OR.
1178 + STRING(1:6).EQ.'NATYPE'.OR. STRING(1:7).EQ.'NAGROUP' .OR.
1179 + STRING(1:5).EQ.'NASET' .OR.
1180 + STRING(1:3).EQ.'PRI' .OR. NSCHECK .GT.0 .OR.
1181 + STRING(1:2).EQ.'OW' .OR. STRING(1:2).EQ.'SM' .OR.
1182 + STRING(1:4).EQ.'GRVP' .OR. STRING(1:5).EQ.'GRV-P' .OR.
1183 + STRING(1:5).EQ.'GRVPI' .OR. STRING(1:6).EQ.'GRV-PI' .OR.
1184 + STRING(1:4).EQ.'ABFK' .OR. STRING(1:3).EQ.'DOG' .OR.
1185C + STRING(1:4).EQ.'DO-G' .OR. STRING(1:4).EQ.'DOPH' .OR.
1186C + STRING(1:5).EQ.'DO-PH' .OR. STRING(1:2).EQ.'DG' .OR.
1187 + STRING(1:2).EQ.'LA' .OR. STRING(1:2).EQ.'GS' .OR.
1188 + STRING(1:3).EQ.'GRV' .OR. STRING(1:3).EQ.'GRS' .OR.
1189 + STRING(1:4).EQ.'WHIT' .OR. STRING(1:3).EQ.'SAS' .OR.
1190 + STRING(1:2).EQ.'AC' ) THEN
1191 TEMP = VALUE
1192 IF(ABS(VALUE).LT.1.E-2) TEMP = IVALUE
1193 ENDIF
1194C Get user's setting for FIRST (option: silent mode => FIRST = .FALSE.)
1195 IF(STRING(1:5).EQ.'FIRST') THEN
1196 IF(TEMP.EQ.0) FIRST = .FALSE.
1197 IF(TEMP.EQ.1) FIRST = .TRUE.
1198 ENDIF
1199C.
1200C. Return, because PDFSET was called ONLY for initialisation
1201 IF(STRING(1:5).EQ.'INIT0') THEN
1202 IFLSET = 0
1203C Print status of PDFLIB
1204 IF(FIRST) THEN
1205 IF(ISTART.EQ.0) CALL PDFVERS
1206 ENDIF
1207 ISTART = 1
1208C. Reset PDFSET parameters
1209 DO K = 1,NCHDIM
1210 PARM(K) = ' '
1211 ENDDO
1212 RETURN
1213 ENDIF
1214C.
1215C Check if PDFLIB format of new version 4.0 or of versions 3.0 or less
1216 IF(STRING(1:4).EQ.'MODE') THEN
1217 NEWVER = .FALSE.
1218 MODE = TEMP
1219 KMODE = 1
1220 ENDIF
1221C Now NEW version 4.0 or bigger
1222C Get particle type
1223 IF(STRING(1:6).EQ.'NPTYPE') THEN
1224 NPTYPE = TEMP
1225 KPTYP = 1
1226 ENDIF
1227C Get Author Group Number
1228 IF(STRING(1:6).EQ.'NGROUP') THEN
1229 NGROUP = TEMP
1230 KGROP = 1
1231 ENDIF
1232C Get PDF Set within Author Group
1233 IF(STRING(1:4).EQ.'NSET' ) THEN
1234 NSET = TEMP
1235 KNSET = 1
1236 ENDIF
1237C Particle Type = Nucleons
1238 IF(STRING(1:2).EQ.'NU' .OR. STRING(1:2).EQ.'PR'
1239 + .OR. NSCHECK.GT.0) THEN
1240 NPTYPE = 1
1241 KPTYP = 1
1242 ENDIF
1243 IF(STRING(1:3).EQ.'PRI') THEN
1244 NPTYPE = 1
1245 KPTYP = 1
1246 NGROUP = 1
1247 KGROP = 1
1248 NSET = TEMP
1249 KNSET = 1
1250 ENDIF
1251 IF(NSCHECK.GE.1 .AND. NSCHECK.LE.5) THEN
1252 NPTYPE = 1
1253 KPTYP = 1
1254 NGROUP = 1
1255 KGROP = 1
1256 NSET = TEMP
1257 KNSET = 1
1258 ENDIF
1259 IF(NSCHECK.EQ.6) THEN
1260 NPTYPE = 1
1261 KPTYP = 1
1262 NGROUP = 2
1263 KGROP = 1
1264 NSET = TEMP
1265 KNSET = 1
1266 ENDIF
1267 IF(NSCHECK.GE.7 .AND. NSCHECK.LE.9) THEN
1268 NPTYPE = 1
1269 KPTYP = 1
1270 NGROUP = 3
1271 KGROP = 1
1272 NSET = TEMP
1273 KNSET = 1
1274 ENDIF
1275 IF(NSCHECK.GE.10 .AND. NSCHECK.LE.11) THEN
1276 NPTYPE = 1
1277 KPTYP = 1
1278 NGROUP = 4
1279 KGROP = 1
1280 NSET = TEMP
1281 KNSET = 1
1282 ENDIF
1283 IF(NSCHECK.EQ.12) THEN
1284 NPTYPE = 1
1285 KPTYP = 1
1286 NGROUP = 5
1287 KGROP = 1
1288 NSET = TEMP
1289 KNSET = 1
1290 ENDIF
1291 IF(NSCHECK.EQ.13) THEN
1292 NPTYPE = 1
1293 KPTYP = 1
1294 NGROUP = 6
1295 KGROP = 1
1296 NSET = TEMP
1297 KNSET = 1
1298 ENDIF
1299 IF(NSCHECK.EQ.14) THEN
1300 NPTYPE = 1
1301 KPTYP = 1
1302 NGROUP = 7
1303 KGROP = 1
1304 NSET = TEMP
1305 KNSET = 1
1306 ENDIF
1307C Particle Type = Pions
1308 IF(STRING(1:2).EQ.'PI') THEN
1309 NPTYPE = 2
1310 KPTYP = 1
1311 ENDIF
1312 IF(STRING(1:2).EQ.'OW' .OR.
1313 + STRING(1:3).EQ.'OWP' .OR. STRING(1:4).EQ.'OW-P' .OR.
1314 + STRING(1:4).EQ.'OWPI' .OR. STRING(1:5).EQ.'OW-PI') THEN
1315 NPTYPE = 2
1316 KPTYP = 1
1317 NGROUP = 1
1318 KGROP = 1
1319 NSET = TEMP
1320 KNSET = 1
1321 ENDIF
1322 IF(STRING(1:2).EQ.'SM' .OR.
1323 + STRING(1:5).EQ.'SMRSP' .OR. STRING(1:6).EQ.'SMRS-P' .OR.
1324 + STRING(1:6).EQ.'SMRSPI' .OR. STRING(1:7).EQ.'SMRS-PI') THEN
1325 NPTYPE = 2
1326 KPTYP = 1
1327 NGROUP = 3
1328 KGROP = 1
1329 NSET = TEMP
1330 KNSET = 1
1331 ENDIF
1332 IF(STRING(1:4).EQ.'GRVP' .OR. STRING(1:5).EQ.'GRV-P' .OR.
1333 + STRING(1:5).EQ.'GRVPI' .OR. STRING(1:6).EQ.'GRV-PI') THEN
1334 NPTYPE = 2
1335 KPTYP = 1
1336 NGROUP = 5
1337 KGROP = 1
1338 NSET = TEMP
1339 KNSET = 1
1340 ENDIF
1341 IF(STRING(1:4).EQ.'ABFK' .OR.
1342 + STRING(1:6).EQ.'ABFKWP' .OR. STRING(1:7).EQ.'ABFKW-P' .OR.
1343 + STRING(1:7).EQ.'ABFKWPI' .OR. STRING(1:8).EQ.'ABFKW-PI') THEN
1344 NPTYPE = 2
1345 KPTYP = 1
1346 NGROUP = 6
1347 KGROP = 1
1348 NSET = TEMP
1349 KNSET = 1
1350 ENDIF
1351C Particle Type = Photons
1352 IF(STRING(1:2).EQ.'PH') THEN
1353 NPTYPE = 3
1354 KPTYP = 1
1355 ENDIF
1356 IF(STRING(1:3).EQ.'DOG' .OR. STRING(1:4).EQ.'DO-G' .OR.
1357 + STRING(1:4).EQ.'DOPH' .OR. STRING(1:5).EQ.'DO-PH') THEN
1358 NPTYPE = 3
1359 KPTYP = 1
1360 NGROUP = 1
1361 KGROP = 1
1362 NSET = TEMP
1363 KNSET = 1
1364 ENDIF
1365 IF(STRING(1:2).EQ.'DG' .OR.
1366 + STRING(1:3).EQ.'DGG' .OR. STRING(1:4).EQ.'DG-G' .OR.
1367 + STRING(1:4).EQ.'DGPH' .OR. STRING(1:5).EQ.'DG-PH') THEN
1368 NPTYPE = 3
1369 KPTYP = 1
1370 NGROUP = 2
1371 KGROP = 1
1372 NSET = TEMP
1373 KNSET = 1
1374 ENDIF
1375 IF(STRING(1:2).EQ.'LA' .OR.
1376 + STRING(1:4).EQ.'LACG' .OR. STRING(1:5).EQ.'LAC-G' .OR.
1377 + STRING(1:5).EQ.'LACPH' .OR. STRING(1:6).EQ.'LAC-PH') THEN
1378 NPTYPE = 3
1379 KPTYP = 1
1380 NGROUP = 3
1381 KGROP = 1
1382 NSET = TEMP
1383 KNSET = 1
1384 ENDIF
1385 IF(STRING(1:2).EQ.'GS' .OR.
1386 + STRING(1:3).EQ.'GSG' .OR. STRING(1:4).EQ.'GS-G' .OR.
1387 + STRING(1:4).EQ.'GSPH' .OR. STRING(1:5).EQ.'GS-PH') THEN
1388 NPTYPE = 3
1389 KPTYP = 1
1390 NGROUP = 4
1391 KGROP = 1
1392 NSET = TEMP
1393 KNSET = 1
1394 ENDIF
1395 IF(STRING(1:4).EQ.'GRVG' .OR. STRING(1:5).EQ.'GRV-G' .OR.
1396 + STRING(1:5).EQ.'GRVPH' .OR. STRING(1:6).EQ.'GRV-PH' .OR.
1397 + STRING(1:5).EQ.'GRSPH' .OR. STRING(1:6).EQ.'GRS-PH') THEN
1398 NPTYPE = 3
1399 KPTYP = 1
1400 NGROUP = 5
1401 KGROP = 1
1402 NSET = TEMP
1403 KNSET = 1
1404 ENDIF
1405 IF(STRING(1:2).EQ.'AC' .OR.
1406 + STRING(1:6).EQ.'ACFGPG' .OR. STRING(1:7).EQ.'ACFGP-G' .OR.
1407 + STRING(1:7).EQ.'ACFGPPH' .OR. STRING(1:8).EQ.'ACFGP-PH') THEN
1408 NPTYPE = 3
1409 KPTYP = 1
1410 NGROUP = 6
1411 KGROP = 1
1412 NSET = TEMP
1413 KNSET = 1
1414 ENDIF
1415 IF(STRING(1:2).EQ.'WI' .OR.
1416 + STRING(1:5).EQ.'WHITG' .OR. STRING(1:6).EQ.'WHIT-G' .OR.
1417 + STRING(1:6).EQ.'WHITPH' .OR. STRING(1:7).EQ.'WHIT-PH') THEN
1418 NPTYPE = 3
1419 KPTYP = 1
1420 NGROUP = 8
1421 KGROP = 1
1422 NSET = TEMP
1423 KNSET = 1
1424 ENDIF
1425 IF(STRING(1:2).EQ.'SA' .OR.
1426 + STRING(1:4).EQ.'SASG' .OR. STRING(1:5).EQ.'SAS-G' .OR.
1427 + STRING(1:5).EQ.'SASPH' .OR. STRING(1:6).EQ.'SAS-PH') THEN
1428 NPTYPE = 3
1429 KPTYP = 1
1430 NGROUP = 9
1431 KGROP = 1
1432 NSET = TEMP
1433 KNSET = 1
1434 ENDIF
1435C
1436C Get particle type for Nuclear Corrections
1437 IF(STRING(1:6).EQ.'NATYPE') THEN
1438 NATYPE = TEMP
1439 KPTYP = 1
1440 NAFLAG = 1
1441 ENDIF
1442C Get Author Group Number for Nuclear Corrections
1443 IF(STRING(1:7).EQ.'NAGROUP') THEN
1444 NAGROUP = TEMP
1445 KGROP = 1
1446 ENDIF
1447C Get PDF Set within Author Group for Nuclear Corrections
1448 IF(STRING(1:5).EQ.'NASET' ) THEN
1449 NASET = TEMP
1450 KNSET = 1
1451 ENDIF
1452C
1453C... set defaults for private structure functions set (Nset = 1)
1454 IF(NPTYPE.EQ.1 .AND. NGROUP.EQ.1 .AND. NSET.EQ.1) THEN
1455 NFL = IFL
1456 LO = LORD
1457 TMAS = TM
1458 QCDL4 = ZEROD
1459 QCDL5 = ZEROMIN
1460 XMIN = ZEROD
1461 XMAX = 0.9999D0
1462 Q2MIN = ZEROD
1463 Q2MAX = 1.000D8
1464 ENDIF
1465 IF(STRING(1:3).EQ.'NFL' ) NFLU = TEMP
1466 IF(STRING(1:2).EQ.'LO' ) LOU = TEMP
1467 IF(STRING(1:4).EQ.'TMAS' ) TMASU = VALUE
1468C.
1469 IF(STRING(1:5).EQ.'QCDL4') QCDL4U = VALUE
1470 IF(STRING(1:5).EQ.'QCDL5') QCDL5U = VALUE
1471C.
1472 IF(STRING(1:4).EQ.'XMIN') XMINU = VALUE
1473 IF(STRING(1:4).EQ.'XMAX') XMAXU = VALUE
1474 IF(STRING(1:5).EQ.'Q2MIN') Q2MINU = VALUE
1475 IF(STRING(1:5).EQ.'Q2MAX') Q2MAXU = VALUE
1476 10 CONTINUE
1477 20 CONTINUE
1478C Print status of PDFLIB
1479 IF(FIRST) THEN
1480 IF(ISTART.EQ.0) CALL PDFVERS
1481 ENDIF
1482 ISTART = 1
1483C Print values
1484 IF(IFLPRT.GE.1) THEN
1485 IF(NV.GT.0) THEN
1486 IF(NEWVER) THEN
1487 WRITE(N6,2000) (PARM(I)(1:6),I=1,NV)
1488 ELSE
1489 WRITE(N6,2000) (PARM(I)(1:4),I=1,NV)
1490 ENDIF
1491 WRITE(N6,2001) (VAL(I),I=1,NV)
1492 ELSE
1493 WRITE(N6,2002)
1494 ENDIF
1495 ELSE
1496 IF(NV.EQ.0) WRITE(N6,2002)
1497 ENDIF
1498 2000 FORMAT(/,' Parm = ',10(1X,A8))
1499 2001 FORMAT(1H ,' Val = ',10(1X,F8.4))
1500 2002 FORMAT(/,' PDFSET : Warning!! NO Parameters read in, defaults
1501 + (Nucleon PDFs) used !!')
1502C.
1503C. check on conflicting parameters in case of mixed versions
1504 IF((KPTYP*KMODE).EQ.1 .OR. (KGROP*KMODE).EQ.1 .OR.
1505 + (KNSET*KMODE).EQ.1) THEN
1506 WRITE(N6,*)
1507 + ' PDFLIB : Warning!! Mixed parameter settings for OLD and NE
1508 +W version of PDFLIB format'
1509 WRITE(N6,*)
1510 + ' Warning : ALL parameters set to NEW default (Nucleon PDFs) !!'
1511 NEWVER = .TRUE.
1512 NPTYPE = LPTYPE
1513 NGROUP = LGROUP
1514 NSET = LNSET
1515 NATYPE = LATYPE
1516 NAGROUP = LAGROUP
1517 NASET = LNASET
1518 ENDIF
1519C.
1520#include "pdf/w50511c.inc"
1521C.
1522C User wants old version (3.0 or less) of PDFLIB format
1523 IF(.NOT.NEWVER) THEN
1524 IF (MODE .GE.0 .AND. MODE .LE.MODEMX) THEN
1525 IF(MODE.EQ.0) THEN
1526 NPTYPE = LPTYPO
1527 NGROUP = LGROPO
1528 NSET = LNSETO
1529 ELSE
1530 NPTYPE = NPTYCR(MODE)
1531 NGROUP = NGROCR(MODE)
1532 NSET = NSETCR(MODE)
1533C Check on validity of parameter values
1534 IF (NPTYPE.LT.0 .OR. NGROUP.LT.0 .OR. NSET.LT.0) THEN
1535 IF(FIRST) THEN
1536 WRITE(N6,*) ' PDFLIB : MODE value INCORRECT, MODE = ',
1537 + MODE
1538 WRITE(N6,*) ' Warning : MODE value set to OLD default !!'
1539 ENDIF
1540 NPTYPE = LPTYPO
1541 NGROUP = LGROPO
1542 NSET = LNSETO
1543 GOTO 30
1544 ENDIF
1545 ENDIF
1546 ELSE
1547 IF(FIRST) THEN
1548 WRITE(N6,*) ' PDFLIB : MODE value INCORRECT, MODE = ',MODE
1549 WRITE(N6,*) ' Warning : MODE value set to OLD default !!'
1550 ENDIF
1551 NPTYPE = LPTYPO
1552 NGROUP = LGROPO
1553 NSET = LNSETO
1554 ENDIF
1555C User wants new version (4.0 or bigger) of PDFLIB format
1556 ELSE
1557C Check on validity of parameter values
1558 IF (NSET.EQ.0) THEN
1559 IF(FIRST) THEN
1560 WRITE(N6,*)
1561 + ' PDFLIB : DEFAULT parameter set (Nucleon PDFs) has been
1562 + chosen !'
1563 WRITE(N6,*)
1564 + ' PDFLIB : Nptype = ',NPTYPE,' Ngroup = ',NGROUP,
1565 + ' Nset = ',NSET
1566 ENDIF
1567 NPTYPE = LPTYPE
1568 NGROUP = LGROUP
1569 NSET = LNSET
1570 ENDIF
1571 IF (NPTYPE.LT.0 .OR. NPTYPE.EQ.0 .OR. NPTYPE.GT.NPTYMX) THEN
1572 IF(FIRST) THEN
1573 WRITE(N6,*)
1574 + ' PDFLIB : PARTICLE TYPE number INCORRECT, Nptype = ',
1575 + NPTYPE,' Ngroup = ',NGROUP,' Nset = ',NSET
1576 WRITE(N6,*)
1577 + ' Warning : ALL parameters set to default (Nucleon PDFs) !!'
1578 ENDIF
1579 NPTYPE = LPTYPE
1580 NGROUP = LGROUP
1581 NSET = LNSET
1582 ENDIF
1583 IF (NGROUP.LT.0 .OR. NGROUP.EQ.0 .OR. NGROUP.GT.NGRMAX) THEN
1584 IF(FIRST) THEN
1585 WRITE(N6,*)
1586 + ' PDFLIB : GROUP number INCORRECT, Nptype = ',
1587 + NPTYPE,' Ngroup = ',NGROUP,' Nset = ',NSET
1588 WRITE(N6,*)
1589 + ' Warning : ALL parameters set to default (Nucleon PDFs) !!'
1590 ENDIF
1591 NPTYPE = LPTYPE
1592 NGROUP = LGROUP
1593 NSET = LNSET
1594 ENDIF
1595 IF (NSET.LT.0 .OR. NSET.EQ.0 .OR. NSET.GT.NPGSMX(NPTYPE,NGROUP))
1596 + THEN
1597 IF(FIRST) THEN
1598 WRITE(N6,*)
1599 + ' PDFLIB : NSET value for PDFs INCORRECT, Nptype = ',
1600 + NPTYPE,' Ngroup = ',NGROUP,' Nset = ',NSET
1601 WRITE(N6,*)
1602 + ' Warning : ALL parameters set to default (Nucleon PDFs) !!'
1603 ENDIF
1604 NPTYPE = LPTYPE
1605 NGROUP = LGROUP
1606 NSET = LNSET
1607 ENDIF
1608C Check on validity of parameter values in case of NPTYPE=4 (wrong)
1609 IF (NPTYPE.EQ.4) THEN
1610 IF(FIRST) THEN
1611 WRITE(N6,*)
1612 + ' PDFLIB : User may want NPDFs with Nuclear Corrections !!'
1613 WRITE(N6,*)
1614 + ' PDFLIB : Nptype = ',NPTYPE,' Ngroup = ',NGROUP,
1615 + ' Nset = ',NSET
1616 WRITE(N6,*)
1617 + ' PDFLIB : Nptype = 4 is not allowed !!'
1618 WRITE(N6,*)
1619 + ' PDFLIB : Please check ALL INPUT parameters for PDFSET and
1620 + STRUCTA !!'
1621 ENDIF
1622C. Reset PDFSET parameters
1623 DO K = 1,NCHDIM
1624 PARM(K) = ' '
1625 ENDDO
1626 RETURN
1627 ENDIF
1628C now check on Nuclear Corrections
1629 IF (NASET.EQ.0) THEN
1630 IF(FIRST) THEN
1631 WRITE(N6,*)
1632 + ' PDFLIB : DEFAULT parameter set has been chosen for Nuclea
1633 +r Corrections!'
1634 WRITE(N6,*)
1635 + ' PDFLIB : NAtype = ',NATYPE,' NAgroup = ',NAGROUP,
1636 + ' NAset = ',NASET
1637 ENDIF
1638 NATYPE = LATYPE
1639 NAGROUP = LAGROUP
1640 NASET = LNASET
1641 ENDIF
1642 IF (NATYPE.LT.0 .OR. NATYPE.EQ.0 .OR. NATYPE.GT.NPTYMX)
1643 +THEN
1644 IF(FIRST) THEN
1645 WRITE(N6,*)
1646 + ' PDFLIB : PARTICLE TYPE number INCORRECT, NAtype = ',
1647 + NATYPE,' NAgroup = ',NAGROUP,' NAset = ',NASET
1648 WRITE(N6,*)
1649 + ' Warning : ALL parameters set to default for Nuclear Correctio
1650 +ns!!'
1651 ENDIF
1652 NATYPE = LATYPE
1653 NAGROUP = LAGROUP
1654 NASET = LNASET
1655 ENDIF
1656 IF (NAGROUP.LT.0 .OR. NAGROUP.EQ.0 .OR. NAGROUP.GT.NGRMAX)
1657 +THEN
1658 IF(FIRST) THEN
1659 WRITE(N6,*)
1660 + ' PDFLIB : NAGROUP number INCORRECT, Nptype = ',
1661 + NATYPE,' NAgroup = ',NAGROUP,' NAset = ',NASET
1662 WRITE(N6,*)
1663 + ' Warning : ALL parameters set to default for Nuclear Correctio
1664 +ns!!'
1665 ENDIF
1666 NATYPE = LATYPE
1667 NAGROUP = LAGROUP
1668 NASET = LNASET
1669 ENDIF
1670 IF (NASET.LT.0 .OR. NASET.EQ.0 .OR.
1671 + NASET.GT.NPGSMX(NATYPE,NAGROUP)) THEN
1672 IF(FIRST) THEN
1673 WRITE(N6,*)
1674 + ' PDFLIB : NASET value for PDFs INCORRECT, NAtype = ',
1675 + NATYPE,' NAgroup = ',NAGROUP,' NAset = ',NASET
1676 WRITE(N6,*)
1677 + ' Warning : ALL parameters set to default for Nuclear Correctio
1678 +ns !!'
1679 ENDIF
1680 NATYPE = LATYPE
1681 NAGROUP = LAGROUP
1682 NASET = LNASET
1683 ENDIF
1684 ENDIF
1685 30 CONTINUE
1686C.
1687C Check on retracted PDF sets
1688 IF (NSETFL(NPTYPE,NGROUP,NSET).EQ.0) THEN
1689 IF(FIRST) THEN
1690 WRITE(N6,*)
1691 + ' PDFLIB : Your selected PDF set has been retracted by the aut
1692 +hors'
1693 WRITE(N6,*)
1694 + ' Warning : ALL parameters set to default (Nucleon PDFs) !!'
1695 ENDIF
1696 NPTYPE = LPTYPE
1697 NGROUP = LGROUP
1698 NSET = LNSET
1699 ENDIF
1700C.
1701C Number of flavours in alpha(s) calculation
1702 IF (ABS(NFL).LT.3 .OR. ABS(NFL).GT.6) THEN
1703 IF(FIRST) THEN
1704 IF(NPTYPE.EQ.1 .AND. NGROUP.EQ.1 .AND. NSET.EQ.1) WRITE(N6,*)
1705 + ' You are calling the DUMMY private SF routine !!'
1706 WRITE(N6,*) ' PDFLIB : NFL value outside range, NFL = ',NFL
1707 WRITE(N6,*) ' Warning : NFL value set to default !!'
1708 ENDIF
1709 NFL = IFL
1710 ENDIF
1711C.
1712C Special settings for Number of Flavours for CTEQ4 and CTEQ55 F3 and F4
1713C
1714 IF (FIRST) THEN
1715 NFL = IFL
1716 IF(NPTYPE.EQ.1 .AND. NGROUP.EQ.4 .AND. NSET.EQ.44) NFL=-3
1717 IF(NPTYPE.EQ.1 .AND. NGROUP.EQ.4 .AND. NSET.EQ.45) NFL=-4
1718 IF(NPTYPE.EQ.1 .AND. NGROUP.EQ.4 .AND. NSET.EQ.51) NFL=-3
1719 IF(NPTYPE.EQ.1 .AND. NGROUP.EQ.4 .AND. NSET.EQ.52) NFL=-4
1720 ENDIF
1721C.
1722C User defined Number of Flavours
1723 IF (NFLU .GT. 0) THEN
1724 NFL = NFLU
1725 IF(FIRST) THEN
1726 IF(NFLU.NE.NFLP) THEN
1727 WRITE(N6,*) ' PDFLIB : NFL value'
1728 WRITE(N6,*) ' Warning : NON standard settings, NFL value = ',
1729 + NFL,' set by user !!'
1730 ENDIF
1731 ENDIF
1732 ENDIF
1733C.
1734C Order of alpha(s) calculation
1735 LO=LORD
1736 IF ((NPTYPE.EQ.1.AND.NGROUP.EQ.1.AND.NSET.GE. 2) .OR.
1737 + (NPTYPE.EQ.1.AND.NGROUP.EQ.2.AND.NSET.LE. 5) .OR.
1738 + (NPTYPE.EQ.1.AND.NGROUP.EQ.3.AND.NSET.GE.72.AND.NSET.LE.76).OR.
1739 + (NPTYPE.EQ.1.AND.NGROUP.EQ.4.AND.NSET.GE.11.AND.NSET.LE.12).OR.
1740 + (NPTYPE.EQ.1.AND.NGROUP.EQ.4.AND.NSET.EQ.17) .OR.
1741 + (NPTYPE.EQ.1.AND.NGROUP.EQ.4.AND.NSET.EQ.23) .OR.
1742 + (NPTYPE.EQ.1.AND.NGROUP.EQ.4.AND.NSET.EQ.29) .OR.
1743 + (NPTYPE.EQ.1.AND.NGROUP.EQ.4.AND.NSET.EQ.32) .OR.
1744 + (NPTYPE.EQ.1.AND.NGROUP.EQ.5.AND.NSET.EQ. 2) .OR.
1745 + (NPTYPE.EQ.1.AND.NGROUP.EQ.5.AND.NSET.GE. 4.AND.NSET.LE.5).OR.
1746 + (NPTYPE.EQ.1.AND.NGROUP.EQ.5.AND.NSET.GE. 8.AND.NSET.LE.9).OR.
1747 + (NPTYPE.EQ.2.AND.NGROUP.EQ.1.AND.NSET.LE. 2) .OR.
1748 + (NPTYPE.EQ.2.AND.NGROUP.EQ.5.AND.NSET.EQ. 2) .OR.
1749 + (NPTYPE.EQ.3.AND.NGROUP.EQ.1.AND.NSET.EQ. 1) .OR.
1750 + (NPTYPE.EQ.3.AND.NGROUP.EQ.2.AND.NSET.LE. 4) .OR.
1751 + (NPTYPE.EQ.3.AND.NGROUP.EQ.3.AND.NSET.LE. 3) .OR.
1752 + (NPTYPE.EQ.3.AND.NGROUP.EQ.4.AND.NSET.GE. 2.AND.NSET.LE.3).OR.
1753 + (NPTYPE.EQ.3.AND.NGROUP.EQ.4.AND.NSET.EQ. 5) .OR.
1754 + (NPTYPE.EQ.3.AND.NGROUP.EQ.5.AND.NSET.GE. 3.AND.NSET.LE.4).OR.
1755 + (NPTYPE.EQ.3.AND.(NGROUP.GE.8.AND.NGROUP.LE.9))) LO=1
1756 IF ((NPTYPE.EQ.1.AND.NGROUP.EQ.4.AND.NSET.EQ. 46)) LO=1
1757 IF ((NPTYPE.EQ.1.AND.NGROUP.EQ.4.AND.NSET.EQ. 55)) LO=1
1758 IF ((NPTYPE.EQ.1.AND.NGROUP.EQ.5.AND.NSET.EQ. 12)) LO=1
1759 IF (LOU .GT. 0) THEN
1760 LO = LOU
1761 IF(FIRST) THEN
1762 IF(LOU.NE.LOP) THEN
1763 WRITE(N6,*) ' PDFLIB : LO value'
1764 WRITE(N6,*) ' Warning : NON standard settings, LO value = ',
1765 + LO,' set by user !!'
1766 ENDIF
1767 ENDIF
1768 ENDIF
1769C.
1770C Top Quark mass
1771 IF(TMAS.LE.ZEROD) THEN
1772 IF(FIRST) THEN
1773 WRITE(N6,*)' PDFLIB : TMAS value zero or negativ, Tmas = ',TMAS
1774 WRITE(N6,*)' Warning : TMAS value set to default !!'
1775 ENDIF
1776 TMAS = TM
1777 ENDIF
1778 IF (TMASU .GT. ZEROMIN) THEN
1779 TMAS = TMASU
1780 IF(FIRST) THEN
1781 IF(TMASU.NE.TMASP) THEN
1782 WRITE(N6,*)' PDFLIB : TMAS value'
1783 WRITE(N6,*) ' Warning : NON standard settings, TMAS value = ',
1784 + TMAS,' set by user !!'
1785 ENDIF
1786 ENDIF
1787 ENDIF
1788C.
1789C. this lambda (NFL=4) is given by the authors and
1790C. NOT used in the calculation
1791 QCDL4 = PQCDL4(NPTYPE,NGROUP,NSET)
1792 IF(NPTYPE.EQ.1 .AND. NGROUP.EQ.1 .AND. NSET.EQ.1) THEN
1793 IF (QCDL5 .LT. ZEROMIN) THEN
1794 IF(FIRST) THEN
1795 WRITE (N6,*) ' PDFLIB : QCD-Lambda5 is not defined !! '
1796 WRITE (N6,*) ' Warning : QCD-Lambda5 set to QCDL5 = 0.0D0'
1797 ENDIF
1798 QCDL4 = ZEROD
1799 QCDL5 = ZEROD
1800C. Reset PDFSET parameters
1801 DO K = 1,NCHDIM
1802 PARM(K) = ' '
1803 ENDDO
1804 RETURN
1805 ENDIF
1806 ENDIF
1807 IF (QCDL4U .GT. ZEROMIN) THEN
1808 QCDL4 = QCDL4U
1809 IF(FIRST) THEN
1810 IF(QCDL4U.NE.QCDL4P) THEN
1811 WRITE (N6,*) ' PDFLIB : QCD-Lambda4 value'
1812 WRITE(N6,*) ' Warning : NON standard settings, QCDL4 value = ',
1813 + QCDL4,' set by user !!'
1814 ENDIF
1815 ENDIF
1816 ENDIF
1817C.
1818C. this is the lambda (NFL=5) used in the calculation
1819 QCDL5 = PQCDL5(NPTYPE,NGROUP,NSET)
1820 IF (QCDL5U .GT. ZEROMIN) THEN
1821 QCDL5 = QCDL5U
1822 IF(FIRST) THEN
1823 IF(QCDL5U.NE.QCDL5P) THEN
1824 WRITE (N6,*) ' PDFLIB : QCD-Lambda5 value'
1825 WRITE(N6,*) ' Warning : NON standard settings, QCDL5 value = ',
1826 + QCDL5,' set by user !!'
1827 ENDIF
1828 ENDIF
1829 ENDIF
1830C.
1831C Minimal X values
1832 XMIN = PXMIN(NPTYPE,NGROUP,NSET)
1833 IF (XMINU .GT. ZEROMIN) THEN
1834 XMIN = XMINU
1835 IF(FIRST) THEN
1836 IF(XMINU.NE.XMINP) THEN
1837 WRITE (N6,*) ' PDFLIB : XMIN value'
1838 WRITE(N6,*) ' Warning : NON standard settings, XMIN value = ',
1839 + XMIN,' set by user !!'
1840 ENDIF
1841 ENDIF
1842 ENDIF
1843C.
1844C Maximal X values
1845 XMAX = PXMAX(NPTYPE,NGROUP,NSET)
1846 IF (XMAXU .GT. ZEROMIN) THEN
1847 XMAX = XMAXU
1848 IF(FIRST) THEN
1849 IF(XMAXU.NE.XMAXP) THEN
1850 WRITE (N6,*) ' PDFLIB : XMAX value'
1851 WRITE(N6,*) ' Warning : NON standard settings, XMAX value = ',
1852 + XMAX,' set by user !!'
1853 ENDIF
1854 ENDIF
1855 ENDIF
1856C.
1857C Minimal Q**2 values
1858 Q2MIN = PQ2MIN(NPTYPE,NGROUP,NSET)
1859 IF (Q2MINU .GT. ZEROMIN) THEN
1860 Q2MIN = Q2MINU
1861 IF(FIRST) THEN
1862 IF(Q2MINU.NE.Q2MINP) THEN
1863 WRITE (N6,*) ' PDFLIB : Q2MIN value'
1864 WRITE(N6,*) ' Warning : NON standard settings, Q2MIN value = ',
1865 + Q2MIN,' set by user !!'
1866 ENDIF
1867 ENDIF
1868 ENDIF
1869C.
1870C Maximal Q**2 values
1871 Q2MAX = PQ2MAX(NPTYPE,NGROUP,NSET)
1872 IF (Q2MAXU .GT. ZEROMIN) THEN
1873 Q2MAX = Q2MAXU
1874 IF(FIRST) THEN
1875 IF(Q2MAXU.NE.Q2MAXP) THEN
1876 WRITE (N6,*) ' PDFLIB : Q2MAX value'
1877 WRITE(N6,*) ' Warning : NON standard settings, Q2MAX value = ',
1878 + Q2MAX,' set by user !!'
1879 ENDIF
1880 ENDIF
1881 ENDIF
1882C.
1883C Save the PDFLIB parameters
1884 IFLPRTP = IFLPRT
1885 NPTYPP = NPTYPE
1886 NGROPP = NGROUP
1887 NSETP = NSET
1888 NFLP = NFL
1889 LOP = LO
1890 TMASP = TMAS
1891 QCDL4P = QCDL4
1892 QCDL5P = QCDL5
1893 XMINP = XMIN
1894 XMAXP = XMAX
1895 Q2MINP = Q2MIN
1896 Q2MAXP = Q2MAX
1897 NATYPP = NATYPE
1898 NAGROPP = NAGROUP
1899 NASETP = NASET
1900C.
1901C... Reset variables to zero for error counting
1902 PDFWGT = ONED
1903 WXMIN = ZEROD
1904 WXMAX = ZEROD
1905 WQ2MIN = ZEROD
1906 WQ2MAX = ZEROD
1907 WTXMIN = ZEROD
1908 WTXMAX = ZEROD
1909 WTQ2MIN = ZEROD
1910 WTQ2MAX = ZEROD
1911C.
1912C print COMMON block values
1913 IF(IFLPRT.GE.2) THEN
1914 WRITE(N6,3000) NPTYPE,NGROUP,NSET,SFNAME(NPTYPE,NGROUP,NSET),
1915 + MODECR(NPTYPE,NGROUP,NSET)
1916 WRITE(N6,3001) NFL,LO,TMAS
1917 WRITE(N6,3002) QCDL4,QCDL5
1918 WRITE(N6,3003) XMIN,XMAX
1919 WRITE(N6,3004) Q2MIN,Q2MAX
1920 IF(NAFLAG.EQ.1) THEN
1921 WRITE(N6,3005) NATYPE,NAGROUP,NASET,
1922 + SFNAME(NATYPE,NAGROUP,NASET)
1923 ENDIF
1924 ENDIF
1925 3000 FORMAT(/,' Nptype = ',I1,' Ngroup = ',I2,' Nset = ',I3,
1926 + ' Name = "',A8,'" CrMode = ',I3)
1927 3001 FORMAT(1H ,' Nfl = ',I2,', LO = ',I1,
1928 + ', Tmas = ',F7.2,' GeV/c**2')
1929 3002 FORMAT(1H ,' QCDL4 = ',F7.4,' GeV, QCDL5 = ',F7.4,' GeV')
1930 3003 FORMAT(1H ,' Xmin = ',1X,E8.2,', Xmax = ',E11.5)
1931 3004 FORMAT(1H ,' Q2min =',F7.3,' (GeV/c)**2, Q2max = ',E8.2,
1932 + ' (GeV/c)**2')
1933 3005 FORMAT(/,' NAtype = ',I1,' NAgroup = ',I2,' NAset = ',I3,
1934 + ' Name = "',A8,'"')
1935C.
1936C. Reset PDFSET parameters
1937 DO K = 1,NCHDIM
1938 PARM(K) = ' '
1939 ENDDO
1940C.
1941 RETURN
1942 END