ALIROOT/pydata.F

   1
   2 C*********************************************************************
   3
   4       BLOCK DATA PYDATA
   5
   6 C...Give sensible default values to all status codes and parameters.
   7       COMMON/PYSUBS/MSEL,MSUB(200),KFIN(2,-40:40),CKIN(200)
   8       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
   9       COMMON/PYINT1/MINT(400),VINT(400)
  10       COMMON/PYINT2/ISET(200),KFPR(200,2),COEF(200,20),ICOL(40,4,2)
  11       COMMON/PYINT3/XSFX(2,-40:40),ISIG(1000,3),SIGH(1000)
  12       COMMON/PYINT4/WIDP(21:40,0:40),WIDE(21:40,0:40),WIDS(21:40,3)
  13       COMMON/PYINT5/NGEN(0:200,3),XSEC(0:200,3)
  14       COMMON/PYINT6/PROC(0:200)
  15       CHARACTER PROC*28
  16       COMMON/PYINT7/SIGT(0:6,0:6,0:5)
  17       SAVE /PYSUBS/,/PYPARS/,/PYINT1/,/PYINT2/,/PYINT3/,/PYINT4/,
  18      &/PYINT5/,/PYINT6/,/PYINT7/
  19
  20 C...Default values for allowed processes and kinematics constraints.
  21       DATA MSEL/1/
  22       DATA MSUB/200*0/
  23       DATA ((KFIN(I,J),J=-40,40),I=1,2)/16*0,4*1,4*0,6*1,5*0,5*1,0,
  24      &5*1,5*0,6*1,4*0,4*1,16*0,16*0,4*1,4*0,6*1,5*0,5*1,0,5*1,5*0,
  25      &6*1,4*0,4*1,16*0/
  26       DATA CKIN/
  27      &   2.0, -1.0,  0.0, -1.0,  1.0,  1.0, -10.,  10., -10.,  10.,
  28      1  -10.,  10., -10.,  10., -10.,  10., -1.0,  1.0, -1.0,  1.0,
  29      2   0.0,  1.0,  0.0,  1.0, -1.0,  1.0, -1.0,  1.0,   0.,   0.,
  30      3   2.0, -1.0,   0.,   0.,  0.0, -1.0,  0.0, -1.0,  4.0, -1.0,
  31      4  12.0, -1.0, 12.0, -1.0, 12.0, -1.0, 12.0, -1.0,   0.,   0.,
  32      5   0.0, -1.0,  0.0, -1.0,  0.0, -1.0,   0.,   0.,   0.,   0.,
  33      6   140*0./
  34
  35 C...Default values for main switches and parameters. Reset information.
  36       DATA (MSTP(I),I=1,100)/
  37      &     3,    1,    2,    0,    0,    0,    0,    0,    0,    0,
  38      1     1,    0,    1,    0,    5,    0,    0,    0,    0,    0,
  39      2     1,    0,    1,    0,    0,    0,    0,    0,    0,    1,
  40      3     1,    2,    0,    1,    0,    2,    1,    5,    2,    0,
  41      4     1,    1,    3,    7,    3,    1,    1,    2,    2,    0,
  42      5     9,    1,    1,    1,    5,    1,    1,    6,    1,    0,
  43      6     1,    3,    2,    2,    1,    1,    2,    0,    0,    0,
  44      7     1,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  45      8     1,    1,  100,    0,    0,    0,    0,    0,    0,    0,
  46      9     1,    4,    1,    2,    0,    0,    0,    0,    0,    0/
  47       DATA (MSTP(I),I=101,200)/
  48      &     3,    1,    0,    0,    0,    0,    0,    0,    0,    0,
  49      1     1,    1,    1,    0,    0,    0,    0,    0,    0,    0,
  50      2     0,    1,    2,    1,    1,   20,    0,    0,   10,    0,
  51      3     0,    4,    0,    1,    0,    0,    0,    0,    0,    0,
  52      4     0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  53      5     0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  54      6     0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  55      7     0,    2,    0,    0,    0,    0,    0,    0,    0,    0,
  56      8     5,  720, 1995,   11,   29,  408,    0,    0,    0,    0,
  57      9     0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
  58       DATA (PARP(I),I=1,100)/
  59      &  0.25,  10.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  60      1    0.,   0.,  1.0, 0.01,  0.6,  1.0,  1.0,   0.,   0.,   0.,
  61      2    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  62      3   1.5,  2.0, 0.075, 1.0,  0.2,   0.,  2.0, 0.70, 0.006,  0.,
  63      4  0.02,  2.0, 0.10, 1000., 2054., 123., 246., 0.,   0.,   0.,
  64      5   1.0,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  65      6  0.25,  1.0, 0.25,  1.0,  2.0, 1E-3,  4.0, 1E-3,   0.,   0.,
  66      7   4.0, 0.25,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  67      8  1.40, 1.55,  0.5,  0.2, 0.33, 0.66,  0.7,  0.5,   0.,   0.,
  68      9  0.44, 0.20,  2.0,  1.0,   0.,  3.0,  1.0, 0.75, 0.44,  2.0/
  69       DATA (PARP(I),I=101,200)/
  70      &   0.5, 0.28,  1.0,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  71      1   2.0,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  72      2   1.0,  0.4,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  73      3  0.01,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  74      4    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  75      5    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  76      6  2.20, 23.6, 18.4, 11.5,   0.,   0.,   0.,   0.,   0.,   0.,
  77      7    0.,   0.,   0.,  1.0,   0.,   0.,   0.,   0.,   0.,   0.,
  78      8    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  79      9    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0./
  80       DATA MSTI/200*0/
  81       DATA PARI/200*0./
  82       DATA MINT/400*0/
  83       DATA VINT/400*0./
  84
  85 C...Constants for the generation of the various processes.
  86       DATA (ISET(I),I=1,100)/
  87      &    1,    1,    1,   -1,    3,   -1,   -1,    3,   -2,    2,
  88      1    2,    2,    2,    2,    2,    2,   -1,    2,    2,    2,
  89      2   -1,    2,    2,    2,    2,    2,   -1,    2,    2,    2,
  90      3    2,   -1,    2,    2,    2,    2,   -1,   -1,   -1,   -1,
  91      4   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,
  92      5   -1,   -1,    2,    2,   -1,   -1,   -1,    2,   -1,   -1,
  93      6   -1,   -1,   -1,   -1,   -1,   -1,   -1,    2,    2,    2,
  94      7    4,    4,    4,   -1,   -1,    4,    4,   -1,   -1,    2,
  95      8    2,    2,    2,    2,    2,    2,    2,    2,    2,   -2,
  96      9    0,    0,    0,    0,    0,    9,   -2,   -2,   -2,   -2/
  97       DATA (ISET(I),I=101,200)/
  98      &   -1,    1,    1,   -2,   -2,   -2,   -2,   -2,   -2,    2,
  99      1    2,    2,    2,    2,    2,   -1,   -1,   -1,   -2,   -2,
 100      2    5,    5,    5,    5,   -2,   -2,   -2,   -2,   -2,   -2,
 101      3    6,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,
 102      4    1,    1,    1,    1,    1,   -2,    1,    1,    1,   -2,
 103      5    1,    1,    1,   -2,   -2,    1,    1,    1,   -2,   -2,
 104      6    2,    2,    2,    2,    2,    2,    2,    2,   -2,   -2,
 105      7    2,    2,    5,    5,   -2,    2,    2,    5,    5,   -2,
 106      8    5,    5,   -2,   -2,   -2,    5,    5,   -2,   -2,   -2,
 107      9   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2/
 108       DATA ((KFPR(I,J),J=1,2),I=1,50)/
 109      &   23,    0,   24,    0,   25,    0,   24,    0,   25,    0,
 110      &   24,    0,   23,    0,   25,    0,    0,    0,    0,    0,
 111      1    0,    0,    0,    0,   21,   21,   21,   22,   21,   23,
 112      1   21,   24,   21,   25,   22,   22,   22,   23,   22,   24,
 113      2   22,   25,   23,   23,   23,   24,   23,   25,   24,   24,
 114      2   24,   25,   25,   25,    0,   21,    0,   22,    0,   23,
 115      3    0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
 116      3    0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
 117      4    0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
 118      4    0,   24,    0,   25,    0,   21,    0,   22,    0,   23/
 119       DATA ((KFPR(I,J),J=1,2),I=51,100)/
 120      5    0,   24,    0,   25,    0,    0,    0,    0,    0,    0,
 121      5    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 122      6    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 123      6    0,    0,    0,    0,   21,   21,   24,   24,   23,   24,
 124      7   23,   23,   24,   24,   23,   24,   23,   25,   22,   22,
 125      7   23,   23,   24,   24,   24,   25,   25,   25,    0,  211,
 126      8    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 127      8  443,   21,10441,   21,20443,   21,  445,   21,    0,    0,
 128      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 129      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
 130       DATA ((KFPR(I,J),J=1,2),I=101,150)/
 131      &   23,    0,   25,    0,   25,    0,    0,    0,    0,    0,
 132      &    0,    0,    0,    0,    0,    0,    0,    0,   22,   25,
 133      1   21,   25,    0,   25,   21,   25,   22,   22,   21,   22,
 134      1   22,   23,   23,   23,   24,   24,    0,    0,    0,    0,
 135      2   25,    6,   25,    6,   25,    0,   25,    0,    0,    0,
 136      2    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 137      3   23,    5,    0,    0,    0,    0,    0,    0,    0,    0,
 138      3    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 139      4   32,    0,   34,    0,   37,    0,   40,    0,   39,    0,
 140      4    0,    0,    7,    0,    8,    0,   38,    0,    0,    0/
 141       DATA ((KFPR(I,J),J=1,2),I=151,200)/
 142      5   35,    0,   35,    0,   35,    0,    0,    0,    0,    0,
 143      5   36,    0,   36,    0,   36,    0,    0,    0,    0,    0,
 144      6    6,   37,   39,    0,   39,   39,   39,   39,   11,    0,
 145      6   11,    0,    0,    7,    0,    8,    0,    0,    0,    0,
 146      7   23,   35,   24,   35,   35,    0,   35,    0,    0,    0,
 147      7   23,   36,   24,   36,   36,    0,   36,    0,    0,    0,
 148      8   35,    6,   35,    6,    0,    0,    0,    0,    0,    0,
 149      8   36,    6,   36,    6,    0,    0,    0,    0,    0,    0,
 150      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 151      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
 152       DATA COEF/4000*0./
 153       DATA (((ICOL(I,J,K),K=1,2),J=1,4),I=1,40)/
 154      1 4,0,3,0,2,0,1,0,3,0,4,0,1,0,2,0,2,0,0,1,4,0,0,3,3,0,0,4,1,0,0,2,
 155      2 3,0,0,4,1,4,3,2,4,0,0,3,4,2,1,3,2,0,4,1,4,0,2,3,4,0,3,4,2,0,1,2,
 156      3 3,2,1,0,1,4,3,0,4,3,3,0,2,1,1,0,3,2,1,4,1,0,0,2,2,4,3,1,2,0,0,1,
 157      4 3,2,1,4,1,4,3,2,4,2,1,3,4,2,1,3,3,4,4,3,1,2,2,1,2,0,3,1,2,0,0,0,
 158      5 4,2,1,0,0,0,1,0,3,0,0,3,1,2,0,0,4,0,0,4,0,0,1,2,2,0,0,1,4,4,3,3,
 159      6 2,2,1,1,4,4,3,3,3,3,4,4,1,1,2,2,3,2,1,3,1,2,0,0,4,2,1,4,0,0,1,2,
 160      7 4,0,0,0,4,0,1,3,0,0,3,0,2,4,3,0,3,4,0,0,1,0,0,1,0,0,3,4,2,0,0,2,
 161      8 3,0,0,0,1,0,0,0,0,0,3,0,2,0,0,0,2,0,3,1,2,0,0,0,3,2,1,0,1,0,0,0,
 162      9 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
 163      & 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0/
 164
 165 C...Character constants: name of processes.
 166       DATA PROC(0)/                    'All included subprocesses   '/
 167       DATA (PROC(I),I=1,20)/
 168      1'f + f~ -> gamma*/Z0         ',  'f + f~'' -> W+/-             ',
 169      2'f + f~ -> H0                ',  'gamma + W+/- -> W+/-        ',
 170      3'Z0 + Z0 -> H0               ',  'Z0 + W+/- -> W+/-           ',
 171      4'                            ',  'W+ + W- -> H0               ',
 172      5'                            ',  'f + f'' -> f + f'' (QFD)      ',
 173      6'f + f'' -> f + f'' (QCD)      ','f + f~ -> f'' + f~''          ',
 174      7'f + f~ -> g + g             ',  'f + f~ -> g + gamma         ',
 175      8'f + f~ -> g + Z0            ',  'f + f~'' -> g + W+/-         ',
 176      9'f + f~ -> g + H0            ',  'f + f~ -> gamma + gamma     ',
 177      &'f + f~ -> gamma + Z0        ',  'f + f~'' -> gamma + W+/-     '/
 178       DATA (PROC(I),I=21,40)/
 179      1'f + f~ -> gamma + H0        ',  'f + f~ -> Z0 + Z0           ',
 180      2'f + f~'' -> Z0 + W+/-        ', 'f + f~ -> Z0 + H0           ',
 181      3'f + f~ -> W+ + W-           ',  'f + f~'' -> W+/- + H0        ',
 182      4'f + f~ -> H0 + H0           ',  'f + g -> f + g              ',
 183      5'f + g -> f + gamma          ',  'f + g -> f + Z0             ',
 184      6'f + g -> f'' + W+/-          ', 'f + g -> f + H0             ',
 185      7'f + gamma -> f + g          ',  'f + gamma -> f + gamma      ',
 186      8'f + gamma -> f + Z0         ',  'f + gamma -> f'' + W+/-      ',
 187      9'f + gamma -> f + H0         ',  'f + Z0 -> f + g             ',
 188      &'f + Z0 -> f + gamma         ',  'f + Z0 -> f + Z0            '/
 189       DATA (PROC(I),I=41,60)/
 190      1'f + Z0 -> f'' + W+/-         ', 'f + Z0 -> f + H0            ',
 191      2'f + W+/- -> f'' + g          ', 'f + W+/- -> f'' + gamma      ',
 192      3'f + W+/- -> f'' + Z0         ', 'f + W+/- -> f'' + W+/-       ',
 193      4'f + W+/- -> f'' + H0         ', 'f + H0 -> f + g             ',
 194      5'f + H0 -> f + gamma         ',  'f + H0 -> f + Z0            ',
 195      6'f + H0 -> f'' + W+/-         ', 'f + H0 -> f + H0            ',
 196      7'g + g -> f + f~             ',  'g + gamma -> f + f~         ',
 197      8'g + Z0 -> f + f~            ',  'g + W+/- -> f + f~''         ',
 198      9'g + H0 -> f + f~            ',  'gamma + gamma -> f + f~     ',
 199      &'gamma + Z0 -> f + f~        ',  'gamma + W+/- -> f + f~''     '/
 200       DATA (PROC(I),I=61,80)/
 201      1'gamma + H0 -> f + f~        ',  'Z0 + Z0 -> f + f~           ',
 202      2'Z0 + W+/- -> f + f~''        ', 'Z0 + H0 -> f + f~           ',
 203      3'W+ + W- -> f + f~           ',  'W+/- + H0 -> f + f~''        ',
 204      4'H0 + H0 -> f + f~           ',  'g + g -> g + g              ',
 205      5'gamma + gamma -> W+ + W-    ',  'gamma + W+/- -> Z0 + W+/-   ',
 206      6'Z0 + Z0 -> Z0 + Z0          ',  'Z0 + Z0 -> W+ + W-          ',
 207      7'Z0 + W+/- -> Z0 + W+/-      ',  'Z0 + Z0 -> Z0 + H0          ',
 208      8'W+ + W- -> gamma + gamma    ',  'W+ + W- -> Z0 + Z0          ',
 209      9'W+/- + W+/- -> W+/- + W+/-  ',  'W+/- + H0 -> W+/- + H0      ',
 210      &'H0 + H0 -> H0 + H0          ',  'q + gamma -> q'' + pi+/-     '/
 211       DATA (PROC(I),I=81,100)/
 212      1'q + q~ -> Q + Q~, massive   ',  'g + g -> Q + Q~, massive    ',
 213      2'f + q -> f'' + Q, massive    ', 'g + gamma -> Q + Q~, massive',
 214      3'gamma + gamma -> F + F~, mas',  'g + g -> J/Psi + g          ',
 215      4'g + g -> chi_0c + g         ',  'g + g -> chi_1c + g         ',
 216      5'g + g -> chi_2c + g         ',  '                            ',
 217      6'Elastic scattering          ',  'Single diffractive (XB)     ',
 218      7'Single diffractive (AX)     ',  'Double  diffractive         ',
 219      8'Low-pT scattering           ',  'Semihard QCD 2 -> 2         ',
 220      9'                            ',  '                            ',
 221      &'                            ',  '                            '/
 222       DATA (PROC(I),I=101,120)/
 223      1'g + g -> gamma*/Z0          ',  'g + g -> H0                 ',
 224      2'gamma + gamma -> H0         ',  '                            ',
 225      3'                            ',  '                            ',
 226      4'                            ',  '                            ',
 227      5'                            ',  'f + f~ -> gamma + H0        ',
 228      6'f + f~ -> g + H0            ',  'q + g -> q + H0             ',
 229      7'g + g -> g + H0             ',  'g + g -> gamma + gamma      ',
 230      8'g + g -> g + gamma          ',  'g + g -> gamma + Z0         ',
 231      9'g + g -> Z0 + Z0            ',  'g + g -> W+ + W-            ',
 232      &'                            ',  '                            '/
 233       DATA (PROC(I),I=121,140)/
 234      1'g + g -> Q + Q~ + H0        ',  'q + q~ -> Q + Q~ + H0       ',
 235      2'f + f'' -> f + f'' + H0       ',
 236      2'f + f'' -> f" + f"'' + H0     ',
 237      3'                            ',  '                            ',
 238      4'                            ',  '                            ',
 239      5'                            ',  '                            ',
 240      6'g + g -> Z0 + q + q~        ',  '                            ',
 241      7'                            ',  '                            ',
 242      8'                            ',  '                            ',
 243      9'                            ',  '                            ',
 244      &'                            ',  '                            '/
 245       DATA (PROC(I),I=141,160)/
 246      1'f + f~ -> gamma*/Z0/Z''0     ', 'f + f~'' -> W''+/-            ',
 247      2'f + f~'' -> H+/-             ', 'f + f~'' -> R                ',
 248      3'q + l -> LQ                 ',  '                            ',
 249      4'd + g -> d*                 ',  'u + g -> u*                 ',
 250      5'g + g -> eta_techni         ',  '                            ',
 251      6'f + f~ -> H''0               ', 'g + g -> H''0                ',
 252      7'gamma + gamma -> H''0        ', '                            ',
 253      8'                            ',  'f + f~ -> A0                ',
 254      9'g + g -> A0                 ',  'gamma + gamma -> A0         ',
 255      &'                            ',  '                            '/
 256       DATA (PROC(I),I=161,180)/
 257      1'f + g -> f'' + H+/-          ', 'q + g -> LQ + l~            ',
 258      2'g + g -> LQ + LQ~           ',  'q + q~ -> LQ + LQ~          ',
 259      3'f + f~ -> f'' + f~'' (gamma/Z)',
 260      3'f +f~'' -> f" + f~"'' (W)     ',
 261      4'q + q'' -> q" + d*           ',  'q + q'' -> q" + u*           ',
 262      5'                            ',  '                            ',
 263      6'f + f~ -> Z0 + H''0          ', 'f + f~'' -> W+/- + H''0       ',
 264      7'f + f'' -> f + f'' + H''0      ',
 265      7'f + f'' -> f" + f"'' + H''0    ',
 266      8'                            ',  'f + f~ -> Z0 + A0           ',
 267      9'f + f~'' -> W+/- + A0        ',
 268      9'f + f'' -> f + f'' + A0       ',
 269      &'f + f'' -> f" + f"'' + A0     ',
 270      &'                            '/
 271       DATA (PROC(I),I=181,200)/
 272      1'g + g -> Q + Q~ + H''0       ',  'q + q~ -> Q + Q~ + H''0      ',
 273      2'                            ',  '                            ',
 274      3'                            ',  'g + g -> Q + Q~ + A0        ',
 275      4'q + q~ -> Q + Q~ + A0       ',  '                            ',
 276      5'                            ',  '                            ',
 277      6'                            ',  '                            ',
 278      7'                            ',  '                            ',
 279      8'                            ',  '                            ',
 280      9'                            ',  '                            ',
 281      &'                            ',  '                            '/
 282
 283 C...Cross sections and slope offsets.
 284       DATA SIGT/294*0./
 285
 286       END