cinteNuc.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cintenuc
subroutine	cintepion
subroutine	cintekaon
subroutine	cintedmes
subroutine	cintedelta
subroutine	cinteheavy
subroutine	cintemuon
subroutine	csetmuonpol (val)
subroutine	cintennb
subroutine	cinteddb
subroutine	cintesigma
subroutine	cintelambda
subroutine	cintegzai
subroutine	cintebomega
subroutine	cinteeta
subroutine	cinteetap
subroutine	cintelambdac
subroutine	cinterho
subroutine	cinteomega
subroutine	cintephi
subroutine	cinteds
subroutine	cintetau

Function/Subroutine Documentation

cintebomega()

subroutine cintebomega

(

)

Definition at line 328 of file cinteNuc.f.

References cbomegadcy(), and chacol().

Referenced by cintenep().

!     *****************
      use modxsecmedia
      implicit none

#include  "Ztrack.h"
#include  "Ztrackv.h"
!

      if(intinfarray(processno)%process .eq. 'decay') then
         call cbomegadcy(movedtrack%p,  pwork, nproduced)
      else
         call chacol(movedtrack%p, targetnucleonno, targetprotonno,
     *     pwork, nproduced)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cinteddb()

subroutine cinteddb

(

)

Definition at line 256 of file cinteNuc.f.

References cddbdc(), and cerrormsg().

Referenced by cintenep().

!     ********************may be used if 1ry is dd~
      implicit none
      
!----      include 'Ztrack.h'
#include  "Ztrack.h"
!----      include 'Ztrackv.h'
#include  "Ztrackv.h"
!
      character*70  msg

      if(intinfarray(processno)%process .eq. 'decay') then
         call cddbdc(movedtrack%p,  pwork, nproduced)
      else
         write(msg, *) ' in cintennb: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintedelta()

subroutine cintedelta

(

)

Definition at line 132 of file cinteNuc.f.

References cdeltadecay().

Referenced by cintenep().

!     ********************
      use modxsecmedia
      implicit none
#include  "Zcode.h"
#include  "Ztrack.h"
#include  "Ztrackv.h"
      if(intinfarray(processno)%process .eq. 'decay') then
         call cdeltadecay(movedtrack%p,  pwork, nproduced)
      else
         write(0,*)
     *     ' Delta resonance collision should not happen '
         stop
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintedmes()

subroutine cintedmes

(

)

Definition at line 98 of file cinteNuc.f.

References cddecay(), ce2p(), cerrormsg(), chacol(), cmkptc(), k0l, and kkaon.

Referenced by cintenep().

!     ********************
      use modxsecmedia
      implicit none

!----      include 'Ztrack.h'
#include  "Zcode.h"
#include  "Ztrack.h"
!----      include 'Ztrackv.h'
#include  "Ztrackv.h"
!
      character*70 msg
      integer icg
      icg =  movedtrack%p%charge
      if(intinfarray(processno)%process .eq. 'decay') then
         call cddecay(movedtrack%p,  pwork, nproduced)
      elseif(intinfarray(processno)%process .eq. 'coll') then   
!            use kaon 
         if(icg .eq. 0) then
            call cmkptc(kkaon, k0l, 0, movedtrack%p)
         else
            call cmkptc(kkaon, 0, icg, movedtrack%p)
         endif
         call ce2p(movedtrack)
         call chacol(movedtrack%p, targetnucleonno, targetprotonno,
     *     pwork, nproduced)
      else
         write(msg, *) ' in cinteDmeson: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cinteds()

subroutine cinteds

(

)

Definition at line 473 of file cinteNuc.f.

References cdsdecay(), and cerrormsg().

Referenced by cintenep().

!     *******************
      implicit none
      
#include  "Ztrack.h"
#include  "Ztrackv.h"
!
      character*80 msg

      if(intinfarray(processno)%process .eq. 'decay') then
         call cdsdecay(movedtrack%p,  pwork, nproduced)
      elseif(intinfarray(processno)%process .eq. 'coll') then   
         nproduced = 0     ! neglect coll.
      else
         write(msg, *) ' in cinteds: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cinteeta()

subroutine cinteeta

(

)

Definition at line 345 of file cinteNuc.f.

References cetadecay(), and chacol().

Referenced by cintenep().

!     *****************
      use modxsecmedia
      implicit none

#include  "Ztrack.h"
#include  "Ztrackv.h"
!

      if(intinfarray(processno)%process .eq. 'decay') then
         call cetadecay(movedtrack%p,  pwork, nproduced)
      else
         call chacol(movedtrack%p, targetnucleonno, targetprotonno,
     *     pwork, nproduced)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cinteetap()

subroutine cinteetap

(

)

Definition at line 362 of file cinteNuc.f.

References cetapdecay().

Referenced by cintenep().

!     *****************
      use modxsecmedia
      implicit none

#include  "Ztrack.h"
#include  "Ztrackv.h"
!

      if(intinfarray(processno)%process .eq. 'decay') then
         call cetapdecay(movedtrack%p,  pwork, nproduced)
      else
!           negelct col.
         nproduced = 0
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintegzai()

subroutine cintegzai

(

)

Definition at line 311 of file cinteNuc.f.

References cgzaidecay(), and chacol().

Referenced by cintenep().

!     *****************
      use modxsecmedia
      implicit none

#include  "Ztrack.h"
#include  "Ztrackv.h"
!

      if(intinfarray(processno)%process .eq. 'decay') then
         call cgzaidecay(movedtrack%p,  pwork, nproduced)
      else
         call chacol(movedtrack%p, targetnucleonno, targetprotonno,
     *     pwork, nproduced)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cinteheavy()

subroutine cinteheavy

(

)

Definition at line 149 of file cinteNuc.f.

References cerrormsg(), and cheavyint().

Referenced by cintenep().

      use modxsecmedia
!     ********************
      implicit none
#include  "Zcode.h"      
#include  "Ztrack.h"
#include  "Ztrackv.h"
#include  "Ztrackp.h"
#include  "Zevhnp.h"
#include  "Zevhnv.h"

!
      character*70 msg

      if(intinfarray(processno)%process .eq. 'coll') then
         call cheavyint(movedtrack%p, targetnucleonno,
     *         targetprotonno,  pwork, nproduced)
      else
         write(msg, *) ' in cinteHeavy: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintekaon()

subroutine cintekaon

(

)

Definition at line 71 of file cinteNuc.f.

References cerrormsg(), chacol(), and ckaondecay().

Referenced by cintenep().

!     ********************
      use modxsecmedia
      implicit none
!----      include 'Ztrack.h'
#include  "Ztrack.h"
!----      include 'Ztrackv.h'
#include  "Ztrackv.h"
!----      include 'Ztrackp.h'
#include  "Ztrackp.h"
!
      character*70 msg

      if(intinfarray(processno)%process .eq. 'coll') then
         call chacol(movedtrack%p, targetnucleonno, targetprotonno,
     *     pwork, nproduced)
      elseif(intinfarray(processno)%process .eq. 'decay') then
         call ckaondecay(movedtrack%p, incmuonpolari, 
     *           pwork, nproduced, muonpolarization)
      else
         write(msg, *) ' in cinteKaon: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintelambda()

subroutine cintelambda

(

)

Definition at line 294 of file cinteNuc.f.

References chacol(), and clambdadcy().

Referenced by cintenep().

!     ****************
      use modxsecmedia
      implicit none

#include  "Ztrack.h"
#include  "Ztrackv.h"
!

      if(intinfarray(processno)%process .eq. 'decay') then
         call clambdadcy(movedtrack%p,  pwork, nproduced)
      else
         call chacol(movedtrack%p, targetnucleonno, targetprotonno,
     *     pwork, nproduced)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintelambdac()

subroutine cintelambdac

(

)

Definition at line 379 of file cinteNuc.f.

References ce2p(), cerrormsg(), chacol(), clambdacdcy(), cmkptc(), and knuc.

Referenced by cintenep().

!     ********************
      use modxsecmedia
      implicit none
      
#include  "Ztrack.h"
#include  "Ztrackv.h"
#include  "Zcode.h"
!
      character*80 msg

      if(intinfarray(processno)%process .eq. 'decay') then
         call clambdacdcy(movedtrack%p,  pwork, nproduced)
      elseif(intinfarray(processno)%process .eq. 'coll') then   
!            use proton or p-bar
         if( movedtrack%p%charge > 0 ) then
            call cmkptc(knuc, -1, 1, movedtrack%p)
         else
            call cmkptc(knuc, 1, -1, movedtrack%p)
         endif
         call ce2p(movedtrack)
         call chacol(movedtrack%p, targetnucleonno, targetprotonno,
     *     pwork, nproduced)
      else
         write(msg, *) ' in cinteLambdac: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintemuon()

subroutine cintemuon

(

)

Definition at line 174 of file cinteNuc.f.

References ccapnu(), cfixtarget(), cgetcaprate(), cmuinte(), cmuneudcy(), rndc(), and t0mu.

Referenced by cintenep().

      use modxsecmedia
!     ********************
      implicit none

!----      include 'Ztrack.h'
#include  "Ztrack.h"
!----      include 'Ztrackv.h'
#include  "Ztrackv.h"
!----      include 'Zlife.h'
#include  "Zlife.h"
#include  "Ztrackp.h"
      real*8  capr, t0cap, u1, u2

!
      if(intinfarray(processno)%process .eq. 'decay') then
!           see if negative muon is being captured
         if(movedtrack%p%charge .eq. -1) then
            if(movedtrack%p%fm%p(4) .le. movedtrack%p%mass*1.001) then
!                   assume stopping muon
               call cgetcaprate( media(mediumno) )  ! v7.640
!               call cmucap(TargetNucleonNo, TargetProtonNo, capr)
               capr = media(mediumno)%xs
               t0cap = 1./capr
               call rndc(u1)
               call rndc(u2)
               if( - log(u1)*t0mu .gt. - log(u2)*t0cap) then
!                       capture
                  intinfarray(processno)%process = "capt"
                  call cfixtarget( media(mediumno) )
                  call ccapnu(targetnucleonno, targetprotonno,
     *               pwork, nproduced)
               else
!!!!              MuonPolarization = 0.   ! stopping mu has no pol.
!                                          keep the value ;it's better.
                  call cmuneudcy(movedtrack%p, muonpolarization, 
     *            eabsorb(1),   pwork, nproduced)
               endif
            else
               call cmuneudcy(movedtrack%p, muonpolarization,
     *            eabsorb(1),   pwork, nproduced)
            endif
         else
            call cmuneudcy(movedtrack%p, muonpolarization,
     *            eabsorb(1),   pwork, nproduced)
         endif
      else
         call cmuinte
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintennb()

subroutine cintennb

(

)

Definition at line 235 of file cinteNuc.f.

References cerrormsg(), and cnnbdc().

Referenced by cintenep().

!     ********************may be used if 1ry is nn~
      implicit none
      
!----      include 'Ztrack.h'
#include  "Ztrack.h"
!----      include 'Ztrackv.h'
#include  "Ztrackv.h"
!
      character*70 msg

      if(intinfarray(processno)%process .eq. 'decay') then
         call cnnbdc(movedtrack%p,  pwork, nproduced)
      else
         write(msg, *) ' in cintennb: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call  cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintenuc()

subroutine cintenuc

(

)

Definition at line 2 of file cinteNuc.f.

References cerrormsg(), and chacol().

Referenced by cintenep().

      use modxsecmedia
      implicit none


#include  "Ztrack.h"
#include  "Ztrackv.h"

!          cinteNuc.  Treat Nucelon interactions with air target.
      character*70 msg
!////////////
      logical show
      common /showshow/ show
!/////////////
      if(intinfarray(processno)%process .eq. 'coll') then
!/////////////
         if(show) then
            write(0,*) 'entering chAcol'
         endif
!//////////////
         call chacol(movedtrack%p, targetnucleonno, targetprotonno,
     *     pwork, nproduced)
!/////////
         if(show) then
            write(0,*) ' aft chAcol: Npro=',nproduced
         endif
!////////////////
      else
         write(msg, *) ' in cinteNuc: process=',
     *                intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call cerrormsg(msg,0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cinteomega()

subroutine cinteomega

(

)

Definition at line 431 of file cinteNuc.f.

References cerrormsg(), and comgdc().

Referenced by cintenep().

!     *******************
      implicit none
      
#include  "Ztrack.h"
#include  "Ztrackv.h"
!
      character*80 msg

      if(intinfarray(processno)%process .eq. 'decay') then
         call comgdc(movedtrack%p,  pwork, nproduced)
      elseif(intinfarray(processno)%process .eq. 'coll') then   
         nproduced = 0     ! neglect coll.
      else
         write(msg, *) ' in cinteomega: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintephi()

subroutine cintephi

(

)

Definition at line 452 of file cinteNuc.f.

References cerrormsg(), and cphidc().

Referenced by cintenep().

!     *******************
      implicit none
      
#include  "Ztrack.h"
#include  "Ztrackv.h"
!
      character*80 msg

      if(intinfarray(processno)%process .eq. 'decay') then
         call cphidc(movedtrack%p,  pwork, nproduced)
      elseif(intinfarray(processno)%process .eq. 'coll') then   
         nproduced = 0     ! neglect coll.
      else
         write(msg, *) ' in cintephi: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintepion()

subroutine cintepion

(

)

Definition at line 37 of file cinteNuc.f.

References cerrormsg(), chacol(), cpi0decay(), and cpimudecay().

Referenced by cintenep().

!     ********************
      use modxsecmedia
      implicit none
      
#include  "Ztrack.h"
#include  "Ztrackv.h"
#include  "Ztrackp.h"

!
      character*70 msg

      if(intinfarray(processno)%process .eq. 'coll') then
         call chacol(movedtrack%p, targetnucleonno, targetprotonno,
     *     pwork, nproduced)

      elseif(intinfarray(processno)%process .eq. 'decay') then
         if(movedtrack%p%charge .eq. 0) then
            call cpi0decay(movedtrack%p, pwork, nproduced)
         else
            call cpimudecay(movedtrack%p, incmuonpolari, 
     *           pwork, nproduced, muonpolarization)
         endif
      else
         write(msg, *) ' in cintePion: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         write(0,*) ' Et=', movedtrack%p%fm%p(4), ' chg=',
     *                movedtrack%p%charge
         call cerrormsg(msg,0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cinterho()

subroutine cinterho

(

)

Definition at line 410 of file cinteNuc.f.

References cerrormsg(), and crhodc().

Referenced by cintenep().

!     *******************
      implicit none
      
#include  "Ztrack.h"
#include  "Ztrackv.h"
!
      character*80 msg

      if(intinfarray(processno)%process .eq. 'decay') then
         call crhodc(movedtrack%p,  pwork, nproduced)
      elseif(intinfarray(processno)%process .eq. 'coll') then   
         nproduced = 0     ! neglect coll.
      else
         write(msg, *) ' in cinterho: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintesigma()

subroutine cintesigma

(

)

Definition at line 277 of file cinteNuc.f.

References chacol(), and csigmadecay().

Referenced by cintenep().

!     *****************
      use modxsecmedia
      implicit none

#include  "Ztrack.h"
#include  "Ztrackv.h"
!

      if(intinfarray(processno)%process .eq. 'decay') then
         call csigmadecay(movedtrack%p,  pwork, nproduced)
      else
         call chacol(movedtrack%p, targetnucleonno, targetprotonno,
     *     pwork, nproduced)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cintetau()

subroutine cintetau

(

)

Definition at line 495 of file cinteNuc.f.

References cerrormsg(), and ctauneudcy().

Referenced by cintenep().

!     *******************
      implicit none
      
#include  "Ztrack.h"
#include  "Ztrackv.h"
!
      character*80 msg

      if(intinfarray(processno)%process .eq. 'decay') then
         call ctauneudcy(movedtrack%p,  pwork, nproduced)
      elseif(intinfarray(processno)%process .eq. 'coll') then   
         nproduced = 0     ! neglect coll.
      else
         write(msg, *) ' in cintetau: process=',
     *           intinfarray(processno)%process,
     *               ' undef. ProcessNo=',processno
         call cerrormsg(msg, 0)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

csetmuonpol()

subroutine csetmuonpol

(

real*8

val

)

Definition at line 224 of file cinteNuc.f.

References d0.

Referenced by cconv_prim_e(), and xbgrun().

!         This may be used when you set muon as a primary for which
!        polarization is not fixed.
!
      implicit none
#include  "Ztrack.h"
#include "Ztrackv.h"
      real*8 val
      muonpolarization =min(1.d0, max(val, -1.d0 ))

Here is the caller graph for this function: