cinteElec.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cinteelec
subroutine	cbrems
subroutine	cknockon (cg)
subroutine	canihi
subroutine	cmbrem

Function/Subroutine Documentation

canihi()

subroutine canihi

(

)

Definition at line 167 of file cinteElec.f.

References canihiea(), ce2p(), cmkptc(), csetdircos(), ctransvectz(), d0, kcasg, kcossn(), and kphoton.

Referenced by cinteelec().

!     *****************
       implicit none
!----       include '../Particle/Zcode.h'
#include  "Zcode.h"
!----       include 'Ztrack.h'
#include  "Ztrack.h"
!----       include 'Ztrackv.h'
#include  "Ztrackv.h"

       real*8 e1, er, cos1, cosr, tmp, sine
       real*8 cs, sn, sinr
       type(track)::atrack
       type(coord)::dc
       type(coord)::dc1
       type(coord)::dcr
!
       call canihiea(movedtrack%p%fm%p(4), e1, er, cos1, cosr)
!               
       tmp=1.d0-cos1*cos1
       if(tmp .lt. 0.d0) then
          tmp=0.d0
          cos1=-1.d0
       endif
       sine=sqrt(tmp)
       call kcossn(cs,sn)
       dc%r(1) = cs*sine
       dc%r(2) = sn*sine
       dc%r(3) = cos1
!             
       call ctransvectz(movedtrack%vec%w, dc, dc1)
!            higher gamma
       atrack = movedtrack
       call csetdircos(dc1, atrack)
       call cmkptc(kphoton, kcasg, 0,  atrack%p)
       atrack%p%fm%p(4) = e1
       call ce2p(atrack)


       nproduced = nproduced + 1
       pwork(nproduced) = atrack%p
!       call cpush(aTrack)
!                low gamma
       tmp=1.d0-cosr*cosr
       if(tmp .lt. 0.d0) then
          tmp=0.d0
          cosr=-1.d0
       endif
       sinr=sqrt(tmp)
       dc%r(1) = -cs*sinr
       dc%r(2) = -sn*sinr
       dc%r(3) = cosr
       call ctransvectz(movedtrack%vec%w, dc, dcr)
       call csetdircos(dcr, atrack)
       atrack%p%fm%p(4) = er
       call ce2p(atrack)


       nproduced = nproduced + 1
       pwork(nproduced) = atrack%p
!       call cpush(aTrack)

Here is the call graph for this function:

Here is the caller graph for this function:

cbrems()

subroutine cbrems

(

)

Definition at line 34 of file cinteElec.f.

References cbremang(), cbremerglpm(), cbremsenergy(), ce2p(), cmkptc(), cos, csetdircos(), ctransvectz(), d0, kcasg, kcossn(), and kphoton.

Referenced by cinteelec().

!        *****************
         implicit none

#include  "Zcode.h"
#include  "Ztrackp.h"
#include  "Ztrack.h"
#include  "Ztrackv.h"


         real*8 beg, den, cvh2den
         real*8 theta, sint, cost, sn, cs
         type(coord)::dc1
         type(coord)::w
         type(track)::atrack

!           get brem gamma energy
         if(lpmeffect .and. 
     *        trackbefmove%p%fm%p(4) .gt. lpmbrememin) then ! use BefMove%p
!            den = cthick2den(TrackBefMove.pos.depth)
            den = cvh2den(trackbefmove%pos%height)  ! this should be better
            call cbremerglpm(trackbefmove%p%fm%p(4), den, beg)
         else
            call cbremsenergy(movedtrack%p%fm%p(4), beg)
         endif

         atrack = movedtrack

!         
         atrack%p%fm%p(4) = movedtrack%p%fm%p(4) -beg
         call ce2p(atrack)
!
!         see           !! note above
         pwork(nproduced+2) = atrack%p


          atrack%p%fm%p(4) = beg
          call cmkptc(kphoton,  kcasg, 0,  atrack%p)
!           brems angle at low energy
         if(  trackbefmove%p%fm%p(4)  .lt. 10. ) then
             call cbremang( trackbefmove%p%fm%p(4), trackbefmove%p%mass, 
     *       beg, targetatomicn, theta)
             ! v7.646; v8.03 revived
!            call cBremAng(1, TrackBefMove%p%fm%p(4), beg, theta)  removed
            if(theta .lt. 0.03d0) then
               sint = theta
               cost = 1.- theta**2 / 2
            else
                sint = sin(theta)
                cost = cos(theta)
             endif
             call kcossn(cs,sn)
             w%r(1) = cs*sint
             w%r(2) = sn*sint
             w%r(3) = cost
             call ctransvectz(movedtrack%vec%w, w, dc1)
             call csetdircos(dc1, atrack)
          endif
          call ce2p(atrack)
          pwork(nproduced+1) = atrack%p
          nproduced = nproduced + 2

Here is the call graph for this function:

Here is the caller graph for this function:

cinteelec()

subroutine cinteelec

(

)

Definition at line 2 of file cinteElec.f.

References canihi(), cbrems(), cerrormsg(), cknockon(), and cmbrem().

Referenced by cinteraction().

       implicit none
!----       include 'Ztrack.h'
#include  "Ztrack.h"
!----       include 'Ztrackv.h'
#include  "Ztrackv.h"
         character*70 msg
!
         if(intinfarray(processno)%process .eq. 'brems') then
            call  cbrems
         elseif(intinfarray(processno)%process .eq. 'mscat') then
            call cknockon(-1)
         elseif(intinfarray(processno)%process .eq. 'bscat') then
            call cknockon(1)
         elseif(intinfarray(processno)%process .eq. 'anihi') then
            call canihi
         elseif(intinfarray(processno)%process .eq. 'mbrem') then
            call cmbrem
         else

              write(msg, *) ' process#=', processno,
     *       'and process for electron=',
     *        intinfarray(processno)%process, ' undef'
              call cerrormsg(msg, 0)
         endif

Here is the call graph for this function:

Here is the caller graph for this function:

cknockon()

subroutine cknockon

(

integer

cg

)

Definition at line 97 of file cinteElec.f.

References cbhabhaea(), ce2p(), cmkptc(), cmollerea(), csetdircos(), ctransvectz(), d0, kcossn(), and kelec.

Referenced by cinteelec().

!      *****************  moller and bhabha scattering
       implicit none
!----       include '../Particle/Zcode.h'
#include  "Zcode.h"
!----       include 'Ztrack.h'
#include  "Ztrack.h"
!----       include 'Ztrackv.h'
#include  "Ztrackv.h"
#include  "Zelemagp.h"
!     
       integer cg
       real*8 e1, er, tmp, cos1, cosr, sine, cs, sn, sinr
       type(coord)::dc
       type(coord)::dc1
       type(coord)::dcr
!
       type(track)::atrack
!
       if(cg .eq. -1) then
          call cmollerea(movedtrack%p%fm%p(4), recoilkinemine,
     *                   e1, er, cos1, cosr)
       else
          call cbhabhaea(movedtrack%p%fm%p(4), recoilkinemine,
     *    e1, er, cos1, cosr)
       endif

       tmp=1.d0-cos1*cos1
       if(tmp .lt. 0.d0) then
          tmp=0.d0
          cos1=1.d0
       endif
       sine=sqrt(tmp)
       call kcossn(cs,sn)
       dc%r(1) = cs*sine
       dc%r(2) = sn*sine
       dc%r(3) = cos1
!           convert angle
       call ctransvectz(movedtrack%vec%w, dc, dc1)
       atrack = movedtrack
       atrack%p%fm%p(4) = e1
       call csetdircos(dc1, atrack)
       call ce2p(atrack)

       pwork(nproduced + 2) = atrack%p
!            knock on electron
       tmp=1.d0-cosr*cosr
       if(tmp .lt. 0.d0) then
          tmp=0.d0
          cosr=1.d0
       endif
       sinr=sqrt(tmp)
       dc%r(1) = -cs*sinr
       dc%r(2) = -sn*sinr
       dc%r(3) = cosr
       call ctransvectz(movedtrack%vec%w, dc, dcr)
       atrack%p%fm%p(4) = er
       call csetdircos(dcr, atrack)
!
       if(cg .eq. 1) then
          call cmkptc(kelec,  0,  -1,  atrack%p)
       endif

       call ce2p(atrack)

       pwork(nproduced + 1) = atrack%p

       nproduced = nproduced + 2

Here is the call graph for this function:

Here is the caller graph for this function:

cmbrem()

subroutine cmbrem

(

)

Definition at line 231 of file cinteElec.f.

References ce2p(), cmbreme(), cmkptc(), kcasg, kphoton, and upsilon().

Referenced by cinteelec().

!          magnetic brems
!        *****************
         implicit none
#include  "Zcode.h"
#include  "Ztrack.h"
#include  "Ztrackv.h"

         real*8 x, beg

         type(track)::atrack

!           get brem gamma energy
         call cmbreme(upsilon, x) 
         beg = x * movedtrack%p%fm%p(4)  ! gamma energy

         atrack = movedtrack
         atrack%p%fm%p(4) = movedtrack%p%fm%p(4) -beg
         call ce2p(atrack)
         pwork( nproduced + 2 ) = atrack%p   
!           stack gamma (probably lower energy than e)
         atrack%p%fm%p(4) = beg
         call cmkptc(kphoton,  kcasg, 0,  atrack%p)
         call ce2p(atrack)
         pwork(nproduced + 1) = atrack%p

         nproduced = nproduced + 2         

Here is the call graph for this function:

Here is the caller graph for this function: