EvtGenSimNoBkg.py

#!/usr/bin/env python3
 

 
"""
<header>
  <output>EvtGenSimNoBkg.root</output>
  <contact>arul.prakash@physik.uni-muenchen.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,
    register_module,
)
from simulation import add_simulation
from beamparameters import add_beamparameters
from validation import statistics_plots, event_timing_plot
 
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 and L1 trigger simulation, no background files
add_simulation(main)
 
# memory profile
main.add_module(register_module("Profile"))
 
# output
output = register_module("RootOutput")
output.param("outputFileName", "../EvtGenSimNoBkg.root")
main.add_module(output)
 
main.add_module('Progress')
process(main, calculateStatistics=True)
 
statistics_plots(
    "EvtGenSimNoBkg_statistics.root",
    contact="arul.prakash@physik.uni-muenchen.de",
    job_desc="a standard simulation job with generic EvtGen events",
    prefix="EvtGenSimNoBkg",
)
event_timing_plot(
    "../EvtGenSimNoBkg.root",
    "EvtGenSimNoBkg_statistics.root",
    contact="arul.prakash@physik.uni-muenchen.de",
    job_desc="a standard simulation job with generic EvtGen events",
    prefix="EvtGenSimNoBkg",
)