seeSkel.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 68 of file seeSkel.f.

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

      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

      write(*,*) ' Zfirst=',zfirst.pos.depth

      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 258 of file seeSkel.f.

References charge, code, erg, o, and subcode.

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

      integer i

      read(from)  np
      do i = 1, np
         read(from) o(i)
         write(*,'(a, 4i3, 1pE11.3)' ) ' HE ',
     *          o(i).where, o(i).code, o(i).subcode, o(i).charge,
     *          o(i).erg
      enddo

chookbgevent()

subroutine chookbgevent

(

)

Definition at line 53 of file seeSkel.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

Here is the call graph for this function:

chookbgrun()

subroutine chookbgrun

(

)

Definition at line 17 of file seeSkel.f.

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

      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

Here is the call graph for this function:

chookeint()

subroutine chookeint

(

integer

never

)

Definition at line 233 of file seeSkel.f.

      implicit none
      integer never             ! input & output
      never = 1

chookenevent()

subroutine chookenevent

(

)

Definition at line 161 of file seeSkel.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 186 of file seeSkel.f.

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

chookgint()

subroutine chookgint

(

integer

never

)

Definition at line 242 of file seeSkel.f.

      implicit none
      integer never   ! input & output
      never = 1

chooknepint()

subroutine chooknepint

(

integer

never

)

Definition at line 251 of file seeSkel.f.

      implicit none
      integer never   ! input & output
      never = 1

chookobs()

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

Definition at line 134 of file seeSkel.f.

References code, and p.

!
!     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
!
!
       write(*,*) 'o ', atrack.p.code, atrack.p.fm.p(4)
!      call xObs(aTrack, id)

chooktrace()

subroutine chooktrace

(

)

Definition at line 198 of file seeSkel.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 274 of file seeSkel.f.

References charge, chookobs(), code, coszenith, erg, mass, o, p, r, subcode, t, 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

      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

Here is the call graph for this function:

cpushinci()

subroutine cpushinci

(

)

Definition at line 306 of file seeSkel.f.

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

      implicit none
#include "../../SkelFlesh/Zprivate.h"
#include "Zmaxdef.h"
#include "Ztrack.h"
#include "Ztrackv.h"
#include "Zstackv.h"
      integer i, remain

      type(track):: atrack
      type(track):: tt(max_stack_size)

      call cinitstack  ! empty the stack

      read(mdev)  nooflowe
      do  i = 1, nooflowe
         read(mdev) atrack
         call cpush(atrack)
      enddo
!           sort stack dscendent order
      call csortstack
      do i = 1, nooflowe
         call cpop(tt(i), remain)
      enddo
      do i = nooflowe, 1, -1
         write(*,* ) 'LE ', tt(i).p.code, tt(i).p.charge, 
     *   tt(i).p.fm.p(4)-tt(i).p.mass, tt(i).p.mass
      enddo
      call cinitstack           ! empty the stack

Here is the call graph for this function: