PDF/spdf/sfmfap.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:30:00  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:38  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13       SUBROUTINE SFMFAP(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(A')  c
  17 c                                                                    c
  18 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  19 C
  20 #include "pdf/impdp.inc"
  21 C
  22       DIMENSION Z(18,8),ZZ(18),F(8)
  23       DATA ((Z(I,K),I=1,18),K=1,4) /
  24      +2.385035, -1.750356,  0.370491,  0.581567,  0.033300,  0.012245,
  25      +3.925280, -0.226187, -0.062425, -0.862021,  1.178930,  0.038682,
  26      +4.894751, -1.377950,  0.280011,  0.061201,  0.749858,  0.092008,
  27
  28      +0.270620, -0.195488,  0.040248,  0.357825,  0.156226, -0.020230,
  29      +4.339434, -0.525551, -0.016963,  6.583414, -0.366604,  1.082098,
  30      +2.174049,  0.263928, -1.336412,  0.101026,  1.140400, -0.026956,
  31
  32      +0.852604, -0.683552,  0.214748, -0.105954, -0.211982,  0.010073,
  33      +4.305152,  0.992946,  0.019575, -1.042052,  0.027041, -0.255378,
  34      +4.441180, -1.800216,  0.129922,  0.639220,  0.833309, -0.341173,
  35
  36      +0.328730, -0.173874,  0.006988, -0.228390,  0.025648, -0.083929,
  37      +9.748830, -0.727386,  0.397301, -1.384145, -0.538920,  0.432677,
  38      +10.440065, -4.280664, 1.752318, -0.192712,  0.860909, -0.191209 /
  39
  40       DATA ((Z(I,K),I=1,18),K=5,8) /
  41      +0.012332, -0.002796,  0.000000, -0.205990,  0.123876,  0.070738,
  42      +10.590572, 0.702136, -0.391293,  0.000000, 16.654520,  8.695223,
  43      +190.873402,-92.774861,-17.858415, -1.078797,0.650973,  0.080318,
  44
  45      +0.042355, -0.025763,  0.006208, -0.117133, -0.090933, -0.031820,
  46      +8.486509, -1.178987,  0.706619, -1.478869,  0.416645,  0.202853,
  47      +12.931922, -4.845551,  0.221202,  0.563591, 0.954011, -0.329131,
  48
  49      +0.001518,  0.018562, -0.005142, -0.067847, -0.212202,  0.017218,
  50      +5.612756,  3.786456, -2.084374, -0.140671, -3.278622,  1.815593,
  51      +6.910856,  0.785627, -3.351047,  1.384535, -0.069658, -0.078327,
  52
  53      +-0.004515, 0.013875, -0.004472, -0.146553, -0.034326, -0.043137,
  54      +2.593083,  3.836694, -0.763188,  1.095127, -2.988716,  1.042826,
  55      +-2.628259,  7.263671, -2.674676, 1.660073,  0.206697, -0.081092/
  56
  57       AL=0.231
  58       Q2=SCALE*SCALE
  59       Q0=2.
  60       ARG=( LOG(SCALE/AL)/ LOG(Q0/AL))
  61       S= LOG(ARG)
  62
  63       DO 10 J=1,8
  64
  65       A=Z(1,J)+Z(2,J)*S+Z(3,J)*S*S
  66       B=Z(4,J)+Z(5,J)*S+Z(6,J)*S*S
  67       C=Z(7,J)+Z(8,J)*S+Z(9,J)*S*S
  68       D=Z(10,J)+Z(11,J)*S+Z(12,J)*S*S
  69       E=Z(13,J)+Z(14,J)*S+Z(15,J)*S*S
  70       H=Z(16,J)+Z(17,J)*S+Z(18,J)*S*S
  71    10 F(J)=A*X**B*(1.-X)**C*(1.+D*SQRT(X)+E*X)*( LOG(1./X))**H
  72       UPV=F(1)
  73       DNV=F(2)
  74       GLU=F(3)
  75       DUSUM=F(4)
  76       DUDIF=F(5)
  77       STR=F(6)
  78       CHM=F(7)
  79       BOT=F(8)
  80       DSEA=0.5D0*(DUSUM+DUDIF)
  81       USEA=0.5D0*(DUSUM-DUDIF)
  82       IF(Q2.LT.2.7D0) CHM=0.D0
  83       IF(Q2.LT.30.D0) BOT=0.D0
  84       IF(CHM.LT.0.D0) CHM=0.D0
  85       IF(BOT.LT.0.D0) BOT=0.D0
  86       RETURN
  87       END