cbremsPath.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cbremspath (Ee, path)
subroutine	cbremsenergy (Ee, Eg)

Function/Subroutine Documentation

cbremsenergy()

subroutine cbremsenergy	(	real*8	Ee,
		real*8	Eg
	)

Definition at line 42 of file cbremsPath.f.

References a, cbrcsampe(), cbrlsampe(), cbrssampe(), d, e, and masele.

Referenced by cbrems().

      implicit none
#include "Zmass.h"
#include "ZbpCnst.h"


      real*8 ee            ! input. e-/e+ energy in GeV
      real*8 eg            ! output. sampled Eg in GeV


 10   continue
      if(ee .le. breemaxs) then
!            Seltzer region
         call cbrssampe(ee,  eg)

      elseif(ee .le. compscre) then
!            screeinig region
         call cbrlsampe(ee,  eg)

      elseif(ee .le. bremeemaxl) then
!            complete screeing region
         call cbrcsampe(ee, eg)

      else
!            LPM region; same as path part. should not come
!         call cBrHSampE(Ee, Eg)
        call cbrcsampe(ee, eg)
      endif
      if( ee - masele <= eg ) goto 10

Here is the call graph for this function:

Here is the caller graph for this function:

cbremspath()

subroutine cbremspath	(	real*8	Ee,
		real*8	path
	)

Definition at line 2 of file cbremsPath.f.

References cbrcsampp(), cbrlsampp(), cbrssampp(), and rndc().

Referenced by csampeintl().

!          generic brems sampling routine
      implicit none
#include "ZbpCnst.h"
!      entry cBrSampP(Ee, prob, path)


      real*8 ee            ! input. Electron/positron energy in GeV
      real*8 prob          ! output. Brems prob. per r%l
      real*8 path          ! ouptut. sampled path in r%l

      real*8 u

      if(ee .le. breemaxs) then
!            Seltzer region
         call cbrssampp(ee,  prob)
      elseif(ee .le. compscre) then
!            screeinig region
         call cbrlsampp(ee,  prob)
      elseif(ee .le. bremeemaxl) then
!            complete screeing region
         call cbrcsampp(ee, prob)
      else
!         LPM region;
!         call cBrHSampP(Ee, prob); this cannot be used due to rho change.
!         should not come here; LPM routine is called
!         from outer place. However, control may come here
!         because of small BremEemaxL ~ 94 GeV. 
!         If you are very serious about LPM, the effect may
!         be seen in 100 GeV region at very small gamma energies
!         However, the effect in air shower development could be
!         seen  at  least > 10^15 eV. Therefore  we use complete
!         screening cross section.
!         
         call cbrcsampp(ee, prob)
      endif
      call rndc(u)
      path = - log(u)/prob

Here is the call graph for this function:

Here is the caller graph for this function: