cBrSSamp.f

Go to the documentation of this file.

      subroutine cbrssampp(Ee, prob)
!          brems prob. sampling by Seltzer data.
      implicit none
#include "ZbpCnst.h"
#include "ZbpTable.h"

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

      real*8 ale


      if(ee .le. breemins) then
         prob= 1.d-40
      else
         ale=log10(ee)
         call kintp3(brtxs,
     *   1, brtxts, brleemins,
     *   brdetxs, ale, prob) 
      endif
      end
!     ************
      subroutine cbrssampe( Eein, Eg)
!     ************
!         brems energy by Seltzer
      implicit none
#include "Zmass.h"
#include "ZbpCnst.h"
#include "ZbpTable.h"

      real*8 Eein,  Eg

      real*8 Ee
      real*8 u, ale, us, ans
      integer count
      logical ok
      
      count = 0
      if(eein-masele .le. bregmins ) then
!         simply avoid loop;  ( Eek < 10 keV)
         call rndc(u)
         eg =( eein-masele ) * u
         return  ! ************
      elseif(eein .lt. breemins) then
!         this and above  happne when electron loses energy
!         after sampling of processes 
         ee =  breemins
      else
         ee = eein
      endif
      ale = log10(ee)

 10   continue
      call rndc(u)
      if(u .gt. bruminsa) then
!          region A

         call k4ptdi(brstsa, 
     *        bruszsa, 
     *        bres,
     *        bruszsa, 
     *        bruminsa,
     *        brleemins,
     *        brdusa,
     *        brdes, u,  ale,  ans)  
         eg= exp( ans*(1.-u))*bregmins
      else
!         region B
         us = u**0.25d0

         call k4ptdi(brstsb, 
     *        bruszsb, 
     *        bres,
     *        bruszsb, 
     *        0.d0,
     *        brleemins,
     *        brdusb,
     *        brdes, us,  ale,  ans)  
         eg = exp(-ans*u)*(ee-masele)
!         Eg = exp(-ans*u)*Ee
      endif
      ok = eein-eg .ge. masele
      if(.not. ok) then
         count  = count + 1
         if(count .gt. 100 ) then
            call rndc(u)
            eg = (eein-masele)*u
            return  !  ******
         endif
         goto 10
      endif
      end