caveInteNuc.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	caveintenuc (pj, tgtMassN, avn)

Function/Subroutine Documentation

caveintenuc()

subroutine caveintenuc	(	type(ptcl)	pj,
		integer	tgtMassN,
		real*8	avn
	)

Definition at line 13 of file caveInteNuc.f.

References caaxsec2(), cavenoofsucc(), and cinelx().

Referenced by csampintenuc().

        implicit none

#include  "Zcode.h"
#include  "Zptcl.h"
#include  "Zheavyp.h"

        type(ptcl):: pj   !  projectile heavy
        real*8  avn  !  average no. of interacting nucleons 
        integer tgtmassn  ! target mass # ( # of nucleons)
        
        real*8  xspp, xspa, xsaa, tgta, a3, sucave, pne, tgtz
        integer ihg
!
        tgta = tgtmassn
        tgtz = tgta*0.4   ! almost dummy
        a3 = tgta**0.333333333
        ihg = charge2heavyg(pj%charge)
        pne = pj%fm%p(4)/heavyg2massn(ihg)

!        call cppXsec(pne, xspp)
!        call cxp2xAXsec(tgtA, xspp, xspA)
        call cinelx(pj, tgta,tgtz, xspa)
!                         this is 0.4d0*tgtA is rather dummy
        call caaxsec2(pj, tgta, tgtz,  xsaa)
        avn = heavyg2massn(ihg) * xspa /xsaa
!             this inlcude successive collision inside the target.
!            Cosmos needs the first collision inside the target
!            so divide this by the average number of successive
!            collisions. However, this theory has no firm basis, 
!            sot that we put a switch do use sucave =1
!              get <Nsuc>
        if(howintnuc .eq. 0) then
           call cavenoofsucc(a3, xspp, sucave)
           avn = avn/sucave
        endif

Here is the call graph for this function:

Here is the caller graph for this function: