Context Navigation

spectrum_init.F @ 114

Last change on this file since 114 was 114, checked in by moretto, 11 years ago
actual version of ESAF at CCin2p3
File size: 1.9 KB

Line
1	SUBROUTINE SPECTRUM_INIT
2	#include "event.inc"
3	#include "shower.inc"
4	#include "detector.inc"
5	CHARACTER*22 CHTAG(1)/'HILLAS PARAMETRIZATION'/
6
7	Emin = 1.e0 ! MeV
8	Emax = 1.e4 ! MeV
9	Ne = 10000
10	Ns = 20
11	Smin = 0.1
12	Smax = 2.0
13	step_e = (Emax-Emin)/(Ne-1)
14	step_s = (Smax-Smin)/(Ns-1)
15	Nint = Ne*Ns
16	i_portion = Nint/10
17	*
18	CALL BOOK_HISTOS('DUMP')
19	IF(SPEC.EQ.1) THEN
20	*
21	WRITE(6,*) ' '
22	WRITE(6,98)
23	*
24	DO k = 1, 1
25	l = 0
26	ID = 1000*k
27	WRITE(6,99) CHTAG(k)
28	DO i = 1, Ne
29	e = Emin + (i-1)*step_e
30	DO j = 1, Ns
31	s = Smin + (j-1)*step_s
32	l = l + 1
33	CALL HFILL(ID,s,e,ENERGY_THRESHOLD(s,e))
34	CALL HFILL(ID+1,s,e,ENERGY_SPECTRUM(s,e))
35	fraction_integral =
36	+ ENERGY_SPECTRUM_INTEGRAL(s,e)
37	CALL HFILL(ID+2,s,e,fraction_integral)
38	ntmp = l/i_portion
39	ntmp = l - ntmp*i_portion
40	IF(ntmp.EQ.0) WRITE(6,100)
41	ENDDO
42	ENDDO
43	WRITE(6,101)
44	WRITE(6,*) ' '
45	ENDDO
46
47	* ... NOW WRITE THE FRACTION OF ELECTRONS WITH
48	* ... THE ENERGY ABOVE THE CHERENKOV THRESHOLD
49
50	SPEC = 0
51	DO k = 1, 1
52	ID = 1000*k
53	DO i = 1, Ne
54	e = Emin + (i-1)*step_e
55	DO j = 1, Ns
56	s = Smin + (j-1)*step_s
57	ch_threshold = CHERENKOV_THRESHOLD(s,e)
58	CALL HFILL(ID+3,s,e,ch_threshold)
59	ENDDO
60	ENDDO
61	ENDDO
62	SPEC = 1
63	ENDIF
64	98 FORMAT(1x,'COMPUTING THE ENERGY DISTRIBUTION INTEGRALS: ')
65	99 FORMAT($,1x,'FOR ',A,11x,'(breath deeply)')
66	100 FORMAT($,1x,'.')
67	101 FORMAT($,1x,'[done]')
68	END

Note: See TracBrowser for help on using the repository browser.