a shorter notation:
(p+len/2).x(), (p+len/2).y(), (p+len/2).z(), p_i.charge, p_i.radius