PDF/spdf/sfmsfg.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:30:27  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:45  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13       SUBROUTINE SFMSFG(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GLU)
  14 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  15 c                                                                    c
  16 c  This program returns the fitted parton distributions for MRS(G)   c
  17 c  ... it is an analytic approximation to the "grid version" and     c
  18 c  gives much smoother distributions. It is, however, much slower.   c
  19 c  Full details of the fit are contained in the preprint "Pinning    c
  20 c  Down the Gluon in the Proton", by A.D. Martin, R.G. Roberts and   c
  21 c  W.J. Stirling, Rutherford preprint RAL-95-021 (1995), to which    c
  22 c  reference should be made if appropriate.                          c
  23 c                                                                    c
  24 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  25 C
  26 #include "pdf/impdp.inc"
  27 C
  28       DIMENSION Z(18,8),ZZ(18),F(8)
  29       DATA ((Z(I,K),I=1,18),K=1,4) /
  30      +1.665297, -0.590753,  0.086807,  0.676836, -0.027729,  0.006186,
  31      +3.340307,  0.260426,  0.024214, -0.787023,  1.101313, -0.302934,
  32      +6.098057, -2.645910,  0.175731,  0.557278,  0.363394, -0.001549,
  33
  34      +0.171268, -0.086429,  0.012191,  0.388998,  0.139088, -0.014565,
  35      +4.049580, -0.495801,  0.071001,  8.829433,  0.246082, -0.319317,
  36      +5.310913, -6.033850,  1.335430,  0.391294,  0.974633, -0.026356,
  37
  38      +1.121933, -0.980891,  0.250001, -0.279071, -0.118949,  0.001370,
  39      +5.733053,  0.962694, -0.537345, -4.338513,  4.055084, -1.249222,
  40      +11.591153, -9.385871,  2.080363, 0.233747,  0.649513, -0.033386,
  41
  42      +0.811684, -0.732591,  0.188398, -0.146012, -0.450745,  0.162278,
  43      +10.716037,  0.423859, -0.332372, -3.347546, 1.476767, -0.164336,
  44      +10.221619, -2.432161,  0.693669, -0.372843, -0.451730, 0.548307 /
  45
  46       DATA ((Z(I,K),I=1,18),K=5,8) /
  47      +0.018819, -0.019057,  0.014581,  0.473385, -0.026701, -0.191575,
  48      +8.942582, -0.562553,  2.098186, -1.302327, 29.018376,-11.868145,
  49      +105.320622,-53.953297,  1.329035, 0.878000, 0.499894, -0.971650,
  50
  51      +0.093684, -0.074091,  0.033660, -0.107765, -0.408746,  0.100063,
  52      +9.275198,  0.692364,  0.110208, -2.787409,  0.041578,  0.488094,
  53      +10.637883, 3.815478, -4.848576,  0.184162,  0.185983, -0.186796,
  54
  55      +0.001026,  0.015587, -0.001414, -0.111311, -0.346509,  0.094542,
  56      +5.764679,  4.958334, -2.102611, -1.341127, -2.951861,  1.855254,
  57      +12.079860, -0.297841, -3.809543, 1.480855, -0.719267,  0.175987,
  58
  59      +-0.001533,  0.004429, -0.000718, -0.097267, -0.221771, 0.056117,
  60      +3.630148,  2.941006, -0.512863, -0.483095,  1.599856, -1.197675,
  61      +8.366443, -4.004534,  1.022508,  2.526478, -0.713666,  0.277209/
  62
  63       AL=0.254
  64       Q2=SCALE*SCALE
  65       Q0=2.
  66       ARG=( LOG(SCALE/AL)/ LOG(Q0/AL))
  67       S= LOG(ARG)
  68
  69       DO 10 J=1,8
  70
  71       A=Z(1,J)+Z(2,J)*S+Z(3,J)*S*S
  72       B=Z(4,J)+Z(5,J)*S+Z(6,J)*S*S
  73       C=Z(7,J)+Z(8,J)*S+Z(9,J)*S*S
  74       D=Z(10,J)+Z(11,J)*S+Z(12,J)*S*S
  75       E=Z(13,J)+Z(14,J)*S+Z(15,J)*S*S
  76       H=Z(16,J)+Z(17,J)*S+Z(18,J)*S*S
  77    10 F(J)=A*X**B*(1.-X)**C*(1.+D*SQRT(X)+E*X)*( LOG(1./X))**H
  78       UPV=F(1)
  79       DNV=F(2)
  80       GLU=F(3)
  81       DUSUM=F(4)
  82       DUDIF=F(5)
  83       STR=F(6)
  84       CHM=F(7)
  85       BOT=F(8)
  86       DSEA=0.5D0*(DUSUM+DUDIF)
  87       USEA=0.5D0*(DUSUM-DUDIF)
  88       IF(Q2.LT.2.7D0) CHM=0.D0
  89       IF(Q2.LT.30.D0) BOT=0.D0
  90       IF(CHM.LT.0.D0) CHM=0.D0
  91       IF(BOT.LT.0.D0) BOT=0.D0
  92       RETURN
  93       END