saut = 1
pos = f[::saut, 3139] #position
sig_part = f[::saut, 0:1999] #signal
ax_part = f[::saut, 3134] #theta
d_part = f[::saut, 3133] #rayon
hydro = f[::saut, 3143]