cmuPrsmp.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cmuprsmpp (Emu, prob, path)
subroutine	cmuprsmpe (Emu, Epair)
subroutine	cmuprsmp0 (a, b, vc, v)

Function/Subroutine Documentation

cmuprsmp0()

subroutine cmuprsmp0	(	real*8	a,
		real*8	b,
		real*8	vc,
		real*8	v
	)

Definition at line 50 of file cmuPrsmp.f.

References d0, false, and rndc().

Referenced by cmuprsmpe().

!     **********************************
      implicit none
!       sample v from  dv/v**(2-b)/(1+av**b)**2
!       This is a good approximation for pair spectrum
!       from muon 
      real*8 a     ! input.
      real*8 b     ! input.
      real*8 vc    ! input.  minimum of v
      real*8 v     ! output. sampled v


      real*8 u, temp


      logical ok

      ok = .false.
      do while (.not. ok)
         if(b .ne. 1.0d0) then
            temp = vc**(b-1)
            call rndc(u)
            v = ((1-temp)*u + temp)**(1./(b-1))
            call rndc(u)
            ok =(u .lt. ((1.+a*temp*vc)/(1.+a*v**b))**2) 
         else
            call rndc(u)
            v = vc**u
            call rndc(u)
            ok =(u .lt. ((1.+a*vc)/(1.+a*v))**2) 
         endif
      enddo

Here is the call graph for this function:

Here is the caller graph for this function:

cmuprsmpe()

subroutine cmuprsmpe	(	real*8	Emu,
		real*8	Epair
	)

Definition at line 32 of file cmuPrsmp.f.

References cmuprsmp0(), mupa, mupb, muqb, and mura.

Referenced by cmuinte().

      implicit none
#include "Zcmuint.h"
      real*8 emu           ! input. muon total energy in GeV
      real*8 epair         ! output. sampled energy loss of muon
                           !        (Epair/Emu > media.cnst.muPrVmin)

      real*8  a, b, v

      a = mupa*emu**(-0.40) + mura
      b = mupb*emu**muqb  +  0.92

      call cmuprsmp0(a, b, muprvmin, v)
      epair = emu*v

Here is the call graph for this function:

Here is the caller graph for this function:

cmuprsmpp()

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

Definition at line 2 of file cmuPrsmp.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 pair creation prob. /X0
      real*8 path            ! output. sampled path in r%l
                             ! if Emu < media.cnst.muPrEmin,  prob=0 
                             !  and path becomes big

      real*8 u, ale

      if(emu .le. mupremin) then
         prob = 0.
      elseif(emu .le. mupremax1) then
         ale = log10(emu)
         call kintp3(muprtx, 1, muprtxt, muprlemin,
     *   muprdetx, ale, prob)
      else
!          prob is const for  v > vmin; Emu > muPrEmax1         
         prob = muprtx(muprtxt)
      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: