cmydecay.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cmydecay (Jdecay, tau, a, nin, nout)

Function/Subroutine Documentation

cmydecay()

subroutine cmydecay	(	integer, dimension(klast), intent(in)	Jdecay,
		real(8), intent(in)	tau,
		type (ptcl), dimension(*)	a,
		integer, intent(in)	nin,
		integer, intent(out)	nout
	)

Definition at line 4 of file cmydecay.f.

References a, b, cbomegadcy(), cddecay(), cerrormsg(), cetadecay(), cgetctau(), cgzaidecay(), ckaondecay(), clambdacdcy(), clambdadcy(), cpi0decay(), csigmadecay(), false, k0s, kbomega, kdmes, keta, kgzai, kkaon, klambda, klambdac, kpion, ksigma, and rndc().

Referenced by __gencol2.f__(), __gencol3.f__(), and __gencol4.f__().

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

      integer,intent(in)::jdecay(klast)  ! if Jdecay(i) !=  0
         !  the ptcl is decayed, if i is one of
         !  kbomega, kdmes, keta, kgzai,  klambda, klambdac, kpion, ksigma
         !  (kpion --> pi0)
      real(8),intent(in)::tau ! life time (sec), above which the partcles
                 ! are  regarded as stable and not made to decay.
                 ! Below this, decay is forced. For the decay product,
                 ! the same is applied.
            ! IF tau =0, random sampling is tried if the
            !  particle decay or not using correct life time and 
            !  distance dist. (see below)
            ! if decay   takes place, decay product is put 
      integer,intent(in):: nin  ! # of ptcls in a
      type(ptcl):: a(*)  ! input/output.  size of a must be large
                            ! say, > (nout - nin) *3 + nin
      integer,intent(out):: nout   !  ptcls in a. nout >= nin.


      type(ptcl):: b(maxn)  ! working
      integer i, j, code, m, n
      real(8):: pol, ctau
      real(8),save::dist=141.2 ! m from col. to LHCf
      real(8) u, gamma, decayl
      logical decay
!//////
#ifdef  SEEDETAIL
      integer nchgin, nchgout
      nchgin = 0
      do i = 1, nin
         if(a(i)%charge /= 0) nchgin= nchgin+1
      enddo
#endif
!//////////////
      n = 0
      do i = 1, nin
         code = a(i)%code
         j = 0
         if( jdecay(code) /=  0 .or. tau == 0. ) then
            call cgetctau(a(i), ctau)
            if(tau == 0. ) then
!                random sample for decay
               gamma = a(i)%fm%p(4)/a(i)%mass
               call rndc(u) 
               decayl = -log(u)*ctau*gamma
               decay = (decayl < dist )
            else
               decay = (ctau/3.e8 < tau)
            endif
            if(.not. decay) then
               n = n + 1
               b(n) = a(i)
            else
               if( code == kpion ) then
                  if( a(i)%charge == 0) then
                     call cpi0decay( a(i),  b(n+1), j)
                  endif
               elseif( code == kkaon ) then
                  if( a(i)%subcode == k0s ) then
                     call ckaondecay(a(i), .false.,  b(n+1), j, pol)
                  endif
               elseif( code .eq. kdmes ) then
                  call cddecay( a(i), b(n+1), j)
               elseif( code .eq. keta ) then
                  call cetadecay( a(i), b(n+1), j)
               elseif( code .eq. kgzai ) then
                  call cgzaidecay( a(i), b(n+1), j ) 
               elseif( code .eq. klambda ) then
                  call clambdadcy( a(i), b(n+1), j ) 
               elseif( code .eq. klambdac ) then
                  call clambdacdcy( a(i), b(n+1), j ) 
               elseif( code .eq. ksigma ) then
                  call csigmadecay( a(i), b(n+1), j ) 
               elseif( code .eq. kbomega ) then
                  call cbomegadcy( a(i), b(n+1), j ) 
               else
                  write(*,*)  ' code =', code
                  call cerrormsg('cmydecay error', 0)
               endif
               if(j > 0 ) then
!/////////////////////
#ifdef SEEDETAIL
                  write(0,*) ' j=',j, ' code=',code, ' chg=',a(i)%charge
#endif
!///////////////
                  n = n + j
               endif
            endif
         else
            n = n + 1
            b(n) = a(i)
         endif
      enddo
!/////////
#ifdef SEEDETAIL
      nchgout = 0
      do i = 1, n
         if(b(i)%charge /= 0) nchgout = nchgout+1
      enddo
      write(0,*) ' chgin out', nchgin, nchgout
#endif
!/////////////////
      do i = 1, n
         a(i) = b(i)
      enddo
      nout = n

Here is the call graph for this function:

Here is the caller graph for this function: