++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
mc_herwig_dijet.tcl
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# input module for generators
module input GenInputManager

mod disable Pythia
mod disable Isajet
mod disable Bgenerator
mod disable FAKE_EVENT

#-->--> Herwig DiJet Production at 1.96 TeV 

mod enable herwig
mod talk herwig
# 1500 is QCD 2 -> 2 hard parton scattering
  Process set 1500
# Following gets 1 line on screen per 100 events
#  report set 100
# seeds for herwig - note there may be TWO random number schemes ***
   Randseed1 set $random1
   Randseed2 set $random2
   Maxer set 5000
#
 Beam
   Beam1  set  P
   Beam2  set  PBAR
   Pbeam1 set 980.
   Pbeam2 set 980.
 exit
 Masses
   top set 175
 exit
#
#  Take default structure function
#  Structure functions, default is MRS 41 (cdf v6_202)
# Struc_Function
#   Autpdf set MRS
#   Modpdf set 41
# CTEQ Structure Functions
#  Autpdf set CTEQ
# - CTEQ 3L
#     Modpdf set  29
# - CTEQ 4M
#    Modpdf set  34
# - CTEQ 4HJ
#    Modpdf set  40
# exit

# Following 6 lines lead to a print (internal to herwig) of 1 event,
#   1 line per entity (so ~500 lines for ttbar), each line has width
#   ~120 characters. If use   prvtx set 0   then VHEP info is not printed
#    and width is ~80.
#   Default is no event printing (maxpr = 0).
 Prints
   help
   maxpr set 2
   prvtx set 0
   show
 exit
#
# Decays
#   help
# Following line sets pizero as stable, so get no pizero-to-gamma.gamma
#    decays. Default is pizero decays.
#   Pizstable set 1
# See manual (Herwig 6 sect 3.4 p 14) for Modbos info - inoperative on ttbar
#   Modbos1 set 4
#   Modbos2 set 4
#   show
# exit
#
# default Clpow leads to no B-Baryons
 Hadronization
    Clpow set 1.26
 exit

 Hards
    Ptmin set $pt_cut
 exit
#
 show
exit