diff --git a/wremnants/combine_theory_helper.py b/wremnants/combine_theory_helper.py
index 7d3c40e45..ff741072a 100644
--- a/wremnants/combine_theory_helper.py
+++ b/wremnants/combine_theory_helper.py
@@ -379,11 +379,6 @@ def add_gamma_np_uncertainties(self):
         var_vals = gamma_vals
         var_names = [f"{gamma_nuisance_name}Down", f"{gamma_nuisance_name}Up"]
 
-        if "Lambda" in self.np_model:
-            Lambda4_nuisance_name = "scetlibNPLambda4"
-            var_vals.extend(["Lambda4.01", "Lambda4.16"])
-            var_names.extend([f"{Lambda4_nuisance_name}Down", f"{Lambda4_nuisance_name}Up"])
-
         logger.debug(f"Adding gamma uncertainties from syst entries {gamma_vals}")
 
 
@@ -437,7 +432,8 @@ def add_resum_scale_uncertainty():
     def add_uncorrelated_np_uncertainties(self):
         np_map = {
             "Lambda2" : ["-0.25", "0.25",],
-            "Delta_Lambda2" : ["-0.02", "0.02",]
+            "Delta_Lambda2" : ["-0.02", "0.02",],
+            "Lambda4" : [".01", ".16"],
         } if "Lambda" in self.np_model else {
             "Omega" : ["0.", "0.8"],
             "Delta_Omega" : ["-0.02", "0.02"],
diff --git a/wremnants/theory_tools.py b/wremnants/theory_tools.py
index ebd63a626..17e31ba5b 100644
--- a/wremnants/theory_tools.py
+++ b/wremnants/theory_tools.py
@@ -624,7 +624,8 @@ def is_flavor_dependent_np(var_label):
             return var_label.startswith("Omega") \
                     or var_label.startswith("Delta_Omega") \
                     or var_label.startswith("Lambda2") \
-                    or var_label.startswith("Delta_Lambda2")
+                    or var_label.startswith("Delta_Lambda2") \
+                    or var_label.startswith("Lambda4")
 
         # special treatment for Lambda2/Omega since they need to be decorrelated in charge and possibly rapidity
         if isinstance(var_axis, hist.axis.StrCategory) and any(is_flavor_dependent_np(var_label) for var_label in var_axis):