cgpXsec.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cgpxsec (a, energy, xs)

Function/Subroutine Documentation

cgpxsec()

subroutine cgpxsec	(	real*8	a,
		real*8	energy,
		real*8	xs
	)

Definition at line 11 of file cgpXsec.f.

References cgpxs1().

Referenced by cgetphotopxs(), cghpath(), cinelx0(), and epgencolbyphoton().

      implicit none
!       a: input. real*8  Target Mass No.
!     energy: input real*8.  gamma energy in GeV
!     xs:  output. real*8.   cross-section in mb
!
      real*8 a, energy, xs

      real*8 loge, pw

!
!           gp x-section,   xs in mb
      call cgpxs1(energy, xs)

      if( a .eq. 1.) then
!          nothing to do
      elseif(energy .lt. 2.) then    
         xs = xs*a*(1.02-.135*energy)
      elseif(energy .lt. 1.e6 ) then
!          before v7.50
!         xs = xs*a*0.75    ! this is close to A**0.91 dependence
!          pw = around 0.91 and decreasing to 0.82 
!             (by Bezrukov and Bugaev, Sov.J. of Nuc. phys.
!              33(5), may, 1981 
         loge = log(energy)
         pw = (-0.00062749*loge + 0.004126)*loge +0.24338+0.6666 
         xs = xs* a**pw
      else
!            extraplation to reach A**2/3 dependence
         pw = 0.18*(energy/1.e6)**(-0.07) +0.6666
         xs = xs* a**pw
      endif

Here is the call graph for this function:

Here is the caller graph for this function: