chookFlesh.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	chookbgrun
subroutine	chookbgevent
subroutine	cbegin1ev (nomore)
subroutine	chookobs (aTrack, id)
subroutine	chookenevent
subroutine	chookenrun
subroutine	chooktrace
subroutine	chookeint (never)
subroutine	chookgint (never)
subroutine	chooknepint (never)
subroutine	cgethes (from)
subroutine	cobshes
subroutine	cpushinci

Function/Subroutine Documentation

cbegin1ev()

subroutine cbegin1ev

(

integer

nomore

)

Definition at line 67 of file chookFlesh.f.

References cgethes(), cinitracking(), cmkincident(), and cobshes().

      implicit none
#include "../../SkelFlesh/Zprivate.h"
#include "Ztrack.h"
#include "Ztrackv.h"
#include "Ztrackp.h"
#include "Zobs.h"
#include "Zobsp.h"
#include "Zobsv.h"
#include "Zcode.h"
#include "Zmanager.h"
#include "Zmanagerp.h"
      
      integer nomore       !  output. 0 still there  are showers
                           !          1 no more skeleton showers to be fleshed
!          event number, primary      

      type(track):: incident, zsave
      type(coord):: angle
      
      integer i
      integer seed(2)
      integer cumnum, num, jeof, fin
      read( mdev, end=1000, err=999 ) cumnum, num, seedsave, zfirst

      eventsintherun = eventsintherun + 1
      eventno = eventno + 1
!                 reset the seed.
      call rnd1r(seedsave)
!         next incident; confirmed to be the same one as preserved one
      call cmkincident(incident, fin)
      if(fin .ne. 0 ) goto 1000
      zsave = zfirst    ! save;  this is reset in next 
      call cinitracking( incident )   
!          set first interaction pos
      zfirst = zsave
!      call cresetTimer(Zfirst)



!          do your own init for a one event here
!      ==========================================================


!      ==========================================================
!

      call cgethes(mdev)  ! get high energy ptlcs
      call cobshes        ! imitate their observation
      nomore = 0
      return

 1000 continue
      nomore = 1
      return
 999  continue
      write(0,*) ' Mdev read err'
      stop 1111

Here is the call graph for this function:

cgethes()

subroutine cgethes

(

integer

from

)

Definition at line 270 of file chookFlesh.f.

References o.

      implicit none
#include "../../SkelFlesh/Zprivate.h"
      integer from

      integer i

      read(from)  np
      do i = 1, np
         read(from) o(i)
      enddo

chookbgevent()

subroutine chookbgevent

(

)

Definition at line 51 of file chookFlesh.f.

References cbegin1ev(), cerrormsg(), and cpushinci().

      implicit none
#include "../../SkelFlesh/Zprivate.h"


      integer nomore
      call cbegin1ev( nomore )
      if( nomore .eq. 1) then
         call cerrormsg('all events have been fleshed', 1)
         stop                   !!!!!!!!!!!!  
      endif
      call cpushinci
 
      call xbgevent
      return

Here is the call graph for this function:

chookbgrun()

subroutine chookbgrun

(

)

Definition at line 15 of file chookFlesh.f.

References cerrormsg(), cgetfname(), copenfw2(), cprintobs(), cprintprim(), cquhookc(), cquhooki(), cwriteparam(), and xbgrun().

      implicit none
#include "Zmanagerp.h"
#include "../../SkelFlesh/Zprivate.h"

      real*8  temp
      character*100 msg
      integer icon
      integer i
      eventno = 0

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

      call cquhooki(1, mdev)      ! get skeleton memo dev #
      call cquhookc(1, msg)       ! get file name for sekelton data
      call cgetfname(msg, mskel)  ! add host name etc if needed
      call copenfw2(mdev, mskel, 2, icon)
      if(icon .ne. 1) then
         call cerrormsg(mskel,1)
         call cerrormsg(' could not be opened',0)
      endif

      call xbgrun
      call ihist                ! instanciate histogram

Here is the call graph for this function:

chookeint()

subroutine chookeint

(

integer

never

)

Definition at line 245 of file chookFlesh.f.

      implicit none
      integer never             ! input & output
      never = 1

chookenevent()

subroutine chookenevent

(

)

Definition at line 173 of file chookFlesh.f.

References cfintracking().

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

      integer i

!       for Job ='newflesh', we must call cfinTracking ourselves.
      call cfintracking
!           end of 1 event; if you need to do some here is
!           the place

      call xenevent

Here is the call graph for this function:

chookenrun()

subroutine chookenrun

(

)

Definition at line 198 of file chookFlesh.f.

References cprintstatus().

      implicit none
      call  cprintstatus   ! if don't like,  comment out

Here is the call graph for this function:

chookgint()

subroutine chookgint

(

integer

never

)

Definition at line 254 of file chookFlesh.f.

      implicit none
      integer never   ! input & output
      never = 1

chooknepint()

subroutine chooknepint

(

integer

never

)

Definition at line 263 of file chookFlesh.f.

      implicit none
      integer never   ! input & output
      never = 1

chookobs()

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

Definition at line 131 of file chookFlesh.f.

References rndc(), and wgt.

!
!     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.  2 ==> reached at an observation level
!                           1 ==> aTrack is going out from
!                                 outer boundery.
!                           2 ==> reached at an observation level
!                           3 ==> reached at inner boundery.
      type(track):: atrack
!
!
      integer n, i
      real*8 eps, u

      if(atrack.wgt .gt. 1.) then
         n=atrack.wgt
         eps = atrack.wgt - n
         call rndc(u)
         if(u .lt. eps) then
            n = n + 1
         endif
      else
         n = 1
      endif
      do i = 1, n 
         call xobs(atrack, id)
      enddo

Here is the call graph for this function:

chooktrace()

subroutine chooktrace

(

)

Definition at line 210 of file chookFlesh.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

cobshes()

subroutine cobshes

(

)

Definition at line 283 of file chookFlesh.f.

References charge, chookobs(), code, coszenith, erg, false, mass, o, p, r, subcode, t, true, x, xyz, and y.

      implicit none
#include "../../SkelFlesh/Zprivate.h"
#include "Ztrack.h"
!
!           memorized high energy showers at the skeleton making
!     time is put into the chookObs as if they are really observed
      type(track):: atrack

      integer i
      logical  heobs            ! if T, currently observing 
      common /zheobs/ heobs     !  particles those obsrved at skeelton making time

      heobs = .true.
      do i = 1, np
         atrack.where =  o(i).where 
         atrack.p.code =  o(i).code 
         atrack.p.subcode = o(i).subcode 
         atrack.p.charge = o(i).charge 
         atrack.t = o(i).atime 
         atrack.p.fm.p(4) = o(i).erg
         atrack.p.mass = o(i).mass 
         atrack.pos.xyz.r(1) = o(i).x 
         atrack.pos.xyz.r(2) = o(i).y 
         atrack.vec.w.r(1) = o(i).wx
         atrack.vec.w.r(2) = o(i).wy
         atrack.vec.w.r(3) = o(i).wz
         atrack.vec.coszenith = o(i).zenith 
         call chookobs(atrack, 2)
      enddo
      heobs = .false.

Here is the call graph for this function:

cpushinci()

subroutine cpushinci

(

)

Definition at line 319 of file chookFlesh.f.

References charge, cinitstack(), code, cpush(), csortstack(), depth, mass, and p.

      implicit none
#include "../../SkelFlesh/Zprivate.h"
#include "Ztrack.h"
#include "Ztrackv.h"
      integer i

      type(track)::atrack

      call cinitstack  ! empty the stack

      read(mdev)  nooflowe
      do  i = 1, nooflowe
         read(mdev) atrack
!           aTrack is already complete track so push it directly.
         call cpush(atrack)
         write(*,'(2i4,g13.4, g13.4)')
     *     atrack.p.code, atrack.p.charge, 
     *     atrack.p.fm.p(4)-atrack.p.mass, atrack.pos.depth*0.1
      enddo
      stop
!           sort stack dscendent order
      call csortstack
      

Here is the call graph for this function: