det2prim.f

Go to the documentation of this file.

!
!     get transformation matrix for converting detector system vector
!     into primary system vector.
!    Since primary system depends on each sampled primary, the 
!    random number by param file must be fixed.
!    This program gets the conversion matrix for the 1 st event
!    specified in the param
!

#include  "BlockData/cblkGene.h"
      implicit none

#include  "Zglobalc.h"
#include  "Ztrack.h"
#include  "Ztrackp.h"
#include  "Ztrackv.h"
#include  "Zobs.h"
#include  "Zobsv.h"




      type(coord):: aa, bb, cc, det, detyaxis
      integer fin
      type(track):: inci

      open(10, file="param")

      call creadparam(10)
      call cbeginrun
      call cmkincident(inci, fin)
      call cinitracking(inci)

      det = obssites(noofsites)%pos%xyz

      write(0, *) ' detector base pos in Exyz'
      write(0,'(3g17.9)')  det%x, det%y, det%z
      write(0,*) 'incident position in Exyz'
      write(0,*) inci%pos%xyz%x, inci%pos%xyz%y,
     *     inci%pos%xyz%z
      call cxyz2det(det, inci%pos%xyz, cc)
      write(0,*) 'incident position in Det'
      write(0,*) cc%x, cc%y, cc%z

      call cvecprod(detzaxis, detxaxis, detyaxis)
      write(0,*) '   '
      write(0,*) ' detector  X,Y,Z axis vecgor in Exyz'
      write(0,'(3g17.9)') detxaxis%x, detxaxis%y, detxaxis%z
      write(0,'(3g17.9)') detyaxis%x, detyaxis%y, detyaxis%z
      write(0,'(3g17.9)') detzaxis%x, detzaxis%y, detzaxis%z
!       we try to convert {X,Y,Z}primary into the one in
!           detector system-->a,b,c.
      call citransvectzx(1, detzaxis, detxaxis, xprimary, aa)
      call citransvectzx(2, detzaxis, detxaxis, yprimary, bb)
      call citransvectzx(2, detzaxis, detxaxis, zprimary, cc)
!      primary etc is the vector  representing the primary system axes
!     in the Detector system. Then, a given vector R in Detector system,
!     its component to the a,b,c is Rx, Ry, Rz in the Primary is primary
!     system.
!
      write(0,*) ' Primary systrem X,Y,Z axis vector in Exyz'
      write(0,'(3g17.9)')  xprimary%x, xprimary%y, xprimary%z
      write(0,'(3g17.9)')  yprimary%x, yprimary%y, yprimary%z
      write(0,'(3g17.9)')  zprimary%x, zprimary%y, zprimary%z

      write(0,*) ' For a given vector in the Detector system'
      write(0,*) ' take scaler product of  the following vector to get'
      write(0,*) ' its x value in  the primary system. For y, z, use '
      write(0,*) ' 2nd, 3rd row'
      write(*,'(3g17.9)')  aa%x, aa%y, aa%z
      write(*,'(3g17.9)')  bb%x, bb%y, bb%z
      write(*,'(3g17.9)')  cc%x, cc%y, cc%z

      end
      subroutine chooktrace
      end
      subroutine chookceren
      end
      subroutine chookcerens
      end
      subroutine chookcerene
      end
      subroutine chookbgrun
      end