twiss parameter calculation added
simplified down to 1 type of particles