B2A104-SimulateAndReconstruct-withoutBeamBkg.py

#!/usr/bin/env python3
 
 
 
 
 
import basf2 as b2
from modularAnalysis import inputMdst
from simulation import add_simulation
from reconstruction import add_reconstruction
from mdst import add_mdst_output
 
# create a path
my_path = b2.create_path()
 
# load input ROOT file
inputMdst(filename=b2.find_file('B2A101-Y4SEventGeneration-evtgen.root'), path=my_path)
 
# simulation
add_simulation(path=my_path)
 
# reconstruction
add_reconstruction(path=my_path)
 
# dump in MDST format
add_mdst_output(path=my_path,
                mc=True,
                filename='B2A101-Y4SEventGeneration-gsim-BKGx0.root')
 
# Show progress of processing
my_path.add_module('ProgressBar')
 
# Process the events
b2.process(my_path)
 
# print out the summary
print(b2.statistics)