cmigdFunc.f File Reference

Go to the source code of this file.

Functions/Subroutines
real *8 function	cmigdg (s)
real *8 function	cmigdpsi (s)
real *8 function	cmigdgzai (s)
real *8 function	cmigdsforbrem (ee, rho, v)
real *8 function	cmigdsforpair (eg, rho, v)
subroutine	csetlpmcnst (s1in, logs1in, vmin, X0in)

Function/Subroutine Documentation

cmigdg()

real*8 function cmigdg

(

real*8

s

)

Definition at line 52 of file cmigdFunc.f.

References d0, and parameter().

      implicit none
      real*8 s  ! input 
!
!        G(s)/( 12pi*s^2/(1+12pi*s^2)) is approximated by polinomials
!
      real*8 xxx1(5)   ! pol coeff. for s < 0.12  G(s) error is 3 x 10^-5 
      real*8 xxx2(7)   ! pol coeff. for s > 0.11       error is 1 x 10^-3

      real*8 pi, const, f
      integer i
      parameter(pi= 3.141592654d0, const = 12.d0 * pi)

          data ( xxx1(i), i=  1,   5)/
     1     0.99995995    ,  -6.2614068    ,   57.629644    ,
     2   -242.27058    ,   450.83728                                    
     * /   
  
          data ( xxx2(i), i=  1,   7)/
     1         0.78283668    ,  -1.3317861    ,   12.312840    ,
     2   -33.848229    ,   44.733673    ,  -29.201775    ,
     3         7.5561750                                                     
     * /   


       if(s .lt. 0.12) then
          f  = 0.       
          do i = 5, 2, -1
             f = (f + xxx1(i))* s
          enddo
          f = f + xxx1(1)
          cmigdg = f * const*s**2/(1. + const*s**2)
       elseif(s .lt. 1.1) then
          f  = 0.       
          do i = 7, 2, -1
             f = (f + xxx2(i))* s
          enddo
          f = f + xxx2(1)
          cmigdg = f * const*s**2/(1. + const*s**2)
       else
          cmigdg = 1. - 0.022/s**4
       endif

Here is the call graph for this function:

cmigdgzai()

real*8 function cmigdgzai

(

real*8

s

)

Definition at line 146 of file cmigdFunc.f.

       implicit none
       real*8 s
       common /cmigdc/ s1, logs1, vm, logvm, smpa1, smpa2, suma,
     *   a2prob, x0
      real*8 vm, logvm, smpa1, smpa2, suma, a2prob, x0
      real*8 s1  ! /1.1185e-4/  ! (Z^0.333/183)^2 ; Z=7.25
      real*8 logs1   !  /-9.0983/  ! ln(s1)


!
       if( s .gt. 1.) then
          cmigdgzai = 1.
       elseif(s .gt. s1) then
          cmigdgzai = log(s)/logs1 + 1.
       else
          cmigdgzai =2.
       endif

cmigdpsi()

real*8 function cmigdpsi

(

real*8

s

)

Definition at line 103 of file cmigdFunc.f.

References d0.

       implicit none
       real*8 s  ! input

       real*8 f
       real*8 xxx1(4)   ! coef. of pol. s<0.3  better than 0.1 % for Psi
       real*8 xxx2(5)   ! //            s> 0.3  //
       integer i

          data ( xxx1(i), i=  1,   4)/
     1    1.00    ,   2.6978063    ,  -9.4242869    ,
     2    11.468298                                                     
     * /   

          data ( xxx2(i), i=  1,   5)/
     1    1.2095058    ,  0.57895055    ,  -1.6531094    ,
     2    1.4846143    , -0.46392960                                    
     * /   

       if(s .lt. 0.3) then
          f = 0.
          do i = 4, 2, -1
             f = (f + xxx1(i)) * s
          enddo
          f = f + xxx1(1)
          cmigdpsi = f * 6.d0 *s /(1.d0 + 6.d0*s)
       elseif(s .lt. 1.2) then
          f = 0.
          do i = 5, 2, -1
             f = (f + xxx2(i)) * s
          enddo
          f = f + xxx2(1)
          cmigdpsi = f * 6.d0 *s /(1.d0 + 6.d0*s)
       else
          cmigdpsi = 1. - 0.012/s**4
       endif

cmigdsforbrem()

real*8 function cmigdsforbrem	(	real*8	ee,
		real*8	rho,
		real*8	v
	)

Definition at line 170 of file cmigdFunc.f.

References d0, e, e3, and masele.

      implicit none
#include "Zmass.h"
!      #include "Zelemagp.h"
!
      real*8 ee     ! input.  Electron energy in GeV.
      real*8 rho    ! input.  air density in kg/m3
      real*8 v      ! input.  Fractional energy of gamma, Eg/Ee
!
      real*8 x, ss, ss2, cmigdgzai, temp
      integer i

      common /cmigdc/ s1, logs1, vm, logvm, smpa1, smpa2, suma,
     *   a2prob, x0
      real*8 vm, logvm, smpa1, smpa2, suma, a2prob, x0
      real*8 s1  ! /1.1185e-4/  ! (Z^0.333/183)^2 ; Z=7.25
      real*8 logs1   !  /-9.0983/  ! ln(s1)

      if(v .eq. 1.) then
         cmigdsforbrem = 1.
      else
         x = ee * rho   ! GeV kg/m3
         ss = 1.
         temp =  masele*  x0/x *100. * v/(1.-v)   !100: X0/rho in cm

         do i = 1, 10
!                normally 3 to 4 iteration is ok
            ss2 = 1.37e3 * sqrt(temp/cmigdgzai(ss))
            if(abs(ss/ss2 -1.d0) .lt. 3.e-3) goto 200
            ss = ss2
         enddo
         ss2 = 1. 
 200     continue
         cmigdsforbrem = ss2
!       
!        write(*, *) i, sngl(ss2), sngl(v), sngl(ee*rho)
!       
      endif

cmigdsforpair()

real*8 function cmigdsforpair	(	real*8	eg,
		real*8	rho,
		real*8	v
	)

Definition at line 214 of file cmigdFunc.f.

References d0, e, e3, and masele.

      implicit none
#include "Zmass.h"
!      #include "Zelemagp.h"
!
      real*8 eg     ! input.  Gamma energy in GeV.
      real*8 rho    ! input.  air density in kg/m3
      real*8 v      ! input.  Fractional energy of e-/e+ Ee/Eg
!
      real*8 x, ss, ss2, cmigdgzai, temp
      integer i
      common /cmigdc/ s1, logs1, vm, logvm, smpa1, smpa2, suma,
     *   a2prob, x0
      real*8 vm, logvm, smpa1, smpa2, suma, a2prob, x0
      real*8 s1  ! /1.1185e-4/  ! (Z^0.333/183)^2 ; Z=7.25
      real*8 logs1   !  /-9.0983/  ! ln(s1)      

      if(v .eq. 1.) then
         cmigdsforpair = 1.
      else
         x = eg * rho   ! GeV kg/m3
         ss = 1.
         temp =  masele*  x0/x *100./v/(1.-v)   !100: X0/rho in cm
         do i = 1, 10
!                normally 3 to 4 iteration is ok
            ss2 = 1.37e3 * sqrt(temp/cmigdgzai(ss))
            if(abs(ss/ss2 -1.d0) .lt. 3.e-3) goto 200
            ss = ss2
         enddo
         ss2 = 1. 
 200     continue
         cmigdsforpair = ss2
!       
!         write(*, *) i, sngl(ss2), sngl(v), sngl(eg*rho)
!       
      endif

csetlpmcnst()

subroutine csetlpmcnst	(	real*8	s1in,
		real*8	logs1in,
		real*8	vmin,
		real*8	X0in
	)

Definition at line 255 of file cmigdFunc.f.

Referenced by cbeginrun().

      implicit none
      real*8 s1in ! Migdals s1 = (Z^0.333/183)^2 
      real*8 logs1in ! log(s1)
      real*8 vmin  ! Egmin/Ee
      real*8 x0in  ! radiation length in kg/m^2. if it is
                   ! g/cm2, multiply it by 0.1/.00001 = 10.0
      common /cmigdc/ s1, logs1, vm, logvm, smpa1, smpa2, suma,
     *   a2prob, x0
      real*8 vm, logvm, smpa1, smpa2, suma, a2prob, x0
      real*8 s1  ! /1.1185e-4/  ! (Z^0.333/183)^2 ; Z=7.25
      real*8 logs1   !  /-9.0983/  ! ln(s1)

      s1 = s1in
      logs1 = logs1in
      vm = vmin
      logvm = log(vm)
      smpa1 = 0.5
      smpa2 = -4.0*logvm/3.0
      suma = smpa1 + smpa2
      a2prob = smpa2/suma
      x0 = x0in  

Here is the caller graph for this function: