PDF/spdf/kmrsgeb.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:28:51  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 KMRSGEB(X,SCALE,MODE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  15 C***************************************************************C
  16 C                                                               C
  17 C     -----  VARIABLE LAMBDA AND GLUONS  ----                   C
  18 C                                                               C
  19 C     NEW VERSIONS !!!! OCTOBER 1990                            C
  20 C     "........................ "  A.D. MARTIN,                 C
  21 C     R.G. ROBERTS AND W.J. STIRLING PREPRINT DTP-90-?? )       C
  22 C                                                               C
  23 C  MODE 1 CORRESPONDS TO  KWIECINSKI,                           C
  24 C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)                        C
  25 C  WITH LAMBDA(4) = 190 MEV, ETAG = 5.10                        C
  26 C  AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2   "B0 FIT"         C
  27 C                                                               C
  28 C  MODE 2 CORRESPONDS TO                                        C
  29 C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)                        C
  30 C  WITH LAMBDA(4) = 135 MEV, ETAG = 4.65                        C
  31 C  AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2   "B135 FIT"       C
  32 C                                                               C
  33 C  MODE 3 CORRESPONDS TO                                        C
  34 C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)                        C
  35 C  WITH LAMBDA(4) = 160 MEV, ETAG = 4.25                        C
  36 C  AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2   "B160 FIT"       C
  37 C                                                               C
  38 C  MODE 4 CORRESPONDS TO                                        C
  39 C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)                        C
  40 C  WITH LAMBDA(4) = 200 MEV, ETAG = 5.65                        C
  41 C  AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2   "B200 FIT"       C
  42 C                                                               C
  43 C  MODE 5 CORRESPONDS TO                                        C
  44 C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)                        C
  45 C  WITH LAMBDA(4) = 235 MEV, ETAG = 5.20                        C
  46 C  AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2   "B235 FIT"       C
  47 C                                                               C
  48 C                                                               C
  49 C             >>>>>>>>  CROSS CHECK  <<<<<<<<                   C
  50 C                                                               C
  51 C    THE FIRST NUMBER IN THE "B0"   GRID IS  .01727             C
  52 C    THE FIRST NUMBER IN THE "B135" GRID IS  .01683             C
  53 C    THE FIRST NUMBER IN THE "B160" GRID IS  .01663             C
  54 C    THE FIRST NUMBER IN THE "B200" GRID IS  .01601             C
  55 C    THE FIRST NUMBER IN THE "B235" GRID IS  .01571             C
  56 C                                                               C
  57 C                                                               C
  58 C                                                               C
  59 C                         -*-                                   C
  60 C                                                               C
  61 C    (NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION        C
  62 C    IS RETURNED I.E. G(X) = GLU/X ETC, AND THAT "SEA"          C
  63 C    IS THE LIGHT QUARK SEA I.E. UBAR(X)=DBAR(X)                C
  64 C    = SEA/X FOR A PROTON.  IF IN DOUBT, CHECK THE              C
  65 C    MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV)          C
  66 C                                                               C
  67 C                         -*-                                   C
  68 C                                                               C
  69 C     (THE RANGE OF APPLICABILITY IS CURRENTLY:                 C
  70 C     10**-5 < X < 1  AND  5 < Q**2 < 1.31 * 10**6              C
  71 C     HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST             C
  72 C     - PROBLEMS, COMMENTS ETC TO WJS@UK.AC.DUR.HEP             C
  73 C                                                               C
  74 C                                                               C
  75 C***************************************************************C
  76 #include "pdf/impdp.inc"
  77       IF(MODE.EQ.1) CALL STRKB0O(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  78       IF(MODE.EQ.2) CALL STRKB02(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  79       IF(MODE.EQ.3) CALL STRKB03(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  80       IF(MODE.EQ.4) CALL STRKB04(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  81       IF(MODE.EQ.5) CALL STRKB05(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  82       RETURN
  83       END