cputTrInfo.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cputtrinfo
subroutine	cwrttrinfo (f, t)
subroutine	ccoordfortr (how, f, t)

Function/Subroutine Documentation

ccoordfortr()

subroutine ccoordfortr	(	integer	how,
		type(coord)	f,
		type(coord)	t
	)

Definition at line 119 of file cputTrInfo.f.

References cerrormsg(), cxyz2det(), and cxyz2prim().

Referenced by chooktrace(), cputcerenkov(), and cputtrinfo().

      implicit none
!     ***************************************
#include  "Ztrack.h"
#include  "Ztrackp.h"
#include  "Ztrackv.h"
#include  "Zobs.h"
#include  "Zobsv.h"
      
      integer how   ! input. see Trace in cputTrInfo.
      type(coord)::f
      type(coord)::t  ! output. transformed 'from to' coord 

      if(how .le. 20) then
!              to 1ry system
         call cxyz2prim(obssites(noofsites)%pos%xyz,
     *                  trackbefmove%pos%xyz, f)
         call cxyz2prim(obssites(noofsites)%pos%xyz,
     *                 movedtrack%pos%xyz, t)
         if(how .ge. 10) then
            f%r(3) = obssites(noofsites)%pos%depth 
     *            - trackbefmove%pos%depth 
            t%r(3) = obssites(noofsites)%pos%depth 
     *            - movedtrack%pos%depth 
         endif
      elseif(how .le. 40) then
!              to detector system
         call cxyz2det(obssites(noofsites)%pos%xyz, 
     *                 trackbefmove%pos%xyz, f)
         call cxyz2det(obssites(noofsites)%pos%xyz, 
     *                 movedtrack%pos%xyz, t)
         if(how .ge. 30) then
            f%r(3) = obssites(noofsites)%pos%depth 
     *            - trackbefmove%pos%depth 
            t%r(3) = obssites(noofsites)%pos%depth 
     *            - movedtrack%pos%depth 
         endif
      elseif(how .le. 60) then
         f = trackbefmove%pos%xyz
         t = movedtrack%pos%xyz
         if(how .ge. 50) then
            f%r(3) = obssites(noofsites)%pos%depth 
     *            - trackbefmove%pos%depth 
            t%r(3) = obssites(noofsites)%pos%depth 
     *            - movedtrack%pos%depth 
         endif
      else
         call cerrormsg('100>= how > 60 for ccoordForTr ', 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cputtrinfo()

subroutine cputtrinfo

(

)

Definition at line 35 of file cputTrInfo.f.

References ccoordfortr(), chooktrace(), cputcerenkov(), and cwrttrinfo().

Referenced by ctracking().

      implicit none

#include  "Ztrack.h"
#include  "Ztrackp.h"
!  #include  "Ztrackv.h"

      
      type(coord)::f
      type(coord)::t

      if(trace .ge. 100 .and. trace .lt. 160) then
         call chooktrace
         return    !      ******************
      endif
      if(trace .gt. 60 .or. trace .gt. 160) then
         call cputcerenkov   ! Cerenkov output
         return    !      *****************
      else
!               convert coord.
         call ccoordfortr(trace, f, t)
         call cwrttrinfo(f, t)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cwrttrinfo()

subroutine cwrttrinfo	(	type(coord)	f,
		type(coord)	t
	)

Definition at line 59 of file cputTrInfo.f.

Referenced by cputtrinfo().

      implicit none

#include  "Ztrack.h"
#include  "Ztrackp.h"
#include  "Ztrackv.h"
#include  "Zobs.h"
#include  "Zobsv.h"
      
      type(coord)::f
      type(coord)::t

      real*8  xxx/-1.d37/, yyy/-1.d36/, zzz/1.d34/
      integer kkk/-1000/, chg/-1000/
      save xxx, yyy, zzz, kkk, chg

      if(kkk .ne. movedtrack%p%code .or. f%r(1) .ne. xxx
     *    .or. f%r(2) .ne. yyy .or. f%r(3) .ne. zzz  .or.
     *     chg .ne. movedtrack%p%charge ) then
         if(xxx .ne. -1.d37) then
            write(tracedev, *) 
            write(tracedev, *) 
         endif
         if(timestructure) then
!!            write(TraceDev, '(3g16.8, i4, g11.4, i4, g16.8)')
            write(tracedev,
     *      '(1p,3g24.16,0p, i4,1p, g24.16,0p, i4, 1p, g24.16)')
     *        f%r(1), f%r(2), f%r(3), trackbefmove%p%code,
     *        trackbefmove%p%fm%p(4)-trackbefmove%p%mass,
     *        trackbefmove%p%charge,  trackbefmove%t
         else
            write(tracedev, '(3g16.8, i4, g11.4, i4)')
     *        f%r(1), f%r(2), f%r(3), trackbefmove%p%code,
     *        trackbefmove%p%fm%p(4)-trackbefmove%p%mass,
     *        trackbefmove%p%charge
         endif
      endif
      if(timestructure) then
!!!         write(TraceDev,  '(3g16.8, i4, g11.4, i4, g16.8)')
            write(tracedev,
     *      '(1p,3g24.16,0p, i4,1p, g24.16,0p, i4, 1p, g24.16)')
     *      t%r(1), t%r(2), t%r(3), movedtrack%p%code,
     *      movedtrack%p%fm%p(4)-movedtrack%p%mass,
     *      movedtrack%p%charge,  movedtrack%t
      else
         write(tracedev,  '(3g16.8, i4, g11.4, i4)')
     *      t%r(1), t%r(2), t%r(3), movedtrack%p%code,
     *      movedtrack%p%fm%p(4)-movedtrack%p%mass,
     *      movedtrack%p%charge
      endif
      xxx = t%r(1)
      yyy = t%r(2)
      zzz = t%r(3)
      kkk = movedtrack%p%code
      chg = movedtrack%p%charge

Here is the caller graph for this function: