HIJING/hipyset1_35/pystfu_hijing.F

   1 * $Id$
   2
   3 C*********************************************************************
   4
   5       SUBROUTINE PYSTFU_HIJING(KF,X,Q2,XPQ,JBT)
   6
   7 C                       *******JBT specifies beam or target of the particle
   8 C...Gives proton and pi+ parton structure functions according to a few
   9 C...different parametrizations. Note that what is coded is x times the
  10 C...probability distribution, i.e. xq(x,Q2) etc.
  11 #include "hiparnt.inc"
  12 #include "hijcrdn.inc"
  13 C                       ********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
  21 C...The following data lines are coefficients needed in the
  22 C...Eichten, Hinchliffe, Lane, Quigg proton structure function
  23 C...parametrizations, see below.
  24 C...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/
  26 C...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/
  53 C...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/
  80 C...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/
 107 C...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/
 134 C...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/
 161 C...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/
 188 C...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/
 215 C...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
 243 C...The following data lines are coefficients needed in the
 244 C...Duke, Owens proton structure function parametrizations, see below.
 245 C...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/
 254 C...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/
 263 C...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/
 272 C...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/
 281 C...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
 291 C...The following data lines are coefficients needed in the
 292 C...Owens pion structure function parametrizations, see below.
 293 C...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/
 302 C...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/
 311 C...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/
 320 C...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
 330 C...Euler's beta function, requires ordinary Gamma function
 331       EULBET(X,Y)=PYGAMM_HIJING(X)*PYGAMM_HIJING(Y)/PYGAMM_HIJING(X+Y)
 332
 333 C...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
 347 C...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
 357 C...Proton structure functions from Eichten, Hinchliffe, Lane, Quigg.
 358 C...Allowed variable range: 5 GeV2 < Q2 < 1E8 GeV2; 1E-4 < x < 1
 359
 360 C...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
 380 C...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
 394 C...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
 402 C...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
 413 C...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
 434 C...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
 458 C...Proton structure functions from Duke, Owens.
 459 C...Allowed variable range: 4 GeV2 < Q2 < approx 1E6 GeV2.
 460
 461 C...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
 471 C...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
 485 C...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
 496 C...Proton structure functions from Diemoz, Ferroni, Longo, Martinelli.
 497 C...These are accessed via PYSTFE_HIJING since the files needed may not always
 498 C...available.
 499       ELSEIF(MSTP(51).GE.11.AND.MSTP(51).LE.13) THEN
 500         CALL PYSTFE_HIJING(2212,X,Q2,XPQ)
 501
 502 C...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
 510 C...Pion structure functions from Owens.
 511 C...Allowed variable range: 4 GeV2 < Q2 < approx 2000 GeV2.
 512
 513 C...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
 524 C...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
 536 C...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
 547 C...Unknown pion parametrization.
 548       ELSE
 549         WRITE(MSTU(11),1200) MSTP(51)
 550       ENDIF
 551
 552 C...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
 568 C...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
 573 C...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
 591 300     CONTINUE
 592 C                       ********consider the nuclear effect on the structure
 593 C                               fucntion which also depends on the impact
 594 C                               parameter of the nuclear reaction
 595
 596 400     CONTINUE
 597 C...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