63 reimplement Module::event()
68 track = pid.getRelatedFrom(
'Tracks')
70 mcpart = track.getRelatedTo(
'MCParticles')
71 if not track
or not mcpart:
73 b2.B2WARNING(
'problems with track <-> mcparticle relations')
75 print(f
'event: {int(event)}, track: {int(track.getArrayIndex())}')
77 pdg = abs(mcpart.getPDG())
78 momentumVec = mcpart.getMomentum()
79 momentum = momentumVec.Mag()
85 selectedpart = Belle2.Const.kaon
88 logl_sel = pid.getLogL(selectedpart, pid_dedx)
89 logl_pi = pid.getLogL(Belle2.Const.pion, pid_dedx)
90 dedx_DLL = logl_pi - logl_sel
92 logl_sel = pid.getLogL(selectedpart, pid_top)
93 logl_pi = pid.getLogL(Belle2.Const.pion, pid_top)
94 top_DLL = logl_pi - logl_sel
96 if pdg == selectedpart.getPDGCode():
97 hist[0].Fill(momentum, dedx_DLL)
98 hist[2].Fill(momentum, top_DLL)
100 hist[1].Fill(momentum, dedx_DLL)
101 hist[3].Fill(momentum, top_DLL)