diff --git a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py
index 349915b..ef2de58 100644
--- a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py
+++ b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py
@@ -136,7 +136,7 @@
 omega_2pi.Inputs = ["Phys/StdAllLoosePions/Particles"]
 omega_2pi.DecayDescriptors = ["[omega(782) -> pi+ pi-]cc"]
 omega_2pi.MotherCut = " (M<1214) & (M>294) & (VFASPF(VCHI2/VDOF) < 12.0) "
-omega_2pi.DaughtersCuts = { "pi+": "PT>0", "pi-": "PT>0" }
+omega_2pi.DaughtersCuts = { "pi+": "(PT>0) & (PROBNNghost<0.3)", "pi-": "(PT>0) & (PROBNNghost<0.3)" }
 
 Lc_star4pi = CombineParticles("Lc_star4pi")
 Lc_star4pi.Inputs = ["Phys/omega_2pi/Particles", "Phys/Lc2PKPiforB2DMuNuX/Particles"]
@@ -228,7 +228,7 @@
     Lc_starpi0pipi.Inputs.append('Phys/%s/Particles' % pc)
 Lc_starpi0pipi.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ pi0 pi+ pi-]cc"]
 Lc_starpi0pipi.MotherCut = "(M<3500) & (M>2300) & (VFASPF(VCHI2/VDOF) < 7.0) "
-Lc_starpi0pipi.DaughtersCuts = { "pi0": "PT>0", "pi+": "PT>0", "pi-": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
+Lc_starpi0pipi.DaughtersCuts = { "pi0": "PT>0", "pi+": "(PT>0) & (PROBNNghost<0.3)", "pi-": "(PT>0) & (PROBNNghost<0.3)", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
 
 Lb_SL_Lcpi0pipi = CombineParticles("Lb_SL_Lcpi0pipi")
 Lb_SL_Lcpi0pipi.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/Lc_starpi0pipi/Particles"]
@@ -288,7 +288,7 @@
 Lc_star3pi.Inputs = ["Phys/StdAllLoosePions/Particles", "Phys/Lc2PKPiforB2DMuNuX/Particles"]
 Lc_star3pi.DecayDescriptors = ["[Sigma_c+ -> Lambda_c+ pi+ pi+ pi-]cc"]
 Lc_star3pi.MotherCut = "(M<3500) & (M>2300) & (VFASPF(VCHI2/VDOF) < 7.0) "
-Lc_star3pi.DaughtersCuts = { "pi+": "PT>0", "pi+": "PT>0", "pi-": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
+Lc_star3pi.DaughtersCuts = { "pi+": "(PT>0)  & (PROBNNghost<0.3)", "pi+": "(PT>0)  & (PROBNNghost<0.3)", "pi-": "(PT>0)  & (PROBNNghost<0.3)", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
 
 Lb_SL_Sc3pi = CombineParticles("Lb_SL_Sc3pi")
 Lb_SL_Sc3pi.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/Lc_star3pi/Particles"]