1 | SUBROUTINE GENERATE_ESAFEVENT(Lmax,KEY_TRIG) |
---|
2 | IMPLICIT NONE |
---|
3 | #include "detector.inc" |
---|
4 | #include "event.inc" |
---|
5 | #include "cwntvars.inc" |
---|
6 | #include "constants.inc" |
---|
7 | |
---|
8 | INTEGER FLUOR,CHER,i |
---|
9 | REAL Lmin,Lmax,th1,ph1,z,time |
---|
10 | LOGICAL KEY_TRIG |
---|
11 | * |
---|
12 | CALL TIMED(time) |
---|
13 | th1 = THETA/PI*180 |
---|
14 | ph1 = PHI/PI*180 |
---|
15 | z = RINT(3) - R_EARTH |
---|
16 | CALL DEVELOP_ESAFSHOWER(Lmax,KEY_TRIG) |
---|
17 | IF(KEY_TRIG) THEN |
---|
18 | CALL TIMED(time) |
---|
19 | WRITE(6,1) time |
---|
20 | WRITE(6,11) |
---|
21 | WRITE(6,12) EINT |
---|
22 | WRITE(6,15) RINT(1) |
---|
23 | WRITE(6,13) TH1 |
---|
24 | WRITE(6,16) RINT(2) |
---|
25 | WRITE(6,14) PH1 |
---|
26 | WRITE(6,17) z |
---|
27 | ENDIF |
---|
28 | 1 FORMAT('SLAST GENERATED EVENT IN:',1x,F6.3,' sec.') |
---|
29 | 11 FORMAT($,'EVENT PARAMETERS:') |
---|
30 | 12 FORMAT($,1x,'ENERGY',2x,E7.1,' eV.') |
---|
31 | 13 FORMAT($,18x,'THETA ',2x,F6.1,' deg.') |
---|
32 | 14 FORMAT($,18x,'PHI ',2x,F6.1,' deg.') |
---|
33 | 15 FORMAT(2x,'Xint ',2x,F9.3,' km.') |
---|
34 | 16 FORMAT(2x,'Yint ',2x,F9.3,' km.') |
---|
35 | 17 FORMAT(2x,'Zint ',2x,F9.3,' km.') |
---|
36 | END |
---|