]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | 1*********************************************************************** |
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 | *********************************************************************** | |
70 | 1*********************************************************************** | |
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 - ******************************** | |
139 | 1*********************************************************************** | |
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 - ******************************** | |
208 | 1*********************************************************************** | |
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 - ******************************** | |
277 | 1*********************************************************************** | |
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 - ******************************** | |
346 | 1*********************************************************************** | |
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 - ******************************** | |
415 | 1*********************************************************************** | |
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 - ******************************** | |
484 | 1*********************************************************************** | |
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 - ******************************** | |
553 | 1*********************************************************************** | |
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 - ******************************** | |
622 | 1*********************************************************************** | |
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 - ******************************** | |
691 | 1*********************************************************************** | |
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 - ******************************** | |
760 | 1*********************************************************************** | |
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 - ******************************** | |
829 | 1*********************************************************************** | |
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 - ******************************** | |
898 | 1*********************************************************************** | |
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 - ******************************** | |
967 | 1*********************************************************************** | |
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 - ******************************** | |
1036 | 1*********************************************************************** | |
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 - ******************************** | |
1105 | 1*********************************************************************** | |
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 - ******************************** | |
1174 | 1*********************************************************************** | |
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 - ******************************** | |
1243 | 1*********************************************************************** | |
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 - ******************************** | |
1312 | 1*********************************************************************** | |
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 - ******************************** | |
1381 | 1*********************************************************************** | |
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 - ******************************** | |
1450 | 1*********************************************************************** | |
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 - ******************************** | |
1519 | 1*********************************************************************** | |
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 - ******************************** | |
1588 | 1*********************************************************************** | |
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 - ******************************** | |
1657 | 1*********************************************************************** | |
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 - ******************************** | |
1726 | 1*********************************************************************** | |
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 - ******************************** | |
1795 | 1*********************************************************************** | |
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 - ******************************** | |
1864 | 1*********************************************************************** | |
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 - ******************************** | |
1933 | 1*********************************************************************** | |
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 - ******************************** | |
2002 | 1*********************************************************************** | |
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 - ******************************** | |
2071 | 1*********************************************************************** | |
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 - ******************************** | |
2140 | 1*********************************************************************** | |
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 - ******************************** | |
2209 | 1*********************************************************************** | |
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 - ******************************** | |
2278 | 1*********************************************************************** | |
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 - ******************************** | |
2347 | 1*********************************************************************** | |
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 - ******************************** | |
2416 | 1*********************************************************************** | |
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 - ******************************** | |
2485 | 1*********************************************************************** | |
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 - ******************************** | |
2554 | 1*********************************************************************** | |
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 - ******************************** | |
2623 | 1*********************************************************************** | |
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 - ******************************** | |
2692 | 1*********************************************************************** | |
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 - ******************************** | |
2761 | 1*********************************************************************** | |
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 - ******************************** |