csampCollInA.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	csampcollina (pj, ia, nc)

Function/Subroutine Documentation

csampcollina()

subroutine csampcollina	(	type(ptcl)	pj,
		integer	ia,
		integer	nc
	)

Definition at line 28 of file csampCollInA.f.

References cinelx(), d0, parameter(), and rndc().

Referenced by chacoladhoc().

        implicit none

#include  "Zcode.h"
#include  "Zptcl.h"
        integer i
        type(ptcl):: pj  !input projectile particle
        integer ia       ! target mass no.
        integer nc       ! output number of collisions.

        real*8 a1, a2, a3, da
        real*8  a, xs
!
        integer  idxa, idxxs

        real*8 u
        real*8  xs1, xs2, dxs
        integer mm, nn, kk
        parameter( mm = 14, nn = 16,  kk= 16)
!                 cumProb,  xsec     A

        real*4 cumprob(mm, nn, kk)




!        parameter (a1 = 4.0**0.3333333, a2 = 208.**0.333333333333,
        parameter(a1 = 1.5874011, a2 = 5.9249921,
     *  da = (a2-a1)/(kk-1) )

!         parameter ( xs1 =log10(15.), xs2 =log10(300.) )
        parameter( xs1 = 1.176091259, xs2 = 2.477121255, 
     *  dxs =( xs2 - xs1)/15.)

#include  "cumProb.h"
!             get cross-section for proton target.        
!        call cxpXsec(pj, xs)
        call cinelx(pj, 1.d0, 1.d0,  xs)
        xs = log10(xs)
!
         a = ia
         a3 = a**0.3333333333
         idxa = (a3- a1)/da + 1
         idxxs = (xs - xs1)/dxs + 1
         if( (a3 - idxa * da - a1) .gt. (idxa*da + da + a1 - a3) ) then
            idxa = idxa  + 1
         endif
         idxa =max(1, min(idxa, kk))

         if( (xs - idxxs * dxs - xs1) .gt.
     *        (idxxs * dxs + dxs + xs1 - xs) ) then
            idxxs = idxxs +1
         endif
         idxxs =max(1, min(idxxs, nn))
         call rndc(u)
         do   i=1,  mm
            if(u .le. cumprob(i, idxxs, idxa) ) then
               nc=i
               goto 100
            endif
         enddo
  100    continue

Here is the call graph for this function:

Here is the caller graph for this function: