cprotonAXsec.f

Go to the documentation of this file.

!         p-A, pi-A, k-A  inelastic cross-seciton  etc.
!
!      real*8 a, e, xs, ex
!      a = 14.5
!      e = 2.e0
!      do while (e .lt. 1.e22/1.e9)
!         ex = e - 0.938
!         call cprotonAXsec(a, ex, xs)
!         write(*, *)  sngl(e),sngl(xs)
!         e = e * 10.**0.2
!      enddo
!      end
      subroutine cprotonaxsec(a, e, xs)
!      a: real*8. input. target mass number
!      e: real*8. input. kinetic energy of the ptcl(in GeV)
!     xs: real*8. output. xsection in mb
      implicit none
      real*8 a, e, xs
!
      real*8 xsxp
!
      call cppxsec(e, xsxp)
      call cxp2xaxsec(a, xsxp, xs)
      end
!     --------------------------------
      subroutine cpbaraxsec(a, e, xs)
      implicit none
      real*8 a, e, xs
!
      real*8 xsxp
!
      call cpbarpxsec(e, xsxp)
      call cxp2xaxsec(a, xsxp, xs)
      end
!     ------------------------------------------
      subroutine cpiminusaxsec(a, e, xs)
      implicit none
      real*8 a, e, xs
!
      real*8 xsxp
      
      call cpiminuspxsec(e, xsxp)
      call cxp2xaxsec(a, xsxp, xs)
      end
!     ------------------------------------------
      subroutine cpiplusaxsec(a, e, xs)
      implicit none
      real*8 a, e, xs
!
      real*8 xsxp
      
      call cpipluspxsec(e, xsxp)
      call cxp2xaxsec(a, xsxp, xs)
      end
!     ------------------------------------------
      subroutine ckminusaxsec(a, e, xs)
      implicit none
      real*8 a, e, xs
!
      real*8 xsxp
      
      call ckminuspxsec(e, xsxp)
      call cxp2xaxsec(a, xsxp, xs)
      end
      subroutine ckplusaxsec(a, e, xs)
      implicit none
      real*8 a, e, xs
!
      real*8 xsxp
      
      call ckpluspxsec(e, xsxp)
      call cxp2xaxsec(a, xsxp, xs)
      end