* * $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