chook.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	chookbgrun
integer function	csighandler (sig, code, context)
subroutine	chookbgevent
subroutine	chookobs (aTrack, id)
subroutine	chookenevent
subroutine	chookenrun
subroutine	chooktrace
subroutine	chookeint (never)
subroutine	chookgint (never)
subroutine	chooknepint (never)

Function/Subroutine Documentation

chookbgevent()

subroutine chookbgevent

(

)

Definition at line 54 of file chook.f.

      call xbgevent

chookbgrun()

subroutine chookbgrun

(

)

Definition at line 10 of file chook.f.

References cprintobs(), cprintprim(), csighandler(), cwriteparam(), and xbgrun().

      implicit none
#include "Zmanagerp.h"
#include "Ztrack.h"
#include "Ztrackv.h"

!
!

!         If you feel writing the parameters on stderr is
!         a bother, comment out the next or
!         use other device than ErrorOut.
!         Also you may comment out all output routines below.
#ifdef sun4
      external csighandler
      integer  ieeer, ieee_handler
      ieeer = ieee_handler('set', 'invalid', csighandler)
#endif

!
!            namelist output
      call cwriteparam(errorout, 0)
!            primary information
      call cprintprim(errorout)
!            observation level information
      call cprintobs(errorout)

      call xbgrun

Here is the call graph for this function:

chookeint()

subroutine chookeint

(

integer

never

)

Definition at line 160 of file chook.f.

            implicit none

#include  "Ztrack.h"
#include  "Ztrackv.h"
!  #include  "Ztrackp.h"
      
      integer never   ! input & output
      
!         don't make never = 1, if you want to get
!         information after an electron made interaction
!         if this is made non zero, this routine will never be called.
!
!   MovedTrack is the electron that made interaction
!   Pwork contains produced particles.
!   Nproduced has the number of particles in Pwork
!   IntInfArray(ProcessNo) contains the type of interaction
!
!        default setting
      never = 1
!
!        IntInfArray(ProcessNo).process will have one of
!       'brems', 'mscat', 'bscat', 'anihi' or 'mbrem'
!

chookenevent()

subroutine chookenevent

(

)

Definition at line 92 of file chook.f.

References cqeventno().

      implicit none
#include "Zmanagerp.h"
#include "Ztrack.h"
#include "Ztrackv.h"
#include "Zobs.h"
#include "Zobsp.h"
#include "Zobsv.h"
      integer  nevent, ntevent
 
      call xenevent

      call cqeventno(nevent, ntevent)
      write(0,*) 'end_of_event #=',ntevent

Here is the call graph for this function:

chookenrun()

subroutine chookenrun

(

)

Definition at line 111 of file chook.f.

      implicit none
      write(0,*) " end of run"

chookgint()

subroutine chookgint

(

integer

never

)

Definition at line 189 of file chook.f.

            implicit none

#include  "Ztrack.h"
#include  "Ztrackv.h"
!  #include  "Ztrackp.h"
      
      integer never   ! input & output
      
!         don't make never = 1, if you want to get
!         information after a gamma ray made interaction
!         if this is made non zero, this routine will never be called.
!
!   MovedTrack is the gamma that made interaction
!   Pwork contains produced particles.
!   Nproduced has the number of particles in Pwork
!   IntInfArray(ProcessNo) contains the type of interaction
!
!        default setting
      never = 1
!         IntInfArray(ProcessNo).process will have one of
!        'pair', 'comp', 'photoe' 'photop' 'mpair'
!       

chooknepint()

subroutine chooknepint

(

integer

never

)

Definition at line 217 of file chook.f.

References cirot3vec(), crot3vec(), cutptcl(), d0, p, and process().

            implicit none


#include  "Zcode.h"
#include  "Ztrack.h"
#include  "Ztrackv.h"

!  #include  "Ztrackp.h"
      
      integer never   ! input & output
      
!         don't make never = 1, if you want to get
!         information after a non-e-g particle  made interaction
!         if this is made 1, this routine will never be called.
!         make never 5 if  the event is not good for you
!         and retry the event generation.
!   MovedTrack is the particle that made interaction
!   Pwork contains produced particles.
!   Nproduced has the number of particles in Pwork
!   IntInfArray(ProcessNo) contains the type of interaction
!
!        default setting
!      never = 1
!
      integer::i
      type(ptcl):: proj
      real(8):: pz 
 
      if(  intinfarray(processno).process /= 'coll' ) then  ! even if this is absent, OK
         return                                             ! ;inside cutptcl same is done
      endif

!       first make Z axis as the incident particle
      do i = 1, nproduced
         call cirot3vec(i, movedtrack.p.fm, pwork(i).fm, pwork(i).fm)
      enddo

!           make projectile move +z direction
!        with the same energy as MovedTrack.p
      proj =  movedtrack.p
      pz =sqrt( dot_product(proj.fm.p(1:3), proj.fm.p(1:3)) )
      proj.fm.p(1:3) =(/0.d0, 0.d0, pz/)
      call cutptcl(proj)

!          again convert Z axis to be the one where
!       MovedTrack is defined
      do i = 1, nproduced
         call crot3vec( movedtrack.p.fm, pwork(i).fm, pwork(i).fm )
      enddo
      never = 0

Here is the call graph for this function:

chookobs()

subroutine chookobs	(	type(track)	aTrack,
		integer	id
	)

Definition at line 64 of file chook.f.

!
!     Note that every real variable is in double  precision so
!     that you may output it in sigle precision to save the memory.
!     In some cases it is essential to put it in sigle (say,
!     for gnuplot).
! 
      implicit none
#include "Zmanagerp.h"
#include "Zcode.h"
#include "Ztrack.h"
#include "Ztrackv.h"
#include  "Zheavyp.h"
#include "Zobs.h"
#include "Zobsp.h"

      integer id  ! input.  1 ==> aTrack is going out from
!                                 outer boundery.
!                           2 ==> reached at an observation level
!                           3 ==> reached at inner boundery.
      type(track):: atrack

      call xobs(atrack, id)

chooktrace()

subroutine chooktrace

(

)

Definition at line 124 of file chook.f.

References height.

            implicit none

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

       real*4 h1,  h2
!
!    Every time a particle is moved in the atmosphere, this routine is called,
!    if trace > 100
!         For a one track segment,
!     TrackBefMove  has  track information at the beginning of the segment.
!     MoveTrack    has   track information at the end of the segment.
!   
!     You can know the  information a track contains in the 
!     chookObs routine. (Note however, no conversion of coordinate
!     has been done.  The values are in the Earth xyz system.)
!     Besides quantities explained there, you can use, for a  given 'track'
!
!     atrack.pos.xyz.x, atrack.pos.xyz.y, atrack.pos.xyz.z    (x,y.z)
!     atrack.pos.radiallen   (distance from the center of the earth)
!     atrack.pos.depth       (vertical depth)
!     atrack.pos.height      (vertical heigth from sea level)  
!

      h1 = trackbefmove.pos.height- obssites(noofsites).pos.height
      h2 = movedtrack.pos.height - obssites(noofsites).pos.height

csighandler()

integer function csighandler	(	integer	sig,
		integer	code,
		integer, dimension(5)	context
	)

Definition at line 41 of file chook.f.

      implicit none
#include "Zmanagerp.h"
      integer sig, code, context(5)
      write(errorout, *)  ' f.p exception content=' , context(4)
!      call abort()