diff --git a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py
index c245b1a..8ee85aa 100644
--- a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py
+++ b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py
@@ -102,7 +102,6 @@
 
 ######################################################################################################
 
-
 ########## Tree for Lc* -> Lc pi pi ##################################################################
 Lc_star = CombineParticles("Lc_star")
 Lc_star.Inputs = ["Phys/StdAllLoosePions/Particles", "Phys/Lc2PKPiforB2DMuNuX/Particles"]
@@ -165,8 +164,12 @@
 
 
 ########## Tree for Lc* -> Lc pi0 ##################################################################
+pi0containers = ['StdLooseMergedPi0', 'StdLoosePi024e', 'StdLoosePi02gee', 'StdLoosePi02gg']
+
 Lc_star1pi0 = CombineParticles("Lc_star1pi0")
 Lc_star1pi0.Inputs = ["Phys/StdLooseMergedPi0/Particles", "Phys/Lc2PKPiforB2DMuNuX/Particles"]
+for pc in pi0containers:
+    Lc_star1pi0.Inputs.append('Phys/%s/Particles' % pc)
 Lc_star1pi0.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ pi0]cc"]
 Lc_star1pi0.MotherCut = "(M<3000) & (M>2300) & (VFASPF(VCHI2/VDOF) < 15.0) "
 Lc_star1pi0.DaughtersCuts = { "pi0": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
@@ -195,7 +198,9 @@
 
 ########## Tree for Lc* -> Lc pi0 pi pi ##################################################################
 Lc_starpi0pipi = CombineParticles("Lc_starpi0pipi")
-Lc_starpi0pipi.Inputs = ["Phys/StdLooseMergedPi0/Particles", "Phys/Lc2PKPiforB2DMuNuX/Particles", "Phys/StdAllLoosePions/Particles"]
+Lc_starpi0pipi.Inputs = ["Phys/Lc2PKPiforB2DMuNuX/Particles", "Phys/StdAllLoosePions/Particles"]
+for pc in pi0containers:
+    Lc_starpi0pipi.Inputs.append('Phys/%s/Particles' % pc)
 Lc_starpi0pipi.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ pi0 pi+ pi-]cc"]
 Lc_starpi0pipi.MotherCut = "(M<3000) & (M>2300) & (VFASPF(VCHI2/VDOF) < 15.0) "
 Lc_starpi0pipi.DaughtersCuts = { "pi0": "PT>0", "pi+": "PT>0", "pi-": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }