5 * Revision 1.1.1.2 1996/10/30 08:28:53 cernlib
8 * Revision 1.1.1.1 1996/04/12 15:29:25 plothow
12 #include "pdf/pilot.h"
13 SUBROUTINE MRSDSHP(X,SCALE,MODE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GL)
14 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
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
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
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
41 C >>>>>>>> CROSS CHECK <<<<<<<< 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
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
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
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
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,
72 IF(MODE.EQ.3)CALL SFMSBD0(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,
74 IF(MODE.EQ.4)CALL SFMSBDM(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,
76 IF(MODE.EQ.5)CALL SFMSBH(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GL)