cBrLSamp.f

Go to the documentation of this file.

      subroutine cbrlsampp(Ee, prob)
      implicit none
#include "ZbpCnst.h"
#include "ZbpTable.h"

      real*8 Ee
      real*8 prob  ! output probability of Brems / X0

      real*8 ale
!        Ee
!           BremEemin  BremEeCS  BremEemaxL  BremEemaxH
!     ----------------------------------------------------
!                |         |            |          |
!     no brems       p.s        c.s        lpm       lpm by
!                                                    rejection
!

      if(ee .le. bremeemin) then
         prob= 1.d-40
      else
         ale=log10(ee)
         call kintp3(brtxl,
     *   1, bremtxtl, bremleemin,
     *   bremdetxl, ale, prob) 
      endif
      end
!     ************
      subroutine cbrlsampe(Ee, Eg)
!     ************
      implicit none
#include "ZbpCnst.h"
#include "ZbpTable.h"

#include "Zmass.h"

      real*8 Ee,  Eg

      real*8 u, ale, us, ans

      call rndc(u)
      ale = log10(ee)
      if(u .gt. bremuminla) then
!          region A

         call k4ptdi(brstla, 
     *        bremuszla, 
     *        bremesize,
     *        bremuszla, 
     *        bremuminla,
     *        bremleemin,
     *        bremdula,
     *        bremdel, u,  ale,  ans)  
         eg= exp( ans*(1.-u))*bremegmin
      else
!         region B
         us = u**0.25d0

         call k4ptdi(brstlb, 
     *        bremuszlb, 
     *        bremesize,
     *        bremuszlb, 
     *        0.d0,
     *        bremleemin,
     *        bremdulb,
     *        bremdel, us,  ale,  ans)  
         eg = exp(-ans*u)*(ee - masele)
      endif
      end