PDF/spdf/ctq5pdf.F

   1 C------------------------------------------------------------------------------
   2 #include "pdf/pilot.h"
   3 C============================================================================
   4 C                CTEQ Parton Distribution Functions: Version 5.0
   5 C                             Nov. 1, 1999
   6 C
   7 C   Ref: "GLOBAL QCD ANALYSIS OF PARTON STRUCTURE OF THE NUCLEON:
   8 C         CTEQ5 PPARTON DISTRIBUTIONS"
   9 C
  10 C  hep-ph/9903282; to be published in Eur. Phys. J. C 1999.
  11 C
  12 C  These PDF's use quadratic interpolation of attached tables. A parametrized
  13 C  version of the same PDF's without external tables is under construction.
  14 C  They will become available later.
  15 C
  16 C   This package contains 7 sets of CTEQ5 PDF's; plus two updated ones.
  17 C   The undated CTEQ5M1 and CTEQHQ1 use an improved evolution code.
  18 C   Both the original and the updated ones fit current data with comparable
  19 C   accuracy.  The CTEQHQ1 set also involve a different choice of scale,
  20 C   hence differs from CTEQHQ slightly more.  It is preferred over CTEQ5HQ.
  21 C
  22 C   Details are:
  23 C ---------------------------------------------------------------------------
  24 C  Iset   PDF        Description       Alpha_s(Mz)  Lam4  Lam5   Table_File
  25 C ---------------------------------------------------------------------------
  26 C   1    CTEQ5M   Standard MSbar scheme   0.118     326   226    cteq5m.tbl
  27 C   2    CTEQ5D   Standard DIS scheme     0.118     326   226    cteq5d.tbl
  28 C   3    CTEQ5L   Leading Order           0.127     192   146    cteq5l.tbl
  29 C   4    CTEQ5HJ  Large-x gluon enhanced  0.118     326   226    cteq5hj.tbl
  30 C   5    CTEQ5HQ  Heavy Quark             0.118     326   226    cteq5hq.tbl
  31 C   6    CTEQ5F3  Nf=3 FixedFlavorNumber  0.106     (Lam3=395)   cteq5f3.tbl
  32 C   7    CTEQ5F4  Nf=4 FixedFlavorNumber  0.112     309   XXX    cteq5f4.tbl
  33 C         --------------------------------------------------------
  34 C   8    CTEQ5M1  Improved CTEQ5M         0.118     326   226    cteq5m1.tbl
  35 C   9    CTEQ5HQ1 Improved CTEQ5HQ        0.118     326   226    ctq5hq1.tbl
  36 C ---------------------------------------------------------------------------
  37 C
  38 C  The available applied range is 10^-5 << x << 1 and 1.0 << Q << 10,000 (GeV).
  39 C   Lam5 (Lam4, Lam3) represents Lambda value (in MeV) for 5 (4,3) flavors.
  40 C   The matching alpha_s between 4 and 5 flavors takes place at Q=4.5 GeV,
  41 C   which is defined as the bottom quark mass, whenever it can be applied.
  42 C
  43 C   The Table_Files are assumed to be in the working directory.
  44 C
  45 C   Before using the PDF, it is necessary to do the initialization by
  46 C       Call SetCtq5(Iset)
  47 C   where Iset is the desired PDF specified in the above table.
  48 C
  49 C   The function Ctq5Pdf (Iparton, X, Q)
  50 C   returns the parton distribution inside the proton for parton [Iparton]
  51 C   at [X] Bjorken_X and scale [Q] (GeV) in PDF set [Iset].
  52 C   Iparton  is the parton label (5, 4, 3, 2, 1, 0, -1, ......, -5)
  53 C                            for (b, c, s, d, u, g, u_bar, ..., b_bar),
  54 C      whereas CTEQ5F3 has, by definition, only 3 flavors and gluon;
  55 C              CTEQ5F4 has only 4 flavors and gluon.
  56 C
  57 C   For detailed information on the parameters used, e.q. quark masses,
  58 C   QCD Lambda, ... etc.,  see info lines at the beginning of the
  59 C   Table_Files.
  60 C
  61 C   These programs, as provided, are in double precision.  By removing the
  62 C   "Implicit Double Precision" lines, they can also be run in single
  63 C   precision.
  64 C
  65 C   If you have detailed questions concerning these CTEQ5 distributions,
  66 C   or if you find problems/bugs using this package, direct inquires to
  67 C   Hung-Liang Lai(lai@phys.nthu.edu.tw) or Wu-Ki Tung(Tung@pa.msu.edu).
  68 C
  69 C===========================================================================
  70 #if defined(CERNLIB_SINGLE)
  71       Function Ctq5Pdf (Iset, Iparton, X, Q)
  72 #endif
  73 #if defined(CERNLIB_DOUBLE)
  74       DOUBLE PRECISION Function Ctq5Pdf (Iset, Iparton, X, Q)
  75 #endif
  76 #include "pdf/impdp.inc"
  77       Data Isetold, Isetmin, Isetmax / -987, 1, 10 /
  78       save
  79
  80 C             If data file not initialized, do so.
  81       If(Iset.ne.Isetold) then
  82          if(Iset.eq.1)  Call Ctq5M  (Iset)
  83          if(Iset.eq.2)  Call Ctq5D  (Iset)
  84          if(Iset.eq.3)  Call Ctq5L  (Iset)
  85          if(Iset.eq.4)  Call Ctq5HJ (Iset)
  86          if(Iset.eq.5)  Call Ctq5HQ (Iset)
  87          if(Iset.eq.6)  Call Ctq5F3 (Iset)
  88          if(Iset.eq.7)  Call Ctq5F4 (Iset)
  89          if(Iset.eq.8)  Call Ctq5M1 (Iset)
  90          if(Iset.eq.9)  Call Ctq5HQ1 (Iset)
  91          Isetold=Iset
  92       Endif
  93 C
  94       Ctq5Pdf = Ctq5PnX (Iparton, X, Q)
  95       if(Ctq5Pdf.lt.0.D0)  Ctq5Pdf = 0.D0
  96 C
  97       Return
  98 C                             ********************
  99       End