B2A102-ccbarEventGeneration.py

#!/usr/bin/env python3
 
 
 
 
 
import basf2 as b2
from modularAnalysis import setupEventInfo
from mdst import add_mdst_output
from generators import add_inclusive_continuum_generator
 
# Defining one path
my_path = b2.create_path()
 
# generation of 100 events according to the specified DECAY table
# e+e- -> ccbar -> D*+ anything
# D*+ -> D0 pi+; D0 -> K- pi+
#
setupEventInfo(100, path=my_path)
add_inclusive_continuum_generator(finalstate="ccbar",
                                  particles=["D*+"],
                                  userdecfile=b2.find_file('analysis/examples/simulations/B2A102-ccbarEventGeneration.dec'),
                                  path=my_path)
 
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)