cgetctau.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cgetctau (proj, ctau)

Function/Subroutine Documentation

cgetctau()

subroutine cgetctau	(	type(ptcl)	proj,
		real*8	ctau
	)

Definition at line 6 of file cgetctau.f.

References c, k0l, k0s, kbomega, kddb, kdmes, kds, keta, ketap, kgzai, kkaon, klambda, klambdac, kmuon, knnb, knuc, komega, kphi, kpion, krho, ksigma, ktau, parameter(), regptcl, t0d0, t0eta, t0gzai0, t0k, t0kl, t0ks, t0lambda, t0lambdac, t0mu, t0pi, t0pi0, t0sigma0, t0sigmam, t0sigmap, and t0tau.

Referenced by cdecayleng(), and cmydecay().

      implicit none

#include  "Zglobalc.h"
#include  "Zlife.h"
#include  "Zptcl.h"
#include  "Zcode.h"
!     **************************************************
      type(ptcl)::proj
      real*8 ctau
!
!
      real*8 tcpi, tck, tcmu, tcks, tckl, tcpi0, tcd0, tcdc, 
     *       tceta, tcsigma0, tcsigmap, tcsigmam, tcgzai0,
     *       tcgzaim, tcbomega, tclambda, tclambdac

      parameter(tcpi=t0pi*c, tck=t0k*c,
     *        tcmu=t0mu*c, tcks=t0ks*c,
     *        tckl=t0kl*c, tcsigma0=t0sigma0*c,
     *        tcsigmap = t0sigmap*c, tcsigmam=t0sigmam*c,
     *        tcgzai0 = t0gzai0 * c, tcgzaim=t0gzaim*c,
     *        tclambda =t0lambda*c, tclambdac=t0lambdac*c,
     *        tcbomega = t0bomega*c,
     *        tcpi0=t0pi0*c, tcd0=t0d0*c,
     *        tcdc =t0dc*c, tceta= t0eta*c)

      real(8),parameter:: tcds=t0ds*c, tctau=t0tau*c, tcxic=t0xic*c,
     *     tcomec0 = t0omec0*c, tcetap= t0etap*c, tcxic0=t0xic0*c
!     *     tcrho = t0rho*c,
!     *     tcomega = t0omega*c, tcphi=t0phi*c 
            
      
      if(proj%code .eq. kpion) then
          if(proj%charge .ne. 0) then
             ctau = tcpi
          else
             ctau = tcpi0
          endif
      elseif(proj%code .eq. kkaon) then
         if(proj%charge .ne. 0) then
            ctau = tck
         elseif(abs(proj%subcode) .eq. k0l) then
            ctau = tckl
         elseif(abs(proj%subcode) .eq. k0s) then
            ctau = tcks
         else
            write(0,*) ' K0 subcode=', proj%subcode, ' stragne'
            stop
         endif
      elseif(proj%code .eq. kmuon) then
         ctau = tcmu
      elseif(proj%code == knuc ) then
         ctau = infty
      elseif(proj%code .eq. knnb) then
         ctau =0.
      elseif(proj%code .eq. kddb) then
         ctau = 0.
      elseif(proj%code .eq. kdmes) then
         if(proj%charge .eq. 0) then
            ctau = tcd0
         else
            ctau = tcdc
         endif
      elseif(proj%code .eq.  ksigma) then
         if(proj%charge .eq. 0) then
            ctau = tcsigma0
         elseif(proj%charge .eq. 1) then
            if(proj%subcode .eq. regptcl) then
               ctau = tcsigmap
            else
               ctau = tcsigmam
            endif
         else
            if(proj%subcode .eq. regptcl) then
               ctau = tcsigmam
            else
               ctau = tcsigmap
            endif
         endif
      elseif(proj%code .eq. kgzai) then
         if(proj%charge .eq. 0) then
            ctau = tcgzai0
         else
            ctau = tcgzaim
         endif
      elseif(proj%code .eq. klambda ) then
         ctau = tclambda
      elseif(proj%code .eq. klambdac) then
         ctau = tclambdac
      elseif(proj%code .eq. kbomega) then
         ctau = tcbomega
      elseif(proj%code .eq. krho) then
         ctau = 0.
      elseif(proj%code .eq. kphi) then
         ctau =0.
      elseif(proj%code .eq. komega) then
         ctau = 0.
      elseif(proj%code .eq. keta ) then
         ctau = tceta
      elseif(proj%code .eq. kds  ) then
         ctau = tcds
      elseif( proj%code .eq. ketap  ) then
         ctau = tcetap
      elseif( proj%code .eq. kdelta  ) then
         ctau = 0
      elseif( proj%code .eq. kxic  ) then
         if( proj%charge == 0) then
            ctau = tcxic0
         else
            ctau = tcxic
         endif
      elseif( proj%code .eq. komec0  ) then
         ctau = tcomec0
      elseif( proj%code .eq. ktau  ) then
         ctau = tctau
      else
         ctau = infty
      endif

Here is the call graph for this function:

Here is the caller graph for this function: