guout.F

Go to the documentation of this file.

      SUBROUTINE guout                                                          
*                                                                               
*     User routine called at the end of each event                              
*                                                                               
#include "geant321/gcflag.inc"
#include "geant321/gckine.inc"
#include "pvolum.inc"
#include "celoss.inc"     
*
* *** drawing
*
#ifndef batch                                                          
      IF (iswit(1).NE.0) THEN
         CALL gdhead(110110,'testem2',0.)
         CALL gdshow(1)
         CALL gdxyz(0)
      END IF
#endif      
*
*
* *** statistic
*                                                                              
      dlc = 0.                                                                  
      drc = 0.                                                                  

*     longitudinal profile
*                                                                                
      DO 2 i = 1,nltot
         sel1(i) = sel1(i) + dedl(i)                                          
         sel2(i) = sel2(i) + dedl(i)**2                                       
         dlc      = dlc      + dedl(i)                                          
         sel1c(i) = sel1c(i) + dlc                                              
         sel2c(i) = sel2c(i) + dlc**2
         
     bin = (float(i)-0.5)*dlx0
         CALL hfill(4,bin,100*dedl(i)/(dlx0*pkine(3)),1.)
         bin = float(i)*dlx0
         CALL hfill(6,bin,100*dlc    /pkine(3),1.)                                                                                                                 
   2  CONTINUE                                                                  

*     radial profile
*                                                                                     
      DO 3 i = 1,nrtot                                                                     
         ser1(i) = ser1(i) + dedr(i)                                          
         ser2(i) = ser2(i) + dedr(i)**2                                       
         drc      = drc      + dedr(i)                                          
         ser1c(i) = ser1c(i) + drc                                              
         ser2c(i) = ser2c(i) + drc**2
     
         bin = (float(i)-0.5)*drx0
         CALL hfill(8,bin,100*dedr(i)/(drx0*pkine(3)),1.)
         bin = float(i)*drx0
         CALL hfill(10,bin,100*drc    /pkine(3),1.)                                                  
   3  CONTINUE                                                                  
*
*     energy deposited and track length
*                                                                                
      eseen = 100.*dlc/pkine(3)                                                 
      CALL hfill(1, eseen,0.,1.)                                                
      CALL hfill(2,strch/x0,0.,1.)                                            
      CALL hfill(3,strne/x0,0.,1.)
*
      strch1 = strch1 + strch
      strch2 = strch2 + strch**2
      strne1 = strne1 + strne
      strne2 = strne2 + strne**2              
*                                                                              
      END