######################################
# generate and simulate Z -> tautau process: 
# File name ztautau_pythia.tcl 
#####################################

module input Pythia
talk Pythia

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

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 ----> Z/gamma production (zero order)

  set_msub -index=1 -value=1

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

# turn on decay Z ----> tau+ + tau-

  set_mdme -channelIndex=186 -decayType=1

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

# select only z production 
  set_mstp -index=43 -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_z_tautau.ybs
   output enable stream1
   output path stream1 AllPath
   output list
   exit

path list 

show 

begin -nev 10

show timer

exit