PairConversions2.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
 
 
 
"""
This script saves e+ or e- from photon conversions into a pair in MCParticles.
 
<header>
  <contact>dorisykim@ssu.ac.kr</contact>
  <description>
      Saves 100 generic BBbar events with EvtGen + secondary e+ or e- from pair conversions created by Geant4 in MCParticles.
      The corresponding secondaryPhysicsProcess ID is 14, which is defined as fGammaConversion in G4EmProcessSubType.h.
      The detector simulation mixed with background, trigger simulation, and standard reconstruction is done.
  </description>
</header>
"""
 
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 triggr simulation
add_simulation(main, bkgfiles=get_background_files())
 
# saving e+ or e- from pair conversions with kinetic energy > 10.0 MeV.
b2.set_module_parameters(main, "FullSim", StorePairConversions=True, PairConversionsEnergyCut=10.0)
 
# reconstruction
add_reconstruction(main)
 
# output
main.add_module("RootOutput", outputFileName="EvtGenSimRecYesPairConversions.root")
 
 
b2.process(main)
 
# Print call statistics
print(b2.statistics)