Belle II Software light-2406-ragdoll
1_gensimrec_1110021010.py
1#!/usr/bin/env python3
2
3
10
11"""
12<header>
13 <output>../1110021010.dst.root</output>
14 <contact>Frank Meier; frank.meier@duke.edu</contact>
15</header>
16"""
17# Generates and reconstructs 1000 [B0 -> rho0 gamma] decay events for nightly
18# validation. This script is NOT AN OFFICIAL DATA PRODUCTION script. Always
19# always refer to the B2P project for the latest official scripts.
20
21import basf2
22from generators import add_evtgen_generator
23from simulation import add_simulation
24from reconstruction import add_reconstruction
25from mdst import add_mdst_output
26
27DECAY_MODE_ID = 1110021010
28OUTPUT_FILENAME = f"../{DECAY_MODE_ID}.dst.root"
29
30basf2.set_random_seed(12345) # for reproducibility
31main = basf2.Path()
32main.add_module('EventInfoSetter', evtNumList=[1000])
33add_evtgen_generator(
34 main, 'signal',
35 basf2.find_file(f'decfiles/dec/{DECAY_MODE_ID}.dec')
36) # signal event generation
37add_simulation(main) # detector simulation without background overlay
38add_reconstruction(main) # reconstruction
39add_mdst_output(main, filename=OUTPUT_FILENAME) # mdst output
40main.add_module('Progress')
41basf2.process(main)
42print(basf2.statistics)