]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:28:53 cernlib | |
6 | * Version 7.04 | |
7 | * | |
8 | * Revision 1.1.1.1 1996/04/12 15:29:25 plothow | |
9 | * Version 7.01 | |
10 | * | |
11 | * | |
12 | #include "pdf/pilot.h" | |
13 | SUBROUTINE MRSDSHP(X,SCALE,MODE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GL) | |
14 | C***************************************************************C | |
15 | C C | |
16 | C C | |
17 | C NEW MRS VERSIONS: S0', D0', D-' (NOVEMBER 1992) C | |
18 | C MODE 1 IS THE 1990 KMRS(B0) SET; C | |
19 | C MODES 2-4 ARE NEW SETS FITTED TO THE RECENT NMC C | |
20 | C AND CCFR ***FINAL*** STRUCTURE FUNCTION DATA. C | |
21 | C DIFFERENCES BETWEEN THESE AND THE APRIL 1992 C | |
22 | C VERSIONS ARE SMALL! THE THREE NEW SETS HAVE C | |
23 | C LAMBDA(MSbar,NF=4) = 230 MeV C | |
24 | C C | |
25 | C MODE 5 is a new (October 1993) preliminary "fit" C | |
26 | C to the recent H1 + ZEUS F2 data at small x, C | |
27 | C otherwise similar to D0', D-'. The exponent of the C | |
28 | C gluon and sea at small x for this new set is -0.3 C | |
29 | C (cf. -0.5 for D-') C | |
30 | C C | |
31 | C THE REFERENCE FOR SETS (2-4) IS: C | |
32 | C A.D. Martin, R.G. Roberts and W.J. Stirling, C | |
33 | C Phys. Lett. B306 (1993) 145 C | |
34 | C C | |
35 | C MODE 1: KMRS(B0) (Lambda(4) = 190 MeV) C | |
36 | C MODE 2: MRS(S0') (updated B0, Lambda(4) = 230 MeV) C | |
37 | C MODE 3: MRS(D0') (... but with ubar not= dbar) C | |
38 | C MODE 4: MRS(D-') (updated B-, ubar not= dbar) C | |
39 | C MODE 5: MRS(H) (a preliminary "HERA" fit, Oct. 1993) C | |
40 | C C | |
41 | C >>>>>>>> CROSS CHECK <<<<<<<< C | |
42 | C C | |
43 | C THE FIRST NUMBER IN THE "1" GRID IS 0.01727 C | |
44 | C THE FIRST NUMBER IN THE "2" GRID IS 0.01617 C | |
45 | C THE FIRST NUMBER IN THE "3" GRID IS 0.00820 C | |
46 | C THE FIRST NUMBER IN THE "4" GRID IS 0.00678 C | |
47 | C THE FIRST NUMBER IN THE "5" GRID IS 0.01399 C | |
48 | C C | |
49 | C NOTE THE EXTRA ARGUMENT IN THIS SUBROUTINE MRSEB, C | |
50 | C TO ALLOW FOR THE POSSIBILITY OF A *** DIFFERENT *** C | |
51 | C UBAR AND DBAR SEA! C | |
52 | C C | |
53 | C -*- C | |
54 | C C | |
55 | C (NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION C | |
56 | C IS RETURNED I.E. G(X) = GLU/X ETC. IF IN DOUBT, CHECK THE C | |
57 | C MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV) C | |
58 | C C | |
59 | C -*- C | |
60 | C C | |
61 | C (THE RANGE OF APPLICABILITY IS CURRENTLY: C | |
62 | C 10**-5 < X < 1 AND 5 < Q**2 < 1.31 * 10**6 C | |
63 | C HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST C | |
64 | C - PROBLEMS, COMMENTS ETC TO WJS@UK.AC.DUR.HEP C | |
65 | C C | |
66 | C C | |
67 | C***************************************************************C | |
68 | #include "pdf/impdp.inc" | |
69 | IF(MODE.EQ.1)CALL STRKB0(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GL) | |
70 | IF(MODE.EQ.2)CALL SFMSBS0(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT, | |
71 | +GL) | |
72 | IF(MODE.EQ.3)CALL SFMSBD0(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT, | |
73 | +GL) | |
74 | IF(MODE.EQ.4)CALL SFMSBDM(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT, | |
75 | +GL) | |
76 | IF(MODE.EQ.5)CALL SFMSBH(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GL) | |
77 | RETURN | |
78 | END |