HIJING/hipyset1_35/pydata_hijing.F

   1 * $Id$
   2
   3 C*********************************************************************
   4
   5       BLOCK DATA PYDATA_HIJING
   6
   7 C...Give sensible default values to all status codes and parameters.
   8 #include "pysubs_hijing.inc"
   9 #include "pypars_hijing.inc"
  10 #include "pyint1_hijing.inc"
  11 #include "pyint2_hijing.inc"
  12 #include "pyint3_hijing.inc"
  13 #include "pyint4_hijing.inc"
  14 #include "pyint5_hijing.inc"
  15 #include "pyint6_hijing.inc"
  16
  17 C...Default values for allowed processes and kinematics constraints.
  18       DATA MSEL/1/
  19       DATA MSUB/200*0/
  20       DATA ((KFIN(I,J),J=-40,40),I=1,2)/40*1,0,80*1,0,40*1/
  21       DATA CKIN/
  22      &   2.0, -1.0,  0.0, -1.0,  1.0,  1.0, -10.,  10., -10.,  10.,
  23      1  -10.,  10., -10.,  10., -10.,  10., -1.0,  1.0, -1.0,  1.0,
  24      2   0.0,  1.0,  0.0,  1.0, -1.0,  1.0, -1.0,  1.0,   0.,   0.,
  25      3   2.0, -1.0,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  26      4   160*0./
  27
  28 C...Default values for main switches and parameters. Reset information.
  29       DATA (MSTP(I),I=1,100)/
  30      &     3,    1,    2,    0,    0,    0,    0,    0,    0,    0,
  31      1     0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  32      2     0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  33      3     1,    2,    0,    0,    0,    2,    0,    0,    0,    0,
  34      4     1,    0,    3,    7,    1,    0,    0,    0,    0,    0,
  35      5     1,    1,   20,    6,    0,    0,    0,    0,    0,    0,
  36      6     1,    2,    2,    2,    1,    0,    0,    0,    0,    0,
  37      7     1,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  38      8     1,    1,  100,    0,    0,    0,    0,    0,    0,    0,
  39      9     1,    4,    0,    0,    0,    0,    0,    0,    0,    0/
  40       DATA (MSTP(I),I=101,200)/
  41      &     1,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  42      1     1,    1,    1,    0,    0,    0,    0,    0,    0,    0,
  43      2     0,    1,    2,    1,    1,   20,    0,    0,    0,    0,
  44      3     0,    4,    0,    1,    0,    0,    0,    0,    0,    0,
  45      4     0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  46      5     0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  47      6     0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  48      7     0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
  49      8     5,    3, 1989,   11,   24,    0,    0,    0,    0,    0,
  50      9     0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
  51       DATA (PARP(I),I=1,100)/
  52      &  0.25,  10.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  53      1    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  54      2    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  55      3   1.5,  2.0, 0.075,  0.,  0.2,   0.,   0.,   0.,   0.,   0.,
  56      4    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  57      5   1.0, 2.26, 1.E4, 1.E-4,  0.,   0.,   0.,   0.,   0.,   0.,
  58      6  0.25,  1.0, 0.25,  1.0,  2.0, 1.E-3, 4.0,   0.,   0.,   0.,
  59      7   4.0,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  60      8   1.6, 1.85,  0.5,  0.2, 0.33, 0.66,  0.7,  0.5,   0.,   0.,
  61      9  0.44, 0.44,  2.0,  1.0,   0.,  3.0,  1.0, 0.75,   0.,   0./
  62       DATA (PARP(I),I=101,200)/
  63      & -0.02,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  64      1   2.0,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  65      2   0.4,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  66      3  0.01,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  67      4    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  68      5    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  69      6    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  70      7    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  71      8    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,
  72      9    0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0.,   0./
  73       DATA MSTI/200*0/
  74       DATA PARI/200*0./
  75       DATA MINT/400*0/
  76       DATA VINT/400*0./
  77
  78 C...Constants for the generation of the various processes.
  79       DATA (ISET(I),I=1,100)/
  80      &    1,    1,    1,   -1,    3,   -1,   -1,    3,   -2,   -2,
  81      1    2,    2,    2,    2,    2,    2,   -1,    2,    2,    2,
  82      2   -1,    2,    2,    2,    2,    2,   -1,    2,    2,    2,
  83      3    2,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,
  84      4   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,   -1,
  85      5   -1,   -1,    2,   -1,   -1,   -1,   -1,   -1,   -1,   -1,
  86      6   -1,   -1,   -1,   -1,   -1,   -1,   -1,    2,   -1,   -1,
  87      7    4,    4,    4,   -1,   -1,    4,    4,   -1,   -1,   -2,
  88      8    2,    2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,
  89      9    0,    0,    0,   -1,    0,    5,   -2,   -2,   -2,   -2/
  90       DATA (ISET(I),I=101,200)/
  91      &   -1,    1,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,
  92      1    2,    2,    2,    2,   -1,   -1,   -1,   -2,   -2,   -2,
  93      2   -1,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,
  94      3   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,
  95      4    1,    1,    1,   -2,   -2,   -2,   -2,   -2,   -2,   -2,
  96      5   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,
  97      6    2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,
  98      7   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,
  99      8   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,
 100      9   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2,   -2/
 101       DATA ((KFPR(I,J),J=1,2),I=1,50)/
 102      &   23,    0,   24,    0,   25,    0,   24,    0,   25,    0,
 103      &   24,    0,   23,    0,   25,    0,    0,    0,    0,    0,
 104      1    0,    0,    0,    0,   21,   21,   21,   22,   21,   23,
 105      1   21,   24,   21,   25,   22,   22,   22,   23,   22,   24,
 106      2   22,   25,   23,   23,   23,   24,   23,   25,   24,   24,
 107      2   24,   25,   25,   25,    0,   21,    0,   22,    0,   23,
 108      3    0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
 109      3    0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
 110      4    0,   24,    0,   25,    0,   21,    0,   22,    0,   23,
 111      4    0,   24,    0,   25,    0,   21,    0,   22,    0,   23/
 112       DATA ((KFPR(I,J),J=1,2),I=51,100)/
 113      5    0,   24,    0,   25,    0,    0,    0,    0,    0,    0,
 114      5    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 115      6    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 116      6    0,    0,    0,    0,   21,   21,   24,   24,   22,   24,
 117      7   23,   23,   24,   24,   23,   24,   23,   25,   22,   22,
 118      7   23,   23,   24,   24,   24,   25,   25,   25,    0,    0,
 119      8    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 120      8    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 121      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 122      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
 123       DATA ((KFPR(I,J),J=1,2),I=101,150)/
 124      &   23,    0,   25,    0,    0,    0,    0,    0,    0,    0,
 125      &    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 126      1   21,   25,    0,   25,   21,   25,   22,   22,   22,   23,
 127      1   23,   23,   24,   24,    0,    0,    0,    0,    0,    0,
 128      2    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 129      2    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 130      3    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 131      3    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 132      4   32,    0,   37,    0,   40,    0,    0,    0,    0,    0,
 133      4    0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
 134       DATA ((KFPR(I,J),J=1,2),I=151,200)/
 135      5    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 136      5    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 137      6    0,   37,    0,    0,    0,    0,    0,    0,    0,    0,
 138      6    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 139      7    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 140      7    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 141      8    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 142      8    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 143      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0,
 144      9    0,    0,    0,    0,    0,    0,    0,    0,    0,    0/
 145       DATA COEF/4000*0./
 146       DATA (((ICOL(I,J,K),K=1,2),J=1,4),I=1,40)/
 147      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,
 148      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,
 149      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,
 150      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,
 151      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,
 152      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,
 153      7 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,
 154      8 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,
 155      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,
 156      & 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/
 157
 158 C...Character constants: name of processes.
 159       DATA PROC(0)/                    'All included subprocesses   '/
 160       DATA (PROC(I),I=1,20)/
 161      1'f + fb -> gamma*/Z0         ',  'f + fb'' -> W+/-             ',
 162      2'f + fb -> H0                ',  'gamma + W+/- -> W+/-        ',
 163      3'Z0 + Z0 -> H0               ',  'Z0 + W+/- -> W+/-           ',
 164      4'                            ',  'W+ + W- -> H0               ',
 165      5'                            ',  '                            ',
 166      6'f + f'' -> f + f''            ','f + fb -> f'' + fb''          ',
 167      7'f + fb -> g + g             ',  'f + fb -> g + gamma         ',
 168      8'f + fb -> g + Z0            ',  'f + fb'' -> g + W+/-         ',
 169      9'f + fb -> g + H0            ',  'f + fb -> gamma + gamma     ',
 170      &'f + fb -> gamma + Z0        ',  'f + fb'' -> gamma + W+/-     '/
 171       DATA (PROC(I),I=21,40)/
 172      1'f + fb -> gamma + H0        ',  'f + fb -> Z0 + Z0           ',
 173      2'f + fb'' -> Z0 + W+/-        ', 'f + fb -> Z0 + H0           ',
 174      3'f + fb -> W+ + W-           ',  'f + fb'' -> W+/- + H0        ',
 175      4'f + fb -> H0 + H0           ',  'f + g -> f + g              ',
 176      5'f + g -> f + gamma          ',  'f + g -> f + Z0             ',
 177      6'f + g -> f'' + W+/-          ', 'f + g -> f + H0             ',
 178      7'f + gamma -> f + g          ',  'f + gamma -> f + gamma      ',
 179      8'f + gamma -> f + Z0         ',  'f + gamma -> f'' + W+/-      ',
 180      9'f + gamma -> f + H0         ',  'f + Z0 -> f + g             ',
 181      &'f + Z0 -> f + gamma         ',  'f + Z0 -> f + Z0            '/
 182       DATA (PROC(I),I=41,60)/
 183      1'f + Z0 -> f'' + W+/-         ', 'f + Z0 -> f + H0            ',
 184      2'f + W+/- -> f'' + g          ', 'f + W+/- -> f'' + gamma      ',
 185      3'f + W+/- -> f'' + Z0         ', 'f + W+/- -> f'' + W+/-       ',
 186      4'f + W+/- -> f'' + H0         ', 'f + H0 -> f + g             ',
 187      5'f + H0 -> f + gamma         ',  'f + H0 -> f + Z0            ',
 188      6'f + H0 -> f'' + W+/-         ', 'f + H0 -> f + H0            ',
 189      7'g + g -> f + fb             ',  'g + gamma -> f + fb         ',
 190      8'g + Z0 -> f + fb            ',  'g + W+/- -> f + fb''         ',
 191      9'g + H0 -> f + fb            ',  'gamma + gamma -> f + fb     ',
 192      &'gamma + Z0 -> f + fb        ',  'gamma + W+/- -> f + fb''     '/
 193       DATA (PROC(I),I=61,80)/
 194      1'gamma + H0 -> f + fb        ',  'Z0 + Z0 -> f + fb           ',
 195      2'Z0 + W+/- -> f + fb''        ', 'Z0 + H0 -> f + fb           ',
 196      3'W+ + W- -> f + fb           ',  'W+/- + H0 -> f + fb''        ',
 197      4'H0 + H0 -> f + fb           ',  'g + g -> g + g              ',
 198      5'gamma + gamma -> W+ + W-    ',  'gamma + W+/- -> gamma + W+/-',
 199      6'Z0 + Z0 -> Z0 + Z0          ',  'Z0 + Z0 -> W+ + W-          ',
 200      7'Z0 + W+/- -> Z0 + W+/-      ',  'Z0 + Z0 -> Z0 + H0          ',
 201      8'W+ + W- -> gamma + gamma    ',  'W+ + W- -> Z0 + Z0          ',
 202      9'W+/- + W+/- -> W+/- + W+/-  ',  'W+/- + H0 -> W+/- + H0      ',
 203      &'H0 + H0 -> H0 + H0          ',  '                            '/
 204       DATA (PROC(I),I=81,100)/
 205      1'q + qb -> Q + QB, massive   ',  'g + g -> Q + QB, massive    ',
 206      2'                            ',  '                            ',
 207      3'                            ',  '                            ',
 208      4'                            ',  '                            ',
 209      5'                            ',  '                            ',
 210      6'Elastic scattering          ',  'Single diffractive          ',
 211      7'Double diffractive          ',  'Central diffractive         ',
 212      8'Low-pT scattering           ',  'Semihard QCD 2 -> 2         ',
 213      9'                            ',  '                            ',
 214      &'                            ',  '                            '/
 215       DATA (PROC(I),I=101,120)/
 216      1'g + g -> gamma*/Z0          ',  'g + g -> H0                 ',
 217      2'                            ',  '                            ',
 218      3'                            ',  '                            ',
 219      4'                            ',  '                            ',
 220      5'                            ',  '                            ',
 221      6'f + fb -> g + H0            ',  'q + g -> q + H0             ',
 222      7'g + g -> g + H0             ',  'g + g -> gamma + gamma      ',
 223      8'g + g -> gamma + Z0         ',  'g + g -> Z0 + Z0            ',
 224      9'g + g -> W+ + W-            ',  '                            ',
 225      &'                            ',  '                            '/
 226       DATA (PROC(I),I=121,140)/
 227      1'g + g -> f + fb + H0        ',  '                            ',
 228      2'                            ',  '                            ',
 229      3'                            ',  '                            ',
 230      4'                            ',  '                            ',
 231      5'                            ',  '                            ',
 232      6'                            ',  '                            ',
 233      7'                            ',  '                            ',
 234      8'                            ',  '                            ',
 235      9'                            ',  '                            ',
 236      &'                            ',  '                            '/
 237       DATA (PROC(I),I=141,160)/
 238      1'f + fb -> gamma*/Z0/Z''0     ', 'f + fb'' -> H+/-             ',
 239      2'f + fb -> R                 ',  '                            ',
 240      3'                            ',  '                            ',
 241      4'                            ',  '                            ',
 242      5'                            ',  '                            ',
 243      6'                            ',  '                            ',
 244      7'                            ',  '                            ',
 245      8'                            ',  '                            ',
 246      9'                            ',  '                            ',
 247      &'                            ',  '                            '/
 248       DATA (PROC(I),I=161,180)/
 249      1'f + g -> f'' + H+/-          ', '                            ',
 250      2'                            ',  '                            ',
 251      3'                            ',  '                            ',
 252      4'                            ',  '                            ',
 253      5'                            ',  '                            ',
 254      6'                            ',  '                            ',
 255      7'                            ',  '                            ',
 256      8'                            ',  '                            ',
 257      9'                            ',  '                            ',
 258      &'                            ',  '                            '/
 259       DATA (PROC(I),I=181,200)/     20*'                            '/
 260
 261       END