cPrLSamp.f

Go to the documentation of this file.

      subroutine cprlsampp(Eg, prob)
      implicit none
#include "ZbpCnst.h"
#include "ZbpTable.h"
      real*8 Eg
      real*8 prob  ! output probability of Pair / X0


      real*8 ale


      if(eg .lt. pairegmin) then
         prob= 1.d-40
      else
         ale=log10(eg)
         call kintp3(prtxl,
     *   1, pairtxtl, pairlegmin,
     *   pairdetxl, ale, prob) 
      endif
      end
!     ************
      subroutine cprlsampe(Eg, Ee)
!     ************
!          samples higher energy pair electron
      implicit none
#include "ZbpCnst.h"
#include "ZbpTable.h"

      real*8 Ee,  Eg

      real*8 u, ale, us, ans, ex

      call rndc(u)
      ale = log10(eg)
      if(u .gt. pairuminla) then
!          region A
         call k4ptdi(prstla, 
     *        pairuszla, 
     *        pairesize,
     *        pairuszla, 
     *        pairuminla,
     *        pairlegmin,
     *        pairdula,
     *        pairdela, u,  ale,  ans)  
         ee = (ans*(1.-u) + 0.5d0)*eg
      else
!         region B
         us = u**0.25d0
         ex = sqrt(ale - pairlegmin)
         
         call k4ptdi(prstlb, 
     *        pairuszlb, 
     *        pairesize,
     *        pairuszlb, 
     *        0.d0,
     *        0.d0,
     *        pairdulb,
     *        pairdelb, us,  ex,  ans)  
         ee = ans* eg
      endif
      end