geant3/src/guout.F

0001
0002       SUBROUTINE GUOUT
0003 *
0004 *     User routine called at the end of each event
0005 *
0006 #include "geant321/gcflag.inc"
0007 #include "calor.inc"
0008 #include "celoss.inc"
0009 #include "histo.inc"
0010 *
0011       CHARACTER*4 unit
0012
0013 *
0014 *     histo: energy deposit per tracking medium
0015       do k=1,NbAbsor
0016         if (histo(k)) call hfill (k,EdepAbs(k)/histUnit(k),0.,1.)
0017       enddo
0018 *
0019 *     cumul energy deposited, charged track length and energy leakage
0020       do k=1,NbAbsor
0021         sumEdep(k) = sumEdep(k) + EdepAbs(k)
0022         sumTrck(k) = sumTrck(k) + TrckAbs(k)
0023 *
0024         su2Edep(k) = su2Edep(k) + EdepAbs(k)**2
0025         su2Trck(k) = su2Trck(k) + TrckAbs(k)**2
0026       enddo
0027 *
0028 * *** print
0029       IF (ISWIT(1).NE.0) THEN
0030         do k=1,NbAbsor
0031           CALL GEVKEV (EdepAbs(k),Eunit,unit)
0032           PRINT 751,k,Eunit,unit,TrckAbs(k)
0033         enddo
0034         print 749
0035       ENDIF
0036 *
0037 #ifndef batch
0038 *
0039 * *** drawing
0040 *
0041       IF (ISWIT(1).NE.0) THEN
0042          CALL GDHEAD (110110,'TestEm3',0.)
0043          CALL GDSHOW (3)
0044          CALL GDXYZ  (0)
0045       END IF
0046 #endif
0047   749 FORMAT(/)
0048   751 FORMAT(1X,'Absorber',I2,' : Energy deposit: ',F8.3,A4,
0049      &                        '   Track length: '  ,F8.3,' cm')
0050 *
0051       END