public void A_Z_Radionuclide(Radionuclide nuc, out int A, out int Z)
{
switch (nuc)
{
case Radionuclide.Ra_226:
A = 226;
Z = 88;
break;
case Radionuclide.Rn_222:
A = 222;
Z = 86;
break;
case Radionuclide.Po_218:
A = 218;
Z = 84;
break;
case Radionuclide.At_218:
A = 218;
Z = 85;
break;
case Radionuclide.Rn_218:
A = 218;
Z = 86;
break;
case Radionuclide.Pb_214:
A = 214;
Z = 82;
break;
case Radionuclide.Bi_214:
A = 214;
Z = 83;
break;
case Radionuclide.Po_214:
A = 214;
Z = 84;
break;
case Radionuclide.Tl_210:
A = 210;
Z = 81;
break;
case Radionuclide.Pb_210:
A = 210;
Z = 82;
break;
default:
A = 0;
Z = 0;
break;
}
}
public double halflife_Radionuclide(Radionuclide nuc)
{
double hl;
//FILL THIS !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
switch (nuc)
{
case Radionuclide.Ra_226:
hl = 1600;
break;
case Radionuclide.Rn_222:
hl = 3;
break;
case Radionuclide.Po_218:
hl = 0;
break;
case Radionuclide.At_218:
hl = 0;
break;
case Radionuclide.Rn_218:
hl = 0;
break;
case Radionuclide.Pb_214:
hl = 0;
break;
case Radionuclide.Bi_214:
hl = 0;
break;
case Radionuclide.Po_214:
hl = 0;
break;
case Radionuclide.Tl_210:
hl = 0;
break;
case Radionuclide.Pb_210:
hl = 0;
break;
default:
hl = 0;
break;
}
return(hl);
}
public Emission_type Emission_Radionuclide(Radionuclide nuc)
{
Emission_type em;
//FILL THIS !!!!!!!!!!!!!!!!!!!!!!!!!!!!!
switch (nuc)
{
case Radionuclide.Ra_226:
em = Emission_type.ALPHA;
break;
case Radionuclide.Rn_222:
em = Emission_type.ALPHA;
break;
case Radionuclide.Po_218:
em = Emission_type.ALPHA;
break;
case Radionuclide.At_218:
em = Emission_type.ALPHA;
break;
case Radionuclide.Rn_218:
em = Emission_type.ALPHA;
break;
case Radionuclide.Pb_214:
em = Emission_type.ALPHA;
break;
case Radionuclide.Bi_214:
em = Emission_type.ALPHA;
break;
case Radionuclide.Po_214:
em = Emission_type.ALPHA;
break;
case Radionuclide.Tl_210:
em = Emission_type.BETA;
break;
case Radionuclide.Pb_210:
em = Emission_type.ALPHA;
break;
default:
em = Emission_type.ALPHA;
break;
}
return(em);
}
public Isotope(Radionuclide nuc, double activ = 0, Radionuclide prog = Radionuclide.Un)
{
ln2 = Math.Log(2);
nuclide = nuc;
isotope_name = Isotope_by_name(nuc);
A_Z_Radionuclide(nuc, out A, out Z);
decay_type = Emission_Radionuclide(nuc);
half_life = halflife_Radionuclide(nuc);
activity = activ;
N = activity * half_life / ln2;
progeny = prog;
}
private string Isotope_by_name(Radionuclide name)
{
string isotope_name = "Unknown";
switch (name)
{
case Radionuclide.Ra_226:
isotope_name = "Radium-226";
break;
case Radionuclide.Rn_222:
isotope_name = "Radon-222";
break;
case Radionuclide.Po_218:
isotope_name = "Polunium-218";
break;
case Radionuclide.At_218:
isotope_name = "Astatium-222";
break;
case Radionuclide.Rn_218:
isotope_name = "Radon-218";
break;
case Radionuclide.Pb_214:
isotope_name = "Plumbum-214";
break;
case Radionuclide.Bi_214:
isotope_name = "Bismutum-214";
break;
case Radionuclide.Po_214:
isotope_name = "Polonium-214";
break;
case Radionuclide.Tl_210:
isotope_name = "Thallium-210";
break;
case Radionuclide.Pb_210:
isotope_name = "Plumbum-210";
break;
default:
break;
}
return(isotope_name);
}
public Isotope(string name, double half_l, int z = 0, int a = 0, double activ = 0, Emission_type type = Emission_type.ALPHA, double a_prob = 1.0, Radionuclide progen = Radionuclide.Un)
{
ln2 = Math.Log(2);
nuclide = Radionuclide.Un;
isotope_name = name;
half_life = half_l;
activity = activ;
N = activity * half_life / ln2; //Check this!
decay_type = type;
Z = z;
A = a;
//Advanced:
progeny = progen;
A_prob = a_prob;
B_prob = 1 - A_prob;
//Dose sim:
excited = false;
}
public Radionuclide Isotope_by_A_Z(int a, int z)
{
Radionuclide nuc = Radionuclide.Un;
switch (z)
{
case 88:
if (226 == a)
{
nuc = Radionuclide.Ra_226;
}
else
{
nuc = Radionuclide.Un;
}
break;
case 86:
if (222 == a)
{
nuc = Radionuclide.Rn_222;
}
else if (218 == a)
{
nuc = Radionuclide.Rn_218;
}
else
{
nuc = Radionuclide.Un;
}
break;
case 85:
if (218 == a)
{
nuc = Radionuclide.At_218;
}
else
{
nuc = Radionuclide.Un;
}
break;
case 84:
if (218 == a)
{
nuc = Radionuclide.Po_218;
}
else if (214 == a)
{
nuc = Radionuclide.Po_214;
}
else
{
nuc = Radionuclide.Un;
}
break;
case 82:
if (214 == a)
{
nuc = Radionuclide.Pb_214;
}
else if (210 == a)
{
nuc = Radionuclide.Pb_210;
}
else
{
nuc = Radionuclide.Un;
}
break;
case 81:
if (210 == a)
{
nuc = Radionuclide.Tl_210;
}
else
{
nuc = Radionuclide.Un;
}
break;
default:
break;
}
return(nuc);
}