PDF/spdf/sfsas24.F

   1 C-----------------------------------------------------------------------
   2 #include "pdf/pilot.h"
   3       SUBROUTINE SFSAS24(DX,DQ2,DP2,IP2,
   4      +                     DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL)
   5 C
   6 C   ********************************************************************
   7 C   *                                                                  *
   8 C   *        Interface to SASset 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 SAS group version 2
  15 C :::::::::::: Lambda = 0.200 GeV, Q**2 = 4 GeV**2 (MS_bar)
  16 C
  17 #include "pdf/expdp.inc"
  18      +          DX,DQ2,DP2,
  19      +          DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL
  20       DIMENSION XPDFGM(-6:6)
  21       REAL X, Q, Q2, P2, F2GAM, XPDFGM
  22       PARAMETER (ISET=4)
  23 #include "pdf/w50511.inc"
  24 C
  25       X  = DX
  26       Q  = SQRT(DQ2)
  27       Q2 = DQ2
  28       P2 = DP2
  29       IP = IP2
  30       IF(IP.LT.0 .OR. IP.GT.7) IP=0
  31 C
  32 C     generate the individual structure fcn calls
  33 C
  34       CALL SASGAM2(ISET,X,Q2,P2,IP,F2GAM,XPDFGM)
  35       UPV = XPDFGM(2)
  36       DUPV = UPV
  37       DNV = XPDFGM(1)
  38       DDNV = DNV
  39       SEAU = XPDFGM(-2)
  40       DSEA = SEAU
  41       SEAD = XPDFGM(-1)
  42       DSEAD = SEAD
  43       STR = XPDFGM(-3)
  44       DSTR = STR
  45       CHM = XPDFGM(-4)
  46       DCHM = CHM
  47       BOT = XPDFGM(-5)
  48       DBOT = BOT
  49       TOP = 0.
  50 C      IF (DSCAL.GT.TMAS) TOP = XPDFGM(6)
  51       DTOP = TOP
  52       GL = XPDFGM(0)
  53       DGL = GL
  54 C
  55       RETURN
  56       END