1 | |
---|
2 | SUBROUTINE GUOUT |
---|
3 | * |
---|
4 | * User routine called at the end of each event |
---|
5 | * |
---|
6 | #include "geant321/gcflag.inc" |
---|
7 | #include "geant321/gckine.inc" |
---|
8 | #include "pvolum.inc" |
---|
9 | #include "celoss.inc" |
---|
10 | * |
---|
11 | * *** drawing |
---|
12 | * |
---|
13 | #ifndef batch |
---|
14 | IF (ISWIT(1).NE.0) THEN |
---|
15 | CALL GDHEAD (110110,'testem2',0.) |
---|
16 | CALL GDSHOW (1) |
---|
17 | CALL GDXYZ (0) |
---|
18 | END IF |
---|
19 | #endif |
---|
20 | * |
---|
21 | * |
---|
22 | * *** statistic |
---|
23 | * |
---|
24 | DLC = 0. |
---|
25 | DRC = 0. |
---|
26 | |
---|
27 | * longitudinal profile |
---|
28 | * |
---|
29 | DO 2 I = 1,NLTOT |
---|
30 | SEL1 (I) = SEL1 (I) + DEDL(I) |
---|
31 | SEL2 (I) = SEL2 (I) + DEDL(I)**2 |
---|
32 | DLC = DLC + DEDL(I) |
---|
33 | SEL1C(I) = SEL1C(I) + DLC |
---|
34 | SEL2C(I) = SEL2C(I) + DLC**2 |
---|
35 | |
---|
36 | BIN = (FLOAT(I)-0.5)*DLX0 |
---|
37 | CALL HFILL(4,BIN,100*DEDL(I)/(DLX0*PKINE(3)),1.) |
---|
38 | BIN = FLOAT(I)*DLX0 |
---|
39 | CALL HFILL(6,BIN,100*DLC /PKINE(3),1.) |
---|
40 | 2 CONTINUE |
---|
41 | |
---|
42 | * radial profile |
---|
43 | * |
---|
44 | DO 3 I = 1,NRTOT |
---|
45 | SER1 (I) = SER1 (I) + DEDR(I) |
---|
46 | SER2 (I) = SER2 (I) + DEDR(I)**2 |
---|
47 | DRC = DRC + DEDR(I) |
---|
48 | SER1C(I) = SER1C(I) + DRC |
---|
49 | SER2C(I) = SER2C(I) + DRC**2 |
---|
50 | |
---|
51 | BIN = (FLOAT(I)-0.5)*DRX0 |
---|
52 | CALL HFILL(8,BIN,100*DEDR(I)/(DRX0*PKINE(3)),1.) |
---|
53 | BIN = FLOAT(I)*DRX0 |
---|
54 | CALL HFILL(10,BIN,100*DRC /PKINE(3),1.) |
---|
55 | 3 CONTINUE |
---|
56 | * |
---|
57 | * energy deposited and track length |
---|
58 | * |
---|
59 | ESEEN = 100.*DLC/PKINE(3) |
---|
60 | CALL HFILL(1, ESEEN,0.,1.) |
---|
61 | CALL HFILL(2,STRCH/X0,0.,1.) |
---|
62 | CALL HFILL(3,STRNE/X0,0.,1.) |
---|
63 | * |
---|
64 | STRCH1 = STRCH1 + STRCH |
---|
65 | STRCH2 = STRCH2 + STRCH**2 |
---|
66 | STRNE1 = STRNE1 + STRNE |
---|
67 | STRNE2 = STRNE2 + STRNE**2 |
---|
68 | * |
---|
69 | END |
---|