modefield Module Reference

Data Types
type	efield

Functions/Subroutines
subroutine	cefieldread
subroutine	cefieldwrite
subroutine	cefieldcheck
subroutine	cefielddet2xyz
subroutine	cefield0 (aTrack, Efout)
subroutine	cefieldcheckr (aTrack, i, out)

Variables
integer, parameter	nmaxefield =5
integer, save	howefield =0
character(len=1), save	defofr ='h'
type(efield), dimension(nmaxefield), save	myef
integer, save	nefield = 0
logical, dimension(nmaxefield), save	checkh = .false.
logical, dimension(nmaxefield), save	checkr = .false.
logical, dimension(nmaxefield), save	checkt = .false.
logical, save	useh =.false.
logical, save	user =.false.
logical, save	uset =.false.
real(8), dimension(3, nmaxefield), save	efxyz

Function/Subroutine Documentation

cefield0()

subroutine modefield::cefield0	(	type (track)	aTrack,
		real(8), dimension(3), intent(out)	Efout
	)

Definition at line 123 of file cEfield0.f.

References c, cefieldcheck(), cefieldcheckr(), cefielddet2xyz(), checkh, checkr, checkt, efxyz, false, myef, nefield, useh, and uset.

Referenced by cgetefield().

!      Electric field is given depending on particle position
!      and/or time.  
!      H          T                      R
!    given   not given --> H is used   R at that H is used. if R is out 
!                                      of range, E=0. else E at H
!                                      if R is not given  E at H
!  not given    given  --> T is used   same as above, use T instead of H.
!    given      given  --> T is used    //
!  not gvien   not given               use first R,.if it exists
!                                      if in range, use E(1) else E=0
!                                      if R absent, use E(1)
      implicit none
#include  "Zglobalc.h"
#include "Ztrack.h"
#include "Zincidentv.h"
      type(track):: atrack  ! input. Track info in E-xyz.
         ! position and "where" info are used.
      real(8),intent(out):: efout(3) ! Elecric field
                       ! in E-xyz system.  V/m.
      integer:: i, hindex, tindex, out
      logical,save:: first=.true.
      real(8)::time
      
      real(8)::r
      real(8)::zero(3)=(/0.,0.,0./)

      if( first ) then   ! convert Ef in det into E-xyz
         call cefieldcheck
         call cefielddet2xyz
         first = .false.
!!!!!!!
!         call cEfieldWrite
!!!!!!!!!
      endif

      hindex = 0
      if(useh) then
         do i = 1, nefield
            if(checkh(i)) then
               if( atrack%pos%height >= myef(i)%H1 .and.
     *              atrack%pos%height <= myef(i)%H2 ) then
                  hindex = i
                  exit
               endif
            endif
         enddo
         if( hindex == 0) then  ! H is out of range
            efout(:) = zero(:)
         else
              !  see if R is in the range
            call cefieldcheckr(atrack, hindex, out)
            if(out == 0 ) then
             !  Efout(:) = myEf(hindex)%Ef(:)  ! bug
               efout(:) = efxyz(:,hindex)
            else
               efout(:) = zero(:)
            endif
         endif
      elseif(uset) then
         tindex = 0
         do i = 1, nefield
            if(checkt(i)) then
               time = atrack%t/c*tonsec
               if( time >= myef(i)%T1 .and.
     *             time <= myef(i)%T2 ) then
                  tindex = i
                  exit
               endif
            endif
         enddo
         if( tindex == 0 ) then  ! T is out of range
            efout(:) = zero(:)
         else
              !  see if R is in the range
            call cefieldcheckr(atrack, tindex, out)
            if(out == 0 ) then
!               Efout(:) = myEf(tindex)%Ef(:)  bug
               efout(:) = efxyz(:,tindex)
            else
               efout(:) = zero(:)
            endif
         endif
      else  ! H,T is not used. 
         i = 1   ! use first data if possible
         if( checkr(i) ) then   ! check first R if any
            call cefieldcheckr(atrack, i, out)
            if(out == 0 ) then
!               Efout(:) = myEf(i)%Ef(:)
               efout(:) = efxyz(:,i)
            else
               efout(:) = zero(:)
            endif
         else  ! no check of R.  assume first Ef is applied
!            Efout(:) = myEf(i)%Ef(:)
            efout(:) = efxyz(:,i)
         endif
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cefieldcheck()

subroutine modefield::cefieldcheck

(

)

Definition at line 53 of file cEfield0.f.

References checkh, checkr, checkt, howefield, myef, nefield, nmaxefield, useh, user, and uset.

Referenced by cefield0().

      implicit none
      integer:: i

      if(howefield /= 1 ) return  !!!!!
      do i =  nmaxefield, 1, -1   ! last non zero Ef 
         if(dot_product(myef(i)%Ef(:),myef(i)%Ef(:))  > 0 ) then
            nefield = i
            exit
         endif
      enddo
      checkh(1:nefield) =
     *     myef(1:nefield)%H1**2+myef(1:nefield)%H2**2 > 0.
      checkr(1:nefield) =
     *     myef(1:nefield)%R1**2+myef(1:nefield)%R2**2 > 0.
      checkt(1:nefield) =
     *     myef(1:nefield)%T1**2+myef(1:nefield)%T2**2 > 0.
      do i = 1, nefield
         useh = useh .or. checkh(i)
         user = user .or. checkr(i)
         uset = uset .or. checkt(i)
      enddo
      useh = useh .and. .not. uset  ! useT has priority
      
      if( nefield >= 1 ) then
              ! check if height is from low to heigh
         do i = 1, nefield
            if( myef(i)%H1 > myef(i)%H2) then
               write(0,*)  i,'-th Efiled height error'
               write(0,*)
     *              'H1=', myef(i)%H1,' >', myef(i)%H2,'=H2'
               stop
            endif
            if( myef(i)%R1 > myef(i)%R2) then
               write(0,*)  i,'-th Efiled Radius error'
               write(0,*)
     *              'R1=', myef(i)%R1, ' >',  myef(i)%R2,'=R2'
               stop
            endif
            if( myef(i)%T1 > myef(i)%T2) then
               write(0,*)  i,'-th Efiled Time error'
               write(0,*)
     *              'T1=', myef(i)%T1, ' >',  myef(i)%T2,'=T2'
               stop
            endif
            if( nefield > i ) then
               if( myef(i+1)%H1 < myef(i)%H2 ) then
                  write(0,*)
     *                 i,'-th Ef height region<=',myef(i)%H2
                  write(0,*)
     *                 ' but is not < next region=',myef(i+1)%H1
                  stop
               endif
            endif                  
         enddo
      endif
      write(0,*) ' # of Efield(s)=',nefield

Here is the caller graph for this function:

cefieldcheckr()

subroutine modefield::cefieldcheckr	(	type (track)	aTrack,
		integer, intent(in)	i,
		integer, intent(out)	out
	)

Definition at line 223 of file cEfield0.f.

References cgethcoredist(), cgetpcoredist(), checkr, defofr, and myef.

Referenced by cefield0().

      implicit none
#include "Ztrack.h"
      type(track):: atrack  ! input
      integer,intent(in):: i ! R or H index
      integer,intent(out)::out ! 0 aTrack is in the range or no check is needed
                   !  else out of range
      real(8):: r

      if( checkr(i) ) then
                ! core distance
         if(defofr == 'p' ) then ! perpedicular distance
            call cgetpcoredist(atrack, r)
         elseif(defofr == 'h') then ! horizontal distance
            call cgethcoredist(atrack, r)
         else
            write(0,*) ' DefofR =',defofr, ' invalid'
            stop
         endif
         if( r >= myef(i)%R1 .and.
     *        r <= myef(i)%R2 ) then
            out = 0
         else                   ! R is out of range
            out = 1
         endif
      else
         out = 0
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cefielddet2xyz()

subroutine modefield::cefielddet2xyz

(

)

Definition at line 113 of file cEfield0.f.

References cdet2xyzd(), efxyz, myef, and nefield.

Referenced by cefield0().

      ! convert Ef in det system into E-xyz system
      ! should be called once before starting tracking
      implicit none
      integer::i
      do i = 1, nefield
         call cdet2xyzd(myef(i)%Ef(:), efxyz(:,i))
      enddo

Here is the call graph for this function:

Here is the caller graph for this function:

cefieldread()

subroutine modefield::cefieldread

(

)

Definition at line 43 of file cEfield0.f.

      implicit none
      read(*, efparam) 

cefieldwrite()

subroutine modefield::cefieldwrite

(

)

Definition at line 48 of file cEfield0.f.

      implicit none
        write(0, efparam) 

Variable Documentation

checkh

logical, dimension(nmaxefield), save checkh = .false.

Definition at line 28 of file cEfield0.f.

Referenced by cefield0(), and cefieldcheck().

        logical,save:: checkh(nmaxefield) = .false.

checkr

logical, dimension(nmaxefield), save checkr = .false.

Definition at line 29 of file cEfield0.f.

Referenced by cefield0(), cefieldcheck(), and cefieldcheckr().

        logical,save:: checkr(nmaxefield) = .false.

checkt

logical, dimension(nmaxefield), save checkt = .false.

Definition at line 30 of file cEfield0.f.

Referenced by cefield0(), and cefieldcheck().

        logical,save:: checkt(nmaxefield) = .false.

defofr

character(len=1), save defofr ='h'

Definition at line 22 of file cEfield0.f.

Referenced by cefieldcheckr().

        character(len=1),save:: defofr ='h'  ! define distance to the shower axis

efxyz

real(8), dimension(3,nmaxefield), save efxyz

Definition at line 37 of file cEfield0.f.

Referenced by cefield0(), and cefielddet2xyz().

        real(8),save::efxyz(3,nmaxefield) !copy of

howefield

integer, save howefield =0

Definition at line 17 of file cEfield0.f.

Referenced by cefieldcheck(), cgetefield(), cmagdef(), and cmaxmovlen().

        integer,save:: howefield=0  

myef

type(efield), dimension(nmaxefield), save myef

Definition at line 25 of file cEfield0.f.

Referenced by cefield0(), cefieldcheck(), cefieldcheckr(), and cefielddet2xyz().

        type(efield),save:: myef(nmaxefield)

nefield

integer, save nefield = 0

Definition at line 26 of file cEfield0.f.

Referenced by cefield0(), cefieldcheck(), and cefielddet2xyz().

        integer,save:: nefield = 0  ! acutal # of E flied 

nmaxefield

integer, parameter nmaxefield =5

Definition at line 14 of file cEfield0.f.

Referenced by cefieldcheck().

        integer,parameter:: nmaxefield=5

useh

logical, save useh =.false.

Definition at line 31 of file cEfield0.f.

Referenced by cefield0(), and cefieldcheck().

        logical,save:: useh=.false.

user

logical, save user =.false.

Definition at line 32 of file cEfield0.f.

Referenced by cefieldcheck().

        logical,save:: user=.false.

uset

logical, save uset =.false.

Definition at line 33 of file cEfield0.f.

Referenced by cefield0(), and cefieldcheck().

        logical,save:: uset=.false.