let
mass p = (p.mz -. mass_proton) *. p.z