PDF/spdf/sfctq56.F

   1 C-----------------------------------------------------------------------
   2 #include "pdf/pilot.h"
   3       SUBROUTINE SFCTQ56(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 the CTEQ collaboration version 5
  15 C :::::::::::: Lambda3 = 0.395 GeV (Nf=3 Fixed flavour number)
  16 C
  17 #include "pdf/impdp.inc"
  18       PARAMETER (ISET=6)
  19 #include "pdf/w50511.inc"
  20 C
  21       X = DX
  22       Q = DSCAL
  23 C
  24 C     generate the individual structure fcn calls
  25 C
  26       UPV = Ctq5Pdf(ISET, 1,X,Q)
  27       DNV = Ctq5Pdf(ISET, 2,X,Q)
  28       SEAU = Ctq5Pdf(ISET,-1,X,Q)
  29       DSEA = SEAU*X
  30       DUPV = (UPV-SEAU)*X
  31       SEAD = Ctq5Pdf(ISET,-2,X,Q)
  32       DSEAD = SEAD*X
  33       DDNV = (DNV-SEAD)*X
  34       STR = Ctq5Pdf(ISET,-3,X,Q)
  35       DSTR = STR*X
  36       CHM = Ctq5Pdf(ISET,-4,X,Q)
  37       DCHM = CHM*X
  38       BOT = Ctq5Pdf(ISET,-5,X,Q)
  39       DBOT = BOT*X
  40       TOP = 0.
  41 C      IF (DSCAL.GT.TMAS) TOP = Ctq5Pdf(ISET,-6,X,Q)
  42       DTOP = TOP*X
  43       GL = Ctq5Pdf(ISET, 0,X,Q)
  44       DGL = GL*X
  45 C
  46       RETURN
  47       END