001 ## wbbgen_test.tcl njets=3 04-Oct-2001
002 ## NOTE: for now only choice for ndns is 101 [which gives pdf set cteq5m];
003 ## MUST execute script generatorMods/examples/wbbgen/wbbgen_setup (or
004 ## similar) to get soft link to cteq5m.tbl .
005 ## But, if you do want another pdf, see wbbgen_i/doc/wbbgen_cdf.txt.
006 ## Following takes ~ 20 mins. to run, produces ~10 (temp) files
007 #############################
008 #show
009 #help
010 path enable AllPath
011 #
012 module input GenInputManager
013 ## Following gets 1 line on screen per 1 event (at hadronization stage)
014 talk GenInputManager
015 report set 1
016 exit
017 #
018 ## set herwig random number here; wbbgen random number in talk Wbbgen below
019 talk RandomGenManager
020 RandomNumberMenu
021 RandomSeedHERWIG1 set 1111110
022 RandomSeedHERWIG2 set 1111119
023 show
024 exit
025 exit
026
027
028 mod enable Wbbgen
029 mod enable QQModule
030
031 talk Wbbgen
032 # help
033 # show
034 #
035 ## Following 5 lines lead to a print (internal to herwig) of maxpr events,
036 ## 1 line per entity (so ~500 lines for ttbar), each line has width
037 ## ~120 characters. If use prvtx set 0 then VHEP info is not printed
038 ## and width is ~80.
039 ## Default is no event printing (maxpr = 0).
040 Prints
041 maxpr set 2
042 prvtx set 0
043 exit
044 #
045 ## default Clpow leads to no B-Baryons
046 Hadronization
047 Clpow set 1.3
048 exit
049 #
050 ## Here is the explicitly wbbgen talkto
051 Wbbuser
052 # help
053 ibeam2 set 1
054 ebeam set 980.
055 bmass set 4.75
056 njets set 3
057 ndns set 101
058 ptjmin set 15
059 ptbmin set 15
060 drjmin set 0.0
061 drbmin set 0.01
062 ptlmin set 0.1
063 etalmax set 9.9
064 drlmin set 0.01
065 ptnmin set 0.1
066 nopt set 200000
067 niter set 1
068 maxev set 200000
069 seed1 set 1111118
070 seed2 set 2222226
071 gridr set 0
072 gridw set 41
073 wbbfile set wbbfile1
074 exit
075 #
076 show
077 exit
078 ###
079 # note default is to NOT call qq, but BEWARE if you do not use qq
080 talk QQModule
081 Decay_Bc set 1
082 Decay_Upsilon set 1
083 Decay_B_Baryons set 1
084 Decay_prompt_charm set 1
085 # show
086 exit
087
088 ###
089 path list
090 #
091 ## To make an output file (wbbgen_test.out), uncomment the following 7 lines:
092 mod output FileOutput
093 module talk FileOutput
094 output create myout wbbgen_test.out
095 output enable myout
096 output path myout AllPath
097 output keepList myout EVCL_StorableBank LRIH_StorableBank HEPG_StorableBank
098 output list
099 exit
100 ###
101 ## Following gets 10 events; change 10 to 1000 for 1000 events, etc.
102 ## But Note: max. no. of events = maxev*effic, which may be < nev set here
103 begin -nev 10
104 show timer
105 exit
106 #########################
107
Send problems or questions to cdfcode@fnal.gov
