source: JEM-EUSO/esaf_lal/tags/v1_r0/esaf/packages/simulation/externals/slast/src/cherenkov/energy_spectrum_integral.F @ 117

Last change on this file since 117 was 117, checked in by moretto, 11 years ago

ESAF version compilable on mac OS
File size: 869 bytes
Line 
1      REAL FUNCTION ENERGY_SPECTRUM_INTEGRAL(S,ET)
2#include "shower.inc"
3#include "detector.inc"
4      EXTERNAL FunINT
5      COMMON/SPEC_INT/AGE,E
6      PARAMETER (zero=0., one=1., NSEG = 1, RELTOL = 1.e-8, ABSTOL = 0.)
7      REAL A,B,S,ET
8      INTEGER i_aux_hbk_spec,ID
9      SAVE i_aux_hbk_spec
10      DATA i_aux_hbk_spec/0/
11*
12      AGE = S
13      E   = ET
14      A = zero
15      B = one
16      ID = 1000*DIST
17*
18c      IF (SPEC.EQ.1) THEN
19         CALL RADAPT(FunINT,A,B,NSEG,RELTOL,ABSTOL,RES,ERR)
20         ENERGY_SPECTRUM_INTEGRAL =  RES
21c      ELSE
22c         ENERGY_SPECTRUM_INTEGRAL = HXY(ID+2,S,ET)
23c      ENDIF
24 1000 CONTINUE
25      END
26
27      REAL FUNCTION FunINT(y)
28      IMPLICIT NONE
29#include "shower.inc"
30      REAL y,S,ET,ENERGY_SPECTRUM_HILLAS,f
31      COMMON/SPEC_INT/S,ET
32      IF(DIST.EQ.1) f = ENERGY_SPECTRUM_HILLAS(S,ET/y)
33      FunINT = ET*f
34      END
Note: See TracBrowser for help on using the repository browser.