COSMOS v7.655  COSMOSv7655
(AirShowerMC)
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().

11  implicit none
12 ! a: input. real*8 Target Mass No.
13 ! energy: input real*8. gamma energy in GeV
14 ! xs: output. real*8. cross-section in mb
15 !
16  real*8 a, energy, xs
17 
18  real*8 loge, pw
19 
20 !
21 ! gp x-section, xs in mb
22  call cgpxs1(energy, xs)
23 
24  if( a .eq. 1.) then
25 ! nothing to do
26  elseif(energy .lt. 2.) then
27  xs = xs*a*(1.02-.135*energy)
28  elseif(energy .lt. 1.e6 ) then
29 ! before v7.50
30 ! xs = xs*a*0.75 ! this is close to A**0.91 dependence
31 ! pw = around 0.91 and decreasing to 0.82
32 ! (by Bezrukov and Bugaev, Sov.J. of Nuc. phys.
33 ! 33(5), may, 1981
34  loge = log(energy)
35  pw = (-0.00062749*loge + 0.004126)*loge +0.24338+0.6666
36  xs = xs* a**pw
37  else
38 ! extraplation to reach A**2/3 dependence
39  pw = 0.18*(energy/1.e6)**(-0.07) +0.6666
40  xs = xs* a**pw
41  endif
! constants thru Cosmos real ! if multiplied to deg radian Torad ! light velocity m sec ! infinty ! kg m2 *Togpcm2 g cm2 ! g cm2 *Tokgpm2 kg m2 ! cm *Tom m ! m *Tocm cm ! g cm3 *Tokgpm3 kg m3 ! kg m3 *Togpcm3 g cm3 ! sec *Tonsec nsec ! Tesla m ! Avogadro *A2deninv ! mfp *n * xs
Definition: Zglobalc.h:18
real(8), save pw
Definition: csoftenPiK.f:36
real(4), save a
Definition: cNRLAtmos.f:20
subroutine cgpxs1(Eg, xs)
Definition: cgpxs1.f:3
Here is the call graph for this function:
Here is the caller graph for this function: