From f4d0da0149b5b2322da78c4a540e3737502bc35c Mon Sep 17 00:00:00 2001
From: mszydagis <mszydagis@albany.edu>
Date: Sun, 8 May 2022 20:42:28 -0400
Subject: [PATCH] Made sure negative excitons warning is only for LXe not gas.
 -MMS

---
 examples/loopNEST.in | 2 +-
 src/NEST.cpp         | 4 ++--
 src/execNEST.cpp     | 2 +-
 3 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/examples/loopNEST.in b/examples/loopNEST.in
index a418b40e..f9d482c5 100644
--- a/examples/loopNEST.in
+++ b/examples/loopNEST.in
@@ -7,7 +7,7 @@
 #Search the word "loopNEST" in src/execNEST.cpp for more info on what each of these input arguments is doing
 
 #First declare ranges for each of the free parameters
-#Here, in several for-loops, using 11 parameters with set ranges
+#Here, in several for-loops, using 11 parameters with set ranges. 3H Run3 best: .00426 43.4 -.524 .775 74 .00989 2.95 .0967 .541 .648 .949. 2e5 skew S1 Hybrid
 for a in 0.001 0.002 0.005 #step size in mm
 do
 for b in 30 40 50 #base for dE/dx power law at high energies
diff --git a/src/NEST.cpp b/src/NEST.cpp
index fcb57ef9..62af3f9d 100644
--- a/src/NEST.cpp
+++ b/src/NEST.cpp
@@ -769,8 +769,8 @@ YieldResult NESTcalc::GetYieldNR(
     Ni = Nq / (1. + NexONi);
     Nex = Nq - Ni;
   }
-
-  if (Nex < 0.)
+  
+  if ( Nex < 0. && density >= 1. )
     cerr << "\nCAUTION: You are approaching the border of NEST's validity for "
             "high-energy (OR, for LOW) NR, or are beyond it, at "
          << energy << " keV." << endl;
diff --git a/src/execNEST.cpp b/src/execNEST.cpp
index a5f4ec4f..d99eb71e 100644
--- a/src/execNEST.cpp
+++ b/src/execNEST.cpp
@@ -105,7 +105,7 @@ int main(int argc, char** argv) {
     type = "MIP";
     dEOdxBasis = true;
     eMin = -TestSpectra::CH3T_spectrum(0.,18.6);
-    //s1CalculationMode = NEST::S1CalculationMode::Parametric;
+    s1CalculationMode = NEST::S1CalculationMode::Parametric;
     eMax = -1.;
     inField = 180.;
     position = "-1";