001 ## Following is an example of running only Herwig (no simulation)
002 ## do something like:
003 ## cdfGen herwig_test.tcl > herwig_test.log
004 #############################
005 #show
006 #help
007 path enable AllPath
008 #
009 module input GenInputManager
010 # Following gets 1 line on screen per 100 events
011 talk GenInputManager
012 report set 100
013 exit
014
015 mod enable herwig
016 mod enable QQModule
017 mod enable TauolaModule
018
019 ## Herwig random number seeds (can also read in, see RestoreRandomStreams)
020 talk RandomGenManager
021 RandomNumberMenu
022 RandomSeedHERWIG1 set 2222222
023 RandomSeedHERWIG2 set 3333333
024 show
025 exit
026 exit
027
028
029
030 talk herwig
031 # help
032 # 1453 is qq' -> W -> tau.nutau
033 # 2300 is SM Higgs + jet prodn.
034 # 2100 is W + jet production
035 # 1500 is QCD 2 -> 2 hard parton scattering
036 Process set 1453
037 #
038 Beam
039 # Default Beam energies are 980. 980.
040 Pbeam1 set 980.
041 Pbeam2 set 980.
042 exit
043 Masses
044 top set 175
045 exit
046 #
047 # Structure functions, default is MRST 1998 LO (v6_4)
048 # Struc_Function
049 # Autpdf set MRS
050 # Modpdf set 41
051 # exit
052 #
053 # Following 6 lines lead to a print (internal to herwig) of 1 event,
054 # 1 line per entity (so ~500 lines for ttbar), each line has width
055 # ~120 characters. If use prvtx set 0 then VHEP info is not printed
056 # and width is ~80.
057 # Default is no event printing (maxpr = 0).
058 Prints
059 help
060 maxpr set 2
061 prvtx set 0
062 show
063 exit
064 #
065 # Decays
066 # help
067 # Following line sets pizero as stable, so get no pizero-to-gamma.gamma
068 # decays. Default is pizero decays.
069 # Pizstable set 1
070 # See manual (Herwig 6 sect 3.4 p 14) for Modbos info - inoperative on ttbar
071 # Modbos1 set 4
072 # Modbos2 set 4
073 # show
074 # exit
075 #
076 # default Clpow leads to no B-Baryons
077 Hadronization
078 Clpow set 1.26
079 exit
080 #
081 show
082 exit
083 ###
084
085 # note default is to NOT call qq, but BEWARE if you do not use qq
086 talk QQModule
087 Decay_Bc set 1
088 Decay_Upsilon set 1
089 Decay_B_Baryons set 1
090 Decay_prompt_charm set 1
091 # show
092 exit
093
094
095 # To make an output file (herwig_tau.out), uncomment the following lines:
096 #
097
098 ##module talk DHOutput
099 ## output create myout herwig_tau.out
100 ## output enable myout
101 ## output path myout AllPath
102 ## output keepList myout EVCL_StorableBank LRIH_StorableBank HEPG_StorableBank
103 ## output list
104 ##exit
105
106 ###
107 # Following gets 10 events; change 10 to 1000 for 1000 events, etc.
108 begin -nev 10
109 show timer
110 exit
111 #########################
112
113
114
Send problems or questions to cdfcode@fnal.gov
