]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:29:46 cernlib | |
6 | * Version 7.04 | |
7 | * | |
8 | * Revision 1.1.1.1 1996/04/12 15:29:35 plothow | |
9 | * Version 7.01 | |
10 | * | |
11 | * | |
12 | #include "pdf/pilot.h" | |
13 | C----------------------------------------------------------------------- | |
14 | SUBROUTINE SFEHLQ1(DX,DSCALE,DUP,DDN,DSEA,DSTR,DCHM,DBOT,DTOP,DGL) | |
15 | C ::::::: EHLQ SET 1 SUPERFAST DOUBLE PRECISION ::::::: | |
16 | C | |
17 | #include "pdf/expdp.inc" | |
18 | + DX,DSCALE,DUP,DDN,DSEA,DSTR,DCHM,DBOT,DTOP,DGL | |
19 | PARAMETER (ALAM=0.2, Q02=5.0) | |
20 | #include "pdf/w50511.inc" | |
21 | DIMENSION XQ(9),TX(6),TT(6),TS(6),NEHLQ(8),CEHLQ(6,6,2,8) | |
22 | C | |
23 | C...The following data lines are coefficients needed in the | |
24 | C...Eichten, Hinchliffe, Lane, Quigg proton structure function | |
25 | C...parametrizations, see below. | |
26 | C...Powers of 1-x in different cases. | |
27 | DATA NEHLQ/3,4,7,5,7,7,7,7/ | |
28 | C...Expansion coefficients for up valence quark distribution. | |
29 | DATA (((CEHLQ(IX,IT,NX,1),IX=1,6),IT=1,6),NX=1,2)/ | |
30 | 1 7.677E-01,-2.087E-01,-3.303E-01,-2.517E-02,-1.570E-02,-1.000E-04, | |
31 | 2-5.326E-01,-2.661E-01, 3.201E-01, 1.192E-01, 2.434E-02, 7.620E-03, | |
32 | 3 2.162E-01, 1.881E-01,-8.375E-02,-6.515E-02,-1.743E-02,-5.040E-03, | |
33 | 4-9.211E-02,-9.952E-02, 1.373E-02, 2.506E-02, 8.770E-03, 2.550E-03, | |
34 | 5 3.670E-02, 4.409E-02, 9.600E-04,-7.960E-03,-3.420E-03,-1.050E-03, | |
35 | 6-1.549E-02,-2.026E-02,-3.060E-03, 2.220E-03, 1.240E-03, 4.100E-04, | |
36 | 1 2.395E-01, 2.905E-01, 9.778E-02, 2.149E-02, 3.440E-03, 5.000E-04, | |
37 | 2 1.751E-02,-6.090E-03,-2.687E-02,-1.916E-02,-7.970E-03,-2.750E-03, | |
38 | 3-5.760E-03,-5.040E-03, 1.080E-03, 2.490E-03, 1.530E-03, 7.500E-04, | |
39 | 4 1.740E-03, 1.960E-03, 3.000E-04,-3.400E-04,-2.900E-04,-1.800E-04, | |
40 | 5-5.300E-04,-6.400E-04,-1.700E-04, 4.000E-05, 6.000E-05, 4.000E-05, | |
41 | 6 1.700E-04, 2.200E-04, 8.000E-05, 1.000E-05,-1.000E-05,-1.000E-05/ | |
42 | C...Expansion coefficients for down valence quark distribution. | |
43 | DATA (((CEHLQ(IX,IT,NX,2),IX=1,6),IT=1,6),NX=1,2)/ | |
44 | 1 3.813E-01,-8.090E-02,-1.634E-01,-2.185E-02,-8.430E-03,-6.200E-04, | |
45 | 2-2.948E-01,-1.435E-01, 1.665E-01, 6.638E-02, 1.473E-02, 4.080E-03, | |
46 | 3 1.252E-01, 1.042E-01,-4.722E-02,-3.683E-02,-1.038E-02,-2.860E-03, | |
47 | 4-5.478E-02,-5.678E-02, 8.900E-03, 1.484E-02, 5.340E-03, 1.520E-03, | |
48 | 5 2.220E-02, 2.567E-02,-3.000E-05,-4.970E-03,-2.160E-03,-6.500E-04, | |
49 | 6-9.530E-03,-1.204E-02,-1.510E-03, 1.510E-03, 8.300E-04, 2.700E-04, | |
50 | 1 1.261E-01, 1.354E-01, 3.958E-02, 8.240E-03, 1.660E-03, 4.500E-04, | |
51 | 2 3.890E-03,-1.159E-02,-1.625E-02,-9.610E-03,-3.710E-03,-1.260E-03, | |
52 | 3-1.910E-03,-5.600E-04, 1.590E-03, 1.590E-03, 8.400E-04, 3.900E-04, | |
53 | 4 6.400E-04, 4.900E-04,-1.500E-04,-2.900E-04,-1.800E-04,-1.000E-04, | |
54 | 5-2.000E-04,-1.900E-04, 0.000E+00, 6.000E-05, 4.000E-05, 3.000E-05, | |
55 | 6 7.000E-05, 8.000E-05, 2.000E-05,-1.000E-05,-1.000E-05,-1.000E-05/ | |
56 | C...Expansion coefficients for up and down sea quark distributions. | |
57 | DATA (((CEHLQ(IX,IT,NX,3),IX=1,6),IT=1,6),NX=1,2)/ | |
58 | 1 6.870E-02,-6.861E-02, 2.973E-02,-5.400E-03, 3.780E-03,-9.700E-04, | |
59 | 2-1.802E-02, 1.400E-04, 6.490E-03,-8.540E-03, 1.220E-03,-1.750E-03, | |
60 | 3-4.650E-03, 1.480E-03,-5.930E-03, 6.000E-04,-1.030E-03,-8.000E-05, | |
61 | 4 6.440E-03, 2.570E-03, 2.830E-03, 1.150E-03, 7.100E-04, 3.300E-04, | |
62 | 5-3.930E-03,-2.540E-03,-1.160E-03,-7.700E-04,-3.600E-04,-1.900E-04, | |
63 | 6 2.340E-03, 1.930E-03, 5.300E-04, 3.700E-04, 1.600E-04, 9.000E-05, | |
64 | 1 1.014E+00,-1.106E+00, 3.374E-01,-7.444E-02, 8.850E-03,-8.700E-04, | |
65 | 2 9.233E-01,-1.285E+00, 4.475E-01,-9.786E-02, 1.419E-02,-1.120E-03, | |
66 | 3 4.888E-02,-1.271E-01, 8.606E-02,-2.608E-02, 4.780E-03,-6.000E-04, | |
67 | 4-2.691E-02, 4.887E-02,-1.771E-02, 1.620E-03, 2.500E-04,-6.000E-05, | |
68 | 5 7.040E-03,-1.113E-02, 1.590E-03, 7.000E-04,-2.000E-04, 0.000E+00, | |
69 | 6-1.710E-03, 2.290E-03, 3.800E-04,-3.500E-04, 4.000E-05, 1.000E-05/ | |
70 | C...Expansion coefficients for gluon distribution. | |
71 | DATA (((CEHLQ(IX,IT,NX,4),IX=1,6),IT=1,6),NX=1,2)/ | |
72 | 1 9.482E-01,-9.578E-01, 1.009E-01,-1.051E-01, 3.456E-02,-3.054E-02, | |
73 | 2-9.627E-01, 5.379E-01, 3.368E-01,-9.525E-02, 1.488E-02,-2.051E-02, | |
74 | 3 4.300E-01,-8.306E-02,-3.372E-01, 4.902E-02,-9.160E-03, 1.041E-02, | |
75 | 4-1.925E-01,-1.790E-02, 2.183E-01, 7.490E-03, 4.140E-03,-1.860E-03, | |
76 | 5 8.183E-02, 1.926E-02,-1.072E-01,-1.944E-02,-2.770E-03,-5.200E-04, | |
77 | 6-3.884E-02,-1.234E-02, 5.410E-02, 1.879E-02, 3.350E-03, 1.040E-03, | |
78 | 1 2.948E+01,-3.902E+01, 1.464E+01,-3.335E+00, 5.054E-01,-5.915E-02, | |
79 | 2 2.559E+01,-3.955E+01, 1.661E+01,-4.299E+00, 6.904E-01,-8.243E-02, | |
80 | 3-1.663E+00, 1.176E+00, 1.118E+00,-7.099E-01, 1.948E-01,-2.404E-02, | |
81 | 4-2.168E-01, 8.170E-01,-7.169E-01, 1.851E-01,-1.924E-02,-3.250E-03, | |
82 | 5 2.088E-01,-4.355E-01, 2.239E-01,-2.446E-02,-3.620E-03, 1.910E-03, | |
83 | 6-9.097E-02, 1.601E-01,-5.681E-02,-2.500E-03, 2.580E-03,-4.700E-04/ | |
84 | C...Expansion coefficients for strange sea quark distribution. | |
85 | DATA (((CEHLQ(IX,IT,NX,5),IX=1,6),IT=1,6),NX=1,2)/ | |
86 | 1 4.968E-02,-4.173E-02, 2.102E-02,-3.270E-03, 3.240E-03,-6.700E-04, | |
87 | 2-6.150E-03,-1.294E-02, 6.740E-03,-6.890E-03, 9.000E-04,-1.510E-03, | |
88 | 3-8.580E-03, 5.050E-03,-4.900E-03,-1.600E-04,-9.400E-04,-1.500E-04, | |
89 | 4 7.840E-03, 1.510E-03, 2.220E-03, 1.400E-03, 7.000E-04, 3.500E-04, | |
90 | 5-4.410E-03,-2.220E-03,-8.900E-04,-8.500E-04,-3.600E-04,-2.000E-04, | |
91 | 6 2.520E-03, 1.840E-03, 4.100E-04, 3.900E-04, 1.600E-04, 9.000E-05, | |
92 | 1 9.235E-01,-1.085E+00, 3.464E-01,-7.210E-02, 9.140E-03,-9.100E-04, | |
93 | 2 9.315E-01,-1.274E+00, 4.512E-01,-9.775E-02, 1.380E-02,-1.310E-03, | |
94 | 3 4.739E-02,-1.296E-01, 8.482E-02,-2.642E-02, 4.760E-03,-5.700E-04, | |
95 | 4-2.653E-02, 4.953E-02,-1.735E-02, 1.750E-03, 2.800E-04,-6.000E-05, | |
96 | 5 6.940E-03,-1.132E-02, 1.480E-03, 6.500E-04,-2.100E-04, 0.000E+00, | |
97 | 6-1.680E-03, 2.340E-03, 4.200E-04,-3.400E-04, 5.000E-05, 1.000E-05/ | |
98 | C...Expansion coefficients for charm sea quark distribution. | |
99 | DATA (((CEHLQ(IX,IT,NX,6),IX=1,6),IT=1,6),NX=1,2)/ | |
100 | 1 9.270E-03,-1.817E-02, 9.590E-03,-6.390E-03, 1.690E-03,-1.540E-03, | |
101 | 2 5.710E-03,-1.188E-02, 6.090E-03,-4.650E-03, 1.240E-03,-1.310E-03, | |
102 | 3-3.960E-03, 7.100E-03,-3.590E-03, 1.840E-03,-3.900E-04, 3.400E-04, | |
103 | 4 1.120E-03,-1.960E-03, 1.120E-03,-4.800E-04, 1.000E-04,-4.000E-05, | |
104 | 5 4.000E-05,-3.000E-05,-1.800E-04, 9.000E-05,-5.000E-05,-2.000E-05, | |
105 | 6-4.200E-04, 7.300E-04,-1.600E-04, 5.000E-05, 5.000E-05, 5.000E-05, | |
106 | 1 8.098E-01,-1.042E+00, 3.398E-01,-6.824E-02, 8.760E-03,-9.000E-04, | |
107 | 2 8.961E-01,-1.217E+00, 4.339E-01,-9.287E-02, 1.304E-02,-1.290E-03, | |
108 | 3 3.058E-02,-1.040E-01, 7.604E-02,-2.415E-02, 4.600E-03,-5.000E-04, | |
109 | 4-2.451E-02, 4.432E-02,-1.651E-02, 1.430E-03, 1.200E-04,-1.000E-04, | |
110 | 5 1.122E-02,-1.457E-02, 2.680E-03, 5.800E-04,-1.200E-04, 3.000E-05, | |
111 | 6-7.730E-03, 7.330E-03,-7.600E-04,-2.400E-04, 1.000E-05, 0.000E+00/ | |
112 | C...Expansion coefficients for bottom sea quark distribution. | |
113 | DATA (((CEHLQ(IX,IT,NX,7),IX=1,6),IT=1,6),NX=1,2)/ | |
114 | 1 9.010E-03,-1.401E-02, 7.150E-03,-4.130E-03, 1.260E-03,-1.040E-03, | |
115 | 2 6.280E-03,-9.320E-03, 4.780E-03,-2.890E-03, 9.100E-04,-8.200E-04, | |
116 | 3-2.930E-03, 4.090E-03,-1.890E-03, 7.600E-04,-2.300E-04, 1.400E-04, | |
117 | 4 3.900E-04,-1.200E-03, 4.400E-04,-2.500E-04, 2.000E-05,-2.000E-05, | |
118 | 5 2.600E-04, 1.400E-04,-8.000E-05, 1.000E-04, 1.000E-05, 1.000E-05, | |
119 | 6-2.600E-04, 3.200E-04, 1.000E-05,-1.000E-05, 1.000E-05,-1.000E-05, | |
120 | 1 8.029E-01,-1.075E+00, 3.792E-01,-7.843E-02, 1.007E-02,-1.090E-03, | |
121 | 2 7.903E-01,-1.099E+00, 4.153E-01,-9.301E-02, 1.317E-02,-1.410E-03, | |
122 | 3-1.704E-02,-1.130E-02, 2.882E-02,-1.341E-02, 3.040E-03,-3.600E-04, | |
123 | 4-7.200E-04, 7.230E-03,-5.160E-03, 1.080E-03,-5.000E-05,-4.000E-05, | |
124 | 5 3.050E-03,-4.610E-03, 1.660E-03,-1.300E-04,-1.000E-05, 1.000E-05, | |
125 | 6-4.360E-03, 5.230E-03,-1.610E-03, 2.000E-04,-2.000E-05, 0.000E+00/ | |
126 | C...Expansion coefficients for top sea quark distribution. | |
127 | DATA (((CEHLQ(IX,IT,NX,8),IX=1,6),IT=1,6),NX=1,2)/ | |
128 | 1 4.410E-03,-7.480E-03, 3.770E-03,-2.580E-03, 7.300E-04,-7.100E-04, | |
129 | 2 3.840E-03,-6.050E-03, 3.030E-03,-2.030E-03, 5.800E-04,-5.900E-04, | |
130 | 3-8.800E-04, 1.660E-03,-7.500E-04, 4.700E-04,-1.000E-04, 1.000E-04, | |
131 | 4-8.000E-05,-1.500E-04, 1.200E-04,-9.000E-05, 3.000E-05, 0.000E+00, | |
132 | 5 1.300E-04,-2.200E-04,-2.000E-05,-2.000E-05,-2.000E-05,-2.000E-05, | |
133 | 6-7.000E-05, 1.900E-04,-4.000E-05, 2.000E-05, 0.000E+00, 0.000E+00, | |
134 | 1 6.623E-01,-9.248E-01, 3.519E-01,-7.930E-02, 1.110E-02,-1.180E-03, | |
135 | 2 6.380E-01,-9.062E-01, 3.582E-01,-8.479E-02, 1.265E-02,-1.390E-03, | |
136 | 3-2.581E-02, 2.125E-02, 4.190E-03,-4.980E-03, 1.490E-03,-2.100E-04, | |
137 | 4 7.100E-04, 5.300E-04,-1.270E-03, 3.900E-04,-5.000E-05,-1.000E-05, | |
138 | 5 3.850E-03,-5.060E-03, 1.860E-03,-3.500E-04, 4.000E-05, 0.000E+00, | |
139 | 6-3.530E-03, 4.460E-03,-1.500E-03, 2.700E-04,-3.000E-05, 0.000E+00/ | |
140 | C | |
141 | C | |
142 | C...Proton structure functions from Eichten, Hinchliffe, Lane, Quigg. | |
143 | C...Allowed variable range: 5 GeV^2 < Q^2 < 1E8 GeV^2; 1E-4 < x < 1 | |
144 | C | |
145 | C...Determine Lamdba and x and t expansion variables. | |
146 | X = DX | |
147 | SCALE = DSCALE | |
148 | KF = ABS(NFL) | |
149 | Q2 =SCALE**2 | |
150 | TMIN=LOG(Q02/ALAM**2) | |
151 | TMAX=LOG(1E8/ALAM**2) | |
152 | T=LOG(MAX(1.,Q2/ALAM**2)) | |
153 | VT=MAX(-1.,MIN(1.,(2.*T-TMAX-TMIN)/(TMAX-TMIN))) | |
154 | NX=1 | |
155 | IF(X.LE.0.1) NX=2 | |
156 | IF(NX.EQ.1) VX=(2.*X-1.1)/0.9 | |
157 | IF(NX.EQ.2) VX=MAX(-1.,(2.*LOG(X)+11.51293)/6.90776) | |
158 | CXS=1. | |
159 | ||
160 | C...Chebyshev polynomials for x and t expansion. | |
161 | TX(1)=1. | |
162 | TX(2)=VX | |
163 | TX(3)=2.*VX**2-1. | |
164 | TX(4)=4.*VX**3-3.*VX | |
165 | TX(5)=8.*VX**4-8.*VX**2+1. | |
166 | TX(6)=16.*VX**5-20.*VX**3+5.*VX | |
167 | TT(1)=1. | |
168 | TT(2)=VT | |
169 | TT(3)=2.*VT**2-1. | |
170 | TT(4)=4.*VT**3-3.*VT | |
171 | TT(5)=8.*VT**4-8.*VT**2+1. | |
172 | TT(6)=16.*VT**5-20.*VT**3+5.*VT | |
173 | ||
174 | C...Calculate structure functions. | |
175 | DO 120 KFL=1,6 | |
176 | XQSUM=0. | |
177 | DO 110 IT=1,6 | |
178 | DO 110 IX=1,6 | |
179 | 110 XQSUM=XQSUM+CEHLQ(IX,IT,NX,KFL)*TX(IX)*TT(IT) | |
180 | 120 XQ(KFL)=XQSUM*(1.-X)**NEHLQ(KFL)*CXS | |
181 | C | |
182 | C...Put into output variables | |
183 | DUP = XQ(1) | |
184 | DDN = XQ(2) | |
185 | DSEA = XQ(3) | |
186 | DGL = XQ(4) | |
187 | DSTR = XQ(5) | |
188 | DCHM = XQ(6) | |
189 | DBOT = 0.D0 | |
190 | DTOP = 0.D0 | |
191 | C | |
192 | IF (KF.GE.5) THEN | |
193 | C | |
194 | C...Special expansion for bottom (threshold effects). | |
195 | BOT = 0. | |
196 | TMIN=8.1905 | |
197 | IF(T.LE.TMIN) GOTO 140 | |
198 | VT=MAX(-1.,MIN(1.,(2.*T-TMAX-TMIN)/(TMAX-TMIN))) | |
199 | TT(1)=1. | |
200 | TT(2)=VT | |
201 | TT(3)=2.*VT**2-1. | |
202 | TT(4)=4.*VT**3-3.*VT | |
203 | TT(5)=8.*VT**4-8.*VT**2+1. | |
204 | TT(6)=16.*VT**5-20.*VT**3+5.*VT | |
205 | XQSUM=0. | |
206 | DO 130 IT=1,6 | |
207 | DO 130 IX=1,6 | |
208 | 130 XQSUM=XQSUM+CEHLQ(IX,IT,NX,7)*TX(IX)*TT(IT) | |
209 | BOT=XQSUM*(1.-X)**NEHLQ(7)*CXS | |
210 | 140 CONTINUE | |
211 | ENDIF | |
212 | DBOT = BOT | |
213 | C-- | |
214 | IF (KF.GE.6) THEN | |
215 | C | |
216 | C...Special expansion for top (threshold effects). | |
217 | TOP = 0. | |
218 | TMIN=11.5528 | |
219 | TMIN=TMIN+2.*LOG(TMAS/30.) | |
220 | TMAX=TMAX+2.*LOG(TMAS/30.) | |
221 | IF(T.LE.TMIN) GOTO 160 | |
222 | VT=MAX(-1.,MIN(1.,(2.*T-TMAX-TMIN)/(TMAX-TMIN))) | |
223 | TT(1)=1. | |
224 | TT(2)=VT | |
225 | TT(3)=2.*VT**2-1. | |
226 | TT(4)=4.*VT**3-3.*VT | |
227 | TT(5)=8.*VT**4-8.*VT**2+1. | |
228 | TT(6)=16.*VT**5-20.*VT**3+5.*VT | |
229 | XQSUM=0. | |
230 | DO 150 IT=1,6 | |
231 | DO 150 IX=1,6 | |
232 | 150 XQSUM=XQSUM+CEHLQ(IX,IT,NX,8)*TX(IX)*TT(IT) | |
233 | TOP=XQSUM*(1.-X)**NEHLQ(8)*CXS | |
234 | 160 CONTINUE | |
235 | ENDIF | |
236 | DTOP = TOP | |
237 | C | |
238 | RETURN | |
239 | END |