best_jmuon() and loop over best_tracks takes unexpected amount of time
The best_jmuon() takes 20 minutes to process one run but this was already reported. Besides it seems that the more i (loop counter) increases the more the processing time for each iteration raises. The script works retrieving relevant parameters for 2000 events in like 5 minutes, but for 10 000 we reach the CPU time limit.
/sps/km3net/users/lcerisy/sun_orca6.py (the script)
/sps/km3net/users/lcerisy/datav6_ORCA_test.jchain.aanet.00009257.root (for instance)
`(...)
best_tracks = ki.tools.best_jmuon(tracks) // --------------------------------------------> takes like 20 minutes.
indices = np.arange(len(best_tracks))
with open("/sps/km3net/users/lcerisy/output_sun_orca6.dat", "a") as output:
for i in indices: //-------------------------------------------------------------> takes hours.
tr = best_tracks[i]
v = (tr.dir_x,tr.dir_y,tr.dir_z)
if type(tr.dir_x) != float or type(tr.dir_y) != float or type(tr.dir_z) != float:
continue
az = azimuth(v) #radian!!!
zen = zenith(v) #radian!!!
time = np.datetime64(f.events[i].t_sec,"s")
dist_to_sun = orca_sun_dist1(az,time,zen,loc)(....)`