###############################################################################
#  this example shows how to generate and simulate in the detector 
#  ppbar --> ttbar process
#
#  CDF2SIM defaults: - path AllPath is enabled
#                    - all the modules are disabled 
#  so one needs to enable the modules he needs
#
# to run cdf2sim in batch mode one needs to add "-b" flag to disable 
# initialization of ROOT graphics :
#                                       cdf2sim.exe input_file.tcl -b
###############################################################################
#
# Select the Generator as Pythia
module input Pythia
#                                   
# Setup Pythia                                
talk Pythia
#

# Make sure the center of mass energy is 1800
   CMEnergy set 2000
   
#
# Set the process  (See Section 9.2, Page 144 in Pythia Manual for MSEL values)
  MSEL set 6  
  commonMenu
#
# Set pt range for 2 -> 2 process (e.g. of two top's)
# See Section 9.2 Page 145 for ckin variable description)
    set_ckin -index=3  -value=5.
    set_ckin -index=4  -value=750.
#    
# mstp Parameter description can be found
#     Section 9.3 Pages 149-162
#                 Pages 210-211
#          Pages 225-227
#
#
#  Sets the top decay before hadronization with W polarization effects                                      
    set_mstp -index=48 -value=2
#    
# Initial state radiation turned on (value=1) or off (value = 0)
    set_mstp -index=61 -value=1

# Final state radiation turned on (value=1) or off(value=0)    
    set_mstp -index=71 -value=1
#    
# Setup the structure functions.                                      
# Tell it to use the pdflib (value=2)
    set_mstp -index=52 -value=2

# Select CTEQ 4L (value=4032) from the pdflib    
    set_mstp -index=51 -value=4032
    
# 
# Calculation of alpha_s at hard interaction first-order(value=1) 
#                         second-order(value=2) 
    set_mstp -index=2 -value=2
#
# K factors in hard cross section: separate K factors for 
#                   color annil and ordinary diag (value=2)    
    set_mstp -index=33 -value=2
#
# Setup Multiple interactions: (see page 225)    
    set_mstp -index=82 -value=2
    set_parp -index=82 -value=1.

# Fragmentaion flag for generic quarks (see Page 269)
    set_mstj -index=11 -value=3
# Fragmentation flag for b quarks from top (See Page 270-271)    
    set_mstj -index=27 -value=2
# Fragmentation parameter  (see page 273)  
    set_parj -index=55 -value=-0.006
    
#    
# Set_pmas command takes just one index ?
    set_pmas -masscode=6  -mass=175
    
# Setup the decay of the W's (Turn off W- decays to quarks)
#
# Note: If you are reading in a decay table above, these
# changes below do not take effect.
#
# First look at the quark decays.   (decayType=0: Channel off )
#                    (decayType=1: Channel on  )
#   (see page 281)            (decayType=2: Channel on for particle/Off for antipart)
#                    (decayType=3: Channel off for particle/on for antipart) 
# W -> ud
  set_mdme -channelIndex=190 -decayType=2
# W -> cd
  set_mdme -channelIndex=191 -decayType=2
# W -> td  (BR = 0.)
  set_mdme -channelIndex=192 -decayType=2
# W -> us
  set_mdme -channelIndex=194 -decayType=2
# W -> cs
  set_mdme -channelIndex=195 -decayType=2
# W -> ts  (BR = 0.)
  set_mdme -channelIndex=196 -decayType=2
# W -> ub
  set_mdme -channelIndex=198 -decayType=2
# W -> cb 
  set_mdme -channelIndex=199 -decayType=2
# W -> tb  (BR = 0.)
  set_mdme -channelIndex=200 -decayType=2
  
# turn on decay W ----> l + nu_l  (Both W+ and W- can go leptonically)
# W -> e nu_e
  set_mdme -channelIndex=206 -decayType=1
# W -> mu nu_mu
  set_mdme -channelIndex=207 -decayType=1
# W -> tau nu_tau
  set_mdme -channelIndex=208 -decayType=1

# Setup the decay of the Z 
# First look at the quark decays.   (decayType=0: Channel off )
#                    (decayType=1: Channel on )
# Z -> d dbar
  set_mdme -channelIndex=174 -decayType=1
# Z -> u ubar
  set_mdme -channelIndex=175 -decayType=1
# Z -> s sbar
  set_mdme -channelIndex=176 -decayType=1
# Z -> c cbar
  set_mdme -channelIndex=177 -decayType=1
# Z -> b bbar
  set_mdme -channelIndex=178 -decayType=1
# Z -> t tbar  (BR = 0.)
  set_mdme -channelIndex=179 -decayType=1
# Z -> e e
  set_mdme -channelIndex=182 -decayType=1
# Z -> nu_e nu_e
  set_mdme -channelIndex=183 -decayType=1
# Z -> mu mu
  set_mdme -channelIndex=184 -decayType=1
# Z -> nu_mu nu_mu
  set_mdme -channelIndex=185 -decayType=1
# Z -> tau tau
  set_mdme -channelIndex=186 -decayType=1
# Z -> nu_tau nu_tau
  set_mdme -channelIndex=187 -decayType=1
  
  exit

# Write out the decay table. This will have the changes made above
# (No real need to do this...other than checking the changes above)
#  decayFileMode set write
#  decayFile set decay.dbt
    
exit
#
#  Setup GEANT3
module enable GEANT3 
talk GEANT3
   DETECTORS OFF PIPE
   source $env(CDFSOFT2_DIR)/RootMods/examples/cdf2sim/tcl/geant3_output_banks.tcl
exit
   
#  Disable plotting package
module disable g3hist

# Show the path
path list
#
# Set up the output
module output FileOutput
talk FileOutput
  output create AA ./results/ttbar_gensim.rt
  output path AA AllPath
  source $env(CDFSOFT2_DIR)/RootMods/examples/cdf2sim/tcl/cdf2sim_AA_keepList.tcl
  output list
  exit
#
# Run this bad boy... 
begin -nev 100

# Show the statistics
show
exit