cmuNsmp.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cmunsmpp (Emu, prob, path)
subroutine	cmunsmpe (Emu, Et)

Function/Subroutine Documentation

cmunsmpe()

subroutine cmunsmpe	(	real*8	Emu,
		real*8	Et
	)

Definition at line 33 of file cmuNsmp.f.

References d0, kpolintp2(), and rndc().

Referenced by cmuinte().

      implicit none
#include "Zcmuint.h"
      real*8 emu           ! input. muon total energy in GeV
      real*8 et            ! output. sampled energy transfer
!

      real*8  ale, u, uu,   v


      real*8 a
      real*8 error

      data a/0.02d0/   ! to change this, you must change creation part



      ale = log10(min(emu,  munemax))


      call rndc(u)
      uu = a*u/(1+a-u)  ! uniform in this variable
         call kpolintp2(0.d0, 1, mundu,
     *        munlemin, 1, munde,
     *        muntbl,  munusize,  munusize, munesize,
     *        5, 3,   uu, ale, v, error)
      et = 10.d0**v * emu

Here is the call graph for this function:

Here is the caller graph for this function:

cmunsmpp()

subroutine cmunsmpp	(	real*8	Emu,
		real*8	prob,
		real*8	path
	)

Definition at line 2 of file cmuNsmp.f.

References kintp3(), and rndc().

Referenced by csampmuint().

      implicit none
#include "Zcmuint.h"

      real*8 emu             ! input.  muon total energy in GeV
      real*8 prob            ! output. muon nuc. int prob. /X0
      real*8 path            ! output. sampled path in r%l
                             ! if Emu < media.cnst.muNEmin,  prob=0 
                             ! and path becomes big

      real*8 u, ale

      if(emu .le. munemin) then
         prob = 0.
      elseif(emu .le. munemax1) then
         ale = log10(emu)
         call kintp3(muntx, 1, muntxt, munlemin,
     *   mundetx, ale, prob)
      else
!        small power dependence
         prob = muntx(muntxt)*
     *      (emu/munemax1)**munpwtx
      endif
      if(prob .gt. 0.) then
         call rndc(u)
         path =- log(u)/prob
      else
         path = 1.d30
      endif

Here is the call graph for this function:

Here is the caller graph for this function: