43 type(
ptcl),
intent(in):: pj
45 type(
ptcl),
intent(out):: a(*)
50 if(u .lt. .3941
d0)
then 54 call c2bdcy(pj, a(1), a(2))
56 elseif(u .lt. .7209
d0 )
then 61 call cnbdcy(3, pj%mass, a, 0, w, icon)
64 call cibst1(i, pj, a(i), a(i))
66 elseif( u .lt. 0.9501
d0 )
then 71 call cnbdcy(3, pj%mass, a, 0, w, icon)
74 call cibst1(i, pj, a(i), a(i))
76 elseif( u .lt. 0.9923
d0)
then 81 call cnbdcy(3, pj%mass, a, 0, w, icon)
84 call cibst1(i, pj, a(i), a(i))
86 elseif( u .lt. 0.99261
d0)
then 91 call cnbdcy(3, pj%mass, a, 0, w, icon)
94 call cibst1(i, pj, a(i), a(i))
96 elseif( u .lt. 0.9926158
d0)
then 101 call c2bdcy(pj, a(1), a(2))
107 call cnbdcy(3, pj%mass, a, 0, w, icon)
110 call cibst1(i, pj, a(i), a(i))
subroutine cibst1(init, p1, p2, po)
max ptcl codes in the kelec
subroutine cnbdcy(n, ecm, p, jw, w, icon)
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
subroutine cetadecay(pj, a, np)
subroutine c2bdcy(p, p1, p2)
subroutine cmkptc(code, subcode, charge, p)
max ptcl codes in the kpion
max ptcl codes in the kmuon