* $Id$ C********************************************************************* BLOCK DATA PYDATA_HIJING C...Give sensible default values to all status codes and parameters. #include "pysubs_hijing.inc" #include "pypars_hijing.inc" #include "pyint1_hijing.inc" #include "pyint2_hijing.inc" #include "pyint3_hijing.inc" #include "pyint4_hijing.inc" #include "pyint5_hijing.inc" #include "pyint6_hijing.inc" C...Default values for allowed processes and kinematics constraints. DATA MSEL/1/ DATA MSUB/200*0/ DATA ((KFIN(I,J),J=-40,40),I=1,2)/40*1,0,80*1,0,40*1/ DATA CKIN/ & 2.0, -1.0, 0.0, -1.0, 1.0, 1.0, -10., 10., -10., 10., 1 -10., 10., -10., 10., -10., 10., -1.0, 1.0, -1.0, 1.0, 2 0.0, 1.0, 0.0, 1.0, -1.0, 1.0, -1.0, 1.0, 0., 0., 3 2.0, -1.0, 0., 0., 0., 0., 0., 0., 0., 0., 4 160*0./ C...Default values for main switches and parameters. Reset information. DATA (MSTP(I),I=1,100)/ & 3, 1, 2, 0, 0, 0, 0, 0, 0, 0, 1 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3 1, 2, 0, 0, 0, 2, 0, 0, 0, 0, 4 1, 0, 3, 7, 1, 0, 0, 0, 0, 0, 5 1, 1, 20, 6, 0, 0, 0, 0, 0, 0, 6 1, 2, 2, 2, 1, 0, 0, 0, 0, 0, 7 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 8 1, 1, 100, 0, 0, 0, 0, 0, 0, 0, 9 1, 4, 0, 0, 0, 0, 0, 0, 0, 0/ DATA (MSTP(I),I=101,200)/ & 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 2 0, 1, 2, 1, 1, 20, 0, 0, 0, 0, 3 0, 4, 0, 1, 0, 0, 0, 0, 0, 0, 4 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 7 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 8 5, 3, 1989, 11, 24, 0, 0, 0, 0, 0, 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/ DATA (PARP(I),I=1,100)/ & 0.25, 10., 0., 0., 0., 0., 0., 0., 0., 0., 1 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 2 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 3 1.5, 2.0, 0.075, 0., 0.2, 0., 0., 0., 0., 0., 4 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 5 1.0, 2.26, 1.E4, 1.E-4, 0., 0., 0., 0., 0., 0., 6 0.25, 1.0, 0.25, 1.0, 2.0, 1.E-3, 4.0, 0., 0., 0., 7 4.0, 0., 0., 0., 0., 0., 0., 0., 0., 0., 8 1.6, 1.85, 0.5, 0.2, 0.33, 0.66, 0.7, 0.5, 0., 0., 9 0.44, 0.44, 2.0, 1.0, 0., 3.0, 1.0, 0.75, 0., 0./ DATA (PARP(I),I=101,200)/ & -0.02, 0., 0., 0., 0., 0., 0., 0., 0., 0., 1 2.0, 0., 0., 0., 0., 0., 0., 0., 0., 0., 2 0.4, 0., 0., 0., 0., 0., 0., 0., 0., 0., 3 0.01, 0., 0., 0., 0., 0., 0., 0., 0., 0., 4 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 5 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 6 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 7 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 8 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 9 0., 0., 0., 0., 0., 0., 0., 0., 0., 0./ DATA MSTI/200*0/ DATA PARI/200*0./ DATA MINT/400*0/ DATA VINT/400*0./ C...Constants for the generation of the various processes. DATA (ISET(I),I=1,100)/ & 1, 1, 1, -1, 3, -1, -1, 3, -2, -2, 1 2, 2, 2, 2, 2, 2, -1, 2, 2, 2, 2 -1, 2, 2, 2, 2, 2, -1, 2, 2, 2, 3 2, -1, -1, -1, -1, -1, -1, -1, -1, -1, 4 -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 5 -1, -1, 2, -1, -1, -1, -1, -1, -1, -1, 6 -1, -1, -1, -1, -1, -1, -1, 2, -1, -1, 7 4, 4, 4, -1, -1, 4, 4, -1, -1, -2, 8 2, 2, -2, -2, -2, -2, -2, -2, -2, -2, 9 0, 0, 0, -1, 0, 5, -2, -2, -2, -2/ DATA (ISET(I),I=101,200)/ & -1, 1, -2, -2, -2, -2, -2, -2, -2, -2, 1 2, 2, 2, 2, -1, -1, -1, -2, -2, -2, 2 -1, -2, -2, -2, -2, -2, -2, -2, -2, -2, 3 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, 4 1, 1, 1, -2, -2, -2, -2, -2, -2, -2, 5 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, 6 2, -2, -2, -2, -2, -2, -2, -2, -2, -2, 7 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, 8 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2, 9 -2, -2, -2, -2, -2, -2, -2, -2, -2, -2/ DATA ((KFPR(I,J),J=1,2),I=1,50)/ & 23, 0, 24, 0, 25, 0, 24, 0, 25, 0, & 24, 0, 23, 0, 25, 0, 0, 0, 0, 0, 1 0, 0, 0, 0, 21, 21, 21, 22, 21, 23, 1 21, 24, 21, 25, 22, 22, 22, 23, 22, 24, 2 22, 25, 23, 23, 23, 24, 23, 25, 24, 24, 2 24, 25, 25, 25, 0, 21, 0, 22, 0, 23, 3 0, 24, 0, 25, 0, 21, 0, 22, 0, 23, 3 0, 24, 0, 25, 0, 21, 0, 22, 0, 23, 4 0, 24, 0, 25, 0, 21, 0, 22, 0, 23, 4 0, 24, 0, 25, 0, 21, 0, 22, 0, 23/ DATA ((KFPR(I,J),J=1,2),I=51,100)/ 5 0, 24, 0, 25, 0, 0, 0, 0, 0, 0, 5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6 0, 0, 0, 0, 21, 21, 24, 24, 22, 24, 7 23, 23, 24, 24, 23, 24, 23, 25, 22, 22, 7 23, 23, 24, 24, 24, 25, 25, 25, 0, 0, 8 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 8 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/ DATA ((KFPR(I,J),J=1,2),I=101,150)/ & 23, 0, 25, 0, 0, 0, 0, 0, 0, 0, & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1 21, 25, 0, 25, 21, 25, 22, 22, 22, 23, 1 23, 23, 24, 24, 0, 0, 0, 0, 0, 0, 2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4 32, 0, 37, 0, 40, 0, 0, 0, 0, 0, 4 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/ DATA ((KFPR(I,J),J=1,2),I=151,200)/ 5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6 0, 37, 0, 0, 0, 0, 0, 0, 0, 0, 6 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 7 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 7 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 8 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 8 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 9 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/ DATA COEF/4000*0./ DATA (((ICOL(I,J,K),K=1,2),J=1,4),I=1,40)/ 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, 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, 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, 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, 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, 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, 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, 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, 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, & 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/ C...Character constants: name of processes. DATA PROC(0)/ 'All included subprocesses '/ DATA (PROC(I),I=1,20)/ 1'f + fb -> gamma*/Z0 ', 'f + fb'' -> W+/- ', 2'f + fb -> H0 ', 'gamma + W+/- -> W+/- ', 3'Z0 + Z0 -> H0 ', 'Z0 + W+/- -> W+/- ', 4' ', 'W+ + W- -> H0 ', 5' ', ' ', 6'f + f'' -> f + f'' ','f + fb -> f'' + fb'' ', 7'f + fb -> g + g ', 'f + fb -> g + gamma ', 8'f + fb -> g + Z0 ', 'f + fb'' -> g + W+/- ', 9'f + fb -> g + H0 ', 'f + fb -> gamma + gamma ', &'f + fb -> gamma + Z0 ', 'f + fb'' -> gamma + W+/- '/ DATA (PROC(I),I=21,40)/ 1'f + fb -> gamma + H0 ', 'f + fb -> Z0 + Z0 ', 2'f + fb'' -> Z0 + W+/- ', 'f + fb -> Z0 + H0 ', 3'f + fb -> W+ + W- ', 'f + fb'' -> W+/- + H0 ', 4'f + fb -> H0 + H0 ', 'f + g -> f + g ', 5'f + g -> f + gamma ', 'f + g -> f + Z0 ', 6'f + g -> f'' + W+/- ', 'f + g -> f + H0 ', 7'f + gamma -> f + g ', 'f + gamma -> f + gamma ', 8'f + gamma -> f + Z0 ', 'f + gamma -> f'' + W+/- ', 9'f + gamma -> f + H0 ', 'f + Z0 -> f + g ', &'f + Z0 -> f + gamma ', 'f + Z0 -> f + Z0 '/ DATA (PROC(I),I=41,60)/ 1'f + Z0 -> f'' + W+/- ', 'f + Z0 -> f + H0 ', 2'f + W+/- -> f'' + g ', 'f + W+/- -> f'' + gamma ', 3'f + W+/- -> f'' + Z0 ', 'f + W+/- -> f'' + W+/- ', 4'f + W+/- -> f'' + H0 ', 'f + H0 -> f + g ', 5'f + H0 -> f + gamma ', 'f + H0 -> f + Z0 ', 6'f + H0 -> f'' + W+/- ', 'f + H0 -> f + H0 ', 7'g + g -> f + fb ', 'g + gamma -> f + fb ', 8'g + Z0 -> f + fb ', 'g + W+/- -> f + fb'' ', 9'g + H0 -> f + fb ', 'gamma + gamma -> f + fb ', &'gamma + Z0 -> f + fb ', 'gamma + W+/- -> f + fb'' '/ DATA (PROC(I),I=61,80)/ 1'gamma + H0 -> f + fb ', 'Z0 + Z0 -> f + fb ', 2'Z0 + W+/- -> f + fb'' ', 'Z0 + H0 -> f + fb ', 3'W+ + W- -> f + fb ', 'W+/- + H0 -> f + fb'' ', 4'H0 + H0 -> f + fb ', 'g + g -> g + g ', 5'gamma + gamma -> W+ + W- ', 'gamma + W+/- -> gamma + W+/-', 6'Z0 + Z0 -> Z0 + Z0 ', 'Z0 + Z0 -> W+ + W- ', 7'Z0 + W+/- -> Z0 + W+/- ', 'Z0 + Z0 -> Z0 + H0 ', 8'W+ + W- -> gamma + gamma ', 'W+ + W- -> Z0 + Z0 ', 9'W+/- + W+/- -> W+/- + W+/- ', 'W+/- + H0 -> W+/- + H0 ', &'H0 + H0 -> H0 + H0 ', ' '/ DATA (PROC(I),I=81,100)/ 1'q + qb -> Q + QB, massive ', 'g + g -> Q + QB, massive ', 2' ', ' ', 3' ', ' ', 4' ', ' ', 5' ', ' ', 6'Elastic scattering ', 'Single diffractive ', 7'Double diffractive ', 'Central diffractive ', 8'Low-pT scattering ', 'Semihard QCD 2 -> 2 ', 9' ', ' ', &' ', ' '/ DATA (PROC(I),I=101,120)/ 1'g + g -> gamma*/Z0 ', 'g + g -> H0 ', 2' ', ' ', 3' ', ' ', 4' ', ' ', 5' ', ' ', 6'f + fb -> g + H0 ', 'q + g -> q + H0 ', 7'g + g -> g + H0 ', 'g + g -> gamma + gamma ', 8'g + g -> gamma + Z0 ', 'g + g -> Z0 + Z0 ', 9'g + g -> W+ + W- ', ' ', &' ', ' '/ DATA (PROC(I),I=121,140)/ 1'g + g -> f + fb + H0 ', ' ', 2' ', ' ', 3' ', ' ', 4' ', ' ', 5' ', ' ', 6' ', ' ', 7' ', ' ', 8' ', ' ', 9' ', ' ', &' ', ' '/ DATA (PROC(I),I=141,160)/ 1'f + fb -> gamma*/Z0/Z''0 ', 'f + fb'' -> H+/- ', 2'f + fb -> R ', ' ', 3' ', ' ', 4' ', ' ', 5' ', ' ', 6' ', ' ', 7' ', ' ', 8' ', ' ', 9' ', ' ', &' ', ' '/ DATA (PROC(I),I=161,180)/ 1'f + g -> f'' + H+/- ', ' ', 2' ', ' ', 3' ', ' ', 4' ', ' ', 5' ', ' ', 6' ', ' ', 7' ', ' ', 8' ', ' ', 9' ', ' ', &' ', ' '/ DATA (PROC(I),I=181,200)/ 20*' '/ END