void KVRangeYanez::AddElementalMaterial(Int_t z, Int_t a)
{
// Adds a material composed of a single isotope of a chemical element.
// If the isotope (a) is not specified, we create a material containing the naturally
// occuring isotopes of the given element, weighted according to their abundance.
// If the mass is given, the material symbol will be "AX" where X is the symbol for the element
// e.g. "48Ca", "124Sn", etc.
// and the material name will be "Xxx-A" where Xxx is the name of the element
// e.g. "Calcium-48", "Tin-124", etc.
// Otherwise, we just use the element symbol and name for naturally-occurring
// mixtures of atomic elements ("Ca", "Calcium", etc.).
KVIonRangeTableMaterial* mat;
if (!a) {
mat = MakeNaturallyOccuringElementMixture(z);
} else {
if (!gNDTManager) {
Error("AddElementalMaterial",
"Nuclear data tables have not been initialised");
return;
}
KVElementDensity* ed = (KVElementDensity*)gNDTManager->GetData(z, a, "ElementDensity");
if (!ed) {
Error("AddElementalMaterial",
"No element found in ElementDensity NDT-table with Z=%d", z);
return;
}
TString state = "solid";
if (ed->IsGas()) state = "gas";
mat = new KVRangeYanezMaterial(this, Form("%s-%d", ed->GetElementName(), a),
Form("%d%s", a, ed->GetElementSymbol()),
state, ed->GetValue(), z, a);
mat->Initialize();
}
CheckMaterialsList();
fMaterials->Add(mat);
mat->ls();
}
KVIonRangeTableMaterial::KVIonRangeTableMaterial(const KVIonRangeTable* tab, const Char_t* name, const Char_t* symbol,
const Char_t* state, Double_t density, Double_t Z, Double_t A)
: KVBase(name, symbol),
fTable(tab),
fState(state),
fComposition(0),
fCompound(kFALSE),
fMixture(kFALSE),
fDens(density),
fZmat(Z),
fAmat(A),
fMoleWt(0),
fDeltaE(0),
fEres(0),
fRange(0),
fStopping(0)
{
// Create new material with given (long) name and symbol
// symbol convention: for elements, use element symbol. for compounds, use chemical formula.
// e.g. silicon: name="Silicon" symbol="Si"
// C4H10: name="Isobutane", symbol="C4H10"
// state="solid", "liquid", "gas"
//
// Density [g/cm**3] must be given for solid compounds/mixtures.
// Densities of atomic elements are known (gNDTManager->GetValue(z, a, "ElementDensity")), they
// will be used automatically, unless a different value is given here.
// Densities of gases are calculated from the molar weight, temperature and pressure.
if (Z > 0 && density < 0) {
KVElementDensity* ed = (KVElementDensity*)gNDTManager->GetData(Z, A, "ElementDensity");
if (!ed) {
Warning("KVIonRangeTableMaterial",
"No element found in density table with Z=%f, density unknown", Z);
}
else
fDens = ed->GetValue();
}
}
KVIonRangeTableMaterial* KVRangeYanez::MakeNaturallyOccuringElementMixture(Int_t z)
{
// create a material containing the naturally occuring isotopes of the given element,
// weighted according to their abundance.
if (!gNDTManager) {
Error("MakeNaturallyOccuringElementMixture",
"Nuclear data tables have not been initialised");
return NULL;
}
KVElementDensity* ed = (KVElementDensity*)gNDTManager->GetData(z, z, "ElementDensity");
if (!ed) {
Error("AddElementalMaterial",
"No element found in ElementDensity NDT-table with Z=%d", z);
return 0x0;
}
TString state = "solid";
if (ed->IsGas()) state = "gas";
KVRangeYanezMaterial* mat =
new KVRangeYanezMaterial(this,
ed->GetElementName(),
ed->GetElementSymbol(),
state, ed->GetValue());
KVNucleus nuc(z);
KVNumberList isotopes = nuc.GetKnownARange();
isotopes.Begin();
while (!isotopes.End()) {
nuc.SetA(isotopes.Next());
Double_t abundance = nuc.GetAbundance() / 100.;
if (abundance > 0.) mat->AddMixtureElement(z, nuc.GetA(), 1, abundance);
}
mat->Initialize();
return (KVIonRangeTableMaterial*)mat;
}