C============================================================================ C CTEQ Parton Distribution Functions: version 6 & 6.5 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 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 hep-ph/0702268, to appear in JHEP C C Ref[7]: "The Charm Content of the Nucleon" C By : J. Pumplin, H.L. Lai, W.K. Tung C hep-ph/0701220, to appear in Phys. Rev. D C 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 C Details about calling convention are: C --------------------------------------------------------------------------- C Iset PDF-set Description Alpha_s(Mz)**Lam4 Lam5 Table_File C =========================================================================== C Standard, "best-fit", sets: Ref.[1] C -------------------------- C 1 CTEQ6M Standard MSbar scheme 0.118 326 226 cteq6m.tbl C 2 CTEQ6D Standard DIS scheme 0.118 326 226 cteq6d.tbl C 3 CTEQ6L Leading Order 0.118** 326** 226 cteq6l.tbl C 4 CTEQ6L1 Leading Order 0.130** 215** 165 cteq6l1.tbl C 200 CTEQ6.1M Ref.[2]: updated CTEQ6M (see below, under "uncertainty" section) 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 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 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 ** 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 10^-6 < x < 1 and 1.3 < Q < 10,000 (GeV). C (For CTEQ6.5S/C series: 10^-7 < x < 1 and 1.3 < Q < 10^5 (GeV)) 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 Pumplin@pa.msu.edu or Tung@pa.msu.edu. C C===========================================================================