[807] | 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 |
---|