cresetMom.f

Go to the documentation of this file.

      subroutine cresetmom(aTrack)
!      Assum the case where the direction cos has been changed
!      and the directionof  vector p should be  changed.
!        (see cresetMombyEDir)
      implicit none

#include  "Ztrack.h"
      type(track)::aTrack  ! input/output.
!
      real*8 p
      call cpxyzp(atrack%p%fm, p)
      atrack%p%fm%p(1) = p * atrack%vec%w%r(1)
      atrack%p%fm%p(2) = p * atrack%vec%w%r(2)
      atrack%p%fm%p(3) = p * atrack%vec%w%r(3)

      end       subroutine cresetmom

      subroutine cresetmombyedir(aTrack)
!          Assume the case where the direction cos is correctly given
!     and E has beeen changed, Vec p should not be used.

!     If E< mass, E is reset to mass. and vec p becomes 0.
      implicit none
#include  "Ztrack.h"
      type(track)::aTrack  ! input/output.
      real(8)::pabs
      pabs = atrack%p%fm%p(4)**2 - atrack%p%mass**2
      if( pabs <= 0. ) then
         atrack%p%fm%p(4) = atrack%p%mass
         pabs = 0.
      else
         pabs = sqrt(pabs)
      endif
      atrack%p%fm%p(1:3) = atrack%vec%w%r(1:3)*pabs
      end      subroutine cresetmombyedir