PDF/spdf/seadis.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:29:19  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:29  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13 C----------------------------------------------------------------------
  14       FUNCTION SEADIS(XX)
  15       DIMENSION
  16      +       SQ2(2:3)
  17       COMMON/W5051I1/
  18      +           GPSI(2:3),GAA(2:3),GPSIA(2:3),GAPSI(2:3),GMINUS(2:3),
  19      1           AL(2:3),BET(2:3),GPLUS(2:3),
  20      2           VNQ02(2:3),SNQ02(2:3),GNQ02(2:3)
  21       COMMON/W5051I2/ Q2L
  22 *- sea quarks' momenta
  23       PARAMETER (Q02=0.5,ALAM2=0.045*0.045)
  24
  25       DATA ZERO/0./,HALF/0.5/,ONE/1./,TWO/2./
  26       DATA PTF/0.25/,PSF/0.75/
  27       DATA IRAZ/0/, Q2OLD/-1./
  28       SAVE SQ2, BSQ2
  29
  30       IF(Q2L.NE.Q2OLD) THEN
  31           Q2OLD = Q2L
  32           S     = LOG( LOG( MAX(Q02,Q2L)/ALAM2)/ LOG( Q02/ALAM2))
  33           DO I=2,3
  34              D1    = SNQ02(I)*EXP(-GPSI(I)*S)
  35              D2PLUS=(ONE-AL(I))*(SNQ02(I)+VNQ02(I))-BET(I)*GNQ02(I)
  36              D2MINU=AL(I)*(SNQ02(I)+VNQ02(I))+BET(I)*GNQ02(I)
  37              D2    =D2PLUS*EXP(-GPLUS (I)*S)    +
  38      1              D2MINU*EXP(-GMINUS(I)*S)    -
  39      2              VNQ02(I)*EXP(-GPSI(I)*S)
  40              SQ2(I)=PSF*D2   + PTF*D1
  41           ENDDO
  42           BSQ2=SQ2(2)/SQ2(3) - TWO
  43       ENDIF
  44       SEADIS=SQ2(2) * (ONE-XX)**BSQ2 * (BSQ2+ONE)
  45       RETURN
  46       END