ugeom.F

Go to the documentation of this file.

      SUBROUTINE ugeom                                                          
*                                                                               
* *** Define user geometry set up                                               
*                                                                               
#include "geant321/gcbank.inc"
#include "pvolum.inc"
*                                                                                                                                       
      dimension zair(2),aair(2),wair(2)
      dimension zh2o(2),ah2o(2),wh2o(2)                                                  
      dimension zbgo(3),abgo(3),wbgo(3)                                         
      dimension zpbwo(3),apbwo(3),wpbwo(3)
      
      dimension par(3)      
                                                  
*                                                                               
* *** Air mixture parameters                                             
      DATA zair/  7.00,  8.00 /                  
      DATA aair/ 14.01, 16.00 /              
      DATA wair/  0.70,  0.30 /
*                                                                               
* *** H2O compound parameters                                             
      DATA zh2o/  1.00,  8.00 /                  
      DATA ah2o/  1.01, 16.00 /              
      DATA wh2o/  2.  ,    1. /                
*                                                                               
* *** BGO compound parameters 
      DATA zbgo/  8.00, 32.00,  83.00 /                                                         
      DATA abgo/ 16.00, 72.59, 208.98 /                                                                                          
      DATA wbgo/ 12.  ,  3.  , 4.     /                                                      
*                                                                               
* *** PbWO4 compound parameters
      DATA zpbwo/  8.00,  74.00,  82.00 /                                       
      DATA apbwo/ 16.00, 183.84, 207.19 /                                                                                                 
      DATA wpbwo/  4.  ,   1.  ,   1.   /                                          
*                                                                               
* *** Defines USER perticular materials
*
      CALL gsmixt( 1,'Air'      , aair , zair , 1.29e-3, 2,wair)
      CALL gsmixt( 2,'Water'    , ah2o , zh2o , 1.0    ,-2,wh2o)
      CALL gsmate( 3,'Ar Liquid', 40.00, 18.  , 1.39   ,14.0 ,84.0,0,0)
      CALL gsmate( 4,'Aluminium', 26.98, 13.  , 2.7    , 8.9 ,37.2,0,0)
      CALL gsmate( 5,'Iron'     , 55.85, 26.  , 7.87   , 1.76,17.1,0,0)
      CALL gsmixt( 6,'BGO'      , abgo , zbgo , 7.1    ,-3,wbgo)
      CALL gsmixt( 7,'PbWO4'    , apbwo, zpbwo, 8.28   ,-3,wpbwo)
      CALL gsmate( 8,'Lead '    ,207.19, 82.  ,11.35   ,0.56,18.5,0,0)
      CALL gsmate( 9,'Tungsten' ,183.85, 74.  ,19.30   ,0.35, 9.6,0,0)
      
      CALL gsmate(10,'Copper'   , 63.55, 29.  , 8.96   , 1.43,17.1,0,0)
      CALL gsmate(11,'Uranium'  ,238.03, 92.  ,18.95   , 0.32,17.1,0,0)
*                                                                               
* *** Defines USER tracking media parameters                                    
      fieldm = 0.0                                                             
      ifield = 0                                                               
      tmaxfd = 10.0                                                             
      stemax = 1000.
      IF (stepmax.gt.0.) stemax = stepmax                                                                  
      deemax = 0.20                                                            
      epsil  = 0.001                                                           
      stmin  = 0.010                                                           
*                                                                               
      CALL gstmed( 1,'Absorber',imat, 0 ,ifield,fieldm,tmaxfd,                
     *               stemax,deemax,epsil,stmin, 0 , 0 )        
*                                                                               
* ***  Defines USER'S VOLUMES                                                   
      jma = lq(jmate-imat)                                                      
      x0 = q(jma + 9)                                                           
      r1     =  nrtot*drx0*x0                                                   
      z1     =  nltot*dlx0*x0*0.5                                               
*                                                                               
      par(1) = 0.                                                               
      par(2) = r1                                                               
      par(3) = z1                                                               
      CALL gsvolu( 'ECAL' , 'TUBE' , 1 , par , 3 , ivol)                     
*                                                                               
* *** Close geometry banks. (obligatory system routine)                         
      CALL ggclos                                                               
**
* *** dessin
      CALL gsatt('*'   ,'SEEN',1)
      CALL gsatt('RING','SEEN',0)
*
      DO ix =1,3
         CALL gdopen(ix)
         scale =   9.5/z1
         paxis =   0.
         saxis =   0.2*z1
         CALL gdrawc('ECAL',ix,0.,10.,10.,scale,scale)
CCC         CALL GDAXIS (PAXIS,PAXIS,PAXIS,SAXIS)
         CALL gdscal( 10., 0.3)
         CALL gdclos
      END DO
*
                                                                               
      END