PDF/spdf/ctq1pd.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:32:14  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.2  1996/10/29 09:30:56  pdflib
   9 * more LINNUX corrections in DATA statements
  10 *
  11 * Revision 1.1.1.1  1996/04/12 15:29:07  plothow
  12 * Version 7.01
  13 *
  14 *
  15 #include "pdf/pilot.h"
  16       FUNCTION Ctq1Pd (ISET, IPRTN, XX, QQ, IRT)
  17
  18 C            Returns xf(x,Q) -- the momentum fraction distribution !!
  19 C            Returns valence and sea rather than combined flavor distr.
  20
  21 C            Iset : PDF set label
  22
  23 C            Iprtn  : Parton label:   2, 1 = d_ and u_ valence
  24 C                                     0 = gluon
  25 C                            -1, ... -6 = u, d, s, c, b, t sea quarks
  26
  27 C            XX  : Bjorken-x
  28 C            QQ  : scale parameter "Q"
  29 C            Irt : Return code
  30
  31 C     IMPLICIT DOUBLE PRECISION (A-H, O-Z)
  32 C+SEQ, IMPDP.
  33       PARAMETER (D0=0D0, D1=1D0, D2=2D0, D3=3D0, D4=4D0, D10=1D1)
  34       PARAMETER (NEX = 5, MXFL = 6, NPN = 3, NST = 30, NEXPT=20)
  35       PARAMETER (NST4 = NST*4)
  36
  37 #include "pdf/expdp.inc"
  38      +       P012, P34, P5, PDFD
  39       REAL PDF
  40
  41       DIMENSION
  42      >   IORD(NST), ISCH(NST), NQRK(NST),ALM(NST)
  43      >   , VLM(4:6,NST), QMS(4:6, NST)
  44      >   , XMN(NST), QMN(NST), QMX(NST), NEXP(NEXPT)
  45      >   , MEX(NST), MPN(NST), EXPN(NEXPT, NST), EXPNOR(NEXPT)
  46
  47 C                                                             CTEQ1M
  48       DATA
  49      >  ISCH(1), IORD(1), NQRK(1), ALM(1) /  1,  2,  6,  .152 /
  50      >  ,(VLM(I,1), I=4,6) / .231,    .152,    .059  /
  51      >  ,(QMS(I,1), I=4,6) / 1.50,   5.00,  180.0 /
  52      >  ,XMN(1), QMN(1), QMX(1) /  1.E-5,  2.00,  1.E3  /
  53      >  ,MEX(1), MPN(1), NEXP(1) /  5, 3, 8  /
  54      >  ,(EXPN(I, 1), I=1,8)
  55      >  / 0.989, 1.00, 1.02, 0.978, 1.10, 0.972, 0.987, 0.846 /
  56 C                                                             CTEQ1MS
  57       DATA
  58      >  ISCH(2), IORD(2), NQRK(2), ALM(2) /  1,  2,  6, .152  /
  59      >  ,(VLM(I,2), I=4,6) / .231,    .152,    .059  /
  60      >  ,(QMS(I,2), I=4,6) / 1.50,   5.00,  180.0 /
  61      >  ,XMN(2), QMN(2), QMX(2) /  1.E-5,  2.00,  1.E3  /
  62      >  ,MEX(2), MPN(2), NEXP(2) /  5, 3, 8  /
  63      >  ,(EXPN(I, 2), I=1,8 )
  64      >  / 0.989, 1.00, 1.02, 0.984, 1.05, 0.891, 0.923, 0.824 /
  65 C                                                             CTEQ1ML
  66       DATA
  67      >  ISCH(3), IORD(3), NQRK(3), ALM(3) /  1,  2,  6, .220  /
  68      >  ,(VLM(I,3), I=4,6) / .322,    .220,     .088  /
  69      >  ,(QMS(I,3), I=4,6) / 1.50,   5.00,  180.0 /
  70      >  ,XMN(3), QMN(3), QMX(3) /  1.E-5,  2.00,  1.E3  /
  71      >  ,MEX(3), MPN(3), NEXP(3) /  5, 3, 8 /
  72      >  ,(EXPN(I, 3), I=1,8 )
  73      >  / 0.985, 1.00, 1.01, 0.977, 1.07, 1.31, 1.19, 1.09 /
  74
  75 C                                                             CTEQ1D
  76       DATA
  77      >  ISCH(4), IORD(4), NQRK(4), ALM(4) /  2,  2,  6, .164  /
  78      >  ,(VLM(I,4), I=4,6) / .247,    .164,    .064  /
  79      >  ,(QMS(I,4), I=4,6) / 1.50,   5.00,  180.0 /
  80      >  ,XMN(4), QMN(4), QMX(4) /  1.E-5,  2.00,  1.E3  /
  81      >  ,MEX(4), MPN(4), NEXP(4) /  5, 3, 8 /
  82      >  ,(EXPN(I, 4), I=1,8 )
  83      >  / 0.983, 1.00, 1.01, 0.975, 0.964, 1.23, 1.00, 1.12 /
  84 C                                                             CTEQ1L
  85       DATA
  86      >  ISCH(5), IORD(5), NQRK(5), ALM(5) /  1,  1,  6, .125  /
  87      >  ,(VLM(I,5), I=4,6) / .168,    .125,     .063   /
  88      >  ,(QMS(I,5), I=4,6) / 1.50,   5.00,  180.0 /
  89      >  ,XMN(5), QMN(5), QMX(5) /  1.E-5,  2.00,  1.E3  /
  90      >  ,MEX(5), MPN(5), NEXP(5) /  5, 3, 8  /
  91      >  ,(EXPN(I, 5), I=1,8 )
  92      >  / 0.982, 1.01, 1.00, 0.972, 0.840, 0.959, 0.930, 0.861 /
  93
  94       DATA IST, LP, QSTO, ALN2 / 0, -10, 1.2345, 0.6931 /
  95
  96       X  = XX
  97       QQ = MAX(QQ,QMN(ISET))
  98       IF(ISET.EQ.IST.AND.IPRTN.EQ.LP.AND.QSTO.EQ.QQ) GOTO 100
  99
 100       IRT = 0
 101       IP = ABS(IPRTN)
 102       IF (IP.GE.5.AND.QQ.LE.QMS(IP, ISET)) THEN
 103         Ctq1Pd = 0.0
 104         RETURN
 105       ENDIF
 106
 107 C if heavy parton, different logarithmic scale
 108
 109       IF(IP.GE.5) THEN
 110         QI = QMS(IP,ISET)
 111       ELSE
 112         QI   = QMN (ISET)
 113       ENDIF
 114
 115       ALAM = ALM (ISET)
 116       STA = LOG(QQ/ALAM)
 117       STB = LOG(QI/ALAM)
 118
 119 C      SBL = LOG(QQ/Alam) / LOG(Qi/Alam)
 120       SBL = STA/STB
 121       SB = LOG (SBL)
 122       SB2 = SB*SB
 123       SB3 = SB2*SB
 124
 125       IFLV = 3 - IPRTN
 126
 127       GOTO (1, 2, 3, 4, 5), ISET
 128
 129     1 GOTO(11,12,13,14,15,16,17,18,19)IFLV
 130 C   ifl =     2
 131    11 A0=0.3636E+01*(1.0 + 0.3122E+00*SB+0.1396E+00*SB2+0.4251E+00*SB3)
 132       A1=0.6930E+00-.2574E-01*SB+0.1047E+00*SB2-.2794E-01*SB3
 133       A2=0.3195E+01+0.4045E+00*SB-.3737E+00*SB2-.1677E+00*SB3
 134       A3=0.1009E+00*(1.0 -.1784E+01*SB+0.6263E+00*SB2+0.7337E-01*SB3)
 135      $  -1.0
 136       A4=0.2910E+00-.2793E+00*SB+0.6155E-01*SB2+0.5150E-02*SB3
 137       A5=0.0000E+00+0.3185E+00*SB+0.1953E+00*SB2+0.4184E-01*SB3
 138       GOTO 100
 139 C   ifl =     1
 140    12 A0=0.2851E+00*(1.0 + 0.3617E+00*SB-.4526E+00*SB2+0.5787E-01*SB3)
 141       A1=0.2690E+00+0.1104E-01*SB+0.1888E-01*SB2-.1031E-01*SB3
 142       A2=0.3766E+01+0.7850E+00*SB-.3053E+00*SB2+0.1822E+00*SB3
 143       A3=0.2865E+02*(1.0 -.9774E+00*SB+0.5958E+00*SB2-.1234E+00*SB3)
 144      $  -1.0
 145       A4=0.8230E+00-.3612E+00*SB+0.5520E-01*SB2+0.1571E-01*SB3
 146       A5=0.0000E+00+0.2145E-01*SB+0.2289E+00*SB2-.4947E-01*SB3
 147       GOTO 100
 148 C   ifl =     0
 149    13 A0=0.2716E+01*(1.0 -.2092E+01*SB+0.1500E+01*SB2-.3703E+00*SB3)
 150       A1=-.3100E-01-.7963E+00*SB+0.1129E+01*SB2-.4191E+00*SB3
 151       A2=0.8015E+01+0.1168E+01*SB-.1625E+01*SB2-.1130E+01*SB3
 152       A3=0.4813E+02*(1.0 -.4951E+00*SB-.8715E+00*SB2+0.5893E+00*SB3)
 153      $  -1.0
 154       A4=0.2773E+01-.6329E+00*SB-.1048E+01*SB2+0.1418E+00*SB3
 155       A5=0.0000E+00+0.5048E+00*SB+0.2390E+01*SB2-.4159E+00*SB3
 156       GOTO 100
 157 C   ifl =    -1
 158    14 A0=0.3085E+00*(1.0 + 0.9422E+00*SB-.2606E+01*SB2+0.1364E+01*SB3)
 159       A1=0.5000E-02-.6433E+00*SB+0.4980E+00*SB2-.1780E+00*SB3
 160       A2=0.7490E+01+0.9112E+00*SB-.2047E+01*SB2+0.1456E+01*SB3
 161       A3=0.1145E-01*(1.0 + 0.4610E+01*SB+0.1699E+01*SB2+0.1296E+00*SB3)
 162      $  -1.0
 163       A4=0.6030E+00-.8081E+00*SB+0.9410E+00*SB2-.4458E+00*SB3
 164       A5=0.0000E+00-.1736E+01*SB+0.2863E+01*SB2-.1268E+01*SB3
 165       GOTO 100
 166 C   ifl =    -2
 167    15 A0=0.1324E+00*(1.0 -.1050E+01*SB+0.4844E+00*SB2-.1043E+00*SB3)
 168       A1=-.1580E+00+0.1672E+00*SB-.4100E+00*SB2+0.1793E+00*SB3
 169       A2=0.8559E+01-.7351E-01*SB+0.5898E+00*SB2-.2655E+00*SB3
 170       A3=0.2378E+02*(1.0 -.1108E+00*SB-.1646E-01*SB2+0.1129E-01*SB3)
 171      $  -1.0
 172       A4=0.1477E+01+0.3312E-01*SB-.2191E+00*SB2+0.9588E-01*SB3
 173       A5=0.0000E+00+0.1850E+01*SB-.1481E+01*SB2+0.6222E+00*SB3
 174       GOTO 100
 175 C   ifl =    -3
 176    16 A0=0.3208E+00*(1.0 -.4755E+00*SB-.4003E+00*SB2+0.2300E+00*SB3)
 177       A1=-.3200E-01-.3357E+00*SB+0.3222E-01*SB2+0.5011E-01*SB3
 178       A2=0.1164E+02+0.1048E+01*SB-.1097E+01*SB2-.4431E+00*SB3
 179       A3=0.5065E+02*(1.0 + 0.2484E+00*SB-.9235E+00*SB2+0.1935E+00*SB3)
 180      $  -1.0
 181       A4=0.3300E+01-.6785E+00*SB+0.5337E+00*SB2-.4035E+00*SB3
 182       A5=0.0000E+00-.2496E+00*SB+0.3903E+00*SB2+0.1392E+00*SB3
 183       GOTO 100
 184 C   ifl =    -4
 185    17 A0=0.7967E-06*(1.0 + 0.1587E+01*SB+0.1812E+02*SB2-.1333E+02*SB3)
 186      $  *SQRT(STA - STB)
 187       A1=0.1096E+01-.1236E+01*SB+0.1014E+02*SB2+0.1940E+01*SB3
 188       A2=0.4366E+00+0.1197E+02*SB-.5471E+00*SB2-.5427E+01*SB3
 189       A3=0.4650E+03*(1.0 + 0.1310E+02*SB-.1918E+02*SB2+0.6791E+01*SB3)
 190      $  -1.0
 191       A4=-.8486E+00+0.7457E+00*SB-.1083E+02*SB2-.1210E+01*SB3
 192       A5=0.3494E+01-.3511E+01*SB-.1766E+01*SB2+0.3442E+01*SB3
 193       GOTO 100
 194 C   ifl =    -5
 195    18 A0=0.1713E-03*(1.0 + 0.2562E+02*SB-.2988E+02*SB2+0.4798E+01*SB3)
 196      $  *SQRT(STA - STB)
 197       A1=-.5276E-01+0.4105E+00*SB-.1079E+01*SB2+0.6278E+00*SB3
 198       A2=0.4515E+01+0.8369E+01*SB-.1192E+02*SB2+0.3403E+01*SB3
 199       A3=0.1756E+01*(1.0 + 0.1325E+02*SB-.2997E+02*SB2+0.1758E+02*SB3)
 200      $  -1.0
 201       A4=0.3557E-01+0.4159E+01*SB-.6947E+01*SB2+0.2982E+01*SB3
 202       A5=0.2551E+01+0.2168E+01*SB-.5119E+01*SB2+0.3739E+01*SB3
 203       GOTO 100
 204 C   ifl =    -6
 205    19 A0=0.7510E-04*(1.0 + 0.2836E+02*SB-.3000E+02*SB2-.2979E+02*SB3)
 206      $  *SQRT(STA - STB)
 207       A1=-.1855E+00+0.4543E+00*SB-.1448E+01*SB2+0.2009E-01*SB3
 208       A2=0.6775E+01-.4210E+01*SB-.1221E+01*SB2+0.1199E+02*SB3
 209       A3=0.1070E+01*(1.0 + 0.8356E+01*SB-.2992E+02*SB2+0.2433E+02*SB3)
 210      $  -1.0
 211       A4=-.4601E-01+0.4248E+01*SB-.1736E+01*SB2+0.1187E+02*SB3
 212       A5=0.2771E+01+0.1382E+01*SB-.4797E+01*SB2+0.1273E+01*SB3
 213       GOTO 100
 214
 215
 216     2 GOTO(21,22,23,24,25,26,27,28,29)IFLV
 217 C                                                             CTEQ1MS
 218 C   ifl =     2
 219    21 A0=0.1828E+01*(1.0 -.8698E+00*SB+0.2906E+00*SB2-.2003E-01*SB3)
 220       A1=0.6060E+00+0.8595E-01*SB-.4934E-01*SB2+0.2221E-01*SB3
 221       A2=0.3454E+01-.3115E+00*SB+0.1321E+01*SB2-.3490E+00*SB3
 222       A3=0.2616E+00*(1.0 -.1670E+01*SB+0.2333E+01*SB2+0.7730E-01*SB3)
 223      $  -1.0
 224       A4=0.8920E+00-.8500E-02*SB+0.4960E+00*SB2-.4045E-01*SB3
 225       A5=0.0000E+00+0.1091E+01*SB-.1613E+00*SB2+0.3773E-01*SB3
 226       GOTO 100
 227 C   ifl =     1
 228    22 A0=0.2885E+00*(1.0 + 0.3388E+00*SB-.4550E+00*SB2+0.6005E-01*SB3)
 229       A1=0.2730E+00+0.1198E-01*SB+0.1880E-01*SB2-.1077E-01*SB3
 230       A2=0.3736E+01+0.7687E+00*SB-.2731E+00*SB2+0.1638E+00*SB3
 231       A3=0.2741E+02*(1.0 -.9585E+00*SB+0.5925E+00*SB2-.1239E+00*SB3)
 232      $  -1.0
 233       A4=0.8040E+00-.3546E+00*SB+0.6123E-01*SB2+0.1086E-01*SB3
 234       A5=0.0000E+00+0.4277E-01*SB+0.2187E+00*SB2-.4646E-01*SB3
 235       GOTO 100
 236 C   ifl =     0
 237    23 A0=0.8416E-01*(1.0 -.1996E+01*SB+0.1903E+01*SB2-.6722E+00*SB3)
 238       A1=-.4790E+00-.5459E+00*SB+0.1638E+01*SB2-.4342E+00*SB3
 239       A2=0.5071E+01+0.1470E+01*SB-.2401E+01*SB2+0.1273E+01*SB3
 240       A3=0.2847E+02*(1.0 + 0.1124E+00*SB-.1338E+01*SB2+0.7115E+00*SB3)
 241      $  -1.0
 242       A4=0.4990E+00-.7208E+00*SB+0.3333E-03*SB2-.2354E+00*SB3
 243       A5=0.0000E+00-.4480E+00*SB+0.3720E+01*SB2-.1838E+01*SB3
 244       GOTO 100
 245 C   ifl =    -1
 246    24 A0=0.4378E+00*(1.0 -.1244E+01*SB+0.3278E+01*SB2-.2098E+01*SB3)
 247       A1=0.3500E-01-.1298E+01*SB+0.1229E+01*SB2-.3665E+00*SB3
 248       A2=0.6781E+01+0.4078E+01*SB-.9711E+00*SB2-.1536E+01*SB3
 249       A3=0.1527E-03*(1.0 + 0.1430E+02*SB+0.3000E+02*SB2+0.2771E+02*SB3)
 250      $  -1.0
 251       A4=0.3060E+00+0.1011E+01*SB-.2045E+01*SB2+0.9422E+00*SB3
 252       A5=0.0000E+00-.3205E+01*SB+0.2683E+01*SB2-.1746E+00*SB3
 253       GOTO 100
 254 C   ifl =    -2
 255    25 A0=0.7413E-01*(1.0 + 0.1291E+01*SB-.2667E+01*SB2+0.1076E+01*SB3)
 256       A1=-.2730E+00-.1206E+00*SB+0.1828E+00*SB2-.1001E+00*SB3
 257       A2=0.7719E+01+0.1537E+01*SB-.6410E+00*SB2-.3920E-01*SB3
 258       A3=0.1799E+02*(1.0 -.1334E+01*SB+0.1916E+01*SB2-.8878E+00*SB3)
 259      $  -1.0
 260       A4=0.1167E+01-.9176E-01*SB+0.5132E+00*SB2-.3460E+00*SB3
 261       A5=0.0000E+00-.5023E+00*SB+0.1951E+01*SB2-.8427E+00*SB3
 262       GOTO 100
 263 C   ifl =    -3
 264    26 A0=0.6551E+00*(1.0 -.5968E-01*SB+0.5621E-02*SB2-.2074E+00*SB3)
 265       A1=0.2800E-01-.1138E+01*SB+0.1178E+01*SB2-.4425E+00*SB3
 266       A2=0.7553E+01+0.3996E+01*SB-.4448E+01*SB2+0.1673E+01*SB3
 267       A3=0.9264E-01*(1.0 -.1760E+01*SB+0.1634E+01*SB2-.4067E+00*SB3)
 268      $  -1.0
 269       A4=0.1970E+00+0.5256E+00*SB-.9775E+00*SB2+0.4488E+00*SB3
 270       A5=0.0000E+00-.3668E+01*SB+0.4757E+01*SB2-.1717E+01*SB3
 271       GOTO 100
 272 C   ifl =    -4
 273    27 A0=0.1486E-03*(1.0 + 0.2107E+01*SB-.1056E+02*SB2+0.1403E+02*SB3)
 274      $  * SQRT(STA - STB)
 275       A1=0.2115E+00-.1702E+01*SB+0.2571E+01*SB2-.1177E+01*SB3
 276       A2=0.3533E+01+0.1367E+01*SB-.3397E+01*SB2+0.6260E+01*SB3
 277       A3=0.1096E+02*(1.0 + 0.9213E+01*SB-.2020E+02*SB2+0.1084E+02*SB3)
 278      $  -1.0
 279       A4=0.7041E+00-.7236E+00*SB+0.2766E-01*SB2+0.7352E+00*SB3
 280       A5=0.3904E+01-.4398E+01*SB+0.7056E+01*SB2-.3722E+01*SB3
 281       GOTO 100
 282 C   ifl =    -5
 283    28 A0=0.1201E-03*(1.0 + 0.5408E+01*SB-.1489E+02*SB2+0.1667E+02*SB3)
 284      $  * SQRT(STA - STB)
 285       A1=0.1420E-01-.1525E+01*SB+0.2408E+01*SB2-.1154E+01*SB3
 286       A2=0.4254E+01+0.2836E+01*SB-.6018E+00*SB2+0.4133E+00*SB3
 287       A3=0.5696E+01*(1.0 + 0.9451E+01*SB-.2029E+02*SB2+0.1033E+02*SB3)
 288      $  -1.0
 289       A4=0.4775E+00-.6695E+00*SB+0.2747E+00*SB2-.1051E+00*SB3
 290       A5=0.3330E+01-.5133E+01*SB+0.6921E+01*SB2-.3283E+01*SB3
 291       GOTO 100
 292 C   ifl =    -6
 293    29 A0=0.7697E-04*(1.0 + 0.2801E+02*SB-.1901E+02*SB2-.2880E+02*SB3)
 294      $  *SQRT(STA - STB)
 295       A1=-.2249E+00+0.4432E+00*SB-.1454E+01*SB2+0.3509E-01*SB3
 296       A2=0.6642E+01-.2702E+01*SB+0.8229E+01*SB2+0.8243E+01*SB3
 297       A3=0.1146E+01*(1.0 + 0.8104E+01*SB-.2998E+02*SB2+0.2812E+02*SB3)
 298      $  -1.0
 299       A4=-.6421E-01+0.4246E+01*SB-.2908E+01*SB2+0.9686E-02*SB3
 300       A5=0.2606E+01+0.1261E+01*SB-.4933E+01*SB2+0.3476E+00*SB3
 301       GOTO 100
 302
 303
 304     3 GOTO(31,32,33,34,35,36,37,38,39)IFLV
 305 C                                                             CTEQ1ML
 306 C   ifl =     2
 307    31 A0=0.3777E+01*(1.0 + 0.6986E+00*SB-.20655E+01*SB2+.10334E+01*SB3)
 308       A1=0.7100E+00+.2880E-01*SB-.7930E-01*SB2+0.5600E-01*SB3
 309       A2=0.3259E+01+0.1508E+01*SB-.3932E+01*SB2+0.20613E+01*SB3
 310       A3=0.1304E+00*(1.0 -.2016E+00*SB-.30015E+01*SB2+0.19118E+01*SB3)
 311      $     -1.0
 312       A4=0.2890E+00-0.4311E+00*SB+0.7387E+00*SB2-.3697E+00*SB3
 313       A5=0.0000E+00+0.4320E+00*SB+0.2449E+00*SB2-0.6670E-01*SB3
 314       GOTO 100
 315 C   ifl =     1
 316    32 A0=0.2780E+00*(1.0 + 0.4355E+00*SB-0.4584E+00*SB2+0.4390E-01*SB3)
 317       A1=0.2760E+00+0.1420E-01*SB+0.1480E-01*SB2-.9800E-02*SB3
 318       A2=0.3710E+01+0.8250E+00*SB-.3581E+00*SB2+0.1978E+00*SB3
 319       A3=0.2928E+02*(1.0 -.10154E+01*SB+0.6037E+00*SB2-.1175E+00*SB3)
 320      $     -1.0
 321       A4=0.8070E+00-.3575E+00*SB+0.4920E-01*SB2+0.1584E-01*SB3
 322       A5=0.0000E+00+0.1860E-01*SB+0.2080E+00*SB2-.450E-01*SB3
 323       GOTO 100
 324 C   ifl =     0
 325    33 A0=0.2924E+01*(1.0 -.18916E+01*SB+0.1191E+01*SB2-.2492E+00*SB3)
 326       A1=0.0000E+00-.9167E+00*SB+0.11147E+01*SB2-.3329E+00*SB3
 327       A2=0.8529E+01+0.7080E+00*SB-.11345E+01*SB2-.10563E+01*SB3
 328       A3=0.1420E+03*(1.0 -.15346E+01*SB+0.7261E+00*SB2-.5730E-01*SB3)
 329      $     -1.0
 330       A4=0.3396E+01-.11541E+01*SB-.8834E+00*SB2+0.2430E+00*SB3
 331       A5=0.0000E+00+0.1645E+00*SB+0.19041E+01*SB2+0.1474E+00*SB3
 332       GOTO 100
 333 C   ifl =    -1
 334    34 A0=0.3471E+00*(1.0- 0.1753E+00*SB-.9189E+00*SB2+0.6211E+00*SB3)
 335       A1=0.1900E-01-.4579E+00*SB+0.2112E+00*SB2-.6180E-01*SB3
 336       A2=0.7301E+01-.17308E+01*SB+.13666E+01*SB2-.6400E-02*SB3
 337       A3=0.1853E-04*(1.0 -.18260E+02*SB-.2872E+02*SB2-.23456E+02*SB3)
 338      $     -1.0
 339       A4=0.4400E+00-.4672E+00*SB+0.6532E+00*SB2-.3222E+00*SB3
 340       A5=0.0000E+00-.4679E+00*SB+0.10741E+01*SB2-.5663E+00*SB3
 341       GOTO 100
 342 C   ifl =    -2
 343    35 A0=0.1702E+00*(1.0 -.1041E+01*SB+0.4064E+00*SB2-.5888E-01*SB3)
 344       A1=-.9300E-01-.4742E-01*SB-.1959E+00*SB2+0.1039E+00*SB3
 345       A2=0.9119E+01-.7331E-01*SB+0.3506E+00*SB2-.2081E+00*SB3
 346       A3=0.2981E+02*(1.0 -.1912E+00*SB-.8947E-02*SB2+0.8805E-02*SB3)
 347      $     -1.0
 348       A4=0.1668E+01-.6678E-02*SB-.2894E+00*SB2+0.1221E+00*SB3
 349       A5=0.0000E+00+0.1245E+01*SB-.7843E+00*SB2+0.3724E+00*SB3
 350       GOTO 100
 351 C   ifl =    -3
 352    36 A0=0.3910E+00*(1.0 -.1103E+01*SB+0.5383E+00*SB2-.1083E+00*SB3)
 353       A1=-.1400E-01-.2471E+00*SB-.8042E-01*SB2+0.7193E-01*SB3
 354       A2=0.9812E+01-.4860E+01*SB+0.5958E+01*SB2-.2342E+01*SB3
 355       A3=0.3749E+00*(1.0 -.3569E+01*SB+0.5456E+01*SB2-.2344E+01*SB3)
 356      $     -1.0
 357       A4=0.4940E+00+0.2772E+00*SB-.2732E+00*SB2+0.6466E-01*SB3
 358       A5=0.0000E+00+0.3927E+00*SB-.3216E+00*SB2+0.2164E+00*SB3
 359       GOTO 100
 360 C   ifl =    -4
 361    37 A0=0.3815E-02*(1.0 + 0.2039E+02*SB-.2834E+02*SB2+0.1070E+02*SB3)
 362      $  * SQRT(STA - STB)
 363       A1=-.2789E-01-.7345E-03*SB-.3251E+00*SB2+0.1946E+00*SB3
 364       A2=0.3223E+01-.4268E+00*SB+0.4387E+01*SB2-.2401E+01*SB3
 365       A3=0.3338E-01*(1.0 -.1163E+02*SB+0.2995E+02*SB2-.1471E+02*SB3)
 366      $     -1.0
 367       A4=0.3646E+00-.5767E+00*SB+0.6088E+00*SB2-.2514E+00*SB3
 368       A5=0.1200E+01+0.2178E+00*SB-.4230E+00*SB2+0.4739E+00*SB3
 369       GOTO 100
 370 C   ifl =    -5
 371    38 A0=0.1666E-02*(1.0 + 0.9518E+01*SB-.4715E+01*SB2-.1060E+01*SB3)
 372      $  * SQRT(STA - STB)
 373       A1=-.1231E+00+0.1656E+00*SB-.5219E+00*SB2+0.2750E+00*SB3
 374       A2=0.3693E+01+0.4922E+01*SB-.1200E+02*SB2+0.7929E+01*SB3
 375       A3=0.1778E+00*(1.0 + 0.3036E+01*SB-.1184E+02*SB2+0.7940E+01*SB3)
 376      $     -1.0
 377       A4=0.5353E+00-.1401E+01*SB+0.1970E+01*SB2-.9405E+00*SB3
 378       A5=0.1590E+01+0.1025E+01*SB-.2318E+01*SB2+0.1380E+01*SB3
 379       GOTO 100
 380 C   ifl =    -6
 381    39 A0=0.4319E-03*(1.0 + 0.1100E+02*SB-.9520E+00*SB2+0.1434E+02*SB3)
 382      $  * SQRT(STA - STB)
 383       A1=-.2512E+00+0.3554E+00*SB-.4120E+00*SB2+0.1328E+00*SB3
 384       A2=0.4764E+01-.3513E+00*SB+0.1199E+02*SB2-.8290E+01*SB3
 385       A3=0.8458E-01*(1.0 + 0.2618E+01*SB+0.4407E+01*SB2+0.2991E+02*SB3)
 386      $    -1.0
 387       A4=0.3991E+00-.1363E+01*SB+0.1526E+01*SB2-.3179E+01*SB3
 388       A5=0.1981E+01+0.1496E+01*SB-.1501E+01*SB2+0.3880E+01*SB3
 389       GOTO 100
 390
 391     4 GOTO(41,42,43,44,45,46,47,48,49)IFLV
 392 C                                                             CTEQ1D
 393 C   ifl =     2
 394    41 A0=0.1634E+01*(1.0 -.8336E+00*SB+0.1640E+00*SB2+0.1530E+00*SB3)
 395       A1=0.5790E+00+0.8587E-01*SB-.6087E-01*SB2+0.1361E-01*SB3
 396       A2=0.2839E+01+0.3720E+00*SB+0.5264E+00*SB2+0.3538E-01*SB3
 397       A3=0.1095E+00*(1.0 -.4830E+00*SB+0.3708E+01*SB2-.6165E+00*SB3)
 398      $  -1.0
 399       A4=0.8010E+00-.1432E+00*SB+0.1442E+01*SB2-.1286E+01*SB3
 400       A5=0.0000E+00+0.1035E+01*SB-.5910E-01*SB2-.1982E+00*SB3
 401       GOTO 100
 402 C   ifl =     1
 403    42 A0=0.3535E+00*(1.0 + 0.4352E+00*SB-.2095E+00*SB2-.8455E-02*SB3)
 404       A1=0.2660E+00-.4096E-03*SB+0.1502E-01*SB2-.1163E-01*SB3
 405       A2=0.3514E+01+0.8219E+00*SB-.2330E+00*SB2+0.1055E+00*SB3
 406       A3=0.2200E+02*(1.0 -.9716E+00*SB+0.4552E+00*SB2-.8202E-01*SB3)
 407      $  -1.0
 408       A4=0.9000E+00-.3207E+00*SB-.4808E-01*SB2+0.3492E-01*SB3
 409       A5=0.0000E+00-.6273E-01*SB+0.1497E+00*SB2-.5683E-01*SB3
 410       GOTO 100
 411 C   ifl =     0
 412    43 A0=0.2743E+01*(1.0 -.2027E+01*SB+0.1517E+01*SB2-.4145E+00*SB3)
 413       A1=0.7000E-02-.9431E+00*SB+0.1231E+01*SB2-.4834E+00*SB3
 414       A2=0.8200E+01+0.1827E+01*SB-.3453E+01*SB2+0.6763E+00*SB3
 415       A3=0.4975E+02*(1.0 -.2329E+00*SB-.1245E+01*SB2+0.7194E+00*SB3)
 416      $  -1.0
 417       A4=0.2387E+01-.4077E+00*SB-.5542E+00*SB2-.9677E-02*SB3
 418       A5=0.0000E+00+0.2702E+00*SB+0.2389E+01*SB2-.8274E+00*SB3
 419       GOTO 100
 420 C   ifl =    -1
 421    44 A0=0.2015E+00*(1.0 -.2133E+00*SB-.6770E+00*SB2+0.5011E+00*SB3)
 422       A1=-.7700E-01-.7104E-01*SB-.3720E+00*SB2+0.2159E+00*SB3
 423       A2=0.8008E+01-.2049E+01*SB+0.1800E+01*SB2-.4660E+00*SB3
 424       A3=0.2923E-05*(1.0 + 0.2327E+02*SB+0.1500E+02*SB2+0.2633E+02*SB3)
 425      $  -1.0
 426       A4=0.9020E+00-.9191E+00*SB+0.1104E+01*SB2-.5863E+00*SB3
 427       A5=0.0000E+00+0.5840E+00*SB-.8720E+00*SB2+0.4234E+00*SB3
 428       GOTO 100
 429 C   ifl =    -2
 430    45 A0=0.9117E-01*(1.0 -.4089E+00*SB-.4361E+00*SB2+0.2512E+00*SB3)
 431       A1=-.2370E+00+0.2492E+00*SB-.3267E+00*SB2+0.1055E+00*SB3
 432       A2=0.8447E+01+0.6009E+00*SB+0.1003E+01*SB2-.1287E+01*SB3
 433       A3=0.3106E+02*(1.0 -.3901E-01*SB+0.1443E+00*SB2-.3433E+00*SB3)
 434      $  -1.0
 435       A4=0.1629E+01+0.7855E-01*SB-.1573E+00*SB2-.8595E-01*SB3
 436       A5=0.0000E+00+0.1558E+01*SB-.6295E+00*SB2+0.1847E+00*SB3
 437       GOTO 100
 438 C   ifl =    -3
 439    46 A0=0.3997E+00*(1.0 -.1046E+01*SB+0.6194E+00*SB2-.1342E+00*SB3)
 440       A1=0.2000E-02-.2544E+00*SB-.1958E+00*SB2+0.1458E+00*SB3
 441       A2=0.9613E+01-.3919E+01*SB+0.9573E+01*SB2-.5623E+01*SB3
 442       A3=0.3620E+00*(1.0 -.1858E+01*SB+0.8312E+01*SB2-.5900E+01*SB3)
 443      $  -1.0
 444       A4=0.3840E+00+0.3572E+00*SB-.1191E+01*SB2+0.7310E+00*SB3
 445       A5=0.0000E+00+0.3351E+00*SB-.7709E+00*SB2+0.4296E+00*SB3
 446       GOTO 100
 447 C   ifl =    -4
 448    47 A0=0.2156E-03*(1.0 + 0.2879E+02*SB-.2310E+02*SB2+0.9812E+01*SB3)
 449      $  * SQRT(STA - STB)
 450       A1=0.9086E-01-.1250E+00*SB-.7373E-01*SB2-.2201E-01*SB3
 451       A2=0.3588E+01+0.4518E+01*SB-.8930E-01*SB2+0.9163E-02*SB3
 452       A3=0.5216E+01*(1.0 + 0.5912E+00*SB-.4111E+00*SB2+0.7330E+00*SB3)
 453      $  -1.0
 454       A4=0.3145E+00+0.1233E+01*SB-.7478E+00*SB2+0.4657E+00*SB3
 455       A5=0.2723E+01-.4110E+00*SB+0.4868E-01*SB2-.3075E+00*SB3
 456       GOTO 100
 457 C   ifl =    -5
 458    48 A0=0.7476E-03*(1.0 + 0.1454E+02*SB-.2509E+02*SB2+0.1184E+02*SB3)
 459      $  * SQRT(STA - STB)
 460       A1=-.1955E-01-.1712E+00*SB-.1686E+00*SB2+0.2339E+00*SB3
 461       A2=0.4616E+01-.6859E+00*SB-.3959E+01*SB2+0.5530E+01*SB3
 462       A3=0.9881E+01*(1.0 -.1239E+02*SB+0.2721E+02*SB2-.1850E+02*SB3)
 463      $  -1.0
 464       A4=0.1200E+02-.1133E+02*SB+0.8138E+01*SB2+0.1199E+02*SB3
 465       A5=0.2226E+01-.5738E+00*SB+0.5239E+00*SB2+0.3825E+00*SB3
 466       GOTO 100
 467 C   ifl =    -6
 468    49 A0=0.8392E-06*(1.0 + 0.1844E+02*SB-.1110E+02*SB2-.2504E+02*SB3)
 469      $  * SQRT(STA - STB)
 470       A1=0.2127E+00-.5602E+00*SB+0.4777E+01*SB2-.1014E+02*SB3
 471       A2=0.1229E+01+0.7495E+01*SB-.5024E+01*SB2-.1200E+02*SB3
 472       A3=0.2868E+02*(1.0 + 0.7634E+01*SB-.2916E+02*SB2+0.2953E+02*SB3)
 473      $  -1.0
 474       A4=0.5970E+00+0.1138E+01*SB-.1439E+01*SB2-.1966E+01*SB3
 475       A5=0.6429E+01-.6673E+00*SB+0.7008E+01*SB2-.1157E+02*SB3
 476       GOTO 100
 477
 478     5 GOTO(51,52,53,54,55,56,57,58,59)IFLV
 479 C                                                             CTEQ1L
 480 C   ifl =     2
 481    51 A0=  1.791*(1.0 -0.449*SB-0.445*SB2+  0.401*SB3)
 482       A1=  0.608+  0.069*SB+  0.005*SB2-0.037*SB3
 483       A2=  3.470-0.375*SB+  2.267*SB2-1.261*SB3
 484       A3=  0.315*(1.0 -2.628*SB+  6.481*SB2-3.834*SB3)-1.0
 485       A4=  1.007-0.732*SB+  1.490*SB2-0.966*SB3
 486       A5=  0.000+  0.741*SB+  0.563*SB2-0.525*SB3
 487       GOTO 100
 488 C   ifl =     1
 489    52 A0=  0.513*(1.0 +   0.032*SB-0.120*SB2+  0.013*SB3)
 490       A1=  0.276+  0.052*SB+  0.000*SB2-0.006*SB3
 491       A2=  3.579+  0.763*SB-0.135*SB2+  0.083*SB3
 492       A3= 17.993*(1.0 -0.725*SB+  0.241*SB2-0.020*SB3)-1.0
 493       A4=  1.120-0.357*SB+  0.008*SB2+  0.028*SB3
 494       A5=  0.000+  0.311*SB+  0.029*SB2-0.010*SB3
 495       GOTO 100
 496 C   ifl =     0
 497    53 A0=  2.710*(1.0 -1.773*SB+  0.970*SB2-0.149*SB3)
 498       A1= -0.010-1.636*SB+  2.087*SB2-0.637*SB3
 499       A2=  7.174+  2.102*SB-2.209*SB2-0.420*SB3
 500       A3= 29.904*(1.0 -0.756*SB-0.506*SB2+  0.605*SB3)-1.0
 501       A4=  2.572-0.437*SB-0.968*SB2+  0.243*SB3
 502       A5=  0.000-1.776*SB+  4.266*SB2-0.335*SB3
 503       GOTO 100
 504 C   ifl =    -1
 505    54 A0=  0.278*(1.0 - 1.022*SB+  0.6457*SB2-0.1824*SB3)
 506       A1=  0.0862*SB-0.8657*SB2+  0.4185*SB3
 507       A2= 11.000-1.2809*SB+ 1.2516*SB2+0.061*SB3
 508       A3= 37.338*(1.0 - 0.9404*SB+  0.2517*SB2+0.1364*SB3)-1.0
 509       A4=  1.960-  0.3385*SB-0.3422*SB2+0.3653*SB3
 510       A5=  0.000+1.424*SB-2.7503*SB2+  1.2226*SB3
 511       GOTO 100
 512 C   ifl =    -2
 513    55 A0=  0.154*(1.0 -0.659*SB+  0.005*SB2+  0.061*SB3)
 514       A1= -0.128+  0.279*SB-0.786*SB2+  0.363*SB3
 515       A2=  8.649+  0.071*SB+  0.351*SB2-0.051*SB3
 516       A3= 43.685*(1.0 -0.603*SB+  0.037*SB2+  0.134*SB3)-1.0
 517       A4=  2.238-0.338*SB-0.199*SB2+  0.157*SB3
 518       A5=  0.000+  1.681*SB-2.068*SB2+  0.975*SB3
 519       GOTO 100
 520 C   ifl =    -3
 521    56 A0=  0.372*(1.0 -1.939*SB+  1.504*SB2-0.440*SB3)
 522       A1=  0.009+  0.610*SB-1.387*SB2+  0.579*SB3
 523       A2= 10.273-4.833*SB+  6.583*SB2-2.633*SB3
 524       A3=  0.160*(1.0 +  10.325*SB-2.027*SB2+  1.571*SB3)-1.0
 525       A4=  0.819-1.660*SB+  1.845*SB2-0.829*SB3
 526       A5=  0.000+  3.558*SB-3.940*SB2+  1.302*SB3
 527       GOTO 100
 528 C   ifl =    -4
 529    57 A0=  (7.5242E-5)*(1.0+22.0905*SB+7.1209*SB2-8.303*SB3)*
 530      $     SQRT(STA - STB)
 531       A1=  0.125-0.3027*SB+0.1564*SB2-0.091*SB3
 532       A2=  2.0388+1.2161*SB+11.5296*SB2-8.0659*SB3
 533       A3=  14.849*(1.0 -2.556*SB+3.5268*SB2-1.6353*SB3)-1.0
 534       A4=  0.3061-0.0901*SB+0.953*SB2-0.4871*SB3
 535       A5=  2.7352+0.1811*SB-0.5167*SB2+0.0543*SB3
 536       GOTO 100
 537 C   ifl =    -5
 538    58 A0=  (3.751E-4)*(1.0 + 21.5993*SB+3.1379*SB2-18.8328*SB3)*
 539      $     SQRT(STA - STB)
 540       A1= -0.0256-0.7717*SB+ 1.1499*SB2-0.5037*SB3
 541       A2=  4.9241+4.0107*SB-4.7012*SB2+0.1097*SB3
 542       A3=  2.842*(1.0 -2.2184*SB+  2.0293*SB2-0.6907*SB3)-1.0
 543       A4=  -0.1352+ 0.8753*SB-1.2626*SB2+  0.667*SB3
 544       A5=  1.5627-0.4917*SB+ 1.5927*SB2-0.351*SB3
 545       GOTO 100
 546 C   ifl =    -6
 547    59 A0=(2.725E-4)*(1.0 +  18.8497*SB-26.5797*SB2-29.0774*SB3)*
 548      $     SQRT(STA - STB)
 549       A1= -0.2204-1.0048*SB+0.9415*SB2-0.4274*SB3
 550       A2=  11.034-9.8362*SB-11.1034*SB2-9.1977*SB3
 551       A3=  2.084*(1.0 -2.881*SB+1.2778*SB2-2.9328*SB3)-1.0
 552       A4= -0.0872+  0.200*SB-1.6187*SB2-1.6058*SB3
 553       A5=  0.8684+4.7047*SB-1.4614*SB2-5.2309*SB3
 554       GOTO 100
 555
 556   100 CONTINUE
 557 C 100  CtqPdf = A0*(x**A1)*((1.-x)**A2)*(1.+A3*(x**A4))
 558 C     $     *(log(1.+1./x))**A5
 559       P012 = A0*(X**A1)*((1.-X)**A2)
 560       P34 = (1.+A3*(X**A4))
 561       P5 = (LOG(1.+1./X))**A5
 562 C
 563       PDFD = P012 * P34 * P5
 564       PDF = PDFD
 565 C
 566       Ctq1Pd = PDF
 567       IF(Ctq1Pd.LT.0.0) Ctq1Pd = 0.0
 568
 569       IST = ISET
 570
 571       LP  = IPRTN
 572       QSTO = QQ
 573
 574       RETURN
 575 C                                  -----------------------
 576       ENTRY Wlamd1 (ISET, IORDER)
 577
 578       IORDER = IORD (ISET)
 579       Wlamd1 = ALM  (ISET)
 580
 581         Ctq1Pd = 0.0
 582       RETURN
 583 C                                  -----------------------
 584       ENTRY PrCtq1
 585      >        (ISET, IORDR, ISCHEME, MXFLV,
 586      >        ALAM4, ALAM5, ALAM6, AMAS4, AMAS5, AMAS6,
 587      >        XMIN, QINI, QMAX, EXPNOR)
 588
 589 C                           Return QCD parameters and Fitting parameters
 590 C                           associated with parton distribution set Iset.
 591 C    Iord    : Order Of Fit
 592 C    Ischeme : (0, 1, 2)  for  (LO, MS-bar-NLO, DIS-NLO) resp.
 593 C    MxFlv   : Maximum number of flavors included
 594 C    Alam_i  : i = 4,5,6  Effective lambda for i-flavors
 595
 596 C    Amas_i  : i = 4,5,6  Mass parameter for flavor i
 597 C    Xmin, Qini, Qmax : self explanary
 598 C    ExpNor(I) : Normalization factor for the experimental data set used in
 599 C                obtaining the best global fit for parton distributions Iset:
 600 C     I = 1,     2,      3,     4,     5,     6,     7,     8
 601 C      BCDMS   NMC90  NMC280  CCFR   E605    WA70   E706   UA6
 602
 603       IORDR  = IORD (ISET)
 604       ISCHEME= ISCH (ISET)
 605       MXFLV  = NQRK (ISET)
 606
 607       ALAM4  = VLM(4,ISET)
 608       ALAM5  = VLM(5,ISET)
 609       ALAM6  = VLM(6,ISET)
 610
 611       AMAS4  = QMS(4,ISET)
 612       AMAS5  = QMS(5,ISET)
 613       AMAS6  = QMS(6,ISET)
 614
 615       XMIN   = XMN  (ISET)
 616       QINI   = QMN  (ISET)
 617       QMAX   = QMX  (ISET)
 618
 619       DO 101 IEXP = 1, NEXP(ISET)
 620         EXPNOR(IEXP) = EXPN(IEXP, ISET)
 621   101 CONTINUE
 622
 623         Ctq1Pd = 0.0
 624       RETURN
 625 C                         *************************
 626       END