------
A material
Example
-------
import dtsa2.jmGen as jmg
# Note: composition came from J. Donovan: Std 160 NBS K-412 mineral glass
l_comps = [epq.Element.O, epq.Element.Si, epq.Element.Mg, epq.Element.Ca,
epq.Element.Fe, epq.Element.Al, epq.Element.Mn, epq.Element.Na]
l_mfs = [0.43597 ,0.21199 ,0.11657 ,0.10899,
0.07742, 0.04906, 0.00077, 0.00043]
mat = jmg.create_material_from_mf(l_comps, l_mfs, 2.66, "K412")
print(mat)
print(mat.toHTMLTable())
"""
import dtsa2.jmGen as jmg
# Note: composition came from J. Donovan: Std 160 NBS K-412 mineral glass
l_comps = [
epq.Element.O, epq.Element.Si, epq.Element.Mg, epq.Element.Ca,
epq.Element.Fe, epq.Element.Al, epq.Element.Mn, epq.Element.Na
]
l_mfs = [
0.43597, 0.21199, 0.11657, 0.10899, 0.07742, 0.04906, 0.00077, 0.00043
]
mat = jmg.create_material_from_mf(l_comps, l_mfs, 2.66, "K412")
print(mat)
print(mat.toHTMLTable())
O 0.05390
== ========
Tot 1.0
N.B. John Donovan notes:
From NBS (NIST), Microanalysis glasses
2 x 2 x 20 mm rod
Traces by EPMA (NBS) by R. Marinenko
**Oxygen by difference**
"""
l_comps = [epq.Element.Al, epq.Element.Mg, epq.Element.P, epq.Element.O]
l_mfs = [0.06470, 0.06650, 0.32980, 0.5390]
k496 = jmg.create_material_from_mf(l_comps, l_mfs, 2.6, "K-496")
# K497
l_comps = [
epq.Element.Al, epq.Element.Mg, epq.Element.P, epq.Element.O,
epq.Element.Pb, epq.Element.Si, epq.Element.Ti, epq.Element.Fe,
epq.Element.Zr, epq.Element.Ce, epq.Element.Ta
]
l_mfs = [
0.045950, 0.055450, 0.211790, 0.680870, 0.000860, 0.000960, 0.000950,
0.000970, 0.000730, 0.001390, 0.000080
]
"""
Al 0.045950
Mg 0.055450
"""
l_comps = [
epq.Element.C, epq.Element.Mn, epq.Element.P, epq.Element.S,
epq.Element.Si, epq.Element.Cu, epq.Element.Ni, epq.Element.Cr,
epq.Element.V, epq.Element.Mo, epq.Element.Co, epq.Element.Pb,
epq.Element.Fe
]
l_mfs = [
0.046, 1.630, 0.020, 0.018, 0.502, 0.169, 12.180, 18.450, 0.047, 2.380,
0.10, 0.001, 64.457
]
srm1155 = jmg.create_material_from_mf(l_comps, l_mfs, 7.97, "SRM1155")
c = material("C", density=2.267)
mn = material("Mn", density=7.43)
p = material("P", density=1.88)
s = material("S", density=2.0)
si = material("Si", density=2.32)
cu = material("Cu", density=8.96)
mo = material("Mo", density=10.2)
co = material("Co", density=8.86)
pb = material("Pb", density=11.34)
fe = material("Fe", density=7.874)
xrts = []
trs = mc3.suggestTransitions(srm1155, e0)
Iron 26 0.1049 0.1049 0.0477566
Barium 56 0.1343 0.1343 0.0248635
-- 19.14 1 1.0 1.0
"""
l_comps = [
epq.Element.O, epq.Element.Al, epq.Element.Si, epq.Element.Ca,
epq.Element.Fe, epq.Element.Ti, epq.Element.Ba
]
# Note add a bit of Ti (impurity) in video not in list
l_mfs = [0.3872, 0.0794, 0.187, 0.1072, 0.1049, 0.0006, 0.1343]
k309 = jmg.create_material_from_mf(l_comps, l_mfs, 2.6, "K309")
si = material("Si", density=2.32)
caf2 = material("CaF2", density=3.18)
fe = material("Fe", density=7.874)
ti = material("Ti", density=4.506)
zn = material("Zn", density=7.14)
baf2 = material("BaF2", density=4.89)
al2o3 = material("Al2O3", density=3.95)
xrts = []
trs = mc3.suggestTransitions(k309, e0)
for tr in trs:
xrts.append(tr)
mat = jmg.create_material_from_mf(l_comps, l_mfs, 2.66, "K412")
print(mat)
print(mat.toHTMLTable())
K309
Density 3 g/cm³
Element Z Mass Frac Norm Mass Atom Frac
------- - --------- ----- ---------
Oxygen 8 0.3872 0.3872 0.615282
Aluminum 13 0.0794 0.0794 0.0748163
Silicon 14 0.187 0.187 0.169278
Calcium 20 0.1072 0.1072 0.0680035
Iron 26 0.1049 0.1049 0.0477566
Barium 56 0.1343 0.1343 0.0248635
"""
import dtsa2.jmGen as jmg
l_comps = [
epq.Element.O, epq.Element.Al, epq.Element.Si, epq.Element.Ca,
epq.Element.Fe, epq.Element.Ba
]
l_mfs = [0.3872, 0.0794, 0.1870, 0.1072, 0.1049, 0.1343]
mat = jmg.create_material_from_mf(l_comps, l_mfs, 3.0, "K309")
print(mat)
print(mat.toHTMLTable())
spcDir = outPath + "/sim-k978"
jmg.ensureDir(spcDir)
os.chdir(spcDir)
l_comps = [ epq.Element.Mg, epq.Element.Al, epq.Element.Si,
epq.Element.K, epq.Element.Ca, epq.Element.Ti,
epq.Element.Rb, epq.Element.Sr, epq.Element.Cs,
epq.Element.La, epq.Element.Eu, epq.Element.Gd,
epq.Element.O]
l_mfs = [ 0.013569, 0.089228, 0.21072,
0.018678, 0.067467, 0.006714,
0.020574, 0.100541, 0.002122,
0.019185, 0.019432, 0.019520,
0.393897]
k978 = jmg.create_material_from_mf(l_comps, l_mfs, 3.929, "K-978")
c = material("C", density=2.1)
layers = [ [c, 2.0e-9],
[k978, 50.0e-6]
]
k978Sim = mc3.multiFilm(layers, det, e0, withPoisson=False, nTraj=nTraj,
dose=dose, sf=True, bf=True, xtraParams={})
sName = "%g-nm-C-on-K978" % tNmC
k978Sim.rename(sName)
k978Sim.setAsStandard(k978)
k978Sim.display()
fi = spcDir + "/"
fi += sName
fi += "-%g-Traj.msa" % (nTraj)