getST2.f File Reference

Go to the source code of this file.

Functions/Subroutines
program	__getst2.f__
subroutine	inteflux (comp, ans)
real *8 function	primdn (eorp)
subroutine	chooktrace
subroutine	chookceren
subroutine	chookcerens
subroutine	chookcerene
subroutine	chookbgrun

Function/Subroutine Documentation

__getst2.f__()

program __getst2.f__

(

)

Definition at line 23 of file getST2.f.

References azmmax, azmmin, cbeginrun(), coszenith, cprintprim(), creadparam(), d0, false, i, inteflux(), rigc, true, zen1, and zen2.

      integer  i

Here is the call graph for this function:

chookbgrun()

subroutine chookbgrun

(

)

Definition at line 191 of file getST2.f.

chookceren()

subroutine chookceren

(

)

Definition at line 185 of file getST2.f.

chookcerene()

subroutine chookcerene

(

)

Definition at line 189 of file getST2.f.

chookcerens()

subroutine chookcerens

(

)

Definition at line 187 of file getST2.f.

chooktrace()

subroutine chooktrace

(

)

Definition at line 183 of file getST2.f.

inteflux()

subroutine inteflux	(	type (component)	comp,
		real*8	ans
	)

Definition at line 102 of file getST2.f.

References cconv_prim_e2(), cmkptc(), d, inteprim2(), kdexpintfb(), kdwhereis(), primdn(), and rigc.

      implicit none
#include "Zglobalc.h"
#include  "Zptcl.h"
#include  "Zprimary.h"
      type(component):: comp
      real*8 ans
      include 'Zflux.h'
      
      real*8 eps, error, ans2, eth, e_or_p
      integer icon


      integer imax
      external primdn
      real*8 primdn
      type(ptcl)::aptcl
      integer j
      data eps/1.d-4/



      compx = comp
      imax = comp%no_of_seg
      if(rigc .eq. 0.) then
!         no rigidy  cut. integrate segmented power functions
         call inteprim2(comp, 1, imax, ans)
      else
         call cmkptc(comp%code, comp%subcode, comp%charge, aptcl)
         eth =sqrt( (rigc*comp%charge)**2 + aptcl%mass**2 )
         call cconv_prim_e2(comp, eth, e_or_p)
         call kdwhereis(e_or_p, comp%no_of_seg+1, comp%energy, 1, j)
         if(j .le. imax) then
!              energy integral; 
            call kdexpintfb(primdn, comp%energy(1), comp%energy(j+1),
     *           eps,  ans, error,  icon)
            ans = ans 
         endif
!             add E> comp.energy(j+1)
         if(j+1 .lt. imax ) then
            call inteprim2(comp, j+1, imax, ans2)
         else
            ans2 = 0.
         endif
         ans = ans + ans2
      endif

Here is the call graph for this function:

primdn()

real*8 function primdn

(

real*8

eorp

)

Definition at line 151 of file getST2.f.

References azmmin, cconv_prim_e(), cprimflux0(), crigcut(), degree, e, primdn(), and zen1.

      implicit none

#include "Zglobalc.h"
#include "Zptcl.h"
#include "Zprimary.h"
!          primary flux at E       
      real*8 eorp

      include 'Zflux.h'



      type(ptcl):: aptcl
      real*8 rig,  prob, flux, zeny

      e = eorp
!          E= E_or_P  must be converted to rigidity
      call cconv_prim_e(compx, e, aptcl)

      rig =sqrt( aptcl%fm%p(4)**2  - aptcl%mass**2)/aptcl%charge
      if(degree) then
         zeny = zen1/torad
      else
         zeny = zen1
      endif
      call crigcut(azmmin, zeny, rig, prob)
      call cprimflux0(compx, e, flux)  ! E = EorP

      primdn = flux * prob 

Here is the call graph for this function: