| 21886bb6 |
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 |
git://git.uio.no / u / mrichter / AliRoot.git / blame