SecondaryParticles2.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
 
 
 
"""
This script saves all the secondary particles in MCParticles.
 
Saves 100 generic BBbar events with EvtGen + all the secondary particles created by Geant4 in MCParticles.
The detector simulation mixed with background, trigger simulation, and standard reconstruction is done.
"""
 
import basf2 as b2
from simulation import add_simulation
from reconstruction import add_reconstruction
from background import get_background_files
 
b2.set_random_seed(12345)
 
 
main = b2.create_path()
 
# specify number of events to be generated
main.add_module('EventInfoSetter', evtNumList=[100], runList=[1], expList=[1])
 
# generate BBbar events
main.add_module('EvtGenInput')
 
# detector and L1 trigger simulation
add_simulation(main, bkgfiles=get_background_files())
 
# saving secondary particles with kinetic energy > 1.0 MeV.
b2.set_module_parameters(main, "FullSim", StoreAllSecondaries=True, SecondariesEnergyCut=1.0)
 
# reconstruction
add_reconstruction(main)
 
# output
main.add_module("RootOutput", outputFileName="EvtGenSimRecYesSecondaries.root")
 
 
b2.process(main)
 
# Print call statistics
print(b2.statistics)