Go to the source code of this file.
|
subroutine | ksamppw (ini, coef, power, node, n, x, fx) |
|
subroutine | ksamppw2 (coef, power, xp, cum, intgv, np, x, fx) |
|
subroutine | ksamppwx (c, p, x0, n, xp, c2) |
|
subroutine | ksamppwmrg (c1, p1, x1, n1, c2, p2, x2, n2, c, p, x, n) |
|
◆ ksamppw()
Definition at line 95 of file ksampPw.f.
References d0, parameter(), and rndc().
Referenced by cmbreme1(), cmbreme2(), and main().
122 real*8 cum(*), intgv(*)
134 if(power(
i) .ne. 1.
d0)
then 135 inte(
i) = coef(
i)/(1.0-power(
i))*
136 * (node(
i+1)**(1.-power(
i)) - node(
i)**(1.-power(
i)))
138 inte(
i) = coef(
i) * log(node(
i+1)/node(
i))
147 ci(
i) = ci(
i) + ci(
i-1)
158 if(u .le. ci(
i))
then 169 if(power(
j) .ne. 1.
d0)
then 170 temp = node(
j+1)**(1.-power(
j)) -
171 * inte(
j)*u/( coef(
j)/(1. - power(
j)) )
173 x =
temp**(1./(1.-power(
j)))
175 x = (node(
j+1)/node(
j))**u * node(
j)
177 fx = coef(
j) *
x **(-power(
j))
183 entry ksamppwq(cum, intgv, np)
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
real(4), dimension(:), allocatable, save temp
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
! 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
◆ ksamppw2()
subroutine ksamppw2 |
( |
real*8, dimension(np) |
coef, |
|
|
real*8, dimension(np) |
power, |
|
|
real*8, dimension(np+1) |
xp, |
|
|
real*8, dimension(np) |
cum, |
|
|
real*8, dimension(np) |
intgv, |
|
|
integer |
np, |
|
|
real*8 |
x, |
|
|
real*8 |
fx |
|
) |
| |
Definition at line 199 of file ksampPw.f.
References rndc().
221 if(u .le. cum(
i))
then 230 temp = xp(
j+1)**(1.-power(
j)) -
231 * intgv(
j)*u/( coef(
j)/(1. - power(
j)) )
233 x =
temp**(1./(1.-power(
j)))
234 fx = coef(
j) *
x**(-power(
j))
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
real(4), dimension(:), allocatable, save temp
! 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
◆ ksamppwmrg()
subroutine ksamppwmrg |
( |
real*8, dimension(n1) |
c1, |
|
|
real*8, dimension(n1) |
p1, |
|
|
real*8, dimension(n1+1) |
x1, |
|
|
integer |
n1, |
|
|
real*8, dimension(n2) |
c2, |
|
|
real*8, dimension(n2) |
p2, |
|
|
real*8, dimension(n2+1) |
x2, |
|
|
integer |
n2, |
|
|
real*8, dimension(*) |
c, |
|
|
real*8, dimension(*) |
p, |
|
|
real*8, dimension(*) |
x, |
|
|
integer |
n |
|
) |
| |
Definition at line 273 of file ksampPw.f.
Referenced by cmbreme1(), and cmbreme2().
297 x(1) = max(
x1(1),
x2(1))
298 do while (
i .le. n1 .and.
j .le. n2)
299 if(
x1(
i+1) .lt.
x2(
j+1) )
then 306 elseif(
x1(
i+1) .eq.
x2(
j+1))
then block data include Zlatfit h c fitting region data x1(1)/0.03/
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
block data include Zlatfit h c fitting region data x2(1)/0.5/data x1(2)/0.3/
! 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
dE dx *! Nuc Int sampling table c
◆ ksamppwx()
subroutine ksamppwx |
( |
real*8, dimension(n) |
c, |
|
|
real*8, dimension(n) |
p, |
|
|
real*8, dimension(n) |
x0, |
|
|
integer |
n, |
|
|
real*8, dimension(n-1) |
xp, |
|
|
real*8, dimension(n) |
c2 |
|
) |
| |
Definition at line 240 of file ksampPw.f.
Referenced by cmbreme1(), cmbreme2(), and main().
261 * (c2(
i)/c2(
i+1) )**(1./(
p(
i) -
p(
i+1)))
********************block data cblkHeavy ********************integer j data *HeavyG2symbol p
dE dx *! Nuc Int sampling table c