]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PDF/spdf/kmrsgeb.F
Separated TOF libraries (base,rec,sim)
[u/mrichter/AliRoot.git] / PDF / spdf / kmrsgeb.F
CommitLineData
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"
13C-----------------------------------------------------------------------
14 SUBROUTINE KMRSGEB(X,SCALE,MODE,UPV,DNV,SEA,STR,CHM,BOT,GL)
15C***************************************************************C
16C C
17C ----- VARIABLE LAMBDA AND GLUONS ---- C
18C C
19C NEW VERSIONS !!!! OCTOBER 1990 C
20C "........................ " A.D. MARTIN, C
21C R.G. ROBERTS AND W.J. STIRLING PREPRINT DTP-90-?? ) C
22C C
23C MODE 1 CORRESPONDS TO KWIECINSKI, C
24C MARTIN, ROBERTS, STIRLING (BCDMS FIT) C
25C WITH LAMBDA(4) = 190 MEV, ETAG = 5.10 C
26C AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2 "B0 FIT" C
27C C
28C MODE 2 CORRESPONDS TO C
29C MARTIN, ROBERTS, STIRLING (BCDMS FIT) C
30C WITH LAMBDA(4) = 135 MEV, ETAG = 4.65 C
31C AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2 "B135 FIT" C
32C C
33C MODE 3 CORRESPONDS TO C
34C MARTIN, ROBERTS, STIRLING (BCDMS FIT) C
35C WITH LAMBDA(4) = 160 MEV, ETAG = 4.25 C
36C AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2 "B160 FIT" C
37C C
38C MODE 4 CORRESPONDS TO C
39C MARTIN, ROBERTS, STIRLING (BCDMS FIT) C
40C WITH LAMBDA(4) = 200 MEV, ETAG = 5.65 C
41C AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2 "B200 FIT" C
42C C
43C MODE 5 CORRESPONDS TO C
44C MARTIN, ROBERTS, STIRLING (BCDMS FIT) C
45C WITH LAMBDA(4) = 235 MEV, ETAG = 5.20 C
46C AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2 "B235 FIT" C
47C C
48C C
49C >>>>>>>> CROSS CHECK <<<<<<<< C
50C C
51C THE FIRST NUMBER IN THE "B0" GRID IS .01727 C
52C THE FIRST NUMBER IN THE "B135" GRID IS .01683 C
53C THE FIRST NUMBER IN THE "B160" GRID IS .01663 C
54C THE FIRST NUMBER IN THE "B200" GRID IS .01601 C
55C THE FIRST NUMBER IN THE "B235" GRID IS .01571 C
56C C
57C C
58C C
59C -*- C
60C C
61C (NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION C
62C IS RETURNED I.E. G(X) = GLU/X ETC, AND THAT "SEA" C
63C IS THE LIGHT QUARK SEA I.E. UBAR(X)=DBAR(X) C
64C = SEA/X FOR A PROTON. IF IN DOUBT, CHECK THE C
65C MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV) C
66C C
67C -*- C
68C C
69C (THE RANGE OF APPLICABILITY IS CURRENTLY: C
70C 10**-5 < X < 1 AND 5 < Q**2 < 1.31 * 10**6 C
71C HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST C
72C - PROBLEMS, COMMENTS ETC TO WJS@UK.AC.DUR.HEP C
73C C
74C C
75C***************************************************************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