small improvements: replace pow() by multiplications, store particle mass in addition to the squared mass in order to avoid sqrt(mass2)