| 1 | |
|---|
| 2 | SUBROUTINE UGLAST |
|---|
| 3 | * |
|---|
| 4 | * Termination routine to print histograms and statistics |
|---|
| 5 | * |
|---|
| 6 | #include "geant321/gcbank.inc" |
|---|
| 7 | #include "geant321/gcflag.inc" |
|---|
| 8 | #include "geant321/gckine.inc" |
|---|
| 9 | #include "pvolum.inc" |
|---|
| 10 | #include "celoss.inc" |
|---|
| 11 | * |
|---|
| 12 | DIMENSION XSEL1(NBIN),XSEL1C(NBIN),XSER1(NBIN),XSER1C(NBIN), |
|---|
| 13 | + XSEL2(NBIN),XSEL2C(NBIN),XSER2(NBIN),XSER2C(NBIN) |
|---|
| 14 | |
|---|
| 15 | DIMENSION EDIST(100),EDISTC(100) |
|---|
| 16 | * |
|---|
| 17 | * *** Normalize and print energy distribution |
|---|
| 18 | XEVENT=IEVENT |
|---|
| 19 | CNORM = 100./(XEVENT*PKINE(3)) |
|---|
| 20 | * |
|---|
| 21 | * *** mean total energy deposit by charged and by neutral |
|---|
| 22 | EDEPCH = CNORM*EDEPCH |
|---|
| 23 | EDEPNE = CNORM*EDEPNE |
|---|
| 24 | * |
|---|
| 25 | * *** longitudinal profile |
|---|
| 26 | DO 2 I = 1,NLTOT |
|---|
| 27 | XSEL1 (I) = CNORM * SEL1 (I) |
|---|
| 28 | XSEL2 (I) = CNORM*SQRT(ABS(XEVENT*SEL2 (I) - SEL1 (I)**2)) |
|---|
| 29 | XSEL1C(I) = CNORM * SEL1C(I) |
|---|
| 30 | XSEL2C(I) = CNORM*SQRT(ABS(XEVENT*SEL2C(I) - SEL1C(I)**2)) |
|---|
| 31 | 2 CONTINUE |
|---|
| 32 | CALL HPAK (5,XSEL2) |
|---|
| 33 | CALL HPAK (7,XSEL2C) |
|---|
| 34 | * |
|---|
| 35 | * *** radial profile |
|---|
| 36 | DO 3 I = 1,NRTOT |
|---|
| 37 | XSER1 (I) = CNORM * SER1 (I) |
|---|
| 38 | XSER2 (I) = CNORM*SQRT(ABS(XEVENT*SER2 (I) - SER1 (I)**2)) |
|---|
| 39 | XSER1C(I) = CNORM * SER1C(I) |
|---|
| 40 | XSER2C(I) = CNORM*SQRT(ABS(XEVENT*SER2C(I) - SER1C(I)**2)) |
|---|
| 41 | 3 CONTINUE |
|---|
| 42 | CALL HPAK ( 9,XSER2) |
|---|
| 43 | CALL HPAK (11,XSER2C) |
|---|
| 44 | * |
|---|
| 45 | * *** total track length |
|---|
| 46 | CNORM = 1./(XEVENT*X0) |
|---|
| 47 | XTRCH1 = CNORM*STRCH1 |
|---|
| 48 | XTRCH2 = CNORM*SQRT(ABS(XEVENT*STRCH2 - STRCH1**2)) |
|---|
| 49 | XTRNE1 = CNORM*STRNE1 |
|---|
| 50 | XTRNE2 = CNORM*SQRT(ABS(XEVENT*STRNE2 - STRNE1**2)) |
|---|
| 51 | * |
|---|
| 52 | * *** Print profiles (under condition iswit(2).gt.0) |
|---|
| 53 | * |
|---|
| 54 | if (iswit(2).gt.0) then |
|---|
| 55 | PRINT 749 |
|---|
| 56 | PRINT 750 |
|---|
| 57 | PRINT 751 |
|---|
| 58 | DO 15 I=1,NLTOT |
|---|
| 59 | B0 = (I-1)*DLX0 |
|---|
| 60 | B1 = I*DLX0 |
|---|
| 61 | PRINT 754,B0,B1,XSEL1(I),XSEL2(I),B1,XSEL1C(I),XSEL2C(I) |
|---|
| 62 | 15 CONTINUE |
|---|
| 63 | |
|---|
| 64 | PRINT 760 |
|---|
| 65 | PRINT 751 |
|---|
| 66 | DO 16 I=1,NRTOT |
|---|
| 67 | B0 = (I-1)*DRX0 |
|---|
| 68 | B1 = I*DRX0 |
|---|
| 69 | PRINT 754,B0,B1,XSER1(I),XSER2(I),B1,XSER1C(I),XSER2C(I) |
|---|
| 70 | 16 CONTINUE |
|---|
| 71 | endif |
|---|
| 72 | * |
|---|
| 73 | * *** normalize histo of energy ditribution of contributing particles |
|---|
| 74 | * and compute cumulative distribution |
|---|
| 75 | SUM = HSUM (21) |
|---|
| 76 | CALL HUNPAK(21,EDIST,'HIST',1) |
|---|
| 77 | EDIST( 1) = EDIST(1)*100/SUM |
|---|
| 78 | EDISTC(1) = EDIST(1) |
|---|
| 79 | DO 17 I=2,100 |
|---|
| 80 | EDIST (I) = EDIST (I)*100/SUM |
|---|
| 81 | EDISTC(I) = EDISTC(I-1) + EDIST(I) |
|---|
| 82 | 17 CONTINUE |
|---|
| 83 | CALL HPAK (21,EDIST) |
|---|
| 84 | CALL HPAK (22,EDISTC) |
|---|
| 85 | * |
|---|
| 86 | * *** print summary |
|---|
| 87 | PRINT 770 |
|---|
| 88 | PRINT 771,XSEL1C(NLTOT),XSEL2C(NLTOT) |
|---|
| 89 | PRINT 774,EDEPCH |
|---|
| 90 | PRINT 775,EDEPNE |
|---|
| 91 | PRINT 772,XTRCH1,XTRCH2 |
|---|
| 92 | PRINT 773,XTRNE1,XTRNE2 |
|---|
| 93 | PRINT 749 |
|---|
| 94 | * |
|---|
| 95 | * *** Save selected histograms |
|---|
| 96 | CALL HRPUT(0,'testem2.hbook','N') |
|---|
| 97 | * |
|---|
| 98 | * *** terminaison |
|---|
| 99 | CALL GLAST |
|---|
| 100 | * |
|---|
| 101 | * *** close HIGZ |
|---|
| 102 | CALL HPLEND |
|---|
| 103 | * |
|---|
| 104 | 749 FORMAT(//) |
|---|
| 105 | 750 FORMAT(15X,'LATERAL PROFILE',35X,'CUMULATIVE LATERAL PROFILE'/) |
|---|
| 106 | 751 FORMAT( 8X,'Bin',12X,' Mean ',5X,' rms', |
|---|
| 107 | * 19X,'Bin', 9X,' Mean ',5X,' rms',/) |
|---|
| 108 | 754 FORMAT( 3X,F5.2,'->',F5.2,' radl: ',F7.2,'% ',F7.2,'%', |
|---|
| 109 | * 13X, '0->',F5.2,' radl: ',F7.2,'% ',F7.2,'%') |
|---|
| 110 | 760 FORMAT(///,15X,'RADIAL PROFILE',35X,'CUMULATIVE RADIAL PROFILE'/) |
|---|
| 111 | 770 FORMAT(/,30X,'SUMMARY',/) |
|---|
| 112 | 771 FORMAT( 25X,'energy deposit : ',F7.2,' % E0 +- ',F7.2,' % E0') |
|---|
| 113 | 772 FORMAT( 25X,'charged traklen: ',F7.2,' radl +- ',F7.2,' radl') |
|---|
| 114 | 773 FORMAT( 25X,'neutral traklen: ',F7.2,' radl +- ',F7.2,' radl') |
|---|
| 115 | 774 FORMAT( 25X,'edep by charged: ',F7.2,' % E0') |
|---|
| 116 | 775 FORMAT( 25X,'edep by neutral: ',F7.2,' % E0') |
|---|
| 117 | * |
|---|
| 118 | END |
|---|
| 119 | |
|---|
