cgetZenith.f

Go to the documentation of this file.

#include "ZsubstRec.h"
!    
      subroutine csetdircos(dc, aTrack)
      implicit none

#include  "Ztrack.h"
      type(track)::aTrack
      type(coord)::dc
!
#ifdef SUBSTREC      
      atrack%vec%w = dc
#else
      atrack%vec%w%r = dc%r
#endif
      call cgetzenith(atrack, atrack%vec%coszenith)
      end
!         get cos of zenith angle from track information
!
      subroutine cgetzenith(aTrack, cosz)
      implicit none

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

!
      cosz = - (atrack%vec%w%r(1) * atrack%pos%xyz%r(1) + 
     *          atrack%vec%w%r(2) * atrack%pos%xyz%r(2) + 
     *          atrack%vec%w%r(3) * atrack%pos%xyz%r(3))/
     *          atrack%pos%radiallen
      if(cosz .gt. 1.d0) then
         cosz = 1.d0
      elseif(cosz .lt. -1.d0) then
         cosz = -1.d0
      endif
           ! or
!     *     sqrt(aTrack.pos.xyz.r(1)**2 +
!     *          aTrack.pos.xyz.r(2)**2 +
!     *          aTrack.pos.xyz.r(3)**2)
      end