gutrev.F

Go to the documentation of this file.

      SUBROUTINE gutrev                                                         
*                                                                               
*     User routine to control tracking of one event                 
*     Called by GRUN                                                            
*                                                                     
#include "celoss.inc"
#include "histo.inc"                                                                       
*
* *** par event initialisations                                                                          
      eevent  = 0.
*
* *** process the event                                                                               
      CALL gtreve
*
* *** returns if no energy deposition
      if (eevent.le.0.) return      
*
* *** histo: total energy deposited
      ih = 2
      if (histo(ih)) call hfill(ih,eevent/histunit(ih),0.,1.)
*
* *** cumul energy deposited
      etot  = etot  + eevent
      etot2 = etot2 + eevent*eevent
*                                                                                           
      END