COSMOS v7.655  COSMOSv7655
(AirShowerMC)
cbremLMPXsec.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine cbremlpmxsec (e, rhoin, xs)
 

Function/Subroutine Documentation

◆ cbremlpmxsec()

subroutine cbremlpmxsec ( real*8  e,
real*8  rhoin,
real*8  xs 
)

Definition at line 29 of file cbremLMPXsec.f.

References e10.

Referenced by csampeintl().

29  implicit none
30 !
31 ! input. e: electron/positron energy. in GeV. must be 1.e6<e<1.e13
32 ! i.e., between 10**15 to 10**22 eV.
33 ! input. rhoin: air density in kg/m**3. must be in 1.3 to 3.162x10^-6
34 ! (approx from 0km to 90 km above sea level.
35 !
36 ! output. xs: brems probability /r.l (number of occurence of brems)
37 !
38 ! Note: Input out side of the mensioned range may be safely
39 ! accepted.
40 !
41 
42  real*8 e, rhoin, xs
43  real*8 x, xxx(6), xlog
44  integer i
45  data ( xxx(i), i= 1, 6)/
46  1 2.6807884 , -0.12901903 , 0.27612910e-01,
47  2 -0.30710641e-02, 0.86706897e-04, -0.67901808e-06
48  * /
49  x = e * rhoin * 1.e-3 ! GeV g/cm3
50  if(x .lt. 30.) then
51  xs = 11.447
52  elseif(x .lt. 1.e10) then
53  xs = 0.
54  xlog = log(x)
55  do i = 6, 2, -1
56  xs = (xs + xxx(i))* xlog
57  enddo
58  xs = xs + xxx(1)
59  xs = exp(xs)
60  else
61  xs =0.04198*(x/1.e10)**(-0.5)
62  endif
dE dx *! Nuc Int sampling table e
Definition: cblkMuInt.h:130
nodes i
! 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
*************************block data cblkTracking *************************implicit none data *ExactThick *Freec *RatioToE0 *MagChgDist *TimeStructure *Truncn *Truncx data *IncMuonPolari *KEminObs *ThinSampling *EthinRatio *Generate *LpmEffect *MagPairEmin e10
Definition: cblkTracking.h:9
! structure defining a particle at production ! Basic idea of what is to be contained in ! the particle structue is that dynamical ones should be included those derivable from the particle code ! is not included ******************************************************type fmom momentum sequence union map real e endmap map real * x
Definition: Zptcl.h:21
Here is the caller graph for this function: