Made sure negative excitons warning is only for LXe not gas. -MMS

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

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;

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";