      real function amass(p)

      implicit none

      integer ndim
      parameter (ndim=4)

      real p(ndim)

      amass = sqrt(p(4)**2-p(1)**2-p(2)**2-p(3)**2)

      return
      end
