ugeom.F

Go to the documentation of this file.

      SUBROUTINE ugeom                                                          
*                                                                               
* *** Define user geometry set up                                               
*                                                                                                                                                                                                                       
*
#include "detector.inc"  
*                                                                               
      dimension par(3)
      
      dimension aair(2),zair(2),wair(2)
      dimension ah2o(2),zh2o(2),wh2o(2)                                        
*                                                                               
* *** Air compound parameters                                                   
      DATA aair/14.01, 16.00/                                            
      DATA zair/ 7.  ,  8.  /                                                     
      DATA wair/ 0.7 ,  0.3 /
*                                                                               
* *** Air compound parameters                                                   
      DATA ah2o/ 1.01, 16.00/                                            
      DATA zh2o/ 1.  ,  8.  /                                                     
      DATA wh2o/ 2.  ,  1.  /                                                           
*
* *** Defines USER perticular materials
      CALL gsmixt( 1,'Air'      , aair ,zair, 1.29e-3, 2   , wair)                                                
      CALL gsmate( 2,'H2 Liquid',  1.01,  1., 0.0708 , 865., 790., 0,0)
      CALL gsmixt( 3,'Water'    , ah2o ,zh2o, 1.0    ,-2   , wh2o)      
      CALL gsmate( 4,'Liquid Ar', 39.95, 18., 1.39   , 14.0, 84.0, 0,0)             
      CALL gsmate( 5,'Aluminium', 26.98, 13., 2.7    ,  8.9, 37.2, 0,0)
      CALL gsmate( 6,'Iron'     , 55.85, 26., 7.87   , 1.76, 17.1, 0,0)
      CALL gsmate( 7,'Tungsten' ,183.85, 74., 19.30  , 0.35, 18.5, 0,0)                    
      CALL gsmate( 8,'Lead'     ,207.19, 82., 11.35  , 0.56, 18.5, 0,0)                    
      CALL gsmate( 9,'Uranium'  ,238.03, 92., 18.95  , 0.32, 12. , 0,0)                    
      CALL gsmate(10,'Germanium', 72.61, 32., 5.323  , 2.30, 16.6, 0,0)                                          
*                                                                               
* *** Defines USER tracking media parameters
      ifield = 0                                                             
      fieldm = 0.                                                               
      tmaxfd = 10.0                                                             
      stemax = 1000.                                                            
      deemax = 1.                                                            
      epsil  = 0.001                                                           
      stmin  = 0.010                                                           
*                                                                               
      CALL gstmed( 1,'Container',imat, 0 ,ifield,fieldm,tmaxfd,                
     *               stemax,deemax,epsil,stmin, 0 , 0 )        
*
*
* *** Geometry                                                                               
      par(1) = boxsize/2.                                                               
      par(2) = boxsize/2.                                                               
      par(3) = boxsize/2.                                                               
      CALL gsvolu('aBox','BOX ',1,par,3,ivol)                     
*                                                                               
* *** Close geometry banks. (obligatory system routine)                         
      CALL ggclos                                                               
*
*
* *** dessin
      CALL gsatt('*','SEEN',1)
*
      DO ix = 1,3
        CALL gdopen(ix)
        scale =   18./boxsize
        paxis =   0.
        saxis =   0.1*boxsize
        CALL gdrawc('aBox',ix,0.,10.,9.3,scale,scale)
        CALL gdaxis(paxis,paxis,paxis,saxis)
        CALL gdscal(10. , 0.3)
        CALL gdclos
      END DO
*
                                                                               
      END