]>
Commit | Line | Data |
---|---|---|
21886bb6 | 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 SFMSFG(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(G) c | |
17 | c ... it is an analytic approximation to the "grid version" and c | |
18 | c gives much smoother distributions. It is, however, much slower. c | |
19 | c Full details of the fit are contained in the preprint "Pinning c | |
20 | c Down the Gluon in the Proton", by A.D. Martin, R.G. Roberts and c | |
21 | c W.J. Stirling, Rutherford preprint RAL-95-021 (1995), to which c | |
22 | c reference should be made if appropriate. c | |
23 | c c | |
24 | cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc | |
25 | C | |
26 | #include "pdf/impdp.inc" | |
27 | C | |
28 | DIMENSION Z(18,8),ZZ(18),F(8) | |
29 | DATA ((Z(I,K),I=1,18),K=1,4) / | |
30 | +1.665297, -0.590753, 0.086807, 0.676836, -0.027729, 0.006186, | |
31 | +3.340307, 0.260426, 0.024214, -0.787023, 1.101313, -0.302934, | |
32 | +6.098057, -2.645910, 0.175731, 0.557278, 0.363394, -0.001549, | |
33 | ||
34 | +0.171268, -0.086429, 0.012191, 0.388998, 0.139088, -0.014565, | |
35 | +4.049580, -0.495801, 0.071001, 8.829433, 0.246082, -0.319317, | |
36 | +5.310913, -6.033850, 1.335430, 0.391294, 0.974633, -0.026356, | |
37 | ||
38 | +1.121933, -0.980891, 0.250001, -0.279071, -0.118949, 0.001370, | |
39 | +5.733053, 0.962694, -0.537345, -4.338513, 4.055084, -1.249222, | |
40 | +11.591153, -9.385871, 2.080363, 0.233747, 0.649513, -0.033386, | |
41 | ||
42 | +0.811684, -0.732591, 0.188398, -0.146012, -0.450745, 0.162278, | |
43 | +10.716037, 0.423859, -0.332372, -3.347546, 1.476767, -0.164336, | |
44 | +10.221619, -2.432161, 0.693669, -0.372843, -0.451730, 0.548307 / | |
45 | ||
46 | DATA ((Z(I,K),I=1,18),K=5,8) / | |
47 | +0.018819, -0.019057, 0.014581, 0.473385, -0.026701, -0.191575, | |
48 | +8.942582, -0.562553, 2.098186, -1.302327, 29.018376,-11.868145, | |
49 | +105.320622,-53.953297, 1.329035, 0.878000, 0.499894, -0.971650, | |
50 | ||
51 | +0.093684, -0.074091, 0.033660, -0.107765, -0.408746, 0.100063, | |
52 | +9.275198, 0.692364, 0.110208, -2.787409, 0.041578, 0.488094, | |
53 | +10.637883, 3.815478, -4.848576, 0.184162, 0.185983, -0.186796, | |
54 | ||
55 | +0.001026, 0.015587, -0.001414, -0.111311, -0.346509, 0.094542, | |
56 | +5.764679, 4.958334, -2.102611, -1.341127, -2.951861, 1.855254, | |
57 | +12.079860, -0.297841, -3.809543, 1.480855, -0.719267, 0.175987, | |
58 | ||
59 | +-0.001533, 0.004429, -0.000718, -0.097267, -0.221771, 0.056117, | |
60 | +3.630148, 2.941006, -0.512863, -0.483095, 1.599856, -1.197675, | |
61 | +8.366443, -4.004534, 1.022508, 2.526478, -0.713666, 0.277209/ | |
62 | ||
63 | AL=0.254 | |
64 | Q2=SCALE*SCALE | |
65 | Q0=2. | |
66 | ARG=( LOG(SCALE/AL)/ LOG(Q0/AL)) | |
67 | S= LOG(ARG) | |
68 | ||
69 | DO 10 J=1,8 | |
70 | ||
71 | A=Z(1,J)+Z(2,J)*S+Z(3,J)*S*S | |
72 | B=Z(4,J)+Z(5,J)*S+Z(6,J)*S*S | |
73 | C=Z(7,J)+Z(8,J)*S+Z(9,J)*S*S | |
74 | D=Z(10,J)+Z(11,J)*S+Z(12,J)*S*S | |
75 | E=Z(13,J)+Z(14,J)*S+Z(15,J)*S*S | |
76 | H=Z(16,J)+Z(17,J)*S+Z(18,J)*S*S | |
77 | 10 F(J)=A*X**B*(1.-X)**C*(1.+D*SQRT(X)+E*X)*( LOG(1./X))**H | |
78 | UPV=F(1) | |
79 | DNV=F(2) | |
80 | GLU=F(3) | |
81 | DUSUM=F(4) | |
82 | DUDIF=F(5) | |
83 | STR=F(6) | |
84 | CHM=F(7) | |
85 | BOT=F(8) | |
86 | DSEA=0.5D0*(DUSUM+DUDIF) | |
87 | USEA=0.5D0*(DUSUM-DUDIF) | |
88 | IF(Q2.LT.2.7D0) CHM=0.D0 | |
89 | IF(Q2.LT.30.D0) BOT=0.D0 | |
90 | IF(CHM.LT.0.D0) CHM=0.D0 | |
91 | IF(BOT.LT.0.D0) BOT=0.D0 | |
92 | RETURN | |
93 | END |