]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
7ef50f50 | 5 | * Revision 1.4 2000/09/18 10:02:36 hristov |
6 | * Makefile added to PDF8 | |
7 | * | |
21886bb6 | 8 | * Revision 1.1.1.2 1996/10/30 08:29:47 cernlib |
9 | * Version 7.04 | |
10 | * | |
11 | * Revision 1.1.1.1 1996/04/12 15:29:35 plothow | |
12 | * Version 7.01 | |
13 | * | |
14 | * | |
15 | #include "pdf/pilot.h" | |
16 | C----------------------------------------------------------------------- | |
17 | SUBROUTINE SFEHLQ2(DX,DSCALE,DUP,DDN,DSEA,DSTR,DCHM,DBOT,DTOP,DGL) | |
18 | C ::::::: EHLQ SET 2 SUPERFAST DOUBLE PRECISION ::::::: | |
19 | C | |
7ef50f50 | 20 | #if defined(CERNLIB_DOUBLE) |
21 | DOUBLE PRECISION | |
22 | #endif | |
23 | #if defined(CERNLIB_SINGLE) | |
24 | REAL | |
25 | #endif | |
21886bb6 | 26 | + DX,DSCALE,DUP,DDN,DSEA,DSTR,DCHM,DBOT,DTOP,DGL |
27 | PARAMETER (ALAM=0.29, Q02=5.0) | |
28 | #include "pdf/w50511.inc" | |
29 | DIMENSION XQ(9),TX(6),TT(6),TS(6),NEHLQ(8),CEHLQ(6,6,2,8) | |
30 | C | |
31 | C...The following data lines are coefficients needed in the | |
32 | C...Eichten, Hinchliffe, Lane, Quigg proton structure function | |
33 | C...parametrizations, see below. | |
34 | C...Powers of 1-x in different cases. | |
35 | DATA NEHLQ/3,4,7,6,7,7,7,7/ | |
36 | C...Expansion coefficients for up valence quark distribution. | |
37 | DATA (((CEHLQ(IX,IT,NX,1),IX=1,6),IT=1,6),NX=1,2)/ | |
38 | 1 7.237E-01,-2.189E-01,-2.995E-01,-1.909E-02,-1.477E-02, 2.500E-04, | |
39 | 2-5.314E-01,-2.425E-01, 3.283E-01, 1.119E-01, 2.223E-02, 7.070E-03, | |
40 | 3 2.289E-01, 1.890E-01,-9.859E-02,-6.900E-02,-1.747E-02,-5.080E-03, | |
41 | 4-1.041E-01,-1.084E-01, 2.108E-02, 2.975E-02, 9.830E-03, 2.830E-03, | |
42 | 5 4.394E-02, 5.116E-02,-1.410E-03,-1.055E-02,-4.230E-03,-1.270E-03, | |
43 | 6-1.991E-02,-2.539E-02,-2.780E-03, 3.430E-03, 1.720E-03, 5.500E-04, | |
44 | 1 2.410E-01, 2.884E-01, 9.369E-02, 1.900E-02, 2.530E-03, 2.400E-04, | |
45 | 2 1.765E-02,-9.220E-03,-3.037E-02,-2.085E-02,-8.440E-03,-2.810E-03, | |
46 | 3-6.450E-03,-5.260E-03, 1.720E-03, 3.110E-03, 1.830E-03, 8.700E-04, | |
47 | 4 2.120E-03, 2.320E-03, 2.600E-04,-4.900E-04,-3.900E-04,-2.300E-04, | |
48 | 5-6.900E-04,-8.200E-04,-2.000E-04, 7.000E-05, 9.000E-05, 6.000E-05, | |
49 | 6 2.400E-04, 3.100E-04, 1.100E-04, 0.000E+00,-2.000E-05,-2.000E-05/ | |
50 | C...Expansion coefficients for down valence quark distribution. | |
51 | DATA (((CEHLQ(IX,IT,NX,2),IX=1,6),IT=1,6),NX=1,2)/ | |
52 | 1 3.578E-01,-8.622E-02,-1.480E-01,-1.840E-02,-7.820E-03,-4.500E-04, | |
53 | 2-2.925E-01,-1.304E-01, 1.696E-01, 6.243E-02, 1.353E-02, 3.750E-03, | |
54 | 3 1.318E-01, 1.041E-01,-5.486E-02,-3.872E-02,-1.038E-02,-2.850E-03, | |
55 | 4-6.162E-02,-6.143E-02, 1.303E-02, 1.740E-02, 5.940E-03, 1.670E-03, | |
56 | 5 2.643E-02, 2.957E-02,-1.490E-03,-6.450E-03,-2.630E-03,-7.700E-04, | |
57 | 6-1.218E-02,-1.497E-02,-1.260E-03, 2.240E-03, 1.120E-03, 3.500E-04, | |
58 | 1 1.263E-01, 1.334E-01, 3.732E-02, 7.070E-03, 1.260E-03, 3.400E-04, | |
59 | 2 3.660E-03,-1.357E-02,-1.795E-02,-1.031E-02,-3.880E-03,-1.280E-03, | |
60 | 3-2.100E-03,-3.600E-04, 2.050E-03, 1.920E-03, 9.800E-04, 4.400E-04, | |
61 | 4 7.700E-04, 5.400E-04,-2.400E-04,-3.900E-04,-2.400E-04,-1.300E-04, | |
62 | 5-2.600E-04,-2.300E-04, 2.000E-05, 9.000E-05, 6.000E-05, 4.000E-05, | |
63 | 6 9.000E-05, 1.000E-04, 2.000E-05,-2.000E-05,-2.000E-05,-1.000E-05/ | |
64 | C...Expansion coefficients for up and down sea quark distributions. | |
65 | DATA (((CEHLQ(IX,IT,NX,3),IX=1,6),IT=1,6),NX=1,2)/ | |
66 | 1 1.008E-01,-7.100E-02, 1.973E-02,-5.710E-03, 2.930E-03,-9.900E-04, | |
67 | 2-5.271E-02,-1.823E-02, 1.792E-02,-6.580E-03, 1.750E-03,-1.550E-03, | |
68 | 3 1.220E-02, 1.763E-02,-8.690E-03,-8.800E-04,-1.160E-03,-2.100E-04, | |
69 | 4-1.190E-03,-7.180E-03, 2.360E-03, 1.890E-03, 7.700E-04, 4.100E-04, | |
70 | 5-9.100E-04, 2.040E-03,-3.100E-04,-1.050E-03,-4.000E-04,-2.400E-04, | |
71 | 6 1.190E-03,-1.700E-04,-2.000E-04, 4.200E-04, 1.700E-04, 1.000E-04, | |
72 | 1 1.081E+00,-1.189E+00, 3.868E-01,-8.617E-02, 1.115E-02,-1.180E-03, | |
73 | 2 9.917E-01,-1.396E+00, 4.998E-01,-1.159E-01, 1.674E-02,-1.720E-03, | |
74 | 3 5.099E-02,-1.338E-01, 9.173E-02,-2.885E-02, 5.890E-03,-6.500E-04, | |
75 | 4-3.178E-02, 5.703E-02,-2.070E-02, 2.440E-03, 1.100E-04,-9.000E-05, | |
76 | 5 8.970E-03,-1.392E-02, 2.050E-03, 6.500E-04,-2.300E-04, 2.000E-05, | |
77 | 6-2.340E-03, 3.010E-03, 5.000E-04,-3.900E-04, 6.000E-05, 1.000E-05/ | |
78 | C...Expansion coefficients for gluon distribution. | |
79 | DATA (((CEHLQ(IX,IT,NX,4),IX=1,6),IT=1,6),NX=1,2)/ | |
80 | 1 2.367E+00, 4.453E-01, 3.660E-01, 9.467E-02, 1.341E-01, 1.661E-02, | |
81 | 2-3.170E+00,-1.795E+00, 3.313E-02,-2.874E-01,-9.827E-02,-7.119E-02, | |
82 | 3 1.823E+00, 1.457E+00,-2.465E-01, 3.739E-02, 6.090E-03, 1.814E-02, | |
83 | 4-1.033E+00,-9.827E-01, 2.136E-01, 1.169E-01, 5.001E-02, 1.684E-02, | |
84 | 5 5.133E-01, 5.259E-01,-1.173E-01,-1.139E-01,-4.988E-02,-2.021E-02, | |
85 | 6-2.881E-01,-3.145E-01, 5.667E-02, 9.161E-02, 4.568E-02, 1.951E-02, | |
86 | 1 3.036E+01,-4.062E+01, 1.578E+01,-3.699E+00, 6.020E-01,-7.031E-02, | |
87 | 2 2.700E+01,-4.167E+01, 1.770E+01,-4.804E+00, 7.862E-01,-1.060E-01, | |
88 | 3-1.909E+00, 1.357E+00, 1.127E+00,-7.181E-01, 2.232E-01,-2.481E-02, | |
89 | 4-2.488E-01, 9.781E-01,-8.127E-01, 2.094E-01,-2.997E-02,-4.710E-03, | |
90 | 5 2.506E-01,-5.427E-01, 2.672E-01,-3.103E-02,-1.800E-03, 2.870E-03, | |
91 | 6-1.128E-01, 2.087E-01,-6.972E-02,-2.480E-03, 2.630E-03,-8.400E-04/ | |
92 | C...Expansion coefficients for strange sea quark distribution. | |
93 | DATA (((CEHLQ(IX,IT,NX,5),IX=1,6),IT=1,6),NX=1,2)/ | |
94 | 1 6.478E-02,-4.537E-02, 1.643E-02,-3.490E-03, 2.710E-03,-6.700E-04, | |
95 | 2-2.223E-02,-2.126E-02, 1.247E-02,-6.290E-03, 1.120E-03,-1.440E-03, | |
96 | 3-1.340E-03, 1.362E-02,-6.130E-03,-7.900E-04,-9.000E-04,-2.000E-04, | |
97 | 4 5.080E-03,-3.610E-03, 1.700E-03, 1.830E-03, 6.800E-04, 4.000E-04, | |
98 | 5-3.580E-03, 6.000E-05,-2.600E-04,-1.050E-03,-3.800E-04,-2.300E-04, | |
99 | 6 2.420E-03, 9.300E-04,-1.000E-04, 4.500E-04, 1.700E-04, 1.100E-04, | |
100 | 1 9.868E-01,-1.171E+00, 3.940E-01,-8.459E-02, 1.124E-02,-1.250E-03, | |
101 | 2 1.001E+00,-1.383E+00, 5.044E-01,-1.152E-01, 1.658E-02,-1.830E-03, | |
102 | 3 4.928E-02,-1.368E-01, 9.021E-02,-2.935E-02, 5.800E-03,-6.600E-04, | |
103 | 4-3.133E-02, 5.785E-02,-2.023E-02, 2.630E-03, 1.600E-04,-8.000E-05, | |
104 | 5 8.840E-03,-1.416E-02, 1.900E-03, 5.800E-04,-2.500E-04, 1.000E-05, | |
105 | 6-2.300E-03, 3.080E-03, 5.500E-04,-3.700E-04, 7.000E-05, 1.000E-05/ | |
106 | C...Expansion coefficients for charm sea quark distribution. | |
107 | DATA (((CEHLQ(IX,IT,NX,6),IX=1,6),IT=1,6),NX=1,2)/ | |
108 | 1 9.980E-03,-1.945E-02, 1.055E-02,-6.870E-03, 1.860E-03,-1.560E-03, | |
109 | 2 5.700E-03,-1.203E-02, 6.250E-03,-4.860E-03, 1.310E-03,-1.370E-03, | |
110 | 3-4.490E-03, 7.990E-03,-4.170E-03, 2.050E-03,-4.400E-04, 3.300E-04, | |
111 | 4 1.470E-03,-2.480E-03, 1.460E-03,-5.700E-04, 1.200E-04,-1.000E-05, | |
112 | 5-9.000E-05, 1.500E-04,-3.200E-04, 1.200E-04,-6.000E-05,-4.000E-05, | |
113 | 6-4.200E-04, 7.600E-04,-1.400E-04, 4.000E-05, 7.000E-05, 5.000E-05, | |
114 | 1 8.698E-01,-1.131E+00, 3.836E-01,-8.111E-02, 1.048E-02,-1.300E-03, | |
115 | 2 9.626E-01,-1.321E+00, 4.854E-01,-1.091E-01, 1.583E-02,-1.700E-03, | |
116 | 3 3.057E-02,-1.088E-01, 8.022E-02,-2.676E-02, 5.590E-03,-5.600E-04, | |
117 | 4-2.845E-02, 5.164E-02,-1.918E-02, 2.210E-03,-4.000E-05,-1.500E-04, | |
118 | 5 1.311E-02,-1.751E-02, 3.310E-03, 5.100E-04,-1.200E-04, 5.000E-05, | |
119 | 6-8.590E-03, 8.380E-03,-9.200E-04,-2.600E-04, 1.000E-05,-1.000E-05/ | |
120 | C...Expansion coefficients for bottom sea quark distribution. | |
121 | DATA (((CEHLQ(IX,IT,NX,7),IX=1,6),IT=1,6),NX=1,2)/ | |
122 | 1 8.980E-03,-1.459E-02, 7.510E-03,-4.410E-03, 1.310E-03,-1.070E-03, | |
123 | 2 5.970E-03,-9.440E-03, 4.800E-03,-3.020E-03, 9.100E-04,-8.500E-04, | |
124 | 3-3.050E-03, 4.440E-03,-2.100E-03, 8.500E-04,-2.400E-04, 1.400E-04, | |
125 | 4 5.300E-04,-1.300E-03, 5.600E-04,-2.700E-04, 3.000E-05,-2.000E-05, | |
126 | 5 2.000E-04, 1.400E-04,-1.100E-04, 1.000E-04, 0.000E+00, 0.000E+00, | |
127 | 6-2.600E-04, 3.200E-04, 0.000E+00,-3.000E-05, 1.000E-05,-1.000E-05, | |
128 | 1 8.672E-01,-1.174E+00, 4.265E-01,-9.252E-02, 1.244E-02,-1.460E-03, | |
129 | 2 8.500E-01,-1.194E+00, 4.630E-01,-1.083E-01, 1.614E-02,-1.830E-03, | |
130 | 3-2.241E-02,-5.630E-03, 2.815E-02,-1.425E-02, 3.520E-03,-4.300E-04, | |
131 | 4-7.300E-04, 8.030E-03,-5.780E-03, 1.380E-03,-1.300E-04,-4.000E-05, | |
132 | 5 3.460E-03,-5.380E-03, 1.960E-03,-2.100E-04, 1.000E-05, 1.000E-05, | |
133 | 6-4.850E-03, 5.950E-03,-1.890E-03, 2.600E-04,-3.000E-05, 0.000E+00/ | |
134 | C...Expansion coefficients for top sea quark distribution. | |
135 | DATA (((CEHLQ(IX,IT,NX,8),IX=1,6),IT=1,6),NX=1,2)/ | |
136 | 1 4.260E-03,-7.530E-03, 3.830E-03,-2.680E-03, 7.600E-04,-7.300E-04, | |
137 | 2 3.640E-03,-6.050E-03, 3.030E-03,-2.090E-03, 5.900E-04,-6.000E-04, | |
138 | 3-9.200E-04, 1.710E-03,-8.200E-04, 5.000E-04,-1.200E-04, 1.000E-04, | |
139 | 4-5.000E-05,-1.600E-04, 1.300E-04,-9.000E-05, 3.000E-05, 0.000E+00, | |
140 | 5 1.300E-04,-2.100E-04,-1.000E-05,-2.000E-05,-2.000E-05,-1.000E-05, | |
141 | 6-8.000E-05, 1.800E-04,-5.000E-05, 2.000E-05, 0.000E+00, 0.000E+00, | |
142 | 1 7.146E-01,-1.007E+00, 3.932E-01,-9.246E-02, 1.366E-02,-1.540E-03, | |
143 | 2 6.856E-01,-9.828E-01, 3.977E-01,-9.795E-02, 1.540E-02,-1.790E-03, | |
144 | 3-3.053E-02, 2.758E-02, 2.150E-03,-4.880E-03, 1.640E-03,-2.500E-04, | |
145 | 4 9.200E-04, 4.200E-04,-1.340E-03, 4.600E-04,-8.000E-05,-1.000E-05, | |
146 | 5 4.230E-03,-5.660E-03, 2.140E-03,-4.300E-04, 6.000E-05, 0.000E+00, | |
147 | 6-3.890E-03, 5.000E-03,-1.740E-03, 3.300E-04,-4.000E-05, 0.000E+00/ | |
148 | C | |
149 | C | |
150 | C...Proton structure functions from Eichten, Hinchliffe, Lane, Quigg. | |
151 | C...Allowed variable range: 5 GeV^2 < Q^2 < 1E8 GeV^2; 1E-4 < x < 1 | |
152 | C | |
153 | C...Determine Lamdba and x and t expansion variables. | |
154 | X = DX | |
155 | SCALE = DSCALE | |
156 | KF = ABS(NFL) | |
157 | Q2 =SCALE*SCALE | |
158 | TMIN=LOG(Q02/ALAM**2) | |
159 | TMAX=LOG(1E8/ALAM**2) | |
160 | T=LOG(MAX(1.,Q2/ALAM**2)) | |
161 | VT=MAX(-1.,MIN(1.,(2.*T-TMAX-TMIN)/(TMAX-TMIN))) | |
162 | NX=1 | |
163 | IF(X.LE.0.1) NX=2 | |
164 | IF(NX.EQ.1) VX=(2.*X-1.1)/0.9 | |
165 | IF(NX.EQ.2) VX=MAX(-1.,(2.*LOG(X)+11.51293)/6.90776) | |
166 | CXS=1. | |
167 | ||
168 | C...Chebyshev polynomials for x and t expansion. | |
169 | TX(1)=1. | |
170 | TX(2)=VX | |
171 | TX(3)=2.*VX**2-1. | |
172 | TX(4)=4.*VX**3-3.*VX | |
173 | TX(5)=8.*VX**4-8.*VX**2+1. | |
174 | TX(6)=16.*VX**5-20.*VX**3+5.*VX | |
175 | TT(1)=1. | |
176 | TT(2)=VT | |
177 | TT(3)=2.*VT**2-1. | |
178 | TT(4)=4.*VT**3-3.*VT | |
179 | TT(5)=8.*VT**4-8.*VT**2+1. | |
180 | TT(6)=16.*VT**5-20.*VT**3+5.*VT | |
181 | ||
182 | C...Calculate structure functions. | |
183 | DO 120 KFL=1,6 | |
184 | XQSUM=0. | |
185 | DO 110 IT=1,6 | |
186 | DO 110 IX=1,6 | |
187 | 110 XQSUM=XQSUM+CEHLQ(IX,IT,NX,KFL)*TX(IX)*TT(IT) | |
188 | 120 XQ(KFL)=XQSUM*(1.-X)**NEHLQ(KFL)*CXS | |
189 | C | |
190 | C...Put into output variables | |
191 | DUP = XQ(1) | |
192 | DDN = XQ(2) | |
193 | DSEA = XQ(3) | |
194 | DGL = XQ(4) | |
195 | DSTR = XQ(5) | |
196 | DCHM = XQ(6) | |
197 | DBOT = 0.D0 | |
198 | DTOP = 0.D0 | |
199 | C | |
200 | IF (KF.GE.5) THEN | |
201 | C | |
202 | C...Special expansion for bottom (threshold effects). | |
203 | BOT = 0. | |
204 | TMIN=7.4474 | |
205 | IF(T.LE.TMIN) GOTO 140 | |
206 | VT=MAX(-1.,MIN(1.,(2.*T-TMAX-TMIN)/(TMAX-TMIN))) | |
207 | TT(1)=1. | |
208 | TT(2)=VT | |
209 | TT(3)=2.*VT**2-1. | |
210 | TT(4)=4.*VT**3-3.*VT | |
211 | TT(5)=8.*VT**4-8.*VT**2+1. | |
212 | TT(6)=16.*VT**5-20.*VT**3+5.*VT | |
213 | XQSUM=0. | |
214 | DO 130 IT=1,6 | |
215 | DO 130 IX=1,6 | |
216 | 130 XQSUM=XQSUM+CEHLQ(IX,IT,NX,7)*TX(IX)*TT(IT) | |
217 | BOT=XQSUM*(1.-X)**NEHLQ(7)*CXS | |
218 | 140 CONTINUE | |
219 | ENDIF | |
220 | DBOT = BOT | |
221 | C-- | |
222 | IF (KF.GE.6) THEN | |
223 | C | |
224 | C...Special expansion for top (threshold effects). | |
225 | TOP = 0. | |
226 | TMIN=10.8097 | |
227 | TMIN=TMIN+2.*LOG(TMAS/30.) | |
228 | TMAX=TMAX+2.*LOG(TMAS/30.) | |
229 | IF(T.LE.TMIN) GOTO 160 | |
230 | VT=MAX(-1.,MIN(1.,(2.*T-TMAX-TMIN)/(TMAX-TMIN))) | |
231 | TT(1)=1. | |
232 | TT(2)=VT | |
233 | TT(3)=2.*VT**2-1. | |
234 | TT(4)=4.*VT**3-3.*VT | |
235 | TT(5)=8.*VT**4-8.*VT**2+1. | |
236 | TT(6)=16.*VT**5-20.*VT**3+5.*VT | |
237 | XQSUM=0. | |
238 | DO 150 IT=1,6 | |
239 | DO 150 IX=1,6 | |
240 | 150 XQSUM=XQSUM+CEHLQ(IX,IT,NX,8)*TX(IX)*TT(IT) | |
241 | TOP=XQSUM*(1.-X)**NEHLQ(8)*CXS | |
242 | 160 CONTINUE | |
243 | ENDIF | |
244 | DTOP = TOP | |
245 | C | |
246 | RETURN | |
247 | END |