PDF/spdf/rxqbar.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:29:17  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       SUBROUTINE RXQBAR (X,QQ,XQBAR)
  14 C.......................................................................
  15 C.
  16 C.    SEA
  17 C.
  18 C.......................................................................
  19 #include "pdf/impdp.inc"
  20       DATA XLAM /0.50D00/, Q0/2.1D0/, Q00/1.8D0/
  21       DATA ONED/1.D0/, THRED/3.D0/
  22       IF(QQ.LT.Q0) QQ=Q0
  23       S=LOG(QQ/XLAM**2)/LOG(Q00/XLAM**2)
  24       S=LOG(S)
  25       AS=0.1467D00-0.1211D00*S+0.0274D00*S**2
  26       ETAS=7.0D00+0.0217D00*S+0.0037D00*S**2
  27       APS=0.1853D00*S-0.0608D00*S**2
  28       ETAPS=9.5041D00+1.0165D00*S-0.1049D00*S**2
  29       BS=0.1682D00*S+0.4473D00*S**2
  30       CS=25.8997D00+3.9572D00*S+1.6331D00*S**2
  31       XQBAR=AS*((ONED-X)**ETAS)+APS*((ONED-X)**ETAPS)+BS*EXP(-CS*X)
  32       RETURN
  33       END