asinfo.f File Reference

Go to the source code of this file.

Functions/Subroutines
real function	sizex (E0in, NN, code, s)
real function	dep2cogdep (E0in, NN, dep)
real function	cogdep2nmu (E0, cogdep)
real function	cogdep2nh (E0, cogdep)
real function	s2nh (E0, s)
real function	cogdep2ne (E0, cogdep)
real function	cogdep2ng (E0, cogdep)
real function	cogdep2s (E0, NN, cogdep)
real function	dep2s (E0, NN, dep)
real function	s2ng (E0, s)
real function	s2ne (E0, s)
real function	s2nmu (E0, s)
real function	s2cogdep (E0, s)

Function/Subroutine Documentation

cogdep2ne()

real function cogdep2ne	(	real	E0,
		real	cogdep
	)

Definition at line 130 of file asinfo.f.

References cogdep2ne(), e, and e10.

!     ************************
      implicit none
      real e0 ! input GeV
      real cogdep ! input dep/cog
      real a, b, c, d
      real eps
      eps = 3.e-3

      a=1.e11
      if(abs(1.- 1e8/e0) .lt. eps ) then
         b=8.1044
         c=7.4068
         d=1.12022
      elseif( abs(1.- 1e9/e0) .lt.eps)  then
         b= 7.278 
         c= 5.063 
         d= 1.555
      elseif( abs(1.-1.e10/e0) .lt. eps) then
         b=6.3661
         c=2.7725
         d=2.2605
      else
         b =  -0.3774*log(e0) +  15.07
         c = -1.0063*log(e0) +  25.93
         d = 0.24760*log(e0) -3.48601
      endif   
      cogdep2ne=a*cogdep**b * exp(-c*cogdep**d)

Here is the call graph for this function:

cogdep2ng()

real function cogdep2ng	(	real	E0,
		real	cogdep
	)

Definition at line 160 of file asinfo.f.

References cogdep2ng(), and e10.

!     ************************
      implicit none
      real e0 ! input GeV
      real cogdep ! input dep/cog
      real a, b, c, d

      a=1.e12
      if(e0 .eq. 1e8) then
         b=8.00  
         c=8.034 
         d=0.941 
      elseif(e0 .eq. 1e9) then
         b=7.71  
         c=5.652 
         d=1.383
      elseif(e0 .eq. 1e10) then
         b= 6.91 
         c= 3.271 
         d= 2.022
      else
         b=-0.23669*log(e0)+12.445
         c=-1.03427*log(e0)+ 27.085
         d=0.2347*log(e0) -3.41583
      endif   
      cogdep2ng=a*cogdep**b * exp(-c*cogdep**d)

Here is the call graph for this function:

cogdep2nh()

real function cogdep2nh	(	real	E0,
		real	cogdep
	)

Definition at line 109 of file asinfo.f.

References cogdep2nh().

!     ***********************
      implicit none
      real e0 ! input GeV
      real cogdep ! input dep/cog
!
      real cogdep2nmu
!        use mu at present
      cogdep2nh= cogdep2nmu(e0,cogdep)

Here is the call graph for this function:

cogdep2nmu()

real function cogdep2nmu	(	real	E0,
		real	cogdep
	)

Definition at line 82 of file asinfo.f.

References cogdep2nmu(), and e10.

!     ************************
      implicit none
      real e0 ! input GeV
      real cogdep ! input dep/cog
      real a, b, c, d
      a=3.2e8 
      if(e0 .eq. 1e8) then
        b=4.8316
        c=5.9692 
        d=0.709
      elseif(e0 .eq. 1e9) then
         b=4.6796
         c=3.7835
         d=1.091
      elseif(e0 .eq. 1e10) then
         b=4.1251
         c=1.8791
         d=1.818
      else
         b = -0.1534*log(e0) + 7.7246
         c = -0.8881*log(e0) + 22.282
         d = 0.24081*log(e0) -3.7845
      endif   
      cogdep2nmu=a*cogdep**b * exp(-c*cogdep**d)

Here is the call graph for this function:

cogdep2s()

real function cogdep2s	(	real	E0,
		integer	NN,
		real	cogdep
	)

Definition at line 188 of file asinfo.f.

References cogdep2s().

!     ***************
!      for all energy      
      implicit none
      real e0   !  input not used now
      integer nn ! input incident nucleon number. not used now
      real cogdep ! input dep/cog
      cogdep2s = 10.28*cogdep**0.0407-9.228

Here is the call graph for this function:

dep2cogdep()

real function dep2cogdep	(	real	E0in,
		integer	NN,
		real	dep
	)

Definition at line 56 of file asinfo.f.

References dep2cogdep(), e, and e10.

!     ************************
      implicit none
      real e0in  ! input primary nucleus total energy in GeV.
      integer nn ! input incident nucleon number. 
      real dep ! input slant depth in g/cm^2

      real e0, eps
      eps = 3.e-3
      if(nn .eq. 0) then
         e0=e0in
      else
         e0=e0in/nn
      endif
!       function value is dep/cog. 
      if(abs(1.0- e0/1.e8) .lt. eps) then
         dep2cogdep=0.001442*dep
      elseif( abs(1.0-e0/1.e9) .lt. eps) then
         dep2cogdep=0.001247*dep
      elseif( abs(1.0-e0/1.e10) .lt. eps ) then
         dep2cogdep=0.001094*dep
      else
         dep2cogdep=(-7.556e-05*log(e0) + 0.002826999)* dep      
      endif

Here is the call graph for this function:

dep2s()

real function dep2s	(	real	E0,
		integer	NN,
		real	dep
	)

Definition at line 198 of file asinfo.f.

References dep2s().

!     ***************
      implicit none 
      real e0  !input GeV
      integer nn ! input  incident nucleon number. not used now
      real dep ! input dep g/cm^2
      
      real dep2cogdep, cogdep2s
      dep2s=cogdep2s(e0,nn, dep2cogdep(e0,nn, dep))

Here is the call graph for this function:

s2cogdep()

real function s2cogdep	(	real	E0,
		real	s
	)

Definition at line 241 of file asinfo.f.

References s2cogdep().

!     *************
      real e0    ! not used now
      real s     ! input. age
      s2cogdep= ( (s+9.228)/10.28 )**(1./0.0407)

Here is the call graph for this function:

s2ne()

real function s2ne	(	real	E0,
		real	s
	)

Definition at line 219 of file asinfo.f.

References s2ne().

!     ***************
      implicit none 
      real e0  !input GeV
      real s    ! input age
      
      real s2cogdep, cogdep2ne

      s2ne=cogdep2ne(e0, s2cogdep(e0, s))

Here is the call graph for this function:

s2ng()

real function s2ng	(	real	E0,
		real	s
	)

Definition at line 209 of file asinfo.f.

References s2ng().

!     ***************
      implicit none 
      real e0  !input total incident energy/n  GeV
      real s ! input  age
      
      real cogdep2ng, s2cogdep
      s2ng=cogdep2ng(e0,s2cogdep(e0,s))

Here is the call graph for this function:

s2nh()

real function s2nh	(	real	E0,
		real	s
	)

Definition at line 120 of file asinfo.f.

References s2nh().

!     ***********************
      implicit none
      real e0 ! input GeV
      real s  ! input age
!        use mu at present
      real s2nmu
      s2nh = s2nmu(e0, s)

Here is the call graph for this function:

s2nmu()

real function s2nmu	(	real	E0,
		real	s
	)

Definition at line 230 of file asinfo.f.

References s2nmu().

!     ***************
      implicit none 
      real e0  !input GeV
      real s    ! input age
      
      real s2cogdep, cogdep2nmu

      s2nmu=cogdep2nmu(e0,s2cogdep(e0, s))

Here is the call graph for this function:

sizex()

real function sizex	(	real	E0in,
		integer	NN,
		integer	code,
		real	s
	)

Definition at line 17 of file asinfo.f.

References sizex().

      implicit none
      real e0in ! input primary total energy in GeV
      integer nn ! input. incident nucleon number 0--> photon, 1-->proton, 4-->He ...56-->Fe
      integer code ! input code
      real s ! input age

      real s2ng, s2ne, s2nmu, s2nh, e0, fac
      if(nn .ge. 1) then
         e0=e0in/nn
         fac=nn
      else
         e0=e0in
         fac =1
      endif
      if(e0 .gt. 1.e11 .or. e0 .lt. 1.e7) then
         write(0,*)
     *   "*******warning: E0/n=",e0, " is out of range ******"
      endif
         
      if(code .eq. 1) then
         sizex=s2ng(e0, s)
         if(nn .eq. 0) fac = 1.2
      elseif(code .eq. 2) then
         sizex=s2ne(e0, s)
         if(nn .eq. 0) fac = 1.2
      elseif(code .eq. 3) then
         sizex=s2nmu(e0, s)
         if(nn .eq. 0) fac = 0.01
      elseif(code .le. 6) then
         sizex=s2nh(e0, s)
         if(nn .eq. 0) fac = 0.01
      else
         write(0,*) ' code error=',code, ' to sizex'
         stop 6776
      endif
      sizex=sizex*fac

Here is the call graph for this function: