Added phase space calculation script