5 * Revision 1.1.1.2 1996/10/30 08:28:49 cernlib
8 * Revision 1.1.1.1 1996/04/12 15:29:24 plothow
12 #include "pdf/pilot.h"
13 C-----------------------------------------------------------------------
14 SUBROUTINE HMRSGEB(X,SCALE,MODE,UPV,DNV,SEA,STR,CHM,BOT,GL)
15 C***************************************************************C
17 C ----- VARIABLE GLUONS ---- C
19 C NEW VERSIONS !!!! JANUARY 1990 (AS DESCRIBED IN C
20 C "PARTON DISTRIBUTIONS ... " P.N. HARRIMAN, A.D. MARTIN, C
21 C R.G. ROBERTS AND W.J. STIRLING PREPRINT DTP-90-04 ) C
23 C MODE 1 CORRESPONDS TO HARRIMAN, C
24 C MARTIN, ROBERTS, STIRLING (EMC FIT) WITH LAMBDA= 100 MEV C
25 C AND XG --> CONSTANT AS X--> 0 AT Q0**2 C
27 C MODE 2 CORRESPONDS TO HARRIMAN, C
28 C MARTIN, ROBERTS, STIRLING (EMC FIT) WITH LAMBDA= 100 MEV C
29 C AND XG --> X**0.5 AS X--> 0 AT Q0**2 C
31 C MODE 3 CORRESPONDS TO HARRIMAN, C
32 C MARTIN, ROBERTS, STIRLING (EMC FIT) WITH LAMBDA= 100 MEV C
33 C AND XG --> X**-0.5 AS X--> 0 AT Q0**2 C
35 C MODE 4 CORRESPONDS TO HARRIMAN, C
36 C MARTIN, ROBERTS, STIRLING (BCDMS FIT) WITH LAMBDA= 190 MEV C
37 C AND XG --> CONSTANT AS X--> 0 AT Q0**2 C
43 C (NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION C
44 C IS RETURNED I.E. G(X) = GLU/X ETC, AND THAT "SEA" C
45 C IS THE LIGHT QUARK SEA I.E. UBAR(X)=DBAR(X) C
46 C = SEA/X FOR A PROTON. IF IN DOUBT, CHECK THE C
47 C MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV) C
51 C (THE RANGE OF APPLICABILITY IS CURRENTLY: C
52 C 10**-5 < X < 1 AND 5 < Q**2 < 1.31 * 10**6 C
53 C HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST C
54 C - PROBLEMS, COMMENTS ETC TO WJS@UK.AC.DUR.HEP C
57 C***************************************************************C
58 #include "pdf/impdp.inc"
59 IF(MODE.EQ.1) CALL STRH3E(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
60 IF(MODE.EQ.2) CALL STRH3EP(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
61 IF(MODE.EQ.3) CALL STRH3EM(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
62 IF(MODE.EQ.4) CALL STRH3B(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)