ccapnu.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	ccapnu (atms, z, a, np)

Function/Subroutine Documentation

ccapnu()

subroutine ccapnu	(	integer	atms,
		integer	z,
		type(ptcl)	a,
		integer	np
	)

Definition at line 2 of file ccapnu.f.

References cmkptc(), d0, kcossn(), kneumu, ksgmrm(), masmu, parameter(), regptcl, and rndc().

Referenced by cintemuon().

       implicit none

#include  "Zptcl.h"
#include  "Zmass.h"
#include  "Zcode.h"
       integer  atms, z
       integer np       ! output.
       type(ptcl):: a  ! output. 
!
!         atms: input.  atomic mass
!            z: input.   charge of the capturer
!         give neu(mu) which comes from captured
!         negative muon.
!        er:   recoil nutron average kinetic energy.
       real*8 er, xm, e, cs, sn, u, cost, x, sint

         parameter(er=15.e-3, xm=er/masmu/(masmu-er) )
!
         if(atms .eq. 1) then
             e=masmu
         else
!                sample x from x*exp(-x)dx with mean xm.
             call ksgmrm(1.00d0, xm, x)
!               energy of neutrino; peaked at around 100 MeV
             e= 1./(1./masmu + x )
         endif
!          make neutrino of muon type
         call cmkptc(kneumu, regptcl, 0, a)
         call kcossn(cs, sn)
         call rndc(u)
         cost=2*u-1.
         sint = sqrt(1.d0 - cost**2)
         a%fm%p(1) = cs*sint*e
         a%fm%p(2) = sn*sint*e
         a%fm%p(3) = cost*e
         a%fm%p(4) = e
         np=1

Here is the call graph for this function:

Here is the caller graph for this function: