PDF/spdf/ctq1opf.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.4  1997/07/13 13:22:34  cernlib
   6 * Import version 7.09
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:07  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13       FUNCTION Ctq1OPf (Iset, Iparton, X, Q, Irt)
  14
  15 C     CTEQ distribution function in a parametrized form.
  16 C     (No data tables are needed.)
  17
  18 C    The returned function value is the PROBABILITY density for a given FLAVOR.
  19
  20 C  !!   A companion function Ctq1Pd (Iset, Iparton, X, Q, Irt) gives the
  21 C  !!   VALENCE and SEA MOMENTUM FRACTION distributions.  See next module.
  22
  23 C     Ref.: "CTEQ Parton Distributions and Flavor Dependence of the Sea Quarks"
  24 C         by: J. Botts, J.G. Morfin, J.F. Owens, J. Qiu, W.K. Tung & H. Weerts
  25 C             MSUHEP-92-27, Fermilab-Pub-92/371, FSU-HEP-92-1225, ISU-NP-92-17
  26
  27 C     Since this is an initial distribution, and there may be updates, it is
  28 C     useful for the authors to maintain a record of the distribution list.
  29 C     Please do not freely distribute this program package; instead, refer any
  30 C     interested colleague to direct their request for a copy to:
  31 C     Botts@msupa.pa.msu.edu  or  Botts@msupa (bitnet)  or  MSUHEP::Botts
  32
  33 C     If you have any questions concerning these distributions, direct inquires
  34 C     to Jim Botts or Wu-Ki Tung (username Tung at same E-mail nodes as above).
  35
  36 C $Header$
  37 C $Log$
  38 C Revision 1.1.1.4  1997/07/13 13:22:34  cernlib
  39 C Import version 7.09
  40 C
  41 C Revision 1.1.1.1  1996/04/12 15:29:07  plothow
  42 C Version 7.01
  43 C
  44 C Revision 1.1  93/02/08  18:35:25  wkt
  45 C Initial revision
  46 C
  47 C
  48 C     This function returns the CTEQ parton distributions f^Iset_Lp/proton
  49 C     where Iset (= 1, 2, ..., 5) is the set label;
  50 C       Name convention for CTEQ distributions:  CTEQnSx  where
  51 C            n : version number                      (currently n = 1)
  52 C            S : factorization scheme label: = [M D L] for [MS-bar DIS LO] resp.
  53 C            x : special characteristics, if any
  54 C                        (e.g. S for singular gluon, L for "LEP lambda value")
  55
  56 C     Lp  is the parton label (6, 5, 4, 3, 2, 1, 0, -1, ......, -6)
  57 C                         for (t, b, c, s, d, u, g, u_bar, ..., t_bar)
  58
  59 C     X, Q are the usual x, Q; Irt is a return error code (not implemented yet).
  60
  61 C --> Iset = 1, 2, 3, 4, 5 correspond to the following CTEQ global fits:
  62 C     cteq1M, cteq1MS, cteq1ML, cteq1D, cteq1L  respectively.
  63
  64 C --> QCD parameters for parton distribution set Iset can be obtained inside
  65 C         the user's program by:
  66 C     Call ParCtqO
  67 C    >        (Iset, Iord, Ischeme, MxFlv,
  68 C    >         Alam4, Alam5, Alam6, Amas4, Amas5, Amas6,
  69 C    >         Xmin, Qini, Qmax, ExpNor)
  70 C     where the arguments are self-explanary -- see details in next module.
  71
  72 C     The range of (x, Q) used in this round of global analysis is, approxi-
  73 C     mately,  0.01 < x < 0.75 ; and 4 GeV^2 < Q^2 < 400 GeV^2.
  74
  75 C      The range of (x, Q) used in the reparametrization of the QCD evolved
  76 C      parton distributions is 10E-5 < x < 1 ; 2 GeV < Q < 1 TeV.  The
  77 C      functional form of this parametrization is:
  78
  79 C      A0 * x^A1 * (1-x)^A2 * (1 + A3 * x^A4) * [log(1+1/x)]^A5
  80
  81 C      with the A'coefficients being smooth functions of Q.
  82 C      Since this function is positive definite and smooth, it provides sensible
  83 C      extrapolations of the parton distributions if they are called beyond
  84 C      the original range in an application.   There is no artificial boundaries
  85 C      or sharp cutoff's.
  86
  87 C     IMPLICIT DOUBLE PRECISION (A-H, O-Z)
  88 C+SEQ, IMPDP.
  89       REAL Ctq1OPd,X,Q
  90 C
  91       Ifl = Iparton
  92       JFL = ABS(Ifl)
  93 C                                                                       Valence
  94       IF (Ifl .LE. 0) THEN
  95         VL = 0.
  96       ELSEIF (Ifl .LE. 2) THEN
  97         VL = Ctq1OPd (Iset, Ifl, X, Q, Irt)
  98       ELSE
  99         VL = 0.
 100       ENDIF
 101 C                                                                         Sea
 102       SEA = Ctq1OPd (Iset, -JFL, X, Q, Irt)
 103 C                                               Full (probability) Distribution
 104       Ctq1OPf = (VL + SEA) / X
 105
 106       Return
 107 C                         *************************
 108       END