cresetDirec.f

Go to the documentation of this file.

!        call this when momentum has been changed to reset the
!     direction information.
!
      subroutine cresetdirec(aTrack)
      implicit none

#include  "Ztrack.h"
      type(track)::aTrack
!
      real*8 pabs

      call cpxyzp(atrack%p%fm, pabs)
      if(pabs .gt.  0) then
         atrack%vec%w%r(1) = atrack%p%fm%p(1)/pabs
         atrack%vec%w%r(2) = atrack%p%fm%p(2)/pabs
         atrack%vec%w%r(3) = atrack%p%fm%p(3)/pabs
      else
         atrack%vec%w%r(1) = 0.
         atrack%vec%w%r(2) = 0.
         atrack%vec%w%r(3) = 1.
      endif
      call cgetzenith(atrack, atrack%vec%coszenith)
      end