COSMOS v7.655  COSMOSv7655
(AirShowerMC)
Functions/Subroutines
kcombi.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine kcombi (n, m, c)
 

Function/Subroutine Documentation

◆ kcombi()

subroutine kcombi ( real * 8  n,
real * 8  m,
real * 8  c 
)

Definition at line 24 of file kcombi.f.

References d0, and parameter().

24  implicit none
25  real * 8 n, m, c
26 
27  real*8 big, sq2pi
28  parameter(big = 20., sq2pi =2.5066282746310)
29  real*8 ep, x, z1, z2, z, tmp, kgamma
30 
31 !
32  ep(x)=(1.d0/288./x + 1.d0/12.)/x + 1.
33 !
34  z1=n+1.
35  z2=m+1.
36  z=n-m+1.
37  if(n .gt. big .and. (n-m) .gt. big .and. m .gt. big) then
38  tmp=exp( (n+.5)*log( z1/z) + m*log(z/z2) + 1.
39  * -.5 *log(z2) )
40  c=tmp/sq2pi * ep(z1)/ep(z2)/ep(z)
41  elseif(n .gt. big .and. (n-m) .gt. big) then
42  tmp=exp(m * (log(z)-1.) + (n+.5)*log(z1/z) )
43  c=tmp/kgamma(z2)*ep(z1)/ep(z)
44  elseif(n .gt. big .and. m .gt. big) then
45  tmp=exp( (n+.5)*log(z1) - (m+.5)*log(z2) + m-n)
46  c=tmp*ep(z1)/kgamma(z)/ep(z2)
47  else
48  c=kgamma(z1)/kgamma(z2)/kgamma(z)
49  endif
parameter
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
z
nodes z
Definition: cstdatmos0.fNew2.h:93
kgamma
real *8 function kgamma(x)
Definition: kgamma.f:25
d0
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5
m
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data m
Definition: ZlatfitBD.h:35
n
integer n
Definition: Zcinippxc.h:1
x
! 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
c
dE dx *! Nuc Int sampling table c
Definition: cblkMuInt.h:130
Here is the call graph for this function: