]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PDF/spdf/kmrseb.F
Possibility to have different binaries in the same tree introduced
[u/mrichter/AliRoot.git] / PDF / spdf / kmrseb.F
CommitLineData
fe4da5cc 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.2 1996/10/30 08:28:50 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 KMRSEB(X,SCALE,MODE,UPV,DNV,SEA,STR,CHM,BOT,GL)
15C***************************************************************C
16C C
17C ----- VARIABLE QUARKS AND GLUONS AT SMALL X ---- C
18C C
19C NEW VERSIONS !!!! JULY 1990 C
20C "........................ " J. KWIECINSKI, A.D. MARTIN, C
21C R.G. ROBERTS AND W.J. STIRLING PREPRINT DTP-90-46 ) C
22C C
23C MODE 1 CORRESPONDS TO HARRIMAN, C
24C MARTIN, ROBERTS, STIRLING (BCDMS FIT) WITH LAMBDA= 190 MEV C
25C WITH SMALL X BEHAVIOUR DETERMINED FROM FIT "HB FIT" C
26C C
27C MODE 2 CORRESPONDS TO KWIECINSKI, C
28C MARTIN, ROBERTS, STIRLING (BCDMS FIT) WITH LAMBDA= 190 MEV C
29C AND XG,XQ --> CONSTANT AS X--> 0 AT Q0**2 "B0 FIT" C
30C C
31C MODE 3 CORRESPONDS TO KWIECINSKI, C
32C MARTIN, ROBERTS, STIRLING (BCDMS FIT) WITH LAMBDA= 190 MEV C
33C AND XG,XQ --> X**-1/2 AS X--> 0 AT Q0**2 "B- FIT" C
34C C
35C MODE 4 CORRESPONDS TO KWIECINSKI, C
36C MARTIN, ROBERTS, STIRLING (BCDMS FIT) WITH LAMBDA= 190 MEV C
37C AND XG,XQ --> X**-1/2 AS X--> 0 AT Q0**2 "B-(5) FIT" C
38C I.E. WITH WEAK (R=5 GEV^-1) SHADOWING INCLUDED C
39C C
40C MODE 5 CORRESPONDS TO KWIECINSKI, C
41C MARTIN, ROBERTS, STIRLING (BCDMS FIT) WITH LAMBDA= 190 MEV C
42C AND XG,XQ --> X**-1/2 AS X--> 0 AT Q0**2 "B-(2) FIT" C
43C I.E. WITH STRONG (R=2 GEV^-1) SHADOWING INCLUDED C
44C C
45C C
46C >>>>>>>> CROSS CHECK <<<<<<<< C
47C C
48C THE FIRST NUMBER IN THE "HB" GRID IS .03058 C
49C THE FIRST NUMBER IN THE "B0" GRID IS .01727 C
50C THE FIRST NUMBER IN THE "B-" GRID IS .01543 C
51C THE FIRST NUMBER IN THE "B-5"GRID IS .0 C
52C THE FIRST NUMBER IN THE "B-2"GRID IS .0 C
53C C
54C C
55C C
56C -*- C
57C C
58C (NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION C
59C IS RETURNED I.E. G(X) = GLU/X ETC, AND THAT "SEA" C
60C IS THE LIGHT QUARK SEA I.E. UBAR(X)=DBAR(X) C
61C = SEA/X FOR A PROTON. IF IN DOUBT, CHECK THE C
62C MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV) C
63C C
64C -*- C
65C C
66C (THE RANGE OF APPLICABILITY IS CURRENTLY: C
67C 10**-5 < X < 1 AND 5 < Q**2 < 1.31 * 10**6 C
68C HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST C
69C - PROBLEMS, COMMENTS ETC TO WJS@UK.AC.DUR.HEP C
70C C
71C C
72C***************************************************************C
73#include "pdf/impdp.inc"
74 IF(MODE.EQ.1) CALL STRH3B(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
75 IF(MODE.EQ.2) CALL STRKB0O(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
76 IF(MODE.EQ.3) CALL STRKBM(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
77 IF(MODE.EQ.5) CALL SFKBMR2(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
78 IF(MODE.EQ.5) CALL SFKBMR5(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GL)
79 RETURN
80 END