Atlfast_RunPythia.py

Go to the documentation of this file.

include( "PartPropSvc/PartPropSvc.py" )
theApp.Dlls += [ "TruthExamples", "Pythia_i" ]
theApp.ExtSvc += ["AtRndmGenSvc" ]

AtRndmGenSvc = Service( "AtRndmGenSvc" )
AtRndmGenSvc.Seeds = [ "PYTHIA 85909879 9707499", "PYTHIA_INIT 87928165 11160529" ]

#theApp.TopAlg += ["Pythia"]


#======================================
#
#  Define the physics process
#  Pass Pythia configuration options
#  through Pythia.PythiaCommand
#
#======================================

Pythia = Algorithm( "Pythia" )

# Switch off all processes
Pythia.PythiaCommand += [ "pysubs msel 0" ]

# Turn on Z0 production
Pythia.PythiaCommand += [ "pysubs msub 1 1" ]

# Turn off Z0/gamma* interference
Pythia.PythiaCommand += [ "pypars mstp 43 2" ]

# Switch off all Z0 decay modes except electrons
Pythia.PythiaCommand += [ "pydat3 mdme 174 1 0", "pydat3 mdme 175 1 0",
                          "pydat3 mdme 176 1 0", "pydat3 mdme 177 1 0",
                          "pydat3 mdme 178 1 0", "pydat3 mdme 179 1 0",
                          "pydat3 mdme 180 1 0", "pydat3 mdme 181 1 0",
                          "pydat3 mdme 182 1 1", "pydat3 mdme 183 1 0",
                          "pydat3 mdme 184 1 0", "pydat3 mdme 185 1 0",
                          "pydat3 mdme 186 1 0", "pydat3 mdme 187 1 0" ]

topSequence += Pythia

---