ugimat.F

Go to the documentation of this file.

      SUBROUTINE ugimat                                                          
*                                                                               
* *** Define user materials                                               
*
#include "geant321/gcbank.inc"                                                                                                                                                                                                                       
*
      dimension aair(2),zair(2),wair(2)
      dimension aco2(2),zco2(2),wco2(2)
      dimension ah2o(2),zh2o(2),wh2o(2)
      dimension ag10(4),zg10(4),wg10(4)
      dimension asci(2),zsci(2),wsci(2)                                              
*      
* *** Air compound parameters          
      DATA aair/14.01, 16.00/
      DATA zair/ 7.  ,  8.  /
      DATA wair/ 0.7 ,  0.3 /
*                                                                               
* *** CO2 compound parameters
      DATA aco2/12.01, 16.00/
      DATA zco2/ 6.  ,  8.  /
      DATA wco2/ 1.  ,  2.  /
*                                                                               
* *** Water compound parameters
      DATA ah2o/ 1.01, 16.00/
      DATA zh2o/ 1.  ,  8.  /
      DATA wh2o/ 2.  ,  1.  /
*                                                                               
* *** G10 compound parameters
      DATA ag10/ 1.01, 12.00, 16.00, 28.00/
      DATA zg10/ 1.  ,  6.  ,  8.  , 14.  /
      DATA wg10/ 3.  ,  3.  ,  2.  ,  1.  /      
*                                                                               
* *** Scintillator compound parameters
      DATA asci/12.01,  1.01/
      DATA zsci/ 6.  ,  1.  /
      DATA wsci/ 9.  , 10.  /                                                                               

*
* *** Defines  materials
      CALL gsmixt( 1,'Air'      , aair ,zair, 1.29e-3, 2   , wair)
      CALL gsmixt( 2,'CO2 gas'  , aco2 ,zco2, 27.0e-3,-2   , wco2)
      CALL gsmate( 3,'H2 Liquid',  1.008,  1., 0.0708 , 865., 790., 0,0)
      CALL gsmixt( 4,'Water'    , ah2o ,zh2o, 1.0    ,-2   , wh2o)
      CALL gsmate( 5,'Liquid Ar', 39.95, 18., 1.39   , 14.0, 84.0, 0,0)
      CALL gsmate( 6,'Aluminium', 26.98, 13., 2.7    ,  8.9, 37.2, 0,0)
      CALL gsmate( 7,'Iron     ', 55.85, 26., 7.87   , 1.76, 17.1, 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,'Silicon  ', 28.09, 14.,  2.33  , 9.36, 45.5, 0,0)
      CALL gsmate(11,'Tungsten ',183.85, 74., 19.30  , 0.35,  9.6, 0,0)
      CALL gsmixt(12,'NemaG10'  , ag10 ,zg10, 1.7    ,-4   , wg10)
      CALL gsmate(13,'Copper   ', 63.55, 29., 8.96   , 1.43, 15.0, 0,0)
      CALL gsmixt(14,'Scintilla', asci ,zsci, 1.032  ,-2   , wsci)
      CALL gsmate(15,'Gold     ',196.97, 79., 19.32  , 0.33,  9.6, 0,0)
*
*     total number of materals
      mtot = 15            
* 
* *** overwrite the computed radlength of some mixture
      jma = lq(jmate-14)
      q(jma+9) = 42.549            
*                                               
*                                                                               
* *** Defines USER tracking media parameters
      ifield = 0                                                             
      fieldm = 0.                                                           
      tmaxfd = 10.0                                                             
      stemax = 1000.                                                            
      deemax = 0.20                                                            
      epsil  = 0.001                                                           
      stmin  = 0.010                                                           
*
      do im =1,mtot                                                                               
        CALL gstmed(im,'tracking medium',im, 0 ,ifield,fieldm,tmaxfd,                
     *                 stemax,deemax,epsil,stmin, 0 , 0 )
      enddo        
*
      END