| 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(5,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(7,BIN,100*DEDR(I)/(DRX0*PKINE(3)),1.) |
|---|
| 53 | BIN = FLOAT(I)*DRX0 |
|---|
| 54 | CALL HFILL(8,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 |
|---|
