Merging the VirtualMC branch to the main development branch (HEAD)
[u/mrichter/AliRoot.git] / HIJING / hipyset1_35 / pystfu_hijing.F
CommitLineData
e74335a4 1* $Id$
2
3C*********************************************************************
4
5 SUBROUTINE PYSTFU_HIJING(KF,X,Q2,XPQ,JBT)
6
7C *******JBT specifies beam or target of the particle
8C...Gives proton and pi+ parton structure functions according to a few
9C...different parametrizations. Note that what is coded is x times the
10C...probability distribution, i.e. xq(x,Q2) etc.
11#include "hiparnt.inc"
12#include "hijcrdn.inc"
13C ********COMMON BLOCK FROM HIJING
14#include "ludat1_hijing.inc"
15#include "ludat2_hijing.inc"
16#include "pypars_hijing.inc"
17#include "pyint1_hijing.inc"
18 DIMENSION XPQ(-6:6),XQ(6),TX(6),TT(6),TS(6),NEHLQ(8,2),
19 &CEHLQ(6,6,2,8,2),CDO(3,6,5,2),COW(3,5,4,2)
20
21C...The following data lines are coefficients needed in the
22C...Eichten, Hinchliffe, Lane, Quigg proton structure function
23C...parametrizations, see below.
24C...Powers of 1-x in different cases.
25 DATA NEHLQ/3,4,7,5,7,7,7,7,3,4,7,6,7,7,7,7/
26C...Expansion coefficients for up valence quark distribution.
27 DATA (((CEHLQ(IX,IT,NX,1,1),IX=1,6),IT=1,6),NX=1,2)/
28 1 7.677E-01,-2.087E-01,-3.303E-01,-2.517E-02,-1.570E-02,-1.000E-04,
29 2-5.326E-01,-2.661E-01, 3.201E-01, 1.192E-01, 2.434E-02, 7.620E-03,
30 3 2.162E-01, 1.881E-01,-8.375E-02,-6.515E-02,-1.743E-02,-5.040E-03,
31 4-9.211E-02,-9.952E-02, 1.373E-02, 2.506E-02, 8.770E-03, 2.550E-03,
32 5 3.670E-02, 4.409E-02, 9.600E-04,-7.960E-03,-3.420E-03,-1.050E-03,
33 6-1.549E-02,-2.026E-02,-3.060E-03, 2.220E-03, 1.240E-03, 4.100E-04,
34 1 2.395E-01, 2.905E-01, 9.778E-02, 2.149E-02, 3.440E-03, 5.000E-04,
35 2 1.751E-02,-6.090E-03,-2.687E-02,-1.916E-02,-7.970E-03,-2.750E-03,
36 3-5.760E-03,-5.040E-03, 1.080E-03, 2.490E-03, 1.530E-03, 7.500E-04,
37 4 1.740E-03, 1.960E-03, 3.000E-04,-3.400E-04,-2.900E-04,-1.800E-04,
38 5-5.300E-04,-6.400E-04,-1.700E-04, 4.000E-05, 6.000E-05, 4.000E-05,
39 6 1.700E-04, 2.200E-04, 8.000E-05, 1.000E-05,-1.000E-05,-1.000E-05/
40 DATA (((CEHLQ(IX,IT,NX,1,2),IX=1,6),IT=1,6),NX=1,2)/
41 1 7.237E-01,-2.189E-01,-2.995E-01,-1.909E-02,-1.477E-02, 2.500E-04,
42 2-5.314E-01,-2.425E-01, 3.283E-01, 1.119E-01, 2.223E-02, 7.070E-03,
43 3 2.289E-01, 1.890E-01,-9.859E-02,-6.900E-02,-1.747E-02,-5.080E-03,
44 4-1.041E-01,-1.084E-01, 2.108E-02, 2.975E-02, 9.830E-03, 2.830E-03,
45 5 4.394E-02, 5.116E-02,-1.410E-03,-1.055E-02,-4.230E-03,-1.270E-03,
46 6-1.991E-02,-2.539E-02,-2.780E-03, 3.430E-03, 1.720E-03, 5.500E-04,
47 1 2.410E-01, 2.884E-01, 9.369E-02, 1.900E-02, 2.530E-03, 2.400E-04,
48 2 1.765E-02,-9.220E-03,-3.037E-02,-2.085E-02,-8.440E-03,-2.810E-03,
49 3-6.450E-03,-5.260E-03, 1.720E-03, 3.110E-03, 1.830E-03, 8.700E-04,
50 4 2.120E-03, 2.320E-03, 2.600E-04,-4.900E-04,-3.900E-04,-2.300E-04,
51 5-6.900E-04,-8.200E-04,-2.000E-04, 7.000E-05, 9.000E-05, 6.000E-05,
52 6 2.400E-04, 3.100E-04, 1.100E-04, 0.000E+00,-2.000E-05,-2.000E-05/
53C...Expansion coefficients for down valence quark distribution.
54 DATA (((CEHLQ(IX,IT,NX,2,1),IX=1,6),IT=1,6),NX=1,2)/
55 1 3.813E-01,-8.090E-02,-1.634E-01,-2.185E-02,-8.430E-03,-6.200E-04,
56 2-2.948E-01,-1.435E-01, 1.665E-01, 6.638E-02, 1.473E-02, 4.080E-03,
57 3 1.252E-01, 1.042E-01,-4.722E-02,-3.683E-02,-1.038E-02,-2.860E-03,
58 4-5.478E-02,-5.678E-02, 8.900E-03, 1.484E-02, 5.340E-03, 1.520E-03,
59 5 2.220E-02, 2.567E-02,-3.000E-05,-4.970E-03,-2.160E-03,-6.500E-04,
60 6-9.530E-03,-1.204E-02,-1.510E-03, 1.510E-03, 8.300E-04, 2.700E-04,
61 1 1.261E-01, 1.354E-01, 3.958E-02, 8.240E-03, 1.660E-03, 4.500E-04,
62 2 3.890E-03,-1.159E-02,-1.625E-02,-9.610E-03,-3.710E-03,-1.260E-03,
63 3-1.910E-03,-5.600E-04, 1.590E-03, 1.590E-03, 8.400E-04, 3.900E-04,
64 4 6.400E-04, 4.900E-04,-1.500E-04,-2.900E-04,-1.800E-04,-1.000E-04,
65 5-2.000E-04,-1.900E-04, 0.000E+00, 6.000E-05, 4.000E-05, 3.000E-05,
66 6 7.000E-05, 8.000E-05, 2.000E-05,-1.000E-05,-1.000E-05,-1.000E-05/
67 DATA (((CEHLQ(IX,IT,NX,2,2),IX=1,6),IT=1,6),NX=1,2)/
68 1 3.578E-01,-8.622E-02,-1.480E-01,-1.840E-02,-7.820E-03,-4.500E-04,
69 2-2.925E-01,-1.304E-01, 1.696E-01, 6.243E-02, 1.353E-02, 3.750E-03,
70 3 1.318E-01, 1.041E-01,-5.486E-02,-3.872E-02,-1.038E-02,-2.850E-03,
71 4-6.162E-02,-6.143E-02, 1.303E-02, 1.740E-02, 5.940E-03, 1.670E-03,
72 5 2.643E-02, 2.957E-02,-1.490E-03,-6.450E-03,-2.630E-03,-7.700E-04,
73 6-1.218E-02,-1.497E-02,-1.260E-03, 2.240E-03, 1.120E-03, 3.500E-04,
74 1 1.263E-01, 1.334E-01, 3.732E-02, 7.070E-03, 1.260E-03, 3.400E-04,
75 2 3.660E-03,-1.357E-02,-1.795E-02,-1.031E-02,-3.880E-03,-1.280E-03,
76 3-2.100E-03,-3.600E-04, 2.050E-03, 1.920E-03, 9.800E-04, 4.400E-04,
77 4 7.700E-04, 5.400E-04,-2.400E-04,-3.900E-04,-2.400E-04,-1.300E-04,
78 5-2.600E-04,-2.300E-04, 2.000E-05, 9.000E-05, 6.000E-05, 4.000E-05,
79 6 9.000E-05, 1.000E-04, 2.000E-05,-2.000E-05,-2.000E-05,-1.000E-05/
80C...Expansion coefficients for up and down sea quark distributions.
81 DATA (((CEHLQ(IX,IT,NX,3,1),IX=1,6),IT=1,6),NX=1,2)/
82 1 6.870E-02,-6.861E-02, 2.973E-02,-5.400E-03, 3.780E-03,-9.700E-04,
83 2-1.802E-02, 1.400E-04, 6.490E-03,-8.540E-03, 1.220E-03,-1.750E-03,
84 3-4.650E-03, 1.480E-03,-5.930E-03, 6.000E-04,-1.030E-03,-8.000E-05,
85 4 6.440E-03, 2.570E-03, 2.830E-03, 1.150E-03, 7.100E-04, 3.300E-04,
86 5-3.930E-03,-2.540E-03,-1.160E-03,-7.700E-04,-3.600E-04,-1.900E-04,
87 6 2.340E-03, 1.930E-03, 5.300E-04, 3.700E-04, 1.600E-04, 9.000E-05,
88 1 1.014E+00,-1.106E+00, 3.374E-01,-7.444E-02, 8.850E-03,-8.700E-04,
89 2 9.233E-01,-1.285E+00, 4.475E-01,-9.786E-02, 1.419E-02,-1.120E-03,
90 3 4.888E-02,-1.271E-01, 8.606E-02,-2.608E-02, 4.780E-03,-6.000E-04,
91 4-2.691E-02, 4.887E-02,-1.771E-02, 1.620E-03, 2.500E-04,-6.000E-05,
92 5 7.040E-03,-1.113E-02, 1.590E-03, 7.000E-04,-2.000E-04, 0.000E+00,
93 6-1.710E-03, 2.290E-03, 3.800E-04,-3.500E-04, 4.000E-05, 1.000E-05/
94 DATA (((CEHLQ(IX,IT,NX,3,2),IX=1,6),IT=1,6),NX=1,2)/
95 1 1.008E-01,-7.100E-02, 1.973E-02,-5.710E-03, 2.930E-03,-9.900E-04,
96 2-5.271E-02,-1.823E-02, 1.792E-02,-6.580E-03, 1.750E-03,-1.550E-03,
97 3 1.220E-02, 1.763E-02,-8.690E-03,-8.800E-04,-1.160E-03,-2.100E-04,
98 4-1.190E-03,-7.180E-03, 2.360E-03, 1.890E-03, 7.700E-04, 4.100E-04,
99 5-9.100E-04, 2.040E-03,-3.100E-04,-1.050E-03,-4.000E-04,-2.400E-04,
100 6 1.190E-03,-1.700E-04,-2.000E-04, 4.200E-04, 1.700E-04, 1.000E-04,
101 1 1.081E+00,-1.189E+00, 3.868E-01,-8.617E-02, 1.115E-02,-1.180E-03,
102 2 9.917E-01,-1.396E+00, 4.998E-01,-1.159E-01, 1.674E-02,-1.720E-03,
103 3 5.099E-02,-1.338E-01, 9.173E-02,-2.885E-02, 5.890E-03,-6.500E-04,
104 4-3.178E-02, 5.703E-02,-2.070E-02, 2.440E-03, 1.100E-04,-9.000E-05,
105 5 8.970E-03,-1.392E-02, 2.050E-03, 6.500E-04,-2.300E-04, 2.000E-05,
106 6-2.340E-03, 3.010E-03, 5.000E-04,-3.900E-04, 6.000E-05, 1.000E-05/
107C...Expansion coefficients for gluon distribution.
108 DATA (((CEHLQ(IX,IT,NX,4,1),IX=1,6),IT=1,6),NX=1,2)/
109 1 9.482E-01,-9.578E-01, 1.009E-01,-1.051E-01, 3.456E-02,-3.054E-02,
110 2-9.627E-01, 5.379E-01, 3.368E-01,-9.525E-02, 1.488E-02,-2.051E-02,
111 3 4.300E-01,-8.306E-02,-3.372E-01, 4.902E-02,-9.160E-03, 1.041E-02,
112 4-1.925E-01,-1.790E-02, 2.183E-01, 7.490E-03, 4.140E-03,-1.860E-03,
113 5 8.183E-02, 1.926E-02,-1.072E-01,-1.944E-02,-2.770E-03,-5.200E-04,
114 6-3.884E-02,-1.234E-02, 5.410E-02, 1.879E-02, 3.350E-03, 1.040E-03,
115 1 2.948E+01,-3.902E+01, 1.464E+01,-3.335E+00, 5.054E-01,-5.915E-02,
116 2 2.559E+01,-3.955E+01, 1.661E+01,-4.299E+00, 6.904E-01,-8.243E-02,
117 3-1.663E+00, 1.176E+00, 1.118E+00,-7.099E-01, 1.948E-01,-2.404E-02,
118 4-2.168E-01, 8.170E-01,-7.169E-01, 1.851E-01,-1.924E-02,-3.250E-03,
119 5 2.088E-01,-4.355E-01, 2.239E-01,-2.446E-02,-3.620E-03, 1.910E-03,
120 6-9.097E-02, 1.601E-01,-5.681E-02,-2.500E-03, 2.580E-03,-4.700E-04/
121 DATA (((CEHLQ(IX,IT,NX,4,2),IX=1,6),IT=1,6),NX=1,2)/
122 1 2.367E+00, 4.453E-01, 3.660E-01, 9.467E-02, 1.341E-01, 1.661E-02,
123 2-3.170E+00,-1.795E+00, 3.313E-02,-2.874E-01,-9.827E-02,-7.119E-02,
124 3 1.823E+00, 1.457E+00,-2.465E-01, 3.739E-02, 6.090E-03, 1.814E-02,
125 4-1.033E+00,-9.827E-01, 2.136E-01, 1.169E-01, 5.001E-02, 1.684E-02,
126 5 5.133E-01, 5.259E-01,-1.173E-01,-1.139E-01,-4.988E-02,-2.021E-02,
127 6-2.881E-01,-3.145E-01, 5.667E-02, 9.161E-02, 4.568E-02, 1.951E-02,
128 1 3.036E+01,-4.062E+01, 1.578E+01,-3.699E+00, 6.020E-01,-7.031E-02,
129 2 2.700E+01,-4.167E+01, 1.770E+01,-4.804E+00, 7.862E-01,-1.060E-01,
130 3-1.909E+00, 1.357E+00, 1.127E+00,-7.181E-01, 2.232E-01,-2.481E-02,
131 4-2.488E-01, 9.781E-01,-8.127E-01, 2.094E-01,-2.997E-02,-4.710E-03,
132 5 2.506E-01,-5.427E-01, 2.672E-01,-3.103E-02,-1.800E-03, 2.870E-03,
133 6-1.128E-01, 2.087E-01,-6.972E-02,-2.480E-03, 2.630E-03,-8.400E-04/
134C...Expansion coefficients for strange sea quark distribution.
135 DATA (((CEHLQ(IX,IT,NX,5,1),IX=1,6),IT=1,6),NX=1,2)/
136 1 4.968E-02,-4.173E-02, 2.102E-02,-3.270E-03, 3.240E-03,-6.700E-04,
137 2-6.150E-03,-1.294E-02, 6.740E-03,-6.890E-03, 9.000E-04,-1.510E-03,
138 3-8.580E-03, 5.050E-03,-4.900E-03,-1.600E-04,-9.400E-04,-1.500E-04,
139 4 7.840E-03, 1.510E-03, 2.220E-03, 1.400E-03, 7.000E-04, 3.500E-04,
140 5-4.410E-03,-2.220E-03,-8.900E-04,-8.500E-04,-3.600E-04,-2.000E-04,
141 6 2.520E-03, 1.840E-03, 4.100E-04, 3.900E-04, 1.600E-04, 9.000E-05,
142 1 9.235E-01,-1.085E+00, 3.464E-01,-7.210E-02, 9.140E-03,-9.100E-04,
143 2 9.315E-01,-1.274E+00, 4.512E-01,-9.775E-02, 1.380E-02,-1.310E-03,
144 3 4.739E-02,-1.296E-01, 8.482E-02,-2.642E-02, 4.760E-03,-5.700E-04,
145 4-2.653E-02, 4.953E-02,-1.735E-02, 1.750E-03, 2.800E-04,-6.000E-05,
146 5 6.940E-03,-1.132E-02, 1.480E-03, 6.500E-04,-2.100E-04, 0.000E+00,
147 6-1.680E-03, 2.340E-03, 4.200E-04,-3.400E-04, 5.000E-05, 1.000E-05/
148 DATA (((CEHLQ(IX,IT,NX,5,2),IX=1,6),IT=1,6),NX=1,2)/
149 1 6.478E-02,-4.537E-02, 1.643E-02,-3.490E-03, 2.710E-03,-6.700E-04,
150 2-2.223E-02,-2.126E-02, 1.247E-02,-6.290E-03, 1.120E-03,-1.440E-03,
151 3-1.340E-03, 1.362E-02,-6.130E-03,-7.900E-04,-9.000E-04,-2.000E-04,
152 4 5.080E-03,-3.610E-03, 1.700E-03, 1.830E-03, 6.800E-04, 4.000E-04,
153 5-3.580E-03, 6.000E-05,-2.600E-04,-1.050E-03,-3.800E-04,-2.300E-04,
154 6 2.420E-03, 9.300E-04,-1.000E-04, 4.500E-04, 1.700E-04, 1.100E-04,
155 1 9.868E-01,-1.171E+00, 3.940E-01,-8.459E-02, 1.124E-02,-1.250E-03,
156 2 1.001E+00,-1.383E+00, 5.044E-01,-1.152E-01, 1.658E-02,-1.830E-03,
157 3 4.928E-02,-1.368E-01, 9.021E-02,-2.935E-02, 5.800E-03,-6.600E-04,
158 4-3.133E-02, 5.785E-02,-2.023E-02, 2.630E-03, 1.600E-04,-8.000E-05,
159 5 8.840E-03,-1.416E-02, 1.900E-03, 5.800E-04,-2.500E-04, 1.000E-05,
160 6-2.300E-03, 3.080E-03, 5.500E-04,-3.700E-04, 7.000E-05, 1.000E-05/
161C...Expansion coefficients for charm sea quark distribution.
162 DATA (((CEHLQ(IX,IT,NX,6,1),IX=1,6),IT=1,6),NX=1,2)/
163 1 9.270E-03,-1.817E-02, 9.590E-03,-6.390E-03, 1.690E-03,-1.540E-03,
164 2 5.710E-03,-1.188E-02, 6.090E-03,-4.650E-03, 1.240E-03,-1.310E-03,
165 3-3.960E-03, 7.100E-03,-3.590E-03, 1.840E-03,-3.900E-04, 3.400E-04,
166 4 1.120E-03,-1.960E-03, 1.120E-03,-4.800E-04, 1.000E-04,-4.000E-05,
167 5 4.000E-05,-3.000E-05,-1.800E-04, 9.000E-05,-5.000E-05,-2.000E-05,
168 6-4.200E-04, 7.300E-04,-1.600E-04, 5.000E-05, 5.000E-05, 5.000E-05,
169 1 8.098E-01,-1.042E+00, 3.398E-01,-6.824E-02, 8.760E-03,-9.000E-04,
170 2 8.961E-01,-1.217E+00, 4.339E-01,-9.287E-02, 1.304E-02,-1.290E-03,
171 3 3.058E-02,-1.040E-01, 7.604E-02,-2.415E-02, 4.600E-03,-5.000E-04,
172 4-2.451E-02, 4.432E-02,-1.651E-02, 1.430E-03, 1.200E-04,-1.000E-04,
173 5 1.122E-02,-1.457E-02, 2.680E-03, 5.800E-04,-1.200E-04, 3.000E-05,
174 6-7.730E-03, 7.330E-03,-7.600E-04,-2.400E-04, 1.000E-05, 0.000E+00/
175 DATA (((CEHLQ(IX,IT,NX,6,2),IX=1,6),IT=1,6),NX=1,2)/
176 1 9.980E-03,-1.945E-02, 1.055E-02,-6.870E-03, 1.860E-03,-1.560E-03,
177 2 5.700E-03,-1.203E-02, 6.250E-03,-4.860E-03, 1.310E-03,-1.370E-03,
178 3-4.490E-03, 7.990E-03,-4.170E-03, 2.050E-03,-4.400E-04, 3.300E-04,
179 4 1.470E-03,-2.480E-03, 1.460E-03,-5.700E-04, 1.200E-04,-1.000E-05,
180 5-9.000E-05, 1.500E-04,-3.200E-04, 1.200E-04,-6.000E-05,-4.000E-05,
181 6-4.200E-04, 7.600E-04,-1.400E-04, 4.000E-05, 7.000E-05, 5.000E-05,
182 1 8.698E-01,-1.131E+00, 3.836E-01,-8.111E-02, 1.048E-02,-1.300E-03,
183 2 9.626E-01,-1.321E+00, 4.854E-01,-1.091E-01, 1.583E-02,-1.700E-03,
184 3 3.057E-02,-1.088E-01, 8.022E-02,-2.676E-02, 5.590E-03,-5.600E-04,
185 4-2.845E-02, 5.164E-02,-1.918E-02, 2.210E-03,-4.000E-05,-1.500E-04,
186 5 1.311E-02,-1.751E-02, 3.310E-03, 5.100E-04,-1.200E-04, 5.000E-05,
187 6-8.590E-03, 8.380E-03,-9.200E-04,-2.600E-04, 1.000E-05,-1.000E-05/
188C...Expansion coefficients for bottom sea quark distribution.
189 DATA (((CEHLQ(IX,IT,NX,7,1),IX=1,6),IT=1,6),NX=1,2)/
190 1 9.010E-03,-1.401E-02, 7.150E-03,-4.130E-03, 1.260E-03,-1.040E-03,
191 2 6.280E-03,-9.320E-03, 4.780E-03,-2.890E-03, 9.100E-04,-8.200E-04,
192 3-2.930E-03, 4.090E-03,-1.890E-03, 7.600E-04,-2.300E-04, 1.400E-04,
193 4 3.900E-04,-1.200E-03, 4.400E-04,-2.500E-04, 2.000E-05,-2.000E-05,
194 5 2.600E-04, 1.400E-04,-8.000E-05, 1.000E-04, 1.000E-05, 1.000E-05,
195 6-2.600E-04, 3.200E-04, 1.000E-05,-1.000E-05, 1.000E-05,-1.000E-05,
196 1 8.029E-01,-1.075E+00, 3.792E-01,-7.843E-02, 1.007E-02,-1.090E-03,
197 2 7.903E-01,-1.099E+00, 4.153E-01,-9.301E-02, 1.317E-02,-1.410E-03,
198 3-1.704E-02,-1.130E-02, 2.882E-02,-1.341E-02, 3.040E-03,-3.600E-04,
199 4-7.200E-04, 7.230E-03,-5.160E-03, 1.080E-03,-5.000E-05,-4.000E-05,
200 5 3.050E-03,-4.610E-03, 1.660E-03,-1.300E-04,-1.000E-05, 1.000E-05,
201 6-4.360E-03, 5.230E-03,-1.610E-03, 2.000E-04,-2.000E-05, 0.000E+00/
202 DATA (((CEHLQ(IX,IT,NX,7,2),IX=1,6),IT=1,6),NX=1,2)/
203 1 8.980E-03,-1.459E-02, 7.510E-03,-4.410E-03, 1.310E-03,-1.070E-03,
204 2 5.970E-03,-9.440E-03, 4.800E-03,-3.020E-03, 9.100E-04,-8.500E-04,
205 3-3.050E-03, 4.440E-03,-2.100E-03, 8.500E-04,-2.400E-04, 1.400E-04,
206 4 5.300E-04,-1.300E-03, 5.600E-04,-2.700E-04, 3.000E-05,-2.000E-05,
207 5 2.000E-04, 1.400E-04,-1.100E-04, 1.000E-04, 0.000E+00, 0.000E+00,
208 6-2.600E-04, 3.200E-04, 0.000E+00,-3.000E-05, 1.000E-05,-1.000E-05,
209 1 8.672E-01,-1.174E+00, 4.265E-01,-9.252E-02, 1.244E-02,-1.460E-03,
210 2 8.500E-01,-1.194E+00, 4.630E-01,-1.083E-01, 1.614E-02,-1.830E-03,
211 3-2.241E-02,-5.630E-03, 2.815E-02,-1.425E-02, 3.520E-03,-4.300E-04,
212 4-7.300E-04, 8.030E-03,-5.780E-03, 1.380E-03,-1.300E-04,-4.000E-05,
213 5 3.460E-03,-5.380E-03, 1.960E-03,-2.100E-04, 1.000E-05, 1.000E-05,
214 6-4.850E-03, 5.950E-03,-1.890E-03, 2.600E-04,-3.000E-05, 0.000E+00/
215C...Expansion coefficients for top sea quark distribution.
216 DATA (((CEHLQ(IX,IT,NX,8,1),IX=1,6),IT=1,6),NX=1,2)/
217 1 4.410E-03,-7.480E-03, 3.770E-03,-2.580E-03, 7.300E-04,-7.100E-04,
218 2 3.840E-03,-6.050E-03, 3.030E-03,-2.030E-03, 5.800E-04,-5.900E-04,
219 3-8.800E-04, 1.660E-03,-7.500E-04, 4.700E-04,-1.000E-04, 1.000E-04,
220 4-8.000E-05,-1.500E-04, 1.200E-04,-9.000E-05, 3.000E-05, 0.000E+00,
221 5 1.300E-04,-2.200E-04,-2.000E-05,-2.000E-05,-2.000E-05,-2.000E-05,
222 6-7.000E-05, 1.900E-04,-4.000E-05, 2.000E-05, 0.000E+00, 0.000E+00,
223 1 6.623E-01,-9.248E-01, 3.519E-01,-7.930E-02, 1.110E-02,-1.180E-03,
224 2 6.380E-01,-9.062E-01, 3.582E-01,-8.479E-02, 1.265E-02,-1.390E-03,
225 3-2.581E-02, 2.125E-02, 4.190E-03,-4.980E-03, 1.490E-03,-2.100E-04,
226 4 7.100E-04, 5.300E-04,-1.270E-03, 3.900E-04,-5.000E-05,-1.000E-05,
227 5 3.850E-03,-5.060E-03, 1.860E-03,-3.500E-04, 4.000E-05, 0.000E+00,
228 6-3.530E-03, 4.460E-03,-1.500E-03, 2.700E-04,-3.000E-05, 0.000E+00/
229 DATA (((CEHLQ(IX,IT,NX,8,2),IX=1,6),IT=1,6),NX=1,2)/
230 1 4.260E-03,-7.530E-03, 3.830E-03,-2.680E-03, 7.600E-04,-7.300E-04,
231 2 3.640E-03,-6.050E-03, 3.030E-03,-2.090E-03, 5.900E-04,-6.000E-04,
232 3-9.200E-04, 1.710E-03,-8.200E-04, 5.000E-04,-1.200E-04, 1.000E-04,
233 4-5.000E-05,-1.600E-04, 1.300E-04,-9.000E-05, 3.000E-05, 0.000E+00,
234 5 1.300E-04,-2.100E-04,-1.000E-05,-2.000E-05,-2.000E-05,-1.000E-05,
235 6-8.000E-05, 1.800E-04,-5.000E-05, 2.000E-05, 0.000E+00, 0.000E+00,
236 1 7.146E-01,-1.007E+00, 3.932E-01,-9.246E-02, 1.366E-02,-1.540E-03,
237 2 6.856E-01,-9.828E-01, 3.977E-01,-9.795E-02, 1.540E-02,-1.790E-03,
238 3-3.053E-02, 2.758E-02, 2.150E-03,-4.880E-03, 1.640E-03,-2.500E-04,
239 4 9.200E-04, 4.200E-04,-1.340E-03, 4.600E-04,-8.000E-05,-1.000E-05,
240 5 4.230E-03,-5.660E-03, 2.140E-03,-4.300E-04, 6.000E-05, 0.000E+00,
241 6-3.890E-03, 5.000E-03,-1.740E-03, 3.300E-04,-4.000E-05, 0.000E+00/
242
243C...The following data lines are coefficients needed in the
244C...Duke, Owens proton structure function parametrizations, see below.
245C...Expansion coefficients for (up+down) valence quark distribution.
246 DATA ((CDO(IP,IS,1,1),IS=1,6),IP=1,3)/
247 1 4.190E-01, 3.460E+00, 4.400E+00, 0.000E+00, 0.000E+00, 0.000E+00,
248 2 4.000E-03, 7.240E-01,-4.860E+00, 0.000E+00, 0.000E+00, 0.000E+00,
249 3-7.000E-03,-6.600E-02, 1.330E+00, 0.000E+00, 0.000E+00, 0.000E+00/
250 DATA ((CDO(IP,IS,1,2),IS=1,6),IP=1,3)/
251 1 3.740E-01, 3.330E+00, 6.030E+00, 0.000E+00, 0.000E+00, 0.000E+00,
252 2 1.400E-02, 7.530E-01,-6.220E+00, 0.000E+00, 0.000E+00, 0.000E+00,
253 3 0.000E+00,-7.600E-02, 1.560E+00, 0.000E+00, 0.000E+00, 0.000E+00/
254C...Expansion coefficients for down valence quark distribution.
255 DATA ((CDO(IP,IS,2,1),IS=1,6),IP=1,3)/
256 1 7.630E-01, 4.000E+00, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00,
257 2-2.370E-01, 6.270E-01,-4.210E-01, 0.000E+00, 0.000E+00, 0.000E+00,
258 3 2.600E-02,-1.900E-02, 3.300E-02, 0.000E+00, 0.000E+00, 0.000E+00/
259 DATA ((CDO(IP,IS,2,2),IS=1,6),IP=1,3)/
260 1 7.610E-01, 3.830E+00, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00,
261 2-2.320E-01, 6.270E-01,-4.180E-01, 0.000E+00, 0.000E+00, 0.000E+00,
262 3 2.300E-02,-1.900E-02, 3.600E-02, 0.000E+00, 0.000E+00, 0.000E+00/
263C...Expansion coefficients for (up+down+strange) sea quark distribution.
264 DATA ((CDO(IP,IS,3,1),IS=1,6),IP=1,3)/
265 1 1.265E+00, 0.000E+00, 8.050E+00, 0.000E+00, 0.000E+00, 0.000E+00,
266 2-1.132E+00,-3.720E-01, 1.590E+00, 6.310E+00,-1.050E+01, 1.470E+01,
267 3 2.930E-01,-2.900E-02,-1.530E-01,-2.730E-01,-3.170E+00, 9.800E+00/
268 DATA ((CDO(IP,IS,3,2),IS=1,6),IP=1,3)/
269 1 1.670E+00, 0.000E+00, 9.150E+00, 0.000E+00, 0.000E+00, 0.000E+00,
270 2-1.920E+00,-2.730E-01, 5.300E-01, 1.570E+01,-1.010E+02, 2.230E+02,
271 3 5.820E-01,-1.640E-01,-7.630E-01,-2.830E+00, 4.470E+01,-1.170E+02/
272C...Expansion coefficients for charm sea quark distribution.
273 DATA ((CDO(IP,IS,4,1),IS=1,6),IP=1,3)/
274 1 0.000E+00,-3.600E-02, 6.350E+00, 0.000E+00, 0.000E+00, 0.000E+00,
275 2 1.350E-01,-2.220E-01, 3.260E+00,-3.030E+00, 1.740E+01,-1.790E+01,
276 3-7.500E-02,-5.800E-02,-9.090E-01, 1.500E+00,-1.130E+01, 1.560E+01/
277 DATA ((CDO(IP,IS,4,2),IS=1,6),IP=1,3)/
278 1 0.000E+00,-1.200E-01, 3.510E+00, 0.000E+00, 0.000E+00, 0.000E+00,
279 2 6.700E-02,-2.330E-01, 3.660E+00,-4.740E-01, 9.500E+00,-1.660E+01,
280 3-3.100E-02,-2.300E-02,-4.530E-01, 3.580E-01,-5.430E+00, 1.550E+01/
281C...Expansion coefficients for gluon distribution.
282 DATA ((CDO(IP,IS,5,1),IS=1,6),IP=1,3)/
283 1 1.560E+00, 0.000E+00, 6.000E+00, 9.000E+00, 0.000E+00, 0.000E+00,
284 2-1.710E+00,-9.490E-01, 1.440E+00,-7.190E+00,-1.650E+01, 1.530E+01,
285 3 6.380E-01, 3.250E-01,-1.050E+00, 2.550E-01, 1.090E+01,-1.010E+01/
286 DATA ((CDO(IP,IS,5,2),IS=1,6),IP=1,3)/
287 1 8.790E-01, 0.000E+00, 4.000E+00, 9.000E+00, 0.000E+00, 0.000E+00,
288 2-9.710E-01,-1.160E+00, 1.230E+00,-5.640E+00,-7.540E+00,-5.960E-01,
289 3 4.340E-01, 4.760E-01,-2.540E-01,-8.170E-01, 5.500E+00, 1.260E-01/
290
291C...The following data lines are coefficients needed in the
292C...Owens pion structure function parametrizations, see below.
293C...Expansion coefficients for up and down valence quark distributions.
294 DATA ((COW(IP,IS,1,1),IS=1,5),IP=1,3)/
295 1 4.0000E-01, 7.0000E-01, 0.0000E+00, 0.0000E+00, 0.0000E+00,
296 2 -6.2120E-02, 6.4780E-01, 0.0000E+00, 0.0000E+00, 0.0000E+00,
297 3 -7.1090E-03, 1.3350E-02, 0.0000E+00, 0.0000E+00, 0.0000E+00/
298 DATA ((COW(IP,IS,1,2),IS=1,5),IP=1,3)/
299 1 4.0000E-01, 6.2800E-01, 0.0000E+00, 0.0000E+00, 0.0000E+00,
300 2 -5.9090E-02, 6.4360E-01, 0.0000E+00, 0.0000E+00, 0.0000E+00,
301 3 -6.5240E-03, 1.4510E-02, 0.0000E+00, 0.0000E+00, 0.0000E+00/
302C...Expansion coefficients for gluon distribution.
303 DATA ((COW(IP,IS,2,1),IS=1,5),IP=1,3)/
304 1 8.8800E-01, 0.0000E+00, 3.1100E+00, 6.0000E+00, 0.0000E+00,
305 2 -1.8020E+00, -1.5760E+00, -1.3170E-01, 2.8010E+00, -1.7280E+01,
306 3 1.8120E+00, 1.2000E+00, 5.0680E-01, -1.2160E+01, 2.0490E+01/
307 DATA ((COW(IP,IS,2,2),IS=1,5),IP=1,3)/
308 1 7.9400E-01, 0.0000E+00, 2.8900E+00, 6.0000E+00, 0.0000E+00,
309 2 -9.1440E-01, -1.2370E+00, 5.9660E-01, -3.6710E+00, -8.1910E+00,
310 3 5.9660E-01, 6.5820E-01, -2.5500E-01, -2.3040E+00, 7.7580E+00/
311C...Expansion coefficients for (up+down+strange) quark sea distribution.
312 DATA ((COW(IP,IS,3,1),IS=1,5),IP=1,3)/
313 1 9.0000E-01, 0.0000E+00, 5.0000E+00, 0.0000E+00, 0.0000E+00,
314 2 -2.4280E-01, -2.1200E-01, 8.6730E-01, 1.2660E+00, 2.3820E+00,
315 3 1.3860E-01, 3.6710E-03, 4.7470E-02, -2.2150E+00, 3.4820E-01/
316 DATA ((COW(IP,IS,3,2),IS=1,5),IP=1,3)/
317 1 9.0000E-01, 0.0000E+00, 5.0000E+00, 0.0000E+00, 0.0000E+00,
318 2 -1.4170E-01, -1.6970E-01, -2.4740E+00, -2.5340E+00, 5.6210E-01,
319 3 -1.7400E-01, -9.6230E-02, 1.5750E+00, 1.3780E+00, -2.7010E-01/
320C...Expansion coefficients for charm quark sea distribution.
321 DATA ((COW(IP,IS,4,1),IS=1,5),IP=1,3)/
322 1 0.0000E+00, -2.2120E-02, 2.8940E+00, 0.0000E+00, 0.0000E+00,
323 2 7.9280E-02, -3.7850E-01, 9.4330E+00, 5.2480E+00, 8.3880E+00,
324 3 -6.1340E-02, -1.0880E-01, -1.0852E+01, -7.1870E+00, -1.1610E+01/
325 DATA ((COW(IP,IS,4,2),IS=1,5),IP=1,3)/
326 1 0.0000E+00, -8.8200E-02, 1.9240E+00, 0.0000E+00, 0.0000E+00,
327 2 6.2290E-02, -2.8920E-01, 2.4240E-01, -4.4630E+00, -8.3670E-01,
328 3 -4.0990E-02, -1.0820E-01, 2.0360E+00, 5.2090E+00, -4.8400E-02/
329
330C...Euler's beta function, requires ordinary Gamma function
331 EULBET(X,Y)=PYGAMM_HIJING(X)*PYGAMM_HIJING(Y)/PYGAMM_HIJING(X+Y)
332
333C...Reset structure functions, check x and hadron flavour.
334 ALAM=0.
335 DO 100 KFL=-6,6
336 100 XPQ(KFL)=0.
337 IF(X.LT.0..OR.X.GT.1.) THEN
338 WRITE(MSTU(11),1000) X
339 RETURN
340 ENDIF
341 KFA=IABS(KF)
342 IF(KFA.NE.211.AND.KFA.NE.2212.AND.KFA.NE.2112) THEN
343 WRITE(MSTU(11),1100) KF
344 RETURN
345 ENDIF
346
347C...Call user-supplied structure function. Select proton/neutron/pion.
348 IF(MSTP(51).EQ.0.OR.MSTP(52).GE.2) THEN
349 KFE=KFA
350 IF(KFA.EQ.2112) KFE=2212
351 CALL PYSTFE_HIJING(KFE,X,Q2,XPQ)
352 GOTO 230
353 ENDIF
354 IF(KFA.EQ.211) GOTO 200
355
356 IF(MSTP(51).EQ.1.OR.MSTP(51).EQ.2) THEN
357C...Proton structure functions from Eichten, Hinchliffe, Lane, Quigg.
358C...Allowed variable range: 5 GeV2 < Q2 < 1E8 GeV2; 1E-4 < x < 1
359
360C...Determine set, Lamdba and x and t expansion variables.
361 NSET=MSTP(51)
362 IF(NSET.EQ.1) ALAM=0.2
363 IF(NSET.EQ.2) ALAM=0.29
364 TMIN=LOG(5./ALAM**2)
365 TMAX=LOG(1E8/ALAM**2)
366 IF(MSTP(52).EQ.0) THEN
367 T=TMIN
368 ELSE
369 T=LOG(Q2/ALAM**2)
370 ENDIF
371 VT=MAX(-1.,MIN(1.,(2.*T-TMAX-TMIN)/(TMAX-TMIN)))
372 NX=1
373 IF(X.LE.0.1) NX=2
374 IF(NX.EQ.1) VX=(2.*X-1.1)/0.9
375 IF(NX.EQ.2) VX=MAX(-1.,(2.*LOG(X)+11.51293)/6.90776)
376 CXS=1.
377 IF(X.LT.1E-4.AND.ABS(PARP(51)-1.).GT.0.01) CXS=
378 & (1E-4/X)**(PARP(51)-1.)
379
380C...Chebyshev polynomials for x and t expansion.
381 TX(1)=1.
382 TX(2)=VX
383 TX(3)=2.*VX**2-1.
384 TX(4)=4.*VX**3-3.*VX
385 TX(5)=8.*VX**4-8.*VX**2+1.
386 TX(6)=16.*VX**5-20.*VX**3+5.*VX
387 TT(1)=1.
388 TT(2)=VT
389 TT(3)=2.*VT**2-1.
390 TT(4)=4.*VT**3-3.*VT
391 TT(5)=8.*VT**4-8.*VT**2+1.
392 TT(6)=16.*VT**5-20.*VT**3+5.*VT
393
394C...Calculate structure functions.
395 DO 120 KFL=1,6
396 XQSUM=0.
397 DO 110 IT=1,6
398 DO 110 IX=1,6
399 110 XQSUM=XQSUM+CEHLQ(IX,IT,NX,KFL,NSET)*TX(IX)*TT(IT)
400 120 XQ(KFL)=XQSUM*(1.-X)**NEHLQ(KFL,NSET)*CXS
401
402C...Put into output array.
403 XPQ(0)=XQ(4)
404 XPQ(1)=XQ(2)+XQ(3)
405 XPQ(2)=XQ(1)+XQ(3)
406 XPQ(3)=XQ(5)
407 XPQ(4)=XQ(6)
408 XPQ(-1)=XQ(3)
409 XPQ(-2)=XQ(3)
410 XPQ(-3)=XQ(5)
411 XPQ(-4)=XQ(6)
412
413C...Special expansion for bottom (threshold effects).
414 IF(MSTP(54).GE.5) THEN
415 IF(NSET.EQ.1) TMIN=8.1905
416 IF(NSET.EQ.2) TMIN=7.4474
417 IF(T.LE.TMIN) GOTO 140
418 VT=MAX(-1.,MIN(1.,(2.*T-TMAX-TMIN)/(TMAX-TMIN)))
419 TT(1)=1.
420 TT(2)=VT
421 TT(3)=2.*VT**2-1.
422 TT(4)=4.*VT**3-3.*VT
423 TT(5)=8.*VT**4-8.*VT**2+1.
424 TT(6)=16.*VT**5-20.*VT**3+5.*VT
425 XQSUM=0.
426 DO 130 IT=1,6
427 DO 130 IX=1,6
428 130 XQSUM=XQSUM+CEHLQ(IX,IT,NX,7,NSET)*TX(IX)*TT(IT)
429 XPQ(5)=XQSUM*(1.-X)**NEHLQ(7,NSET)
430 XPQ(-5)=XPQ(5)
431 140 CONTINUE
432 ENDIF
433
434C...Special expansion for top (threshold effects).
435 IF(MSTP(54).GE.6) THEN
436 IF(NSET.EQ.1) TMIN=11.5528
437 IF(NSET.EQ.2) TMIN=10.8097
438 TMIN=TMIN+2.*LOG(PMAS(6,1)/30.)
439 TMAX=TMAX+2.*LOG(PMAS(6,1)/30.)
440 IF(T.LE.TMIN) GOTO 160
441 VT=MAX(-1.,MIN(1.,(2.*T-TMAX-TMIN)/(TMAX-TMIN)))
442 TT(1)=1.
443 TT(2)=VT
444 TT(3)=2.*VT**2-1.
445 TT(4)=4.*VT**3-3.*VT
446 TT(5)=8.*VT**4-8.*VT**2+1.
447 TT(6)=16.*VT**5-20.*VT**3+5.*VT
448 XQSUM=0.
449 DO 150 IT=1,6
450 DO 150 IX=1,6
451 150 XQSUM=XQSUM+CEHLQ(IX,IT,NX,8,NSET)*TX(IX)*TT(IT)
452 XPQ(6)=XQSUM*(1.-X)**NEHLQ(8,NSET)
453 XPQ(-6)=XPQ(6)
454 160 CONTINUE
455 ENDIF
456
457 ELSEIF(MSTP(51).EQ.3.OR.MSTP(51).EQ.4) THEN
458C...Proton structure functions from Duke, Owens.
459C...Allowed variable range: 4 GeV2 < Q2 < approx 1E6 GeV2.
460
461C...Determine set, Lambda and s expansion parameter.
462 NSET=MSTP(51)-2
463 IF(NSET.EQ.1) ALAM=0.2
464 IF(NSET.EQ.2) ALAM=0.4
465 IF(MSTP(52).LE.0) THEN
466 SD=0.
467 ELSE
468 SD=LOG(LOG(MAX(Q2,4.)/ALAM**2)/LOG(4./ALAM**2))
469 ENDIF
470
471C...Calculate structure functions.
472 DO 180 KFL=1,5
473 DO 170 IS=1,6
474 170 TS(IS)=CDO(1,IS,KFL,NSET)+CDO(2,IS,KFL,NSET)*SD+
475 & CDO(3,IS,KFL,NSET)*SD**2
476 IF(KFL.LE.2) THEN
477 XQ(KFL)=X**TS(1)*(1.-X)**TS(2)*(1.+TS(3)*X)/(EULBET(TS(1),
478 & TS(2)+1.)*(1.+TS(3)*TS(1)/(TS(1)+TS(2)+1.)))
479 ELSE
480 XQ(KFL)=TS(1)*X**TS(2)*(1.-X)**TS(3)*(1.+TS(4)*X+TS(5)*X**2+
481 & TS(6)*X**3)
482 ENDIF
483 180 CONTINUE
484
485C...Put into output arrays.
486 XPQ(0)=XQ(5)
487 XPQ(1)=XQ(2)+XQ(3)/6.
488 XPQ(2)=3.*XQ(1)-XQ(2)+XQ(3)/6.
489 XPQ(3)=XQ(3)/6.
490 XPQ(4)=XQ(4)
491 XPQ(-1)=XQ(3)/6.
492 XPQ(-2)=XQ(3)/6.
493 XPQ(-3)=XQ(3)/6.
494 XPQ(-4)=XQ(4)
495
496C...Proton structure functions from Diemoz, Ferroni, Longo, Martinelli.
497C...These are accessed via PYSTFE_HIJING since the files needed may not always
498C...available.
499 ELSEIF(MSTP(51).GE.11.AND.MSTP(51).LE.13) THEN
500 CALL PYSTFE_HIJING(2212,X,Q2,XPQ)
501
502C...Unknown proton parametrization.
503 ELSE
504 WRITE(MSTU(11),1200) MSTP(51)
505 ENDIF
506 GOTO 230
507
508 200 IF((MSTP(51).GE.1.AND.MSTP(51).LE.4).OR.
509 &(MSTP(51).GE.11.AND.MSTP(51).LE.13)) THEN
510C...Pion structure functions from Owens.
511C...Allowed variable range: 4 GeV2 < Q2 < approx 2000 GeV2.
512
513C...Determine set, Lambda and s expansion variable.
514 NSET=1
515 IF(MSTP(51).EQ.2.OR.MSTP(51).EQ.4.OR.MSTP(51).EQ.13) NSET=2
516 IF(NSET.EQ.1) ALAM=0.2
517 IF(NSET.EQ.2) ALAM=0.4
518 IF(MSTP(52).LE.0) THEN
519 SD=0.
520 ELSE
521 SD=LOG(LOG(MAX(Q2,4.)/ALAM**2)/LOG(4./ALAM**2))
522 ENDIF
523
524C...Calculate structure functions.
525 DO 220 KFL=1,4
526 DO 210 IS=1,5
527 210 TS(IS)=COW(1,IS,KFL,NSET)+COW(2,IS,KFL,NSET)*SD+
528 & COW(3,IS,KFL,NSET)*SD**2
529 IF(KFL.EQ.1) THEN
530 XQ(KFL)=X**TS(1)*(1.-X)**TS(2)/EULBET(TS(1),TS(2)+1.)
531 ELSE
532 XQ(KFL)=TS(1)*X**TS(2)*(1.-X)**TS(3)*(1.+TS(4)*X+TS(5)*X**2)
533 ENDIF
534 220 CONTINUE
535
536C...Put into output arrays.
537 XPQ(0)=XQ(2)
538 XPQ(1)=XQ(3)/6.
539 XPQ(2)=XQ(1)+XQ(3)/6.
540 XPQ(3)=XQ(3)/6.
541 XPQ(4)=XQ(4)
542 XPQ(-1)=XQ(1)+XQ(3)/6.
543 XPQ(-2)=XQ(3)/6.
544 XPQ(-3)=XQ(3)/6.
545 XPQ(-4)=XQ(4)
546
547C...Unknown pion parametrization.
548 ELSE
549 WRITE(MSTU(11),1200) MSTP(51)
550 ENDIF
551
552C...Isospin conjugation for neutron, charge conjugation for antipart.
553 230 IF(KFA.EQ.2112) THEN
554 XPS=XPQ(1)
555 XPQ(1)=XPQ(2)
556 XPQ(2)=XPS
557 XPS=XPQ(-1)
558 XPQ(-1)=XPQ(-2)
559 XPQ(-2)=XPS
560 ENDIF
561 IF(KF.LT.0) THEN
562 DO 240 KFL=1,4
563 XPS=XPQ(KFL)
564 XPQ(KFL)=XPQ(-KFL)
565 240 XPQ(-KFL)=XPS
566 ENDIF
567
568C...Check positivity and reset above maximum allowed flavour.
569 DO 250 KFL=-6,6
570 XPQ(KFL)=MAX(0.,XPQ(KFL))
571 250 IF(IABS(KFL).GT.MSTP(54)) XPQ(KFL)=0.
572
573C...consider nuclear effect on the structure function
574 IF((JBT.NE.1.AND.JBT.NE.2).OR.IHPR2(6).EQ.0
575 & .OR.IHNT2(16).EQ.1) GO TO 400
576 ATNM=IHNT2(2*JBT-1)
577 IF(ATNM.LE.1.0) GO TO 400
578 IF(JBT.EQ.1) THEN
579 BBR2=(YP(1,IHNT2(11))**2+YP(2,IHNT2(11))**2)/1.44/ATNM**0.66666
580 ELSEIF(JBT.EQ.2) THEN
581 BBR2=(YT(1,IHNT2(12))**2+YT(2,IHNT2(12))**2)/1.44/ATNM**0.66666
582 ENDIF
583 BBR2=MIN(1.0,BBR2)
584 ABX=(ATNM**0.33333333-1.0)
585 APX=HIPR1(6)*4.0/3.0*ABX*SQRT(1.0-BBR2)
586 AAX=1.192*ALOG(ATNM)**0.1666666
587 RRX=AAX*(X**3-1.2*X**2+0.21*X)+1.0
588 & -(APX-1.079*ABX*SQRT(X)/ALOG(ATNM+1.0))*EXP(-X**2.0/0.01)
589 DO 300 KFL=-6,6
590 XPQ(KFL)=XPQ(KFL)*RRX
591300 CONTINUE
592C ********consider the nuclear effect on the structure
593C fucntion which also depends on the impact
594C parameter of the nuclear reaction
595
596400 CONTINUE
597C...Formats for error printouts.
598 1000 FORMAT(' Error: x value outside physical range, x =',1P,E12.3)
599 1100 FORMAT(' Error: illegal particle code for structure function,',
600 &' KF =',I5)
601 1200 FORMAT(' Error: bad value of parameter MSTP(51) in PYSTFU_HIJING,'
602 $ ,' MSTP(51) =',I5)
603
604 RETURN
605 END