cinteNEP.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cintenep
subroutine	csetknock

Function/Subroutine Documentation

cintenep()

subroutine cintenep

(

)

Definition at line 4 of file cinteNEP.f.

References cerrormsg(), cintebomega(), cinteddb(), cintedelta(), cintedmes(), cinteds(), cinteeta(), cinteetap(), cintegzai(), cinteheavy(), cintekaon(), cintelambda(), cintelambdac(), cintemuon(), cintennb(), cintenuc(), cinteomega(), cintephi(), cintepion(), cinterho(), cintesigma(), cintetau(), csetknock(), kalfa, kbomega, kddb, kdmes, kds, keta, ketap, kgnuc, kgzai, khvymax, kkaon, klambda, klambdac, kmuon, knnb, knuc, komega, kphi, kpion, krho, ksigma, ktau, and ktriton.

Referenced by cinteraction().

      implicit none

#include  "Zcode.h"
#include  "Ztrack.h"
#include  "Ztrackv.h"
!
!
!          used to judge if user hook should be called 
!          after MovedTrack interacted.
!
!
      character*80 msg
!
!   **  ptcl stacking is done in each subroutine; should be changed
!           (except for hadronic interactions)
!////////////
!      logical deb
!      common /cccdeb/deb
!/////////////

      
!/////////
!      if(deb) then
!         write(*,*) 'NEP:proc=',  IntInfArray(ProcessNo).process
!      endif
!////////////

      if(intinfarray(processno)%process .eq. 'knock') then
         call csetknock
      else
         if(movedtrack%p%code .eq. kpion) then
            call cintepion
         elseif(movedtrack%p%code .eq. kkaon) then
            call cintekaon
         elseif(movedtrack%p%code .eq. knuc) then
            call cintenuc
         elseif(movedtrack%p%code .eq. kmuon) then
            call cintemuon
         elseif( movedtrack%p%code .eq. kgnuc ) then
            call cinteheavy
         elseif(movedtrack%p%code .ge. kalfa .and.
     *           movedtrack%p%code .le. khvymax  ) then
!cc            elseif(MovedTrack.p.code .ge. kdeut .and.
            call cinteheavy
         elseif(movedtrack%p%code .eq. ktriton) then
            call cinteheavy
         elseif(movedtrack%p%code .eq. kdmes ) then
            call cintedmes
         elseif(movedtrack%p%code .eq. knnb) then
            call cintennb
         elseif(movedtrack%p%code .eq. kddb) then
            call cinteddb
         elseif(movedtrack%p%code .eq. keta) then
            call cinteeta
         elseif(movedtrack%p%code .eq. kgzai) then
            call cintegzai
         elseif(movedtrack%p%code .eq. kbomega) then
            call cintebomega
         elseif(movedtrack%p%code .eq. klambda) then
            call cintelambda
         elseif(movedtrack%p%code .eq. klambdac) then
            call cintelambdac
         elseif(movedtrack%p%code .eq. ksigma) then
            call cintesigma
         elseif(movedtrack%p%code .eq. krho) then
            call cinterho
         elseif(movedtrack%p%code .eq. komega) then
            call cinteomega
         elseif(movedtrack%p%code .eq. kphi) then
            call cintephi
         elseif(movedtrack%p%code .eq. kds ) then
            call cinteds
         elseif(movedtrack%p%code .eq. ketap ) then
            call cinteetap
         elseif(movedtrack%p%code .eq. kdelta ) then
            call cintedelta
         elseif(movedtrack%p%code .eq. ktau ) then
            call cintetau
         else
            write(0,*) ' E=', movedtrack%p%fm%p(4)
            write(msg, *) ' ptcl =', movedtrack%p%code,
     *        ' interaction=', 
     *         intinfarray(processno)%process,
     *        ' should not occure'
            call cerrormsg(msg, 0)
         endif
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

csetknock()

subroutine csetknock

(

)

Definition at line 93 of file cinteNEP.f.

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

Referenced by cintenep().

      implicit none
#include  "Zcode.h"
#include  "Ztrack.h"
#include  "Ztrackv.h"
!

       real*8 e1, er, tmp, cos1, cosr,  cs, sn, sinr, sine
       type(coord)::dc
       type(coord)::dc1
       type(coord)::dcr
!
       type(track)::atrack

       integer idx   ! 
 
       atrack = movedtrack
!                                surv  elec  surv   elec
       call  cknockea(atrack%p,  e1,   er,  cos1,   cosr)
!        we must put recoil electron into Pwork first and
!        then survival high energy ptcl.  This order will
!        be reversed when particles in Pwork are stacked.
!      (this inversion is to manage nuclear interactions where
!       the leading particle is put last in Pwork).
!        Therefore the low energy electron is first extracted.
!  Otherwise, stack overflow may take place for inclinde muons
!     (<= v5.33)
!         idx is to resolve the problem with a minumum effort.
!
!          survival paticle angle negligible always so
!          you may  put dc.r = (0, 0, 1)
       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%p%fm%p(4) = e1
       call csetdircos(dc1, atrack)
       call ce2p(atrack)
!           Nproduced = Nproduced + 1 
       idx = nproduced + 2
!           Pwork(Nproduced) = aTrack.p
       pwork(idx) = 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)
!
       call cmkptc(kelec,  0,  -1,  atrack%p)

       call ce2p(atrack)
!            Nproduced = Nproduced + 1
       idx =  nproduced + 1
       pwork(idx) = atrack%p

       nproduced = nproduced + 2

Here is the call graph for this function:

Here is the caller graph for this function: