spdf/mrsdshp.F

   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