def pure(cls, z, absorption_energy_eV=None):
"""
Returns the material for the specified pure element.
:arg z: atomic number
:type z: :class:`int`
:arg absorption_energy_eV: absorption energy of the different particles
:type absorption_energy_eV: :class:`dict`
"""
name = ep.name(z)
composition = {z: 1.0}
density_kg_m3 = ep.mass_density_kg_m3(z)
return cls(composition, name, density_kg_m3, absorption_energy_eV)
def kanaya_okayama(composition, energy):
"""
Returns the electron range (in meters).
:arg composition: composition in weight fraction.
The composition is specified by a dictionary.
The keys are the atomic numbers and the values are the weight fractions
between ]0.0, 1.0].
:type composition: :class:`dict`
:arg energy: beam energy in eV
"""
r = 0.0;
for z, fraction in composition.items():
dr = (0.0276 * (ep.atomic_mass_kg_mol(z) * 1000.0) * (energy / 1000.0) ** 1.67) / \
(z ** 0.89 * (ep.mass_density_kg_m3(z) / 1000.0))
r += fraction / (dr * 1e-6)
return 1.0 / r;
def calculate_density(composition):
"""
Returns an estimate density from the composition using the pure element
density and this equation.
.. math::
\\frac{1}{\\rho} = \\sum{\\frac{1}{\\rho_i}}
:arg composition: composition in weight fraction.
The composition is specified by a dictionary.
The key are atomic numbers and the values weight fractions.
No wildcard are accepted.
:type composition: :class:`dict`
"""
density = 0.0
for z, fraction in composition.items():
density += fraction / ep.mass_density_kg_m3(z)
return 1.0 / density
def testmass_density(self):
self.assertAlmostEqual(7.1900, ep.mass_density_kg_m3(24) / 1000.0)
def testmass_density(self):
self.assertAlmostEqual(7.1900, ep.mass_density_kg_m3(24) / 1000.0)
self.assertRaises(ValueError, self.ep.mass_density_kg_m3, 85)
self.assertRaises(ValueError, self.ep.mass_density_kg_m3, 87)
self.assertRaises(ValueError, self.ep.mass_density_kg_m3, 118)
def testmass_density(self):
self.assertAlmostEqual(7.1900, ep.mass_density_kg_m3(24) / 1000.0)
def testmass_density(self):
self.assertAlmostEqual(7.1900, ep.mass_density_kg_m3(24) / 1000.0)
self.assertRaises(ValueError, self.ep.mass_density_kg_m3, 85)
self.assertRaises(ValueError, self.ep.mass_density_kg_m3, 87)
self.assertRaises(ValueError, self.ep.mass_density_kg_m3, 118)
