PDF/spdf/spctq51.F

   1 #include "pdf/pilot.h"
   2 C-----------------------------------------------------------------------
   3       SUBROUTINE SPCTQ51(DX,DSCAL,
   4      +                     DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL)
   5 C
   6 C   ********************************************************************
   7 C   *                                                                  *
   8 C   *        Interface to CTEQset of structure functions               *
   9 C   *                                                                  *
  10 C   *        Author:    H. Plothow-Besch (CERN-PPE)                    *
  11 C   *                                                                  *
  12 C   ********************************************************************
  13 C
  14 C :::::::::::: Structure functions from CETQ5M1 Parametrization
  15 C :::::::::::: Lambda = 0.326 GeV (MSbar) (CTEQ5M1)
  16 C
  17 #include "pdf/impdp.inc"
  18 C
  19       PARAMETER (ISET=1)
  20 #include "pdf/w50511.inc"
  21 C
  22       X = DX
  23       Q = DSCAL
  24 C
  25 C     generate the individual structure fcn calls
  26 C
  27       UPV = Ctq5Mi(1,X,Q)
  28       DNV = Ctq5Mi(2,X,Q)
  29       SEAU = Ctq5Mi(-1,X,Q)
  30       DSEA = SEAU * X
  31       DUPV = (UPV-SEAU) * X
  32       SEAD = Ctq5Mi(-2,X,Q)
  33       DSEAD = SEAD * X
  34       DDNV = (DNV-SEAD) * X
  35       STR = Ctq5Mi(-3,X,Q)
  36       DSTR = STR * X
  37       CHM = Ctq5Mi(-4,X,Q)
  38       DCHM = CHM * X
  39       BOT = Ctq5Mi(-5,X,Q)
  40       DBOT = BOT * X
  41       TOP = 0.
  42       IF (DSCAL.GT.TMAS) TOP = Ctq5Mi(-6,X,QL)
  43       DTOP = TOP * X
  44       GL = Ctq5Mi(0,X,Q)
  45       DGL = GL * X
  46 C
  47       RETURN
  48       END