EvtGenSimNoBkg.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
 
"""
<header>
  <output>EvtGenSimNoBkg.root</output>
  <contact>Software team b2soft@mail.desy.de</contact>
  <cacheable/>
  <description>This steering file produces 1000 generic BBbar events with EvtGen
  and runs the detector simulation without mixing in background.</description>
</header>
"""
 
from basf2 import set_random_seed, create_path, process, statistics, \
    register_module
from simulation import add_simulation
from L1trigger import add_tsim
from beamparameters import add_beamparameters
 
set_random_seed(12345)
 
main = create_path()
 
# specify number of events to be generated
eventinfosetter = register_module('EventInfoSetter')
eventinfosetter.param('evtNumList', [1000])
eventinfosetter.param('runList', [1])
eventinfosetter.param('expList', [0])
main.add_module(eventinfosetter)
 
# beam parameters
beamparameters = add_beamparameters(main, "Y4S")
# beamparameters.param("smearVertex", False)
 
# generate BBbar events
evtgeninput = register_module('EvtGenInput')
main.add_module(evtgeninput)
 
# detector simulation, no background files
add_simulation(main)
 
# trigger simulation
add_tsim(main)
 
# memory profile
main.add_module(register_module('Profile'))
 
# output
output = register_module('RootOutput')
output.param('outputFileName', '../EvtGenSimNoBkg.root')
main.add_module(output)
 
process(main)
 
# Print call statistics
print(statistics)
 
from validation import *
statistics_plots(
    'EvtGenSimNoBkg_statistics.root',
    contact='Software team b2soft@mail.desy.de',
    job_desc='a standard simulation job with generic EvtGen events',
    prefix='EvtGenSimNoBkg'
)
event_timing_plot(
    '../EvtGenSimNoBkg.root', 'EvtGenSimNoBkg_statistics.root',
    contact='Software team b2soft@mail.desy.de',
    job_desc='a standard simulation job with generic EvtGen events',
    prefix='EvtGenSimNoBkg'
)