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 124 of file chook.f.

References cqincident(), cqinirn(), cqnoofprim(), and cqprime().

      implicit none
#include "Ztrack.h"

    type(track):: atrack
    type(coord):: angleatobs
!
      integer no,seed(2)
      real*8 pore
!      write(*, *) ' bigin event generation'
      call cqinirn(seed)
      call cqincident(atrack,angleatobs)
      call cqnoofprim(no)
      call cqprime(pore)

Here is the call graph for this function:

chookbgrun()

subroutine chookbgrun

(

)

Definition at line 79 of file chook.f.

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

      implicit none
#include "Zmanagerp.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)
!     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
!      open component disk file, and get some over all
!     informaiton such as number of component showers
     
      call chookbgrunas 
      

Here is the call graph for this function:

chookeint()

subroutine chookeint

(

integer

never

)

Definition at line 297 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 202 of file chook.f.

References cqfirstid(), and depth.

      implicit none
#include "Ztrack.h"
#include "Ztrackv.h"
#include "Zobs.h"
#include "Zobsp.h"
#include "Zobsv.h"
      include "ZASbyH.h"



!      type(track):: inci
!      type(coord):: angle
      integer i
      real*8 fdepth

!

!      open(59,file='cosmos.dat',access='append')
!      write(59,*) (-1,i=1,8)
      if(observeas) then
         call cqfirstid(fdepth)
!         call cqIncident(inci, angle)
!         bsin = cgetBsin(inci.p, Mag)*1.e4
!                   electron size in B approx.
!         write(59, *) (ASObsSites(i).esize, i=1, NoOfASSites)
!                   size weighted age
!         write(*, *) (ASObsSites(i).age,   i=1, NoOfASSites)
         do i = 1, noofassites
             write(*, *) 's ',  sngl(asobssites(i).pos.depth/10.),
     *           sngl(asobssites(i).esize),
     *           sngl(asobssites(i).age)
!     *          ,sngl(fdepth)
!     *          ,sngl(bsin)
         enddo
         write(*, *)
!   close(59)
      endif

 

Here is the call graph for this function:

chookenrun()

subroutine chookenrun

(

)

Definition at line 249 of file chook.f.

      implicit none

chookgint()

subroutine chookgint

(

integer

never

)

Definition at line 326 of file chook.f.

            implicit none

#include  "Ztrack.h"
#include  "Ztrackv.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 354 of file chook.f.

            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 non zero, this routine will never be called.
!
!   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
!
!        IntInfArray(ProcessNo).process  will have
!             'col' or 'decay'
!   %%%%%%%%%%%%%%%%%%%%%%%%%%
!         replace interacting nucleons by c.s. give never = 2
!     so that no further development by the interacting nucleons.
!
      call chooknepinta
      never = 2

chookobs()

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

Definition at line 147 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 "Zcode.h"
#include "Ztrack.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

!
!     For id =2, you need not output the z value, because it is always
!     0 (within the computational accuracy).
!
      if(id .eq. 2) then
!            output typical quantities.
!        write(*, *) 
!     *  aTrack.where,   !  observation level. integer*2.  1 is highest.
!     *  aTrack.p.code,    !  ptcl code.  integer*2.
!     *  aTrack.p.charge,  !  charge,  integer*2 
!     *  sngl(aTrack.t), !  relateive arrival time in nsec (NOT sec).
!                        !  if TimeStructure is F, nonsense.
!     *  sngl(aTrack.p.fm.e)  ! total energy in GeV.
!     *  sngl(aTrack.pos.xyz.x), sngl(aTrack.pos.xyz.y),  !  x, y in m
!     *  sngl(aTrack.vec.w.x),  ! direc. cos.x in the current detector system.
!     *  sngl(aTrack.vec.w.y),  ! direc. cos.y
!     *  sngl(aTrack.vec.w.z),  ! direc. cos.z
!     *  sngl(aTrack.vec.coszenith) ! cos of zenith angle
!         if(aTrack.p.code .eq. kelec) then
!            write(*, *) aTrack.where
!         endif
!      endif
!         you may need in some case other information such as
!       aTrack.p.subcode   ! sub code of the particle integer*2
!       aTrack.p.mass      ! mass 
!       aTrack.wgt         ! weight of the particle (may not be 1. if
!                           ! ThinSampling =T)
!       aTrack.p.fm.x      ! momentum x component.  Note. Momentum is
!                            given in the  Earth xyz system.

!       aTrack.p.fm.y      !          y
!       aTrack.p.fm.z      !          z

      endif

chooktrace()

subroutine chooktrace

(

)

Definition at line 261 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 > 60. 
!         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 110 of file chook.f.

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