C============================================================================ C CTEQ Parton Distribution Functions: version 6.0-6.6 C April 10, 2002, v6.01 C February 23, 2003, v6.1 C August 6, 2003, v6.11 C December 12, 2004, v6.12 C December 4, 2006, v6.5 (CTEQ6.5M series added) C March 23, 2007, v6.51 (CTEQ6.5S/C series added) C April 24, 2007, v6.52 (minor improvement) C March 30, 2008, v6.6 C March 23, 2010, v6.6as (AS seies added) C C Ref[1]: "New Generation of Parton Distributions with Uncertainties from Global QCD Analysis" C By: J. Pumplin, D.R. Stump, J.Huston, H.L. Lai, P. Nadolsky, W.K. Tung C JHEP 0207:012(2002), hep-ph/0201195 C C Ref[2]: "Inclusive Jet Production, Parton Distributions, and the Search for New Physics" C By : D. Stump, J. Huston, J. Pumplin, W.K. Tung, H.L. Lai, S. Kuhlmann, J. Owens C JHEP 0310:046(2003), hep-ph/0303013 C C Ref[3]: "Neutrino dimuon Production and Strangeness Asymmetry of the Nucleon" C By: F. Olness, J. Pumplin, S. Stump, J. Huston, P. Nadolsky, H.L. Lai, S. Kretzer, J.F. Owens, W.K. Tung C Eur. Phys. J. C40:145(2005), hep-ph/0312323 C C Ref[4]: "CTEQ6 Parton Distributions with Heavy Quark Mass Effects" C By: S. Kretzer, H.L. Lai, F. Olness, W.K. Tung C Phys. Rev. D69:114005(2004), hep-ph/0307022 C C Ref[5]: "Heavy Quark Mass Effects in Deep Inelastic Scattering and Global QCD Analysis" C By : W.K. Tung, H.L. Lai, A. Belyaev, J. Pumplin, D. Stump, C.-P. Yuan C JHEP 0702:053(2007), hep-ph/0611254 C C Ref[6]: "The Strange Parton Distribution of Nucleon: Global Analysis and Applications" C By : H.L. Lai, P. Nadolsky, J. Pumplin, D. Stump, W.K. Tung, C.-P. Yuan C JHEP 0704:089,2007, hep-ph/0702268 C C Ref[7]: "The Charm Content of the Nucleon" C By : J. Pumplin, H.L. Lai, W.K. Tung C Phys.Rev.D75:054029,2007, hep-ph/0701220 C Ref[8]: "Implications of CTEQ global analysis for collider observables" C By : P. M. Nadolsky, H.-L. Lai, Q.-H. Cao, J. Huston, J. Pumplin, D. R. Stump, W.-K. Tung, C.-P. Yuan C arXiv:0802.0007 [hep-ph], submitted to Phys. Rev. D. C C Ref[9]: TBA C This package contains C (1) 4 standard sets of CTEQ6 PDF's (CTEQ6M, CTEQ6D, CTEQ6L, CTEQ6L1) ; C (2) 40 up/down sets (with respect to CTEQ6M) for uncertainty studies from Ref[1]; C (3) updated version of the above: CTEQ6.1M and its 40 up/down eigenvector sets from Ref[2]. C (4) 5 special sets for strangeness study from Ref[3]. C (5) 1 special set for heavy quark study from Ref[4]. C (6) CTEQ6.5M and its 40 up/down eigenvector sets from Ref[5]. C (7) 8 sets of PDFs resulting from the strangeness study, Ref[6]. C (8) 7 sets of PDFs resulting from the charm study, Ref[7]. C (9) CTEQ6.6M and its 44 up/down eigenvector sets from Ref[8]. C (10) Fits with nonperturbative charm from the study in Ref[8]. C (11) Fits with alternative values of the strong coupling strength from the study in Ref[9]. C Details about the calling convention are: C -------------------------------------------------------------------------------- C Iset PDF-set Description Alpha_s(Mz)**Lam4 Lam5 Table_File Ref C ================================================================================ C Standard, "best-fit", sets: C -------------------------- C 1 CTEQ6M Standard MSbar scheme 0.118 326 226 cteq6m.tbl [1] C 2 CTEQ6D Standard DIS scheme 0.118 326 226 cteq6d.tbl [1] C 3 CTEQ6L Leading Order 0.118** 326** 226 cteq6l.tbl [1] C 4 CTEQ6L1 Leading Order 0.130** 215** 165 cteq6l1.tbl [1] C 200 CTEQ6.1M: updated CTEQ6M (see below, under "uncertainty" section) [2] C 400 CTEQ6.6M; the 2008 set (see below, under "uncertainty" section) [8] C C -------------------------- C Special sets with nonperturbative charm at Q_0=1.3 GeV from Ref [8] C -------------------------- C 450 CTEQ6.6C1 BHPS model for IC 0.118 326 226 ctq66.c1.pds C 451 CTEQ6.6C2 BHPS model for IC 0.118 326 226 ctq66.c2.pds C 452 CTEQ6.6C3 Sea-like model 0.118 326 226 ctq66.c3.pds C 453 CTEQ6.6C4 Sea-like model 0.118 326 226 ctq66.c4.pds C Momentum Fraction carried by c+cbar=2c at Q0=1.3 GeV: C Iset: 451 452 453 454 C Mom. frac: 0.01 0.035 0.01 0.035 C -------------------------- C Special CTEQ6.6AS sets with alternative values of strong coupling strength [9] C -------------------------- C 470 CTEQ6.6AS1 0.116 ctq66.as1.pds C 471 CTEQ6.6AS2 0.117 ctq66.as2.pds C 472 CTEQ6.6AS3 0.119 ctq66.as3.pds C 473 CTEQ6.6AS4 0.120 ctq66.as4.pds C -------------------------- C Special sets for strangeness study: Ref.[3] C -------------------------- C 11 CTEQ6A Class A 0.118 326 226 cteq6sa.pds C 12 CTEQ6B Class B 0.118 326 226 cteq6sb.pds C 13 CTEQ6C Class C 0.118 326 226 cteq6sc.pds C 14 CTEQ6B+ Large [S-] 0.118 326 226 cteq6sb+.pds C 15 CTEQ6B- Negative [S-] 0.118 326 226 cteq6sb-.pds C -------------------------- C Special set for Heavy Quark study: Ref.[4] C -------------------------- C 21 CTEQ6HQ 0.118 326 226 cteq6hq.pds C -------------------------- C Released sets for strangeness study: Ref.[6] C -------------------------- s=sbr C 30 CTEQ6.5S0 Best-fit 0.118 326 226 ctq65.s+0.pds C 31 CTEQ6.5S1 Low s+ 0.118 326 226 ctq65.s+1.pds C 32 CTEQ6.5S2 High s+ 0.118 326 226 ctq65.s+2.pds C 33 CTEQ6.5S3 Alt Low s+ 0.118 326 226 ctq65.s+3.pds C 34 CTEQ6.5S4 Alt High s+ 0.118 326 226 ctq65.s+4.pds C -------------------------- s!=sbr C strangeness asymmetry _s- C 35 CTEQ6.5S-0 Best-fit 0.0014 0.118 326 226 ctq65.s-0.pds C 36 CTEQ6.5S-1 Low -0.0010 0.118 326 226 ctq65.s-1.pds C 37 CTEQ6.5S-2 High 0.0050 0.118 326 226 ctq65.s-2.pds C -------------------------- C Released sets for charm study: Ref.[7] C -------------------------- C 40 CTEQ6.5C0 no intrinsic charm 0.118 326 226 ctq65.c0.pds C 41 CTEQ6.5C1 BHPS model for IC 0.118 326 226 ctq65.c1.pds C 42 CTEQ6.5C2 BHPS model for IC 0.118 326 226 ctq65.c2.pds C 43 CTEQ6.5C3 Meson cloud model 0.118 326 226 ctq65.c3.pds C 44 CTEQ6.5C4 Meson cloud model 0.118 326 226 ctq65.c4.pds C 45 CTEQ6.5C5 Sea-like model 0.118 326 226 ctq65.c5.pds C 46 CTEQ6.5C6 Sea-like model 0.118 326 226 ctq65.c6.pds C C Momentum Fraction carried by c,cbar at Q0=1.3 GeV: C Iset:charm ,cbar | Iset:charm ,cbar | Iset:charm ,cbar C 41: 0.002857,0.002857 | 43: 0.003755,0.004817 | 45: 0.005714,0.005714 C 42: 0.010000,0.010000 | 44: 0.007259,0.009312 | 46: 0.012285,0.012285 C C ============================================================================ C For uncertainty calculations using eigenvectors of the Hessian: C --------------------------------------------------------------- C central + 40 up/down sets along 20 eigenvector directions C ----------------------------- C Original version, Ref[1]: central fit: CTEQ6M (=CTEQ6M.00) C ----------------------- C 1xx CTEQ6M.xx +/- sets 0.118 326 226 cteq6m1xx.tbl C where xx = 01-40: 01/02 corresponds to +/- for the 1st eigenvector, ... etc. C e.g. 100 is CTEQ6M.00 (=CTEQ6M), C 101/102 are CTEQ6M.01/02, +/- sets of 1st eigenvector, ... etc. C ==================================================================== C Updated version, Ref[2]: central fit: CTEQ6.1M (=CTEQ61.00) C ----------------------- C 2xx CTEQ61.xx +/- sets 0.118 326 226 ctq61.xx.tbl C where xx = 01-40: 01/02 corresponds to +/- for the 1st eigenvector, ... etc. C e.g. 200 is CTEQ61.00 (=CTEQ6.1M), C 201/202 are CTEQ61.01/02, +/- sets of 1st eigenvector, ... etc. C ==================================================================== C Version with mass effects, Ref[5]: central fit: CTEQ6.5M (=CTEQ65.00) C ----------------------- C 3xx CTEQ65.xx +/- sets 0.118 326 226 ctq65.xx.pds C where xx = 01-40: 01/02 corresponds to +/- for the 1st eigenvector, ... etc. C e.g. 300 is CTEQ65.00 (=CTEQ6.5M), C 301/302 are CTEQ65.01/02, +/- sets of 1st eigenvector, ... etc. C ==================================================================== C Version with mass effects and free strangeness, Ref[8]: C central fit: CTEQ6.6M (=CTEQ66.00) C ----------------------- C 4xx CTEQ66.xx +/- sets 0.118 326 226 ctq66.xx.pds C where xx = 01-44: 01/02 corresponds to +/- for the 1st eigenvector, ... etc. C e.g. 400 is CTEQ66.00 (=CTEQ6.6M), C 401/402 are CTEQ66.01/02, +/- sets of 1st eigenvector, ... etc. C =========================================================================== C ** ALL fits are obtained by using the same coupling strength C \alpha_s(Mz)=0.118 and the NLO running \alpha_s formula, except CTEQ6L1 C which uses the LO running \alpha_s and its value determined from the fit. C For the LO fits, the evolution of the PDF and the hard cross sections are C calculated at LO. More detailed discussions are given in the references. C C The table grids are generated for C * 10^-8 < x < 1 and 1.3 < Q < 10^5 (GeV) for CTEQ6.6 series; C * 10^-7 < x < 1 and 1.3 < Q < 10^5 (GeV) for CTEQ6.5S/C series; C * 10^-6 < x < 1 and 1.3 < Q < 10,000 (GeV) for CTEQ6, CTEQ6.1 series; C C PDF values outside of the above range are returned using extrapolation. C Lam5 (Lam4) represents Lambda value (in MeV) for 5 (4) flavors. C The matching alpha_s between 4 and 5 flavors takes place at Q=4.5 GeV, C which is defined as the bottom quark mass, whenever it can be applied. C C The Table_Files are assumed to be in the working directory. C C Before using the PDF, it is necessary to do the initialization by C Call SetCtq6(Iset) C where Iset is the desired PDF specified in the above table. C C The function Ctq6Pdf (Iparton, X, Q) C returns the parton distribution inside the proton for parton [Iparton] C at [X] Bjorken_X and scale [Q] (GeV) in PDF set [Iset]. C Iparton is the parton label (5, 4, 3, 2, 1, 0, -1, ......, -5) C for (b, c, s, d, u, g, u_bar, ..., b_bar), C C For detailed information on the parameters used, e.q. quark masses, C QCD Lambda, ... etc., see info lines at the beginning of the C Table_Files. C C These programs, as provided, are in double precision. By removing the C "Implicit Double Precision" lines, they can also be run in single C precision. C C If you have detailed questions concerning these CTEQ6 distributions, C or if you find problems/bugs using this package, direct inquires to C nadolsky@physics.smu.edu, pumplin@pa.msu.edu or hllai@tmue.edu.tw. C C===========================================================================