# Noise_Ev=f["RealRes_Noise_Layer"+str(lay+1)][ev][0]
Lay_Ev = Lay_Ev + Noise_Ev
Topo_Ev = gT["TopoClusters" +
str(lay + 1)][ev][0][:layers[lay], :layers[lay]]
for X in range(layers[lay]):
for Y in range(layers[lay]):
if Topo_Ev[X][Y] != 0:
total.append([
lay, X, Y, Lay_Ev[X][Y], Topo_Ev[X][Y],
np.sqrt((X - Traj_Mark[lay][0]) *
(X - Traj_Mark[lay][0]) +
(Y - Traj_Mark[lay][1]) *
(Y - Traj_Mark[lay][1]))
])
if total == []: total.append([0, 0, 0, 0, 0])
LFI = Track.ComputeLFI(total, Traj_Mark)
max_value = max(total, key=itemgetter(4))
dist = 0
mx = int(max_value[4])
en = []
closest = []
for m in range(1, mx + 1):
average_X = 0
average_Y = 0
energy = 0
tot = 0
totA = 0
for cell in total:
if m != cell[4]: continue
if cell[0] == LFI:
average_X = average_X + cell[1]
