Changed the Pb-206 charge loss to be default 1.0 to not double-count …

…it again: problem pointed out by Nishat. -MMS
NESTCollaboration · Nov 11, 2021 · 49d50a8 · 49d50a8
1 parent 229aebd
commit 49d50a8
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/NEST.cpp b/src/NEST.cpp
@@ -9,7 +9,7 @@
         // 0.03)
 //#define InfraredNR 7.
 
-#define ChargeLoss 0.77  // LUX Run03 is the default (wall BG model)
+#define ChargeLoss 1.000  // 0.77 is LUX Run03 default (wall BG model)
 
 using namespace std;
 using namespace NEST;
@@ -649,7 +649,7 @@ YieldResult NESTcalc::GetYieldIon(
   double NexONi =
       alpha +
       0.00178 *
-          pow(atomNum, 1.587);  // Wq_eV=13.7; ThomasImel=0.05; NexONi=0.05;
+          pow(atomNum, 1.587);
   if (ValidityTests::nearlyEqual(A1, 206.) &&
       ValidityTests::nearlyEqual(Z1,
                                  82.)) {  // Pb-206 (from Po-210 alpha decay).