cobservation.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cobservation (id)

Function/Subroutine Documentation

cobservation()

subroutine cobservation

(

integer

id

)

Definition at line 4 of file cobservation.f.

References c, chookeabsorbb(), chookobs(), cxyz2det(), cxyz2detd(), cxyz2prim(), cxyz2primd(), horizontal, and perpendicular.

Referenced by ctracking().

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

      integer id !  input. 1 ==> the ptcl goes out from BorderHeightH
                 !         2 ==> the ptcl is observed at an observation level
                 !         3 ==> the ptcl reached at BorderHeightL
                 !          calles below happen only if DETAILED_TRACKING == 1
                 !         4 ==> the ptcl is going to make an interaction
                 !         5 ==> the ptcl is going to die
                 !         6 ==> the ptcl is being discarded because path >
                 !               limit  
                 !         7 ==> the ptcl is being discarded by angle limit
                 !               However, ondimentional mode  will not 
                 !               call this.
                 !         8 ==> the ptcl moved a step. 

      type(track)::otrack

      integer  iwhere

#if LABELING > 0
      if( movedtrack%vec%coszenith .lt. 0. ) then
!             upgoing.       
         if(movedtrack%where .eq. 1 .or. movedtrack%where .eq. 2 .and.
     *        id .eq. 2) then
!c      if(MovedTrack.where .le. 31 .and. id .eq. 2) then
            if(movedtrack%where .eq. 1) then
               movedtrack%info = movedtrack%info + 1
            else
               movedtrack%info = movedtrack%info + 65336
            endif
!c            iwhere=MovedTrack.where
!               set bit for where-th bit.
!            call ksetbit( MovedTrack.info,  iwhere )
         endif
      endif
#endif
      if(id .eq. 2) then
         otrack = movedtrack
!         no coord. conv. is made unless the ptcl is at observation
!         level.
! 
!           convert time into ns
         if(timestructure .and. abs(obsplane) .ge. 1
     *                    .and. abs(obsplane) .le. 2)  then
            otrack%t = ( otrack%t  - obssites(otrack%where)%minitime )
     *              /c *tonsec
         endif

!          convert coordinate and direction cos. 
!          if ObsPlane < 0, no conversion. 
         if(obsplane .eq. horizontal) then

            call cxyz2det(obssites(otrack%where)%pos%xyz, 
     *                 otrack%pos%xyz, otrack%pos%xyz)
!       
            call cxyz2detd(otrack%vec%w, otrack%vec%w)

         elseif(obsplane .eq. perpendicular) then

            call cxyz2prim(obssites(otrack%where)%pos%xyz, 
     *                 otrack%pos%xyz, otrack%pos%xyz)
            call cxyz2primd(otrack%vec%w, otrack%vec%w)
!///////////
!            write(0,*) ' bef chookobs'
!            write(0,*) 'xyz=', oTrack.pos.xyz
!            write(0,*) 's=', oTrack.vec.w
!////////////
         endif
!          call user routine      
         call chookobs(otrack, id)
         if(eabsorb(1) .ne. 0 ) then
!                 2nd int arg. for normal observation 
!                 point.
            call chookeabsorbb(otrack,2)
         endif
      else
!          The user may change the energy to discard the
!          particle
         call chookobs(movedtrack,id)
      endif

Here is the call graph for this function:

Here is the caller graph for this function: