spdf/hmrs3eb.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:28:47  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 HMRS3EB(X,SCALE,MODE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  15 C***************************************************************C
  16 C                                                               C
  17 C                                                               C
  18 C     NEW VERSIONS !!!! JANUARY 1990  (AS DESCRIBED IN          C
  19 C     "PARTON DISTRIBUTIONS ... " P.N. HARRIMAN, A.D. MARTIN,   C
  20 C     R.G. ROBERTS AND W.J. STIRLING PREPRINT DTP-90-04 )       C
  21 C                                                               C
  22 C         ********* DEBUGGED APRIL 1990********                 C
  23 C                                                               C
  24 C         ****** NOW DOWN TO X=10**-5 *********                 C
  25 C                                                               C
  26 C  MODE 1 CORRESPONDS TO  HARRIMAN,                             C
  27 C  MARTIN, ROBERTS, STIRLING (EMC FIT)    WITH LAMBDA= 100 MEV  C
  28 C                                                               C
  29 C  MODE 2  CORRESPONDS TO HARRIMAN,                             C
  30 C  MARTIN, ROBERTS, STIRLING (BCDMS FIT)  WITH LAMBDA= 190 MEV  C
  31 C                                                               C
  32 C             >>>>>>>>  CROSS CHECK  <<<<<<<<                   C
  33 C                                                               C
  34 C    THE FIRST NUMBER IN THE "E" GRID IS  .01969                C
  35 C    THE FIRST NUMBER IN THE "B" GRID IS  .03058                C
  36 C                                                               C
  37 C                                                               C
  38 C                         -*-                                   C
  39 C                                                               C
  40 C    (NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION        C
  41 C    IS RETURNED I.E. G(X) = GLU/X ETC, AND THAT "SEA"          C
  42 C    IS THE LIGHT QUARK SEA I.E. UBAR(X)=DBAR(X)=               C
  43 C      SEA/X FOR A PROTON.  IF IN DOUBT, CHECK THE              C
  44 C    MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV)          C
  45 C                                                               C
  46 C                         -*-                                   C
  47 C                                                               C
  48 C     (THE RANGE OF APPLICABILITY IS CURRENTLY:                 C
  49 C     10**-5 < X < 1  AND  5 < Q**2 < 1.31 * 10**6              C
  50 C     HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST             C
  51 C     - PROBLEMS, COMMENTS ETC TO WJS@UK.AC.DUR.HEP             C
  52 C                                                               C
  53 C                                                               C
  54 C***************************************************************C
  55 #include "pdf/impdp.inc"
  56       IF(MODE.EQ.1) CALL STRH3E(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  57       IF(MODE.EQ.2) CALL STRH3B(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
  58       RETURN
  59       END