]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:28:52 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 | SUBROUTINE MRSDSHD(X,SCALE,MODE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GL) | |
14 | C***************************************************************C | |
15 | C C | |
16 | C C | |
17 | C MRS *** DIS SCHEME ** VERSIONS: S0',D0',D-' (FEB 1993) C | |
18 | C NEW SET MRS(H) ADDED (DECEMBER 1993) C | |
19 | C C | |
20 | C MODE 1 IS THE MRS(S0'-DIS) SET C | |
21 | C MODE 2 IS THE MRS(D0'-DIS) SET C | |
22 | C MODE 3 IS THE MRS(D-'-DIS) SET C | |
23 | C MODE 4 IS THE MRS(H-DIS) SET C | |
24 | C ALL SETS ARE FITTED TO THE RECENT NMC C | |
25 | C AND CCFR STRUCTURE FUNCTION DATA, C | |
26 | C AND HAVE (D0',D-',H) UBAR < DBAR IN THE SEA. C | |
27 | C BOTH THESE NEW SETS HAVE C | |
28 | C LAMBDA(MSbar,NF=4) = 230 MeV C | |
29 | C C | |
30 | C THE REFERENCE IS: A.D. Martin, R.G. Roberts and C | |
31 | C W.J. Stirling, RAL preprint RAL-93-077 (1993) C | |
32 | C C | |
33 | C C | |
34 | C >>>>>>>> CROSS CHECK <<<<<<<< C | |
35 | C C | |
36 | C THE FIRST NUMBER IN THE "1" GRID IS 0.01733 C | |
37 | C THE FIRST NUMBER IN THE "2" GRID IS 0.00930 C | |
38 | C THE FIRST NUMBER IN THE "3" GRID IS 0.00783 C | |
39 | C THE FIRST NUMBER IN THE "4" GRID IS 0.01508 C | |
40 | C C | |
41 | C -*- C | |
42 | C C | |
43 | C (NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION C | |
44 | C IS RETURNED I.E. G(X) = GLU/X ETC. IF IN DOUBT, CHECK THE C | |
45 | C MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV) C | |
46 | C C | |
47 | C -*- C | |
48 | C C | |
49 | C (THE RANGE OF APPLICABILITY IS CURRENTLY: C | |
50 | C 10**-5 < X < 1 AND 5 < Q**2 < 1.31 * 10**6 C | |
51 | C HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST C | |
52 | C - PROBLEMS, COMMENTS ETC TO WJS@UK.AC.DUR.HEP C | |
53 | C C | |
54 | C C | |
55 | C***************************************************************C | |
56 | #include "pdf/impdp.inc" | |
57 | IF(MODE.EQ.1) CALL SFDISS0(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT, | |
58 | +GL) | |
59 | IF(MODE.EQ.2) CALL SFDISD0(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT, | |
60 | +GL) | |
61 | IF(MODE.EQ.3) CALL SFDISDM(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT, | |
62 | +GL) | |
63 | IF(MODE.EQ.4)CALL SFDISH(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GL) | |
64 | RETURN | |
65 | END |