| 1 |
|
|---|
| 2 | SUBROUTINE UGLAST
|
|---|
| 3 | *
|
|---|
| 4 | * Termination routine to print histograms and statistics
|
|---|
| 5 | *
|
|---|
| 6 | #include "geant321/gcflag.inc"
|
|---|
| 7 | #include "geant321/gckine.inc"
|
|---|
| 8 | #include "geant321/gctrak.inc"
|
|---|
| 9 | #include "calor.inc"
|
|---|
| 10 | #include "celoss.inc"
|
|---|
| 11 | #include "histo.inc"
|
|---|
| 12 | *
|
|---|
| 13 | CHARACTER*20 parnam,matnam
|
|---|
| 14 | CHARACTER*4 unit1,unit2
|
|---|
| 15 | *
|
|---|
| 16 | * *** run summary
|
|---|
| 17 | PRINT 749
|
|---|
| 18 | call GFPART (Ikine,parnam,ityp,dum,duc,dut,udu,idu)
|
|---|
| 19 | call GFMATE (Imat ,matnam,dua,duz,dud,dur,dui,udu,idu)
|
|---|
| 20 | CALL GEVKEV (pkine(1),Ekine,unit1)
|
|---|
| 21 | PRINT 750, Ievent,parnam,Ekine,unit1,sizeX,matnam
|
|---|
| 22 | PRINT 749
|
|---|
| 23 | *
|
|---|
| 24 | fnorm = 1./Ievent
|
|---|
| 25 | *
|
|---|
| 26 | * *** total energy deposit
|
|---|
| 27 | Etot = fnorm*Etot
|
|---|
| 28 | Etot2 = fnorm*Etot2
|
|---|
| 29 | rms = sqrt(abs(Etot2-Etot**2))
|
|---|
| 30 | *
|
|---|
| 31 | CALL GEVKEV (Etot,aveE,unit1)
|
|---|
| 32 | CALL GEVKEV (rms ,rmsE ,unit2)
|
|---|
| 33 | PRINT 752, aveE,unit1,rmsE,unit2
|
|---|
| 34 | *
|
|---|
| 35 | * *** track length of primary track
|
|---|
| 36 | TrLen = fnorm*TrLen
|
|---|
| 37 | TrLen2 = fnorm*TrLen2
|
|---|
| 38 | rms = sqrt(abs(TrLen2-TrLen**2))
|
|---|
| 39 | PRINT 753, TrLen,rms
|
|---|
| 40 | *
|
|---|
| 41 | * *** compare with csda range from table
|
|---|
| 42 | call GFTMAT(Imat,ikine,'RANG',1,pkine(1),csdaRang,dum,ixst)
|
|---|
| 43 | PRINT 763, csdaRang
|
|---|
| 44 | *
|
|---|
| 45 | * *** projected range of primary track
|
|---|
| 46 | XTrLen = fnorm*XTrLen
|
|---|
| 47 | XTrLen2 = fnorm*XTrLen2
|
|---|
| 48 | rms = sqrt(abs(XTrLen2-XTrLen**2))
|
|---|
| 49 | PRINT 754, XTrLen,rms
|
|---|
| 50 | *
|
|---|
| 51 | * *** nb of steps and step size of primary track
|
|---|
| 52 | fNbStep = fnorm*NbStep
|
|---|
| 53 | fNbStep2 = fnorm*NbStep2
|
|---|
| 54 | rms1 = sqrt(abs(fNbStep2-fNbStep**2))
|
|---|
| 55 | *
|
|---|
| 56 | StepLen = fnorm*StepLen
|
|---|
| 57 | StepLen2 = fnorm*StepLen2
|
|---|
| 58 | rms2 = sqrt(abs(StepLen2-StepLen**2))
|
|---|
| 59 | PRINT 755, fNbStep,rms1,StepLen,rms2
|
|---|
| 60 | *
|
|---|
| 61 | * *** transmission coefficients
|
|---|
| 62 | absorbed = fnorm*100*kcoef(1)
|
|---|
| 63 | transmit = fnorm*100*kcoef(2)
|
|---|
| 64 | reflected = fnorm*100*kcoef(3)
|
|---|
| 65 | PRINT 756, absorbed, transmit, reflected
|
|---|
| 66 | *
|
|---|
| 67 | PRINT 749
|
|---|
| 68 | *
|
|---|
| 69 | * *** geant termination
|
|---|
| 70 | CALL GLAST
|
|---|
| 71 | *
|
|---|
| 72 | * *** close HIGZ
|
|---|
| 73 | CALL HPLEND
|
|---|
| 74 | *
|
|---|
| 75 | * *** Save histo
|
|---|
| 76 | if (iswit(2).eq.1) then
|
|---|
| 77 | write(6,99) fileName
|
|---|
| 78 | 99 format(/,15x,'histogram file --> Name: ',A25)
|
|---|
| 79 | CALL HRPUT(0,fileName,'N')
|
|---|
| 80 | endif
|
|---|
| 81 | *
|
|---|
| 82 | * *** formats
|
|---|
| 83 | 749 FORMAT(/, 60(1H-),/)
|
|---|
| 84 | 750 FORMAT( 1X,'The run is',I8,1X,A10,' of',F8.2,A4,' through ',
|
|---|
| 85 | & E12.4,' cm of ',A10)
|
|---|
| 86 | 752 FORMAT( 1X,'Total Energy deposited = ',
|
|---|
| 87 | & F8.3, A4,' +- ', F8.3, A4)
|
|---|
| 88 | 753 FORMAT(/,1X,'Track length of primary track = ',
|
|---|
| 89 | & E12.5,' cm +- ', E12.5,' cm')
|
|---|
| 90 | 763 FORMAT( 1X,'csda Range from table = ',
|
|---|
| 91 | & E12.5,' cm')
|
|---|
| 92 | 754 FORMAT(/,1X,'Projected range = ',
|
|---|
| 93 | & E12.5,' cm +- ', E12.5,' cm')
|
|---|
| 94 | 755 FORMAT(/,1X,'Nb of steps of primary track = ',
|
|---|
| 95 | & F6.2,' +- ', F6.2,
|
|---|
| 96 | & 3X,'Step size = ', E12.5,' cm +- ', E12.5,' cm')
|
|---|
| 97 | 756 FORMAT(/,1X,'absorbed = ',F6.2,' %',3X,'transmit = ', F6.2,' %',
|
|---|
| 98 | & 3X,'reflected = ', F6.2,' %')
|
|---|
| 99 | *
|
|---|
| 100 | END
|
|---|
