PDF version 8 added

[u/mrichter/AliRoot.git] / PDF / spdf / sfmsfa.F

diff --git a/PDF/spdf/sfmsfa.F b/PDF/spdf/sfmsfa.F
new file mode 100644 (file)
index 0000000..7fe57c0
--- /dev/null
+++ b/PDF/spdf/sfmsfa.F
@@ -0,0 +1,87 @@
+*
+* $Id$
+*
+* $Log$
+* Revision 1.1.1.2  1996/10/30 08:30:27  cernlib
+* Version 7.04
+*
+* Revision 1.1.1.1  1996/04/12 15:29:45  plothow
+* Version 7.01
+*
+*
+#include "pdf/pilot.h"
+      SUBROUTINE SFMSFA(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GLU)
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+c                                                                    c
+c  This program returns the fitted parton distributions for MRS(A)   c
+c      *** Reference: DTP/94/34 (1994) ***                           c
+c                                                                    c
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+#include "pdf/impdp.inc"
+      DIMENSION Z(18,8),F(8)
+      DATA ((Z(I,K),I=1,18),K=1,4)/
+     +2.029386, -1.892760,  0.478114,  0.538000,  0.044742,  0.058742,
+     +3.956000, -0.136659, -0.513630, -0.389400,  1.267936,  0.134397,
+     +5.132000, -0.235874,  0.939374, -0.008433,  0.833784,  0.376702,
+
+     +0.295255, -0.169615,  0.028711,  0.330000,  0.174503, -0.045713,
+     +4.706000, -0.536820,  0.114806,  5.029000, -0.322071,  0.949196,
+     +5.561000, -2.161758, -1.334992, -0.002353,  1.134808, -0.156988,
+
+     +0.775400,-0.517656,0.159759,-0.300000,0.049683,-0.098598,
+     +5.300000,0.233708,0.282415,0.000000,-1.100596,0.052043,
+     +5.200000,-3.232076,0.712254,0.003932,1.625446,-0.699441,
+
+     +0.160750, -0.104311,  0.017157, -0.300000,  0.052451, -0.065844,
+     +9.272000, -0.576575,  0.109982, -1.153000,  1.173228,  0.169775,
+     +15.600000, -6.054861,  0.267608, -0.002159,  0.919959, -0.153501/
+
+      DATA ((Z(I,K),I=1,18),K=5,8)/
+     +0.099280, -0.042174,  0.008380,  0.400000, -0.070009,  0.067714,
+     +9.270000,  0.606514, -0.190625,  0.000000,  7.334181, -3.106807,
+     +25.000000,-20.291985,  4.165453, -0.005437,  0.080749,  0.256036,
+
+     +0.040260, -0.021695,  0.004481, -0.300000,  0.024073, -0.042711,
+     +9.272000, -1.998587,  1.019307, -1.153000,  2.932868, -1.591815,
+     +15.600000,-14.772277,  5.885859, -0.003502,  1.155235, -0.269042,
+
+     +0.001753,  0.017854, -0.004894, -0.195549, -0.180373,  0.047210,
+     +6.592858,  1.960693, -1.082506,  0.121744, -0.775791,  0.007846,
+     +8.249182, -7.279327,  1.702087,  0.928112, -0.034017,  0.084953,
+
+     +-0.004324,  0.012668, -0.001675, -0.276236,  0.090440, -0.112318,
+     +2.647901,  3.239016, -0.534956,  1.390272, -2.629704,  0.646575,
+     +-3.178775,  4.968992, -1.683463,  1.055388,  0.781798, -0.491059/
+C
+      AL=0.23
+      Q2=SCALE*SCALE
+      Q0=2.
+      ARG=( LOG(SCALE/AL)/ LOG(Q0/AL))
+      S= LOG(ARG)
+C
+      DO 10 J=1,8
+
+      A=Z(1,J)+Z(2,J)*S+Z(3,J)*S*S
+      B=Z(4,J)+Z(5,J)*S+Z(6,J)*S*S
+      C=Z(7,J)+Z(8,J)*S+Z(9,J)*S*S
+      D=Z(10,J)+Z(11,J)*S+Z(12,J)*S*S
+      E=Z(13,J)+Z(14,J)*S+Z(15,J)*S*S
+      H=Z(16,J)+Z(17,J)*S+Z(18,J)*S*S
+   10 F(J)=A*X**B*(1.D0-X)**C*(1.D0+D*SQRT(X)+E*X)*( LOG(1.D0/X))**H
+      UPV=F(1)
+      DNV=F(2)
+      GLU=F(3)
+      DUSUM=F(4)
+      DUDIF=F(5)
+      STR=F(6)
+      CHM=F(7)
+      BOT=F(8)
+      DSEA=0.5D0*(DUSUM+DUDIF)
+      USEA=0.5D0*(DUSUM-DUDIF)
+      IF(Q2.LT.2.7D0) CHM=0.D0
+      IF(Q2.LT.30.D0) BOT=0.D0
+      IF(CHM.LT.0.D0) CHM=0.D0
+      IF(BOT.LT.0.D0) BOT=0.D0
+      RETURN
+      END
+