clambdacDcy.f File Reference

Go to the source code of this file.

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

Function/Subroutine Documentation

clambdacdcy()

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

Definition at line 8 of file clambdacDcy.f.

References antip, cibst1(), cmkptc(), cnbdcy(), klambda, kmuon, kneumu, kpion, regptcl, and rndc().

Referenced by cintelambdac(), 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
       real*8 u, w
       integer icon
!           lambdac --> lamuda +  mu+ +  neumu  (2 %)
!
       integer:: i

        call rndc(u)
        if(u .lt. .02) then
           if(pj%subcode .eq. regptcl) then
              call cmkptc(klambda, regptcl, 0, a(1))
              call cmkptc(kmuon,   antip,  1,  a(2))
              call cmkptc(kneumu,  regptcl, 0, a(3))
           else
              call cmkptc(klambda, regptcl, 0, a(1))
              call cmkptc(kmuon,   regptcl, -1, a(2))
              call cmkptc(kneumu,  antip, 0, a(3))
           endif
           call cnbdcy(3, pj%mass, a, 0, w, icon)
           np=3
        else
!         many modes with small branching ratios
!         and  no major decay mode but
!           p+anything or n+anything is about 50 %,50%
!        so we use Lambda0 + anything for simplicity
!        and for anything we take pi+,pi0
           call cmkptc(klambda, regptcl, 0, a(1))
           call cmkptc(kpion,   regptcl, 1, a(2))
           call cmkptc(kpion,   regptcl, 0, a(3))
           call cnbdcy(3, pj%mass, a, 0, w, icon)
           np=3
        endif
        do   i=1, np
           call cibst1(i, pj, a(i), a(i))
        enddo

Here is the call graph for this function:

Here is the caller graph for this function: