git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / PDF / dpdf / pdfdoc.doc
CommitLineData
fe4da5cc 11***********************************************************************
2 * *
3 * CERN *
4 * *
5 * W5051 PDFLIB 1997.07.02 *
6 * *
7 * --------------------------------------------------------------- *
8 * *
9 * PDFLIB: Nucleon, Pion and Photon Parton Density Functions *
10 * and alpha(s) Calculations *
11 * *
12 * User's Manual *
13 * *
14 * Version 7.09 *
15 * *
16 * *
17 * H. Plothow-Besch/CERN-PPE *
18 * *
19 * CH-1211 Geneva 23 *
20 * E-mail address: plothow@mail.cern.ch *
21 * Tel.: +41 - 22 - 767 39 80 *
22 * *
23 * A complete manual in PostScript format exists on a separate file *
24 * Requests can be addressed to the CERN Program Library Office *
25 * *
26 * Copyright by H. Plothow-Besch *
27 * *
28 * Abstract *
29 * *
30 * This document describes an integrated package of Parton Den- *
31 * sity Functions called PDFLIB which has been added to the CERN *
32 * Program Library Pool W999 and is labelled as W5051. In this pack- *
33 * age all the different sets of parton density functions of the *
34 * NUCLEON, the PION and the PHOTON which are available today (about *
35 * 100 different sets) have been put together. All these sets have *
36 * been combined in a consistent way such that they all have similar *
37 * calling sequences and no external data files have to be read in *
38 * anymore. A default set has been prepared, although those prefering *
39 * their own set or wanting to test a new one may do so within the *
40 * package. The selection mode, the program parameters as well as *
41 * the possibilities to modify the defaults and to control errors *
42 * occured during execution are described. *
43 * *
44 * The package offers as well a program to calculate the strong *
45 * coupling constant alpha(s) to first or second order (default). *
46 * The correct Lambda(QCD) associated to the selected set of struc- *
47 * ture functions and the number of allowed flavours with respect to *
48 * the given Q**2 is automatically used in the caculation. For those *
49 * structure functions where the evolution has been performed to *
50 * leading order only, an alpha(s) value fixed to first order is re- *
51 * turned. *
52 * *
53 * The package is written in standard FORTRAN-77 and is avai- *
54 * lable on all CERN system areas. *
55 * *
56 * *
57 * *
58 * Submitter: H. Plothow-Besch/CERN-PPE *
59 * Language: FORTRAN 77 *
60 * Library: PDFLIB *
61 * *
62 * *
63 * *
64 * *
65 * *
66 * *
67 * *
68 * *
69 ***********************************************************************
701***********************************************************************
71 * *
72 * *
73 * _______________________________________________________________ *
74 * | | *
75 * | PDFLIB - User's Manual | *
76 * | | *
77 * | CERN Computer Program Library entry W5051 | *
78 * | | *
79 * | Copyright CERN, Geneva 1990, 1991, 1992, 1993, 1994, 1995 | *
80 * | Copyright CERN, Geneva 1996 | *
81 * | | *
82 * | Copyright and any other appropriate legal protection of | *
83 * | these computer programs and associated documentation re- | *
84 * | served in all countries of the world. | *
85 * | | *
86 * | These programs or documentation may not be reproduced | *
87 * | and/or redistributed by any method without prior written | *
88 * | consent of the Director-General of CERN or his delegate. | *
89 * | | *
90 * | Permission for the scientific usage of any programs des- | *
91 * | cribed herein is granted apriori to those scientific in- | *
92 * | stitutes associated with the CERN experimental program or | *
93 * | with whom CERN has concluded a scientific collaboration | *
94 * | agreement. | *
95 * | | *
96 * | Commercial utilisation requires explicit a priori per- | *
97 * | mission from CERN and will be subjected to payment of a | *
98 * | licence fee. | *
99 * | | *
100 * | These statements are without prejudice to any rights that | *
101 * | third parties may hold in these programs. | *
102 * | | *
103 * | | *
104 * | Requests for information should be addressed to: | *
105 * | | *
106 * | CERN Computer Program Library Office | *
107 * | CERN-CN Division | *
108 * | CH-1211 Geneva 23 | *
109 * | Switzerland | *
110 * | Tel. +41 22 767 4951 | *
111 * | Fax. +41 22 767 7155 | *
112 * | Electronic Mail address: | *
113 * | | *
114 * | Internet: cernlib@cern.ch | *
115 * |_____________________________________________________________| *
116 * *
117 * *
118 * *
119 * *
120 * *
121 * *
122 * *
123 * *
124 * *
125 * *
126 * *
127 * *
128 * *
129 * *
130 * *
131 * *
132 * *
133 * *
134 * *
135 * *
136 * *
137 * *
138 ******************************** - 1 - ********************************
1391***********************************************************************
140 * *
141 * Table of Contents *
142 * page *
143 * *
144 * PROGRAM SUMMARY .......................................... 3 *
145 * *
146 * PDFLIB : New version 7.09 ................................ 4 *
147 * *
148 * 0. Major Differences compared to earlier Versions ........... 4 *
149 * 0.1 The actual Version 7.09 wrt Version 6.00 ............ 4 *
150 * 0.2 Version 6.00 wrt Version 5.00 ....................... 5 *
151 * 0.3 Version 5.00 wrt Version 4.00 ....................... 6 *
152 * 0.4 Version 4.00 wrt Version 3.00 ....................... 6 *
153 * 0.5 Version 3.00 wrt Version 2.00 ....................... 6 *
154 * 0.6 Version 2.00 wrt Version 1.00 ....................... 7 *
155 * *
156 * 1. The Library of Parton Densities: PDFLIB ................. 10 *
157 * 1.1 Introduction ....................................... 10 *
158 * 1.2 Structure Functions of the Nucleon ................. 12 *
159 * 1.3 Structure Functions of the Pion .................... 12 *
160 * 1.4 Structure Functions of the Photon .................. 12 *
161 * 1.3 The Calling Sequence: PDG Format ................... 14 *
162 * 1.4 The Library PDFLIB ................................. 14 *
163 * 1.5 Calculations of the Strong Coupling Constant alpha(s)15 *
164 * *
165 * 2. Applications of the Parton Densities Package PDFLIB ..... 10 *
166 * 2.1 Parameter Setting and Selection of Sets ............ 10 *
167 * 2.2 The Calling Sequence: PDFLIB Format ................ 12 *
168 * 2.3 The Calling Sequence: PDG Format ................... 14 *
169 * 2.4 Error Handling ..................................... 14 *
170 * 2.5 Features ........................................... 15 *
171 * 2.6 How to find PDFLIB ................................. 18 *
172 * *
173 * 3. PION and PHOTON Parton Densities 18 *
174 * *
175 * 4. alpha(s) Calculation 18 *
176 * *
177 * 5. Conclusions 19 *
178 * *
179 * 6. Acknowledgements 21 *
180 * *
181 * 7. References 32 *
182 * *
183 * Appendix 22 *
184 * A : Example of how to use PDFLIB ........................ 22 *
185 * B : List of subprograms and COMMON blocks in PDFLIB *
186 * which are relevant to the user ...................... 25 *
187 * C : Complete list of subprograms in PDFLIB .............. 26 *
188 * *
189 * List of Tables 34 *
190 * *
191 * *
192 * *
193 * *
194 * *
195 * *
196 * *
197 * *
198 * *
199 * *
200 * *
201 * *
202 * *
203 * *
204 * *
205 * *
206 * *
207 ******************************** - 2 - ********************************
2081***********************************************************************
209 * *
210 * PROGRAM SUMMARY *
211 * *
212 * *
213 * Titel of Program : PDFLIB, version 7.09 *
214 * *
215 * Catalogue number : W5051 in the CERN Computer Program Library *
216 * *
217 * Program obtainable from : CERN Computer Program Library Office *
218 * *
219 * Computer for which the program is designed and others on which it *
220 * has been tested: *
221 * Computer : IBM 3090, VAX, IBM-RISC, DECS, SUN, APOLLO, HP-UX, *
222 * Silicon Graphics, NEXT, CDC, CRAY, IBM-PC *
223 * *
224 * Operating systems : VM/CMS, VAX/VMS, Unix, Windows NT *
225 * *
226 * Programming language used : FORTRAN 77, program is available in *
227 * the CVS format *
228 * *
229 * Number of lines in distributed progam, *
230 * including test data, etc : 192 000 *
231 * *
232 * Memory required, *
233 * including testprograms and some documentation : 25 Mbytes *
234 * (Source code: 10 Mbytes, Library: 10 Mbytes) *
235 * *
236 * No. of bits in a word : 64 *
237 * *
238 * Peripherals used : line printer (optional) *
239 * *
240 * Keywords : parton density functions, structure functions, W5051, *
241 * nucleon-, proton-, pion-, photon density parametrisa- *
242 * tions, Altarelli-Parisi equation, strong coupling *
243 * constant alpha(s) *
244 * *
245 * Nature of physical problem : *
246 * All theoretical calculations of cross sections involving the *
247 * hadronic structure of the incoming particles, i.e. protons, pions *
248 * or photons, use the parton density at a given longitudinal momen- *
249 * tum X of the parton and at a given momentum transfer Q**2. There *
250 * is a variety of parametrisations available, but each of them has *
251 * its own structure for input and/or output. *
252 * *
253 * Method of solution : *
254 * Collection of all available parton density parametrisations and *
255 * application in an identical format for all of them, regardless of *
256 * the author's convention or of the particle type. *
257 * *
258 * References : *
259 * [1] H. Plothow-Besch, 'PDFLIB: Structure Functions and alpha(s) *
260 * Calculations', User's Manual - Version 1.00, W5051 PDFLIB, *
261 * 1991.03.21, CERN-PPE. *
262 * [2] H. Plothow-Besch, 'PDFLIB: a library of all available parton *
263 * density functions of the nucleon, the pion and the photon and *
264 * the corresponding alpha(s) calculations', *
265 * Comp. Phys. Comm. 75 (1993) 396-416. *
266 * [3] H. Plothow-Besch, 'The Parton Distribution Function Library', *
267 * Int. J. Mod. Phys. A10 (1995) 2901-2910. *
268 * *
269 * *
270 * *
271 * *
272 * *
273 * *
274 * *
275 * *
276 ******************************** - 3 - ********************************
2771***********************************************************************
278 * *
279 * 0. Major Differences compared to earlier Versions *
280 * *
281 * *
282 * 0.1 The actual Version 7.09 wrt Version 6.00: *
283 * *
284 * *
285 * The following major changes compared to version 6.00 of PDFLIB *
286 * have been implemented: *
287 * *
288 * Several NUCLEON structure function sets have been added: *
289 * *
290 * - the NLL sets A of Martin Roberts and Stirling (MRS-A lQ) for *
291 * low Q**2 of May 1995 in the DIS and the MS_bar renormalisation *
292 * schemes, as well as the NLL sets A' for 6 different values of *
293 * alpha_s (MRS-Ap-150, MRS-Ap-201, MRS-Ap-266, MRS-Ap-344, *
294 * MRS-Ap-435, MRS-Ap-542), the J and J' sets (MRS-J and MRS-Jp) *
295 * of May 1996 which include the high p_T-jet data of the TEVA- *
296 * TRON, and the latest MRS sets R1 to R4 (MRS-R1, MRS-R2, *
297 * MRS-R3, MRS-R4) of May 1996 with grid data as well as in its *
298 * parametrized form (MRS-R1-F) and (MRS-R2-F); *
299 * - version 4 of the CTEQ group of June 1996 as NLL parametrisa- *
300 * tion in the $\overline{MS}$ as well as in the DIS renormalisa- *
301 * tion scheme and as LO (CTEQ4M, CTEQ4D, CTEQ4L), the CTEQ4M set *
302 * for 5 different values of alpha_s (CTEQ4A1, CTEQ4A2, CTEQ4A3, *
303 * CTEQ4A4, CTEQ4A5), a set which includes the high $p_T$-jet *
304 * data of the TEVATRON (CTEQ4HJ), and a set for low Q**2 *
305 * (CTEQ4lQ); *
306 * - the new sets of Gl"uck, Reya, Stratmann and Vogelsang as LO *
307 * and of August 1995 as NLL in the MS_bar renormalisation scheme *
308 * of polarized parton distributions in a 'standard' and a *
309 * 'valence' fit version (GRSVpolS, GRSVpolV). A special COMMON *
310 * block (COMMON /GRSVA1/ A1P, A1N, A1D) transmits directly the *
311 * spin asymmetries of the proton, the neutron and the deuteron. *
312 * *
313 * A NEW GROUP of authors of PHOTON structure function sets has been *
314 * added: *
315 * *
316 * - the 4 old LO sets of Schuler and Sjostrand of March 1995, as *
317 * well as the 4 LO sets of their revised version of May 1996 *
318 * (SaS-G-1D, SaS-G-1M, SaS-G-2D, SaS-G-2MD). *
319 * *
320 * Several new PHOTON structure function sets have been added: *
321 * *
322 * - the recent LO set of Gl"uck, Reya and Stratmann of February *
323 * 1996 (GRS-G); *
324 * - the recent LO and HO sets of Gordon and Storrow of August 1996 *
325 * (GS-G-LO, GS-G-HO). *
326 * *
327 * To take into account the need for more input parameters for parton *
328 * densities of the photon a NEW SUBROUTINE STRUCTP has been added *
329 * which should be called instead of STRUCTM in case of applications *
330 * of the photon. *
331 * *
332 * Several CORRECTIONS to existing structure function sets have been *
333 * performed: *
334 * *
335 * - to account for heavy quark contributions (c and b) in the *
336 * nucleon parton densities of Gl"uck, Reya and Vogt (GRV94), *
337 * now automatically the OLD GRV parton densities of 1992 are *
338 * called for c and b quarks as advised in their paper; *
339 * - corrections for LINNUX-systems suggested by V. Balashov have *
340 * been taken into account. *
341 * *
342 * A BUG in the data tables for the parton densities of the photon of *
343 * the set of Aurenche, Fontannaz and Guillet (AFG-G) has been fixed. *
344 * *
345 ******************************** - 4 - ********************************
3461***********************************************************************
347 * *
348 * More details can be found on the next pages of this updated *
349 * User's Manual. *
350 * *
351 * The response from the Physics community has been very satis- *
352 * factory, which is encouraging us to keep up-to-date with the *
353 * latest developments in this area. The author will be pleased to *
354 * learn about new parton density functions, as well as to receive *
355 * suggestions to improve both, the usage and the documentation. *
356 * Please send your comments to plothow@mail.cern.ch. *
357 * *
358 * *
359 * *
360 * *
361 * 0.2 Version 6.00 wrt Version 5.00: *
362 * *
363 * *
364 * The following major changes compared to version 5.00 of PDFLIB *
365 * have been implemented: *
366 * *
367 * Several NUCLEON structure function sets have been added: *
368 * *
369 * - the NLL sets A' and G of Martin, Roberts and Stirling (MRS-Ap) *
370 * and (MRS-G) of February 1995 with grid data as well as in its *
371 * parametrized version (MRS-Ap-F) and (MRS-G-Fit), *
372 * - the version 3 of the CTEQ group of October 1994 as NLL parame- *
373 * trisation in the MS_bar as well as in the DIS renormalisation *
374 * scheme and as LO (CTEQ3M, CTEQ3D, CTEQ3L), *
375 * - the new sets of GRV of November 1994 as NLL in the MS_bar as *
376 * well as in the DIS renormalisation scheme and as LO GRV *
377 * (GRV94-MS, GRV94DI, GRV94-LO). *
378 * *
379 * A NEW group of PHOTON structure function sets have been added: *
380 * *
381 * - the six NLL sets of the WHIT group (WHIT1-G to WHIT6-G). *
382 * *
383 * Please note that as new the default set the MRS set (G) has been *
384 * chosen (NPTYPE = 1, NGROUP = 3, NSET = 41). *
385 * *
386 * *
387 * *
388 * *
389 * 0.3 Version 5.00 wrt Version 4.00: *
390 * *
391 * *
392 * The following major changes compared to version 4.00 of PDFLIB *
393 * have been implemented: *
394 * *
395 * NUCLEON structure function sets have been added: *
396 * *
397 * - the NLL set A of Martin, Roberts and Stirling (MRS-A) and the *
398 * parametrized version (MRS-A-Fit) *
399 * PHOTON structure function sets have been added: *
400 * - the NLL set of Aurenche, Fontannaz and Guillet (AFG-G). *
401 * *
402 * *
403 * *
404 * *
405 * 0.4 Version 4.00 wrt Version 3.00: *
406 * *
407 * *
408 * The following major changes compared to version 3.00 of PDFLIB *
409 * have been implemented: *
410 * *
411 * NUCLEON structure function sets have been updated or added: *
412 * *
413 * *
414 ******************************** - 5 - ********************************
4151***********************************************************************
416 * *
417 * - the updated NLL sets of Martin, Roberts and Stirling (MRS) of *
418 * November 1992, the sets S0', D0' and D-' in the MSbar as well *
419 * as the sets S0', D0' and D-' in the DIS renormalisation *
420 * scheme. Please note that the previous S0, D0 and D- sets have *
421 * been retracted by the authors. The new MRS(H) set in both *
422 * schemes have also been added, *
423 * - the NLL and LO sets of the CTEQ collaboration (CTEQ) in its *
424 * parametrised form of January 1993. These are the sets 1M, 1MS *
425 * and 1ML in the MSbar, the set 1D in the DIS renormalisation *
426 * scheme and the LO set 1L, *
427 * - the new NLL sets A and B of Berger and Meng (BM) of *
428 * February 1993 with special gluon distributions, and finally *
429 * - the new NLL and LO sets of the CTEQ collaboration (CTEQ) in *
430 * its parametrised form of August 1993. These are the sets 2M, *
431 * 2MS, 2MF and 2ML in the MSbar, the set 2D in the DIS renorma- *
432 * lisation scheme and the LO set 2L. The revised CTEQ2 version, *
433 * CTEQ2p, of November 1993 is recommended by the authors. *
434 * *
435 * Please note that the default set has been put to the GRV set HO *
436 * (NPTYPE = 1, NGROUP = 5, NSET = 3 / former MODE = 72). *
437 * *
438 * The MOST DRASTICAL change of the package has been, however, the *
439 * NEW FORMAT of PARAMETER SETTINGS. Instead of using only ONE para- *
440 * meter (MODE) to select a set of parton density functions (PDF), *
441 * each PDF set is now identified by THREE parameters. These are *
442 * NPTYPE, NGROUP and NSET, for PARTICLE TYPE (Nucleons, Pions or *
443 * Photons) AUTHOR GROUP (f.ex. MRS, CTEQ, GRV, etc) and the PARAME- *
444 * TRISATION SET within the group. Therefore the user should provide *
445 * THREE parameters (NPTYPE, NGROUP, NSET) to the subroutine PDFSET, *
446 * each time a PDF set shall be selected. But to make life nice and *
447 * easy, it is also possible to select a PDF set by transferring only *
448 * ONE parameter to PDFSET as before: this parameter should then *
449 * identify the NAME of the author group in an unequivocal way; the *
450 * selection of the PDF set within the author group is made by the *
451 * corresponding VALUE value (f.ex. PARM(1) = 'MRS' with VALUE(1) = *
452 * 29.0D0 would select the 29th set of the Nucleon PDFs of MRS, and *
453 * f.ex. PARM(1) = 'GRVph' with VALUE(1) = 1.0D0 would select the *
454 * first set of the Photon PDFs of GRV). As a suggestion, the NAME for*
455 * each PDF set can be found in the NEW COMMON block /W505110/SFNAME, *
456 * where SFNAME is a CHARACTER*8 array with dimension SFNAME(NPTYMX, *
457 * NGRMAX,NSETMX), and NPTYMX, NGRMAX, NSETMX are defined by the *
458 * parameter sequence, +SEQ, W5051P2. In case the user wants to se- *
459 * lect the DEFAULT, it is sufficient to transfer as the ONLY parame- *
460 * ter PARM(1) = 'Nset' with VALUE(1) = 0.0D0. This change of format *
461 * had been necessary to take into account the inflationary growing *
462 * number of PDFs until now and to foresee already the possibility of *
463 * adding new PDFs in a more flexible way. Please note that in case *
464 * of any INCONSISTANCY of parameter settings, the DEFAULT setting is *
465 * always chosen. *
466 * *
467 * WARNING : Please also note that the internal COMMON block /W50511/ *
468 * has been modified to /W50511/ NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS. *
469 * *
470 * A new feature, PDFLIB in quiet mode, has been implemented. If the *
471 * user does not want to see any output from PDFLIB, not even the *
472 * version number, he/she must set the parameter PARM(1)='first' with *
473 * value VAL(1)=0.D0. Under this condition any output is suppressed. *
474 * *
475 * Full backward compatibility with previous versions is ensured !! *
476 * This means that old programs using PDF sets up to version 3.00 of *
477 * PDFLIB should run without ANY change. To ensure this, four new *
478 * COMMON blocks have been added: *
479 * *
480 * - /W50519/ NEWVER : to flag the NEW/OLD version with the *
481 * LOGICAL variable NEWVER, *
482 * *
483 ******************************** - 6 - ********************************
4841***********************************************************************
485 * *
486 * - /W505120/ NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX) : *
487 * to provide the actual maximum number of parametrisation sets *
488 * per particle type and author group, and to flag the retracted *
489 * PDF sets (NSETFL(I,J,K) = 0), *
490 * - /W505121/ NPTYCR(MODEMX),NGROMX(MODEMX),NSETCR(MODEMAX) : *
491 * to provide the cross reference to the NPTYPE, NGROUP and NSET *
492 * value (NEW format) for a given MODE value (OLD format), *
493 * and *
494 * - /W505122/ MODECR(NPTYMX,NGRMAX,NSETMX) : to provide the *
495 * cross reference to the MODE value (OLD version) for a given *
496 * NPTYPE, NGROUP and NSET value (NEW format). In case of the *
497 * above added PDF sets NO MODE value is available; in that case *
498 * MODECR(I,J,K) = -1 is returned. *
499 * *
500 * *
501 * *
502 * *
503 * 0.5 Version 3.00 wrt Version 2.00 *
504 * *
505 * *
506 * The following major changes compared to version 2.00 of PDFLIB *
507 * had been implemented: *
508 * *
509 * NUCLEON structure function sets had been updated or added: *
510 * *
511 * - the NLL set of Aurenche et al. (ABFOW), *
512 * - the updated NLL set of Kwiecinski, Martin, Roberts and *
513 * Stirling (KMRS) set B0 with u_bar not equal to d_bar *
514 * of April 1992, *
515 * - the new NLL sets of Martin, Roberts and Stirling (MRS) *
516 * set S0, D0 and D- with u_bar not equal to d_bar of April 1992. *
517 * *
518 * Please note that the default set had been put to the MRS set S0 *
519 * (MODE = 45). *
520 * *
521 * PION structure function sets had been added: *
522 * *
523 * - the new LO and NLL sets of Gl"uck, Reya and Vogt (GRV-P) *
524 * of August 1991, *
525 * - the new sets 2 and 3 of Aurenche et al. (ABFKW-P). *
526 * *
527 * Please also note that the ABFKW set had changed MODE number from *
528 * 110 -> 180, as well as the SMRS sets from 120ff -> 130ff, and *
529 * that earlier inconsistencies in the parton density definitions *
530 * had been corrected. *
531 * *
532 * For the first time also PHOTON structure function sets had been *
533 * added to the package. These sets were the following: *
534 * *
535 * - the old Duke and Owens (DO-G) LO and NLL sets, *
536 * - the old LO sets 1 to 4 of Drees and Grassie (DG-G), *
537 * - the new LO sets 1 to 3 of Abramowicz, Charchula and Levy *
538 * (LAC-G), *
539 * - the new LO sets 1 and 2 and the NLL set of Gordon and Storrow *
540 * (GS-G) of July 1991, *
541 * - the new LO and NLL sets, as well as the leading term of the *
542 * NLL calculation of Gl"uck, Reya and Vogt (GRV-G) *
543 * of December 1991, *
544 * - the new NLL sets of Aurenche et al. (ACFGP-G) of April 1992. *
545 * *
546 * Some improvements in the usage of the package had been made: *
547 * *
548 * - more protections had been implemented and inconsistencies *
549 * had been corrected to make the package more robust and *
550 * available on ALL system areas, and finally *
551 * *
552 ******************************** - 7 - ********************************
5531***********************************************************************
554 * *
555 * - the subroutine STRUCTM had been added. STRUCTM is basically *
556 * the same routine as the old steering subroutine STRUCTF apart *
557 * from ONE additional parameter in the argument list which had *
558 * been added to transfer separately the parton densities for *
559 * u_bar and d_bar. The argument list of the subroutine STRUCTF *
560 * was as before; in case of STRUCTF sea = (u_bar + d_bar)/2 is *
561 * transferred. The subroutine PFTOPDG, which is the sterring *
562 * routine in PDG format, was unchanged; this routine always *
563 * transferred u_bar and d_bar separately. *
564 * *
565 * *
566 * *
567 * *
568 * 0.6 Version 2.00 wrt Version 1.00 *
569 * *
570 * *
571 * The following major changes compared to version 1.00 of PDFLIB *
572 * had been implemented: *
573 * *
574 * New parton density functions of the NUCLEON had been added: *
575 * *
576 * - the old LO set of Buras and Gaemers of 1978 (BEBC) for com- *
577 * pleteness, *
578 * - all the LO sets of Diemoz, Ferroni, Longo and Martinelli (DFLM)*
579 * including their parametrisation for the NLL average fit, *
580 * - the new LO set 1.1 of Owens (DO 1.1) of June 1991, *
581 * and finally *
582 * - the new LO and NLL sets of Gl"uck, Reya and Vogt (GRV) of *
583 * May 1991. *
584 * *
585 * PION structure function sets had been added to the package: *
586 * *
587 * - the old Owens (OW-P) sets 1 and 2, *
588 * - the newer Aurenche et al. (ABFKW-P) set and *
589 * - the new Sutton, Martin, Roberts and Stirling (SMRS-P) sets 1, *
590 * 2 and 3. *
591 * *
592 * Some improvements in the usage of the package had been made: *
593 * *
594 * - COMMON block names had been changed to avoid interference *
595 * to those defined by the user (/IMODE/ -> /W50511/, *
596 * /ALAMB/ -> /W50512/), *
597 * - three new COMMON blocks had been added: *
598 * a) /W50513/ XMIN,XMAX,Q2MIN,Q2MAX : *
599 * to provide the minimum and maximum values of X and Q**2 for *
600 * each set of structure functions, *
601 * b) /W50510/ IFLPRT : *
602 * to provide a print flag at different levels to print *
603 * b.1) the variables of the three COMMON blocks *
604 * /W50511/, /W50512/ and /W50513/ at job initialisation *
605 * b.2) error messages during job execution *
606 * c) /W50514W/ PDFWGT : *
607 * to take into account the event weight for error counting *
608 * - a subroutine PDFSTA had been added to allow at job termination *
609 * for a summary how often X and Q**2 limits have been exceeded, *
610 * - the MODE number of some structure functions had been changed *
611 * to keep consistensy, *
612 * - the number of parameters which can be set through the sub- *
613 * routine PDFSET had been extended to 20, *
614 * - more protections had been implemented and inconsistencies *
615 * had been corrected to make the package more robust, *
616 * and finally *
617 * - a subroutine PFTOPDG had been provided as an interface for *
618 * those users who prefer the flavour code convention of the PDG *
619 * instead of the PDFLIB format. *
620 * *
621 ******************************** - 8 - ********************************
6221***********************************************************************
623 * *
624 * 1. The Library of Parton Densities: PDFLIB *
625 * *
626 * *
627 * 1.1. Introduction *
628 * *
629 * *
630 * Perturbative QCD describes the $Q^2$ dependence of the parton *
631 * densities through the solution of the Altarelli-Parisi equations. *
632 * Its predictive power is given by the universality of these densi- *
633 * ties. Usually these densities are extracted from precision mea- *
634 * surements in fixed-target experiments of a given process at some *
635 * scale - recently they have also been extracted from the HERA ep- *
636 * collision experiments -, and then used to perform calculations for *
637 * different reactions in an extrapolated range of $Q^2$ provided *
638 * that the elementary cross sections at the parton level are known. *
639 * To be able to predict cross sections for a wide range of hard *
640 * scattering processes, it is of fundamental importance to have a *
641 * reliable and precise set of parton density distributions, together *
642 * with a precise value of the strong coupling constant alpha_s. *
643 * Therefore the momentum distributions of quarks and gluons are key *
644 * ingredients for the lepton-lepton, lepton-hadron and hadron-hadron *
645 * collider phenomenology, whenever the hadronic structure of nucle- *
646 * ons, pions or photons is involved. Apart from the obvious applica- *
647 * tions mentioned above, the study of parton distributions provides *
648 * also information about flavour dependence, isospin symmetry, par- *
649 * tons in nuclei, sum rules, polarized parton distibutions, etc.. *
650 * For any application the parton densities should therefore be *
651 * sufficiently well known, and the different sets should be easily *
652 * accessible. *
653 * *
654 * *
655 * *
656 * *
657 * 1.2. Structure Functions of the Nucleon *
658 * *
659 * *
660 * The number of sets of parton density functions of the NUCLEON *
661 * available today is about 70 and is still growing. The older sets *
662 * are leading order (LO) evolutions of the Altarelli-Parisi equa- *
663 * tions only, like *
664 * *
665 * - Buras and Gaemers (BEBC) [1], *
666 * - Owens and Reya (OR) [2], *
667 * - Baier, Engels and Petersson (BEP) [3], *
668 * - Gl"uck, Hoffmann and Reya (GHR) [4], *
669 * - Duke and Owens (DO) sets 1 and 2 [5], *
670 * - Eichten, Hinchliffe, Lane and Quigg (EHLQ) sets 1 and 2 [6], *
671 * - Diemoz, Ferroni, Longo and Martinelli (DFLM) [7] with *
672 * Lambda_4(QCD) = 200 MeV for soft and hard valence quark and *
673 * soft and hard gluon distributions], *
674 * *
675 * and *
676 * *
677 * - the set 1.1 of Owens (DO 1.1) [8], replacing the DO [5] sets 1 *
678 * and 2. *
679 * *
680 * Most of the sets, essentially the more recent ones, are *
681 * next-to-leading order evolutions, like *
682 * *
683 * - Diemoz, Ferroni, Longo and Martinelli (DFLM) with *
684 * Lambda_4(QCD) = 160, 260, 360 MeV [7], *
685 * - Martin, Roberts and Stirling (MRS) sets 1 to 3, sets E and B *
686 * and sets E' and B' [9], *
687 * - Harriman, Martin, Roberts and Stirling (HMRS) sets E, E+, E- *
688 * and B [10], *
689 * *
690 ******************************** - 9 - ********************************
6911***********************************************************************
692 * *
693 * - Kwiecinski, Martin, Roberts and Stirling (KMRS) [11] sets B0 *
694 * and B-, and for the B- set with different radii for shadowing, *
695 * - Martin, Roberts and Stirling (MRS-B Lambda) [11] for the B0 *
696 * set also with different values of Lambda_4(QCD) = 135, 160, *
697 * 200 and 235 MeV, *
698 * - Martin, Roberts and Stirling (MRS) [12] sets S0, D0 and D- *
699 * as well as the upgraded KMRS set B0, all having u_bar not *
700 * equal to d_bar, *
701 * - Morfin and Tung (MT) sets 1 to 6 [13], *
702 * - Gl"uck, Reya and Vogt (GRV) [14] including their new para- *
703 * metrisation from August 1991, *
704 * - the set of Aurenche et al. (ABFOW) [15], *
705 * - the upgraded Martin, Roberts and Stirling (MRS) [16] sets *
706 * S0', D0' and D-' as well as the more recent MRS(H) set, all in *
707 * both, the MS_bar and the DIS renormalisation schemes, *
708 * - the upgraded Martin, Roberts and Stirling (MRS) [17] sets *
709 * S0', D0' and D-' in the DIS renormalisation scheme and MRS(H), *
710 * - the Berger and Meng (BM) [17] sets 1 and 2 with special gluon *
711 * distributions, *
712 * - the two first versions of the CTEQ collaboration (CTEQ1 and *
713 * CTEQ2) [18], *
714 * - the set of Martin, Roberts and Stirling (MRS-A) [19] in form *
715 * of a grid and in form of a parametrisation, *
716 * - the two sets of Martin, Roberts and Stirling (MRS-A' and *
717 * MRS-G) [20] in form of a grid and in form of a parametrisation,*
718 * - the third versions of the CTEQ collaboration (CTEQ3) [21], *
719 * - the 1994 set of Gl\"{u}ck, Reya and Vogt (GRV) [22], *
720 * - the recent sets of Martin, Roberts and Stirling for low Q^2 *
721 * (MRS-A lQ) [23] and the MRS-A' sets for 6 different values of *
722 * alpha_s (MRS-Ap-150 to MRS-Ap-542) [24], *
723 * - the more recent sets of Martin, Roberts and Stirling (MRS-J *
724 * and MRS-Jp) [25] and the latest sets R1 to R4 (MRS-R1 to *
725 * MRS-R4) [26] in form of a grid and in form of a parametrisa- *
726 * tion, *
727 * - the fourth versions of the CTEQ collaboration (CTEQ4) [27], *
728 * *
729 * and finally *
730 * *
731 * - the LO and NLL sets of Gl\"{u}ck, Reya, Stratmann and *
732 * Vogelsang (GRSV) [28] of polarized parton densities, *
733 * *
734 * The different NLL sets use either the DIS [44] or the MS_bar *
735 * renormalisation scheme. Some of the sets are available in both *
736 * renormalisation schemes. All sets are made using four flavours in *
737 * the initial state (NF = 4). The value of the QCD scale factor, *
738 * Lambda_4(QCD), ranges from 45 to 500 MeV. *
739 * *
740 * A summary of all parton density functions of the NUCLEON avai- *
741 * lable in the package can be found in Table 1. *
742 * *
743 * *
744 * *
745 * *
746 * *
747 * *
748 * *
749 * *
750 * *
751 * *
752 * *
753 * *
754 * *
755 * *
756 * *
757 * *
758 * *
759 ******************************* - 10 - ********************************
7601***********************************************************************
761 * *
762 * 1.3. Structure Functions of the Pion *
763 * *
764 * *
765 * There is also a limited number of parton density functions of *
766 * the PION available. These are the old set to leading order evolu- *
767 * tion of *
768 * *
769 * - Owens (OW-P) [29], sets 1 and 2, *
770 * *
771 * and the next-to-leading order evolutions of *
772 * *
773 * - Aurenche et al. (ABFKW-P) [30], sets 1, 2 and 3, *
774 * - Sutton, Martin, Roberts and Stirling (SMRS-P) [31], sets 1, 2 *
775 * and 3 and *
776 * - Gl"uck, Reya and Vogt (GRV-P) [32], the LO and NLL sets. *
777 * *
778 * A summary of all parton density functions of the PION availble *
779 * in the package can be found in Table 2. *
780 * *
781 * *
782 * *
783 * *
784 * 1.4. Structure Functions of the Photon *
785 * *
786 * *
787 * In addition there are also a few parton density functions of *
788 * the PHOTON available. These are the old sets to leading order evo- *
789 * lution of *
790 * *
791 * - Duke and Owens (DO-G) (and the NLL) [33] using an asymptotic *
792 * solution of the Altarelli-Parisi equation, *
793 * - Drees and Grassie (DG-G) sets 1 to 4 [34] using a full solu- *
794 * tion of the Altarelli-Parisi equation, *
795 * - Abramowicz, Charchula and Levy (LAC-G) [35] sets 1 to 3 using *
796 * a full solution of the Altarelli-Parisi equation, *
797 * *
798 * and of *
799 * *
800 * - Watanabe, Hagiwara, Izubuchi and Tanaka (WHIT1 to WHIT6) [36], *
801 * *
802 * and the next-to-leading order evolution of *
803 * *
804 * - Gordon and Storrow (GS-G) [37], the LO sets 1 and 2 and the *
805 * NLL set, *
806 * - Gl"uck, Reya and Vogt (GRV-G) [38], the LO, NLL and the *
807 * leading term of the NLL sets, *
808 * - Aurenche et al. (ACFGP-G) [39] sets 1 and 2, without and with *
809 * massive charm, *
810 * *
811 * and finally *
812 * *
813 * _ the LO and NLL sets of Aurenche, Fontannaz and Guillet, *
814 * (AFG-G) [40]. *
815 * *
816 * Recently new sets of of the parton densities of the photon have *
817 * been come out: *
818 * *
819 * - Schuler and Sjostrand (SaS-G) [41], *
820 * - the new sets of Gl"uck, Reya and Stratmann (GRG-G) [42], *
821 * and *
822 * - the new sets LO and NLL of Gordon and Storrow (GS-G) [43]. *
823 * *
824 * A summary of all parton density functions of the PHOTON avai- *
825 * lable in the package can be found in Table 3. *
826 * *
827 * *
828 ******************************* - 11 - ********************************
8291***********************************************************************
830 * *
831 * 1.5 The Library PDFLIB *
832 * *
833 * *
834 * We have put together all these different sets of parton density *
835 * functions in ONE SINGLE PACKAGE [45]. This library PDFLIB can be *
836 * found in and obtained from the CERN Computer Program Library under *
837 * the registration number W5051. *
838 * *
839 * Please notice that the use of the PDFLIB library does not pre- *
840 * vent to consult the publications of the different authers, which *
841 * can be found in the bibliography of this manual, for details about *
842 * each parton density function. *
843 * *
844 * We have modified the sets such that NO external data files for *
845 * the grids have to be read in anymore. All these structure function *
846 * sets have been combined in a consistent way such that they all *
847 * appear in an IDENTICAL STRUCTURE to the user. The selection is *
848 * made via THREE parameters, NPTYPE, NGROUP and NSET, which identify *
849 * a PDF set by its PARTICLE TYPE (Nucleon: NPTYPE = 1, *
850 * Pion: NPTYPE = 2, Photon: NPTYPE = 3), its AUTHOR GROUP (f.ex. *
851 * MRS: NGROUP = 3, CTEQ: NGROUP = 4, GRV: NGROUP = 5) and its PDF *
852 * SET WITHIN THE GROUP (NSET = 1,2,3 etc.). These parameters should *
853 * be set with a call to subroutine PDFSET at the initialization *
854 * phase. A simple SUBROUTINE call, which is IDENTICAL for all appli- *
855 * cations (nucleons, pions, photons) returns the parton densities *
856 * for all partons (u, d, s, c, b, t, gluon and their antiquarks) at *
857 * a given X value (where X is the fraction of the longitudinal *
858 * momentum carried by the parton) and the Q-scale SCALE (in GeV). *
859 * A DEFAULT set has been prepared, although those preferring their *
860 * OWN PRIVATE set or wanting to test a NEW one may do so within the *
861 * package. Error control can be obtained by setting a print flag to *
862 * obtain output either during execution and/or as a summary at job *
863 * termination via a call to the subroutine PDFSTA. *
864 * *
865 * The source code of the different parton density sets has been *
866 * modified in such a way that the library can be used on all the *
867 * different computer systems known today. Full backward compatibi- *
868 * lity will always be assured. PDFLIB has been run under different *
869 * operating systems like VM/CMS, VAX/VMS and ULTRIX without pro- *
870 * blems. The library has been tested on the different system areas *
871 * to ensure identical results within the given machine precision. *
872 * The program has been run and tested on the following computers: *
873 * IBM 3090, VAX, IBM-RISC, DECS, SUN, APOLLO, HP-UX, Silicon *
874 * Graphics, CDC, CRAY, IBM-PC. The programming language is *
875 * FORTRAN 77. The actual source code manager system used is CVS. *
876 * The required memory storage including testprograms and some docu- *
877 * mentation is about 25 Mbytes. The latest edition (July 1997) *
878 * is version 7.09 and can be requested from the CERN Computer *
879 * Program Library. To run PDFLIB a link to the CERN library *
880 * (PACKLIB, MATHLIB and KERNLIB) is required. *
881 * *
882 * *
883 * *
884 * 1.6. Calculation of the Stroing Coupling constant alpha(s) *
885 * *
886 * *
887 * Within the SAME package a program is also provided to calculate *
888 * the strong coupling constant alpha(s) to second order (default) *
889 * or to first order (by user's choice). The correct Lambda(QCD) as- *
890 * sociated to the selected set of structure functions and the number *
891 * of allowed flavours with respect to the given Q**2 is automati- *
892 * cally used in the calculation. For those structure functions where *
893 * the evolution has been performed to leading order only, the *
894 * alpha(s) value to first order is returned. *
895 * *
896 * *
897 ******************************* - 12 - ********************************
8981***********************************************************************
899 * *
900 * 2. Application of the Parton Density Package PDFLIB *
901 * *
902 * *
903 * In the following the use of and the access to the PDFLIB pack- *
904 * age is described. *
905 * *
906 * *
907 * *
908 * 2.1 Parameter Setting and Selection of Sets *
909 * *
910 * *
911 * To access his/her preferred set of structure functions the user *
912 * must define the THREE parameters, NPTYPE, NGROUP and NSET, which *
913 * identify each set of parton density functions via a call to the *
914 * subroutine *
915 * *
916 * CALL PDFSET(PARM,VALUE) *
917 * *
918 * at the initialization phase of his/her MAIN program, or use the *
919 * package as a 'black box' with the default values (see below). An *
920 * example of the application of PDFLIB is given in the appendix. The *
921 * arguments PARM and VALUE are VECTORS of DIMENSION 20 and have the *
922 * following meaning *
923 * *
924 * PARM(I) = character*20 variable, which defines in any order *
925 * the variables *
926 * 'NPTYPE', 'NGROUP', 'NSET', 'MODE', 'INIT0', *
927 * 'NFL', 'LO', 'TMAS', 'QCDL4', 'QCDL5' and *
928 * 'XMIN', 'XMAX', 'Q2MIN', 'Q2MAX' *
929 * VALUE(I) = the corresponding numerical value of the variable *
930 * PARM(I) *
931 * (TMAS, QCDL4, QCDL5, XMIN, XMAX, Q2MIN and Q2MAX *
932 * are DOUBLE PRECISION variables), *
933 * *
934 * where *
935 * *
936 * NPTYPE = number of particle type ranging from 1 to 3 *
937 * (Nucleons: NPTYPE = 1, Pions: NPTYPE = 2 and *
938 * Photons: NPTYPE = 3) *
939* (Default: NPTYPE = 1) *
940 * NGROUP = number of author group ranging from 1 to 9 *
941 * (Default: NGROUP = 3) *
942 * NSET = number of a selected structure function set within *
943 * the author group ranging from 1 to 58 *
944 * (Default: NSET = 41) *
945 * (MODE = number of a selected structure function set ranging *
946 * from 0 to 281 - OLD format !!) *
947 * (Default: MODE = 45) *
948 * INIT0 = in case of PARM(1) = 'INIT0' PDFSET fills as the *
949 * only action the COMMON blocks /W505110/, /W505120/ *
950 * /W505121/ and /W505122/ *
951 * NFL = desired number of flavours in the alpha(s) calcula- *
952 * tion ranging from 3 to 6 *
953 * (Default: NFL = 5) *
954 * LO = order of alpha(s) calculation; if LO = 1, alpha(s) *
955 * to first order only *
956 * (Default: LO = 2) *
957 * TMAS = the user defined value of the top-quark mass in *
958 * GeV/c**2 (optional) *
959 * (Default: TMAS = 180.0D0) *
960 * QCDL4 = QCD scale, Lambda_4(QCD), in GeV for four flavours *
961 * QCDL5 = QCD scale, Lambda_5(QCD), in GeV for five flavours *
962 * corresponding to QCDL4 *
963 * *
964 * and *
965 * *
966 ******************************* - 13 - ********************************
9671***********************************************************************
968 * *
969 * XMIN = minimum allowed X value *
970 * XMAX = maximum allowed X value *
971 * Q2MIN = minimum allowed Q**2 value (in (GeV/c)**2) *
972 * Q2MAX = maximum allowed Q**2 value (in (GeV/c)**2) *
973 * *
974 * for each set of structure functions. *
975 * *
976 * Please also note that PDFSET can be called as often as the user *
977 * likes. In order to redefine the parameters to select other sets of *
978 * structure functions which the user wants to investigate, it might *
979 * be necessary to call PDFSET in an alternating way, but it is al- *
980 * ways mandatory to transfer either the THREE parameters, NPTYPE, *
981 * NGROUP, NSET, or the 'NAME' of the author group with their corres- *
982 * ponding VALUE values. The subroutine PDFSET fills the internal *
983 * COMMON blocks *
984 * *
985 * COMMON/W50511/ NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS *
986 * COMMON/W50512/ QCDL4,QCDL5 *
987 * COMMON/W50513/ XMIN,XMAX,Q2MIN,Q2MAX *
988 * *
989 * at the time the routine is called. Please note that in case of a *
990 * multiple call to PDFSET with fewer parameters redefined than in a *
991 * preceeding call, always the last parameters are kept in memory. *
992 * All variables of the three COMMON blocks can be automatically *
993 * printed at job initialisation by setting the print flag IFLPRT in *
994 * the COMMON block *
995 * *
996 * COMMON/W50510/ IFLPRT *
997 * *
998 * to IFLPRT = 2. *
999 * *
1000 * The program is protected against calculations of Q^2 values *
1001 * below Q2MIN (in this case Q^2 = Q2MIN) and of X values in unphysi- *
1002 * cal regions (X < 0 or X > 1). The execution of the program is *
1003 * stopped in the later case. To control how often the X or Q^2 ran- *
1004 * ges have been exceeded during execution of the user's program a *
1005 * call to the subroutine *
1006 * *
1007 * CALL PDFSTA *
1008 * *
1009 * at the termination phase of his/her MAIN program allows to print a *
1010 * summary of these errors. The print flag in the COMMON/W50510/IFLPRT*
1011 * may be set to IFLPRT = 3 to print an error message EACH TIME a *
1012 * limit has been exceeded during job execution. Please note that *
1013 * this may produce an ENORMOUS AMOUNT OF OUTPUT ! *
1014 * *
1015 * *
1016 * *
1017 * *
1018 * *
1019 * 2.2 The Calling Sequence: PDFLIB Format *
1020 * *
1021 * *
1022 * The main steering routine for a set of structure functions is *
1023 * accessed as follows: *
1024 * *
1025 * CALL STRUCTM(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL) *
1026 * *
1027 * Please note that all variables are defined as DOUBLE PRECISION. *
1028 * The user has to provide the following INPUTs: *
1029 * *
1030 * X = X value of parton *
1031 * SCALE = QCD scale in GeV *
1032 * *
1033 * The subroutine STRUCTM returns the following OUTPUT: *
1034 * *
1035 ******************************* - 14 - ********************************
10361***********************************************************************
1037 * *
1038 * UPV = up valence quark *
1039 * DNV = down valence quark *
1040 * USEA = sea (up_bar) *
1041 * DSEA = sea (down_bar) *
1042 * STR = strange quark *
1043 * CHM = charm quark *
1044 * BOT = bottom quark *
1045 * TOP = top quark *
1046 * GL = gluon *
1047 * *
1048 * In case up_bar is not given separately from down_bar it is set *
1049 * USEA = DSEA. *
1050 * *
1051 * The recommended set of structure functions from the different *
1052 * authors is put in bold characters. If NSET is set to zero (or if *
1053 * any of the parameters, NPTYPE, NGROUP or NSET, is undefined) ONE *
1054 * DEFAULT SET is selected which is always the Nucleon parton densi- *
1055 * ties, the set of MRS set G (NPTYPE = 1, NGROUP = 3, NSET = 41), *
1056 * also in case of NPTYPE = Pion or Photon. *
1057 * *
1058 * Please note that in any of the calling sequences for the *
1059 * nucleon, the pion and the photon it is always returned *
1060 * X * parton distribution function ! *
1061 * *
1062 * *
1063 * *
1064 * *
1065 * 2.3 The Calling Sequence: PDG Format *
1066 * *
1067 * *
1068 * To allow the use of the package with the flavour code conven- *
1069 * tion of the Particle Data Group (PDG) an interface has been writ- *
1070 * ten which translates the PDFLIB format (see Section 2.2) into the *
1071 * PDG format. Instead of calling the subroutine STRUCTM the user *
1072 * preferring the flavour code convention of the PDG accesses the *
1073 * package via a call to the subroutine *
1074 * *
1075 * CALL PFTOPDG(X,SCALE,DXPDF). *
1076 * *
1077 * All variables are defined as DOUBLE PRECISION as before, and *
1078 * X and SCALE are INPUTs provided by the user having the same *
1079 * meaning as before. *
1080 * *
1081 * The subroutine PFTOPDG outputs a vector DXPDF(-6:6) of which *
1082 * its variables have the following meaning: *
1083 * *
1084 * DXPDF(0) = gluon *
1085 * DXPDF(1) = sum of down valence + down_bar quark *
1086 * DXPDF(2) = sum of up valence + up_bar quark *
1087 * DXPDF(3) = strange quark *
1088 * DXPDF(4) = charm quark *
1089 * DXPDF(5) = bottom quark *
1090 * DXPDF(6) = top quark *
1091 * and DXPDF(-1) to DXPDF(-6) are the corresponding antiquarks. *
1092 * *
1093 * DXPDF(-1) = DXPDF(-2) stands for down_bar = up_bar or *
1094 * (down_bar + up_bar)/2 in almost all parametrisations. In case *
1095 * where u_bar is not equal to d_bar DXPDF(-1) = d_bar and *
1096 * DXPDF(-2) = u_bar. In all sets is DXPDF(3) = DXPDF(-3), *
1097 * DXPDF(4) = DXPDF(-4), DXPDF(5) = DXPDF(-5) and *
1098 * DXPDF(6) = DXPDF(-6) so far. *
1099 * *
1100 * *
1101 * *
1102 * *
1103 * *
1104 ******************************* - 15 - ********************************
11051***********************************************************************
1106 * *
1107 * 2.4 Error Handling *
1108 * *
1109 * *
1110 * A note of caution should be addressed here. All structure *
1111 * functions are limited in X and Q**2 which range for most of the *
1112 * sets from *
1113 * *
1114 * 4 - 10 < Q**2 < 10**6 - 10**8 (GeV/c)**2 *
1115 * 10**{-5} - 10**{-4} < X < 1. *
1116 * *
1117 * If the user wants to extend the X or Q**2 ranges, it is pos- *
1118 * sible for some of the structure function sets, but the result *
1119 * should be looked at with great caution. If in doubt, please check *
1120 * with the authors of the structure function set. The program is *
1121 * protected against calculations of Q**2 values below Q2MIN (in this *
1122 * case Q**2 = Q2MIN) and of X values in unphysical regions (X < 0 or *
1123 * X > 1). The execution of the program is stopped in the later case. *
1124 * *
1125 * A COMMON block can be accessed *
1126 * *
1127 * COMMON/W50513/ XMIN,XMAX,Q2MIN,Q2MAX *
1128 * *
1129 * where *
1130 * *
1131 * XMIN = minimum allowed X value *
1132 * XMAX = maximum allowed X value *
1133 * Q2MIN = minimum allowed Q**2 value (in (GeV/c)**2) *
1134 * Q2MAX = maximum allowed Q**2 value (in (GeV/c)**2) *
1135 * *
1136 * is given for each set of structure functions. To control how often *
1137 * the X or Q**2 ranges have been exceeded during execution of the *
1138 * user's program a call to the subroutine *
1139 * *
1140 * CALL PDFSTA *
1141 * *
1142 * at the termination phase of his/her MAIN program allows to print *
1143 * a summary of these errors. If the COMMON block *
1144 * *
1145 * COMMON/W50514W/ PDFWGT *
1146 * *
1147 * where *
1148 * *
1149 * PDFWGT = weight *
1150 * *
1151 * is filled by the user on an event-to-event basis the statistics *
1152 * in the error summary takes into account this weight. *
1153 * *
1154 * The print flag in the COMMON/W50510/IFLPRT may be set to *
1155 * IFLPRT = 3 to print an error message EACH TIME a limit is exeeded *
1156 * during job execution. Please note that setting this flag may pro- *
1157 * duce an ENORMOUS AMOUNT of output! *
1158 * *
1159 * *
1160 * *
1161 * *
1162 * 2.5 Features *
1163 * *
1164 * *
1165 * If NPTYPE = 1, NGROUP = 1 and NSET = 1, the user has the possi- *
1166 * bility to use his/her own private set of structure functions *
1167 * (i.e. CDHS effective structure functions could be inseted here, *
1168 * or the user could implement a new set of structure functions for *
1169 * testing purposes). To insert the private set the DUMMY subroutine *
1170 * *
1171 * SUBROUTINE STRPRIV(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL) *
1172 * *
1173 ******************************* - 16 - ********************************
11741***********************************************************************
1175 * *
1176 * should be activated by the user. All variables have to be defined *
1177 * as DOUBLE PRECISION. The variables have the same meaning as des- *
1178 * cribed for the subroutine STRUCTM in Section 2.2 (for explanations *
1179 * see there). *
1180 * *
1181 * For the calculation of the strong coupling constant alpha(s) *
1182 * (see Section 3) with the user's private set of structure functions *
1183 * the user should provide a value for the variables QCDL4 and QCDL5 *
1184 * via the subroutine PDFSET at the initialization phase. In fact, *
1185 * because only the parameter QCDL5 is used in the alpha(s) calcul- *
1186 * ation, only this definition is mandatory. *
1187 * *
1188 * *
1189 * *
1190 * *
1191 * 2.6. How to find PDFLIB *
1192 * *
1193 * *
1194 * At CERN the PDFLIB library and corresponding documentation are *
1195 * available on all systems in the standard areas. *
1196 * *
1197 * System Library *
1198 * *
1199 * VXCERN CERN:[NEW.LIB]PDFLIB.OLB *
1200 * Unix systems /cern/new/lib/libpdflib.a *
1201 * *
1202 * *
1203 * At CERN on Unix systems the package may be accessed via the *
1204 * 'cernlib' command: *
1205 * *
1206 * \ 'cernlib pdflib packlib mathlib kernlib' \ *
1207 * *
1208 * All necessary information about PDFLIB can also be found on the *
1209 * World Wide Web via the path 'CERN.Computing.CERNLIB.documentation' *
1210 * or directly as *
1211 * *
1212 * for the flat version : *
1213 * http://wwwcn.cern.ch/asdoc/pdflib/pdflib.html *
1214 * or for the PostScript version : *
1215 * http://wwwcn.cern.ch/asdoc/psdir/pdflib.ps.gz *
1216 * *
1217 * The source code of PDFLIB may be obtained from asisftp.cern.ch *
1218 * via anonymous ftp with: *
1219 * *
1220 * get /cernlib/share/new/src/mclibs/pdf or *
1221 * get /cernlib/share/new/src/mclibs/pdf.tar or *
1222 * get /cernlib/share/new/src/mclibs/pdf.tar.gz *
1223 * *
1224 * even, if the last two files are not visible via the command 'ls'. *
1225 * *
1226 * *
1227 * Whenever you use this program, please give reference to the *
1228 * following papers: *
1229 * *
1230 * H. Plothow-Besch, 'The Parton Distribution Function Library', *
1231 * Int. J. Mod. Phys. A10 (1995) 2901-2920. *
1232 * *
1233 * and/or *
1234 * *
1235 * H. Plothow-Besch, 'PDFLIB: A Library of all available Parton *
1236 * Density Functions of the Nucleon, the Pion and the Photon and *
1237 * the corresponding alpha_s Calculations', *
1238 * Comp. Phys. Comm. 75 (1993) 396-416. *
1239 * *
1240 * and/or *
1241 * *
1242 ******************************* - 17 - ********************************
12431***********************************************************************
1244 * *
1245 * H. Plothow-Besch, 'PDFLIB: Structure Functions and alpha_s *
1246 * Calculations', *
1247 * User's Manual - Version 1.00, W5051 PDFLIB, *
1248 * 1991.03.21, CERN-PPE. *
1249 * *
1250 * and/or *
1251 * *
1252 * H. Plothow-Besch, 'PDFLIB: Nucleon, Pion and Photon Parton *
1253 * Density Functions and alpha_s Calculations', *
1254 * Users's Manual - Version 7.09, W5051 PDFLIB, *
1255 * 1997.07.02, CERN-PPE. *
1256 * *
1257 * *
1258 * *
1259 * *
1260 * 3. PION and PHOTON Parton Densities *
1261 * *
1262 * *
1263 * The calling sequences to the PION and the PHOTON sets of parton *
1264 * densities is kept identical to those described in Sections 2.2 and *
1265 * 2.3 for the nucleon structure functions. *
1266 * *
1267 * In case of the PION structure functions it is set and returned *
1268 * UPV = DNV, SEA = STR (USEA = DSEA = STR) and TOP = 0 *
1269 * (DXPDF(1) = DXPDF(2), DXPDF(3) = DXPDF(-1) = DXPDF(-2) and *
1270 * DXPDF(6) = 0). *
1271 * *
1272 * In case of the PHOTON structure functions it is set and *
1273 * returned UPV = USEA, DNV = DSEA and TOP = 0 *
1274 * (DXPDF(1) = 2 * DXPDF(-1), DXPDF(2) = 2 * DXPDF(-2), as a result *
1275 * of the definitions - see Section 2.3 -, and DXPDF(6) = 0), while *
1276 * all the other quarks are set identical to their antiquark densi- *
1277 * ties. Note that the alpha_(QED) has to be taken care of by the *
1278 * user ! *
1279 * *
1280 * To take into account the need for more input parameters for *
1281 * parton densities of the photon a NEW SUBROUTINE STRUCTP has been *
1282 * added which should be called instead of STRUCTM in case of appli- *
1283 * cations of the photon. Therefore, the main steering routine for a *
1284 * set of structure functions of the photon should now be accessed *
1285 * as follows: *
1286 * *
1287 * CALL STRUCTP(X,Q2,P2,IP2,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL) *
1288 * *
1289 * Please note that all variables (INPUT and OUTPUT) are defined as *
1290 * DOUBLE PRECISION. The user has now to provide the following *
1291 * INPUTs: *
1292 * *
1293 * X = x value of parton *
1294 * Q2 = QCD scale in GeV**2 *
1295 * P2 = virtuality of the photon in GeV**2, *
1296 * should be = 0 for an on-shell photon *
1297 * IP2 = parameter to evaluate off-shell anomalous component *
1298 * *
1299 * The OUTPUT parameters are defined as for STRUCTM (see *
1300 * Section 2.2). *
1301 * *
1302 * *
1303 * *
1304 * *
1305 * *
1306 * *
1307 * *
1308 * *
1309 * *
1310 * *
1311 ******************************* - 18 - ********************************
13121***********************************************************************
1313 * *
1314 * 4. alpha(s) Calculation *
1315 * *
1316 * *
1317 * Within the same package a program is provided to calculate the *
1318 * strong coupling constant alpha(s) to second order as a function *
1319 * of Lambda(QCD) of five flavours and the desired number of fla- *
1320 * vours (NFL) for the selected set of structure functions, which *
1321 * fixes Lambda(QCD). The formula on which the calculations are *
1322 * based upon can be found in Ref. [46]. The same three parameters, *
1323 * NPTYPE, NGROUP and NSET, which select a structure function set, *
1324 * is used to steer the calculation of the alpha(s) value at a given *
1325 * scale from the Lambda(QCD) defined in the selected structure *
1326 * function. *
1327 * *
1328 * The value of alpha(s) is matched at the thresholds q = m_q. *
1329 * When invoked with NFL < 0, it chooses NFL as the number of fla- *
1330 * vours for which the masses are less then q. For the quark masses *
1331 * where thresholds are changed the following values have been used: *
1332 * m_charm = 1.5 GeV/c**2, m_bottom = 4.75 GeV/c**2 and *
1333 * m_top = 180 GeV/c**2. *
1334 * *
1335 * The alpha(s) value can be obtained via a call to: *
1336 * FUNCTION ALPHAS2(SCALE) *
1337 * *
1338 * The user has to provide the following INPUT: *
1339 * *
1340 * SCALE = QCD scale in GeV *
1341 * *
1342 * The function ALPHAS2 returns the following OUTPUT: *
1343 * *
1344 * ALPHAS2 = alpha strong to second order, *
1345 * *
1346 * if LO not equal to one. For those structure functions, for which *
1347 * the evolution is done to leading order only, alpha(s) to first *
1348 * order is returned. *
1349 * *
1350 * The same subroutine *
1351 * *
1352 * SUBROUTINE PDFSET(PARM,VALUE) *
1353 * *
1354 * as described in Section 2.1 should be used to change the default *
1355 * values of the variables NPTYPE, NGROUP, NSET, NFL, LO, TMAS, *
1356 * QCDL4, QCDL5, XMIN, XMAX, Q2MIN and Q2MAX. *
1357 * *
1358 * The internal COMMON block *
1359 * *
1360 * COMMON/W50512/ QCDL4,QCDL5 *
1361 * *
1362 * where *
1363 * *
1364 * QCDL4 = QCD scale, Lambda_4(QCD), in GeV for four flavours *
1365 * QCDL5 = QCD scale, Lambda_5(QCD), in GeV for five flavours *
1366 * corresponding to QCDL4 *
1367 * *
1368 * provides the actual value of Lambda_4(QCD) and Lambda_5(QCD) used *
1369 * in the alpha(s) calculation for four and five flavours, respectiv- *
1370 * ly, for each set of structure functions. *
1371 * *
1372 * *
1373 * *
1374 * *
1375 * *
1376 * *
1377 * *
1378 * *
1379 * *
1380 ******************************* - 19 - ********************************
13811***********************************************************************
1382 * *
1383 * 5. Conclusions *
1384 * *
1385 * *
1386 * In the Long Write-up as an example, the NUCLEON structure *
1387 * function distributions at the scale Q**2 = m_W**2 for up valence *
1388 * lence quark is shown in Figs. 1 as a function of the parton X. *
1389 * The average X values for CERN and FNAL as well as for LHC and SSC *
1390 * energies at that scale are indicated. It can be seen from that *
1391 * figures that there is a large spread in shape for the different *
1392 * sets leading to different results in a cross section calculation. *
1393 * It should be pointed out that the spread from the different para- *
1394 * metrisations in the SAME renormalisation scheme is MUCH LARGER *
1395 * than the spread of parametrisations in different renormalisation *
1396 * schemes. This is valid for ALL parton densities, and in particular *
1397 * at small X values ! *
1398 * *
1399 * A comparison to recent data is shown in Fig. 2, where the ratio *
1400 * F_2(n)/F_2(p) at Q**2 = m_W**2 is displayed. Note that the older *
1401 * sets of structure functions, namely the set of Owens and Reya [2] *
1402 * and the set of Baier et al. [3] have only been implemented for *
1403 * completeness. These sets should not be used for cross section cal- *
1404 * culations of any hard processes at high energies anymore. From *
1405 * Fig. 2 we conclude that also the other older sets of structure *
1406 * functions, namely the two sets of Duke and Owens [5] and the two *
1407 * sets of Eichten et al. [6] should be used with some care, because *
1408 * they do not fit the recent low energy deep-inelastic lepton-nuc- *
1409 * leon data from NMC and BCDMS [47]. This is not surprising because *
1410 * these data were not available when the sets have been made. The *
1411 * new set of DO 1.1 superseeds the old DO sets 1 and 2 but still *
1412 * gives very limited results in the low X region, even though data *
1413 * are available there. We conclude that the recent sets of parton *
1414 * densities could be preferred for all theoretical predictions in- *
1415 * volving structure functions, because they fit the present nucleon *
1416 * data best. *
1417 * *
1418 * In Fig. 3 the F_2(p) distribution of the recent data from Ref. *
1419 * [48] at Q**2 = 5 GeV**2, together with the recent data at Q**2 = *
1420 * 15 GeV**2 from the H1 and the ZEUS experiments [49] at the ep- *
1421 * collider HERA is shown as a function of X. Overlayed on the same *
1422 * Figure are a few theoretical predictions of more recent parton *
1423 * parametrisations [12-16] at Q**2 = 15 GeV**2. Please note that *
1424 * most of these parametrisations are NOT TRUE predictions because *
1425 * the data are fully or partially used in the fits. From Figs. 2 *
1426 * and 3 we conclude that the recent sets of parton densities, namely *
1427 * the GRV set HO [14], but also the old MT set B2 [13] - both sets *
1428 * are true predictions -, and the updated MRS set (H) [16], could *
1429 * be preferred for all theoretical predictions involving structure *
1430 * functions, because they fit the present data best which, at the *
1431 * moment,are still suffering from large uncertainties. *
1432 * *
1433 * Figure 4 shows a comparison of F_2(\gamma)/alpha_(QED) at *
1434 * Q**2 = 100 GeV**2 of the different PHOTON structure function sets *
1435 * as a function of X. The charm quark density has been taken into *
1436 * account. The theoretical predictions are compared with experimen- *
1437 * tal data from the JADE experiment [50] at the same Q**2 value. Be- *
1438 * cause of large experimental uncertainties no distinction between *
1439 * the different parton density sets can be made yet. *
1440 * *
1441 * The response from the Physics community has been very satisfac- *
1442 * tory, which is encouraging us to keep up-to-date with the latest *
1443 * developments in this area. Therefore, please CONTINUE SENDING YOUR *
1444 * FEEDBACK ON THE USAGE AND POSSIBLE IMPROVEMENTS TO US. Coming new *
1445 * sets of structure functions should be easily implemented in the *
1446 * package. Authors of {\bf new sets} are kindly asked to PROVIDE US *
1447 * with the RELEVANT INFORMATION. PLEASE RETURN ANY PROBLEMS, QUES- *
1448 * *
1449 ******************************* - 20 - ********************************
14501***********************************************************************
1451 * *
1452 * TIONS, SUGGESTIONS FOR IMPROVEMENTS TO THE AUTHOR OF THE PACKAGE *
1453 * (e-mail address: plothow@mail.cern.ch). *
1454 * *
1455 * *
1456 * *
1457 * *
1458 * *
1459 * *
1460 * 5. Acknowledgements *
1461 * *
1462 * *
1463 * The author of the package would like to thank *
1464 * M.A. Marquina/CERN-CN for help and support making this package *
1465 * available to the Physics community as part of the CERN Program *
1466 * Library. In addition we give many thanks to G. Folger/CERN-CN and *
1467 * I. McLaren/CERN-CN for thoroughly testing out the package on the *
1468 * different system areas. *
1469 * *
1470 * We also would like to thank J.M. Marraffino/FNAL for the help- *
1471 * ful comments and the correction set which he supplied. *
1472 * *
1473 * And finally we acknowledge K. Charchula/DESY who provided the *
1474 * source code of some of the photon structure functions which had *
1475 * previously been implemented in the package. *
1476 * *
1477 * *
1478 * *
1479 * *
1480 * *
1481 * *
1482 * *
1483 * *
1484 * *
1485 * *
1486 * *
1487 * *
1488 * *
1489 * *
1490 * *
1491 * *
1492 * *
1493 * *
1494 * *
1495 * *
1496 * *
1497 * *
1498 * *
1499 * *
1500 * *
1501 * *
1502 * *
1503 * *
1504 * *
1505 * *
1506 * *
1507 * *
1508 * *
1509 * *
1510 * *
1511 * *
1512 * *
1513 * *
1514 * *
1515 * *
1516 * *
1517 * *
1518 ******************************* - 21 - ********************************
15191***********************************************************************
1520 * *
1521 * References *
1522 * *
1523 * *
1524 * [1] A.J.Buras and K.J.F. Gaemers : Nucl. Phys. B132 (1978) 249; *
1525 * BEBC WA 59 Collaboration, K. Varnell et al., Z. Phys. C36 (1987) 1.
1526 * *
1527 * [2] J.F.Owens and E. Reya : Phys. Rev. D17 (1978) 3003. *
1528 * *
1529 * [3] R. Baier, J. Engels and B. Petersson : Z. Phys. C2 (1979) 265. *
1530 * *
1531 * [4] Gl"uck, E. Hoffman and E. Reya : Z. Phys. C13 (1982) 119. *
1532 * *
1533 * [5] D.W. Duke and J.F.Owens : Phys. Rev. D30 (1984) 49. *
1534 * *
1535 * [6] E. Eichten, I Hinchliffe, K. Lane and C. Quigg : *
1536 * Rev. Mod. Phys. 56 (1984) 579 and Rev. Mod. Phys. 58 (1985) 1065.
1537 * *
1538 * [7] M. Diemoz, F. Ferroni, E. Longo and G. Martinelli : *
1539 * Z. Phys. C39 (1988) 21. *
1540 * *
1541 * [8] J.F.Owens : Phys. Lett. 266B (1991) 126. *
1542 * *
1543 * [9] A.D. Martin, R.G. Roberts and W.J. Stirling : Phys.Rev. D37 (1988) 1161;
1544 * Phys. Lett. 206B (1988) 327 and Mod. Phys. Lett. A4 (1989) 1135.
1545 * *
1546 * [10] P.N. Harriman, A.D. Martin, R.G. Roberts and W.J. Stirling : *
1547 * Phys. Rev. D42 (1990) 798 and Phys. Lett. 243B (1990) 421. *
1548 * *
1549 * [11] J. Kwiecinski, A.D. Martin, R.G. Roberts and W.J. Stirling : *
1550 * Phys. Rev. D42 (1990) 3645; *
1551 * A.D. Martin, R.G. Roberts and W.J. Stirling : *
1552 * Phys. Rev. D43 (1991) 3648. *
1553 * *
1554 * [12] A.D. Martin, R.G. Roberts and W.J. Stirling : *
1555 * Phys. Rev. D47 (1993) 867. *
1556 * *
1557 * [13] J. Morfin and W.K. Tung : Z. Phys. C52 (1991) 13. *
1558 * *
1559 * [14] M. Gl"uck, E. Reya and A. Vogt : Z. Phys. C48 (1990) 471; *
1560 * Z. Phys. C53 (1992) 127 and Phys. Lett. 306B (1993) 391. *
1561 * *
1562 * [15] P. Aurenche et al. : Phys. Rev. D39 (1989) 3275. *
1563 * *
1564 * [16] A.D. Martin, R.G. Roberts and W.J. Stirling : *
1565 * Phys. Lett. 306B (1993) 147 and Erratum-ibid. 309B (1993) 492; *
1566 * Durham Preprint, DTP/93/86 (1993). *
1567 * *
1568 * [17] E.L. Berger and R. Meng : Phys. Lett. 304B (1993) 318; *
1569 * CERN-TH 6739/92; ANL-HEP-CP-92-108 and *
1570 * E.L. Berger, R. Meng and J. Qiu : ANL-HEP-CP-92-79. *
1571 * *
1572 * [18] CTEQ collaboration, J. Botts et al. : Phys. Lett. 304B (1993) 159.
1573 * *
1574 * [19] A.D. Martin, R.G. Roberts and W.J. Stirling : *
1575 * Phys. Rev. D50 (1994) 6734. *
1576 * *
1577 * [20] A.D. Martin, R.G. Roberts and W.J. Stirling : *
1578 * Phys. Lett. 354B (1995) 155; *
1579 * A.D. Martin, R.G. Roberts and W.J. Stirling : *
1580 * RAL Preprint, RAL/95-021 (1995). *
1581 * *
1582 * [21] CTEQ collaboration, MSUHEP-41024. *
1583 * *
1584 * [22] M. Gl"uck, E. Reya and A. Vogt : Z. Phys. C67 (1995) 433. *
1585 * *
1586 * *
1587 ******************************* - 22 - ********************************
15881***********************************************************************
1589 * *
1590 * [23] A.D. Martin, R.G. Roberts and W.J. Stirling : *
1591 * Phys. Rev. D51 (1995) 4756. *
1592 * *
1593 * [24] A.D. Martin, R.G. Roberts and W.J. Stirling : *
1594 * Phys. Lett. 356B (1995) 89. *
1595 * *
1596 * [25] E.W.N. Glover, A.D. Martin, R.G. Roberts and W.J. Stirling : *
1597 * Phys. Lett. 381B (1996) 353. *
1598 * *
1599 * [26] A.D. Martin, R.G. Roberts and W.J. Stirling : *
1600 * Durham Preprint, DTP/96/44 (1996). *
1601 * *
1602 * [27] CTEQ collaboration, MSUHEP-60416, CTEQ-604. *
1603 * *
1604 * [28] M. Gl"uck, E. Reya, M. Stratmann and A. Vogelsang : *
1605 * Phys. Lett. 359B (1995) 201; *
1606 * Phys. Rev. D53 (1995) 4775. *
1607 * *
1608 * [29] J.F.Owens : Phys. Rev. D30 (1984) 943. *
1609 * *
1610 * [30] P. Aurenche et al. : Phys. Lett. 233B (1989) 517. *
1611 * *
1612 * [31] P.J. Sutton, A.D. Martin, R.G. Roberts and W.J. Stirling : *
1613 * Phys. Rev. D45 (1992) 2349. *
1614 * *
1615 * [32] M. Gl"uck, E. Reya and A. Vogt : Z. Phys. C53 (1992) 651. *
1616 * *
1617 * [33] D.W. Duke and J.F.Owens : Phys. Rev. D26 (1982) 1600. *
1618 * *
1619 * [34] M. Drees and K. Grassie : Z. Phys. C28 (1985) 451. *
1620 * *
1621 * [35] H. Abramowicz, K. Charchula and A. Levy : *
1622 * Phys. Lett. 269B (1991) 458. *
1623 * *
1624 * [36] K. Hagiwara et al. : Phys. Rev. D51 (1995) 3197. *
1625 * *
1626 * [37] L.E. Gordon and J.K. Storrow : Z. Phys. C56 (1992) 307. *
1627 * *
1628 * [38] M. Gl"uck, E. Reya and A. Vogt : Phys. Rev. D46 (1992) 1973 *
1629 * and Phys. Rev. D45 (1992) 3986. *
1630 * *
1631 * [39] P. Aurenche et al. : Z. Phys. C56 (1992) 589. *
1632 * *
1633 * [40] P. Aurenche M. Fontannaz, J.Ph. Guillet : Z. Phys. C64 (1994) 621.
1634 * *
1635 * [41] G.A. Schuler and T. Sjostrand : Z. Phys. C68 (1995) 607 and *
1636 * Phys. Lett. 276B (1996) 193. *
1637 * *
1638 * [42] M. Gl"uck, E. Reya, M. Stratmann : Phys. Rev. D51 (1995) 3220. *
1639 * Phys. Rev. D53 (1995) 4775. *
1640 * *
1641 * [43] L.E. Gordon and J.K. Storrow : ANL-HEP-PR-96-33 and MC-TH-96-16.
1642 * *
1643 * [44] G. Altarelli, R.K. Ellis and G. Martinelli : *
1644 * Nucl. Phys. B143 (1978) 521 and B157 (1979) 461. *
1645 * *
1646 * [45] H. Plothow-Besch, 'PDFLIB: Structure Functions and alpha(s) *
1647 * Calculations', User's Manual - Version 1.00, W5051 PDFLIB, *
1648 * 1991.03.21, CERN-PPE. *
1649 * H. Plothow-Besch, 'PDFLIB: a library of all available parton *
1650 * density functions of the nucleon, the pion and the photon *
1651 * and the corresponding alpha(s) calculations', *
1652 * Comp. Phys. Comm. 75 (1993) 396-416; *
1653 * H. Plothow-Besch, 'The Parton Distribution Function Library', *
1654 * Int. J. Mod. Phys. A10 (1995) 2901. *
1655 * *
1656 ******************************* - 23 - ********************************
16571***********************************************************************
1658 * *
1659 * [46] Review of Particle Properties, June 1992, *
1660 * Phys. Rev. D45 (1992), Part. II, p. III.54. *
1661 * *
1662 * [47] BCDMS Collaboration, A.C. Benvenuti et al. : *
1663 * Phys. Lett. 237B (1990) 592 and Phys. Lett. 237B (1990) 599; *
1664 * EMC Collaboration, M. Arneodo et al. : Nucl. Phys. B333 (1990) 1;
1665 * NMC Collaboration, D. Allasia et al. : Phys. Lett. 249B (1990) 366;
1666 * P. Amaudruz et al. : Nucl. Phys. B371 (1992) 3.
1667 * *
1668 * [48] NMC Collaboration, P. Amaudruz et al. : Phys. Lett. 295B (1992) 159.
1669 * *
1670 * [36] H1 Collaboration, I. Abt et al. : Nucl. Phys. B407 (1993) 515; *
1671 * ZEUS Collaboration, M. Derrick et al. : Phys. Lett. 316B (1993) 412.
1672 * *
1673 * [50] JADE Collaboration, W. Bartel et al. : Z. Phys. C24 (1984) 231.*
1674 * *
1675 * [51] PAW, CERN Program Library Entry Q121, Long Write-up; *
1676 * R. Brun, O. Couet, C. Vandoni and P Zanarini; *
1677 * see also to References therein. *
1678 * *
1679 * *
1680 * *
1681 * *
1682 * *
1683 * *
1684 * *
1685 * *
1686 * *
1687 * *
1688 * *
1689 * *
1690 * *
1691 * *
1692 * *
1693 * *
1694 * *
1695 * *
1696 * *
1697 * *
1698 * *
1699 * *
1700 * *
1701 * *
1702 * *
1703 * *
1704 * *
1705 * *
1706 * *
1707 * *
1708 * *
1709 * *
1710 * *
1711 * *
1712 * *
1713 * *
1714 * *
1715 * *
1716 * *
1717 * *
1718 * *
1719 * *
1720 * *
1721 * *
1722 * *
1723 * *
1724 * *
1725 ******************************* - 24 - ********************************
17261***********************************************************************
1727 * *
1728 * Appendix A *
1729 * *
1730 * *
1731 * As an example how to use the package, the FORTRAN code to ob- *
1732 * tain Fig. 1 of the Long Write-up which shows the parton density *
1733 * function distributions for the up valence quark as a function of *
1734 * the parton X for all sets of structure functions is given below: *
1735 * *
1736 * PROGRAM PDFUPV *
1737 * C define the maximum number of PDF sets, the name and the cross reference
1738 * PARAMETER (NPTYMX = 3, NGRMAX = 9, NSETMX = 58) *
1739 * COMMON /W505120/ NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX)
1740 * CHARACTER*8 SFNAME(NPTYMX,NGRMAX,NSETMX) *
1741 * COMMON /W505110/ SFNAME *
1742 * COMMON /W505122/ MODECR(NPTYMX,NGRMAX,NSETMX) *
1743 * C define HBOOK settings *
1744 * PARAMETER (NHBMEM = 500000) *
1745 * COMMON /PAWC/ HMEM(NHBMEM) *
1746 * PARAMETER (NB=900, ID=100) *
1747 * C define DOUBLE PRECISION variables for calling sequence to STRUCTM*
1748 * DOUBLE PRECISION DX,DSCALE *
1749 * DOUBLE PRECISION DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL *
1750 * DOUBLE PRECISION ALF,ALPHAS2 *
1751 * REAL X, SCALE, UPV, DNV, USEA, DSEA, STR, CHM, BOT, TOP, GL *
1752 * COMMON/W50516/ FIRST *
1753 * LOGICAL FIRST *
1754 * CHARACTER*20 PARM(20) *
1755 * DOUBLE PRECISION VAL(20) *
1756 * DATA SCALE/80.140/ *
1757 * DATA X00/0.002/,DX0/0.001/,XLOW/0.0015/,XUP/0.9015/ *
1758 * C *
1759 * CALL HLIMIT(NHBMEM) *
1760 * DSCALE=SCALE *
1761 * C first call to PDFSET to initialize COMMON/W505120/,/W505120/ and /W505122/
1762 * PARM(1) = 'Init0' *
1763 * VAL(1) = 0.D0 *
1764 * CALL PDFSET(PARM,VAL) *
1765 * C loop over all existing sets of Nucleon structure functions (SF) *
1766 * NPTYPE = 1 *
1767 * NHB = 0 *
1768 * DO 20 IGR = 1,NGRMAX *
1769 * IF(NPGSMX(NPTYPE,IGR).EQ.0) GOTO 20 *
1770 * DO ISET = 1,NPGSMX(NPTYPE,IGR) *
1771 * C book histograms for each set of SF separately *
1772 * NHB = NHB + 1 *
1773 * CALL HBOOK1(ID+NHB,'U Valence quark ',NB,XLOW,XUP,0.) *
1774 * C force label printing for each set of SF (not only the 1st) *
1775 * FIRST = .TRUE. *
1776 * C define and set parameters *
1777 * PARM(1) = 'Nptype' *
1778 * VAL(1) = NPTYPE *
1779 * PARM(2) = 'Ngroup' *
1780 * VAL(2) = IGR *
1781 * PARM(3) = 'Nset' *
1782 * VAL(3) = ISET *
1783 * CALL PDFSET(PARM,VAL) *
1784 * C loop over all X bins *
1785 * DO 10 I=1,NB *
1786 * X = X00 + (I-1)*DX0 *
1787 * DX = X *
1788 * IF(X.LT.XLOW .OR. X.GT.XUP) GOTO 10 *
1789 * CALL STRUCTM(DX,DSCALE, *
1790 * + DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DTOP,DGL) *
1791 * UPV=DUPV *
1792 * IF(X.GT.0.499 .AND. X.LE.0.500) WRITE(6,1000) X,SCALE,UPV *
1793 * *
1794 ******************************* - 25 - ********************************
17951***********************************************************************
1796 * *
1797 * Example of Fig. 1, cont'd *
1798 * *
1799 * 1000 FORMAT(/,' X= ',F6.4,' Q= ',F6.3,' UPV= ',8.4) *
1800 * CALL HF1(ID+NHB,X,UPV) *
1801 * 10 CONTINUE *
1802 * C get alpha(s) for selected set of SF at Q = SCALE *
1803 * ALF = ALPHAS2(DSCALE) *
1804 * WRITE(6,4000) NPTYPE,IGR,ISET,SFNAME(NPTYPE,IGR,ISET), *
1805 * + MODECR(NPTYPE,IGR,ISET),NSF *
1806 * WRITE(6,4001) DSCALE,ALF *
1807 * 4000 FORMAT(/,' Nptype = ',I1,' Ngroup = ',I1,' Nset = ',I2, *
1808 * + ' Name = "'A8,'" CrMode = ',I3,' HBId = ',I3) *
1809 * 4001 FORMAT(1H ,' SCALE = ',F8.4,' alpha(s) = ',F6.4) *
1810 * C get error summary for each set of SF *
1811 * CALL PDFSTA *
1812 * ENDDO *
1813 * 20 CONTINUE *
1814 * C *
1815 * STOP *
1816 * END *
1817 * *
1818 * *
1819 * *
1820 * *
1821 * *
1822 * This code will produce the following output: *
1823 * *
1824 * *
1825 * 1. HBOOK Output : *
1826 * *
1827 * one histogram for each set of structure functions. The histo- *
1828 * gram output has then been accessed, and each histogram has been *
1829 * superimposed on one single frame using the PAW package [51] to *
1830 * obtain Fig. 1 of the Long Write-up. *
1831 * *
1832 * *
1833 * 2. Print Output (Example is given for NPTYPE = 1, *
1834 * NGROUP = 3, NSET = 41 only : *
1835 * *
1836 * ___________________________________________________________________ *
1837 *| |*
1838 *| 1***** PDFLIB Version: 7.09 Released on 970702 at 16.05 in |*
1839 *| the CERN Computer Program Library W5051 ***** |*
1840 *| ***** Library compiled on 970702 at 16.05 ***** |*
1841 *| |*
1842 *| Nucleon PDFs: Ngroup = 3, Nset = 41, for MRS Set (G) Structure |*
1843 *| ------------------------------------------------------------------|*
1844 *| Functions |*
1845 *| --------- |*
1846 *| |*
1847 *| X= 0.5000 Q= 80.140 UPV= 0.1176 |*
1848 *| |*
1849 *| Nptype = 1 Ngroup = 3 Nset = 41 Name = "MRS-G " CrMode = -1 HBId = 60
1850 *| Scale = 80.1400 alpha(s) = 0.1154 |*
1851 *| |*
1852 *| PDFLIB : Summary from PDFSTA |*
1853 *| Nptype = 1 Ngroup = 3 Nset = 41 Name = "MRS-G " CrMode = -1|*
1854 *| Nfl = -5 LO = 2 Tmas = 180.00 GeV/c**2 |*
1855 *| QCDL4 = 0.2550 GeV, QCDL5 = 0.1708 GeV |*
1856 *| Xmin = 0.10E-04, Xmax = 0.99999E+00, Q2min = 5.000 (GeV/c)**2,
1857 *| Q2max = 0.27E+10 (GeV/c)**2 |*
1858 *| |*
1859 *| PDFSTA: NO errors occured |*
1860 *|___________________________________________________________________|*
1861 * *
1862 * *
1863 ******************************* - 26 - ********************************
18641***********************************************************************
1865 * *
1866 * Appendix B *
1867 * *
1868 * *
1869 * A list of subroutines, functions, COMMON blocks and parameters *
1870 * of RELEVANCE to the user which are used in PDFLIB is given below: *
1871 * *
1872 * *
1873 * List of relevant subprograms and COMMON blocks in PDFLIB: *
1874 * *
1875 * _________________________________________________________________ *
1876 * | | *
1877 * | (S = Subroutine, F = Function, C = COMMON, P = Parameter) | *
1878 * | | *
1879 * | ------------------------------------------------------------- | *
1880 * | | *
1881 * | S PDFSET : to set all parameters for PDFLIB | *
1882 * | S STRUCTF : to access PDFLIB with output in PDFLIB format | *
1883 * | S STRUCTM : to access PDFLIB with output in PDFLIB format | *
1884 * | (with u_bar not equal to d_bar) | *
1885 * | S STRUCTP : to access PDFLIB with output in PDFLIB format | *
1886 * | for PHOTON structure functions | *
1887 * | S PFTOPDG : to access PDFLIB with output in PDG format | *
1888 * | S PDFSTA : to print summary of parameters and error statistics
1889 * | F ALPHAS2 : to calculate alpha(s) to second order QCD | *
1890 * | | *
1891 * | C W50510 : IFLPRT - to set print flag | *
1892 * | C W50511 : NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS | *
1893 * | - to select a structure function set, to | *
1894 * | define number of flavours, flag for leading | *
1895 * | order and value for the mass of the top quark | *
1896 * | C W50512 : QCDL4,QCDL5 - contains Lambda_4(QCD) and | *
1897 * | Lambda_5(QCD) (in GeV) | *
1898 * | C W50513 : XMIN,XMAX,Q2MIN,Q2MAX - contains minimum and | *
1899 * | maximum values of X and Q**2 | *
1900 * | C W50516 : FIRST - to force printing name of structure | *
1901 * | function set in case of multiple choices | *
1902 * | C W50519 : NEWVER - flags NEW/OLD version of PDFLIB format| *
1903 * | C W505110 : SFNAME - contains NAME (CHARACTER*8 array with | *
1904 * | dimension SFNAME(NPTYMX,NGRMAX,NSETMX)) for | *
1905 * | each set of PDFs | *
1906 * | C W505120 : NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX)*
1907 * | - contains maximum number of stucture functions| *
1908 * | per particle type and author group, and FLAG | *
1909 * | for retracted PDF set (NSETFL(I,J,K) = 0), | *
1910 * | C W505121 : NPTYCR(MODEMX),NGROCR(MODEMX),NSETCR(MODEMX) | *
1911 * | - cross reference to the NPTYPE, NGROUP and | *
1912 * | NSET value (NEW format) for a given MODE | *
1913 * | value (OLD format) | *
1914 * | C W505122 : MODECR(NPTYMX,NGRMAX,NSETMX) - cross reference | *
1915 * | to the MODE value (OLD format) for a given | *
1916 * | NPTYPE, NGROUP and NSET value (NEW format) | *
1917 * | P W5051P2 : MODEMX - maximum number of structure function | *
1918 * | sets in the OLD PDFLIB format | *
1919 * | : NPTYMX,NGRMAX,NSETMX - maximum number of | *
1920 * | structure function sets in the NEW PDFLIB | *
1921 * | format for particle type, author group number | *
1922 * | and number of sets within an author group | *
1923 * | P W5051P7 : L6 - to define logical print unit number | *
1924 * |_______________________________________________________________| *
1925 * *
1926 * *
1927 * *
1928 * *
1929 * *
1930 * *
1931 * *
1932 ******************************* - 27 - ********************************
19331***********************************************************************
1934 * *
1935 * Appendix C *
1936 * *
1937 * *
1938 * A complete list of subroutines, functions, COMMON blocks and *
1939 * parameters which are used in PDFLIB is given below: *
1940 * *
1941 * *
1942 * Complete list of subprograms and COMMON blocks in PDFLIB: *
1943 * (in alphabetical order) *
1944 * *
1945 * _________________________________________________________________ *
1946 * | | *
1947 * | (S = Subroutine, F = Function, C = COMMON, P = Parameter) | *
1948 * | | *
1949 * | ------------------------------------------------------------- | *
1950 * | | *
1951 * | S ABFKW1 : ABFKW SF set 1 in the pion (NP=2,NG=6,NS=1) | *
1952 * | S ABFKW2 : ABFKW SF set 2 in the pion (NP=2,NG=6,NS=2) | *
1953 * | S ABFKW3 : ABFKW SF set 3 in the pion (NP=2,NG=6,NS=3) | *
1954 * | S ACFGP1 : ACFGP SF set 1 in the photon (NP=3,NG=6,NS=1) | *
1955 * | S ACFGP2 : ACFGP SF set 1 in the photon (mc) (NP=3,NG=6,NS=2)*
1956 * | F AFCPLU : to get PD with massive charm for ACFGP | *
1957 * | F AFGETFV : to interpolate PDs in grid for ACFGP | *
1958 * | F AFGIN2 : to interpolate PDs in grid for ACFGP | *
1959 * | F AFGINT : to interpolate PDs in grid for ACFGP | *
1960 * | F AFPOLIN : to interpolate PDs in grid for ACFGP | *
1961 * | F AFRATIN : to interpolate PDs in grid for ACFGP | *
1962 * | F ALPHAS2 : to calculate alpha(s) to second order QCD | *
1963 * | F AUGETFV : to interpolate logarithmically PDs in grid | *
1964 * | S AURGAM : to get the PDs in the photon for ACFGP-G | *
1965 * | S AURPI1 : to get the PDs in the pion for ABFKW-P set 1 | *
1966 * | S AURPI2 : to get the PDs in the pion for ABFKW-P set 2 | *
1967 * | S AURPI3 : to get the PDs in the pion for ABFKW-P set 3 | *
1968 * | S AURPR : to get the PDs in the nucleon for ABFOW | *
1969 * | F BBETA : DOUBLE PRECISION beta function using DGAMMA | *
1970 * | F BETA : SINGLE PRECISION beta function using GAMMA | *
1971 * | S BXF3 : BEP SF set for up and down valence quarks | *
1972 * | S BXGLUE : BEP SF set for gluons | *
1973 * | S BXQBAR : BEP SF set for antiquarks | *
1974 * | F CTQ1OPD : CTEQ SF sets 1M, 1MS, 1ML, 1D and 1L | *
1975 * | F CTQ1OPF : Description of CTEQ sets 1M, 1MS, 1ML, 1D, 1L | *
1976 * | F CTQ1PD : CTEQ SF sets 1M, 1MS, 1ML, 1D and 1L | *
1977 * | F CTQ1PF : Description of CTEQ sets 1M, 1MS, 1ML, 1D, 1L | *
1978 * | F CTQ2OPD : CTEQ SF sets 2M, 2MS, 2ML, 2D and 2L | *
1979 * | F CTQ2OPF : Description of CTEQ sets 2M, 2MS, 2ML, 2D, 2L | *
1980 * | S CTQ2OPS : CTEQ SF sets 2M, 2MS, 2ML, 2D and 2L | *
1981 * | F CTQ2PPD : CTEQ SF sets 2pM, 2pMS, 2pML, 2pD and 2pL | *
1982 * | F CTQ2PPF : Description of CTEQ sets 2pM, 2pMS, 2pML, 2pD, 2pL*
1983 * | S CTQ2PPS : CTEQ SF sets 2pM, 2pMS, 2pML, 2pD and 2pL | *
1984 * | F CTQ3PD : CTEQ SF sets 3L, 3M and 3D | *
1985 * | S CTQ3PDS : CTEQ SF sets 3L, 3M and 3D | *
1986 * | F CTQ3PF : Description of CTEQ sets 3L, 3M and 3D | *
1987 * | F CTQ4Fn : Description of CTEQ sets 4L, 4D, 4A1-A5, 4HJ and 4lQ
1988 * | F CTQ4PnX : CTEQ SF sets 4L, 4D, 4A1-A5, 4HJ and 4lQ | *
1989 * | S CTQ4RA1 : Parameters of CTEQ set 4 A1 | *
1990 * | S CTQ4RA2 : Parameters of CTEQ set 4 A2 | *
1991 * | S CTQ4RA4 : Parameters of CTEQ set 4 A4 | *
1992 * | S CTQ4RA5 : Parameters of CTEQ set 4 A5 | *
1993 * | S CTQ4RD : Parameters of CTEQ set 4 D | *
1994 * | S CTQ4RHJ : Parameters of CTEQ set 4 HJ | *
1995 * | S CTQ4RL : Parameters of CTEQ set 4 L | *
1996 * | S CTQ4RM : Parameters of CTEQ set 4 M | *
1997 * | S CTQ4RlQ : Parameters of CTEQ set 4 lQ | *
1998 * | F DBFINT : to interpolate linearly PDs in grid (double precision)
1999 * | S DFLM1 : interface to DFLM 160 (NP=1,NG=2,NS=7) | *
2000 * |_______________________________________________________________| *
2001 ******************************* - 28 - ********************************
20021***********************************************************************
2003 * _________________________________________________________________ *
2004 * | Complete list of subprograms, cont'd | *
2005 * | | *
2006 * | S DFLM2 : interface to DFLM 260 (NP=1,NG=2,NS=8) | *
2007 * | S DFLM3 : interface to DFLM 360 (NP=1,NG=2,NS=9) | *
2008 * | S DFLM4 : interface to DFLM soft valence quarks (NP=1,NG=2,NS=1)
2009 * | S DFLM5 : interface to DFLM hard valence quarks (NP=1,NG=2,NS=2)
2010 * | S DFLM6 : interface to DFLM soft gluons (NP=1,NG=2,NS=3) | *
2011 * | S DFLM7 : interface to DFLM hard gluons (NP=1,NG=2,NS=4) | *
2012 * | S DFLM8 : interface to DFLM LO average fit (NP=1,NG=2,NS=5) *
2013 * | S DFLM9 : interface to DFLM NLL average fit (NP=1,NG=2,NS=6)*
2014 * | S DGPHO1 : DG SF set 1 in the photon (NP=3,NG=2,NS=1) | *
2015 * | S DGPHO2 : DG SF set 2 in the photon (NP=3,NG=2,NS=2) | *
2016 * | S DGPHO3 : DG SF set 3 in the photon (NP=3,NG=2,NS=3) | *
2017 * | S DGPHO4 : DG SF set 4 in the photon (NP=3,NG=2,NS=4) | *
2018 * | S DOPHO1 : DO SF set 1 in the photon (NP=3,NG=1,NS=1) | *
2019 * | S DOPHO2 : DO SF set 2 in the photon (NP=3,NG=1,NS=2) | *
2020 * | F DPOLIN : to interpolate PDs in grid (DOUBLE PRECISION) | *
2021 * | S FXAVER : DFLM SF set of LO central average fit | *
2022 * | S FXG160 : DFLM SF set of NLL with QCDL4 = 160 | *
2023 * | S FXG260 : DFLM SF set of NLL with QCDL4 = 260 | *
2024 * | S FXG360 : DFLM SF set of NLL with QCDL4 = 360 | *
2025 * | S FXGHAR : DFLM SF set of LO hard gluon | *
2026 * | S FXGSOF : DFLM SF set of LO soft gluon | *
2027 * | S FXNLLA : DFLM SF set of NLL central average fit | *
2028 * | S FXVHAR : DFLM SF set of LO hard valence | *
2029 * | S FXVSOF : DFLM SF set of LO soft valence | *
2030 * | F GAMFUN : to calculate GAMMA funct. in either SP or DP prec.*
2031 * | F GAMMADO : to calculate GAMMA funct. in SP (Hastings) | *
2032 * | F GHRDV : GHR SF set for down quarks | *
2033 * | F GHRGL : GHR SF set for gluons | *
2034 * | F GHRUV : GHR SF set for up quarks | *
2035 * | F GHRXI : GHR SF set for antiquarks | *
2036 * | F GHRXS : GHR SF set for strange quarks | *
2037 * | F GR94FV : GRV94 SF set for up + down valence quarks + gluons*
2038 * | F GR94FW : GRV94 SF set for antiquarks | *
2039 * | F GR94FWS : GRV94 SF set for strange quarks | *
2040 * | F GRSF1 : GRS SF set part1 in the photon | *
2041 * | F GRSF2 : GRS SF set part2 in the photon | *
2042 * | S GRSGALO : GRS SF set of LO in the photon (NP=3,NG=5,NS=4)| *
2043 * | S GRSVH1 : GRSV SF set of pol. NLL standard (NP=1,NG=5,NS=10)*
2044 * | S GRSVH2 : GRSV SF set of pol. NLL valence (NP=1,NG=5,NS=11) *
2045 * | S GRSVL1 : GRSV SF set of pol. LO standard (NP=1,NG=5,NS=8) *
2046 * | S GRSVL2 : GRSV SF set of pol. LO valence (NP=1,NG=5,NS=9)| *
2047 * | S GRV94DI : GRV94 SF set of DIS in the nucleon (NP=1,NG=5,NS=7)
2048 * | S GRV94HO : GRV94 SF set of MS_bar in the nucleon (NP=1,NG=5,NS=6)
2049 * | S GRV94LO : GRV94 SF set of LO in the nucleon (NP=1,NG=5,NS=5)*
2050 * | F GRVFGP : GRV SF set for gluons in the pion | *
2051 * | F GRVFQBP : GRV SF set for s, c and b quarks in the pion | *
2052 * | F GRVFV : GRV SF set for up + down valence quarks | *
2053 * | F GRVFVP : GRV SF set for valence quarks in the pion | *
2054 * | F GRVFW : GRV SF set for antiquarks and gluons | *
2055 * | F GRVFWS : GRV SF set for strange, charm and bottom quarks| *
2056 * | S GRVGAH0 : GRV SF set lead. terms of NLL in the photon (NP=3,NG=5,NS=1
2057 * | S GRVGAHO : GRV SF set of NLL in the photon (NP=3,NG=5,NS=2) *
2058 * | S GRVGALO : GRV SF set of LO in the photon (NP=3,NG=5,NS=3)| *
2059 * | F GRVGF : GRV SF set for u and d quarks and gluons in the photon
2060 * | F GRVGFS : GRV SF set for s, c and b quarks in the photon | *
2061 * | S GRVHO : GRV SF set of NLL (NP=1,NG=5,NS=3) | *
2062 * | S GRVLO : GRV SF set of LO (NP=1,NG=5,NS=4) | *
2063 * | S GRVPIHO : GRV SF set of NLL in the pion (NP=2,NG=5,NS=1) | *
2064 * | S GRVPILO : GRV SF set of LO in the pion (NP=2,NG=5,NS=2) | *
2065 * | S GSXCOR : to get X coordinates for SFGSHO + LO sets 1, 2 | *
2066 * | S GS2XCOR : to get X coordinates for SFGSHO + LO 96 | *
2067 * | S GS96HO : GS-96 SF set of HO in the photon (NP=3,NG=4,NS=4) *
2068 * | S GS96LO : GS-96 SF set of LO in the photon (NP=3,NG=4,NS=5) *
2069 * |_______________________________________________________________| *
2070 ******************************* - 29 - ********************************
20711***********************************************************************
2072 * _________________________________________________________________ *
2073 * | Complete list of subprograms, cont'd | *
2074 * | | *
2075 * | S HMRS1EB : Description of HMRS sets E, B (1.90-retracted) | *
2076 * | S HMRS2EB : Description of HMRS sets E, B (3.90) | *
2077 * | S HMRS3EB : Description of HMRS sets E, B (4.90) | *
2078 * | S HMRSGEB : Description of HMRS sets E, E+, E- and B | *
2079 * | S KMRSEB : Description of KMRS sets E and B | *
2080 * | S KMRSGEB : Description of KMRS B0 with L=135,160,200,235 | *
2081 * | S MRSDSHD : Description of MRS sets D0', S0', D-', (H) in DIS *
2082 * | S MRSDSHP : Description of MRS sets D0', S0', D-', (H) | *
2083 * | S MRSEB : Description of MRS sets E and B | *
2084 * | S MRSEBP : Description of MRS sets E' and B' | *
2085 * | S NEWDO1 : DO SF set 1.1 (NP=1,NG=1,NS=10) | *
2086 * | S PDFSET : to set all parameters for PDFLIB | *
2087 * | S PDFSTA : to print summary of param. and error statistics| *
2088 * | S PDFVERS : to get date and time of PDFLIB Library release | *
2089 * | F PDXMT : Description of MT sets 1 to 11 | *
2090 * | S PDZXMT : MT SF sets 1 to 11 | *
2091 * | S PFTOPDG : to access PDFLIB with output in PDG format | *
2092 * | S PHLAC1 : interface to LAC-G set 1 in the photon | *
2093 * | S PHLAC2 : interface to LAC-G set 2 in the photon | *
2094 * | S PHLAC3 : interface to LAC-G set 3 in the photon | *
2095 * | S RXDV : OR SF set for down valence quarks | *
2096 * | S RXGLUE : OR SF set for gluons | *
2097 * | S RXQBAR : OR SF set for antiquarks | *
2098 * | S RXUD : OR SF set for the sum of up + down quarks | *
2099 * | S SASAN1 : | *
2100 * | S SASANO : | *
2101 * | S SASBEH : | *
2102 * | S SASDIR : | *
2103 * | S SASGAM1 : | *
2104 * | S SASGAM2 : | *
2105 * | S SASVM1 : | *
2106 * | S SASVMD : | *
2107 * | S SEADIS : BEBC SF set for antiquarks | *
2108 * | S SFABFOW : ABFOW SF set (NP=1,NG=6,NS=1) | *
2109 * | S SFAFG1 : AFG SF set of NLL for Photons (NP=3,NG=6,NS=3) | *
2110 * | S SFBMA : BM SF set A (NP=1,NG=7,NS=1) | *
2111 * | S SFBMB : BM SF set B (NP=1,NG=7,NS=2) | *
2112 * | S SFCTQ11 : CTEQ SF set 1L (LO) (NP=1,NG=4,NS=12) | *
2113 * | S SFCTQ12 : CTEQ SF set 1M (NP=1,NG=4,NS=13) | *
2114 * | S SFCTQ13 : CTEQ SF set 1MS (NP=1,NG=4,NS=14) | *
2115 * | S SFCTQ14 : CTEQ SF set 1ML (NP=1,NG=4,NS=15) | *
2116 * | S SFCTQ15 : CTEQ SF set 1D (NP=1,NG=4,NS=16) | *
2117 * | S SFCTQ21 : CTEQ SF set 2L (LO) (NP=1,NG=4,NS=17) | *
2118 * | S SFCTQ22 : CTEQ SF set 2M (NP=1,NG=4,NS=18) | *
2119 * | S SFCTQ23 : CTEQ SF set 2MS (NP=1,NG=4,NS=19) | *
2120 * | S SFCTQ24 : CTEQ SF set 2MF (NP=1,NG=4,NS=20) | *
2121 * | S SFCTQ25 : CTEQ SF set 2ML (NP=1,NG=4,NS=21) | *
2122 * | S SFCTQ26 : CTEQ SF set 2D (NP=1,NG=4,NS=22) | *
2123 * | S SFCTQ31 : CTEQ SF set 3L (LO) (NP=1,NG=4,NS=29) | *
2124 * | S SFCTQ32 : CTEQ SF set 3M (NP=1,NG=4,NS=30) | *
2125 * | S SFCTQ33 : CTEQ SF set 3D (NP=1,NG=4,NS=31) | *
2126 * | S SFCTQ41 : CTEQ SF set 4L (LO) (NP=1,NG=4,NS=32) | *
2127 * | S SFCTQ42 : CTEQ SF set 4D (NP=1,NG=4,NS=33) | *
2128 * | S SFCTQ43 : CTEQ SF set 4M (NP=1,NG=4,NS=34) | *
2129 * | S SFCTQ44 : CTEQ SF set 4A1 (NP=1,NG=4,NS=35) | *
2130 * | S SFCTQ45 : CTEQ SF set 4A2 (NP=1,NG=4,NS=36) | *
2131 * | S SFCTQ46 : CTEQ SF set 4A4 (NP=1,NG=4,NS=38) | *
2132 * | S SFCTQ47 : CTEQ SF set 4A5 (NP=1,NG=4,NS=39) | *
2133 * | S SFCTQ48 : CTEQ SF set 4HJ (NP=1,NG=4,NS=40) | *
2134 * | S SFCTQ49 : CTEQ SF set 4lQ (NP=1,NG=4,NS=41) | *
2135 * | S SFDISA : MRS SF set A (DIS) (NP=1,NG=3,NS=44) | *
2136 * | S SFDISD0 : MRS SF set D0' (DIS) (NP=1,NG=3,NS=33) | *
2137 * | S SFDISDM : MRS SF set D-' (DIS) (NP=1,NG=3,NS=34) | *
2138 * |_______________________________________________________________| *
2139 ******************************* - 30 - ********************************
21401***********************************************************************
2141 * _________________________________________________________________ *
2142 * | Complete list of subprograms, cont'd | *
2143 * | | *
2144 * | S SFDISH : MRS SF set (H) (NP=1,NG=3,NS=36) | *
2145 * | S SFDISS0 : MRS SF set S0' (DIS) (NP=1,NG=3,NS=32) | *
2146 * | S SFDSQA : MRS SF set A (DIS) low Q**2 (NP=1,NG=3,NS=44) | *
2147 * | S SFEHLQ1 : EHLQ SF set 1 (NP=1,NG=1,NS=8) | *
2148 * | S SFEHLQ2 : EHLQ SF set 2 (NP=1,NG=1,NS=9) | *
2149 * | S SFGRVH : old GRV SF set of NLL (NP=1,NG=5,NS=1) | *
2150 * | S SFGRVL : old GRV SF set of LO (NP=1,NG=5,NS=2) | *
2151 * | S SFGSHO : GS SF set of NLL in the photon (NP=3,NG=4,NS=1)| *
2152 * | S SFGSLO1 : GS SF set 1 of LO in the photon (NP=3,NG=4,NS=2) *
2153 * | S SFGSLO2 : GS SF set 2 of LO in the photon (NP=3,NG=4,NS=3) *
2154 * | S SFKBMR2 : KMRS SF set B- with R2 shadowing (NP=1,NG=3,NS=19)*
2155 * | S SFKBMR5 : KMRS SF set B- with R5 shadowing (NP=1,NG=3,NS=20)*
2156 * | S SFLACG1 : LAC SF set 1 in the photon (NP=3,NG=3,NS=1) | *
2157 * | S SFLACG2 : LAC SF set 2 in the photon (NP=3,NG=3,NS=2) | *
2158 * | S SFLACG3 : LAC SF set 3 in the photon (NP=3,NG=3,NS=3) | *
2159 * | S SFMFAP : MRS SF set Ap parametrisation (NP=1,NG=3,NS=40)| *
2160 * | S SFMFR1 : MRS SF set R1 parametrisation (NP=1,NG=3,NS=57)| *
2161 * | S SFMFR2 : MRS SF set R2 parametrisation (NP=1,NG=3,NS=58)| *
2162 * | S SFMRSD0 : MRS SF set D0 (NP=1,NG=3,NS=27) | *
2163 * | S SFMRSDM : MRS SF set D- (NP=1,NG=3,NS=28) | *
2164 * | S SFMRSG : MRS SF set G (NP=1,NG=3,NS=41) | *
2165 * | S SFMRSJ : MRS SF set J (NP=1,NG=3,NS=51) | *
2166 * | S SFMRSS0 : MRS SF set S0 (NP=1,NG=3,NS=26) | *
2167 * | S SFMS105 : MRS SF set Ap (150) (NP=1,NG=3,NS=45) | *
2168 * | S SFMS110 : MRS SF set Ap (201) (NP=1,NG=3,NS=46) | *
2169 * | S SFMS115 : MRS SF set Ap (266) (NP=1,NG=3,NS=47) | *
2170 * | S SFMS120 : MRS SF set Ap (344) (NP=1,NG=3,NS=48) | *
2171 * | S SFMS125 : MRS SF set Ap (435) (NP=1,NG=3,NS=49) | *
2172 * | S SFMS130 : MRS SF set Ap (542) (NP=1,NG=3,NS=50) | *
2173 * | S SFMSAP : MRS SF set Ap (NP=1,NG=3,NS=39) | *
2174 * | S SFMSBA : MRS SF set A (NP=1,NG=3,NS=37) | *
2175 * | S SFMSBB0 : MRS SF set B0' (NP=1,NG=3,NS= | *
2176 * | S SFMSBD0 : MRS SF set D0' (NP=1,NG=3.NS=30) | *
2177 * | S SFMSBDM : MRS SF set D-' (NP=1,NG=3,NS=31) | *
2178 * | S SFMSBH : MRS SF set (H) (NP=1,NG=3,NS=35) | *
2179 * | S SFMSBS0 : MRS SF set S0' (NP=1,NG=3,NS=29) | *
2180 * | S SFMSFA : MRS SF set A parametrisation (NP=1,NG=3,NS=38) | *
2181 * | S SFMSFG : MRS SF set G parametrisation (NP=1,NG=3,NS=42) | *
2182 * | S SFMSJP : MRS SF set Jp (NP=1,NG=3,NS=52) | *
2183 * | S SFMSQA : MRS SF set A low Q**2 (NP=1,NG=3,NS=43) | *
2184 * | S SFMSR1 : MRS SF set R1 (NP=1,NG=3,NS=53) | *
2185 * | S SFMSR2 : MRS SF set R2 (NP=1,NG=3,NS=54) | *
2186 * | S SFMSR3 : MRS SF set R3 (NP=1,NG=3,NS=55) | *
2187 * | S SFMSR4 : MRS SF set R4 (NP=1,NG=3,NS=56) | *
2188 * | S SFMT10 : interface to MT set 6 (MS_bar) | *
2189 * | S SFMT11 : interface to MT set LO | *
2190 * | S SFMTU1 : interface to MT set S1 (DIS) | *
2191 * | S SFMTU2 : interface to MT set B1 (DIS) | *
2192 * | S SFMTU3 : interface to MT set B2 (DIS) | *
2193 * | S SFMTU4 : interface to MT set E1 (DIS) | *
2194 * | S SFMTU5 : interface to MT set 6 (u_bar + d_bar)/2 (DIS) | *
2195 * | S SFMTU6 : interface to MT set S1 (MS_bar) | *
2196 * | S SFMTU7 : interface to MT set B1 (MS_bar) | *
2197 * | S SFMTU8 : interface to MT set B2 (MS_bar) | *
2198 * | S SFMTU9 : interface to MT set E1 (MS_bar) | *
2199 * | S SFSAS11 : SaS SF set 1D (V1) (NP=3,NG=9,NS=1) | *
2200 * | S SFSAS12 : SaS SF set 1M (V1) (NP=3,NG=9,NS=2) | *
2201 * | S SFSAS13 : SaS SF set 2D (V1) (NP=3,NG=9,NS=3) | *
2202 * | S SFSAS14 : SaS SF set 2M (V1) (NP=3,NG=9,NS=4) | *
2203 * | S SFSAS21 : SaS SF set 1D (V2) (NP=3,NG=9,NS=5) | *
2204 * | S SFSAS22 : SaS SF set 1M (V2) (NP=3,NG=9,NS=6) | *
2205 * | S SFSAS23 : SaS SF set 2D (V2) (NP=3,NG=9,NS=7) | *
2206 * | S SFSAS24 : SaS SF set 2M (V2) (NP=3,NG=9,NS=8) | *
2207 * |_______________________________________________________________| *
2208 ******************************* - 31 - ********************************
22091***********************************************************************
2210 * _________________________________________________________________ *
2211 * | Complete list of subprograms, cont'd | *
2212 * | | *
2213 * | S SFWHI1 : WHIT SF set 1 (NP=3,NG=8,NS=1) | *
2214 * | S SFWHI2 : WHIT SF set 2 (NP=3,NG=8,NS=2) | *
2215 * | S SFWHI3 : WHIT SF set 3 (NP=3,NG=8,NS=3) | *
2216 * | S SFWHI4 : WHIT SF set 4 (NP=3,NG=8,NS=4) | *
2217 * | S SFWHI5 : WHIT SF set 5 (NP=3,NG=8,NS=5) | *
2218 * | S SFWHI6 : WHIT SF set 6 (NP=3,NG=8,NS=6) | *
2219 * | S SMRSP31 : SMRS SF set 1 in the pion (NP=2,NG=3,NS=1) | *
2220 * | S SMRSP32 : SMRS SF set 2 in the pion (NP=2,NG=3,NS=2) | *
2221 * | S SMRSP33 : SMRS SF set 3 in the pion (NP=2,NG=3,NS=3) | *
2222 * | S SMRSPI : Description of SMRS sets 1, 2, 3 in the pion | *
2223 * | S SRCTQ21 : CTEQ SF set 2pL (LO) (NP=1,NG=4,NS=23) | *
2224 * | S SRCTQ22 : CTEQ SF set 2pM (NP=1,NG=4,NS=24) | *
2225 * | S SRCTQ23 : CTEQ SF set 2pMS (NP=1,NG=4,NS=25) | *
2226 * | S SRCTQ24 : CTEQ SF set 2pMF (NP=1,NG=4,NS=26) | *
2227 * | S SRCTQ25 : CTEQ SF set 2pML (NP=1,NG=4,NS=27) | *
2228 * | S SRCTQ26 : CTEQ SF set 2pD (NP=1,NG=4,NS=28) | *
2229 * | S STRBEBC : BEBC SF set (NP=1,NG=1,NS=2) | *
2230 * | S STRBEP : BEP SF set (NP=1,NG=1,NS=4) | *
2231 * | S STRCBP : MRS SF set B' (NP=1,NG=3,NS=7) | *
2232 * | S STRCEP : MRS SF set E' (NP=1,NG=3,NS=6) | *
2233 * | S STRDO1 : DO SF set 1 (NP=1,NG=1,NS=6) | *
2234 * | S STRDO2 : DO SF set 2 (NP=1,NG=1,NS=7) | *
2235 * | S STRGHR : GHR SF set (NP=1,NG=1,NS=5) | *
2236 * | S STRH1B : HMRS SF set B (NP=1,NG=3,NS=9 - retracted) | *
2237 * | S STRH1E : HMRS SF set E (NP=1,NG=3,NS=8 - retracted) | *
2238 * | S STRH2B : HMRS SF set B (NP=1,NG=3,NS=11) | *
2239 * | S STRH2E : HMRS SF set E (NP=1,NG=3,NS=10) | *
2240 * | S STRH3B : HMRS SF set B (NP=1,NG=3,NS=15) | *
2241 * | S STRH3B1 : HMRS SF set B - L=100 (NP=1,NG=3,NS=16) | *
2242 * | S STRH3B3 : HMRS SF set B - L=300 (NP=1,NG=3,NS=17) | *
2243 * | S STRH3E : HMRS SF set E (NP=1,NG=3,NS=14) | *
2244 * | S STRH3EM : HMRS SF set E- (NP=1,NG=3,NS=13) | *
2245 * | S STRH3EP : HMRS SF set E+ (NP=1,NG=3,NS=12) | *
2246 * | S STRKB0 : KMRS SF set B0 (NP=1,NG=3,NS=21), updated version *
2247 * | S STRKB02 : KMRS SF set B0 - L=135 (NP=1,NG=3,NS=22) | *
2248 * | S STRKB03 : KMRS SF set B0 - L=160 (NP=1,NG=3,NS=23) | *
2249 * | S STRKB04 : KMRS SF set B0 - L=200 (NP=1,NG=3,NS=24) | *
2250 * | S STRKB05 : KMRS SF set B0 - L=235 (NP=1,NG=3,NS=25) | *
2251 * | S STRKB0O : KMRS SF set B0 (old version) | *
2252 * | S STRKBM : KMRS SF set B- (NP=1,NG=3,NS=18) | *
2253 * | S STRMRS1 : MRS SF set 1 (NP=1,NG=3,NS=1) | *
2254 * | S STRMRS2 : MRS SF set 2 (NP=1,NG=3,NS=2) | *
2255 * | S STRMRS3 : MRS SF set 3 (NP=1,NG=3,NS=3) | *
2256 * | S STROWP1 : OW SF set 1 in the pion (NP=2,NG=1,NS=1) | *
2257 * | S STROWP2 : OW SF set 2 in the pion (NP=2,NG=1,NS=2) | *
2258 * | S STRPRIV : Dummy subroutine for a new set of SF | *
2259 * | S STRUCB : MRS SF set B (NP=1,NG=3,NS=4) | *
2260 * | S STRUCE : MRS SF set E (NP=1,NG=3,NS=5) | *
2261 * | S STRUCOR : OR SF set (NP=1,NG=1,NS=3) | *
2262 * | S STRUCTF : to access PDFLIB with output in PDFLIB format | *
2263 * | S STRUCTM : to access PDFLIB with output in PDFLIB format | *
2264 * | S STRUCTP : to access PDFLIB with output in PDFLIB format for PHOTONS
2265 * | TESTGAM : main program to test PDFLIB for photons | *
2266 * | TESTGAMO : main program to test PDFLIB for photons (OLD) | *
2267 * | TESTPDF : main program to test PDFLIB in PDFLIB format | *
2268 * | TESTPDFO : main program to test PDFLIB in PDFLIB format (OLD)*
2269 * | TESTPDG : main program to test PDFLIB in PDG format | *
2270 * | TESTPDGO : main program to test PDFLIB in PDG format (OLD)| *
2271 * | TMAINUPV : main program to show how to use the package | *
2272 * | TMAINUPVO: main program to show how to use the package (OLD) *
2273 * | F V3 : BEBC SF set to calculate X * VALDIS | *
2274 * | S VALDIS : BEBC SF set for up + down valence quarks | *
2275 * | S VLAMBD : Entry in PDZXMT to get Lambda value for MT sets| *
2276 * |_______________________________________________________________| *
2277 ******************************* - 32 - ********************************
22781***********************************************************************
2279 * _________________________________________________________________ *
2280 * | Complete list of subprograms, cont'd | *
2281 * | | *
2282 * | S WATE32 : defines 32 point Gaussian quadrature routine | *
2283 * | S WHIT1G : WHIT SF set 1 for gluons | *
2284 * | S WHIT1Q : WHIT SF set 1 for quarks | *
2285 * | S WHIT2G : WHIT SF set 2 for gluons | *
2286 * | S WHIT2Q : WHIT SF set 2 for quarks | *
2287 * | S WHIT3G : WHIT SF set 3 for gluons | *
2288 * | S WHIT3Q : WHIT SF set 3 for quarks | *
2289 * | S WHIT4G : WHIT SF set 4 for gluons | *
2290 * | S WHIT4Q : WHIT SF set 4 for quarks | *
2291 * | S WHIT5G : WHIT SF set 5 for gluons | *
2292 * | S WHIT5Q : WHIT SF set 5 for quarks | *
2293 * | S WHIT6G : WHIT SF set 6 for gluons | *
2294 * | S WHIT6Q : WHIT SF set 6 for quarks | *
2295 * | | *
2296 * | | *
2297 * | | *
2298 * | | *
2299 * | | *
2300 * | | *
2301 * | | *
2302 * | | *
2303 * | | *
2304 * | | *
2305 * | | *
2306 * | | *
2307 * | | *
2308 * | | *
2309 * | | *
2310 * | | *
2311 * | | *
2312 * | | *
2313 * | | *
2314 * | | *
2315 * | | *
2316 * | | *
2317 * | | *
2318 * | | *
2319 * | | *
2320 * | | *
2321 * | | *
2322 * | | *
2323 * | | *
2324 * | | *
2325 * | | *
2326 * | | *
2327 * | | *
2328 * | | *
2329 * | | *
2330 * | | *
2331 * | | *
2332 * | | *
2333 * | | *
2334 * | | *
2335 * | | *
2336 * | | *
2337 * | | *
2338 * | | *
2339 * | | *
2340 * | | *
2341 * | | *
2342 * | | *
2343 * | | *
2344 * | | *
2345 * |_______________________________________________________________| *
2346 ******************************* - 33 - ********************************
23471***********************************************************************
2348 * _________________________________________________________________ *
2349 * | Complete list of subprograms, cont'd | *
2350 * | | *
2351 * | | *
2352 * | C W50510 : IFLPRT - to set print flag | *
2353 * | C W50511 : NPTYPE,NGROUP,NSET,MODE,NFL,LO,TMAS - | *
2354 * | to select structure function set, to define | *
2355 * | number of flavours, flag for leading order | *
2356 * | and value for the mass of the top quark | *
2357 * | C W50512 : QCDL4,QCDL5 - contains Lambda_4(QCD) and | *
2358 * | Lambda_5(QCD) (in GeV) | *
2359 * | C W50513 : XMIN,XMAX.Q2MIN,Q2MAX - contains minimum | *
2360 * | and maximum values of X and Q**2 | *
2361 * | C W50514 : WXMIN,WXMAX,WQ2MIN,WQ2MAX,WTXMIN,WTXMAX,WQ2MIN,| *
2362 * | WQ2MAX - error counting for fractional and | *
2363 * | total number of X and Q**2 | *
2364 * | C W50514W : PDFWGT - to transmit user filled event weight | *
2365 * | C W50515 : IFLSET,IFLSTA - to transmit flag for call to | *
2366 * | subroutine PDFSET or PDFSTA | *
2367 * | C W50516 : FIRST - to force printing name of structure | *
2368 * | function set in case of multiple choices | *
2369 * | C W50517 : N6 - contains logical print unit number | *
2370 * | C W50518 : save of COMMONs W50510, W50511, W50512, W50513 | *
2371 * | C W50519 : NEWVER - to flag NEW or OLD version of PDFLIB | *
2372 * | format | *
2373 * | C W505110 : SFNAME - contains NAME (CHARACTER*8 array with | *
2374 * | dimension SFNAME(NPTYMX,NGRMAX,NSETMX)) for | *
2375 * | each set of PDFs | *
2376 * | C W505120 : NPGSMX(NPTYMX,NGRMAX),NSETFL(NPTYMX,NGRMAX,NSETMX)*
2377 * | - to provide the actual maximum number of stuc-| *
2378 * | ture function sets per particle type and author| *
2379 * | group, NSETFL(I,J,K) = 0 for retracted PDF set | *
2380 * | C W505121 : NPTYCR(MODEMX),NGROCR(MODEMX),NSETCR(MODEMX) | *
2381 * | - to provide the cross reference to the | *
2382 * | NPTYPE, NGROUP and NSET value (NEW format) for | *
2383 * | a given MODE value (OLD format) | *
2384 * | C W505122 : MODECR(NPTYMX,NGRMAX,NSETMX) - to provide the | *
2385 * | cross reference to the MODE value (OLD format) | *
2386 * | for a given NPTYPE, NGROUP and NSET value | *
2387 * | (NEW format) | *
2388 * | P W5051P0 : LFLPRT - to define print level (D=0) | *
2389 * | P W5051P1 : LPTYPE,LGROUP,LNSET,IFL,LORD,TM - to define | *
2390 * | the defaults for a PDF set given by the | *
2391 * | particle type, author group and number of set | *
2392 * | within an author group, defaults for number of | *
2393 * | flavours, leading order flag and mass of the | *
2394 * | top quark (D=1, 5, 3, 5, 2, 100.D0) | *
2395 * | P W5051P2 : MODEMX - maximum number of structure function | *
2396 * | sets in the OLD PDFLIB format | *
2397 * | : NPTYMX,NGRMAX,NSETMX - maximum number of | *
2398 * | structure function sets in the NEW PDFLIB | *
2399 * | format for particle type, author group number | *
2400 * | and number of sets within an author group | *
2401 * | P W5051P2 : MODEMX - to define maximum number of SF sets | *
2402 * | in PDFLIB (OLD format) (D=281) | *
2403 * | P W5051P7 : L6 - to define logical print unit number (D=6) | *
2404 * | P W5051P7 : L6 - to define logical print unit number | *
2405 * | W50511C : Comments only - complete list of available SFs | *
2406 * | | *
2407 * | | *
2408 * | | *
2409 * | | *
2410 * | | *
2411 * | | *
2412 * | | *
2413 * | | *
2414 * |_______________________________________________________________| *
2415 ******************************* - 34 - ********************************
24161***********************************************************************
2417 * *
2418 * *
2419 * *
2420 * *
2421 * *
2422 * *
2423 * *
2424 * *
2425 * *
2426 * *
2427 * *
2428 * List of Tables *
2429 * *
2430 * *
2431 * 1 List of available sets of NUCLEON structure functions in *
2432 * PDFLIB version 7.09 ...................................... 36 *
2433 * *
2434 * *
2435 * 2 List of available sets of PION structure functions in PDFLIB *
2436 * version 7.09 ............................................. 39 *
2437 * *
2438 * *
2439 * 3 List of available sets of PHOTON structure functions in PDFLIB *
2440 * version 7.09 ............................................. 40 *
2441 * *
2442 * *
2443 * *
2444 * *
2445 * *
2446 * *
2447 * *
2448 * *
2449 * *
2450 * *
2451 * *
2452 * *
2453 * *
2454 * *
2455 * *
2456 * *
2457 * *
2458 * *
2459 * *
2460 * *
2461 * *
2462 * *
2463 * *
2464 * *
2465 * *
2466 * *
2467 * *
2468 * *
2469 * *
2470 * *
2471 * *
2472 * *
2473 * *
2474 * *
2475 * *
2476 * *
2477 * *
2478 * *
2479 * *
2480 * *
2481 * *
2482 * *
2483 * *
2484 ******************************* - 35 - ********************************
24851***********************************************************************
2486 * *
2487 * Table 1 *
2488 * *
2489 * Available sets of NUCLEON structure functions in PDFLIB *
2490 * *
2491 * ------------------------------------------------------------------- *
2492 * *
2493 * Nptype *
2494 * Ngroup Q**2 Name of Set Scheme Reference *
2495 * Nset min *
2496 * QCDL4 *
2497 * ------------------------------------------------------------------- *
2498 * 0 200 0.3 ->GRV HO MS_bar NLL Default *
2499 * ------------------------------------------------------------------- *
2500 * ------------------------------------------------------------------- *
2501 * 1 1 1 PRIVATE to user *
2502 * 1 1 2 45 0.5 BEBC LO [Nucl.Phys. B132 (1978) 249]
2503 * 1 1 3 500 1.8 OR LO [Phys.Rev. D17 (1978) 3003]
2504 * 1 1 4 470 5 BEP LO [Z. Phys. C2 (1979) 265]
2505 * 1 1 5 400 4 GHR LO [Z. Phys. C13 (1982) 119]
2506 * 1 1 6 200 4 DO Set 1 LO [Phys.Rev. D30 (1984) 49]
2507 * 1 1 7 400 4 DO Set 2 LO *
2508 * 1 1 8 200 5 EHLQ Set 1 LO [Rev.Mod.Phys. 56 (1984) 579;
2509 * 1 1 9 290 5 EHLQ Set 2 LO Rev.Mod.Phys. 58 (1985) 1065]
2510 * 1 1 10 177 4 DO Set 1.1 LO [Phys.Lett.266B (1991) 3275]
2511 * ------------------------------------------------------------------- *
2512 * 1 2 1 200 10 DFLM soft valence LO [Z. Phys. C39 (1988) 21
2513 * 1 2 2 200 10 DFLM hard valence LO *
2514 * 1 2 3 200 10 DFLM soft gluon LO *
2515 * 1 2 4 200 10 DFLM hard gluon LO *
2516 * 1 2 5 200 10 DFLM central av. LO *
2517 * 1 2 6 300 10 DFLM central av. DIS NLL *
2518 * 1 2 7 160 10 DFLM 160 DIS NLL [unpublished] *
2519 * 1 2 8 260 10 ->DFLM 260 DIS NLL *
2520 * 1 2 9 360 10 DFLM 360 DIS NLL *
2521 * ------------------------------------------------------------------- *
2522 * 1 3 1 107 5 MRS 1 MS_bar NLL [Phys.Rev. D37 (1988) 1181]
2523 * 1 3 2 250 5 MRS 2 MS_bar NLL *
2524 * 1 3 3 178 5 MRS 3 MS_bar NLL *
2525 * 1 3 4 91 5 MRS E MS_bar NLL [Phys.Lett. 206B (1988) 327]
2526 * 1 3 5 228 5 MRS B MS_bar NLL *
2527 * 1 3 6 91 5 MRS E' MS_bar NLL [Mod.Phys. Lett.A4 (1989) 1135
2528 * 1 3 7 228 5 MRS B' MS_bar NLL *
2529 * 1 3 8 100 5 HMRS E (1.90) MS_bar NLL [retracted] *
2530 * 1 3 9 190 5 HMRS B (1.90) MS_bar NLL *
2531 * 1 3 10 100 5 HMRS E (3.90) MS_bar NLL *
2532 * 1 3 11 190 5 HMRS B (3.90) MS_bar NLL *
2533 * 1 3 12 100 5 HMRS E+(4.90) MS_bar NLL [Phys.Lett. 243B (1990) 421]
2534 * 1 3 13 100 5 HMRS E-(4.90) MS_bar NLL *
2535 * 1 3 14 100 5 HMRS E (4.90) MS_bar NLL [Phys.Rev. D4 (1990) 798]
2536 * 1 3 15 190 5 HMRS B (4.90) MS_bar NLL *
2537 * 1 3 16 100 5 HMRS B (8.90) MS_bar NLL [unpublished] *
2538 * 1 3 17 300 5 HMRS B (8.90) MS_bar NLL *
2539 * 1 3 18 190 5 KMRS B-(7.90) MS_bar NLL [Phys.Rev. D42 (1990) 3645]
2540 * 1 3 19 190 5 KMRS B-R2 Shadow (7.90) MS_bar NLL *
2541 * 1 3 20 190 5 KMRS B-R5 Shadow (7.90) MS_bar NLL *
2542 * 1 3 21 190 5 KMRS B0-190 (7.90) MS_bar NLL (as updated 4.92) *
2543 * 1 3 22 135 5 MRS B0-135 (10.90) MS_bar NLL [Phys.Rev. D43 (1991) 3648]
2544 * 1 3 23 160 5 MRS B0-160 (10.90) MS_bar NLL *
2545 * 1 3 24 200 5 MRS B0-200 (10.90) MS_bar NLL *
2546 * 1 3 25 235 5 MRS B0-235 (10.90) MS_bar NLL *
2547 * 1 3 26 215 5 MRS S0 (4.92) MS_bar NLL [retracted] *
2548 * 1 3 27 215 5 MRS D0 (4.92) MS_bar NLL *
2549 * 1 3 28 215 5 MRS D- (4.92) MS_bar NLL *
2550 * 1 3 29 230 5 MRS S0' (11.92) MS_bar NLL [Phys.Lett. 306B (1993)145
2551 * 1 3 30 230 5 MRS D0' (11.92) MS_bar NLL *
2552 * *
2553 ******************************* - 36 - ********************************
25541***********************************************************************
2555 * *
2556 * Table 1 (cont'd) *
2557 * *
2558 * ------------------------------------------------------------------- *
2559 * *
2560 * Nptype *
2561 * Ngroup Q**2 Name of Set Scheme Reference *
2562 * Nset min *
2563 * QCDL4 *
2564 * ------------------------------------------------------------------- *
2565 * 1 3 31 230 5 MRS D-' (11.92) MS_bar NLL *
2566 * 1 3 32 230 5 MRS S0' (2.93) DIS NLL [RAL-92-078 (1992)]
2567 * 1 3 33 230 5 MRS D0' (2.93) DIS NLL *
2568 * 1 3 34 230 5 MRS D-' (2.93) DIS NLL *
2569 * 1 3 35 230 5 MRS (H) (11.93) MS_bar NLL *
2570 * 1 3 36 230 5 MRS (H) (11.93) DIS NLL *
2571 * 1 3 37 230 5 MRS (A) (5.94 MS_bar NLL [RAL-94-055 (1994)]
2572 * 1 3 38 230 5 MRS (A) Fit (5.94) MS_bar NLL *
2573 * 1 3 39 230 5 MRS (Ap) (2.95 MS_bar NLL [RAL-94-104 (1994)]
2574 * 1 3 40 231 5 MRS (Ap) Fit(2.95) MS_bar NLL *
2575 * 1 3 41 255 5 ->MRS (G) (2.95 MS_bar NLL [RAL-95-021 (1995)]
2576 * 1 3 42 255 5 MRS (G) Fit (2.95) MS_bar NLL *
2577 * 1 3 43 230 .625 MRS Set (A) lowQ^2 MS_bar NLL [Phys.Rev. D51 (1995)4756]
2578 * 1 3 44 230 .625 MRS Set (A) lowQ^2 DIS NLL *
2579 * 1 3 45 150 5 MRS Set (Ap-L150) MS_bar NLL [Phys.Lett. 356B (1995)89]
2580 * 1 3 46 201 5 MRS Set (Ap-L201) MS_bar NLL *
2581 * 1 3 47 266 5 MRS Set (Ap-L266) MS_bar NLL *
2582 * 1 3 48 344 5 MRS Set (Ap-L344) MS_bar NLL *
2583 * 1 3 49 435 5 MRS Set (Ap-L435) MS_bar NLL *
2584 * 1 3 50 542 5 MRS Set (Ap-L542) MS_bar NLL *
2585 * 1 3 51 344 5 MRS Set (J) MS_bar NLL [DTP/96/22 (1996)]
2586 * 1 3 52 507 5 MRS Set (J') MS_bar NLL *
2587 * 1 3 53 241 5 MRS Set (R1) MS_bar NLL [DTP/96/44 (1996)]
2588 * 1 3 54 344 5 MRS Set (R2) MS_bar NLL *
2589 * 1 3 55 241 5 MRS Set (R3) MS_bar NLL *
2590 * 1 3 56 344 5 MRS Set (R4) MS_bar NLL *
2591 * 1 3 57 241 5 MRS Set (R1) MS_bar NLL Fit *
2592 * 1 3 58 344 5 MRS Set (R2) MS_bar NLL Fit *
2593 * ------------------------------------------------------------------- *
2594 * 1 4 1 212 4 MT S1 DIS NLL [Z. Phys. C52 (1991) 13]
2595 * 1 4 2 194 4 MT B1 DIS NLL *
2596 * 1 4 3 191 4 MT B2 DIS NLL *
2597 * 1 4 4 155 4 MT E1 DIS NLL *
2598 * 1 4 5 237 4 MT 6 (1/2s) DIS NLL *
2599 * 1 4 6 212 4 MT S1 MS_bar NLL *
2600 * 1 4 7 194 4 MT B1 MS_bar NLL *
2601 * 1 4 8 191 4 MT B2 MS_bar NLL *
2602 * 1 4 9 155 4 MT E1 MS_bar NLL *
2603 * 1 4 10 237 4 MT 6 (1/2s) MS_bar NLL *
2604 * 1 4 11 144 4 MT LO LO *
2605 * 1 4 12 168 4 CTEQ 1L LO [retracted] *
2606 * 1 4 13 231 4 CTEQ 1M MS_bar NLL *
2607 * 1 4 14 231 4 CTEQ 1MS MS_bar NLL *
2608 * 1 4 15 322 4 CTEQ 1ML MS_bar NLL *
2609 * 1 4 16 247 4 CTEQ 1D DIS NLL *
2610 * 1 4 17 190 4 CTEQ 2L LO [retracted] *
2611 * 1 4 18 213 4 CTEQ 2M MS_bar NLL *
2612 * 1 4 19 208 4 CTEQ 2MS MS_bar NLL *
2613 * 1 4 20 208 4 CTEQ 2MF MS_bar NLL *
2614 * 1 4 21 322 4 CTEQ 2ML MS_bar NLL *
2615 * 1 4 22 235 4 CTEQ 2D DIS NLL *
2616 * 1 4 23 190 4 CTEQ 2pL LO [to be published]*
2617 * 1 4 24 213 4 CTEQ 2pM MS_bar NLL *
2618 * 1 4 25 208 4 CTEQ 2pMS MS_bar NLL *
2619 * 1 4 26 208 4 CTEQ 2pMF MS_bar NLL *
2620 * ------------------------------------------------------------------- *
2621 * *
2622 ******************************* - 37 - ********************************
26231***********************************************************************
2624 * *
2625 * Table 1 (cont'd) *
2626 * *
2627 * ------------------------------------------------------------------- *
2628 * *
2629 * Nptype *
2630 * Ngroup Q**2 Name of Set Scheme Reference *
2631 * Nset min *
2632 * QCDL4 *
2633 * ------------------------------------------------------------------- *
2634 * 1 4 27 322 4 CTEQ 2pML MS_bar NLL *
2635 * 1 4 28 235 4 CTEQ 2pD DIS NLL *
2636 * 1 4 29 177 4 CTEQ 3L LO [MSU-HEP/41024] *
2637 * 1 4 30 239 4 ->CTEQ 3M MS_bar NLL *
2638 * 1 4 31 247 4 CTEQ 3D DIS NLL *
2639 * 1 4 32 235 2.56 CTEQ 4LO LO [MSUHEP-60426, CTEQ-604]
2640 * 1 4 33 296 2.56 CTEQ 4D DIS NLL *
2641 * 1 4 34 296 2.56 CTEQ 4M MS_bar NLL *
2642 * 1 4 35 213 2.56 CTEQ 4A1 MS_bar NLL *
2643 * 1 4 36 253 2.56 CTEQ 4A2 MS_bar NLL *
2644 * 1 4 37 296 2.56 CTEQ 4A3 MS_bar NLL *
2645 * 1 4 38 344 2.56 CTEQ 4A4 MS_bar NLL *
2646 * 1 4 39 399 2.56 CTEQ 4A5 MS_bar NLL *
2647 * 1 4 40 302 2.56 CTEQ 4HJ MS_bar NLL *
2648 * 1 4 41 259 0.49 CTEQ 4lQ MS_bar NLL *
2649 * ------------------------------------------------------------------- *
2650 * 1 5 1 160 0.2 old GRV HO MS_bar NLL [retracted] *
2651 * 1 5 2 220 0.2 old GRV LO LO *
2652 * 1 5 3 200 0.3 GRV HO MS_bar NLL [Z. Phys. C53 (1992) 127]
2653 * 1 5 4 200 0.25 GRV LO LO *
2654 * 1 5 5 200 0.4 GRV 94 LO LO [DO-TH 94/24 (1994)]
2655 * 1 5 6 200 0.4 ->GRV 94 HO MS_bar NLL *
2656 * 1 5 7 200 0.4 GRV 94 HO DIS NLL *
2657 * 1 5 8 200 0.4 GRSV pol. LO standard LO [Phys.Lett. 359B (1995)201
2658 * 1 5 9 200 0.4 GRSV pol. LO valence LO *
2659 * 1 5 10 200 0.4 GRSV pol. HO standard LO [Phys.Rev. D53 (1996) 4775
2660 * 1 5 11 200 0.4 GRSV pol. HO valence LO *
2661 * ------------------------------------------------------------------- *
2662 * 1 6 1 230 2 ABFOW MS_bar NLL [Phys.Rev D39 (1989) 3275]
2663 * ------------------------------------------------------------------- *
2664 * 1 7 1 254 5 BM A MS_bar NLL [Phys.Lett. 304B (1993)318
2665 * 1 7 2 254 5 BM B MS_bar NLL *
2666 * ------------------------------------------------------------------- *
2667 * *
2668 * *
2669 * *
2670 * *
2671 * *
2672 * *
2673 * *
2674 * *
2675 * *
2676 * *
2677 * *
2678 * *
2679 * *
2680 * *
2681 * *
2682 * *
2683 * *
2684 * *
2685 * *
2686 * *
2687 * *
2688 * *
2689 * *
2690 * *
2691 ******************************* - 38 - ********************************
26921***********************************************************************
2693 * *
2694 * *
2695 * *
2696 * *
2697 * *
2698 * *
2699 * *
2700 * *
2701 * Table 2 *
2702 * *
2703 * Available sets of PION structure functions in PDFLIB *
2704 * *
2705 * ------------------------------------------------------------------- *
2706 * *
2707 * Nptype *
2708 * Ngroup Q**2 Name of Set Scheme Reference *
2709 * Nset min *
2710 * QCDL4 *
2711 * ------------------------------------------------------------------- *
2712 * 2 1 1 200 4 OW-P Set 1 LO [Phys.Rev. D30 (1984) 943]
2713 * 2 1 2 400 4 OW-P Set 2 LO *
2714 * ------------------------------------------------------------------- *
2715 * 2 3 1 190 5 SMRS-P 1 MS_bar NLL [Phys.Rev. D45 (1992) 2349]
2716 * 2 3 2 190 5 ->SMRS-P 2 MS_bar NLL *
2717 * 2 3 3 190 5 SMRS-P 3 MS_bar NLL *
2718 * ------------------------------------------------------------------- *
2719 * 2 5 1 200 0.3 ->GRV-P HO MS_bar NLL [Z.Phys. C53 (1992) 651]
2720 * 2 5 2 200 0.25 GRV-P LO LO *
2721 * ------------------------------------------------------------------- *
2722 * 2 6 1 231 2 ->ABFKW-P 1 MS_bar NLL [Phys.Lett.233B (189) 517]
2723 * 2 6 2 181 2 ABFKW-P 2 MS_bar NLL *
2724 * 2 6 3 281 2 ABFKW-P 3 MS_bar NLL *
2725 * ------------------------------------------------------------------- *
2726 * *
2727 * *
2728 * *
2729 * *
2730 * *
2731 * *
2732 * *
2733 * *
2734 * *
2735 * *
2736 * *
2737 * *
2738 * *
2739 * *
2740 * *
2741 * *
2742 * *
2743 * *
2744 * *
2745 * *
2746 * *
2747 * *
2748 * *
2749 * *
2750 * *
2751 * *
2752 * *
2753 * *
2754 * *
2755 * *
2756 * *
2757 * *
2758 * *
2759 * *
2760 ******************************* - 39 - ********************************
27611***********************************************************************
2762 * *
2763 * *
2764 * *
2765 * *
2766 * Table 3 *
2767 * *
2768 * Available sets of PHOTON structure functions in PDFLIB *
2769 * *
2770 * ------------------------------------------------------------------- *
2771 * *
2772 * Nptype *
2773 * Ngroup Q**2 Name of Set Scheme Reference *
2774 * Nset min *
2775 * QCDL4 *
2776 * ------------------------------------------------------------------- *
2777 * 3 1 1 380 10 DO-G LO LO [Phys.Rev. D26 (1982) 1600]
2778 * 3 1 2 440 10 ->DO-G HO MS_bar NLL *
2779 * ------------------------------------------------------------------- *
2780 * 3 2 1 400 1 ->DG-G Set 1 LO [Z. Phys. C28 (1985) 451]
2781 * 3 2 2 400 1 DG-G Set 2 LO *
2782 * 3 2 3 400 20 DG-G Set 3 LO *
2783 * 3 2 4 400 200 DG-G Set 4 LO *
2784 * ------------------------------------------------------------------- *
2785 * 3 3 1 200 5 LAC-G 1 LO [Phys.Lett.269B (1991) 458]
2786 * 3 3 2 200 5 ->LAC-G 2 LO *
2787 * 3 3 3 200 5 LAC-G 3 LO *
2788 * ------------------------------------------------------------------- *
2789 * 3 4 1 200 5.3 GS-G HO MS_bar NLL [Z.Phys. C56 (1992) 307]
2790 * 3 4 2 200 5.3 GS-G LO set 1 LO *
2791 * 3 4 3 200 5.3 GS-G LO set 2 LO *
2792 * 3 4 4 200 5.3 ->GS-G-96 HO MS_bar NLL [ANL-HEP-PR-96-33] *
2793 * 3 4 5 200 5.3 GS-G-96 LO LO *
2794 * ------------------------------------------------------------------- *
2795 * 3 5 1 200 0.3 GRV-G LHO DIS*) NLL [Phys.Rev. D46 (1992) 1973;
2796 * 3 5 2 200 0.3 ->GRV-G HO DIS*) NLL Phys.Rev. D45 (1992) 3986]
2797 * 3 5 3 200 0.25 GRV-G LO LO *
2798 * 3 5 4 200 0.6 GRS-G LO LO [Phys. Rev. D51 (1995) 3220
2799 * ------------------------------------------------------------------- *
2800 * 3 6 1 200 2 ACFGP-G HO MS_bar NLL [Z.Phys. C56 (1992) 589]
2801 * 3 6 2 200 2 ->ACFGP-G HO-mc MS_bar NLL *
2802 * 3 6 3 200 2 ->AFG-G HO MS_bar NLL *
2803 * ------------------------------------------------------------------- *
2804 * 3 8 1 400 4 WHIT-G 1 LO [KEK preprint 93-160]*
2805 * 3 8 2 400 4 ->WHIT-G 2 LO *
2806 * 3 8 3 400 4 WHIT-G 3 LO *
2807 * 3 8 4 400 4 WHIT-G 4 LO *
2808 * 3 8 5 400 4 WHIT-G 5 LO *
2809 * 3 8 6 400 4 WHIT-G 6 LO *
2810 * ------------------------------------------------------------------- *
2811 * 3 9 1 200 0.36 SAS-G Set 1D (V1) LO [Z. Phys. C68 (1995) 607]
2812 * 3 9 2 200 0.36 SAS-G Set 1M (V1) LO *
2813 * 3 9 3 200 4 SAS-G Set 2D (V1) LO *
2814 * 3 9 4 200 4 SAS-G Set 2M (V1) LO *
2815 * 3 9 5 200 0.36 SAS-G Set 1D (V2) LO [CERN-TH/96-04] *
2816 * 3 9 6 200 0.36 SAS-G Set 1M (V2) LO *
2817 * 3 9 7 200 4 SAS-G Set 2D (V2) LO *
2818 * 3 9 8 200 4 SAS-G Set 2M (V2) LO *
2819 * ------------------------------------------------------------------- *
2820 * *
2821 * *
2822 * *
2823 * *) not standard, please consult references. *
2824 * *
2825 * *
2826 * *
2827 * *
2828 * *
2829 ******************************* - 40 - ********************************