C slightly modified: FIO 10 may 99: IRT -> IRT2. DEC compiler problem. C============================================================================ C CTEQ4 Parton Distribution Functions in Parametrized Form C C Preliminary version: Aug. 27, 1998 C C Ref[1]: "IMPROVED PARTON DISTRIBUTIONS FROM GLOBAL ANALYSIS OF RECENT DEEP C INELASTIC SCATTERING AND INCLUSIVE JET DATA" C By: H.L. Lai, J. Huston, S. Kuhlmann, F. Olness, J. Owens, D. Soper C W.K. Tung, H. Weerts C Phys. Rev. D55, 1280 (1997) C C Ref[2]: "CHARM PRODUCTION AND PARTON DISTRIBUTIONS" C By: H.L. Lai and W.K. Tung C Z. Phys. C74, 463 (1997) C These parametrizations were obtained by Brendan Foster. C C The implementation of the parametrization is the same as for CTEQ3 C and CTEQ2 parton distributions. C ------------------------------------------- C This file contains three versions of the same Cteq4 parton distributions: C C Two "front-end" subprograms: C FUNCTION Ctq4Pd (Iset, Iparton, X, Q, Irt) C returns the PROBABILITY density for a GIVEN flavor; C SUBROUTINE Ctq4Pds(Iset, Pdf, XX, QQ, Irt) C returns an array of MOMENTUM densities for ALL flavors; C One lower-level subprogram: C FUNCTION Ctq4df (Iset, Iprtn, XX, QQ, Irt) C returns the MOMENTUM density of a GIVEN valence or sea distribution. C One supplementary function to return the QCD lambda parameter C concerning these distributions is also included (see below). C Although DOUBLE PRECISION is used, conversion to SINGLE PRECISION C is straightforward by removing the C Implicit Double Precision statements. C Since this is an initial distribution of the parametrized version, it C is useful for the authors to know any possible problems. Please send C correspondences to: Tung@pa.msu.edu. C ------------------------------------------- C Detailed instructions follow. C Name convention for CTEQ distributions: CTEQnSx where C n : version number (currently n = 3) C S : factorization scheme label: = [M L D] for [MS-bar LO DIS] c resp. C x : special characteristics, if any C Explanation of functional arguments: C Iset is the set label; C This package contains 9 of the 13 sets of CTEQ4 PDF's. C which were made available previously in table form: C --------------------------------------------------------------------------- C Iset PDF Description Alpha_s(Mz) Lam4 Lam5 C --------------------------------------------------------------------------- C Ref[1] C 1 CTEQ4M Standard MSbar scheme 0.116 298 202 C 2 CTEQ4D Standard DIS scheme 0.116 298 202 C 3 CTEQ4L Leading Order 0.132 236 181 C 4 CTEQ4A1 Alpha_s series 0.110 215 140 C 5 CTEQ4A2 Alpha_s series 0.113 254 169 C XXXX [ 6 not included: CTEQ4A3 == CTEQ4M ] C 7 CTEQ4A4 Alpha_s series 0.119 346 239 C 8 CTEQ4A5 Alpha_s series 0.122 401 282 C XXXX [ 9 CTEQ4HJ not yet parametrized ] C XXXX [10 CTEQ4LQ not yet parametrized ] C --------------------------------------------------------------------------- C Ref[2] C 11 CTEQ4HQ Heavy Quark 0.116 298 202 C ----------------------------------------------------------------------------- C Iprtn is the parton label (6, 5, 4, 3, 2, 1, 0, -1, ......, -6) C for (t, b, c, s, d, u, g, u_bar, ..., t_bar) C *** WARNING: We use the parton label 2 as D-quark, and 1 as U-quark which C might be different with your labels. C X, Q are the usual x, Q; C Irt is a return error code (see individual modules for explanation). C C --------------------------------------------- C Since the QCD Lambda value for the various sets are needed more often than C the other parameters in most applications, a special function C WLAMD4 (ISET, IORDER, NEFF) is provided C which returns the lambda value for Neff = 4,5,6 effective flavors as well as C the order these values pertain to. C ---------------------------------------------- C The range of (x, Q) used in this round of global analysis is, approxi- C mately, 0.01 < x < 0.75 ; and 4 GeV^2 < Q^2 < 400 GeV^2 for fixed target C experiments and 0.0001 < x < 0.1 from HERA data. C The range of (x, Q) used in the reparametrization of the QCD evolved C parton distributions is 10E-6 < x < 1 ; 1.6 GeV < Q < 10 TeV. The C functional form of this parametrization is: C A0 * x^A1 * (1-x)^A2 * (1 + A3 * x^A4) * [log(1+1/x)]^A5 C with the A'coefficients being smooth functions of Q. For heavy quarks, C a threshold factor is applied to A0 which simulates the proper Q-dependence C of the QCD evolution in that region according to the renormalization C scheme defined in Collins-Tung, Nucl. Phys. B278, 934 (1986). C Since this function is positive definite and smooth, it provides sensible C extrapolations of the parton distributions if they are called beyond C the original range in an application. There is no artificial boundaries C or sharp cutoff's. C ******************************************************** FUNCTION CTQ4PD (ISET, IPARTON, X, Q, IRT2) IMPLICIT DOUBLE PRECISION (A-H, O-Z) EXTERNAL CTQ4DF IFL = IPARTON KFL = ABS(IPARTON) MKFL=-KFL IRT=IRT2 IF (IFL.EQ.1 .OR. IFL.EQ.2) THEN VL = CTQ4DF (ISET, IFL, X, Q, IRT) ELSE VL = 0. ENDIF SEA = CTQ4DF (ISET, MKFL, X, Q, IRT) CTQ4PD = (VL + SEA) / X RETURN C ******************************************************** END SUBROUTINE CTQ4PDS (ISET, PDF, X, Q, IRT) IMPLICIT DOUBLE PRECISION (A-H, O-Z) DIMENSION PDF (-6:6) IRT = 0 DO 10 IFL= -6,2 IF (IFL.LE.0) THEN PDF(IFL) = CTQ4DF (ISET,IFL,X,Q,IRT1) PDF(-IFL) = PDF(IFL) ELSE PDF(IFL) = CTQ4DF(ISET,IFL,X,Q,IRT1)+PDF(-IFL) ENDIF IRT = IRT + IRT1 10 CONTINUE RETURN C ********************************************************* END FUNCTION CTQ4DF (ISET, IFL, X, QQ, IRT2) IMPLICIT DOUBLE PRECISION (A-H, O-Z) PARAMETER (ILF=3, IEX=5, NSET=15) DIMENSION > A(Nset, 0:ILF, 0:IEX, -6:2), AF(0:IEX), > IORD(NSET), QMN(NSET), ALAM(NSET), > VLM(4:6,NSET), QMS(4:6) SAVE LSET, LFLV, QSTO SAVE SB, SB2, SB3 SAVE A, AF DATA LSET, LFLV, QSTO / 0, -10, 1.2345 / DATA (QMS(I), I=4,6) / 1.60, 5.00, 180.0 / C ======================== C DATA FOR CTEQ4 SETS C CTEQ4M / CTEQ4A3 DATA > Iord(1), Alam(1), Qmn(1) / 2, .202, 1.60 / > (Vlm(I,1), I=4,6) / .239, .158, .063 / C CTEQ4D DATA > Iord(2), Alam(2), Qmn(2) / 2, .202, 1.60 / > (Vlm(I,2), I=4,6) / .239, .158, .063 / C CTEQ4L DATA > Iord(3), Alam(3), Qmn(3) / 1, .181, 1.60 / > (Vlm(I,3), I=4,6) / .239, .158, .063 / C CTEQ4A1 DATA > Iord(4), Alam(4), Qmn(4) / 2, .140, 1.60 / > (Vlm(I,4), I=4,6) / .239, .158, .063 / C CTEQ4A2 DATA > Iord(5), Alam(5), Qmn(5) / 2, .169, 1.60 / > (Vlm(I,5), I=4,6) / .239, .158, .063 / C CTEQ4A4 DATA > Iord(7), Alam(7), Qmn(7) / 2, .239, 1.60 / > (Vlm(I,7), I=4,6) / .239, .158, .063 / C CTEQ4A5 DATA > Iord(8), Alam(8), Qmn(8) / 2, .282, 1.60 / > (Vlm(I,8), I=4,6) / .239, .158, .063 / C CTEQ4LQ DATA > Iord(10), Alam(10), Qmn(10) / 2, .261, 0.49 / > (Vlm(I,10), I=4,6) / .239, .158, .063 / C CTEQ4HQ DATA > Iord(11), Alam(11), Qmn(11) / 2, .202, 1.60 / > (Vlm(I,11), I=4,6) / .239, .158, .063 / C ======================== C --- Iset = 1, cteq4m --- C A0 Data (A( 1,0,0,I-7),I=1,9) / > 2.150E-01, 1.740E+00, 1.107E+00, -2.763E+00, -2.012E+00, > -1.405E+00, 1.110E-01, 4.112E-01, -2.204E-01/ Data (A( 1,1,0,I-7),I=1,9) / > -8.204E+00, -2.821E+00, -3.377E+00, -2.287E+00, 3.192E-01, > -5.055E-01, -2.356E+00, 5.052E-02, -2.326E-01/ Data (A( 1,2,0,I-7),I=1,9) / > 7.449E+00, -6.476E+00, -4.713E+00, 1.452E+00, -3.564E+00, > 1.620E+00, 7.451E-01, -1.631E-01, -3.938E-01/ Data (A( 1,3,0,I-7),I=1,9) / > -9.321E+00, 4.562E+00, 2.149E+00, -5.296E-01, 1.829E+00, > -9.652E-01, -2.974E-01, 1.462E-02, 1.821E-01/ C A1 Data (A( 1,0,1,I-7),I=1,9) / > -4.515E-01, -2.743E-01, -9.939E-02, -1.501E-01, -1.524E-01, > -2.689E-01, -2.049E-01, 4.866E-01, 5.096E-01/ Data (A( 1,1,1,I-7),I=1,9) / > 2.268E+00, 1.058E+00, 5.319E-02, -1.398E-01, -2.840E-01, > 8.835E-02, 1.047E-01, 4.373E-02, -4.806E-02/ Data (A( 1,2,1,I-7),I=1,9) / > -4.659E+00, -2.093E+00, -4.036E-02, 1.453E-01, 3.579E-01, > -3.156E-01, -1.856E-01, 6.490E-05, 1.149E-01/ Data (A( 1,3,1,I-7),I=1,9) / > -6.002E+00, 1.057E+00, -3.235E-02, -6.800E-02, -1.649E-01, > 1.416E-01, 6.805E-02, -3.232E-03, -3.050E-02/ C A2 Data (A( 1,0,2,I-7),I=1,9) / > 9.744E+00, 6.499E+00, 6.438E+00, 8.000E+00, 7.726E+00, > 8.729E+00, 4.675E+00, 3.689E+00, 4.248E+00/ Data (A( 1,1,2,I-7),I=1,9) / > -7.684E+00, -1.323E+01, -6.494E+00, -2.230E+00, 1.763E+00, > -3.874E+00, 8.076E-01, 8.567E-01, 8.756E-01/ Data (A( 1,2,2,I-7),I=1,9) / > -6.184E+00, 2.988E+01, 5.910E+00, 2.765E+00, -3.766E+00, > 5.774E+00, 5.451E-01, -3.615E-01, -4.974E-01/ Data (A( 1,3,2,I-7),I=1,9) / > 3.000E+01, -1.204E+01, -9.781E-01, -1.128E+00, 1.972E+00, > -2.638E+00, -3.544E-01, 1.638E-01, 2.315E-01/ C A3 Data (A( 1,0,3,I-7),I=1,9) / > 5.519E+00, 2.281E+00, 6.448E+00, 5.910E+00, 9.341E+00, > 2.520E-01, 4.279E+00, 5.384E+00, 1.907E+00/ Data (A( 1,1,3,I-7),I=1,9) / > 2.331E+01, -1.332E+01, -1.886E+01, 3.300E+00, -1.905E+00, > -2.065E+00, 2.144E+00, -3.247E+00, -1.276E+00/ Data (A( 1,2,3,I-7),I=1,9) / > -1.517E+01, 3.000E+01, 2.608E+01, -2.245E+00, 6.458E+00, > 1.589E+00, -1.869E+00, 1.626E-01, 6.981E-01/ Data (A( 1,3,3,I-7),I=1,9) / > -1.997E+01, 4.412E-01, -1.103E+01, -2.476E-01, -3.418E+00, > -4.926E-01, 2.315E-01, 3.590E-01, -2.576E-01/ C A4 Data (A( 1,0,4,I-7),I=1,9) / > 9.042E-01, 1.997E-05, 1.467E+00, 9.820E-01, 1.080E+00, > 8.014E-02, 1.502E+00, 8.451E-01, 3.204E-01/ Data (A( 1,1,4,I-7),I=1,9) / > 2.833E+00, 3.019E+00, -3.765E+00, -3.735E-02, 4.856E-01, > 4.690E-02, -5.526E-01, -3.161E-02, -3.242E-01/ Data (A( 1,2,4,I-7),I=1,9) / > 2.155E+00, -3.112E+00, 4.400E+00, 1.113E-01, -1.151E+00, > -1.050E-01, 3.402E-01, -1.409E-01, 1.900E-01/ Data (A( 1,3,4,I-7),I=1,9) / > -7.187E+00, 2.021E+00, -1.561E+00, -8.972E-02, 5.864E-01, > 7.088E-02, -1.250E-01, 7.219E-02, -8.404E-02/ C A5 Data (A( 1,0,5,I-7),I=1,9) / > 2.153E+00, 1.518E+00, 1.366E+00, -9.964E-03, -5.589E-02, > -8.690E-01, 1.085E-02, -1.106E-01, -5.631E-02/ Data (A( 1,1,5,I-7),I=1,9) / > 3.748E+00, 4.848E+00, 2.313E+00, 2.361E+00, -2.147E-01, > 2.093E+00, 3.154E+00, 4.732E-01, 1.527E-01/ Data (A( 1,2,5,I-7),I=1,9) / > -8.744E+00, -7.348E+00, -8.032E-01, -1.208E+00, 2.777E+00, > -2.503E+00, -1.977E+00, -1.479E-03, 6.039E-01/ Data (A( 1,3,5,I-7),I=1,9) / > -6.564E+00, 2.685E+00, -2.474E-01, 3.168E-01, -1.452E+00, > 1.053E+00, 6.615E-01, 2.632E-03, -2.291E-01/ C --- Iset = 2, cteq4d --- C A0 Data (A( 2,0,0,I-7),I=1,9) / > 2.150E-01, 1.703E+00, 8.201E-01, -3.231E+00, -2.470E+00, > -1.668E+00, 8.275E-01, 4.716E-01, -2.794E-01/ Data (A( 2,1,0,I-7),I=1,9) / > -8.204E+00, -2.788E+00, -4.418E+00, -1.403E+00, 9.576E-01, > 6.169E-03, -3.391E+00, 2.347E-02, -2.033E-01/ Data (A( 2,2,0,I-7),I=1,9) / > 7.449E+00, -6.452E+00, 1.858E+00, 1.144E+00, -3.477E+00, > 1.595E+00, 1.508E+00, -1.860E-01, -4.064E-01/ Data (A( 2,3,0,I-7),I=1,9) / > -9.321E+00, 4.686E+00, -1.846E+00, -7.860E-01, 1.458E+00, > -1.121E+00, -5.686E-01, 1.714E-03, 1.672E-01/ C A1 Data (A( 2,0,1,I-7),I=1,9) / > -4.515E-01, -2.327E-01, -1.154E-02, -1.925E-01, -2.144E-01, > -3.335E-01, -9.517E-02, 4.781E-01, 4.910E-01/ Data (A( 2,1,1,I-7),I=1,9) / > 2.268E+00, 1.013E+00, -1.372E-01, 9.436E-02, -1.098E-01, > 2.790E-01, -1.731E-01, 4.124E-02, -4.667E-02/ Data (A( 2,2,1,I-7),I=1,9) / > -4.659E+00, -2.047E+00, -3.328E-01, -2.820E-03, 3.170E-01, > -4.136E-01, 6.932E-02, -4.196E-03, 1.154E-01/ Data (A( 2,3,1,I-7),I=1,9) / > -6.002E+00, 1.030E+00, 2.635E-01, -5.136E-02, -1.886E-01, > 1.347E-01, -3.498E-03, -5.172E-03, -2.578E-02/ C A2 Data (A( 2,0,2,I-7),I=1,9) / > 9.744E+00, 6.653E+00, 3.000E+00, 7.167E+00, 7.073E+00, > 7.914E+00, 7.487E+00, 3.539E+00, 3.848E+00/ Data (A( 2,1,2,I-7),I=1,9) / > -7.684E+00, -1.150E+01, 1.258E+00, 1.727E-01, 2.682E+00, > -3.094E+00, -7.178E-01, 8.654E-01, 8.521E-01/ Data (A( 2,2,2,I-7),I=1,9) / > -6.184E+00, 2.979E+01, 5.597E+00, 5.920E-01, -4.443E+00, > 6.929E+00, 8.142E-01, -3.217E-01, -6.070E-01/ Data (A( 2,3,2,I-7),I=1,9) / > 3.000E+01, -1.211E+01, -4.084E+00, -8.347E-01, 1.970E+00, > -4.726E+00, -6.482E-01, 1.616E-01, 2.939E-01/ C A3 Data (A( 2,0,3,I-7),I=1,9) / > 5.519E+00, 2.262E+00, -3.772E+00, 7.437E+00, 1.279E+01, > 5.134E-02, 2.937E+01, 5.777E+00, 1.564E+00/ Data (A( 2,1,3,I-7),I=1,9) / > 2.331E+01, -1.210E+01, 1.128E+01, 3.670E+00, -2.974E+00, > -1.752E+00, -1.698E+01, -2.787E+00, -1.191E+00/ Data (A( 2,2,3,I-7),I=1,9) / > -1.517E+01, 2.997E+01, -1.063E+01, -1.457E+00, 4.117E+00, > 1.314E+00, -8.189E+00, 2.432E-01, 6.901E-01/ Data (A( 2,3,3,I-7),I=1,9) / > -1.997E+01, 1.371E+01, 2.545E+00, -1.993E+00, -2.072E+00, > -5.133E-01, 5.971E+00, 3.276E-01, -2.346E-01/ C A4 Data (A( 2,0,4,I-7),I=1,9) / > 9.042E-01, 3.247E-06, 2.543E+00, 1.105E+00, 1.204E+00, > 1.198E-01, 2.131E+00, 1.126E+00, 3.304E-01/ Data (A( 2,1,4,I-7),I=1,9) / > 2.833E+00, 3.187E+00, -3.697E+00, -4.589E-01, 4.337E-01, > 5.657E-02, -9.004E-01, 5.454E-02, -3.291E-01/ Data (A( 2,2,4,I-7),I=1,9) / > 2.155E+00, -3.206E+00, -4.036E+00, 1.262E+00, -7.363E-01, > -1.725E-01, 4.549E-02, -1.452E-01, 1.942E-01/ Data (A( 2,3,4,I-7),I=1,9) / > -7.187E+00, 1.948E+00, 6.889E+00, -8.539E-01, 2.512E-01, > 1.383E-01, 4.671E-02, 4.105E-02, -8.966E-02/ C A5 Data (A( 2,0,5,I-7),I=1,9) / > 2.153E+00, 1.562E+00, 1.669E+00, 1.272E-01, 3.706E-03, > -8.628E-01, -1.614E-01, -7.141E-02, -1.622E-02/ Data (A( 2,1,5,I-7),I=1,9) / > 3.748E+00, 4.863E+00, -2.675E-01, 2.357E+00, -1.890E-01, > 2.056E+00, 3.243E+00, 4.471E-01, 1.398E-01/ Data (A( 2,2,5,I-7),I=1,9) / > -8.744E+00, -7.453E+00, -1.222E+00, -9.743E-01, 2.957E+00, > -2.318E+00, -2.062E+00, -6.167E-03, 6.155E-01/ Data (A( 2,3,5,I-7),I=1,9) / > -6.564E+00, 2.694E+00, 1.302E+00, 2.347E-01, -1.501E+00, > 8.901E-01, 7.712E-01, 3.041E-03, -2.061E-01/ C --- Iset = 3, cteq4l --- C A0 Data (A( 3,0,0,I-7),I=1,9) / > 2.150E-01, 1.768E+00, 1.146E+00, -2.969E+00, -2.167E+00, > -1.731E+00, -2.050E-01, 3.252E-01, -2.751E-01/ Data (A( 3,1,0,I-7),I=1,9) / > -8.204E+00, -3.115E+00, -3.623E+00, -2.435E+00, 2.815E-01, > -2.647E-01, -2.362E+00, -4.450E-03, -2.182E-01/ Data (A( 3,2,0,I-7),I=1,9) / > 7.449E+00, -6.112E+00, -4.675E+00, 1.664E+00, -3.375E+00, > 1.685E+00, 6.991E-01, -1.563E-01, -3.949E-01/ Data (A( 3,3,0,I-7),I=1,9) / > -9.321E+00, 3.890E+00, 1.818E+00, -5.607E-01, 1.770E+00, > -9.514E-01, -3.397E-01, 2.772E-02, 1.538E-01/ C A1 Data (A( 3,0,1,I-7),I=1,9) / > -4.515E-01, -3.179E-01, -1.307E-01, -2.033E-01, -2.081E-01, > -3.134E-01, -2.909E-01, 4.246E-01, 4.422E-01/ Data (A( 3,1,1,I-7),I=1,9) / > 2.268E+00, 1.078E+00, -3.601E-02, -9.540E-02, -2.513E-01, > 1.098E-01, 1.342E-01, 4.566E-02, -4.939E-02/ Data (A( 3,2,1,I-7),I=1,9) / > -4.659E+00, -2.095E+00, -6.512E-03, 4.800E-02, 2.997E-01, > -3.482E-01, -1.979E-01, -4.276E-03, 1.096E-01/ Data (A( 3,3,1,I-7),I=1,9) / > -6.002E+00, 1.025E+00, -2.844E-03, -3.024E-02, -1.549E-01, > 1.352E-01, 7.481E-02, 7.373E-04, -4.024E-02/ C A2 Data (A( 3,0,2,I-7),I=1,9) / > 9.744E+00, 6.011E+00, 5.358E+00, 6.558E+00, 6.596E+00, > 7.529E+00, 3.649E+00, 3.444E+00, 3.894E+00/ Data (A( 3,1,2,I-7),I=1,9) / > -7.684E+00, -1.406E+01, -5.571E+00, -1.336E+00, 1.809E+00, > -2.841E+00, 1.131E+00, 9.319E-01, 1.029E+00/ Data (A( 3,2,2,I-7),I=1,9) / > -6.184E+00, 2.999E+01, 6.393E+00, 1.239E+00, -3.618E+00, > 5.460E+00, 7.021E-01, -2.827E-01, -4.502E-01/ Data (A( 3,3,2,I-7),I=1,9) / > 3.000E+01, -1.212E+01, -1.931E+00, -1.167E-01, 2.039E+00, > -2.872E+00, -4.563E-01, 1.078E-01, 1.674E-01/ C A3 Data (A( 3,0,3,I-7),I=1,9) / > 5.519E+00, 2.641E+00, 5.130E+00, 4.132E+00, 7.851E+00, > 3.224E-01, 1.774E+00, 6.122E+00, 1.771E+00/ Data (A( 3,1,3,I-7),I=1,9) / > 2.331E+01, -1.482E+01, -1.405E+01, 2.802E+00, -2.067E+00, > -2.069E+00, 2.216E+00, -3.213E+00, -1.183E+00/ Data (A( 3,2,3,I-7),I=1,9) / > -1.517E+01, 2.992E+01, 2.103E+01, -3.721E+00, 4.454E+00, > 1.360E+00, -1.603E+00, 1.893E-02, 7.180E-01/ Data (A( 3,3,3,I-7),I=1,9) / > -1.997E+01, -3.961E+00, -9.137E+00, 1.748E+00, -1.704E+00, > -3.925E-01, 7.639E-01, 2.731E-01, -3.024E-01/ C A4 Data (A( 3,0,4,I-7),I=1,9) / > 9.042E-01, 1.562E-04, 1.246E+00, 9.104E-01, 1.083E+00, > 1.748E-02, 1.716E+00, 1.101E+00, 5.414E-01/ Data (A( 3,1,4,I-7),I=1,9) / > 2.833E+00, 3.161E+00, -3.117E+00, -5.168E-01, 2.725E-01, > 1.695E-01, -1.079E+00, -9.618E-02, -3.572E-01/ Data (A( 3,2,4,I-7),I=1,9) / > 2.155E+00, -3.800E+00, 3.556E+00, 8.017E-01, -9.142E-01, > -1.740E-01, 6.207E-01, -1.656E-01, 1.688E-01/ Data (A( 3,3,4,I-7),I=1,9) / > -7.187E+00, 2.379E+00, -1.214E+00, -2.963E-01, 5.363E-01, > 8.281E-02, -1.559E-01, 8.969E-02, -8.020E-02/ C A5 Data (A( 3,0,5,I-7),I=1,9) / > 2.153E+00, 1.525E+00, 1.294E+00, -2.811E-02, -9.895E-02, > -8.407E-01, 8.477E-02, -1.328E-01, -4.547E-02/ Data (A( 3,1,5,I-7),I=1,9) / > 3.748E+00, 4.916E+00, 2.048E+00, 2.154E+00, -3.882E-01, > 1.995E+00, 3.134E+00, 4.198E-01, 1.228E-01/ Data (A( 3,2,5,I-7),I=1,9) / > -8.744E+00, -7.441E+00, -5.208E-01, -9.996E-01, 2.888E+00, > -2.518E+00, -1.915E+00, 1.369E-02, 5.897E-01/ Data (A( 3,3,5,I-7),I=1,9) / > -6.564E+00, 2.854E+00, -1.391E-01, 1.943E-01, -1.548E+00, > 9.592E-01, 6.638E-01, 4.426E-03, -2.397E-01/ C --- Iset = 4, cteq4a1 --- C A0 Data (A( 4,0,0,I-7),I=1,9) / > 2.150E-01, 1.610E+00, 1.071E+00, -2.881E+00, -2.069E+00, > -1.449E+00, 1.716E-01, 4.394E-01, -2.320E-01/ Data (A( 4,1,0,I-7),I=1,9) / > -8.204E+00, -2.930E+00, -3.695E+00, -2.022E+00, 5.780E-01, > -5.779E-02, -2.521E+00, 3.225E-02, -2.020E-01/ Data (A( 4,2,0,I-7),I=1,9) / > 7.449E+00, -5.052E+00, -3.805E+00, 1.357E+00, -3.462E+00, > 1.626E+00, 7.350E-01, -1.760E-01, -4.004E-01/ Data (A( 4,3,0,I-7),I=1,9) / > -9.321E+00, 2.890E+00, 1.803E+00, -5.611E-01, 1.714E+00, > -1.049E+00, -2.854E-01, 5.502E-03, 1.769E-01/ C A1 Data (A( 4,0,1,I-7),I=1,9) / > -4.515E-01, -2.434E-01, -5.541E-02, -1.741E-01, -1.894E-01, > -3.055E-01, -1.896E-01, 4.879E-01, 4.951E-01/ Data (A( 4,1,1,I-7),I=1,9) / > 2.268E+00, 9.846E-01, -9.211E-03, -3.842E-02, -2.094E-01, > 2.042E-01, 7.888E-02, 4.389E-02, -4.584E-02/ Data (A( 4,2,1,I-7),I=1,9) / > -4.659E+00, -2.235E+00, -6.473E-02, 5.031E-02, 3.148E-01, > -3.953E-01, -1.683E-01, -1.614E-03, 1.114E-01/ Data (A( 4,3,1,I-7),I=1,9) / > -6.002E+00, 1.299E+00, -6.827E-04, -4.552E-02, -1.693E-01, > 1.528E-01, 6.943E-02, -1.697E-04, -2.871E-02/ C A2 Data (A( 4,0,2,I-7),I=1,9) / > 9.744E+00, 6.901E+00, 5.572E+00, 7.914E+00, 7.734E+00, > 8.652E+00, 4.907E+00, 3.768E+00, 4.272E+00/ Data (A( 4,1,2,I-7),I=1,9) / > -7.684E+00, -1.999E+01, -1.604E+00, -1.460E+00, 2.436E+00, > -2.999E+00, 6.776E-01, 8.357E-01, 9.176E-01/ Data (A( 4,2,2,I-7),I=1,9) / > -6.184E+00, 2.627E+01, 9.623E-01, 1.645E+00, -4.533E+00, > 4.505E+00, 7.132E-01, -2.912E-01, -5.329E-01/ Data (A( 4,3,2,I-7),I=1,9) / > 3.000E+01, -4.919E+00, 6.764E-01, -6.002E-01, 2.292E+00, > -2.061E+00, -4.890E-01, 1.309E-01, 2.484E-01/ C A3 Data (A( 4,0,3,I-7),I=1,9) / > 5.519E+00, 2.382E+00, 7.197E+00, 6.455E+00, 1.016E+01, > 2.025E-01, 3.369E+00, 5.288E+00, 1.732E+00/ Data (A( 4,1,3,I-7),I=1,9) / > 2.331E+01, -1.557E+01, -1.989E+01, 1.586E+00, -1.873E+00, > -2.172E+00, 2.698E+00, -2.942E+00, -1.206E+00/ Data (A( 4,2,3,I-7),I=1,9) / > -1.517E+01, 1.228E+01, 2.708E+01, -1.206E+00, 4.455E+00, > 1.614E+00, -1.397E+00, 1.205E-01, 6.584E-01/ Data (A( 4,3,3,I-7),I=1,9) / > -1.997E+01, 3.571E+00, -1.158E+01, -3.500E-01, -2.653E+00, > -4.543E-01, -1.655E-01, 3.775E-01, -2.417E-01/ C A4 Data (A( 4,0,4,I-7),I=1,9) / > 9.042E-01, 7.455E-03, 1.417E+00, 1.017E+00, 1.133E+00, > 8.214E-02, 1.244E+00, 8.882E-01, 3.263E-01/ Data (A( 4,1,4,I-7),I=1,9) / > 2.833E+00, 5.030E+00, -3.415E+00, -3.612E-01, 3.413E-01, > 3.082E-02, -2.106E-01, -2.322E-02, -3.329E-01/ Data (A( 4,2,4,I-7),I=1,9) / > 2.155E+00, -7.241E+00, 4.046E+00, 6.949E-01, -7.593E-01, > -1.013E-01, 5.235E-02, -1.188E-01, 1.884E-01/ Data (A( 4,3,4,I-7),I=1,9) / > -7.187E+00, 4.054E+00, -1.418E+00, -3.777E-01, 3.570E-01, > 6.457E-02, -3.007E-02, 6.322E-02, -8.242E-02/ C A5 Data (A( 4,0,5,I-7),I=1,9) / > 2.153E+00, 1.554E+00, 1.597E+00, 4.035E-02, -8.048E-02, > -8.867E-01, 3.002E-02, -8.075E-02, -4.117E-02/ Data (A( 4,1,5,I-7),I=1,9) / > 3.748E+00, 3.947E+00, 1.812E+00, 2.282E+00, -3.147E-01, > 2.067E+00, 3.211E+00, 4.585E-01, 1.184E-01/ Data (A( 4,2,5,I-7),I=1,9) / > -8.744E+00, -6.952E+00, -8.193E-01, -1.093E+00, 2.859E+00, > -2.380E+00, -1.986E+00, 1.683E-02, 6.080E-01/ Data (A( 4,3,5,I-7),I=1,9) / > -6.564E+00, 3.536E+00, -1.154E-01, 2.481E-01, -1.526E+00, > 9.488E-01, 6.907E-01, 1.300E-02, -2.204E-01/ C --- Iset = 5, cteq4a2 --- C A0 Data (A( 5,0,0,I-7),I=1,9) / > 2.150E-01, 1.412E+00, 6.372E-01, -2.776E+00, -2.012E+00, > -1.417E+00, 1.197E-01, 4.277E-01, -2.204E-01/ Data (A( 5,1,0,I-7),I=1,9) / > -8.204E+00, -3.319E+00, -6.721E+00, -2.223E+00, 3.813E-01, > -3.865E-01, -2.384E+00, 5.576E-02, -2.348E-01/ Data (A( 5,2,0,I-7),I=1,9) / > 7.449E+00, -5.621E+00, 2.675E+00, 1.410E+00, -3.527E+00, > 1.634E+00, 7.471E-01, -1.650E-01, -3.940E-01/ Data (A( 5,3,0,I-7),I=1,9) / > -9.321E+00, 4.265E+00, -8.999E-01, -5.310E-01, 1.826E+00, > -1.033E+00, -2.942E-01, 1.057E-02, 1.839E-01/ C A1 Data (A( 5,0,1,I-7),I=1,9) / > -4.515E-01, -2.442E-01, -5.022E-02, -1.604E-01, -1.659E-01, > -2.770E-01, -2.090E-01, 4.882E-01, 5.059E-01/ Data (A( 5,1,1,I-7),I=1,9) / > 2.268E+00, 9.312E-01, 9.016E-03, -1.145E-01, -2.657E-01, > 1.157E-01, 1.059E-01, 4.444E-02, -4.490E-02/ Data (A( 5,2,1,I-7),I=1,9) / > -4.659E+00, -2.003E+00, -9.800E-02, 1.199E-01, 3.384E-01, > -3.320E-01, -1.870E-01, -6.880E-04, 1.147E-01/ Data (A( 5,3,1,I-7),I=1,9) / > -6.002E+00, 1.076E+00, -3.698E-02, -6.169E-02, -1.615E-01, > 1.413E-01, 6.859E-02, -3.954E-03, -3.134E-02/ C A2 Data (A( 5,0,2,I-7),I=1,9) / > 9.744E+00, 6.008E+00, 3.225E+00, 8.016E+00, 7.696E+00, > 8.675E+00, 4.776E+00, 3.724E+00, 4.305E+00/ Data (A( 5,1,2,I-7),I=1,9) / > -7.684E+00, -1.164E+01, 2.695E+00, -2.156E+00, 2.139E+00, > -3.424E+00, 7.940E-01, 8.600E-01, 8.621E-01/ Data (A( 5,2,2,I-7),I=1,9) / > -6.184E+00, 2.995E+01, 8.608E-01, 2.667E+00, -4.195E+00, > 4.931E+00, 6.018E-01, -3.557E-01, -4.792E-01/ Data (A( 5,3,2,I-7),I=1,9) / > 3.000E+01, -1.389E+01, 4.124E-02, -1.090E+00, 2.169E+00, > -2.177E+00, -4.011E-01, 1.579E-01, 2.271E-01/ C A3 Data (A( 5,0,3,I-7),I=1,9) / > 5.519E+00, 1.240E+00, 5.484E+00, 6.011E+00, 9.120E+00, > 2.203E-01, 4.469E+00, 5.296E+00, 1.878E+00/ Data (A( 5,1,3,I-7),I=1,9) / > 2.331E+01, -1.022E+01, -1.305E+01, 2.662E+00, -1.154E+00, > -2.075E+00, 2.397E+00, -3.169E+00, -1.273E+00/ Data (A( 5,2,3,I-7),I=1,9) / > -1.517E+01, 2.991E+01, 2.157E+01, -1.516E+00, 5.863E+00, > 1.596E+00, -1.952E+00, 1.361E-01, 6.958E-01/ Data (A( 5,3,3,I-7),I=1,9) / > -1.997E+01, -4.268E+00, -8.948E+00, -5.133E-01, -3.569E+00, > -4.785E-01, 8.726E-02, 3.491E-01, -2.590E-01/ C A4 Data (A( 5,0,4,I-7),I=1,9) / > 9.042E-01, 4.137E-04, 5.503E-01, 9.957E-01, 1.090E+00, > 7.702E-02, 1.506E+00, 8.544E-01, 3.265E-01/ Data (A( 5,1,4,I-7),I=1,9) / > 2.833E+00, 3.953E+00, -1.564E+00, -1.112E-01, 4.315E-01, > 4.298E-02, -5.256E-01, -1.954E-02, -3.313E-01/ Data (A( 5,2,4,I-7),I=1,9) / > 2.155E+00, -6.077E+00, 4.401E+00, 2.337E-01, -9.890E-01, > -1.052E-01, 3.321E-01, -1.422E-01, 1.836E-01/ Data (A( 5,3,4,I-7),I=1,9) / > -7.187E+00, 4.288E+00, -2.061E+00, -1.481E-01, 4.993E-01, > 7.263E-02, -1.324E-01, 6.559E-02, -7.926E-02/ C A5 Data (A( 5,0,5,I-7),I=1,9) / > 2.153E+00, 1.469E+00, 2.387E+00, -9.174E-03, -7.001E-02, > -8.526E-01, 2.224E-02, -1.018E-01, -4.756E-02/ Data (A( 5,1,5,I-7),I=1,9) / > 3.748E+00, 4.880E+00, 6.115E-01, 2.336E+00, -2.662E-01, > 2.086E+00, 3.143E+00, 4.673E-01, 1.490E-01/ Data (A( 5,2,5,I-7),I=1,9) / > -8.744E+00, -7.750E+00, -1.337E+00, -1.184E+00, 2.799E+00, > -2.476E+00, -1.971E+00, -1.450E-03, 6.077E-01/ Data (A( 5,3,5,I-7),I=1,9) / > -6.564E+00, 3.124E+00, 2.935E-01, 3.084E-01, -1.484E+00, > 1.031E+00, 6.637E-01, 2.274E-03, -2.322E-01/ C --- Iset = 7, cteq4a4 --- C A0 Data (A( 7,0,0,I-7),I=1,9) / > 2.150E-01, 1.710E+00, 9.927E-01, -2.656E+00, -1.914E+00, > -1.314E+00, 4.729E-02, 3.663E-01, -1.965E-01/ Data (A( 7,1,0,I-7),I=1,9) / > -8.204E+00, -2.550E+00, -4.272E+00, -2.446E+00, 2.434E-01, > -1.175E+00, -2.305E+00, 6.992E-02, -1.984E-01/ Data (A( 7,2,0,I-7),I=1,9) / > 7.449E+00, -6.712E+00, -2.841E+00, 1.396E+00, -3.679E+00, > 1.340E+00, 8.093E-01, -1.604E-01, -4.215E-01/ Data (A( 7,3,0,I-7),I=1,9) / > -9.321E+00, 3.288E+00, 1.611E+00, -3.984E-01, 1.884E+00, > -7.693E-01, -3.394E-01, 4.125E-02, 1.678E-01/ C A1 Data (A( 7,0,1,I-7),I=1,9) / > -4.515E-01, -3.744E-01, -5.302E-02, -1.264E-01, -1.233E-01, > -2.313E-01, -2.143E-01, 4.859E-01, 5.237E-01/ Data (A( 7,1,1,I-7),I=1,9) / > 2.268E+00, 1.430E+00, 2.993E-02, -2.060E-01, -3.311E-01, > -3.545E-03, 1.227E-01, 3.453E-02, -4.949E-02/ Data (A( 7,2,1,I-7),I=1,9) / > -4.659E+00, -2.073E+00, -4.987E-02, 1.648E-01, 3.704E-01, > -2.577E-01, -1.958E-01, -7.167E-04, 1.065E-01/ Data (A( 7,3,1,I-7),I=1,9) / > -6.002E+00, 8.063E-01, -5.234E-02, -5.811E-02, -1.605E-01, > 1.270E-01, 6.892E-02, 2.652E-03, -2.637E-02/ C A2 Data (A( 7,0,2,I-7),I=1,9) / > 9.744E+00, 8.013E+00, 4.971E+00, 7.870E+00, 7.585E+00, > 8.960E+00, 4.576E+00, 3.599E+00, 4.248E+00/ Data (A( 7,1,2,I-7),I=1,9) / > -7.684E+00, -1.069E+01, -2.358E+00, -2.345E+00, 2.012E+00, > -4.774E+00, 7.356E-01, 9.408E-01, 1.013E+00/ Data (A( 7,2,2,I-7),I=1,9) / > -6.184E+00, -1.192E+00, 3.356E+00, 2.929E+00, -4.416E+00, > 7.114E+00, 7.747E-01, -3.030E-01, -6.110E-01/ Data (A( 7,3,2,I-7),I=1,9) / > 3.000E+01, 8.309E+00, 2.911E-01, -1.125E+00, 2.449E+00, > -3.147E+00, -4.799E-01, 7.649E-02, 2.244E-01/ C A3 Data (A( 7,0,3,I-7),I=1,9) / > 5.519E+00, 3.914E+00, 6.677E+00, 4.443E+00, 7.261E+00, > 2.790E-01, 4.503E+00, 5.256E+00, 2.025E+00/ Data (A( 7,1,3,I-7),I=1,9) / > 2.331E+01, -1.570E+01, -1.714E+01, 3.485E+00, -1.408E+00, > -1.752E+00, 2.116E+00, -2.677E+00, -1.183E+00/ Data (A( 7,2,3,I-7),I=1,9) / > -1.517E+01, 7.831E+00, 2.061E+01, -7.123E-01, 8.703E+00, > 1.788E+00, -2.310E+00, -7.359E-02, 6.051E-01/ Data (A( 7,3,3,I-7),I=1,9) / > -1.997E+01, 4.017E+00, -4.575E+00, -1.137E+00, -3.987E+00, > -6.858E-01, 4.937E-01, 1.381E-01, -2.329E-01/ C A4 Data (A( 7,0,4,I-7),I=1,9) / > 9.042E-01, 5.888E-02, 1.292E+00, 8.854E-01, 9.847E-01, > 7.867E-02, 1.569E+00, 7.687E-01, 2.982E-01/ Data (A( 7,1,4,I-7),I=1,9) / > 2.833E+00, 4.588E+00, -2.864E+00, 1.092E-01, 5.946E-01, > -1.766E-02, -6.157E-01, -2.775E-02, -3.089E-01/ Data (A( 7,2,4,I-7),I=1,9) / > 2.155E+00, 1.731E+00, 3.642E+00, -5.463E-02, -1.150E+00, > -6.089E-02, 3.466E-01, -8.542E-02, 2.137E-01/ Data (A( 7,3,4,I-7),I=1,9) / > -7.187E+00, -6.281E+00, -1.056E+00, -1.309E-02, 5.479E-01, > 1.284E-01, -1.181E-01, 2.788E-02, -7.193E-02/ C A5 Data (A( 7,0,5,I-7),I=1,9) / > 2.153E+00, 1.066E+00, 1.836E+00, -3.921E-02, -6.622E-02, > -8.516E-01, 1.048E-02, -1.377E-01, -8.910E-02/ Data (A( 7,1,5,I-7),I=1,9) / > 3.748E+00, 5.853E+00, 1.720E+00, 2.394E+00, -1.667E-01, > 2.271E+00, 3.129E+00, 4.443E-01, 1.231E-01/ Data (A( 7,2,5,I-7),I=1,9) / > -8.744E+00, -5.448E+00, -8.893E-01, -1.296E+00, 2.710E+00, > -2.613E+00, -1.938E+00, -3.045E-03, 6.352E-01/ Data (A( 7,3,5,I-7),I=1,9) / > -6.564E+00, 9.772E-01, -2.502E-01, 3.353E-01, -1.415E+00, > 1.103E+00, 6.310E-01, 8.174E-03, -1.934E-01/ C --- Iset = 8, cteq4a5 --- C A0 Data (A( 8,0,0,I-7),I=1,9) / > 2.150E-01, 1.477E+00, 9.355E-01, -2.717E+00, -1.907E+00, > -1.087E+00, -2.732E-01, 3.323E-01, -1.573E-01/ Data (A( 8,1,0,I-7),I=1,9) / > -8.204E+00, -3.761E+00, -4.462E+00, -3.068E+00, -3.691E-01, > -1.442E+00, -2.021E+00, 4.780E-02, -2.413E-01/ Data (A( 8,2,0,I-7),I=1,9) / > 7.449E+00, -5.198E+00, -2.163E+00, 1.156E+00, -3.769E+00, > 1.160E+00, 9.483E-01, -1.772E-01, -3.961E-01/ Data (A( 8,3,0,I-7),I=1,9) / > -9.321E+00, 2.876E+00, 1.166E+00, 1.278E-01, 1.908E+00, > -6.175E-01, -4.658E-01, -3.998E-04, 2.020E-01/ C A1 Data (A( 8,0,1,I-7),I=1,9) / > -4.515E-01, -2.991E-01, -1.210E-01, -6.356E-02, -6.407E-02, > -2.173E-01, -2.723E-01, 4.929E-01, 5.472E-01/ Data (A( 8,1,1,I-7),I=1,9) / > 2.268E+00, 1.012E+00, 7.550E-02, -6.364E-01, -5.592E-01, > -9.241E-02, 2.045E-01, 3.604E-02, -4.053E-02/ Data (A( 8,2,1,I-7),I=1,9) / > -4.659E+00, -1.829E+00, -6.397E-02, 4.547E-01, 5.402E-01, > -2.653E-01, -2.259E-01, -2.897E-03, 1.116E-01/ Data (A( 8,3,1,I-7),I=1,9) / > -6.002E+00, 9.120E-01, -4.241E-02, -1.222E-01, -1.712E-01, > 1.693E-01, 6.348E-02, -3.005E-03, -3.549E-02/ C A2 Data (A( 8,0,2,I-7),I=1,9) / > 9.744E+00, 6.450E+00, 5.356E+00, 7.735E+00, 6.773E+00, > 1.001E+01, 4.561E+00, 3.561E+00, 4.212E+00/ Data (A( 8,1,2,I-7),I=1,9) / > -7.684E+00, -1.386E+01, -2.857E+00, -4.688E+00, 2.222E+00, > -5.131E+00, 7.469E-01, 9.155E-01, 9.560E-01/ Data (A( 8,2,2,I-7),I=1,9) / > -6.184E+00, 2.991E+01, 3.134E+00, 3.242E+00, -4.347E+00, > 5.097E+00, 5.081E-01, -4.114E-01, -6.051E-01/ Data (A( 8,3,2,I-7),I=1,9) / > 3.000E+01, -1.664E+01, 2.038E-01, -4.521E-01, 1.932E+00, > -1.689E+00, -3.263E-01, 1.287E-01, 2.425E-01/ C A3 Data (A( 8,0,3,I-7),I=1,9) / > 5.519E+00, 1.863E+00, 6.824E+00, 3.364E+00, 3.076E+00, > 5.175E-01, 6.560E+00, 5.635E+00, 2.065E+00/ Data (A( 8,1,3,I-7),I=1,9) / > 2.331E+01, -8.006E+00, -1.602E+01, -2.125E-01, 3.115E+00, > -2.111E+00, -9.584E-01, -3.303E+00, -1.287E+00/ Data (A( 8,2,3,I-7),I=1,9) / > -1.517E+01, 2.543E+01, 1.902E+01, -7.670E-01, 9.041E+00, > 1.821E+00, -2.198E+00, 1.162E-01, 6.407E-01/ Data (A( 8,3,3,I-7),I=1,9) / > -1.997E+01, -1.788E+01, -5.294E+00, -1.001E-01, -5.534E+00, > -5.069E-01, 1.102E+00, 2.762E-01, -2.714E-01/ C A4 Data (A( 8,0,4,I-7),I=1,9) / > 9.042E-01, 1.240E-02, 1.208E+00, 6.314E-01, 5.906E-01, > 9.945E-02, 1.463E+00, 7.404E-01, 2.787E-01/ Data (A( 8,1,4,I-7),I=1,9) / > 2.833E+00, 3.141E+00, -2.311E+00, 2.764E-02, 1.110E+00, > -3.865E-02, -6.371E-01, -5.847E-02, -3.454E-01/ Data (A( 8,2,4,I-7),I=1,9) / > 2.155E+00, -5.401E+00, 3.006E+00, -1.049E+00, -1.692E+00, > -1.359E-01, 5.236E-01, -1.499E-01, 2.322E-01/ Data (A( 8,3,4,I-7),I=1,9) / > -7.187E+00, 3.023E+00, -8.972E-01, 5.384E-01, 6.988E-01, > 8.255E-02, -2.219E-01, 3.268E-02, -5.544E-02/ C A5 Data (A( 8,0,5,I-7),I=1,9) / > 2.153E+00, 1.524E+00, 1.652E+00, 5.648E-02, -7.828E-03, > -1.086E+00, 3.816E-02, -1.507E-01, -9.609E-02/ Data (A( 8,1,5,I-7),I=1,9) / > 3.748E+00, 5.347E+00, 1.690E+00, 2.110E+00, -1.283E-01, > 2.401E+00, 3.061E+00, 4.810E-01, 1.761E-01/ Data (A( 8,2,5,I-7),I=1,9) / > -8.744E+00, -7.479E+00, -8.646E-01, -2.083E+00, 2.646E+00, > -2.623E+00, -1.898E+00, -2.853E-03, 6.331E-01/ Data (A( 8,3,5,I-7),I=1,9) / > -6.564E+00, 3.218E+00, -1.659E-01, 7.016E-01, -1.190E+00, > 1.170E+00, 5.987E-01, -4.004E-03, -2.064E-01/ C --- Iset = 11, cteq4hq --- C A0 Data (A(11,0,0,I-7),I=1,9) / > 2.150E-01, 1.811E+00, 1.031E+00, -2.809E+00, -2.021E+00, > -1.374E+00, 6.892E-03, 3.977E-01, -1.897E-01/ Data (A(11,1,0,I-7),I=1,9) / > -8.204E+00, -3.033E+00, -4.563E+00, -2.176E+00, 4.190E-01, > -2.590E-01, -2.339E+00, 3.239E-02, -2.336E-01/ Data (A(11,2,0,I-7),I=1,9) / > 7.449E+00, -4.856E+00, -3.614E+00, 1.295E+00, -3.527E+00, > 1.586E+00, 8.438E-01, -1.727E-01, -4.009E-01/ Data (A(11,3,0,I-7),I=1,9) / > -9.321E+00, 2.487E+00, 1.900E+00, -4.228E-01, 1.806E+00, > -9.920E-01, -3.614E-01, 1.408E-02, 1.704E-01/ C A1 Data (A(11,0,1,I-7),I=1,9) / > -4.515E-01, -2.931E-01, -2.177E-02, -1.676E-01, -1.769E-01, > -2.893E-01, -2.238E-01, 4.897E-01, 5.160E-01/ Data (A(11,1,1,I-7),I=1,9) / > 2.268E+00, 1.007E+00, -3.209E-02, -1.005E-01, -2.528E-01, > 1.499E-01, 1.258E-01, 4.155E-02, -4.962E-02/ Data (A(11,2,1,I-7),I=1,9) / > -4.659E+00, -2.107E+00, -4.628E-02, 8.924E-02, 3.345E-01, > -3.635E-01, -1.859E-01, 7.056E-04, 1.152E-01/ Data (A(11,3,1,I-7),I=1,9) / > -6.002E+00, 1.218E+00, -2.528E-03, -4.699E-02, -1.632E-01, > 1.524E-01, 6.140E-02, 5.708E-04, -3.143E-02/ C A2 Data (A(11,0,2,I-7),I=1,9) / > 9.744E+00, 6.670E+00, 5.201E+00, 7.950E+00, 7.664E+00, > 8.847E+00, 4.534E+00, 3.677E+00, 4.288E+00/ Data (A(11,1,2,I-7),I=1,9) / > -7.684E+00, -1.432E+01, -4.503E+00, -2.350E+00, 2.330E+00, > -3.517E+00, 6.804E-01, 8.597E-01, 9.503E-01/ Data (A(11,2,2,I-7),I=1,9) / > -6.184E+00, 2.986E+01, 5.298E+00, 2.846E+00, -4.544E+00, > 4.550E+00, 9.484E-01, -3.528E-01, -6.182E-01/ Data (A(11,3,2,I-7),I=1,9) / > 3.000E+01, -1.426E+01, -1.227E+00, -1.111E+00, 2.355E+00, > -1.869E+00, -5.950E-01, 1.508E-01, 2.729E-01/ C A3 Data (A(11,0,3,I-7),I=1,9) / > 5.519E+00, 2.512E+00, 8.346E+00, 5.386E+00, 8.413E+00, > 1.845E-01, 4.516E+00, 5.324E+00, 1.981E+00/ Data (A(11,1,3,I-7),I=1,9) / > 2.331E+01, -1.249E+01, -1.776E+01, 1.879E+00, -1.395E+00, > -2.146E+00, 2.514E+00, -3.100E+00, -1.267E+00/ Data (A(11,2,3,I-7),I=1,9) / > -1.517E+01, 2.276E+01, 2.335E+01, -1.023E-01, 6.549E+00, > 1.653E+00, -2.655E+00, 2.154E-01, 7.043E-01/ Data (A(11,3,3,I-7),I=1,9) / > -1.997E+01, -6.234E+00, -1.010E+01, -1.196E+00, -3.786E+00, > -4.753E-01, 4.773E-01, 3.073E-01, -2.621E-01/ C A4 Data (A(11,0,4,I-7),I=1,9) / > 9.042E-01, 6.567E-04, 1.076E+00, 9.574E-01, 1.074E+00, > 8.488E-02, 1.564E+00, 8.212E-01, 3.222E-01/ Data (A(11,1,4,I-7),I=1,9) / > 2.833E+00, 2.478E+00, -1.963E+00, -9.318E-02, 4.307E-01, > 2.164E-02, -5.140E-01, -4.187E-02, -3.172E-01/ Data (A(11,2,4,I-7),I=1,9) / > 2.155E+00, -3.172E+00, 2.047E+00, 2.075E-01, -9.315E-01, > -9.271E-02, 1.660E-01, -1.297E-01, 1.995E-01/ Data (A(11,3,4,I-7),I=1,9) / > -7.187E+00, 2.353E+00, -5.826E-01, -1.252E-01, 4.568E-01, > 6.032E-02, -4.152E-02, 7.040E-02, -8.115E-02/ C A5 Data (A(11,0,5,I-7),I=1,9) / > 2.153E+00, 1.584E+00, 2.169E+00, 2.561E-02, -6.093E-02, > -8.609E-01, 4.342E-02, -1.030E-01, -7.603E-02/ Data (A(11,1,5,I-7),I=1,9) / > 3.748E+00, 4.066E+00, 1.637E+00, 2.308E+00, -2.808E-01, > 2.085E+00, 3.133E+00, 4.741E-01, 1.600E-01/ Data (A(11,2,5,I-7),I=1,9) / > -8.744E+00, -7.020E+00, -8.614E-01, -1.210E+00, 2.802E+00, > -2.444E+00, -1.949E+00, 1.103E-02, 6.071E-01/ Data (A(11,3,5,I-7),I=1,9) / > -6.564E+00, 3.731E+00, -7.693E-02, 2.984E-01, -1.485E+00, > 1.004E+00, 6.416E-01, 1.481E-02, -2.140E-01/ C ======================== IRT2 = 0 IF (ISET.EQ.LSET .AND. QSTO.EQ.QQ) THEN IF (IFL.EQ.LFLV) GOTO 100 IF (IFL.GE.-3 .AND. LFLV.GE.-3) GOTO 501 ENDIF C return 0 if below threshold C or setup threshold Q => Qi IF (IFL.LE.-4) THEN KFL=ABS(IFL) IF (QQ .LE. QMS(KFL)) THEN CTQ4DF = 0.0 RETURN ENDIF QI = QMS(KFL) ELSE QI = QMN(ISET) ENDIF C setup Q values SBL = LOG(QQ/ALAM(ISET)) / LOG(QI/ALAM(ISET)) SB = LOG(SBL) SB2 = SB*SB SB3 = SB2*SB C evaluate exponents 501 DO I=0,IEX AF(I) = A(Iset, 0,I,IFL) + A(Iset, 1,I,IFL) *SB > + A(Iset, 2,I,IFL) *SB2 + A(Iset, 3,I,IFL) *SB3 ENDDO C heavy quark normalization factor has distinct fn form IF (IFL.LE.-4) THEN zz = >A(Iset,1,0,IFL) + A(Iset,2,0,IFL) *SB + A(Iset,3,0,IFL) *SB2 AF(0) = SB **A(Iset, 0,0,IFL) * EXP(zz) ELSE AF(0) = EXP(AF(0)) ENDIF C Evaluate PDF function 100 CTQ4DF = AF(0) * X **AF(1) * (1.-X) **AF(2) > * (1D0 + AF(3) * X **AF(4)) * (LOG(1. + 1./X)) **AF(5) IF (CTQ4DF.LT.0.0) THEN CTQ4DF = 0.0 IRT2 = 1 ENDIF LSET = ISET LFLV = IFL QSTO = QQ RETURN C ENTRY WLAMD4 (ISET, IORDER, NEFF) C IORDER = IORD (ISET) C WLAMD4 = VLM (NEFF,ISET) C RETURN c ************************************ END