cdedxElec.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cdedxelec (eini, ic, dedt)

Function/Subroutine Documentation

cdedxelec()

subroutine cdedxelec	(	real*8	eini,
		integer	ic,
		real*8	dedt
	)

Definition at line 36 of file cdedxElec.f.

References stha, sthb, sthc, sthsa, sthx0, sthx1, w0, and wlg0.

       implicit none
!
!
#include  "ZdedxAir.h"
!
      real*8 eini, dedt
      integer ic
!      external cblkdedx

      real*8 emass, emassg
      real*8 ek, ein, e, gi, beta2, x, a, b, c, cb, x0, x1
      real*8 dltx, wm, wlg
!
      real*8 peake/2.3d-7/, peak/0.858d0/
      data emass/0.511d0/, emassg/0.511d-3/


!
       if(jdef .eq. 0) then
          write(*, *) ' cdedxEleci must be called'
          stop 9999
       endif

       ek=eini-emassg
       ein=eini

       if( ek .lt. peake) then
              dedt=peak*sqrt(ek/peake)
       else
!                  energy in mev unit
              e=ein*1000.d0
              gi=emass/e
              beta2= 1.d0 - gi**2
!                  x=log10(p/mc)
              x=log10(e**2 - 0.2611d0)/2.0d0 + 0.29d0
              a=stha
              b=sthb
              c=sthc
              cb=-c
              x0=sthx0
              x1=sthx1
              if(x .lt. x0) then
!                    4.605x - dlt
                  dltx=4.605d0*x
              elseif(x .lt. x1) then
                  dltx=cb - (x1-x)**3 * sthsa
              else
                  dltx=cb
              endif
              wm=e- emass
!              if(Knckon) then
!                Now we don't use Knckon, since big wm can do
!               the equivalent
                   if(wm .gt. w0) then
                       wlg=wlg0
                   else
                       wlg=log(wm)
                   endif
!              else
!                   wlg=log(wm)
!              endif
              if(ic .eq. -1) then
                  dedt=a/beta2 * (b + 1.12d0 + wlg-beta2 +dltx)
              else
                  dedt=a/beta2 *( b+0.693d0 + wlg -2*beta2 + dltx)
              endif
!                            convert it to gev/(g/cm^2)
              dedt=dedt *1.d-3
          endif