22.05.2014 Views

CERN Program Library Long Writeup W5013 - CERNLIB ...

CERN Program Library Long Writeup W5013 - CERNLIB ...

CERN Program Library Long Writeup W5013 - CERNLIB ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Geant 3.16 GEANT User’s Guide CONS110<br />

Origin : R.Brun, M.Maire, J.Allison Submitted: 01.06.83<br />

Revision : Revised: 08.12.93<br />

Documentation : M.Maire<br />

Mixtures and Compounds<br />

CALL GSMIXT<br />

(IMATE,NAMATE,A,Z,DENS,NLMAT,WMAT*)<br />

Defines mixture or compound IMATE as composed by NLMAT materials defined via the arrays A, Z and WMAT.<br />

Mixtures of compounds can also be defined.<br />

IMATE (INTEGER) user material (mixture) number;<br />

NAMATE (CHARACTER*20) mixture name;<br />

A<br />

(REAL) array of atomic weights;<br />

Z<br />

(REAL) array of atomic numbers;<br />

DENS (REAL) density in g cm −3 ;<br />

NLMAT (INTEGER) number of elements in the mixture;<br />

> 0 WMAT contains the proportion by weights of each material in the mixture;<br />

< 0 WMAT contains the proportion by number of atoms of each kind, the content of WMAT<br />

in output is changed to contain the relative weights;<br />

Method<br />

For a compound (NLMAT < 0), the molecular weight and charge are:<br />

A mol = ∑ i<br />

n i A i<br />

Z mol = ∑ i<br />

n i Z i<br />

where n i = WMAT(I) is the number of atoms of the i th component of the molecule. In this case the proportion<br />

by weight is:<br />

p i = n iA i<br />

A mol<br />

where p i = WMAT(I) in output. From the relative weights, GSMIXT works out the effective atomic weight<br />

and atomic number according to the following formulas:<br />

A eff = ∑ i<br />

p i A i<br />

Z eff = ∑ i<br />

p i Z i<br />

which are stored in the JMATE data structure [CONS199] together with the radiation length.<br />

The radiation length is computed according to the EGS manual [12], formula 268 (37), for an element:<br />

[<br />

1<br />

=4αr0 2 ρX N 1<br />

Av<br />

0 A Z(Z + ξ(Z)) ln 183 ]<br />

Z 1/3 − F c(Z)<br />

where<br />

X 0 radiation length (in cm);<br />

ρ density (in g cm −3 );<br />

51 CONS110 – 1

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!