Commit cf8ee507 authored by LEFEBVREJP email's avatar LEFEBVREJP email
Browse files

Updating radixams to solve nan result.

parent ba3d338c
Pipeline #7748 passed with stages
in 19 minutes and 24 seconds
...@@ -14,7 +14,7 @@ int main(int argc, char **argv) ...@@ -14,7 +14,7 @@ int main(int argc, char **argv)
{ {
CommandLine command(argc, argv); CommandLine command(argc, argv);
command.declareArgument("altitude" command.declareArgument("altitude"
, "Altitude at which to calcuate dose rate (0 for sea level)"); , "Altitude at which to calcuate dose rate (0 for sea level) in m");
command.declareArgument("activity" command.declareArgument("activity"
, "Activity of the area in Ci/m^2"); , "Activity of the area in Ci/m^2");
command.declareArgument("energy" command.declareArgument("energy"
...@@ -61,14 +61,14 @@ int main(int argc, char **argv) ...@@ -61,14 +61,14 @@ int main(int argc, char **argv)
auto mut = radix::gammaRayAttenuationInAir(e); auto mut = radix::gammaRayAttenuationInAir(e);
// mu scalars // mu scalars
auto mutg = -mut * 1.225; auto mutg = mut * 1.225;
// auto muag = -mua * 1.225; // auto muag = -mua * 1.225;
auto mutz = -mut * sigma; auto mutz = mut * sigma;
// auto muaz = -mua * sigma; // auto muaz = -mua * sigma;
// exponential integral values and their ratio // exponential integral values and their ratio
auto EIg = radix::exponentialIntegral(mutg); auto EIg = 1.0 / radix::exponentialIntegral(mutg);
auto EIz = radix::exponentialIntegral(mutz * z); auto EIz = 1.0 / radix::exponentialIntegral(mutz * z);
auto EIRatio = EIz / EIg; auto EIRatio = EIz / EIg;
// dose rate at sea level/altitude // dose rate at sea level/altitude
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment