001 #######################################################################
002 # Dave Waters, 29.6.2001 #
003 # ====================== #
004 # tcl file for generating and simulating di-boson events using PYTHIA #
005 #######################################################################
006 mod enable Pythia
007 mod talk Pythia
008 PythiaMenu
009 msel set 0
010 cmEnergy set 1960
011 commonMenu
012 ########################################################
013 # Set the process code. (PYTHIA manual, pp.113-116)
014 # =================================================
015 # WW :
016 # ----
017 set_msub -index=25 -value=1
018 # WZ :
019 # ----
020 # set_msub -index=23 -value=1
021 # ZZ :
022 # ----
023 # set_msub -index=22 -value=1
024 ########################################################
025 # Set the mass range for the di-bosons :
026 # ======================================
027 set_ckin -index=41 -value=0.0
028 set_ckin -index=42 -value=-1
029 set_ckin -index=43 -value=0.0
030 set_ckin -index=44 -value=-1
031 ########################################################
032 #
033 # s-hat is the scale for the hard process
034 #
035 set_mstp -index=32 -value=4
036 #
037 ########################################################
038 # Set the following to avoid unnecessary tasks
039 # (for example, for generator level studies) :
040 # ============================================
041 # No initial-state QCD or QED radiation :
042 ## set_mstp -index=61 -value=0
043 # No final-state QCD or QED radiation :
044 ## set_mstp -index=71 -value=0
045 # No multiple interactions :
046 ## set_mstp -index=81 -value=0
047 # Switch off fragmentation and decay
048 # (might well crash the detector simulation, so you
049 # probably want to disable SimulationControlMod
050 # in this case).
051 ## set_mstp -index=111 -value=0
052 ########################################################
053 # Select certain decy modes only
054 # ==============================
055 #
056 # Simply comment out all mdme settings to generate
057 # events with the full compliment of Standard Model
058 # decay modes. Note that PYTHIA still uses full
059 # Standard Model widths, even in the absence of
060 # certain decay channels.
061 #
062 # Decay table for Z0 :
063 # --------------------
064 # 174 1 32 0.154211 d dbar
065 # 175 1 32 0.119589 u ubar
066 # 176 1 32 0.154208 s sbar
067 # 177 1 32 0.119452 c cbar
068 # 178 1 32 0.152331 b bbar
069 # 179 1 32 0.000000 t tbar
070 # 180 -1 32 0.000000 b' b'bar
071 # 181 -1 32 0.000000 t' t'bar
072 # 182 1 0 0.033480 e- e+
073 # 183 1 0 0.066615 nu_e nu_ebar
074 # 184 1 0 0.033480 mu- mu+
075 # 185 1 0 0.066615 nu_mu nu_mubar
076 # 186 1 0 0.033404 tau- tau+
077 # 187 1 0 0.066615 nu_tau nu_taubar
078 # 188 -1 0 0.000000 tau'- tau'+
079 # 189 -1 0 0.000000 nu'_tau nu'_taubar
080 #
081 # Turn off all SM decays :
082 # ========================
083 ## set_mdme -channelIndex=174 -decayType=0
084 ## set_mdme -channelIndex=175 -decayType=0
085 ## set_mdme -channelIndex=176 -decayType=0
086 ## set_mdme -channelIndex=177 -decayType=0
087 ## set_mdme -channelIndex=178 -decayType=0
088 ## set_mdme -channelIndex=179 -decayType=0
089 ## set_mdme -channelIndex=182 -decayType=0
090 ## set_mdme -channelIndex=183 -decayType=0
091 ## set_mdme -channelIndex=184 -decayType=0
092 ## set_mdme -channelIndex=185 -decayType=0
093 ## set_mdme -channelIndex=186 -decayType=0
094 ## set_mdme -channelIndex=187 -decayType=0
095 # Enable Z->e+e- only :
096 # =====================
097 ## set_mdme -channelIndex=182 -decayType=1
098
099 # Decay table for W+/- :
100 # ----------------------
101 # 190 1 32 0.321822 dbar u
102 # 191 1 32 0.016518 dbar c
103 # 192 1 32 0.000000 dbar t
104 # 193 -1 32 0.000000 dbar t'
105 # 194 1 32 0.016525 sbar u
106 # 195 1 32 0.321097 sbar c
107 # 196 1 32 0.000000 sbar t
108 # 197 -1 32 0.000000 sbar t'
109 # 198 1 32 0.000010 bbar u
110 # 199 1 32 0.000592 bbar c
111 # 200 1 32 0.000000 bbar t
112 # 201 -1 32 0.000000 bbar t'
113 # 202 -1 32 0.000000 b'bar u
114 # 203 -1 32 0.000000 b'bar c
115 # 204 -1 32 0.000000 b'bar t
116 # 205 -1 32 0.000000 b'bar t'
117 # 206 1 0 0.107838 e+ nu_e
118 # 207 1 0 0.107838 mu+ nu_mu
119 # 208 1 0 0.107759 tau+ nu_tau
120 # 209 -1 0 0.000000 tau'+ nu'_tau
121 #
122 # Turn off all SM decays :
123 # ========================
124 set_mdme -channelIndex=190 -decayType=0
125 set_mdme -channelIndex=191 -decayType=0
126 set_mdme -channelIndex=192 -decayType=0
127 set_mdme -channelIndex=194 -decayType=0
128 set_mdme -channelIndex=195 -decayType=0
129 set_mdme -channelIndex=196 -decayType=0
130 set_mdme -channelIndex=198 -decayType=0
131 set_mdme -channelIndex=199 -decayType=0
132 set_mdme -channelIndex=200 -decayType=0
133 set_mdme -channelIndex=206 -decayType=1
134 set_mdme -channelIndex=207 -decayType=1
135 set_mdme -channelIndex=208 -decayType=1
136 # Enable W->ev only :
137 # ===================
138 ## set_mdme -channelIndex=206 -decayType=1
139 ########################################################
140 exit
141 exit
142 exit
143 source $env(PROJECT_DIR)/mcProduction/tcl/mc_Pythia_pdf_CTEQ5L.tcl
144 source $env(PROJECT_DIR)/mcProduction/tcl/mc_Pythia_underlying_event_A.tcl
Send problems or questions to cdfcode@fnal.gov
