PDF/spdf/mrsdshd.F

   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