doxygen/v8.037/SkelFleshMu_2chookSkel_8f_source.html

 ! This one is for getting
 ! Xcm,  eta, y and ptcl code,  at the first collision,
 ! which generated observed muons.c
 !
 #include "cmain.f"
 #include "chookHybAS.f"
 #include "ctemplCeren.f"
 !
 !  *************************************** hook for Beginning of a Run
 !  * At this moment, all (system-level) initialization for this run
 !  * has been ended.  After this routine is executed, the system goes into the
 !  * event creation loop.
 !  *
       subroutine chookbgrun
       implicit none
 #include "Zmanagerp.h"
 #include "Ztrack.h"
 #include "Zprivate.h"

       character*10  append
       character*100 msg
       logical  ex, opn
       integer klena

 !            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 cquhooki(2, wdev)      ! get working disk dev #


       call cquhooki(3, ngmin)     ! get Nh min
       call cquhooki(4, nhmin)     ! get Ng min
       call cquhooki(5, where)     ! where to check

       call cquhookr(1, sumegmin)      ! sum E min
       call cquhookr(2, sumehmin)


       call cquhookc(1, msg)       ! get file name for sekelton memo
       call cgetfname(msg, mskel)  ! add host name etc if needed
       call cquhookc(2, msg)       ! get file name for working
       call cgetfname(msg, wskel)  ! add host name etc if neededn
       call cquhookc(3, append)    ! append data, if Mskel already exists

       write(msg, *) 'Skeleton is  judged at obs.pos=', Where
       call cerrormsg(msg, 1)
       write(msg, *) ' Ngmin=',ngmin,
      *     ' SumEGmin=',sumegmin/1000.,'TeV'
       call cerrormsg(msg, 1)
       write(msg, *) ' Nhmin=',nhmin,
      *    ' SumEHmin=',sumehmin/1000.,'TeV'
       call cerrormsg(msg, 1)

 !
       inquire(file=mskel(1:klena(mskel)), opened=opn, exist=ex)
       if(opn) then
          call cerrormsg(mskel, 1)
          call cerrormsg(' already opened: starange', 0)
       elseif(ex .and.
      *      append(1:klena(append)) .eq. 'append' ) then
          open(mdev, file=mskel, form='unformatted',status='old')
          call cerrormsg('skeleton node info. will be appended', 1)
 !           skip to the end of file
          do while( .true. )
             read(mdev, end=100)
          enddo
       else
          if(ex .and.
      *      append(1:klena(append)) .ne. 'append' ) then
             call cerrormsg(
      *           'Old skeleton node info. file exists', 1)
             call cerrormsg(
      *           'but node info. will NOT be appended', 1)
          endif
          open(mdev, file=mskel(1:klena(mskel)), form='unformatted',
      *       status='unknown')
       endif

  100  continue

       open(wdev, file=wskel(1:klena(wskel)), form='unformatted',
      *      status='unknown' )

       accepted = 0   ! counter;  accepted as skeleton
 !
       end


 !     *********************************** hook for Beginning of  1 event
 !     *  All system-level initialization for 1 event generation has been
 !     *  eneded at this moment.
 !     *  After this is executed, event generation starts.
 !     *
       subroutine chookbgevent
       implicit none

 #include "Ztrack.h"
 #include "Ztrackv.h"
 #include "Ztrackp.h"
 #include "Zobs.h"
 #include "Zobsp.h"
 #include "Zobsv.h"
 #include "Zcode.h"
 #include "Zprivate.h"


       type(coord )::angle


       integer  EventNo
       integer i, j
       integer seed(2)
       real*8 svEasWait, svEthin(2), kepn
       save

       bgevent = .true.
       np = 0
       call cqincident( incident, angle)
       kepn = incident.p.fm.p(4)
       if( incident.p.code .eq. kgnuc ) then
          kepn = kepn/ incident.p.subcode
       endif
       ethresh = kepn * waitratio

       sveaswait = easwait       ! for safety save
       svethin(1) = ethin(1)           ! //
       svethin(2) = ethin(2)           ! //
       call csetemin(generate2, keminobs2, cutneg, cuteg)
       easwait = sveaswait       ! restore
       ethin(1) = svethin(1)
       ethin(2) = svethin(2)

       rewind  wdev

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

       eventno = eventno + 1
       do i = 1, noofsites
          write(*, 1000)
      *    sngl(obssites(i).pos.depth),
      *    eventno,
      *    incident.p.code,
      *    incident.p.subcode,
      *    incident.p.charge,
      *    incident.p.fm.e,
      *    -angle.r(1),
      *    -angle.r(2),
      *    -angle.r(3)
  1000       format(f10.3,i9,3i4,e15.5,3(1x,f12.8))
       enddo


 !     ===================================================
       end


 !     ************************************ hook for observation
 !     *  One particel information is brought here by the system.
 !     *  All information of the particle is in aTrack
 !     *
       subroutine chookobs(aTrack, id)
 !
 !     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 "Zprivate.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
 !
 !     For id =2, you need not output the z value, because it is always
 !     0 (within the computational accuracy).
 !
 !     ===================================================
       common /testcos/sumg, ng(20), nth, eventno
       real*8 sumg
       integer ng, nth, EventNo

 !     ===================================================
       real*8 user
       real*4 userE(2)
       integer*2 userI(4)
       equivalence(user, usere(1)), (user, useri(1))
       save


       if( id .eq. 2 .and. atrack.p.code .ne. kneumu .and.
      *   atrack.p.code .ne. kneue) then
          np = np + 1
          if( np .gt. npmax) then
             call cerrormsg(
      *      '# of particles >NpMax in observation', 0)
          endif
          o(np).where = atrack.where
          o(np).code = atrack.p.code
          o(np).subcode = atrack.p.subcode
          o(np).charge = atrack.p.charge
          o(np).atime = atrack.t
          o(np).erg = atrack.p.fm.p(4)
          o(np).mass = atrack.p.mass
          o(np).x = atrack.pos.xyz.r(1)
          o(np).y = atrack.pos.xyz.r(2)
          o(np).wx =atrack.vec.w.r(1)
          o(np).wy =atrack.vec.w.r(2)
          o(np).wz =atrack.vec.w.r(3)
          o(np).zenith = atrack.vec.coszenith
          o(np).wgt = atrack.wgt
          o(np).user =  atrack.user
 !     ===================================================
          if(o(np).code .eq. 3) then
             user=atrack.user
             write(*, 959)
 !     *  o(Np).where,
 !     *  o(Np).code,
 !     *  o(Np).charge,
      *  ( o(np).erg - o(np).mass ),
 !     *  ( o(Np).x ), ( o(Np).y ) ,
 !     *  ( o(Np).wx ),
 !     *  ( o(Np).wy ),
 !     *  ( o(Np).wz ),
 !     *  ( o(Np).zenith ),
      *    useri(1)/1000., useri(2)/1000., usere(2)

 ! 959    format(3i3, g14.3,2f16.6,4(1x,f12.8), i3, 1pE11.3)
  959    format(3g14.3,e11.3)
       endif

 !     ===================================================
       endif
 !     ===================================================

 !     ===================================================
       end

 !    *********************************** hook for end of 1 vent
 !    * At this moment, 1 event generation has been ended.
 !    *
       subroutine chookenevent

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

       integer i
       integer seed(2)
       real  sumeg, sumeh
       logical memorize

       integer ng, nh

       ng = 0
       nh = 0

       sumeg = 0
       sumeh =0
 !         count sum Eg etc
       do i = 1, np
          if(o(i).where .eq. where) then
             if(o(i).code .le. kelec) then
                ng = ng + 1
                sumeg = sumeg + o(i).erg
             elseif( ( o(i).code .ge. kpion .and.
      *                o(i).code .le. knuc )
      *             .or.
      *               o(i).code .eq. kgnuc) then
                nh = nh + 1
                sumeh = sumeh + o(i).erg
             endif
          endif
       enddo

 !     ===================================================
       memorize =(ng .ge. ngmin .and. sumeg .ge. sumegmin) .or.
      *  ( nh .ge. nhmin .and. sumeh .ge. sumehmin)

 !/////
       write(0,*) ' memo=', memorize,
      *   ' ng=',ng, ' seg=',sumeg, ' nh=',nh,' seh=',sumeh
 !//////


 !     ===================================================
       if(memorize) then
          accepted = accepted + 1
          call cwriteseed        !  SeedFile
 !          flag for end of 1 event on working disk
          write(wdev)  -1, p
          call cmemorize(wdev, mdev)     !  reocord this event
       endif
       end


 !     ********************************* hook for end of a run
 !     *  all events have been created or time lacks
 !     *
       subroutine chookenrun
       implicit none
 #include "Ztrack.h"
 #include "Zprivate.h"
       character*100 msg
 !     =========================================


 !     =========================================
       call cerrormsg('++++++++++++', 1)
       write(msg, '(i8, a)') accepted,
      *       ' events are memorized as skeleton'
       call cerrormsg(msg, 1)
       call cerrormsg('their  seeds are also memorized', 1)
       call cerrormsg('-----------',1)
       call  cprintstatus   ! if don't like,  comment out
       end
 !     ********************************* hook for trace
 !     *  This is called only when trace > 100
 !     *  User should manage the trace information here.
 !     *  If you use this, you may need some output for trace
 !     *  at the beginning of 1 event generatio and at the end of  1 event
 !     *  generation so that you can identfy each event.
 !     *
 !     *
       subroutine chooktrace
             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

       end
 !     ********************* this is the hook called when
 !       an electron made an interaction.
 !
       subroutine chookeint(never)
       implicit none

 #include  "Ztrack.h"
 #include  "Ztrackv.h"
 #include  "Zcode.h"
 #include  "Zprivate.h"
 !  #include  "Ztrackp.h"
       integer never             ! input & output
       integer nlow
 !
 !    If Job = 'newskel', input  "never" may be < 0,  and MovedTrack
 !    may not be an electron.
 !      never = -1 :  MovedTrack K.E becomes < KEminObs due to
 !                    energy loss during traveling.
 !      never = -2 :  The same as above, but the particle crosses an
 !                    observation depth, so the calling point to this
 !                    routine is different from the never = -1 case.
 !      never = -3 :  K.E >= KEminiObs.  The ptcl is observed at the
 !                    current deepest Obs. level. but at the flesh time
 !                    the deepest level will be more deep so that
 !                    this must be memorized.
 !
 !         For above cases, the product is only MovedTrack and
 !         no particle is in PWork.
 !  Otherwise,
 !   MovedTrack is the electron that made interaction
 !   Pwork contains produced particles (normally gamma, but mayb  e).
 !   Nproduced has the number of particles in Pwork.
 !   IntInfArray(ProcessNo) contains the type of interaction
 !        IntInfArray(ProcessNo).process will have one of
 !       'brems', 'mscat', 'bscat'  'anihi' or 'mbrem'
 !
       if(never .lt. 0) then
          nproduced = 1
 !           IBM user must modify next
          pwork(1) = movedtrack.p
       endif

 !         high energy parent at node might be used
 !        for hybrid AS generation if it is in some
 !        energy region.
       if( movedtrack.p.code .eq. kelec ) then
          if( movedtrack.asflag .eq. 0 ) then
             if( movedtrack.p.fm.p(4) .lt. ethresh ) then
                movedtrack.asflag = -1
             endif
          endif
       endif

       call cmemonode(wdev, never, 2)
       if(movedtrack.asflag .eq. -1) then
 !          decendent should not be used to  generae A.S
          movedtrack.asflag  = -2
       endif

  10   continue
       never = 0
       end

 !     ********************* this is the hook called when
 !       a gamma ray made an interaction.
 !
       subroutine chookgint(never)
       implicit none

 #include  "Ztrack.h"
 #include  "Ztrackv.h"
 #include  "Zprivate.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
 !         IntInfArray(ProcessNo).process will have one of
 !        'pair', 'comp', 'photoe' 'photop' 'mpair'
       call cmemonode(wdev, 1, 1)
       never = 0
       end

 !     ********************* this is the hook called when
 !       non e-g particle made an interaction.
 !
       subroutine chooknepint(never)
       implicit none
 #include  "Zmaxdef.h"
 #include  "Zcode.h"
 #include  "Ztrack.h"
 #include  "Ztrackv.h"
 #include  "Ztrackp.h"
 #include  "Zstackv.h"
 #include  "Zprivate.h"

       integer never   ! input & output

       type(coord):: angle
       type(ptcl)::  aptcl
       type(track):: aTrack
       type(ptcl):: tgt, pj
       real*8 user
       real*4 userE(2)
       integer*2 userI(4)
       equivalence(user, usere(1)), (user, useri(1))
       real*8 y, eta
       integer i, icon
       logical yes/.false./
       integer stackpos
       save

 !            never < 0 never happen
       if( bgevent ) then
 !               just  after collision or before 1st col.
 !               (knockon)
          yes = ( movedtrack.p.code .eq. knuc .or.
      *          movedtrack.p.code .eq. kgnuc ) .and.
      *          intinfarray(processno).process .eq. 'coll'
          if(yes)  then
 !               just after collision of nucleon or heavy primary

 !                copy particle property
             incip = movedtrack.p
 !               userE(2) is used to store  Xcm of a particle
 !               but the last bit is set if parent is pi/K
 !               else  it is reset to 0
 !               userI(1) = eta *1000 userI(2)= y*1000.
 !
             bgevent = .false.
             useri(1)=6
             usere(2) = 1.
             movedtrack.user =  user
             call cmkptc(6, 0, 1, tgt)
             tgt.fm.p(1) = 0.
             tgt.fm.p(2) = 0.
             tgt.fm.p(3) = 0.
             tgt.fm.p(4) = tgt.mass
             pj = incip
             pjcm = pj
             call cirot3vec(1, incip, incip, pj)
             call cgeqm(pj, tgt, cmsptcl, icon)
             call cbst1(1, cmsptcl, pj, pjcm)
 !////////////////
 !            write(*,*) ' pjcm=',
 !     *           pjcm.fm.p(1), pjcm.fm.p(2), pjcm.fm.p(3)
 !////////////
 !              first int. product is in stack
             call cgetcurrentstackpos(stackpos)

             do i = stackpos-nstacked+1, stackpos
                aptcl=stack(i).p
                call cirot3vec(1, incip, stack(i).p, aptcl)
                call cbst1(1, cmsptcl, aptcl, aptcl)
                call cyeta(aptcl, y, eta)
                useri(1)=eta*1000.
                useri(2)=y * 1000.
 !                    xcm
                usere(2)= aptcl.fm.p(3)/pjcm.fm.p(3)
 !/////////////////
 !            if(aptcl.code .le. 5 .and. userE(2) .gt. -1.0) then
 !               write(*,'(a,3i3, 4g11.3)')  'p: ',  aptcl.code,
 !     *         aptcl.subcode, aptcl.charge,
 !     *          eta, y, userE(2), Stack(i).p.fm.p(4)
 !            endif
 !/////////////////
 !                set ptcl code by using last bit
                if(aptcl.code .le. 5) then
                   call ksetbit(usere(2), 1)
                else
                   call krsetbit(usere(2), 1)
                endif
 !                    store in stack
                stack(i).user=user

             enddo
          endif
       endif
 !         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
 !
       if(never .eq. 0) then
          call cmemonode(wdev, 1, 3)
       endif
       end

       subroutine cmemonode(dev, flag, neporeg)
       implicit none
 #include  "Zcode.h"
 #include  "Ztrack.h"
 #include  "Ztrackv.h"
 #include  "Ztrackp.h"
 #include  "Zprivate.h"
 !
       integer dev  ! input. fortran logical dev. number for data write
       integer flag ! input. If this is -3, high energy particles must
                    !       also be memorized.
       integer neporeg ! input 1; from g 2 from e; 3 from NEP
 !
 !
 !       memorize nodal information at interaction.
 !
 !        \
 !         \  projectile    = MovedTrack
 !         /|\  produced particles.: Pwork(i), i=1~Nproduced
 !        / | \
 !
 !   From projectile, the track information is extracted
 !   From produced particles, only those of K.E< KEminObs is
 !   extracted and information needed for further tracking is
 !   obtaned and memorized. The position information is common
 !   to all the children.
 !
 !   Track   information;   pos =
 !                                xyz = r(1~3), sys
 !                                radiallen, depth, height, colheight
 !                            t
 !                          vec  =
 !                                 w = r(1~3), sys
 !                                 coszenith
 !                          wgt
 !                         where
 !                        asflag
 !
 !      Among these, we don't memorize
 !         sys which is always 'xyz'
 !       radiallen: can be computed from height
 !         vec; children knows  their direction
 !         wgt: normally not needed, it should be 1 for skeleton making
 !              So thinsamping must not be used when making skeleton.
 !      asflag: should be always F, (assume for skeleton making, hybrid
 !              air shower is not requested)
 !
 !   write  track info. of projectile
 !
       integer i, nlow
       real*8 ke
       real*8  y, eta
       real*8 user
       real*4 userE(2)
       integer*2 userI(4)
       equivalence(user, usere(1)), (user, useri(1))
       type(ptcl):: aptcl
       save

 !
 !
       nlow = 0
       do i = 1, nproduced
          ke = pwork(i).fm.p(4) - pwork(i).mass
          if( ( pwork(i).code .le. kelec .and.
      *       ke .ge.  cuteg ) .or.
      *       ( pwork(i).code .gt. kelec .and.
      *       ke .ge.  cutneg ) ) then
 !              count low energy ptcls
             if(flag .ne. -3) then
                if( ke .lt.  keminobs) then
                   nlow = nlow + 1
                endif
             else
 !               all ptcl must be memorized
                nlow = nlow + 1
             endif
          endif
       enddo
 !         next is harmful for B-approx air shower
 !         because, parent is not memorized and cannot
 !         generate AS at flesh time at proper positon
 !         asflag == -1 must apear in flesh
 !      if(nlow .eq. 0 )  return !  *************
 !       if there is no return here, the size of skeleton data
 !       becomes 3.2 times larger

       if(nlow .eq. 0  .and. movedtrack.asflag .ne. -1) return


       p.posx = movedtrack.pos.xyz.r(1)
       p.posy = movedtrack.pos.xyz.r(2)
       p.posz = movedtrack.pos.xyz.r(3)
       p.depth = movedtrack.pos.depth
       p.height = movedtrack.pos.height
       if( movedtrack.pos.colheight .gt. 1.e36) then
          p.colheight = 1.e36    ! no col. yet.
       else
          p.colheight =  movedtrack.pos.colheight
       endif
       p.atime = movedtrack.t
       p.where = movedtrack.where
       p.coszenith = movedtrack.vec.coszenith
       p.code = movedtrack.p.code
       p.erg  = movedtrack.p.fm.p(4)
       p.asflag = movedtrack.asflag
       p.user =  movedtrack.user

 !c
 !c     *           p.posx, p.posy, p.posz, p.depth, p.height,
 !c     *           p.colHeight, p.atime, p.where
 !c
 !       write particle info
 !           p(1~4)
 !           mass
 !           code
 !           subcode
 !           charge
 !


       write(dev) nlow, p

       do i = 1, nproduced
          ke = pwork(i).fm.p(4) - pwork(i).mass
          if( ( pwork(i).code .le. kelec .and.
      *       ke .ge.  cuteg ) .or.
      *       ( pwork(i).code .gt. kelec .and.
      *       ke .ge.  cutneg ) ) then
             if(flag .eq. -3 .or.  ke .lt. keminobs) then
                c.code = pwork(i).code
                c.subcode =  pwork(i).subcode
                c.charge =  pwork(i).charge
                c.fm(1) = pwork(i).fm.p(1)
                c.fm(2) = pwork(i).fm.p(2)
                c.fm(3) = pwork(i).fm.p(3)
                c.fm(4) = pwork(i).fm.p(4)
                c.mass = pwork(i).mass
                if(bgevent .and. neporeg .eq. 3 ) then
                   aptcl = pwork(i)
                   call cirot3vec(1, incip,  pwork(i),  aptcl)
                   call cbst1(1, cmsptcl, aptcl, aptcl)
                   call cyeta(aptcl, y, eta)
                   useri(1)=eta*1000.
                   useri(2)=y * 1000.
 !                    xcm
                   usere(2)= aptcl.fm.p(3)/pjcm.fm.p(3)
 !                     set ptcl code by using last bit
                   if(pwork(i).code .le. 5) then
                      call ksetbit(usere(2), 1)
                   else
                      call krsetbit(usere(2), 1)
                   endif
                   c.user = user
 !////////////////////
 !                  if(userE(2) .gt. -1. .and.
 !     *              Pwork(i).code .le. 5) then
 !                     write(*,'(a,3i3, 4g11.3)')
 !     *                 'l: ', Pwork(i).code, Pwork(i).subcode,
 !     *                    Pwork(i).charge,
 !     *                    eta, y, userE(2), Pwork(i).fm.p(4)
 !                  endif
 !/////////////////
                else
                   c.user = p.user
                endif
                write(dev) c
             endif
          endif
       enddo
       end

 !
       subroutine cmemorize(from, to)
       implicit none
 #include  "Ztrack.h"
 #include  "Ztrackv.h"

       integer from      !  working disk file
       integer to        !  permanent disk file where skeleton is sotered


       integer num, cumnum, irevent(2)
 !      type(track):: incident
 !      type(coord):: angle


       rewind from
 !          need not memorize, can be generated at flesh
 !      call cqIncident(incident, angle)

       call cqeventno(num, cumnum)

       call cqinirn(irevent)   ! seed of the event

       write(to) cumnum, num, irevent,  zfirst
       call cputhes(to)   ! put high energy showers.
       write(0,*) ' first Z=',zfirst.pos.depth*0.1, ' g/cm2',
      *  zfirst.pos.height, ' m'

       call cputnodinfo(from, to)  ! put nordal info.

       end

       subroutine cputhes(to)
       implicit none
 #include "Ztrack.h"
 #include "Zprivate.h"
       integer to
 !
 !        write high energy sekeleton data into 'to'
 !
       integer i

       write(to) np

       do i = 1, np
          write(to) o(i)
       enddo

       end

       subroutine cputnodinfo(from, to)
       implicit none
 #include "Ztrack.h"
 #include "Zprivate.h"

       integer from, to

       integer i, nlow

       nlow = 1
       do while ( nlow .ge. 0 )
          read(from) nlow,  p
          write(to)  nlow,  p
          do i = 1, nlow
             read(from) c
             write(to) c
          enddo
       enddo

       end


cgetfname
subroutine cgetfname(fnin, fn)
Definition: copenf.f:275

depth
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos depth
Definition: ZavoidUnionMap.h:1

cerrormsg
subroutine cerrormsg(msg, needrtn)
Definition: cerrorMsg.f:4

cprintobs
subroutine cprintobs(io)
Definition: cprintObs.f:2

cqincident
subroutine cqincident(incident, AngleAtObs)
Definition: cmkIncident.f:486

krsetbit
subroutine krsetbit(i, bit)
Definition: kmanbit.f:53

cbst1
subroutine cbst1(init, p1, p2, po)
Definition: cbst1.f:54

chookgint
subroutine chookgint(never)
Definition: chook.f:191

cputhes
subroutine cputhes(to)
Definition: chookSkel.f:630

cputnodinfo
subroutine cputnodinfo(from, to)
Definition: chookSkel.f:647

e
dE dx *! Nuc Int sampling table e
Definition: cblkMuInt.h:130

ng
others if is ng
Definition: cblkManager.h:9

colheight
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos *Zfirst pos colheight
Definition: ZavoidUnionMap.h:1

r
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz r
Definition: ZavoidUnionMap.h:1

process
! parameters for Elemag process(-> ---------------------------------------------- real *8 RecoilKineMinE !2 Recoil Kinetic Min Energy above which the recoil(=knock-on process) ! is treated. Below this energy, the effect is included as continuous ! energy loss. Used only if KnockOnRatio $>$ 1. ! If this is 0 or if KnockOnRatio=1, KEminObs(gamma)=KEminObs(elec) is used. ! See also KnockOnRatio. real *8 KnockOnRatio !2 KnockOnRatio *KEminoObs is used instead of RecoilKineMinE if KnockOnRatio $< $1. real *8 X0 !2 Radiation length in kg/m$^2$ for air. Normally the user should not touch this. real *8 Ecrit !2 Critical energy in GeV. \newline ! Employed only when calculating air shower size in the hybrid ! air shower generation. The value would be dependent on the ! experimental purpose. The default value, 81 MeV, is bit too ! small in many applications(The air shower size is overestimated). ! Comparisons of sizes by the hybrid method and by the full Monte ! Carlo tell that \newline ! $N_e$(full 3-D M.C) $< N_e$(hybrid AS with $E_c=81$ MeV) $< N_e$(full 1-D M.C) ! $ {\ \lower-1.2pt\vbox{\hbox{\rlap{$<$}\lower5pt\vbox{\hbox{$\sim$}}}}\ } ! N_e$(hybrid AS with $E_c={76}$ MeV) at around shower maximum. ! Hybrid AS is always essentially 1-D. logical Knockon !2 Obsolete. Don 't use this. See RecoilKineMinE ! and KnockonRatio. real *8 AnihiE !2 If E(positron) $<$ AnihiE, annihilation is considered. real *8 Es !2 Modified scattering constant. 19.3d-3 GeV real *8 MaxComptonE !2 Above this energy, Compton scattering is neglected. real *8 MaxPhotoE !2 Above this energy, photoelectric effect is neglected. real *8 MinPhotoProdE !1 Below this energy, no photo-prod of hadron. See also PhotoProd. logical PhotoProd !1 Switch. if .false., no photo prod. of hadron is considered at all. ! See also MinPhotoProdE, HowPhotoP real *8 Excom1 !2(GeV). If photon energy is<=Excom1, use XCOM data for ! compton/p.e/coherent scattering(must be< 100 GeV). real *8 Excom2 !2(GeV). If photon energy is<=Excom2, use XCOM data for ! pair creation cross-section.(must be< 100 GeV). integer Moliere !2 2$\rightarrow$ use Moliere scat.\newline ! 0$\rightarrow$ use Gaussian scattrign. \newline ! 1$\rightarrow$ use Moli\`ere scattering for non-electrons \newline ! 2$\rightarrow$ use Moli\`ere scattering for all charged ! particles. But treatment is not so rigorous as case of 3. ! \newline ! 3$\rightarrow$ use rigorus Moliere scattering. Diff. from 2 is verysmall. May be some effect in the ! core region. integer ALateCor !2 1$\rightarrow$ angular and lateral correlation is taken into account when Moliere=0 .\newline ! t$\rightarrow$ Use angular-lateral correlation by Gaussian ! approximation. No effect is seen if path length is short. !<-) ---------------------------------------------- common/Zelemagc/RecoilKineMinE

kgnuc
max ptcl codes in the kgnuc
Definition: Zcode.h:2

cmemorize
subroutine cmemorize(from, to)
Definition: chookSkel.f:599

xyz
latitude latitude this system is used *****************************************************************! type coord sequence union map real z z in m endmap xyz map real ! latitude in deg is to the north ! longitude in deg is to the east *h ! height in m endmap llh map real ! polar angle ! azimuthal angle *radius ! radial distance endmap sph endunion character *sys ! which system xyz
Definition: Zcoord.h:25

chookenrun
subroutine chookenrun
Definition: chook.f:147

track
Definition: Ztrack.h:44

kelec
max ptcl codes in the kelec
Definition: Zcode.h:2

chooknepint
subroutine chooknepint(never)
Definition: chook.f:219

true
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon true
Definition: cblkElemag.h:7

cprintprim
subroutine cprintprim(out)
Definition: cprintPrim.f:3

kneue
max ptcl codes in the kneue
Definition: Zcode.h:2

nth
!onst int nth
Definition: Zprivate.h:2

chookenevent
subroutine chookenevent
Definition: chook.f:116

cgeqm
subroutine cgeqm(p1, p2, q, icon)
Definition: cgeqm.f:2

cwriteparam
subroutine cwriteparam(io, force)
Definition: cwriteParam.f:4

cirot3vec
subroutine cirot3vec(init, p1, p2, po)
Definition: crot3vec.f:97

csetemin
subroutine csetemin(gen, eminob, emin, emCas)
Definition: ciniTracking.f:242

p
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
Definition: cblkHeavy.h:7

ksetbit
subroutine ksetbit(i, bit)
Definition: kmanbit.f:36

cwriteseed
subroutine cwriteseed
Definition: cwriteSeed.f:15

cquhookr
subroutine cquhookr(i, rv)
Definition: cqUHookr.f:2

asflag
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos *Zfirst pos *Zfirst Zfirst vec w *Zfirst vec w Zfirst vec *Zfirst Zfirst Zfirst asflag
Definition: ZavoidUnionMap.h:1

code
*Zfirst p fm *Zfirst p Zfirst p code
Definition: ZavoidUnionMap.h:1

chooktrace
subroutine chooktrace
Definition: chook.f:275

cmemonode
subroutine cmemonode(dev, flag)
Definition: chookSkel.f:467

knuc
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code knuc
Definition: cblkHeavy.h:7

o
struct ob o[NpMax]
Definition: Zprivate.h:34

kneumu
max ptcl codes in the kneumu
Definition: Zcode.h:2

height
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos height
Definition: ZavoidUnionMap.h:1

y
latitude latitude this system is used *****************************************************************! type coord sequence union map real y
Definition: Zcoord.h:25

cqeventno
subroutine cqeventno(num, cumnum)
Definition: cqEventNo.f:3

coszenith
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos *Zfirst pos *Zfirst Zfirst vec w *Zfirst vec w Zfirst vec coszenith
Definition: ZavoidUnionMap.h:1

erg
float erg[maxp]
Definition: Zprivate.h:7

charge
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p charge
Definition: ZavoidUnionMap.h:1

chookobs
subroutine chookobs(aTrack, id)
Definition: chook.f:59

chookbgevent
subroutine chookbgevent
Definition: chook.f:39

cquhooki
subroutine cquhooki(i, iv)
Definition: cqUHookr.f:15

cmkptc
subroutine cmkptc(code, subcode, charge, p)
Definition: cmkptc.f:15

cqinirn
subroutine cqinirn(ir)
Definition: cwriteSeed.f:4

t
nodes t
Definition: cstdatmos0.fNew2.h:84

where
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos *Zfirst pos *Zfirst Zfirst vec w *Zfirst vec w Zfirst vec *Zfirst Zfirst where
Definition: ZavoidUnionMap.h:1

cquhookc
subroutine cquhookc(i, cv)
Definition: cqUHookr.f:28

ptcl
Definition: Zptcl.h:75

coord
Definition: Zcoord.h:43

mass
*Zfirst p fm *Zfirst p mass
Definition: ZavoidUnionMap.h:1

cprintstatus
subroutine cprintstatus
Definition: cendRun.f:35

chookbgrun
subroutine chookbgrun
Definition: chook.f:15

kpion
max ptcl codes in the kpion
Definition: Zcode.h:2

cyeta
subroutine cyeta(p, y, eta)
Definition: cyeta.f:5

wgt
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos *Zfirst pos *Zfirst Zfirst vec w *Zfirst vec w Zfirst vec *Zfirst wgt
Definition: ZavoidUnionMap.h:1

cgetcurrentstackpos
subroutine cgetcurrentstackpos(stackpos)
Definition: cstack.f:84

chookeint
subroutine chookeint(never)
Definition: chook.f:162

x
! structure defining a particle at production ! Basic idea of what is to be contained in ! the particle structue is that dynamical ones should be included those derivable from the particle code ! is not included ******************************************************type fmom momentum sequence union map real e endmap map real * x
Definition: Zptcl.h:21

subcode
*Zfirst p fm *Zfirst p Zfirst p Zfirst p subcode
Definition: ZavoidUnionMap.h:1

c
dE dx *! Nuc Int sampling table c
Definition: cblkMuInt.h:130