PDF/spdf/smrspi.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:30:53  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:48  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13       SUBROUTINE SMRSPI(X,SCALE,MODE,UPV,DNV,SEA,STR,CHM,BOT,GLU)
  14 C***************************************************************C
  15 C                                                               C
  16 C    PION DISTRIBUTION FUNCTIONS EXTRACTED FROM DRELL-YAN AND   C
  17 C    PROMPT PHOTON PION BEAM DATA, USING (NUCLEAR TARGET        C
  18 C    CORRECTED) HMRS(B) DISTRIBUTIONS FOR THE NUCLEON.          C
  19 C    MAIN UNCERTAINTY IS IN SEA, HENCE FOLLOWING THREE SETS...  C
  20 C    MODE=2 IS THE "BEST FIT" SET.                              C
  21 C                                                               C
  22 C    MODE 1 CORRESPONDS TO A 10% MOMENTUM SEA    (0.11213)      C
  23 C    MODE 2 CORRESPONDS TO A 15% MOMENTUM SEA    (0.16119)      C
  24 C    MODE 3 CORRESPONDS TO A 20% MOMENTUM SEA    (0.23785)      C
  25 C                                                               C
  26 C    (NUMBER IN BRACKETS IS LAST ENTRY IN FIRST ROW OF GRID)    C
  27 C                                                               C
  28 C                         -*-                                   C
  29 C                                                               C
  30 C    REFERENCE: A.D. MARTIN. R.G. ROBERTS. W.J. STIRLING        C
  31 C    AND P.J. SUTTON, PREPRINT RAL-91-058 (1991).               C
  32 C                         -*-                                   C
  33 C                                                               C
  34 C    FOR THESE SETS....                                         C
  35 C    * LAMBDA(MSBAR,NF=4) = 190 MEV                             C
  36 C    * UV=DV, STR=USEA=DSEA, BTM=0                              C
  37 C                                                               C
  38 C    NOTE THAT X TIMES THE PARTON DISTRIBUTION FUNCTION         C
  39 C    IS RETURNED I.E. G(X) = GLU/X ETC, AND THAT "SEA"          C
  40 C    IS THE QUARK SEA / 6   I.E. UBAR(X)=DBAR(X)                C
  41 C    = SEA/X FOR A PION. IF IN DOUBT, CHECK THE                 C
  42 C    MOMENTUM SUM RULE! NOTE ALSO THAT SCALE=Q IN GEV           C
  43 C                                                               C
  44 C                         -*-                                   C
  45 C                                                               C
  46 C     THE RANGE OF APPLICABILITY IS CURRENTLY:                  C
  47 C     10**-5 < X < 1  AND  5 < Q**2 < 1.31 * 10**6              C
  48 C     HIGHER Q**2 VALUES CAN BE SUPPLIED ON REQUEST             C
  49 C     - PROBLEMS, COMMENTS ETC TO WJS@UK.AC.DUR.HEP             C
  50 C                                                               C
  51 C                                                               C
  52 C***************************************************************C
  53 #include "pdf/impdp.inc"
  54       IF(MODE.EQ.1) CALL SMRSP31(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GLU)
  55       IF(MODE.EQ.2) CALL SMRSP32(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GLU)
  56       IF(MODE.EQ.3) CALL SMRSP33(X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GLU)
  57       RETURN
  58       END