cpiMuDecay.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cpimudecay (pj, mupol, a, np, polari)

Function/Subroutine Documentation

cpimudecay()

subroutine cpimudecay	(	type(ptcl)	pj,
		logical	mupol,
		type(ptcl), dimension(*)	a,
		integer	np,
		real*8	polari
	)

Definition at line 10 of file cpiMuDecay.f.

References c2bdcy(), cmkptc(), cpimupolari(), kmuon, and kneumu.

Referenced by cintepion().

      implicit none
!----      include '../../Zptcl.h'
#include  "Zptcl.h"
#include  "Zcode.h"
      
      type(ptcl):: pj    !  input. pion
      logical mupol    !  if T, polarization is taken into account
      type(ptcl):: a(*)   ! to store decay product
      integer np     !    no. of ptcls stored in a 
      real*8   polari
      integer subcode, charge
!
!           pi --> mu+neumu (muon should be set last)
      subcode =  -pj%charge
      call cmkptc(kneumu,  subcode, 0,  a(1))
      charge =  pj%charge
      call cmkptc(kmuon,  0,  charge,   a(2))
      if(pj%fm%p(4) <= pj%mass) then
!                stopping pion.
!          In  some case, Ppi /=0 seems to result in.
!          reset Ppi
         pj%fm%p(1:3) = 0.
         pj%fm%p(4) = pj%mass
      endif
      call c2bdcy(pj, a(1), a(2))
!               set polarization of muon
      if(mupol) then
         call cpimupolari(pj, a(2), polari)
      else
         polari=0.
      endif
      np=2

Here is the call graph for this function:

Here is the caller graph for this function: