csigmaDecay.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	csigmadecay (pj, a, np)
subroutine	csigmapdcy (pj, a, np)
subroutine	csigma0dcy (pj, a, np)
subroutine	csigmamdcy (pj, a, np)

Function/Subroutine Documentation

csigma0dcy()

subroutine csigma0dcy	(	type(ptcl)	pj,
		type(ptcl), dimension(*)	a,
		integer	np
	)

Definition at line 60 of file csigmaDecay.f.

References c2bdcy(), cmkptc(), klambda, and kphoton.

Referenced by csigmadecay().

!
!         sigma0 decay
!         1) Sigma---->Lamda gamma  100 %
      implicit none
#include  "Zptcl.h"
#include  "Zcode.h"
      integer np
      type(ptcl):: pj
      type(ptcl):: a(*)
      integer subcode

      subcode = pj%subcode
      call cmkptc(klambda, subcode, 0, a(1) )
      call cmkptc(kphoton, 0, 0, a(2) )
      call c2bdcy(pj, a(1), a(2))
      np=2

Here is the call graph for this function:

Here is the caller graph for this function:

csigmadecay()

subroutine csigmadecay	(	type(ptcl)	pj,
		type(ptcl), dimension(*)	a,
		integer	np
	)

Definition at line 8 of file csigmaDecay.f.

References csigma0dcy(), csigmamdcy(), and csigmapdcy().

Referenced by cintesigma(), and cmydecay().

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

       integer np               !output. no. of ptcls produced
       type(ptcl):: pj         ! input. kaon
       type(ptcl):: a(*)      ! output. produced ptcls
!
!
      if(pj%charge .eq. 1) then
!           Sigma +
         call csigmapdcy(pj, a, np)
      elseif(pj%charge .eq. 0) then
!           sigma 0
         call csigma0dcy(pj, a, np)
      elseif(pj%charge .eq. -1) then
!           sigma -1         
         call csigmamdcy(pj, a, np)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

csigmamdcy()

subroutine csigmamdcy	(	type(ptcl)	pj,
		type(ptcl), dimension(*)	a,
		integer	np
	)

Definition at line 78 of file csigmaDecay.f.

References c2bdcy(), cmkptc(), knuc, and kpion.

Referenced by csigmadecay().

!
!         sigma- decay
!         1) Sigma---->n pi-   100 %
      implicit none
#include  "Zptcl.h"
#include  "Zcode.h"      
      integer np
      type(ptcl):: pj
      type(ptcl):: a(*)
      integer subcode, charge

      subcode = pj%subcode
      charge = pj%subcode
      
!          n+pi-
         call cmkptc(kpion, 0, charge, a(1))
         call cmkptc(knuc, subcode, 0 , a(2))
         call c2bdcy(pj, a(1), a(2))
         np=2

Here is the call graph for this function:

Here is the caller graph for this function:

csigmapdcy()

subroutine csigmapdcy	(	type(ptcl)	pj,
		type(ptcl), dimension(*)	a,
		integer	np
	)

Definition at line 29 of file csigmaDecay.f.

References c2bdcy(), cmkptc(), knuc, kpion, regptcl, and rndc().

Referenced by csigmadecay().

!
!         sigma+ decay
!         1) Sigma---->p pi0      (51.6%0
!         2)      ---->n pi+       48.4 
      implicit none
#include  "Zptcl.h"
#include  "Zcode.h"
      integer np
      type(ptcl):: pj
      type(ptcl):: a(*)
      integer subcode, charge

      real*8 u
      call rndc(u)
      subcode = pj%subcode
      charge = -pj%subcode
      if(u  .lt. .516) then
!          p+pi0  or conjugae       
         call cmkptc(kpion, 0, 0, a(1))
         call cmkptc(knuc, subcode, charge, a(2))
         call c2bdcy(pj, a(1), a(2))
         np=2
      else
!          n pi+
         call cmkptc(kpion, regptcl, -charge, a(1))
         call cmkptc(knuc,  subcode,  0, a(2))
         call c2bdcy(pj, a(1), a(2))
         np=2
      endif

Here is the call graph for this function:

Here is the caller graph for this function: