PDF/spdf/kmrseb.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:28:50  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:24  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13 C-----------------------------------------------------------------------
  14       SUBROUTINE KMRSEB(X,SCALE,MODE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  15 C***************************************************************C
  16 C                                                               C
  17 C     -----  VARIABLE QUARKS AND GLUONS AT SMALL X ----         C
  18 C                                                               C
  19 C     NEW VERSIONS !!!! JULY 1990                               C
  20 C     "........................ " J. KWIECINSKI, A.D. MARTIN,   C
  21 C     R.G. ROBERTS AND W.J. STIRLING PREPRINT DTP-90-46 )       C
  22 C                                                               C
  23 C  MODE 1 CORRESPONDS TO  HARRIMAN,                             C
  24 C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)  WITH LAMBDA= 190 MEV  C
  25 C  WITH SMALL X BEHAVIOUR DETERMINED FROM FIT  "HB FIT"         C
  26 C                                                               C
  27 C  MODE 2 CORRESPONDS TO  KWIECINSKI,                           C
  28 C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)  WITH LAMBDA= 190 MEV  C
  29 C  AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2   "B0 FIT"         C
  30 C                                                               C
  31 C  MODE 3 CORRESPONDS TO  KWIECINSKI,                           C
  32 C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)  WITH LAMBDA= 190 MEV  C
  33 C  AND XG,XQ --> X**-1/2 AS X--> 0 AT Q0**2    "B- FIT"         C
  34 C                                                               C
  35 C  MODE 4 CORRESPONDS TO  KWIECINSKI,                           C
  36 C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)  WITH LAMBDA= 190 MEV  C
  37 C  AND XG,XQ --> X**-1/2 AS X--> 0 AT Q0**2    "B-(5) FIT"      C
  38 C  I.E. WITH WEAK (R=5 GEV^-1) SHADOWING INCLUDED               C
  39 C                                                               C
  40 C  MODE 5 CORRESPONDS TO  KWIECINSKI,                           C
  41 C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)  WITH LAMBDA= 190 MEV  C
  42 C  AND XG,XQ --> X**-1/2 AS X--> 0 AT Q0**2    "B-(2) FIT"      C
  43 C  I.E. WITH STRONG  (R=2 GEV^-1) SHADOWING INCLUDED            C
  44 C                                                               C
  45 C                                                               C
  46 C             >>>>>>>>  CROSS CHECK  <<<<<<<<                   C
  47 C                                                               C
  48 C    THE FIRST NUMBER IN THE "HB" GRID IS  .03058               C
  49 C    THE FIRST NUMBER IN THE "B0" GRID IS  .01727               C
  50 C    THE FIRST NUMBER IN THE "B-" GRID IS  .01543               C
  51 C    THE FIRST NUMBER IN THE "B-5"GRID IS  .0                   C
  52 C    THE FIRST NUMBER IN THE "B-2"GRID IS  .0                   C
  53 C                                                               C
  54 C                                                               C
  55 C                                                               C
  56 C                         -*-                                   C
  57 C                                                               C
  58 C    (NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION        C
  59 C    IS RETURNED I.E. G(X) = GLU/X ETC, AND THAT "SEA"          C
  60 C    IS THE LIGHT QUARK SEA I.E. UBAR(X)=DBAR(X)                C
  61 C    = SEA/X FOR A PROTON.  IF IN DOUBT, CHECK THE              C
  62 C    MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV)          C
  63 C                                                               C
  64 C                         -*-                                   C
  65 C                                                               C
  66 C     (THE RANGE OF APPLICABILITY IS CURRENTLY:                 C
  67 C     10**-5 < X < 1  AND  5 < Q**2 < 1.31 * 10**6              C
  68 C     HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST             C
  69 C     - PROBLEMS, COMMENTS ETC TO WJS@UK.AC.DUR.HEP             C
  70 C                                                               C
  71 C                                                               C
  72 C***************************************************************C
  73 #include "pdf/impdp.inc"
  74       IF(MODE.EQ.1) CALL STRH3B(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  75       IF(MODE.EQ.2) CALL STRKB0O(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  76       IF(MODE.EQ.3) CALL STRKBM(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  77       IF(MODE.EQ.5) CALL SFKBMR2(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  78       IF(MODE.EQ.5) CALL SFKBMR5(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  79       RETURN
  80       END