| fe4da5cc |
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 SFMSFA(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 *** Reference: DTP/94/34 (1994) *** c |
| 18 | c c |
| 19 | cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
| 20 | #include "pdf/impdp.inc" |
| 21 | DIMENSION Z(18,8),F(8) |
| 22 | DATA ((Z(I,K),I=1,18),K=1,4)/ |
| 23 | +2.029386, -1.892760, 0.478114, 0.538000, 0.044742, 0.058742, |
| 24 | +3.956000, -0.136659, -0.513630, -0.389400, 1.267936, 0.134397, |
| 25 | +5.132000, -0.235874, 0.939374, -0.008433, 0.833784, 0.376702, |
| 26 | |
| 27 | +0.295255, -0.169615, 0.028711, 0.330000, 0.174503, -0.045713, |
| 28 | +4.706000, -0.536820, 0.114806, 5.029000, -0.322071, 0.949196, |
| 29 | +5.561000, -2.161758, -1.334992, -0.002353, 1.134808, -0.156988, |
| 30 | |
| 31 | +0.775400,-0.517656,0.159759,-0.300000,0.049683,-0.098598, |
| 32 | +5.300000,0.233708,0.282415,0.000000,-1.100596,0.052043, |
| 33 | +5.200000,-3.232076,0.712254,0.003932,1.625446,-0.699441, |
| 34 | |
| 35 | +0.160750, -0.104311, 0.017157, -0.300000, 0.052451, -0.065844, |
| 36 | +9.272000, -0.576575, 0.109982, -1.153000, 1.173228, 0.169775, |
| 37 | +15.600000, -6.054861, 0.267608, -0.002159, 0.919959, -0.153501/ |
| 38 | |
| 39 | DATA ((Z(I,K),I=1,18),K=5,8)/ |
| 40 | +0.099280, -0.042174, 0.008380, 0.400000, -0.070009, 0.067714, |
| 41 | +9.270000, 0.606514, -0.190625, 0.000000, 7.334181, -3.106807, |
| 42 | +25.000000,-20.291985, 4.165453, -0.005437, 0.080749, 0.256036, |
| 43 | |
| 44 | +0.040260, -0.021695, 0.004481, -0.300000, 0.024073, -0.042711, |
| 45 | +9.272000, -1.998587, 1.019307, -1.153000, 2.932868, -1.591815, |
| 46 | +15.600000,-14.772277, 5.885859, -0.003502, 1.155235, -0.269042, |
| 47 | |
| 48 | +0.001753, 0.017854, -0.004894, -0.195549, -0.180373, 0.047210, |
| 49 | +6.592858, 1.960693, -1.082506, 0.121744, -0.775791, 0.007846, |
| 50 | +8.249182, -7.279327, 1.702087, 0.928112, -0.034017, 0.084953, |
| 51 | |
| 52 | +-0.004324, 0.012668, -0.001675, -0.276236, 0.090440, -0.112318, |
| 53 | +2.647901, 3.239016, -0.534956, 1.390272, -2.629704, 0.646575, |
| 54 | +-3.178775, 4.968992, -1.683463, 1.055388, 0.781798, -0.491059/ |
| 55 | C |
| 56 | AL=0.23 |
| 57 | Q2=SCALE*SCALE |
| 58 | Q0=2. |
| 59 | ARG=( LOG(SCALE/AL)/ LOG(Q0/AL)) |
| 60 | S= LOG(ARG) |
| 61 | C |
| 62 | DO 10 J=1,8 |
| 63 | |
| 64 | A=Z(1,J)+Z(2,J)*S+Z(3,J)*S*S |
| 65 | B=Z(4,J)+Z(5,J)*S+Z(6,J)*S*S |
| 66 | C=Z(7,J)+Z(8,J)*S+Z(9,J)*S*S |
| 67 | D=Z(10,J)+Z(11,J)*S+Z(12,J)*S*S |
| 68 | E=Z(13,J)+Z(14,J)*S+Z(15,J)*S*S |
| 69 | H=Z(16,J)+Z(17,J)*S+Z(18,J)*S*S |
| 70 | 10 F(J)=A*X**B*(1.D0-X)**C*(1.D0+D*SQRT(X)+E*X)*( LOG(1.D0/X))**H |
| 71 | UPV=F(1) |
| 72 | DNV=F(2) |
| 73 | GLU=F(3) |
| 74 | DUSUM=F(4) |
| 75 | DUDIF=F(5) |
| 76 | STR=F(6) |
| 77 | CHM=F(7) |
| 78 | BOT=F(8) |
| 79 | DSEA=0.5D0*(DUSUM+DUDIF) |
| 80 | USEA=0.5D0*(DUSUM-DUDIF) |
| 81 | IF(Q2.LT.2.7D0) CHM=0.D0 |
| 82 | IF(Q2.LT.30.D0) BOT=0.D0 |
| 83 | IF(CHM.LT.0.D0) CHM=0.D0 |
| 84 | IF(BOT.LT.0.D0) BOT=0.D0 |
| 85 | RETURN |
| 86 | END |
| 87 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame