B2A103-SimulateAndReconstruct.py

#!/usr/bin/env python3
 
 
 
 
 
import basf2 as b2
from mdst import add_mdst_output
from simulation import add_simulation
from reconstruction import add_reconstruction
from modularAnalysis import inputMdst
 
# We previously generated events but we want to simulate them with the
# currently best values, not whatever was valid back then. So we have to
# disable reuse of the same globaltag from when the events were generated
b2.conditions.disable_globaltag_replay()
 
# create path
my_path = b2.create_path()
 
# load input ROOT file
inputMdst(filename=b2.find_file('B2pi0D_D2hh_D2hhh_B2munu_evtgen.root', 'examples', False),
          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='B2pi0D_D2hh_D2hhh_B2munu.root')
 
# Show progress of processing
my_path.add_module('ProgressBar')
 
# Process the events
b2.process(my_path)
 
# print out the summary
print(b2.statistics)