######################################
# generate and simulate W + H process: 
# File name wh_pythia.tcl 
#####################################

module input Pythia
talk Pythia

# 
# set up Pythia parameters to run f + fbar ----> WH
#

MSEL set 0

   call_qq set t
   output set HEPG
 qq
   help 
   Decay_Bc set t 
   Decay_Upsilon set t 
   Decay_B_Baryons set t 
   Decay_prompt_charm set t 
   show 
 exit 

commonMenu
# select f + fbar ----> WH production

  set_msub -index=26 -value=1

# turn off decay W ----> jets: it's for decay table 6.129 
  
   for {set i 190} {$i<210} {incr i} { 
    set_mdme -channelIndex=$i -decayType=0   
   } 

# turn on decay W ----> l + nu_l

  set_mdme -channelIndex=206 -decayType=1
  set_mdme -channelIndex=207 -decayType=1
  set_mdme -channelIndex=208 -decayType=1


# set mass of the Higgs boson to 110 GeV

  set_pmas -masscode=25 -mass=110.0

  show_pmas
# turn on decay H ----> b bbar 
  
   for {set i 210} {$i<289} {incr i} { 
    set_mdme -channelIndex=$i -decayType=0 
   } 

  set_mdme -channelIndex=214 -decayType=1

# pdf function MRSG (group 3 set 41) 
  set_mstp -index=51 -value=3041 
  set_mstp -index=52 -value=2

  exit
exit

module enable GEANT3
# CALORIMETRY CT_TrackingModule JETCLU 

talk GEANT3
  GFLASH ON 
  DETECTORS OFF PIPE 
  BANKS HEPG OBSP OBSV 
  BANKS MSVX MISL MVTX ISLD SIXD 
  BANKS COTD MCOT 
  BANKS CMUD CMPD CMXD CMUE CMPE CMXE
  TOWE 1 
  TOWD 1 
  TOWT 0 
  show 
exit 

# configure output
module output FileOutput

talk FileOutput
   output create stream1 pythia_wh_lnubb.ybs
   output enable stream1
   output path stream1 AllPath
   output list
   exit

path list 

show 

begin -nev 10

show timer

exit