cDeltaDecay.f File Reference

Go to the source code of this file.

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

Function/Subroutine Documentation

cdeltadecay()

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

Definition at line 8 of file cDeltaDecay.f.

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

Referenced by cintedelta().

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

       type(ptcl),intent(in):: pj       !  Delta
       type(ptcl),intent(inout):: a(*)    !  produced ptcls
       integer,intent(out):: np   ! no. of ptcls produced =2

       real(8):: u
       integer subcode, chargepi, chargen
!     Delta+  --> p pi0    or n pi+
!     Delta-  --> n pi-    
!     Delta0  --> n pi0   or  p pi-

!        Ng mode neglected       
       call rndc(u)

       if( pj%charge > 0 ) then
!              Delta+  --> p pi0    or n pi+
          if( u < 0.5 ) then
                !   p pi0 
             chargen = 1
             chargepi = 0
          else
             !   n pi+
             chargen = 0
             chargepi = 1
          endif
          call cmkptc(kpion, regptcl, chargepi, a(1))
          call cmkptc(knuc,  regptcl, chargen, a(2))
       elseif( pj%charge == 0 ) then
          !     Delta0  --> n pi0   or  p pi-
          if( u < 0.5 ) then
                !   n pi0 
             chargen = 0
             chargepi = 0
             subcode = regptcl
          else
             !   p pi-
             chargen = 1
             chargepi = -1
             subcode = antip
          endif
          call cmkptc(kpion, subcode, chargepi, a(1))
          call cmkptc(knuc,  regptcl, chargen, a(2))
       else
!     Delta-  --> n pi-              
          call cmkptc(kpion, antip, -1, a(1))
          call cmkptc(knuc,  regptcl, 0, a(2))
       endif
       

       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: