PDF/spdf/pdzxmt.F

   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 / NST*3, NST*2, 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