getST.f File Reference

Go to the source code of this file.

Functions/Subroutines
program	__getst.f__
subroutine	inteflux (comp, ans)
real *8 function	primdn (eorp)
real *8 function	funczen (zenx)
real *8 function	funcazim (azm)
subroutine	chooktrace
subroutine	chookceren
subroutine	chookcerens
subroutine	chookcerene
subroutine	chookbgrun
subroutine	csetcosdeg (cosin, degin)

Function/Subroutine Documentation

__getst.f__()

program __getst.f__

(

)

Definition at line 20 of file getST.f.

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

      integer  i

Here is the call graph for this function:

chookbgrun()

subroutine chookbgrun

(

)

Definition at line 329 of file getST.f.

chookceren()

subroutine chookceren

(

)

Definition at line 323 of file getST.f.

chookcerene()

subroutine chookcerene

(

)

Definition at line 327 of file getST.f.

chookcerens()

subroutine chookcerens

(

)

Definition at line 325 of file getST.f.

chooktrace()

subroutine chooktrace

(

)

Definition at line 321 of file getST.f.

csetcosdeg()

subroutine csetcosdeg	(	logical	cosin,
		logical	degin
	)

Definition at line 332 of file getST.f.

References degree.

      implicit none
      logical cosin, degin
#include "Zptcl.h"
#include "Zprimary.h"

      include 'Zflux.h'
      
      cosfactor = cosin
      degree = degin

funcazim()

real*8 function funcazim

(

real*8

azm

)

Definition at line 286 of file getST.f.

References cconv_prim_e(), cprimflux0(), crigcut(), d0, degree, e, and zen.

Referenced by funcazim().

      implicit none
#include "Zglobalc.h"
#include "Zptcl.h"
#include "Zprimary.h"
      include 'Zflux.h'

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

!          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 = zen/torad
      else
         zeny = zen
      endif
      call crigcut(azm, zeny, rig, prob)
!/////////////
!      prob = 1.
!/////////////
      call cprimflux0(compx, e, flux)  ! E = EorP
      if(cosfactor .eq. 0 ) then
         funcazim = flux * prob 
      elseif(cosfactor .eq. 1) then
         funcazim = flux * prob * zen
      elseif(cosfactor .eq. 2) then
         funcazim = flux * prob *(1.0d0+ zen)/2.d0
      else
         funcazim = flux * prob *
     *   (zen + ratio*sqrt(1.-zen**2)) / sqrt(1.+ratio**2)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

funczen()

real*8 function funczen

(

real*8

zenx

)

Definition at line 266 of file getST.f.

References azmmax, azmmin, d, kdexpintf(), and zen.

      implicit none
#include "Zptcl.h"
#include "Zprimary.h"

      real*8 zenx

      include 'Zflux.h'

      integer icon
      real*8 eps, ans, error
      real*8 funcazim
      external funcazim
      data  eps/1.d-5/

      zen = zenx
      call kdexpintf(funcazim, azmmin, azmmax, eps, ans, error, icon)
      funczen = ans

Here is the call graph for this function:

inteflux()

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

Definition at line 150 of file getST.f.

References azmmax, azmmin, cconv_prim_e2(), cmkptc(), d, d0, inteprim2(), kdexpintfb(), kdwhereis(), primdn(), rigc, zen1, and zen2.

Referenced by __getst.f__(), __getst2.f__(), and __getst3.f__().

      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)
         if(cosfactor .eq. 0) then
           if(abs(ans/comp%inte_value-1.d0) .gt. 1.d-3 ) then
              write(*,*) ' ans=',ans, ' internal integral=',
     *                    comp%inte_value
              stop 9999
           endif
           ans = ans * (zen2-zen1)* abs(azmmax - azmmin)*torad
         elseif(cosfactor .eq. 1) then
!             take into account the horizontal  area.
            ans = comp%inte_value * (zen2**2- zen1**2)/2 *
     *                   abs(azmmax - azmmin)*torad
         elseif(cosfactor .eq. 2) then
!              
            ans = comp%inte_value
     *           *(zen2-zen1)/2.d0* (1.d0 + (zen1 + zen2)/2)*
     *           abs(azmmax - azmmin)*torad
         elseif(cosfactor .eq. 3) then
!                  cos + rat*sin
            ans = comp%inte_value *
     *       (  (zen2-zen1)/2.0d0 + ratio* (
     *           acos(zen1)-acos(zen2) +
     *           0.5*( zen2*sqrt(1.-zen2**2) -
     *                 zen1*sqrt(1.-zen1**2) )) )
     *           /sqrt(1.+ratio**2)
     *          * abs(azmmax - azmmin)*torad
         else 
            write(*,*) ' error of cosin'
            stop 9999
         endif   
      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 * torad
         endif
!             add E> comp.energy(j+1)
         if(j+1 .lt. imax ) then
            call inteprim2(comp, j+1, imax, ans2)
            if(cosfactor .eq. 0 ) then
               ans2 = ans2 * (zen2- zen1)*abs(azmmax - azmmin)*torad
            elseif(cosfactor .eq. 1) then
               ans2 = ans2 * (zen2**2- zen1**2)/2 *
     *                  abs(azmmax - azmmin)*torad
            elseif(cosfactor .eq. 2) then
!              fai is from 0 to 2pi; 
               ans2 =
     *          ans2 *(zen2-zen1)/2.d0* (1.d0 + (zen1 + zen2)/2.d0)
            else
                ans2 = ans2 *
     *          (  (zen2-zen1)/2.0d0 + ratio* (
     *           acos(zen1)-acos(zen2) +
     *           0.5*( zen2*sqrt(1.-zen2**2) -
     *                 zen1*sqrt(1.-zen1**2) )) )
     *           /sqrt(1.+ratio**2)
     *           *  abs(azmmax - azmmin)*torad
            endif
         else
            ans2 = 0.
         endif
         ans = ans + ans2
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

primdn()

real*8 function primdn

(

real*8

eorp

)

Definition at line 246 of file getST.f.

References e, k16pgaussleg(), zen1, and zen2.

Referenced by inteflux(), and primdn().

      implicit none

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

      include 'Zflux.h'

      real*8 funczen, ans
      external funczen

      e = eorp
!         integral on  zenith 
      call k16pgaussleg(funczen, zen1, zen2, 16, ans)
      primdn = ans

Here is the call graph for this function:

Here is the caller graph for this function: