cetapDecay.f File Reference

Go to the source code of this file.

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

Function/Subroutine Documentation

cetapdecay()

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

Definition at line 33 of file cetapDecay.f.

References c2bdcy(), cibst1(), cmkptc(), cnbdcy(), d0, keta, komega, kphoton, kpion, and rndc().

Referenced by cinteetap().

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

       integer np               !output. no. of ptcls produced
       type(ptcl),intent(in):: pj 
       type(ptcl),intent(out):: a(*)
       integer  i, icon
       real*8 u, w
!
       call rndc(u)
       if(u .lt. 0.426d0) then
!           pi+ pi- eta
          call cmkptc(kpion, -1, 1, a(1))
          call cmkptc(kpion, 1, -1, a(2))
          call cmkptc(keta, 0, 0,  a(3) )
          call cnbdcy(3, pj%mass, a, 0, w, icon)
          np=3
       elseif(u .lt. .715d0 ) then
!            pi+ pi- gamma
          call cmkptc(kpion, -1, 1, a(1))
          call cmkptc(kpion, 1, -1, a(2))
          call cmkptc(kphoton, 0, 0, a(3))
          call cnbdcy(3, pj%mass, a, 0, w, icon)
          np = 3
       elseif( u .lt. 0.943d0 ) then
!           pi0 pi0 eta
          call cmkptc(kpion, 0, 0, a(1))
          call cmkptc(kpion, 0, 0, a(2))
          call cmkptc(keta, 0, 0, a(3))
          call cnbdcy(3, pj%mass, a, 0, w, icon)
          np = 3
       elseif( u .lt. 0.965d0) then
!     g g
          call cmkptc(kphoton, 0, 0, a(1))
          call cmkptc(kphoton, 0, 0, a(2))
          call c2bdcy(pj, a(1), a(2))
          np = 2
       else
!         omega gamma  
          call cmkptc(komega, 0, 0, a(1))
          call cmkptc(kphoton, 0,0, a(2))
          call c2bdcy(pj, a(1), a(2))
          np = 2
       endif
       if( np > 2 ) then
          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: