cdedx2atLE.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cdedx2atle (eini, amassi, dedt)
subroutine	cdedx2atlea (eini, amassi, dedt)

Function/Subroutine Documentation

cdedx2atle()

subroutine cdedx2atle	(	real*8	eini,
		real*8	amassi,
		real*8	dedt
	)

Definition at line 36 of file cdedx2atLE.f.

References cdedx2atlea(), d, and d0.

       implicit none
!
!
#include  "ZdedxAir.h"
      real*8 eini, amassi, dedt
      real*8  ek,  peak2, peak2e
!
      ek=eini-amassi
      if(ek .lt. 1.d-3) then    ! < 1MeV
         peak2e = 1.4d-6* (amassi/100.d-3)**1.6d0
         if(ek .lt. peak2e ) then
            if(ek .gt. 0.) then
               peak2 = 0.6d0*(ek/1.5d-5)**(-0.15d0)
               dedt=peak2*sqrt(ek/peak2e)
            else
               dedt = 0.
            endif
         else
            call cdedx2atlea(eini, amassi, dedt)
         endif
      else
         call cdedx2atlea(eini, amassi, dedt)
      endif
      if(dedt .lt. 0.) then
         dedt = 1.d-7
      endif

Here is the call graph for this function:

cdedx2atlea()

subroutine cdedx2atlea	(	real*8	eini,
		real*8	amassi,
		real*8	dedt
	)

Definition at line 65 of file cdedx2atLE.f.

References d, d0, parameter(), stha, sthb, sthc, sthsa, sthx0, sthx1, w0, and wlg0.

Referenced by cdedx2atle().

       implicit none
!
!
#include  "ZdedxAir.h"

!
      real*8 eini, amassi, dedt

      real*8 emass, emass2
      parameter(emass=.511d0, emass2=emass**2)
      real*8  e, amass, ein, gi, beta2, x, a, b, c, cb, x0
      real*8 dltx, p2, wm, wlg,  x1


       if(jdef .eq. 0) then
          write(*, *) ' cdedxEleci must be called'
          stop 9999
       endif
       ein=eini
!
!                  energy in mev unit
          e=ein*1000.d0
          amass=amassi*1000.d0
          gi=amass/e
          if(gi .ge. 1.d0) then
              dedt=.01d0
          else
            beta2= 1.d0 - gi**2
!                x=log10(p/mc)
            x=log10((e/amass)**2 - 1.d0) / 2.d0
            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
            p2=e**2 -amass**2
            wm=2*emass*p2/( amass**2+ emass2+ emass*e*2)
            if(wm .gt. w0) then
               wlg = wlg0
            else
               wlg=log(wm)
            endif
            dedt=a/beta2 *( b+0.693d0+wlg -2.d0*beta2 + dltx)
        endif
        dedt=dedt *1.d-3

Here is the call graph for this function:

Here is the caller graph for this function: