[u/mrichter/AliRoot.git] / PDF / spdf / pdzxmt.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.2  1996/10/30 08:29:05  cernlib
* Version 7.04
*
* Revision 1.2  1996/10/28 15:04:20  pdflib
* PDFLIB new version 7.04
*
* Revision 1.1.1.1  1996/04/12 15:29:26  plothow
* Version 7.01
*
*
#include "pdf/pilot.h"
      FUNCTION PDZXMT (IST, LP, XX, QQ, IRT)
C
#include "pdf/w50511.inc"
#include "pdf/expdp.inc"
     +       AIEX, FXIEX, PDFD, ONED
      REAL PDF
C
      PARAMETER (D0=0.0, D1=1.0, D2=2.0, D3=3.0, D4=4.0, D10=10.0)
      PARAMETER (NEX = 3, MXFL = 6, NPN = 2, NST = 11)
C
      DIMENSION
     1 AC(0:NEX, 0:NPN, -MXFL:2, NST), A(0:NEX), T(0:NPN), FX(0:NEX),
     1 ALM(NST), Q02(NST), MEX(NST), MPN(NST), MQRK(NST), IORD(NST)
C
      DATA ONED/1.D0/
C
      DATA MEX, MPN, MQRK / NST*3, NST*2, NST*6 /
C                                          Set S1:    PDF parameters fro
      DATA IORD(1)/2/, ALM(1)/0.212/, Q02(1)/4.00/
      DATA (((AC(IEX,IPN,IFL,1), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
     > /  1.34,  0.15,
     >    5.30, -1.96, -0.57,  0.16,  0.43,  1.08, -0.08, -0.02,
     >    0.06, -0.03,  1.62,  0.11,  3.68, -1.94, -0.33,  0.14,
     >    0.53,  0.87, -0.10, -0.01,  0.03,  0.02,  1.88, -0.33,
     >    7.52, -1.34, -2.78,  0.10, -1.13,  2.92,  0.13, -0.04,
     >    0.04, -0.49, -0.99, -0.33,  8.53, -1.55, -1.54,  0.03,
     >   -1.08,  2.02,  0.10, -0.03,  0.39, -0.39, -0.99, -0.33,
     >    8.53, -1.55, -1.54,  0.03, -1.08,  2.02,  0.10, -0.03,
     >    0.39, -0.39, -0.99, -0.33,  8.53, -1.55, -1.54,  0.03,
     >   -1.08,  2.02,  0.10, -0.03,  0.39, -0.39, -3.98, -0.15,
     >    7.46,  0.35,  0.72, -0.06,  0.96,  0.89, -0.63,  0.00,
     >   -0.30, -0.04, -6.28, -0.18,  6.56,  0.65,  2.62,  0.02,
     >    1.40,  1.13, -1.18, -0.03, -0.38, -0.16,-13.08, -0.40,
     >   15.35, -0.43,  8.54,  0.31,-11.83,  3.18, -2.70, -0.12,
     >    4.16, -0.82 /
C                                           Set B1:    PDF parameters fr
      DATA IORD(2)/2/, ALM(2)/0.194/, Q02(2)/4.00/
      DATA (((AC(IEX,IPN,IFL,2), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
     > /  1.30,  0.19,
     >    5.24, -1.81, -0.57,  0.15,  0.44,  1.06, -0.09, -0.02,
     >    0.05, -0.02,  1.59,  0.14,  3.65, -1.81, -0.34,  0.13,
     >    0.53,  0.86, -0.10, -0.01,  0.03,  0.02,  1.48, -0.14,
     >    6.75, -0.50, -2.49, -0.11, -0.54,  2.13,  0.04,  0.03,
     >   -0.15, -0.24, -1.08, -0.13,  8.40, -0.88, -1.33, -0.21,
     >   -0.51,  1.18, -0.03,  0.06,  0.07, -0.05, -1.08, -0.13,
     >    8.39, -0.88, -1.33, -0.21, -0.50,  1.18, -0.03,  0.06,
     >    0.07, -0.05, -1.08, -0.13,  8.39, -0.88, -1.33, -0.21,
     >   -0.50,  1.18, -0.03,  0.06,  0.07, -0.05, -4.22, -0.02,
     >    7.29,  0.90,  0.88, -0.17,  1.08,  0.50, -0.69,  0.03,
     >   -0.39,  0.08, -6.42, -0.09,  6.47,  1.03,  2.67, -0.03,
     >    1.39,  1.00, -1.21, -0.02, -0.42, -0.14,-12.92, -0.36,
     >   15.74, -0.30,  8.33,  0.32,-12.73,  3.35, -2.68, -0.13,
     >    4.51, -0.91 /
C                                             Set B2:  PDF parameters fr
      DATA IORD(3)/2/, ALM(3)/0.191/, Q02(3)/4.00/
      DATA (((AC(IEX,IPN,IFL,3), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
     > /  1.38,  0.18,
     >    5.40, -1.91, -0.59,  0.16,  0.42,  1.11, -0.08, -0.02,
     >    0.06, -0.03,  1.64,  0.09,  3.74, -2.02, -0.33,  0.14,
     >    0.54,  0.88, -0.10, -0.01,  0.03,  0.02,  1.52, -0.72,
     >    7.75, -2.18, -2.71,  0.45, -1.56,  3.75,  0.15, -0.15,
     >    0.16, -0.76, -0.85, -0.82,  9.19, -2.76, -1.43,  0.35,
     >   -0.92,  2.56, -0.03, -0.09,  0.12, -0.40, -0.85, -0.82,
     >    9.19, -2.76, -1.43,  0.35, -0.92,  2.56, -0.03, -0.10,
     >    0.12, -0.40, -0.85, -0.82,  9.19, -2.76, -1.43,  0.35,
     >   -0.92,  2.56, -0.03, -0.10,  0.12, -0.40, -3.74, -0.58,
     >    9.63, -1.09,  0.21,  0.24, -1.13,  2.10, -0.50, -0.07,
     >    0.25, -0.33, -6.07, -0.52,  8.33, -0.52,  2.33,  0.22,
     >    0.28,  1.91, -1.15, -0.07, -0.28, -0.31,-12.08, -0.73,
     >   21.14, -1.92,  7.31,  0.54,-19.17,  4.59, -2.35, -0.18,
     >    6.64, -1.25 /
C                                             Set E1:  PDF parameters fr
      DATA IORD(4)/2/, ALM(4)/0.155/, Q02(4)/4.00/
      DATA (((AC(IEX,IPN,IFL,4), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
     > /  1.43,  0.16,
     >    6.17, -1.94, -0.65,  0.16,  0.43,  1.12, -0.08, -0.02,
     >    0.06, -0.02,  1.69,  0.11,  3.69, -1.99, -0.33,  0.14,
     >    0.54,  0.90, -0.11, -0.01,  0.03,  0.02,  2.11, -0.33,
     >    7.93, -1.51, -3.01,  0.10, -1.40,  3.14,  0.18, -0.04,
     >    0.09, -0.55, -0.84, -0.32,  8.96, -1.70, -1.65,  0.02,
     >   -1.24,  2.15,  0.12, -0.03,  0.45, -0.43, -0.84, -0.32,
     >    8.96, -1.70, -1.65,  0.02, -1.24,  2.15,  0.12, -0.03,
     >    0.45, -0.43, -0.84, -0.32,  8.96, -1.70, -1.65,  0.02,
     >   -1.24,  2.15,  0.12, -0.03,  0.45, -0.43, -3.87, -0.15,
     >    7.83,  0.21,  0.85, -0.07,  1.00,  0.93, -0.73,  0.00,
     >   -0.36, -0.03, -6.09, -0.17,  6.75,  0.54,  2.81,  0.01,
     >    1.74,  1.15, -1.34, -0.03, -0.56, -0.16,-12.56, -0.38,
     >   14.62, -0.41,  8.69,  0.30,-11.27,  3.19, -2.93, -0.12,
     >    4.29, -0.87 /
C                                   1/2 strange sea;  PDF parameters fro
      DATA IORD(5)/2/, ALM(5)/0.237/, Q02(5)/4.00/
      DATA (((AC(IEX,IPN,IFL,5), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
     > /  1.42,  0.16,
     >    5.40, -1.99, -0.59,  0.17,  0.41,  1.12, -0.08, -0.02,
     >    0.06, -0.03,  1.68,  0.08,  3.75, -2.09, -0.33,  0.15,
     >    0.53,  0.89, -0.10, -0.01,  0.03,  0.02,  0.90, -0.17,
     >    5.27, -0.20, -1.86, -0.10,  0.43,  1.67, -0.09,  0.02,
     >   -0.26, -0.14, -1.48, -0.13,  7.83, -0.38, -0.89, -0.19,
     >   -0.06,  0.68, -0.12,  0.04,  0.01,  0.05, -1.48, -0.13,
     >    7.83, -0.38, -0.89, -0.19, -0.05,  0.68, -0.13,  0.04,
     >    0.00,  0.05, -2.26, -0.15,  7.47, -0.23, -0.90, -0.10,
     >   -0.61,  1.22, -0.06,  0.01,  0.28, -0.16, -4.68, -0.06,
     >    5.55,  1.13,  0.92, -0.12,  1.16,  0.50, -0.62,  0.01,
     >   -0.26,  0.03, -6.83, -0.12,  5.24,  1.19,  2.68, -0.01,
     >    1.14,  0.93, -1.13, -0.03, -0.24, -0.13,-14.41, -0.28,
     >   11.48,  0.65,  9.65,  0.15, -7.50,  1.99, -2.98, -0.06,
     >    2.54, -0.43 /
C                                      Set S1-MSbar:  PDF parameters fro
      DATA IORD(6)/2/, ALM(6)/0.212/, Q02(6)/4.0/
      DATA (((AC(IEX,IPN,IFL,6), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
     > /  1.75,  0.11,
     >    6.20, -2.35, -1.02,  0.26, -0.41,  1.68,  0.05, -0.06,
     >    0.29, -0.24,  2.03,  0.06,  4.43, -2.35, -0.78,  0.24,
     >   -0.18,  1.52,  0.03, -0.04,  0.22, -0.19,  1.09, -0.24,
     >    5.97, -0.64, -2.41,  0.08, -0.90,  2.71, -0.12,  0.02,
     >   -0.35, -0.20, -0.14, -0.49, 10.24, -2.57, -1.98,  0.02,
     >   -1.43,  2.32,  0.23, -0.02,  0.44, -0.47, -0.14, -0.49,
     >   10.24, -2.57, -1.98,  0.02, -1.44,  2.32,  0.23, -0.02,
     >    0.45, -0.47, -0.15, -0.49, 10.23, -2.57, -1.98,  0.02,
     >   -1.44,  2.32,  0.23, -0.02,  0.45, -0.47, -2.36, -0.49,
     >    9.00, -1.74, -1.42,  0.44, -0.46,  3.93,  0.21, -0.22,
     >    0.29, -1.34, -2.19, -1.07, 11.30, -4.85, -3.86,  1.56,
     >   -7.20, 10.51,  1.57, -0.73,  3.85, -4.36,-24.77,  7.52,
     >  -99.51, 36.02,-23.00,  0.48,-16.45, 16.51, 34.44, -6.26,
     >   97.19,-40.40 /
C                                      Set B1-MSbar:  PDF parameters fro
      DATA IORD(7)/2/, ALM(7)/0.194/, Q02(7)/4.0/
      DATA (((AC(IEX,IPN,IFL,7), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
     > /  1.66,  0.18,
     >    6.04, -2.09, -0.94,  0.18, -0.25,  1.42,  0.03, -0.03,
     >    0.23, -0.14,  2.00,  0.09,  4.40, -2.24, -0.81,  0.24,
     >   -0.20,  1.53,  0.05, -0.05,  0.25, -0.23,  0.92, -0.07,
     >    5.79, -0.01, -2.28, -0.16, -0.68,  1.93, -0.07,  0.06,
     >   -0.23, -0.11, -0.60, -0.13,  9.31, -1.18, -1.76, -0.27,
     >   -0.94,  1.31,  0.13,  0.09,  0.21, -0.10, -0.60, -0.13,
     >    9.31, -1.18, -1.76, -0.27, -0.94,  1.31,  0.13,  0.09,
     >    0.21, -0.10, -0.60, -0.13,  9.31, -1.18, -1.76, -0.27,
     >   -0.94,  1.31,  0.14,  0.09,  0.21, -0.10, -2.94, -0.19,
     >    7.94, -0.46, -1.12,  0.16, -0.05,  2.93,  0.15, -0.13,
     >    0.27, -1.05, -2.95, -0.62,  9.97, -3.00, -3.21,  0.99,
     >   -6.33,  8.42,  1.38, -0.51,  3.71, -3.61, -3.88, -0.78,
     >    3.80, -2.37, -1.59, -0.07,  2.13,  0.48, -0.05,  0.40,
     >    0.96,  2.30 /
C                                      Set B2-MSbar:  PDF parameters fro
      DATA IORD(8)/2/, ALM(8)/0.191/, Q02(8)/4.0/
      DATA (((AC(IEX,IPN,IFL,8), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
     > /  1.77,  0.13,
     >    6.28, -2.30, -0.98,  0.23, -0.34,  1.60,  0.03, -0.04,
     >    0.26, -0.18,  2.04,  0.03,  4.48, -2.47, -0.75,  0.26,
     >   -0.15,  1.52,  0.02, -0.05,  0.21, -0.19,  0.74, -0.59,
     >    6.31, -1.37, -2.44,  0.42, -1.62,  3.56,  0.07, -0.15,
     >    0.18, -0.77, -0.43, -0.86, 10.16, -3.14, -1.96,  0.43,
     >   -1.91,  3.14,  0.20, -0.14,  0.53, -0.68, -0.43, -0.86,
     >   10.16, -3.14, -1.96,  0.43, -1.91,  3.14,  0.20, -0.14,
     >    0.53, -0.68, -0.43, -0.86, 10.16, -3.14, -1.96,  0.43,
     >   -1.91,  3.14,  0.20, -0.14,  0.53, -0.68, -3.07, -0.66,
     >    8.57, -1.68, -1.03,  0.45, -0.32,  3.48,  0.04, -0.17,
     >    0.17, -0.98, -4.44, -0.68,  6.90, -1.82, -0.13,  0.50,
     >    1.46,  4.11, -0.23, -0.18, -0.53, -1.16, -7.03, -1.13,
     >    8.56, -5.47,  1.10,  1.07, -2.33,  9.08, -0.41, -0.35,
     >    0.87, -2.66 /
C                                      Set E1-MSbar:  PDF parameters fro
      DATA IORD(9)/2/, ALM(9)/0.155/, Q02(9)/4.0/
      DATA (((AC(IEX,IPN,IFL,9), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
     > /  1.79,  0.12,
     >    7.03, -2.29, -1.05,  0.24, -0.38,  1.63,  0.03, -0.04,
     >    0.27, -0.18,  2.12,  0.02,  4.46, -2.57, -0.85,  0.32,
     >   -0.28,  1.82,  0.07, -0.08,  0.29, -0.33,  1.58, -0.28,
     >    6.84, -1.08, -2.68,  0.05, -0.93,  2.76,  0.01,  0.00,
     >   -0.26, -0.32, -0.10, -0.43, 10.43, -2.49, -2.29,  0.09,
     >   -2.14,  2.80,  0.35, -0.06,  0.73, -0.67, -0.10, -0.43,
     >   10.43, -2.49, -2.29,  0.09, -2.14,  2.80,  0.35, -0.06,
     >    0.73, -0.67, -0.11, -0.43, 10.43, -2.48, -2.29,  0.09,
     >   -2.14,  2.80,  0.35, -0.06,  0.73, -0.67, -2.53, -0.35,
     >    8.67, -1.24, -1.16,  0.26, -0.10,  3.26,  0.12, -0.15,
     >    0.27, -1.06, -3.91, -0.44,  6.85, -1.56, -0.19,  0.38,
     >    2.15,  4.07, -0.24, -0.17, -0.74, -1.24, -6.57, -0.90,
     >    7.27, -5.07,  1.15,  0.95, -0.28,  9.02, -0.48, -0.33,
     >    0.28, -2.75 /
C                         (MSbar)   1/2 strange sea;  PDF parameters fro
      DATA IORD(10)/2/, ALM(10)/0.237/, Q02(10)/4.00/
      DATA (((AC(IEX,IPN,IFL,10), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
     > /  1.84,  0.12,
     >    6.34, -2.40, -0.97,  0.22, -0.34,  1.53,  0.03, -0.04,
     >    0.25, -0.16,  2.08,  0.02,  4.53, -2.51, -0.66,  0.19,
     >   -0.04,  1.24, -0.02, -0.01,  0.15, -0.05,  0.31, -0.10,
     >    4.18,  0.34, -1.84, -0.10,  0.05,  1.64, -0.06,  0.01,
     >   -0.12, -0.16, -1.13, -0.15,  8.43, -0.64, -1.26, -0.16,
     >   -0.39,  1.01, -0.01,  0.03,  0.05, -0.06, -1.13, -0.15,
     >    8.43, -0.64, -1.26, -0.16, -0.39,  1.01, -0.01,  0.03,
     >    0.05, -0.06, -1.82, -0.18,  7.94, -0.56, -1.40, -0.06,
     >   -0.82,  1.65,  0.09, -0.01,  0.30, -0.31, -3.69, -0.15,
     >    5.72,  0.26, -0.47,  0.04,  0.93,  1.85, -0.10, -0.05,
     >   -0.11, -0.50, -5.06, -0.25,  4.42, -0.14,  0.39,  0.16,
     >    2.38,  2.72, -0.35, -0.08, -0.63, -0.75, -9.92, -0.38,
     >   -1.27, -1.60,  4.60,  0.24,  9.17,  4.40, -1.53, -0.08,
     >   -2.88, -1.08 /
C                                             Set B0:  PDF parameters fr
      DATA IORD(11)/1/, ALM(11)/0.144/, Q02(11)/4.00/
      DATA (((AC(IEX,IPN,IFL,11), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
     > /  1.38,  0.16,
     >    5.40, -1.97, -0.62,  0.19,  0.59,  1.24, -0.10, -0.02,
     >    0.03, -0.05,  1.67,  0.08,  3.75, -2.09, -0.33,  0.17,
     >    0.70,  0.98, -0.13, -0.01,  0.00,  0.02,  1.52, -0.25,
     >    7.01, -0.79, -3.17, -0.01, -0.90,  2.90,  0.25,  0.00,
     >   -0.08, -0.54, -0.81, -0.07,  9.19, -0.89, -1.13, -0.46,
     >    0.35,  0.33, -0.26,  0.16, -0.49,  0.40, -0.81, -0.07,
     >    9.19, -0.89, -1.13, -0.46,  0.35,  0.33, -0.26,  0.16,
     >   -0.49,  0.40, -0.81, -0.07,  9.19, -0.89, -1.13, -0.46,
     >    0.35,  0.33, -0.26,  0.16, -0.49,  0.40, -3.62, -0.06,
     >    8.30,  0.16,  0.03, -0.21, -0.60,  1.26, -0.48,  0.05,
     >    0.25, -0.15, -6.16, -0.11,  6.49,  0.71,  2.37, -0.05,
     >    1.28,  1.37, -1.24, -0.02, -0.41, -0.26,-12.68, -0.35,
     >   14.87, -0.17,  8.36,  0.28,-12.56,  3.39, -2.89, -0.12,
     >    4.75, -0.96 /

      IRT = 0
      IFL = LP
C
      PDZXMT = 0.
      IBLP = ABS(LP)
      IF(IBLP .GE. 4 .AND. QQ .LE. 1.5) RETURN
      IF(IBLP .GE. 5 .AND. QQ .LE. 5.0) RETURN
C
      TM = TMAS
      X  = XX
      Q0 = SQRT (Q02(IST))
      ALAM = ALM(IST)
      SQ = LOG ( LOG( MAX(Q0,QQ) /ALAM) / LOG(Q0/ALAM) )
      IF(IBLP .EQ. 6) THEN
        SQ = LOG ( LOG( MAX(Q0,QQ) /ALAM * (100./TM)) / LOG(Q0/ALAM))
      ENDIF
C
      FX(0) = EXP(D1)
      FX(1) = X
      FX(2) = 1.- X
      FX(3) = LOG (1.+ 1./X)

      PDFD = ONED
      DO 20 IEX = 0, MEX(IST)
        A(IEX) = AC(IEX, 0, IFL, IST)
        SQN=SQ
        DO 21 IPN = 1, MPN(IST)
          A(IEX) = A(IEX) + AC(IEX, IPN, IFL, IST) * SQN
          SQN    = SQN*SQ
21      CONTINUE
        FXIEX = FX(IEX)
        AIEX  = A(IEX)
        PDFD = PDFD * FXIEX **(AIEX)
   20 CONTINUE
C
C     PDZXMT = PDF / X
      PDF = PDFD
      IF(PDF.GT.0.) PDZXMT = PDF
C
      RETURN
C
      ENTRY VLAMBD (ISET, IORDER)
C
      IORDER = IORD (ISET)
      VLAMBD = ALM  (ISET)
C
      RETURN
C                         *************************
      END
Commit	Line	Data
21886bb6	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.2 1996/10/30 08:29:05 cernlib
	6	* Version 7.04
	7	*
	8	* Revision 1.2 1996/10/28 15:04:20 pdflib
	9	* PDFLIB new version 7.04
	10	*
	11	* Revision 1.1.1.1 1996/04/12 15:29:26 plothow
	12	* Version 7.01
	13	*
	14	*
	15	#include "pdf/pilot.h"
	16	FUNCTION PDZXMT (IST, LP, XX, QQ, IRT)
	17	C
	18	#include "pdf/w50511.inc"
	19	#include "pdf/expdp.inc"
	20	+ AIEX, FXIEX, PDFD, ONED
	21	REAL PDF
	22	C
	23	PARAMETER (D0=0.0, D1=1.0, D2=2.0, D3=3.0, D4=4.0, D10=10.0)
	24	PARAMETER (NEX = 3, MXFL = 6, NPN = 2, NST = 11)
	25	C
	26	DIMENSION
	27	1 AC(0:NEX, 0:NPN, -MXFL:2, NST), A(0:NEX), T(0:NPN), FX(0:NEX),
	28	1 ALM(NST), Q02(NST), MEX(NST), MPN(NST), MQRK(NST), IORD(NST)
	29	C
	30	DATA ONED/1.D0/
	31	C
	32	DATA MEX, MPN, MQRK / NST3, NST2, NST*6 /
	33	C Set S1: PDF parameters fro
	34	DATA IORD(1)/2/, ALM(1)/0.212/, Q02(1)/4.00/
	35	DATA (((AC(IEX,IPN,IFL,1), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
	36	> / 1.34, 0.15,
	37	> 5.30, -1.96, -0.57, 0.16, 0.43, 1.08, -0.08, -0.02,
	38	> 0.06, -0.03, 1.62, 0.11, 3.68, -1.94, -0.33, 0.14,
	39	> 0.53, 0.87, -0.10, -0.01, 0.03, 0.02, 1.88, -0.33,
	40	> 7.52, -1.34, -2.78, 0.10, -1.13, 2.92, 0.13, -0.04,
	41	> 0.04, -0.49, -0.99, -0.33, 8.53, -1.55, -1.54, 0.03,
	42	> -1.08, 2.02, 0.10, -0.03, 0.39, -0.39, -0.99, -0.33,
	43	> 8.53, -1.55, -1.54, 0.03, -1.08, 2.02, 0.10, -0.03,
	44	> 0.39, -0.39, -0.99, -0.33, 8.53, -1.55, -1.54, 0.03,
	45	> -1.08, 2.02, 0.10, -0.03, 0.39, -0.39, -3.98, -0.15,
	46	> 7.46, 0.35, 0.72, -0.06, 0.96, 0.89, -0.63, 0.00,
	47	> -0.30, -0.04, -6.28, -0.18, 6.56, 0.65, 2.62, 0.02,
	48	> 1.40, 1.13, -1.18, -0.03, -0.38, -0.16,-13.08, -0.40,
	49	> 15.35, -0.43, 8.54, 0.31,-11.83, 3.18, -2.70, -0.12,
	50	> 4.16, -0.82 /
	51	C Set B1: PDF parameters fr
	52	DATA IORD(2)/2/, ALM(2)/0.194/, Q02(2)/4.00/
	53	DATA (((AC(IEX,IPN,IFL,2), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
	54	> / 1.30, 0.19,
	55	> 5.24, -1.81, -0.57, 0.15, 0.44, 1.06, -0.09, -0.02,
	56	> 0.05, -0.02, 1.59, 0.14, 3.65, -1.81, -0.34, 0.13,
	57	> 0.53, 0.86, -0.10, -0.01, 0.03, 0.02, 1.48, -0.14,
	58	> 6.75, -0.50, -2.49, -0.11, -0.54, 2.13, 0.04, 0.03,
	59	> -0.15, -0.24, -1.08, -0.13, 8.40, -0.88, -1.33, -0.21,
	60	> -0.51, 1.18, -0.03, 0.06, 0.07, -0.05, -1.08, -0.13,
	61	> 8.39, -0.88, -1.33, -0.21, -0.50, 1.18, -0.03, 0.06,
	62	> 0.07, -0.05, -1.08, -0.13, 8.39, -0.88, -1.33, -0.21,
	63	> -0.50, 1.18, -0.03, 0.06, 0.07, -0.05, -4.22, -0.02,
	64	> 7.29, 0.90, 0.88, -0.17, 1.08, 0.50, -0.69, 0.03,
65	> -0.39, 0.08, -6.42, -0.09, 6.47, 1.03, 2.67, -0.03,
66	> 1.39, 1.00, -1.21, -0.02, -0.42, -0.14,-12.92, -0.36,
67	> 15.74, -0.30, 8.33, 0.32,-12.73, 3.35, -2.68, -0.13,
68	> 4.51, -0.91 /
69	C Set B2: PDF parameters fr
70	DATA IORD(3)/2/, ALM(3)/0.191/, Q02(3)/4.00/
71	DATA (((AC(IEX,IPN,IFL,3), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
72	> / 1.38, 0.18,
73	> 5.40, -1.91, -0.59, 0.16, 0.42, 1.11, -0.08, -0.02,
74	> 0.06, -0.03, 1.64, 0.09, 3.74, -2.02, -0.33, 0.14,
75	> 0.54, 0.88, -0.10, -0.01, 0.03, 0.02, 1.52, -0.72,
76	> 7.75, -2.18, -2.71, 0.45, -1.56, 3.75, 0.15, -0.15,
77	> 0.16, -0.76, -0.85, -0.82, 9.19, -2.76, -1.43, 0.35,
78	> -0.92, 2.56, -0.03, -0.09, 0.12, -0.40, -0.85, -0.82,
79	> 9.19, -2.76, -1.43, 0.35, -0.92, 2.56, -0.03, -0.10,
80	> 0.12, -0.40, -0.85, -0.82, 9.19, -2.76, -1.43, 0.35,
81	> -0.92, 2.56, -0.03, -0.10, 0.12, -0.40, -3.74, -0.58,
82	> 9.63, -1.09, 0.21, 0.24, -1.13, 2.10, -0.50, -0.07,
83	> 0.25, -0.33, -6.07, -0.52, 8.33, -0.52, 2.33, 0.22,
84	> 0.28, 1.91, -1.15, -0.07, -0.28, -0.31,-12.08, -0.73,
85	> 21.14, -1.92, 7.31, 0.54,-19.17, 4.59, -2.35, -0.18,
86	> 6.64, -1.25 /
87	C Set E1: PDF parameters fr
88	DATA IORD(4)/2/, ALM(4)/0.155/, Q02(4)/4.00/
89	DATA (((AC(IEX,IPN,IFL,4), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
90	> / 1.43, 0.16,
91	> 6.17, -1.94, -0.65, 0.16, 0.43, 1.12, -0.08, -0.02,
92	> 0.06, -0.02, 1.69, 0.11, 3.69, -1.99, -0.33, 0.14,
93	> 0.54, 0.90, -0.11, -0.01, 0.03, 0.02, 2.11, -0.33,
94	> 7.93, -1.51, -3.01, 0.10, -1.40, 3.14, 0.18, -0.04,
95	> 0.09, -0.55, -0.84, -0.32, 8.96, -1.70, -1.65, 0.02,
96	> -1.24, 2.15, 0.12, -0.03, 0.45, -0.43, -0.84, -0.32,
97	> 8.96, -1.70, -1.65, 0.02, -1.24, 2.15, 0.12, -0.03,
98	> 0.45, -0.43, -0.84, -0.32, 8.96, -1.70, -1.65, 0.02,
99	> -1.24, 2.15, 0.12, -0.03, 0.45, -0.43, -3.87, -0.15,
100	> 7.83, 0.21, 0.85, -0.07, 1.00, 0.93, -0.73, 0.00,
101	> -0.36, -0.03, -6.09, -0.17, 6.75, 0.54, 2.81, 0.01,
102	> 1.74, 1.15, -1.34, -0.03, -0.56, -0.16,-12.56, -0.38,
103	> 14.62, -0.41, 8.69, 0.30,-11.27, 3.19, -2.93, -0.12,
104	> 4.29, -0.87 /
105	C 1/2 strange sea; PDF parameters fro
106	DATA IORD(5)/2/, ALM(5)/0.237/, Q02(5)/4.00/
107	DATA (((AC(IEX,IPN,IFL,5), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
108	> / 1.42, 0.16,
109	> 5.40, -1.99, -0.59, 0.17, 0.41, 1.12, -0.08, -0.02,
110	> 0.06, -0.03, 1.68, 0.08, 3.75, -2.09, -0.33, 0.15,
111	> 0.53, 0.89, -0.10, -0.01, 0.03, 0.02, 0.90, -0.17,
112	> 5.27, -0.20, -1.86, -0.10, 0.43, 1.67, -0.09, 0.02,
113	> -0.26, -0.14, -1.48, -0.13, 7.83, -0.38, -0.89, -0.19,
114	> -0.06, 0.68, -0.12, 0.04, 0.01, 0.05, -1.48, -0.13,
115	> 7.83, -0.38, -0.89, -0.19, -0.05, 0.68, -0.13, 0.04,
116	> 0.00, 0.05, -2.26, -0.15, 7.47, -0.23, -0.90, -0.10,
117	> -0.61, 1.22, -0.06, 0.01, 0.28, -0.16, -4.68, -0.06,
118	> 5.55, 1.13, 0.92, -0.12, 1.16, 0.50, -0.62, 0.01,
119	> -0.26, 0.03, -6.83, -0.12, 5.24, 1.19, 2.68, -0.01,
120	> 1.14, 0.93, -1.13, -0.03, -0.24, -0.13,-14.41, -0.28,
121	> 11.48, 0.65, 9.65, 0.15, -7.50, 1.99, -2.98, -0.06,
122	> 2.54, -0.43 /
123	C Set S1-MSbar: PDF parameters fro
124	DATA IORD(6)/2/, ALM(6)/0.212/, Q02(6)/4.0/
125	DATA (((AC(IEX,IPN,IFL,6), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
126	> / 1.75, 0.11,
127	> 6.20, -2.35, -1.02, 0.26, -0.41, 1.68, 0.05, -0.06,
128	> 0.29, -0.24, 2.03, 0.06, 4.43, -2.35, -0.78, 0.24,
129	> -0.18, 1.52, 0.03, -0.04, 0.22, -0.19, 1.09, -0.24,
130	> 5.97, -0.64, -2.41, 0.08, -0.90, 2.71, -0.12, 0.02,
131	> -0.35, -0.20, -0.14, -0.49, 10.24, -2.57, -1.98, 0.02,
132	> -1.43, 2.32, 0.23, -0.02, 0.44, -0.47, -0.14, -0.49,
133	> 10.24, -2.57, -1.98, 0.02, -1.44, 2.32, 0.23, -0.02,
134	> 0.45, -0.47, -0.15, -0.49, 10.23, -2.57, -1.98, 0.02,
135	> -1.44, 2.32, 0.23, -0.02, 0.45, -0.47, -2.36, -0.49,
136	> 9.00, -1.74, -1.42, 0.44, -0.46, 3.93, 0.21, -0.22,
137	> 0.29, -1.34, -2.19, -1.07, 11.30, -4.85, -3.86, 1.56,
138	> -7.20, 10.51, 1.57, -0.73, 3.85, -4.36,-24.77, 7.52,
139	> -99.51, 36.02,-23.00, 0.48,-16.45, 16.51, 34.44, -6.26,
140	> 97.19,-40.40 /
141	C Set B1-MSbar: PDF parameters fro
142	DATA IORD(7)/2/, ALM(7)/0.194/, Q02(7)/4.0/
143	DATA (((AC(IEX,IPN,IFL,7), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
144	> / 1.66, 0.18,
145	> 6.04, -2.09, -0.94, 0.18, -0.25, 1.42, 0.03, -0.03,
146	> 0.23, -0.14, 2.00, 0.09, 4.40, -2.24, -0.81, 0.24,
147	> -0.20, 1.53, 0.05, -0.05, 0.25, -0.23, 0.92, -0.07,
148	> 5.79, -0.01, -2.28, -0.16, -0.68, 1.93, -0.07, 0.06,
149	> -0.23, -0.11, -0.60, -0.13, 9.31, -1.18, -1.76, -0.27,
150	> -0.94, 1.31, 0.13, 0.09, 0.21, -0.10, -0.60, -0.13,
151	> 9.31, -1.18, -1.76, -0.27, -0.94, 1.31, 0.13, 0.09,
152	> 0.21, -0.10, -0.60, -0.13, 9.31, -1.18, -1.76, -0.27,
153	> -0.94, 1.31, 0.14, 0.09, 0.21, -0.10, -2.94, -0.19,
154	> 7.94, -0.46, -1.12, 0.16, -0.05, 2.93, 0.15, -0.13,
155	> 0.27, -1.05, -2.95, -0.62, 9.97, -3.00, -3.21, 0.99,
156	> -6.33, 8.42, 1.38, -0.51, 3.71, -3.61, -3.88, -0.78,
157	> 3.80, -2.37, -1.59, -0.07, 2.13, 0.48, -0.05, 0.40,
158	> 0.96, 2.30 /
159	C Set B2-MSbar: PDF parameters fro
160	DATA IORD(8)/2/, ALM(8)/0.191/, Q02(8)/4.0/
161	DATA (((AC(IEX,IPN,IFL,8), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
162	> / 1.77, 0.13,
163	> 6.28, -2.30, -0.98, 0.23, -0.34, 1.60, 0.03, -0.04,
164	> 0.26, -0.18, 2.04, 0.03, 4.48, -2.47, -0.75, 0.26,
165	> -0.15, 1.52, 0.02, -0.05, 0.21, -0.19, 0.74, -0.59,
166	> 6.31, -1.37, -2.44, 0.42, -1.62, 3.56, 0.07, -0.15,
167	> 0.18, -0.77, -0.43, -0.86, 10.16, -3.14, -1.96, 0.43,
168	> -1.91, 3.14, 0.20, -0.14, 0.53, -0.68, -0.43, -0.86,
169	> 10.16, -3.14, -1.96, 0.43, -1.91, 3.14, 0.20, -0.14,
170	> 0.53, -0.68, -0.43, -0.86, 10.16, -3.14, -1.96, 0.43,
171	> -1.91, 3.14, 0.20, -0.14, 0.53, -0.68, -3.07, -0.66,
172	> 8.57, -1.68, -1.03, 0.45, -0.32, 3.48, 0.04, -0.17,
173	> 0.17, -0.98, -4.44, -0.68, 6.90, -1.82, -0.13, 0.50,
174	> 1.46, 4.11, -0.23, -0.18, -0.53, -1.16, -7.03, -1.13,
175	> 8.56, -5.47, 1.10, 1.07, -2.33, 9.08, -0.41, -0.35,
176	> 0.87, -2.66 /
177	C Set E1-MSbar: PDF parameters fro
178	DATA IORD(9)/2/, ALM(9)/0.155/, Q02(9)/4.0/
179	DATA (((AC(IEX,IPN,IFL,9), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
180	> / 1.79, 0.12,
181	> 7.03, -2.29, -1.05, 0.24, -0.38, 1.63, 0.03, -0.04,
182	> 0.27, -0.18, 2.12, 0.02, 4.46, -2.57, -0.85, 0.32,
183	> -0.28, 1.82, 0.07, -0.08, 0.29, -0.33, 1.58, -0.28,
184	> 6.84, -1.08, -2.68, 0.05, -0.93, 2.76, 0.01, 0.00,
185	> -0.26, -0.32, -0.10, -0.43, 10.43, -2.49, -2.29, 0.09,
186	> -2.14, 2.80, 0.35, -0.06, 0.73, -0.67, -0.10, -0.43,
187	> 10.43, -2.49, -2.29, 0.09, -2.14, 2.80, 0.35, -0.06,
188	> 0.73, -0.67, -0.11, -0.43, 10.43, -2.48, -2.29, 0.09,
189	> -2.14, 2.80, 0.35, -0.06, 0.73, -0.67, -2.53, -0.35,
190	> 8.67, -1.24, -1.16, 0.26, -0.10, 3.26, 0.12, -0.15,
191	> 0.27, -1.06, -3.91, -0.44, 6.85, -1.56, -0.19, 0.38,
192	> 2.15, 4.07, -0.24, -0.17, -0.74, -1.24, -6.57, -0.90,
193	> 7.27, -5.07, 1.15, 0.95, -0.28, 9.02, -0.48, -0.33,
194	> 0.28, -2.75 /
195	C (MSbar) 1/2 strange sea; PDF parameters fro
196	DATA IORD(10)/2/, ALM(10)/0.237/, Q02(10)/4.00/
197	DATA (((AC(IEX,IPN,IFL,10), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
198	> / 1.84, 0.12,
199	> 6.34, -2.40, -0.97, 0.22, -0.34, 1.53, 0.03, -0.04,
200	> 0.25, -0.16, 2.08, 0.02, 4.53, -2.51, -0.66, 0.19,
201	> -0.04, 1.24, -0.02, -0.01, 0.15, -0.05, 0.31, -0.10,
202	> 4.18, 0.34, -1.84, -0.10, 0.05, 1.64, -0.06, 0.01,
203	> -0.12, -0.16, -1.13, -0.15, 8.43, -0.64, -1.26, -0.16,
204	> -0.39, 1.01, -0.01, 0.03, 0.05, -0.06, -1.13, -0.15,
205	> 8.43, -0.64, -1.26, -0.16, -0.39, 1.01, -0.01, 0.03,
206	> 0.05, -0.06, -1.82, -0.18, 7.94, -0.56, -1.40, -0.06,
207	> -0.82, 1.65, 0.09, -0.01, 0.30, -0.31, -3.69, -0.15,
208	> 5.72, 0.26, -0.47, 0.04, 0.93, 1.85, -0.10, -0.05,
209	> -0.11, -0.50, -5.06, -0.25, 4.42, -0.14, 0.39, 0.16,
210	> 2.38, 2.72, -0.35, -0.08, -0.63, -0.75, -9.92, -0.38,
211	> -1.27, -1.60, 4.60, 0.24, 9.17, 4.40, -1.53, -0.08,
212	> -2.88, -1.08 /
213	C Set B0: PDF parameters fr
214	DATA IORD(11)/1/, ALM(11)/0.144/, Q02(11)/4.00/
215	DATA (((AC(IEX,IPN,IFL,11), IEX=0,3), IPN=0,2), IFL=2,-6,-1)
216	> / 1.38, 0.16,
217	> 5.40, -1.97, -0.62, 0.19, 0.59, 1.24, -0.10, -0.02,
218	> 0.03, -0.05, 1.67, 0.08, 3.75, -2.09, -0.33, 0.17,
219	> 0.70, 0.98, -0.13, -0.01, 0.00, 0.02, 1.52, -0.25,
220	> 7.01, -0.79, -3.17, -0.01, -0.90, 2.90, 0.25, 0.00,
221	> -0.08, -0.54, -0.81, -0.07, 9.19, -0.89, -1.13, -0.46,
222	> 0.35, 0.33, -0.26, 0.16, -0.49, 0.40, -0.81, -0.07,
223	> 9.19, -0.89, -1.13, -0.46, 0.35, 0.33, -0.26, 0.16,
224	> -0.49, 0.40, -0.81, -0.07, 9.19, -0.89, -1.13, -0.46,
225	> 0.35, 0.33, -0.26, 0.16, -0.49, 0.40, -3.62, -0.06,
226	> 8.30, 0.16, 0.03, -0.21, -0.60, 1.26, -0.48, 0.05,
227	> 0.25, -0.15, -6.16, -0.11, 6.49, 0.71, 2.37, -0.05,
228	> 1.28, 1.37, -1.24, -0.02, -0.41, -0.26,-12.68, -0.35,
229	> 14.87, -0.17, 8.36, 0.28,-12.56, 3.39, -2.89, -0.12,
230	> 4.75, -0.96 /
231
232	IRT = 0
233	IFL = LP
234	C
235	PDZXMT = 0.
236	IBLP = ABS(LP)
237	IF(IBLP .GE. 4 .AND. QQ .LE. 1.5) RETURN
238	IF(IBLP .GE. 5 .AND. QQ .LE. 5.0) RETURN
239	C
240	TM = TMAS
241	X = XX
242	Q0 = SQRT (Q02(IST))
243	ALAM = ALM(IST)
244	SQ = LOG ( LOG( MAX(Q0,QQ) /ALAM) / LOG(Q0/ALAM) )
245	IF(IBLP .EQ. 6) THEN
246	SQ = LOG ( LOG( MAX(Q0,QQ) /ALAM * (100./TM)) / LOG(Q0/ALAM))
247	ENDIF
248	C
249	FX(0) = EXP(D1)
250	FX(1) = X
251	FX(2) = 1.- X
252	FX(3) = LOG (1.+ 1./X)
253
254	PDFD = ONED
255	DO 20 IEX = 0, MEX(IST)
256	A(IEX) = AC(IEX, 0, IFL, IST)
257	SQN=SQ
258	DO 21 IPN = 1, MPN(IST)
259	A(IEX) = A(IEX) + AC(IEX, IPN, IFL, IST) * SQN
260	SQN = SQN*SQ
261	21 CONTINUE
262	FXIEX = FX(IEX)
263	AIEX = A(IEX)
264	PDFD = PDFD * FXIEX **(AIEX)
265	20 CONTINUE
266	C
267	C PDZXMT = PDF / X
268	PDF = PDFD
269	IF(PDF.GT.0.) PDZXMT = PDF
270	C
271	RETURN
272	C
273	ENTRY VLAMBD (ISET, IORDER)
274	C
275	IORDER = IORD (ISET)
276	VLAMBD = ALM (ISET)
277	C
278	RETURN
279	C *************************
280	END