]>
Commit | Line | Data |
---|---|---|
21886bb6 | 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 |