001 #####################################################################
002 # Dave Waters, 19.06.2001 #
003 # ======================= #
004 # Example .tcl file for generating and simulating WGRAD events #
005 # using cdfGen #
006 # To use with cdfSim add GeometryManager, SimInitManager and #
007 # SimulationControlMod setup #
008 # #
009 # 02 Oct 2001 lena: updated for cdfGen use #
010 # 25 Jul 2002 dwaters : use standard PythiaMenu #
011 # #
012 #####################################################################
013
014 path enable AllPath
015 creator set NSIM
016
017 module enable GenInputManager
018 mod input GenInputManager
019 talk GenInputManager
020 report set 10
021 first_event set 1
022 show
023 exit
024
025 mod disable Pythia
026 mod disable Isajet
027 mod disable herwig
028 mod disable Bgenerator
029 mod disable FAKE_EVENT
030 mod disable WGRAD
031 mod disable Wbbgen
032 mod disable VecbosModule
033 mod disable VecunwModule
034 mod disable Vecher
035
036 mod disable TauolaModule
037 mod disable QQModule
038 mod disable EvtGen
039
040 module enable GenOutputManager
041
042 mod enable WGRAD
043 mod talk WGRAD
044 debug set f
045 # usePythia = f : WGRAD final state only - lepton, neutrino, (photon).
046 # usePythia = t : WGRAD+PYTHIA - including parton showers and underlying event.
047 usePythia set t
048 cardsMenu
049 # 1: p pbar collider 2: pp collider :
050 BeamType set 1
051 # 1: W+ production, 2: W- production :
052 W_Charge set 1
053 # 1: constant, 2: s-dependent W width :
054 widthTreatment set 2
055 # 1: electron, 2: muon :
056 DecayLepton set 1
057 # 1: initial, 2: final, 4: all :
058 QED_RadType set 2
059 # W mass (GeV) :
060 W_Mass set 80.4
061 # W width (GeV) :
062 W_Width set 2.1
063 # center of mass energy (GeV) :
064 CME set 1960.0
065 # 0:LO only; 1:NLO-LO ; 2:NLO :
066 Order set 2
067 # 0:mu=MW, 1: mu=.5*MW, 2: mu=2*MW. :
068 Scale set 0
069 # final state: 1: full, 2: a la Berends et al. :
070 fsrType set 1
071 # final state: without coll. cut (0), with (1) :
072 fsrCollinearCut set 0
073 # delta_s (photon energy cut) :
074 PhotonEnergyCut set 0.01
075 # delta_c (collinear cut) :
076 PhotonCollinearCut set 0.001
077 # 0: no s.,no r., 1: with s,no r., 2: with s.,with r. :
078 SmearingRecombination set 0
079 # ptmin and ymax for lepton :
080 Lepton_ptmin set 5.0
081 Lepton_ymax set 5.0
082 # ptmin and ymax for neutrino :
083 Neutrino_ptmin set 5.0
084 Neutrino_ymax set 100.0
085 # ptmin and yman for photon :
086 Photon_ptmin set 0.0
087 Photon_ymax set 100.0
088 # low, high transverse mass range :
089 mT_low set 0.0
090 mT_high set 2000.0
091 # pdflib group :
092 pdflibGroup set 3
093 # pdflib choice :
094 pdflibSet set 54
095 # itmax and ncall for grid and precision calculation :
096 itmaxGrid set 1
097 ncallGrid set 5000
098 itmaxPrecision set 2
099 ncallPrecision set 10000
100 # Use "RandomNumberMenu" below.
101 # random number seed :
102 # randomSeed set 8910
103 exit
104 PythiaMenu
105 commonMenu
106 # Fragmentation and decay : 0 - off, 1 - on (default) :
107 # -----------------------------------------------------
108 # set_mstp -index=111 -value=1
109 #
110 # QCD and QED showers. QED showers must be turned off
111 # if WGRAD is run in NLO mode, to avoid double counting.
112 # 0 - no showers, 1 - QCD only, 2 - QCD+QED (default) :
113 # -------------------------------------------------------
114 set_mstj -index=41 -value=1
115 #
116 # Master switch for initial state showering :
117 # 0 - no showering, 1 - QCD only, 2 - QCD+QED (default) :
118 # -------------------------------------------------------
119 set_mstp -index=61 -value=1
120 #
121 # Master switch for final state showering :
122 # 0 - off, 1 - on (default) :
123 # ---------------------------
124 set_mstp -index=71 -value=1
125 #
126 # Random number seed :
127 # --------------------
128 # set_mrpy -index=1 -value=280273
129 # Multiple interactions : 0 - off, 1 - on (default) :
130 # ---------------------------------------------------
131 set_mstp -index=81 -value=0
132 #
133 # Intrinsic kT : 0 - none, 1 - Gaussian (default), 2 - exponential :
134 # ------------------------------------------------------------------
135 # set_mstp -index=91 -value=1
136 #
137 # kT Gaussian width (GeV) :
138 # -------------------------
139 # set_parp -index=91 -value=4.0
140 #
141 # kT exponential width (GeV) :
142 # ----------------------------
143 # set_parp -index=92 -value=1.63
144 #
145 # Upper bound on intrinsic kT (GeV) :
146 # -----------------------------------
147 # set_parp -index=93 -value=20.0
148 exit
149 exit
150 RandomNumberMenu
151 RandomSeed1 set 1234
152 RandomSeed2 set 5678
153 exit
154 exit
155
156 mod output FileOutput
157 talk FileOutput
158 output create main_stream cdfSim.root
159 output path main_stream AllPath
160 output keepList main_stream \
161 LRIH_StorableBank EVCL_StorableBank HEPG_StorableBank
162 exit
163
164 module enable HepRootManager
165 mod talk HepRootManager
166 verbose set true
167 histfile set wgrad.hist
168 createHistoFile set f
169 exit
170
171 show
172
173 path list
174
175 action on "Name Action"
176
177 begin -nev 10
178 show timer
179
180 exit
181
182
183
184
Send problems or questions to cdfcode@fnal.gov
