cppXsec0.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cppxsec0 (e, xs)
real *8 function	ctotppx (p, logp)
real *8 function	celappx (p, logp)
real *8 function	cxsec (p, logp, a, b, n, c, d)
subroutine	cpbarpxsec0 (e, xs)
real *8 function	ctotpbpx (p, logp)
real *8 function	celapbpx (p, logp)
subroutine	cpiminuspxsec0 (e, xs)
real *8 function	ctotpimpx (p, logp)
real *8 function	celapimpx (p, logp)
subroutine	cpipluspxsec0 (e, xs)
real *8 function	ctotpippx (p, logp)
real *8 function	celapippx (p, logp)
subroutine	ckminuspxsec0 (e, xs)
real *8 function	ctotkmpx (p, logp)
real *8 function	celakmpx (p, logp)
subroutine	ckpluspxsec0 (e, xs)
real *8 function	ctotkppx (p, logp)
real *8 function	celakppx (p, logp)
subroutine	cpimairreacx (e, xs)
subroutine	cpipairreacx (e, xs)
subroutine	cpimairabsx (e, xs)
subroutine	cpipairabsx (e, xs)

Function/Subroutine Documentation

celakmpx()

real*8 function celakmpx	(	real*8	p,
		real*8	logp
	)

Definition at line 303 of file cppXsec0.f.

      implicit none
        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/7.3/, b/0./, n/0./, c/0.29/, d/-2.4/
        celakmpx = cxsec(p, logp, a,  b, n, c, d)

celakppx()

real*8 function celakppx	(	real*8	p,
		real*8	logp
	)

Definition at line 360 of file cppXsec0.f.

      implicit none
        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/5.0/, b/8.1/, n/-1.8/, c/0.16/, d/-1.3/
        celakppx = cxsec(p, logp, a,  b, n, c, d)

celapbpx()

real*8 function celapbpx	(	real*8	p,
		real*8	logp
	)

Definition at line 128 of file cppXsec0.f.

      implicit none
        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/10.2/, b/52.7/, n/-1.16/, c/0.125/, d/-1.28/
        celapbpx = cxsec(p, logp, a,  b, n, c, d)

celapimpx()

real*8 function celapimpx	(	real*8	p,
		real*8	logp
	)

Definition at line 186 of file cppXsec0.f.

      implicit none
        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/1.76/, b/11.2/, n/-0.64/, c/0.043/, d/0./
        celapimpx = cxsec(p, logp, a,  b, n, c, d)

celapippx()

real*8 function celapippx	(	real*8	p,
		real*8	logp
	)

Definition at line 242 of file cppXsec0.f.

      implicit none
        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/0./, b/11.4/, n/-0.4/, c/0.079/, d/0./
        celapippx = cxsec(p, logp, a,  b, n, c, d)

celappx()

real*8 function celappx	(	real*8	p,
		real*8	logp
	)

Definition at line 72 of file cppXsec0.f.

      implicit none
        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/11.9/, b/26.9/, n/-1.21/, c/0.169/, d/-1.85/
        celappx = cxsec(p, logp, a, b, n,  c, d)

ckminuspxsec0()

subroutine ckminuspxsec0	(	real*8	e,
		real*8	xs
	)

Definition at line 263 of file cppXsec0.f.

References d, and kfrge().

        implicit none
        real*8 e, xs
        real*8 ee
        integer i, icon, l

        real*8 lp
        real*8 ctotkmpx, celakmpx
        real*8 ea(23), xsa(23)

!
      data (ea(i),i=   1,  23)/
     1 0.040, 0.070, 0.080, 0.085, 0.090, 0.100, 0.120, 0.130, 0.150,
     2 0.170, 0.200, 0.250, 0.300, 0.400, 0.500, 0.600, 0.640, 1.500,
     3 5.000,10.000,20.000,50.000, 100.0001/
      data (xsa(i),i=   1,  23)/
     1 1.d-5,30.000,38.000,39.800,40.000,40.200,40.000,39.000,35.000,
     230.000,22.500,17.500,18.000,20.000,22.500,26.000,27.000,25.000,
     321.000,20.000,19.000,18.500,18.0000/
!
       if(e .le. 40.e-3)then
           xs=1.e-5
       elseif(e .lt. 100.) then
           ee = e
           call kfrge(ea, 1, 23, ee, l, icon)
           xs=(xsa(l)-xsa(l-1))/(ea(l)-ea(l-1))*(ee-ea(l-1)) +
     *         xsa(l-1)
       else
           lp = log(e)   !  e = p here.
           xs = ctotkmpx(e, lp) - celakmpx(e, lp)
       endif

Here is the call graph for this function:

ckpluspxsec0()

subroutine ckpluspxsec0	(	real*8	e,
		real*8	xs
	)

Definition at line 324 of file cppXsec0.f.

References d, and kfrge().

        implicit none
        real*8 e, xs
        integer i, icon, l
        real*8 ea(13), xsa(13), ee
        real*8 lp
        real*8 ctotkppx, celakppx


      data (ea(i),i=   1,  13)/
     1 0.100, 0.200, 0.300, 0.500, 1.000, 2.000, 3.000, 4.000, 5.000,
     210.000,20.000,50.000, 100.0001 /
      data (xsa(i),i=   1,  13)/
     1 1.d-5, 2.500, 5.000, 7.500,10.000,12.000,13.000,13.500,14.000,
     2 15.000,15.300,15.800,16.200/

        if(e .le. 100.e-3)then
           xs=1.d-5
        elseif(e .lt. 100.) then
           ee = e
           call kfrge(ea, 1, 13, ee, l, icon)
           xs=(xsa(l)-xsa(l-1))/(ea(l)-ea(l-1))*(ee-ea(l-1)) +
     *         xsa(l-1)
        else
           lp = log(e)   !  e = p here.
           xs = ctotkppx(e, lp) - celakppx(e, lp)
        endif

Here is the call graph for this function:

cpbarpxsec0()

subroutine cpbarpxsec0	(	real*8	e,
		real*8	xs
	)

Definition at line 104 of file cppXsec0.f.

        implicit none
        real*8 e, xs
        real*8 lp
        real*8 ctotpbpx, celapbpx


!               due to the formula given in ptcl data 88. p121

        if(e .lt. 100.) then
           xs =  (e/100)**(-.55)*3 +32
        else
           lp = log(e)   !  e = p here.
           xs =  ctotpbpx(e, lp) - celapbpx(e, lp)
        endif

cpimairabsx()

subroutine cpimairabsx	(	real*8	e,
		real*8	xs
	)

Definition at line 446 of file cppXsec0.f.

References d0, d3, and kfrge().

        implicit none
        real*8 e, xs
        integer i, icon, l
        real*8 ea(13), xsa(13), ee
!           mev
      data (ea(i),i=   1,  13)/
     1   0.,  10.,  20.,  30.,  40.,  50.,  70., 100., 140., 180., 200.,
     2 300., 500.001/
      data (xsa(i),i=   1,  13)/
     1   0.,  91., 136., 159., 163., 172., 184., 193., 195., 193., 184.,
     2 115.,   0./
           ee=min(500.d0, e*1.d3)
           call kfrge(ea, 1, 13, ee, l, icon)
           xs=(xsa(l)-xsa(l-1))/(ea(l)-ea(l-1))*(ee-ea(l-1)) +
     *         xsa(l-1)

Here is the call graph for this function:

cpimairreacx()

subroutine cpimairreacx	(	real*8	e,
		real*8	xs
	)

Definition at line 381 of file cppXsec0.f.

References d0, d3, and kfrge().

        implicit none
        real*8 e, xs
        integer i, icon, l
        real*8 ea(23), xsa(23), ee


!            MeV
      data (ea(i),i=   1,  23)/
     1   0.,  10.,  20.,  30.,  40.,  50.,  70.,  80., 100., 120., 140.,
     2 160., 180., 200., 240., 260., 280., 300., 340., 360., 400., 450.,
     3 500.0001/
      data (xsa(i),i=   1,  23)/
     1   0.,  92., 138., 172., 205., 230., 356., 414., 452., 460., 462.,
     2 463., 461., 460., 416., 391., 359., 328., 276., 241., 224., 218.,
     3 207./
           ee=min(500.d0, e*1.d3)
           call kfrge(ea, 1, 23, ee, l, icon)
           xs=(xsa(l)-xsa(l-1))/(ea(l)-ea(l-1))*(ee-ea(l-1)) +
     *         xsa(l-1)

Here is the call graph for this function:

cpiminuspxsec0()

subroutine cpiminuspxsec0	(	real*8	e,
		real*8	xs
	)

Definition at line 148 of file cppXsec0.f.

References d, and kfrge().

        implicit none
        real*8 e, xs
        integer i, l, icon
        real*8 lp
        real*8 ctotpimpx, celapimpx

        real*8 ea(25), xsa(25), ee

      data (ea(i),i=   1,  25)/
     1 0.010, 0.020, 0.030, 0.060, 0.100, 0.200, 0.300, 0.400, 0.450,
     2 0.500, 0.550, 0.670, 0.870, 1.000, 1.200, 1.350, 1.500, 2.000,
     3 2.800, 3.000, 5.000, 7.000,10.000,50.000, 100.001 /
      data (xsa(i),i=   1,  25)/
     1 0.000, 4.000, 6.000,10.000,15.000,15.300,16.500,18.000,19.500,
     222.500,30.000,26.000,33.500,32.000,25.000,23.500,25.000,27.000,
     326.000,24.500,23.500,22.500,22.500,21.500,21.000/
        if(e .le. 10.d-3) then
           xs=1.d-5
        elseif(e .lt. 100.) then
           ee= e
           call kfrge(ea, 1, 25, ee, l, icon)
           xs=(xsa(l)-xsa(l-1))/(ea(l)-ea(l-1))*(ee-ea(l-1)) +
     *         xsa(l-1)
        else
           lp = log(e)   !  e = p here.
           xs = ctotpimpx(e, lp) - celapimpx(e, lp)
        endif

Here is the call graph for this function:

cpipairabsx()

subroutine cpipairabsx	(	real *8	e,
		real *8	xs
	)

Definition at line 474 of file cppXsec0.f.

References d0, d3, and kfrge().

        implicit none
        real *8 e, xs
!
        integer i, icon, l
        real * 8 ea(13), xsa(13), ee
!           mev
      data (ea(i),i=   1,  13)/
     1   0.,  10.,  20.,  30.,  40.,  50.,  70., 100., 140., 180., 200.,
     2 300., 500.001/
      data (xsa(i),i=   1,  13)/
     1   0.,  47.,  92., 115., 129., 140., 161., 178., 184., 184., 182.,
     2 115.,   0./
           ee=min(500.d0, e*1.d3)
           call kfrge(ea, 1, 13, ee, l, icon)
           xs=(xsa(l)-xsa(l-1))/(ea(l)-ea(l-1))*(ee-ea(l-1)) +
     *         xsa(l-1)

Here is the call graph for this function:

cpipairreacx()

subroutine cpipairreacx	(	real*8	e,
		real*8	xs
	)

Definition at line 415 of file cppXsec0.f.

References d0, d3, and kfrge().

        implicit none
        real*8 e, xs

        integer i, icon, l
        real*8 ea(24), xsa(24), ee
!           mev
      data (ea(i),i=   1,  24)/
     1   0.,  10.,  20.,  30.,  40.,  50.,  70.,  80.,  90., 100., 120.,
     2 140., 160., 180., 200., 240., 260., 280., 300., 340., 360., 400.,
     3 450., 500.0001/
      data (xsa(i),i=   1,  24)/
     1   0.,  57.,  94., 126., 155., 195., 316., 379., 408., 416., 437.,
     2 443., 448., 443., 437., 393., 370., 345., 328., 253., 230., 207.,
     3 190., 190./
           ee=min(500.d0, e*1.d3)
           call kfrge(ea, 1, 24, ee, l, icon)
           xs=(xsa(l)-xsa(l-1))/(ea(l)-ea(l-1))*(ee-ea(l-1)) +
     *         xsa(l-1)

Here is the call graph for this function:

cpipluspxsec0()

subroutine cpipluspxsec0	(	real*8	e,
		real*8	xs
	)

Definition at line 207 of file cppXsec0.f.

References d, and kfrge().

        implicit none 
        real*8 e, xs
        
        real*8 ea(18), xsa(18), ee
        integer i, icon, l
        real*8 lp, ctotpippx, celapippx
!
      data (ea(i),i=   1,  18)/
     1 0.100, 0.300, 0.500, 0.600, 0.700, 0.800, 0.900, 1.000, 1.100,
     2 1.200, 1.300, 1.500, 5.000,10.000,30.000,50.000,70.000, 100.001/
      data (xsa(i),i=   1,  18)/
     1 1.d-5, 4.500, 7.500, 8.000,10.500,12.000,13.500,16.000,18.000,
     2 19.500,19.800,20.000,20.200,19.500,19.200,19.200,19.500,20.000/
!
        if(e .lt. 100.d-3) then
           xs=1.d-5
        elseif(e .lt. 100.) then
           ee= e
           call kfrge(ea, 1, 18, ee, l, icon)
           xs=(xsa(l)-xsa(l-1))/(ea(l)-ea(l-1))*(ee-ea(l-1)) +
     *         xsa(l-1)
        else
           lp = log(e)   !  e = p here.
           xs = ctotpippx(e, lp) - celapippx(e, lp)
        endif

Here is the call graph for this function:

cppxsec0()

subroutine cppxsec0	(	real*8	e,
		real*8	xs
	)

Definition at line 30 of file cppXsec0.f.

References d, d0, and kfrge().

        implicit none
        real*8 e, xs
!
        real*8 ee
        integer i, icon, l
        real*8 lp
        real*8 ctotppx, celappx

        real*8 ea(21), xsa(21)
!           ea in GeV
      data (ea(i),i=   1,  21)/
     1 0.301, 0.350, 0.400, 0.450, 0.500, 0.550, 0.600, 0.650, 0.700,
     2 0.800, 0.900, 1.000, 1.500, 2.000, 5.000, 7.000,10.000,20.000,
     3 40.000,50.000, 100.001/
      data (xsa(i),i=   1,  21)/
     1 0.000, 0.500, 1.000, 2.500, 5.000, 7.350,11.000,15.000,18.500,
     222.000,23.000,23.500,25.000,25.500,28.000,29.000,29.500,30.000,
     330.000,29.500,31.000/
      if(e .le. 0.)  then
         xs = 0.

      elseif(e .lt. .30101d0) then
         xs=1.d-5
      elseif(e .lt. 100.) then
         ee = e
         call kfrge(ea, 1, 21, ee, l, icon)
         xs=(xsa(l)-xsa(l-1))/(ea(l)-ea(l-1))*(ee-ea(l-1)) +
     *        xsa(l-1)
      else
         lp = log(e)            ! e = p      
         xs = ctotppx(e, lp) - celappx(e, lp)
      endif

Here is the call graph for this function:

ctotkmpx()

real*8 function ctotkmpx	(	real*8	p,
		real*8	logp
	)

Definition at line 295 of file cppXsec0.f.

      implicit none
        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/32.1/, b/0./, n/0./, c/0.66/, d/-5.6/
        ctotkmpx = cxsec(p, logp, a,  b, n, c, d)

ctotkppx()

real*8 function ctotkppx	(	real*8	p,
		real*8	logp
	)

Definition at line 352 of file cppXsec0.f.

      implicit none
        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/18.1/, b/0./, n/0./, c/0.26/, d/-1.0/
        ctotkppx = cxsec(p, logp, a,  b, n, c, d)

ctotpbpx()

real*8 function ctotpbpx	(	real*8	p,
		real*8	logp
	)

Definition at line 120 of file cppXsec0.f.

      implicit none
        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/38.4/, b/77.6/, n/-0.64/, c/0.26/,  d/-1.2/
        ctotpbpx = cxsec(p, logp, a,  b, n, c, d)

ctotpimpx()

real*8 function ctotpimpx	(	real*8	p,
		real*8	logp
	)

Definition at line 177 of file cppXsec0.f.

      implicit none

        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/33./, b/14./, n/-1.36/, c/0.456/, d/-4.03/
        ctotpimpx = cxsec(p, logp, a,  b, n, c, d)

ctotpippx()

real*8 function ctotpippx	(	real*8	p,
		real*8	logp
	)

Definition at line 234 of file cppXsec0.f.

      implicit none
        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/16.4/, b/19.3/, n/-0.42/, c/0.19/, d/0./
        ctotpippx = cxsec(p, logp, a,  b, n, c, d)

ctotppx()

real*8 function ctotppx	(	real*8	p,
		real*8	logp
	)

Definition at line 64 of file cppXsec0.f.

      implicit none
        real*8 p, logp
        real*8 a, b, n,  c, d
        real*8 cxsec
        data a/48./, b/0./, n/0./, c/.522/, d/-4.51/
        ctotppx = cxsec(p, logp, a,  b, n, c, d)

cxsec()

real*8 function cxsec	(	real*8	p,
		real*8	logp,
		real*8	a,
		real*8	b,
		real*8	n,
		real*8	c,
		real*8	d
	)

Definition at line 83 of file cppXsec0.f.

!        elastic/inelastic  high energy  cross-secion formula.
!        a + bp**n + c*log(p)**2 + d*log(p) 
!        p is momentum in GeV/c
      implicit none
      real*8 p, logp, a, b, n, c, d
!
      cxsec = (c*logp + d)*logp + b*p**n + a