############################################################################
# this .tcl file generates W-gluon single-top events using vesrion 6.129   #
# of Pythia.  It is a modified version of the Drell-Yan tcl file (by J.    #
# Nachmann and J. Goldstein) which used an example .tcl file taken         #
# from the OSU web site which can be found at:                             #
#                                                                          #
#     http://www.physics.ohio-state.edu/~cdf/Offline/CDFOffline.html       #
#                                                                          #
# to run CDF2SIM on it:                                                    #
#                       cdf2sim.exe  filename.tcl                          #
#                                                                          #
# to run cdf2sim in batch mode one needs to add "-b" flag to disable       #
# initialization of ROOT graphics                                          #
#                                                                          #
#                       cdf2sim.exe   filename.tcl -b                      #
#                                                                          #
#############################                                              #
# Author: Pierre Savard     #                                              #
############################################################################ 
# configure input module

module input Pythia
talk Pythia
  run_number set $run
  report     set $report_freq
  rnIoMode   set READ_THEN_SAVE
  rnSaveFreq set 1

#  Listlast set 100

# select complete user control of processes

  MSEL set 36
  commonMenu

# turn on desired processes (for now a partial list, -start with only the
# [2->1, tree] process, [ffbar->gamma*/Z], and later (?) add [2->2, tree]
# processes and (?) [2->1, one-loop] processes, etc...)
#

    set_msub -index=83 -value=1    

# set pt range --NOTE: not clear to me what setting ckin(3) = 10.
# will actually do since we are selecting a 2->1 tree process and all of 
# the ckin() explanations seem to refer to 2->2 and 2->1->2 processes, but
# I will set it this way anyway.....

    set_ckin -index=3  -value=10.
    set_ckin -index=4  -value=750.


# also set eta ranges
    set_ckin -index=13  -value=-4.
    set_ckin -index=14  -value=4.
    set_ckin -index=15  -value=-4.
    set_ckin -index=16  -value=4.


# implement fix (recommended by Pasha M./Stan Thompson) to fix underlying
# event parameters regarding fragmentation

    set_parp -index=81 -value=1.4
    set_parp -index=82 -value=1.55
    set_parp -index=90 -value=0.

# set initial/final state radiation

    set_mstp -index=61 -value=1
    set_mstp -index=71 -value=1

# NOTE: mstp(43) has default=3 for sub-process 1, so don't need to
# set it to 3 in this file...I will do it anyway since later other 
# sub-processes may be included which need this to be set...

    set_mstp -index=43 -value=3

# choose structure functions

    set_mstp -index=52 -value=1
    set_mstp -index=51 -value=10

# tune underlying event parameters

    set_mstp -index=2 -value=2
    set_mstp -index=33 -value=2
    set_mstp -index=82 -value=2
#    set_parp -index=82 -value=1. --have to comment this out bec. of fix above

    set_mstj -index=11 -value=3
    set_mstj -index=27 -value=1
    set_parj -index=55 -value=-0.006

# specify top mass [set_pmas command takes just one index ??]

    set_pmas -masscode=6 -mass=175.
   
# **note that the following lines reference the updated decay table for
#   version 6.129 of Pythia; it is worth checking this regularly
#   since the decay table changes from time to time**

# NOTE: the W stuff is not needed, but I don't see any harm in
# leaving the settings in anyway...comment them out or do what you will
# setup the decay of W's



    set_mdme -channelIndex=190 -decayType=0
    set_mdme -channelIndex=191 -decayType=0
    set_mdme -channelIndex=192 -decayType=0
    set_mdme -channelIndex=193 -decayType=0
    set_mdme -channelIndex=194 -decayType=0
    set_mdme -channelIndex=195 -decayType=0
    set_mdme -channelIndex=196 -decayType=0
    set_mdme -channelIndex=197 -decayType=0
    set_mdme -channelIndex=198 -decayType=0
    set_mdme -channelIndex=199 -decayType=0
    set_mdme -channelIndex=200 -decayType=0
    set_mdme -channelIndex=201 -decayType=0
    set_mdme -channelIndex=202 -decayType=0
    set_mdme -channelIndex=203 -decayType=0
    set_mdme -channelIndex=204 -decayType=0
    set_mdme -channelIndex=205 -decayType=0  
#
    set_mdme -channelIndex=206 -decayType=1
    set_mdme -channelIndex=207 -decayType=1
    set_mdme -channelIndex=208 -decayType=1
    set_mdme -channelIndex=209 -decayType=0

 
# setup the decay of gammas (all quark decays are turned off...)

    set_mdme -channelIndex=162 -decayType=0    
    set_mdme -channelIndex=163 -decayType=0
    set_mdme -channelIndex=164 -decayType=0
    set_mdme -channelIndex=165 -decayType=0
    set_mdme -channelIndex=166 -decayType=0
    set_mdme -channelIndex=167 -decayType=0
    set_mdme -channelIndex=168 -decayType=-1
    set_mdme -channelIndex=169 -decayType=-1
    set_mdme -channelIndex=170 -decayType=1
    set_mdme -channelIndex=171 -decayType=1
    set_mdme -channelIndex=172 -decayType=1
    set_mdme -channelIndex=173 -decayType=-1

    exit
  exit
#-----------------------------------------------------------------------
#  end of the singletop_wgluon.tcl
#-----------------------------------------------------------------------