5 import simulation
as sim
6 from ROOT
import Belle2
8 basf2.set_random_seed(321)
9 basf2.set_log_level(basf2.LogLevel.INFO)
15 module which checks if two collections of EKLMDigits and BKLMDigits are
21 Use some digit information to sort the digits
22 Returns a python-list containing the sorted digits
26 py_list = [x
for x
in unsortedPyStoreArray]
42 load original digits and the packed/unpacked ones, sort and compare them
55 if not len(klm_digits_sorted) == len(klm_digits_unpacked_sorted):
56 basf2.B2FATAL(
"KLMDigits: size not equal after packing and unpacking")
58 for i
in range(len(klm_digits_sorted)):
60 digit = klm_digits_sorted[i]
61 digit_unpacked = klm_digits_unpacked_sorted[i]
64 assert digit.getSubdetector() == digit_unpacked.getSubdetector()
65 assert digit.getSection() == digit_unpacked.getSection()
66 assert digit.getSector() == digit_unpacked.getSector()
67 assert digit.getLayer() == digit_unpacked.getLayer()
68 assert digit.getPlane() == digit_unpacked.getPlane()
69 assert digit.getStrip() == digit_unpacked.getStrip()
70 assert digit.inRPC() == digit_unpacked.inRPC()
71 assert digit.getCharge() == digit_unpacked.getCharge()
72 assert digit.getCTime() == digit_unpacked.getCTime()
73 assert digit.getTDC() == digit_unpacked.getTDC()
74 assert digit.getFitStatus() == digit_unpacked.getFitStatus()
77 main = basf2.create_path()
79 eventinfosetter = basf2.register_module(
'EventInfoSetter')
80 eventinfosetter.param(
'evtNumList', 50)
81 main.add_module(eventinfosetter)
83 particlegun = basf2.register_module(
'ParticleGun')
84 particlegun.param(
'pdgCodes', [13, -13])
85 particlegun.param(
'nTracks', 10)
86 particlegun.param(
'momentumParams', [0.5, 4.0])
87 main.add_module(particlegun)
89 sim.add_simulation(main, components=[
'KLM'])
90 basf2.set_module_parameters(main, type=
'Geometry', useDB=
False, components=[
'KLM'])
92 klm_packer = basf2.register_module(
'KLMPacker')
93 main.add_module(klm_packer)
95 unpacker = basf2.register_module(
'KLMUnpacker')
96 unpacker.param(
'outputKLMDigitsName',
'KLMDigitsUnpacked')
97 main.add_module(unpacker)
101 progress = basf2.register_module(
'Progress')
102 main.add_module(progress)
105 print(basf2.statistics)