cetaDecay.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cetadecay (pj, a, np)

Function/Subroutine Documentation

cetadecay()

subroutine cetadecay	(	type(ptcl), intent(in)	pj,
		type(ptcl), dimension(*), intent(out)	a,
		integer	np
	)

Definition at line 37 of file cetaDecay.f.

References c2bdcy(), cibst1(), cmkptc(), cnbdcy(), d0, kelec, kmuon, kphoton, kpion, and rndc().

Referenced by cinteeta(), and cmydecay().

       implicit none
#include  "Zptcl.h"
#include  "Zcode.h"

       integer np               !output. no. of ptcls produced
!     record /ptcl/ pj         ! input. eta
       type(ptcl),intent(in):: pj 
!     record /ptcl/ a(*)      ! output. produced ptcls
       type(ptcl),intent(out):: a(*)
       integer  i, icon
       real*8 u, w
!
       call rndc(u)
       if(u .lt. .3941d0) then
!           gg
          call cmkptc(kphoton, 0, 0, a(1))
          call cmkptc(kphoton, 0, 0, a(2))
          call c2bdcy(pj, a(1), a(2))
          np=2
       elseif(u .lt. .7209d0 ) then
!           3 pi0
          call cmkptc(kpion, 0, 0, a(1))
          call cmkptc(kpion, 0, 0, a(2))
          call cmkptc(kpion, 0, 0, a(3))
          call cnbdcy(3, pj%mass, a, 0, w, icon)
          np = 3
          do   i=1, np
             call cibst1(i, pj, a(i), a(i))
          enddo
       elseif( u .lt. 0.9501d0 ) then
!           pi+ pi- pi0
          call cmkptc(kpion, 0, 1, a(1))
          call cmkptc(kpion, 0, -1, a(2))
          call cmkptc(kpion, 0, 0, a(3))
          call cnbdcy(3, pj%mass, a, 0, w, icon)
          np = 3
          do   i=1, np
             call cibst1(i, pj, a(i), a(i))
          enddo
       elseif( u .lt. 0.9923d0) then
!     pi+ pi- + g
          call cmkptc(kpion, 0, 1, a(1))
          call cmkptc(kpion, 0, -1, a(2))
          call cmkptc(kphoton, 0, 0, a(3))
          call cnbdcy(3, pj%mass, a, 0, w, icon)
          np = 3
          do   i=1, np
             call cibst1(i, pj, a(i), a(i))
          enddo
       elseif(  u .lt. 0.99261d0) then
!     mu+ + mu- + g
          call cmkptc(kmuon, 0, 1, a(1))
          call cmkptc(kmuon, 0, -1, a(2))
          call cmkptc(kphoton, 0, 0, a(3))
          call cnbdcy(3, pj%mass, a, 0, w, icon)
          np = 3
          do   i=1, np
             call cibst1(i, pj, a(i), a(i))
          enddo
       elseif( u .lt. 0.9926158d0) then
!     mu+ + mu- 
          call cmkptc(kmuon, 0, 1, a(1))
          call cmkptc(kmuon, 0, -1, a(2))
          np = 2
          call c2bdcy(pj, a(1), a(2))
       else
! all rest : regards e+e-g  ~ 6.9e-3
          call cmkptc(kelec, 0, 1, a(1))
          call cmkptc(kelec, 0, -1, a(2))
          call cmkptc(kphoton, 0, 0, a(3))
          call cnbdcy(3, pj%mass, a, 0, w, icon)
          np = 3
          do   i=1, np
             call cibst1(i, pj, a(i), a(i))
          enddo
      endif

Here is the call graph for this function:

Here is the caller graph for this function: