PDF/spdf/pdzxmt.F

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