cprimAcceptance.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cprimacceptance (comp, e_or_p, rigth, cos1, cos2, fai1, fai2, prob)
real *8 function	crigfunc1 (zen)
real *8 function	crigfunc2 (phi)
subroutine	csetcosdeg (cosin, degin)

Function/Subroutine Documentation

cprimacceptance()

subroutine cprimacceptance	(	type (component)	comp,
		real*8	e_or_p,
		real*8	rigth,
		real*8	cos1,
		real*8	cos2,
		real*8	fai1,
		real*8	fai2,
		real*8	prob
	)

Definition at line 3 of file cprimAcceptance.f.

References cconv_prim_e(), crigfunc1(), d, and kdexpintf().

Referenced by __getacceptance.f__(), and cprimflux().

#include  "Zglobalc.h"
#include  "Zptcl.h"
#include  "Zprimary.h"

      real*8 azmmin, azmmax, rig, cosx
      logical degree, cosfactor
      common /zpirmflux/ azmmin, azmmax, rig, cosx, 
     * degree, cosfactor


      type(component)::comp  ! input primary component
      real*8 e_or_p     ! input.  E or p as given in primary file
      real*8 rigth      ! input.  threshold rigidty below which geomagneic
                        !         effect appears.  make it 0 if
                        !         no rigidity cut.
      real*8 cos1, cos2 ! input.  cos zenith range (cos1 < cos2)
      real*8 fai1, fai2 ! input.  azimuthal angle range (fai1< fai2) deg.
      real*8 prob       ! output. average prob. that the primary can come

      type(ptcl):: aptcl
!

      external crigfunc1
      real*8 crigfunc1
      real*8  error, ans

      integer icon
      data eps/1.d-4/



      azmmin = fai1
      azmmax = fai2
      call cconv_prim_e(comp, e_or_p, aptcl)
!         rigidity
      rig =
     *   sqrt( aptcl%fm%p(4)**2 - aptcl%mass**2 )/aptcl%charge

      if(rig .lt. rigth) then
         call kdexpintf(crigfunc1, cos1, cos2, eps, ans, error, icon)
         prob= ans/( (fai2-fai1) * (cos2 - cos1 ))   ! fai is in deg. o.k
      else
         prob = 1.0
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

crigfunc1()

real*8 function crigfunc1

(

real*8

zen

)

Definition at line 50 of file cprimAcceptance.f.

References crigfunc2(), and k16pgaussleg().

Referenced by cprimacceptance().

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

      real*8 azmmin, azmmax, rig, cosx
      logical degree, cosfactor
      common /zpirmflux/ azmmin, azmmax, rig, cosx, 
     * degree, cosfactor

      real*8 zen

      real*8 ans


      external crigfunc2
      real*8 crigfunc2

!        integrate over phi = azmmin, azmmax
      cosx = zen
      call k16pgaussleg(crigfunc2, azmmin, azmmax, 16, ans)
      crigfunc1 = ans 
      if(cosfactor) then
         crigfunc1 = crigfunc1 * cosx
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

crigfunc2()

real*8 function crigfunc2

(

real*8

phi

)

Definition at line 78 of file cprimAcceptance.f.

References crigcut().

Referenced by crigfunc1().

      implicit none
#include "Zglobalc.h"
#include "Zptcl.h"
#include "Zprimary.h"
      real*8 phi   ! iput in degree
!         for fixed cos and phi; get  R(theta, fai, rig)

      real*8 azmmin, azmmax, rig, cosx
      logical degree, cosfactor
      common /zpirmflux/ azmmin, azmmax, rig, cosx, 
     * degree, cosfactor

      real*8 temp
!
      real*8 prob
!
      if(degree) then
         temp = acos(cosx)/torad
      else
         temp = cosx
      endif
      call crigcut(phi,  temp,  rig, prob) 
      crigfunc2 = prob

Here is the call graph for this function:

Here is the caller graph for this function:

csetcosdeg()

subroutine csetcosdeg	(	logical	cosin,
		logical	degin
	)

Definition at line 104 of file cprimAcceptance.f.

Referenced by __getacceptance.f__(), __getst.f__(), and __showspec.f__().

      implicit none
      logical cosin, degin
      real*8 azmmin, azmmax, rig, cosx
      logical degree, cosfactor
      common /zpirmflux/ azmmin, azmmax, rig, cosx, 
     * degree, cosfactor
      
      cosfactor = cosin
      degree = degin

Here is the caller graph for this function: