]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | C----------------------------------------------------------------------- |
2 | #include "pdf/pilot.h" | |
3 | C============================================================================ | |
4 | C CTEQ Parton Distribution Functions: Version 4.6 | |
5 | C June 21, 1996 | |
6 | C Modified: 10/17/96, 1/7/97, 1/15/97 | |
7 | C 2/17/97, 2/21/97 | |
8 | C Last Modified on April 2, 1997 | |
9 | C | |
10 | C Ref[1]: "IMPROVED PARTON DISTRIBUTIONS FROM GLOBAL ANALYSIS OF RECENT DEEP | |
11 | C INELASTIC SCATTERING AND INCLUSIVE JET DATA" | |
12 | C By: H.L. Lai, J. Huston, S. Kuhlmann, F. Olness, J. Owens, D. Soper | |
13 | C W.K. Tung, H. Weerts | |
14 | C Phys. Rev. D55, 1280 (1997) | |
15 | C | |
16 | C Ref[2]: "CHARM PRODUCTION AND PARTON DISTRIBUTIONS" | |
17 | C By: H.L. Lai and W.K. Tung | |
18 | C MSU-HEP-61222, CTEQ-622, e-Print Archive: hep-ph/9701256 | |
19 | C to appear in Z. Phys. | |
20 | C | |
21 | C This package contains 13 sets of CTEQ4 PDF's. Details are: | |
22 | C --------------------------------------------------------------------------- | |
23 | C Iset PDF Description Alpha_s(Mz) Lam4 Lam5 Table_File | |
24 | C --------------------------------------------------------------------------- | |
25 | C Ref[1] | |
26 | C 1 CTEQ4M Standard MSbar scheme 0.116 298 202 cteq4m.tbl | |
27 | C 2 CTEQ4D Standard DIS scheme 0.116 298 202 cteq4d.tbl | |
28 | C 3 CTEQ4L Leading Order 0.132 236 181 cteq4l.tbl | |
29 | C 4 CTEQ4A1 Alpha_s series 0.110 215 140 cteq4a1.tbl | |
30 | C 5 CTEQ4A2 Alpha_s series 0.113 254 169 cteq4a2.tbl | |
31 | C 6 CTEQ4A3 ( same as CTEQ4M ) | |
32 | C 7 CTEQ4A4 Alpha_s series 0.119 346 239 cteq4a4.tbl | |
33 | C 8 CTEQ4A5 Alpha_s series 0.122 401 282 cteq4a5.tbl | |
34 | C 9 CTEQ4HJ High Jet 0.116 303 206 cteq4hj.tbl | |
35 | C 10 CTEQ4LQ Low Q0 0.114 261 174 cteq4lq.tbl | |
36 | C --------------------------------------------------------------------------- | |
37 | C Ref[2] | |
38 | C 11 CTEQ4HQ Heavy Quark 0.116 298 202 cteq4hq.tbl | |
39 | C 12 CTEQ4HQ1 Heavy Quark:Q0=1,Mc=1.3 0.116 298 202 cteq4hq1.tbl | |
40 | C (Improved version of CTEQ4HQ, recommended) | |
41 | C 13 CTEQ4F3 Nf=3 FixedFlavorNumber 0.106 (Lam3=385) cteq4f3.tbl | |
42 | C 14 CTEQ4F4 Nf=4 FixedFlavorNumber 0.111 292 XXX cteq4f4.tbl | |
43 | C --------------------------------------------------------------------------- | |
44 | C | |
45 | C The available applied range is 10^-5 < x < 1 and 1.6 < Q < 10,000 (GeV) | |
46 | C except CTEQ4LQ(4HQ1) for which Q starts at a lower value of 0.7(1.0) GeV. | |
47 | C Lam5 (Lam4, Lam3) represents Lambda value (in MeV) for 5 (4,3) flavors. | |
48 | C The matching alpha_s between 4 and 5 flavors takes place at Q=5.0 GeV, | |
49 | C which is defined as the bottom quark mass, whenever it can be applied. | |
50 | C | |
51 | C The Table_Files are assumed to be in the working directory. | |
52 | C | |
53 | C Before using the PDF, it is necessary to do the initialization by | |
54 | C Call SetCtq4(Iset) | |
55 | C where Iset is the desired PDF specified in the above table. | |
56 | C | |
57 | C The function Ctq4Pdf (Iparton, X, Q) | |
58 | C returns the parton distribution inside the proton for parton [Iparton] | |
59 | C at [X] Bjorken_X and scale [Q] (GeV) in PDF set [Iset]. | |
60 | C Iparton is the parton label (5, 4, 3, 2, 1, 0, -1, ......, -5) | |
61 | C for (b, c, s, d, u, g, u_bar, ..., b_bar), | |
62 | C whereas CTEQ4F3 has, by definition, only 3 flavors and gluon; | |
63 | C CTEQ4F4 has only 4 flavors and gluon. | |
64 | C | |
65 | C For detailed information on the parameters used, e.q. quark masses, | |
66 | C QCD Lambda, ... etc., see info lines at the beginning of the | |
67 | C Table_Files. | |
68 | C | |
69 | C These programs, as provided, are in double precision. By removing the | |
70 | C "Implicit Double Precision" lines, they can also be run in single | |
71 | C precision. | |
72 | C | |
73 | C If you have detailed questions concerning these CTEQ4 distributions, | |
74 | C or if you find problems/bugs using this package, direct inquires to | |
75 | C Hung-Liang Lai(Lai_H@pa.msu.edu) or Wu-Ki Tung(Tung@pa.msu.edu). | |
76 | C | |
77 | C=========================================================================== | |
78 | #if defined(CERNLIB_SINGLE) | |
79 | Function Ctq4Fn (Iset, Iparton, X, Q) | |
80 | #endif | |
81 | #if defined(CERNLIB_DOUBLE) | |
82 | DOUBLE PRECISION FUNCTION Ctq4Fn (Iset, Iparton, X, Q) | |
83 | #endif | |
84 | #include "pdf/impdp.inc" | |
85 | Data Isetold, Isetmin, Isetmax / -987, 1, 10 / | |
86 | save | |
87 | ||
88 | C If data file not initialized, do so. | |
89 | If(Iset.ne.Isetold) then | |
90 | if(Iset.eq.1) Call Ctq4RM (Iset) | |
91 | if(Iset.eq.2) Call Ctq4RD (Iset) | |
92 | if(Iset.eq.3) Call Ctq4RL (Iset) | |
93 | if(Iset.eq.4) Call Ctq4RA1 (Iset) | |
94 | if(Iset.eq.5) Call Ctq4RA2 (Iset) | |
95 | if(Iset.eq.6) Call Ctq4RM (Iset) | |
96 | if(Iset.eq.7) Call Ctq4RA4 (Iset) | |
97 | if(Iset.eq.8) Call Ctq4RA5 (Iset) | |
98 | if(Iset.eq.9) Call Ctq4RHJ (Iset) | |
99 | if(Iset.eq.10) Call Ctq4RlQ (Iset) | |
100 | if(Iset.eq.11) Call Ctq4RHQ (Iset) | |
101 | if(Iset.eq.12) Call Ctq4R_HQ1 (Iset) | |
102 | if(Iset.eq.13) Call Ctq4RF3 (Iset) | |
103 | if(Iset.eq.14) Call Ctq4RF4 (Iset) | |
104 | Isetold=Iset | |
105 | Endif | |
106 | C | |
107 | Ctq4Fn = Ctq4PnX (Iparton, X, Q) | |
108 | if(Ctq4Fn.lt.0.D0) Ctq4Fn = 0.D0 | |
109 | C | |
110 | Return | |
111 | C ******************** | |
112 | End |