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

Go to the source code of this file.

Functions/Subroutines

real *8 function cmbremf1 (up, v)
 
real *8 function cmbremf11 (z)
 

Function/Subroutine Documentation

◆ cmbremf1()

real*8 function cmbremf1 ( real*8  up,
real*8  v 
)

Definition at line 6 of file cmBremF1.f.

References d0.

6  implicit none
7  real*8 up ! input. upsilon
8  real*8 v ! input. faractional energy.
9  real*8 zeta2, cmbremf11
10 
11  zeta2 = v/(1.0d0-v) /3.0d0/up * 2.0d0
12  cmbremf1 = (1.0d0-v)/v * cmbremf11(zeta2)
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data h g *is for param c g data up(2, 1)/7.0d0/
real *8 function cmbremf11(z)
Definition: cmBremF1.f:36
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5
real *8 function cmbremf1(up, v)
Definition: cmBremF1.f:6

◆ cmbremf11()

real*8 function cmbremf11 ( real*8  z)

Definition at line 36 of file cmBremF1.f.

References e, and parameter().

36  implicit none
37  real*8 z
38  integer pol
39  parameter(pol = 13)
40  real*8 coef(pol + 1)
41  real*8 x, y
42  integer i
43  data coef/
44  * -0.4285957, -0.6851291, -0.4785896, -0.1615990, -0.4218448e-01,
45  * -0.8632375e-02, -0.1389391e-02, -0.1832289e-03,
46  * -0.2284316e-04, -0.3119965e-05,
47  * -0.4083642e-06, -0.3924065e-07, -0.2244432e-08, -0.5629939e-10/
48  save coef
49 !
50  if(z .lt. 0.001) then
51  cmbremf11 = 2.1495 * z**0.3333333333 - 1.8138* z
52  elseif(z .gt. 10.) then
53  cmbremf11 = 1.2533 * sqrt(z) * exp(-z)
54  else
55  y = coef(pol + 1)
56  x = log(z)
57  do i = pol, 1, -1
58  y = y * x + coef(i)
59  enddo
60  cmbremf11 = exp(y)
61  endif
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
nodes z
dE dx *! Nuc Int sampling table e
Definition: cblkMuInt.h:130
nodes i
real *8 function cmbremf11(z)
Definition: cmBremF1.f:36
latitude latitude this system is used *****************************************************************! type coord sequence union map real y
Definition: Zcoord.h:25
! 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 call graph for this function: