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 82 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, 
#if defined (KEKB) || (KEKA)
#include "ZavoidUnionMap.h"
!////////////
!      write(0,*) ' rank=',mpirank,' ir=',SeedSave,' Zfirst=',
!     *   Zfirst.pos.depth*0.1, 'g/cm2'
!//////////////
#else
     * zfirst
#endif
      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 297 of file chookFlesh.f.

References o.

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

      integer i

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

chookbgevent()

subroutine chookbgevent

(

)

Definition at line 65 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
      call xihist                ! instanciate histogram
      call xclearnrfai           ! clear Nrfai region

Here is the call graph for this function:

chookbgrun()

subroutine chookbgrun

(

)

Definition at line 19 of file chookFlesh.f.

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

      implicit none
#include "Zmaxdef.h"
#include "Zmanagerp.h"
#include "../../SkelFlesh/Zprivate.h"
#include "Ztrack.h"
#if defined (DOMPI)
#include "mpif.h"
#include "Zmpi.h"
#else
      integer mpirank
      parameter(mpirank=0)
#endif
      real*8  temp
      character*100 msg
      integer icon
      integer i
      eventno = 0

      if(mpirank .eq. 0) then
!            namelist output
         call cwriteparam(errorout, 0)
!            primary information
         call cprintprim(errorout)
!            observation level information
         call cprintobs(errorout)
      endif

      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 272 of file chookFlesh.f.

      implicit none
      integer never             ! input & output
      never = 1

chookenevent()

subroutine chookenevent

(

)

Definition at line 200 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 225 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 281 of file chookFlesh.f.

      implicit none
      integer never   ! input & output
      never = 1

chooknepint()

subroutine chooknepint

(

integer

never

)

Definition at line 290 of file chookFlesh.f.

      implicit none
      integer never   ! input & output
      never = 1

chookobs()

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

Definition at line 153 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"
#include "Zprivate4.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( keepweight ) then
         n = 1
      else
         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
      endif

      do i = 1, n 
         call xobs(atrack, id)
      enddo

Here is the call graph for this function:

chooktrace()

subroutine chooktrace

(

)

Definition at line 237 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 312 of file chookFlesh.f.

References charge, chookobs(), code, coszenith, erg, false, mass, o, p, r, subcode, t, true, wgt, 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 
         atrack.wgt = 1.
         call chookobs(atrack, 2)
      enddo
      heobs = .false.

Here is the call graph for this function:

cpushinci()

subroutine cpushinci

(

)

Definition at line 349 of file chookFlesh.f.

References cinitstack(), cpush(), and csortstack().

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

      type(track)::atrack
      type(track)::zsave
      call cinitstack  ! empty the stack

      read(mdev)  nooflowe
!/////////
!      write(0,*) ' rank=',mpirank,' NoOfLowE=',NoOfLowE
!//////////
      do  i = 1, nooflowe
#if defined (KEKA) || (KEKB)
         zsave = zfirst  ! save; next destroys Zfirst
         read(mdev) 
#include "ZavoidUnionMap.h"
         atrack = zfirst
         zfirst = zsave   ! restore
#else
         read(mdev) atrack
#endif
!           aTrack is already complete track so push it directly.
         call cpush(atrack)
      enddo
!           sort stack dscendent order
      call csortstack
      

Here is the call graph for this function: