getAcceptance.f

Go to the documentation of this file.

#include  "BlockData/cblkGene.h"
      implicit none
#include "Zglobalc.h"
#include "Zmanagerp.h"
#include "ZrigCut.h"
#include "Zptcl.h"
#include "Zobs.h"
#include "Zobsp.h"
#include "Zprimary.h"
#include "Zprimaryc.h"
#include "Zprimaryv.h"
#include "Zincidentp.h"
#ifdef NEXT486
#define IMAG_P dimag
#else
#define IMAG_P imag
#endif


      integer  i, j

      real*8  prob
      real*8  rigc
      real*8 E, azmmin, azmmax, dE, zen1, zen2

      call creadparam(5)
      call cbeginrun
      call cprintprim(errorout)
      write(errorout,*) ' XaxisFromSouth=',xaxisfromsouth
      zen2 = imag_p(coszenith)
      zen1 = real(coszenith)
      azmmax =imag_p(azimuth) + xaxisfromsouth
      azmmin =real(Azimuth) + XaxisFromSouth
      de=10.**0.05
      if(cutofffile .eq. ' ') then
         rigc = 0.
         call csetcosdeg(.false., .false.)
      else
         if(abs(latit-latitofsite) .gt. 5. .or. 
     *        abs(longi-longitofsite) .gt. 5. ) then
            write(errorout, *) 
     *     'check LatitOfSite and LongitOfSite; they differ'//
     *     ' from those in CutoffFile'
            stop 999
         endif
         rigc = 100.
         if(zenvalue .eq. 'deg') then
            call csetcosdeg(.false., .true.)
         else
            call csetcosdeg(.false., .false.)
         endif
      endif
      do i = 1, prim%no_of_comps
         e = prim%each(i)%energy(1)
         do j = 1, 10000
            call cprimacceptance(prim%each(i), e, rigc, zen1, zen2,
     *       azmmin, azmmax,  prob)
            write(*,*) prim%each(i)%label, sngl(e), sngl(prob)
            e = e*de
            if(e .gt. prim%each(i)%energy(prim%each(i)%no_of_seg+1))
     *          goto 100            
         enddo
 100     continue
         write(*,*)
      enddo
      end
      subroutine chooktrace
      end
      subroutine chookceren
      end
      subroutine chookcerens
      end
      subroutine chookcerene
      end
      subroutine chookbgrun
      end