######################################
# generate and simulate Z + H -> nunu bb process: 
# File name zh_pythia.tcl 
#####################################

module input Pythia
talk Pythia

  run_number set $run
  report     set $report_freq
  rnIoMode   set READ_THEN_SAVE
  rnSaveFreq set 1
  call_qq    set t
# 
# set up Pythia parameters to run f + fbar ----> ZH
#

 MSEL set 0

 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 ----> ZH production

  set_msub -index=24 -value=1

# turn off  Z0 decay
  
   for {set i 174} {$i<190} {incr i} { 
    set_mdme -channelIndex=$i -decayType=0
   } 

# turn on decay Z0 ----> nu + nu_l

  set_mdme -channelIndex=183 -decayType=1
  set_mdme -channelIndex=185 -decayType=1
  set_mdme -channelIndex=187 -decayType=1


# set mass of the Higgs boson to 110 GeV

  set_pmas -masscode=25 -value=110.0

# 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 g3hist

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_ZH_nunubb.ybs
   output enable stream1
   output path stream1 AllPath
   output list
   exit

path list 

# define ROOT histogram file
#module enable g3hist

#talk g3hist
#   histfile pythia_WH_nunubb.root
#   exit

# start processing the events

show
begin -nev 10

show timer 

exit