| 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
|
|---|
