mass number matching (idaholab#300)

snschune · Aug 9, 2018 · 9c1b144 · 9c1b144
1 parent 748e618
commit 9c1b144
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 2 deletions.
diff --git a/include/userobjects/PKAGeneratorBase.h b/include/userobjects/PKAGeneratorBase.h
@@ -47,6 +47,9 @@ class PKAGeneratorBase : public DiscreteElementUserObject
 
   /// finds the right ion tag; -1 means that the nuclide is not tracked, otherwise the index in the rasterizer nuclide vector must be retrieved
   int ionTag(const std::vector<Real> & rasterizer_Z, const std::vector<Real> & rasterizer_m, Real Z, Real m) const;
+
+  /// tolerance for mass number matching during tagging
+  const Real _mass_number_tolerance;
 };
 
 #endif // PKAGENERATORBASE_H
diff --git a/src/userobjects/PKAGeneratorBase.C b/src/userobjects/PKAGeneratorBase.C
@@ -14,12 +14,14 @@ template<>
 InputParameters validParams<PKAGeneratorBase>()
 {
   InputParameters params = validParams<DiscreteElementUserObject>();
+  params.addParam<Real>("mass_number_tolerance", 0.5, "Tolerance for matching mass numbers of recoils with rasterizer isotopes.");
   params.addClassDescription("PKA generator user object base class.\n Takes pdf and samples PKAs due to various interactions.");
   return params;
 }
 
 PKAGeneratorBase::PKAGeneratorBase(const InputParameters & parameters) :
-    DiscreteElementUserObject(parameters)
+    DiscreteElementUserObject(parameters),
+    _mass_number_tolerance(getParam<Real>("mass_number_tolerance"))
 {
   setRandomResetFrequency(EXEC_TIMESTEP_END);
 }
@@ -45,7 +47,7 @@ PKAGeneratorBase::ionTag(const std::vector<Real> & rasterizer_Z, const std::vect
   }
   else if (count > 1)
     for (unsigned int j = 0; j < rasterizer_Z.size(); ++j)
-      if (rasterizer_Z[j] == Z && rasterizer_m[j] == m)
+      if (rasterizer_Z[j] == Z && std::abs(rasterizer_m[j] - m) < _mass_number_tolerance)
         return j;
   return -1;
 }