| 21886bb6 |
1 | #include "pdf/pilot.h" |
| 2 | SUBROUTINE SFMFR2(X,Q2,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GLU) |
| 3 | cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
| 4 | c c |
| 5 | c This program returns the fitted parton distributions for MRS(R2) c |
| 6 | c c |
| 7 | c Note:: The gluon fit below Q2=2 GeV2 is not reliable c |
| 8 | c c |
| 9 | cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
| 10 | C |
| 11 | #include "pdf/impdp.inc" |
| 12 | C |
| 13 | #include "pdf/expdp.inc" |
| 14 | + F(8) |
| 15 | REAL Z(18,8) |
| 16 | DATA Z/ |
| 17 | .1.840938, -0.613912, 0.069060, 0.664453, -0.047312, 0.008339, |
| 18 | .3.199642, 0.299148, -0.005252, -1.176342, 1.452607, -0.242739, |
| 19 | .7.607104, -3.009800, 0.110586, 0.305600, 0.324526, 0.004361, |
| 20 | |
| 21 | .0.116204, -0.036926, 0.005132, 0.347559, 0.044734, 0.025428, |
| 22 | .3.776571, 0.129517, -0.106180, 4.583171, 6.631731, -2.496486, |
| 23 | .29.786677,-19.577891, 3.120062, 0.410322, 0.467855, 0.101884, |
| 24 | |
| 25 | .9.333184, -8.070131, 2.146747, -0.089195, -0.436328, 0.122834, |
| 26 | .6.948321, -0.174262, 0.921887, -4.043493, 2.187837, -1.088141, |
| 27 | .7.196254, -4.912953, 2.464480, -1.359569, 0.705893, -0.012063, |
| 28 | |
| 29 | .0.083723, 0.102184, -0.032540, -0.069562, -0.107065, -0.042555, |
| 30 | .6.848496, 2.149966, -0.330963, 2.365100, -5.254575, 1.846434, |
| 31 | .11.934685, -0.271707, -1.758501, 0.554814, 0.050133, -0.160565, |
| 32 | |
| 33 | .-0.221E-04, 0.120E-04,-0.234E-05,-0.603E+00, 0.251E+00,-0.146E-01, |
| 34 | . 0.108E+02,-0.722E+00, 0.876E-02,-0.774E+04,-0.783E+03, 0.152E+04, |
| 35 | .-0.125E+06, 0.519E+05,-0.507E+04,-0.244E+01, 0.146E+01, 0.237E-01, |
| 36 | |
| 37 | .0.285E-01, 0.481E-01,-0.128E-01,-0.686E-01,-0.139E+00,-0.255E-01, |
| 38 | .0.694E+01, 0.204E+01,-0.290E+00, 0.105E+01,-0.445E+01, 0.154E+01, |
| 39 | .0.988E+01, 0.519E+00,-0.159E+01, 0.447E+00, 0.123E+00,-0.158E+00, |
| 40 | |
| 41 | .0.536E-02,-0.149E-01, 0.337E-01,-0.136E+00, 0.397E-01,-0.114E+00, |
| 42 | .0.459E+01, 0.333E+01,-0.115E+00,-0.309E+01, 0.189E+00,-0.109E+00, |
| 43 | .0.182E+02,-0.122E+02, 0.286E+01,-0.323E+00, 0.207E+01,-0.108E+01, |
| 44 | |
| 45 | .-0.111E-01, 0.134E-01,-0.146E-02,-0.130E+00,-0.157E-01,-0.418E-01, |
| 46 | . 0.276E+01, 0.277E+01,-0.697E-01, 0.124E+01,-0.255E+01, 0.466E+00, |
| 47 | .-0.125E+01, 0.581E+01,-0.154E+01, 0.148E+01, 0.353E+00,-0.196E+00/ |
| 48 | DATA ZEROD/0.D0/,PONED/0.1D0/,ONED/1.D0/,ONEDO/1.1D0/,TWOD/2.D0/ |
| 49 | C |
| 50 | UPV= ZEROD |
| 51 | DNV= ZEROD |
| 52 | USEA=ZEROD |
| 53 | DSEA=ZEROD |
| 54 | STR= ZEROD |
| 55 | CHM= ZEROD |
| 56 | BOT= ZEROD |
| 57 | GLU= ZEROD |
| 58 | IF(X.GE.ONED) RETURN |
| 59 | C |
| 60 | AL=0.344D0 |
| 61 | AL2=AL*AL |
| 62 | Q02=ONED |
| 63 | ARG=( LOG(Q2/AL2)/ LOG(Q02/AL2)) |
| 64 | S= LOG(ARG) |
| 65 | C |
| 66 | DO 10 J=1,8 |
| 67 | C |
| 68 | A=Z(1,J)+Z(2,J)*S+Z(3,J)*S*S |
| 69 | B=Z(4,J)+Z(5,J)*S+Z(6,J)*S*S |
| 70 | C=Z(7,J)+Z(8,J)*S+Z(9,J)*S*S |
| 71 | D=Z(10,J)+Z(11,J)*S+Z(12,J)*S*S |
| 72 | E=Z(13,J)+Z(14,J)*S+Z(15,J)*S*S |
| 73 | H=Z(16,J)+Z(17,J)*S+Z(18,J)*S*S |
| 74 | 10 F(J)=A*X**B*(ONED-X)**C*(ONED+D*SQRT(X)+E*X)*( LOG(ONED/X))**H |
| 75 | UPV=F(1) |
| 76 | DNV=F(2) |
| 77 | GLU=F(3) |
| 78 | DUSUM=F(4) |
| 79 | DUDIF=F(5) |
| 80 | STR=F(6) |
| 81 | CHM=F(7) |
| 82 | BOT=F(8) |
| 83 | DSEA=0.5D0*(DUSUM+DUDIF) |
| 84 | USEA=0.5D0*(DUSUM-DUDIF) |
| 85 | IF(Q2.LT.2.7D0) CHM=ZEROD |
| 86 | IF(Q2.LT.30.D0) BOT=ZEROD |
| 87 | IF(CHM.LT.ZEROD) CHM=ZEROD |
| 88 | IF(BOT.LT.ZEROD) BOT=ZEROD |
| 89 | RETURN |
| 90 | END |
git://git.uio.no / u / mrichter / AliRoot.git / blame