chook.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	time (xxx)
subroutine	chookbgrun
subroutine	chookbgevent
subroutine	chookobs (aTrack, id)
subroutine	chookenevent
subroutine	chookenrun
subroutine	chookeint (never)
subroutine	chookgint (never)
subroutine	chooknepint (never)
subroutine	chooktrace

Function/Subroutine Documentation

chookbgevent()

subroutine chookbgevent

(

)

Definition at line 39 of file chook.f.

References cquhookr(), cut, nc, sume, and sumepi0.

Referenced by ceventloop().

      implicit none
#include "Ztrack.h"
#include "Zobs.h"
#include "Zobsv.h"
      include "Zprivate.h"

      integer i, j
      sume = 0
      sumepi0 = 0
      nc = 0
      ncpi0= 0 
      call cquhookr(1, cut)

Here is the call graph for this function:

Here is the caller graph for this function:

chookbgrun()

subroutine chookbgrun

(

)

Definition at line 15 of file chook.f.

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

Referenced by cbeginrun().

      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.

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

Here is the call graph for this function:

Here is the caller graph for this function:

chookeint()

subroutine chookeint

(

integer

never

)

Definition at line 162 of file chook.f.

Referenced by cinteraction(), and ctracking().

            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',or  'anihi'
!

Here is the caller graph for this function:

chookenevent()

subroutine chookenevent

(

)

Definition at line 116 of file chook.f.

Referenced by ceventloop().

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

      integer i, j
!      do i = 1, NoOfSites
!         do j = 1, 3
!            write(*, *) i, j, nc(j, i)
!         enddo
!      enddo
      if(observeas) then
!                   electron size in B approx.
!         write(*,*) 'e ', sngl(sume), nc, sngl(sumepi0), ncpi0
         do i = 1, noofassites
            write(*, *) 's ',  sngl( asdepthlist(i)*0.1 ),
     *       sngl(asobssites(i).esize), sngl(asobssites(i).age)
         enddo
      endif
!      write(*,*)

Here is the caller graph for this function:

chookenrun()

subroutine chookenrun

(

)

Definition at line 147 of file chook.f.

Referenced by cendrun().

      implicit none
#include "Ztrack.h"
#include "Ztrackp.h"

Here is the caller graph for this function:

chookgint()

subroutine chookgint

(

integer

never

)

Definition at line 191 of file chook.f.

Referenced by cinteraction().

            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' or 'photop'
!       

Here is the caller graph for this function:

chooknepint()

subroutine chooknepint

(

integer

never

)

Definition at line 219 of file chook.f.

References charge, code, cut, knuc, kpion, mass, nc, p, process(), sume, and sumepi0.

Referenced by cinteraction().

            implicit none

#include  "Ztrack.h"
#include  "Ztrackv.h"
#include  "Zcode.h"
      include "Zprivate.h"
      integer never   ! input & output
      integer i

!         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
!
!        default setting is 1
      never = 1

!
!        IntInfArray(ProcessNo).process  will have
!             'coll' or 'decay'
      if(cut .lt. 1.) then
         if( intinfarray(processno).process   .eq. 'coll' ) then
            do i = 1, nproduced
!              for pions, kaons put E=mass if X > cut (for cut>0)
               if( pwork(i).code .lt. knuc ) then
                  if( pwork(i).fm.p(4) / movedtrack.p.fm.p(4)
     *                 .gt. cut) then
                     sume = sume +  pwork(i).fm.p(4)
                     nc = nc +1
                     if( pwork(i).code .eq. kpion .and.
     *                   pwork(i).charge .eq. 0 ) then
                        sumepi0 = sumepi0 + pwork(i).fm.p(4)
                        ncpi0= ncpi0 + 1
                     endif
                     pwork(i).fm.p(4) = pwork(i).mass*1.1
                  endif
               endif
            enddo
         endif
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

chookobs()

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

Definition at line 59 of file chook.f.

Referenced by cobservation(), and cobshes().

!
!     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 "Ztrack.h"
#include  "Zcode.h"
      include "Zprivate.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 .and.  aTrack.vec.coszenith .gt. 0 ) then
!         if(aTrack.p.code .le. 3) then
!            nc(aTrack.p.code, aTrack.where) = 
!     *        nc(aTrack.p.code, aTrack.where) +1
!         endif
!      endif
!            output typical quantities.
!        write(*, '(3i5, ) 
!     *  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.p(4))  ! total energy in GeV.
!     *  sngl(aTrack.pos.xyz.r(1)), sngl(aTrack.pos.xyz.r(2)), !  x, y in m
!     *  sngl(aTrack.vec.w.r(1)),  ! direc. cos.x in the current detector system.
!     *  sngl(aTrack.vec.w.r(2)),  ! direc. cos.y
!     *  sngl(aTrack.vec.w.r(3)),  ! direc. cos.z
!     *  sngl(aTrack.vec.coszenith) ! cos of zenith angle
!      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.p(1)      ! momentum x component.  Note. Momentum is
!                            given in the  Earth xyz system.

!       aTrack.p.fm.p(2)      !          y
!       aTrack.p.fm.p(3)      !          z

Here is the caller graph for this function:

chooktrace()

subroutine chooktrace

(

)

Definition at line 275 of file chook.f.

References ccoordfortr(), charge, code, p, and r.

Referenced by cputtrinfo().

            implicit none

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


!      h1 = TrackBefMove.pos.height- ObsSites(NoOfSites).pos.height
!      h2 = MovedTrack.pos.height - ObsSites(NoOfSites).pos.height
      
      type(coord)::  f, t



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

      if( movedtrack.p.charge .eq. 0 ) return      
!               convert coord.
      call ccoordfortr( 21,  f,  t )
      
 
      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, *) 
            write(*, *) 
            write(*, *) 
         endif
!         write(TraceDev, '(4g16.8, i3,i3)')
         write(*, '(4g12.4, i3,i3)')
     *        f.r(1), f.r(2), f.r(3), trackbefmove.t, 
     *        trackbefmove.p.code, trackbefmove.p.charge
         
      endif
!      write(TraceDev, '(4g16.8,i3,i3)') 
      write(*, '(4g12.4,i3,i3)') 
     *        t.r(1), t.r(2), t.r(3), movedtrack.t,
     *        movedtrack.p.code,   movedtrack.p.charge
         
      xxx = t.r(1)
      yyy = t.r(2)
      zzz = t.r(3)
      kkk = movedtrack.p.code
      chg = movedtrack.p.charge

Here is the call graph for this function:

Here is the caller graph for this function:

time()

subroutine time

(

integer

xxx

)

Definition at line 5 of file chook.f.

Referenced by __mknrfai.f__(), __mknrfaifromdat.f__(), __reduceeachsize.f__(), __reducesize.f__(), kdcmjd(), kmjd(), kmjdym(), ksided(), and ksidet().

      integer xxx
      xxx = 1

Here is the caller graph for this function: