cinteNEP.f

Go to the documentation of this file.

!            treat interaction of MovedTrack for non e,g
!
      subroutine cintenep
      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
      end     
      subroutine csetknock
      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
      end