B2A102-ccbarEventGeneration.py

#!/usr/bin/env python3
 
 
 
import basf2 as b2
import modularAnalysis as ma
import reconstruction as re
import generators as ge
 
# Defining one path
my_path = b2.create_path()
 
# generation of 100 events according to the specified DECAY table
# e+e- -> ccbar -> D*+ anthing
# D*+ -> D0 pi+; D0 -> K- pi+
#
ma.setupEventInfo(100, path=my_path)
ge.add_inclusive_continuum_generator(finalstate="ccbar",
                                     particles=["D*+"],
                                     userdecfile=b2.find_file('analysis/examples/simulations/B2A102-ccbarEventGeneration.dec'),
                                     path=my_path)
 
re.add_mdst_output(mc=True, filename='B2A102-ccbarEventGeneration-kkmc.root', path=my_path)
 
# Process the events
b2.process(my_path)
 
# print out the summary
print(b2.statistics)