diff --git a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py
index 9eeca5b..d3274ec 100644
--- a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py
+++ b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_moreStates.py
@@ -103,16 +103,20 @@
 ######################################################################################################
 
 ########## Tree for Lc* -> Lc pi pi ##################################################################
+LcStarMassCut = "(M>2543) & (M<2675)"
+
 Lc_star = CombineParticles("Lc_star")
 Lc_star.Inputs = ["Phys/StdAllLoosePions/Particles", "Phys/Lc2PKPiforB2DMuNuX/Particles"]
 Lc_star.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ pi+ pi-]cc"]
-Lc_star.MotherCut = "(M<3500) & (M>2300) & (VFASPF(VCHI2/VDOF) < 7.0) "
+Lc_star.MotherCut = LcStarMassCut + " & (VFASPF(VCHI2/VDOF) < 7.0) "
 Lc_star.DaughtersCuts = { "pi+": "PT>0", "pi-": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
 
+LbMassCut = "(M<6500)"
+
 Lb_SL = CombineParticles("Lb_SL")
 Lb_SL.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/Lc_star/Particles"]
 Lb_SL.DecayDescriptors = [ "[Lambda_b0 -> Lambda_c(2595)+ mu-]cc" ]
-Lb_SL.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
+Lb_SL.MotherCut = LbMassCut + " & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
 
 decay0 = "[Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^pi+ ^pi-) ^mu- ]CC"
 
@@ -131,50 +135,50 @@
     "K_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) pi+ pi-) mu- ]CC"
 }
 
-########## Tree for Lc** -> (Lc* -> Lc pi pi) pi pi ##################################################
-omega_2pi = CombineParticles("omega_2pi")
-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) & (PROBNNghost<0.3)", "pi-": "(PT>0) & (PROBNNghost<0.3)" }
-omega_2pi.CombinationCut = "ACUTDOCA(1.0*mm,'')"
+# ########## Tree for Lc** -> (Lc* -> Lc pi pi) pi pi ##################################################
+# omega_2pi = CombineParticles("omega_2pi")
+# 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) & (PROBNNghost<0.3)", "pi-": "(PT>0) & (PROBNNghost<0.3)" }
+# omega_2pi.CombinationCut = "ACUTDOCA(1.0*mm,'')"
 
-Lc_star4pi = CombineParticles("Lc_star4pi")
-Lc_star4pi.Inputs = ["Phys/omega_2pi/Particles", "Phys/Lc2PKPiforB2DMuNuX/Particles"]
-## Too much combinatorics this way
-# Lc_star4pi.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ pi+ pi- pi+ pi-]cc"]
-Lc_star4pi.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ omega(782) omega(782)]cc"]
-Lc_star4pi.MotherCut = "(M<3500) & (M>2300) & (VFASPF(VCHI2/VDOF) < 7.0) "
-#Lc_star4pi.DaughtersCuts = { "pi+": "PT>0", "pi-": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
-Lc_star4pi.DaughtersCuts = { "omega(782)": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
-Lc_star4pi.CombinationCut = "ACUTDOCA(1.0*mm,'')"
+# Lc_star4pi = CombineParticles("Lc_star4pi")
+# Lc_star4pi.Inputs = ["Phys/omega_2pi/Particles", "Phys/Lc2PKPiforB2DMuNuX/Particles"]
+# ## Too much combinatorics this way
+# # Lc_star4pi.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ pi+ pi- pi+ pi-]cc"]
+# Lc_star4pi.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ omega(782) omega(782)]cc"]
+# Lc_star4pi.MotherCut = "(M<3500) & (M>2300) & (VFASPF(VCHI2/VDOF) < 7.0) "
+# #Lc_star4pi.DaughtersCuts = { "pi+": "PT>0", "pi-": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
+# Lc_star4pi.DaughtersCuts = { "omega(782)": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
+# Lc_star4pi.CombinationCut = "ACUTDOCA(1.0*mm,'')"
 
-Lb_SL_Lc4pi = CombineParticles("Lb_SL_Lc4pi")
-Lb_SL_Lc4pi.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/Lc_star4pi/Particles"]
-Lb_SL_Lc4pi.DecayDescriptors = [ "[Lambda_b0 -> Lambda_c(2595)+ mu-]cc" ]
-Lb_SL_Lc4pi.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
+# Lb_SL_Lc4pi = CombineParticles("Lb_SL_Lc4pi")
+# Lb_SL_Lc4pi.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/Lc_star4pi/Particles"]
+# Lb_SL_Lc4pi.DecayDescriptors = [ "[Lambda_b0 -> Lambda_c(2595)+ mu-]cc" ]
+# Lb_SL_Lc4pi.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
 
-#decay0_Lc4pi = "[Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^pi+ ^pi- ^pi+ ^pi-) ^mu- ]CC"
-decay0_Lc4pi = "[Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^(omega(782) -> ^pi+ ^pi-) ^(omega(782) -> ^pi+ ^pi-)) ^mu- ]CC"
+# #decay0_Lc4pi = "[Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^pi+ ^pi- ^pi+ ^pi-) ^mu- ]CC"
+# decay0_Lc4pi = "[Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^(omega(782) -> ^pi+ ^pi-) ^(omega(782) -> ^pi+ ^pi-)) ^mu- ]CC"
 
-Lb_SL_Lc4pi_OfflineTree = DecayTreeTuple("Lb_SL_Lc4pi_OfflineTree")
-Lb_SL_Lc4pi_OfflineTree.Inputs = ["Phys/Lb_SL_Lc4pi/Particles"]
-Lb_SL_Lc4pi_OfflineTree.Decay = decay0_Lc4pi
-Lb_SL_Lc4pi_OfflineTree.Branches = {
-    "Lambda_b0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
-    "Lambda_c_STAR": "[Lambda_b0 -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
-    "mu": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) ^mu- ]CC",
-    "Lambda_c": "[Lambda_b0 -> (Lambda_c(2595)+ -> ^(Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
-    "pi_plus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> ^pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
-    "pi_minus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ ^pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
-    "pi_plus2": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> ^pi+ pi-)) mu- ]CC",
-    "pi_minus2": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ ^pi-)) mu- ]CC",
-    "omega": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) ^(omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
-    "omega2": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) ^(omega(782) -> pi+ pi-)) mu- ]CC",
-    "p": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- ^p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
-    "pi_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ ^pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
-    "K_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC"
-}
+# Lb_SL_Lc4pi_OfflineTree = DecayTreeTuple("Lb_SL_Lc4pi_OfflineTree")
+# Lb_SL_Lc4pi_OfflineTree.Inputs = ["Phys/Lb_SL_Lc4pi/Particles"]
+# Lb_SL_Lc4pi_OfflineTree.Decay = decay0_Lc4pi
+# Lb_SL_Lc4pi_OfflineTree.Branches = {
+#     "Lambda_b0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
+#     "Lambda_c_STAR": "[Lambda_b0 -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
+#     "mu": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) ^mu- ]CC",
+#     "Lambda_c": "[Lambda_b0 -> (Lambda_c(2595)+ -> ^(Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
+#     "pi_plus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> ^pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
+#     "pi_minus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ ^pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
+#     "pi_plus2": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> ^pi+ pi-)) mu- ]CC",
+#     "pi_minus2": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ ^pi-)) mu- ]CC",
+#     "omega": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) ^(omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
+#     "omega2": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) (omega(782) -> pi+ pi-) ^(omega(782) -> pi+ pi-)) mu- ]CC",
+#     "p": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- ^p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
+#     "pi_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ ^pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC",
+#     "K_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) (omega(782) -> pi+ pi-) (omega(782) -> pi+ pi-)) mu- ]CC"
+# }
 #Lb_SL_Lc4pi_OfflineTree.Branches = {
 #    "Lambda_b0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi- pi+ pi-) mu- ]CC",
 #    "Lambda_c_STAR": "[Lambda_b0 -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi- pi+ pi-) mu- ]CC",
@@ -223,38 +227,38 @@
 # }
 
 
-########## Tree for Lc** -> (Lc* -> Lc pi pi) pi0 ####################################################
-Lc_starpi0pipi = CombineParticles("Lc_starpi0pipi")
-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<3500) & (M>2300) & (VFASPF(VCHI2/VDOF) < 7.0) "
-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" }
-Lc_starpi0pipi.CombinationCut = "ACUTDOCA(1.0*mm,'')"
+# ########## Tree for Lc** -> (Lc* -> Lc pi pi) pi0 ####################################################
+# Lc_starpi0pipi = CombineParticles("Lc_starpi0pipi")
+# 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<3500) & (M>2300) & (VFASPF(VCHI2/VDOF) < 7.0) "
+# 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" }
+# Lc_starpi0pipi.CombinationCut = "ACUTDOCA(1.0*mm,'')"
 
-Lb_SL_Lcpi0pipi = CombineParticles("Lb_SL_Lcpi0pipi")
-Lb_SL_Lcpi0pipi.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/Lc_starpi0pipi/Particles"]
-Lb_SL_Lcpi0pipi.DecayDescriptors = [ "[Lambda_b0 -> Lambda_c(2595)+ mu-]cc" ]
-Lb_SL_Lcpi0pipi.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
+# Lb_SL_Lcpi0pipi = CombineParticles("Lb_SL_Lcpi0pipi")
+# Lb_SL_Lcpi0pipi.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/Lc_starpi0pipi/Particles"]
+# Lb_SL_Lcpi0pipi.DecayDescriptors = [ "[Lambda_b0 -> Lambda_c(2595)+ mu-]cc" ]
+# Lb_SL_Lcpi0pipi.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
 
-decay0_pi0pipi = "[Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^pi0 ^pi+ ^pi-) ^mu- ]CC"
+# decay0_pi0pipi = "[Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^pi0 ^pi+ ^pi-) ^mu- ]CC"
 
-Lb_SL_Lcpi0pipi_OfflineTree = DecayTreeTuple("Lb_SL_Lcpi0pipi_OfflineTree")
-Lb_SL_Lcpi0pipi_OfflineTree.Inputs = ["Phys/Lb_SL_Lcpi0pipi/Particles"]
-Lb_SL_Lcpi0pipi_OfflineTree.Decay = decay0_pi0pipi
-Lb_SL_Lcpi0pipi_OfflineTree.Branches = {
-    "Lambda_b0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi0 pi+ pi-) mu- ]CC",
-    "Lambda_c_STAR": "[Lambda_b0 -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi0 pi+ pi-) mu- ]CC",
-    "mu": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi0 pi+ pi-) ^mu- ]CC",
-    "Lambda_c": "[Lambda_b0 -> (Lambda_c(2595)+ -> ^(Lambda_c+ -> K- p+ pi+ ) pi0 pi+ pi-) mu- ]CC",
-    "pi0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) ^pi0 pi+ pi-) mu- ]CC",
-    "pi_plus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi0 ^pi+ pi-) mu- ]CC",
-    "pi_minus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi0 pi+ ^pi-) mu- ]CC",
-    "p": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- ^p+ pi+ ) pi0 pi+ pi-) mu- ]CC",
-    "pi_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ ^pi+ ) pi0 pi+ pi-) mu- ]CC",
-    "K_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) pi0 pi+ pi-) mu- ]CC",
-}
+# Lb_SL_Lcpi0pipi_OfflineTree = DecayTreeTuple("Lb_SL_Lcpi0pipi_OfflineTree")
+# Lb_SL_Lcpi0pipi_OfflineTree.Inputs = ["Phys/Lb_SL_Lcpi0pipi/Particles"]
+# Lb_SL_Lcpi0pipi_OfflineTree.Decay = decay0_pi0pipi
+# Lb_SL_Lcpi0pipi_OfflineTree.Branches = {
+#     "Lambda_b0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi0 pi+ pi-) mu- ]CC",
+#     "Lambda_c_STAR": "[Lambda_b0 -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi0 pi+ pi-) mu- ]CC",
+#     "mu": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi0 pi+ pi-) ^mu- ]CC",
+#     "Lambda_c": "[Lambda_b0 -> (Lambda_c(2595)+ -> ^(Lambda_c+ -> K- p+ pi+ ) pi0 pi+ pi-) mu- ]CC",
+#     "pi0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) ^pi0 pi+ pi-) mu- ]CC",
+#     "pi_plus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi0 ^pi+ pi-) mu- ]CC",
+#     "pi_minus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi0 pi+ ^pi-) mu- ]CC",
+#     "p": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- ^p+ pi+ ) pi0 pi+ pi-) mu- ]CC",
+#     "pi_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ ^pi+ ) pi0 pi+ pi-) mu- ]CC",
+#     "K_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) pi0 pi+ pi-) mu- ]CC",
+# }
 
 
 ########## Tree for Sc*++ -> Lc pi+ ##################################################################
@@ -286,68 +290,71 @@
 # }
 
 
-########## Tree for Sc*++ -> (Lc* -> Lc pi+ pi-) pi+ ##################################################################
-Lc_star3pi = CombineParticles("Lc_star3pi")
-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)  & (PROBNNghost<0.3)", "pi+": "(PT>0)  & (PROBNNghost<0.3)", "pi-": "(PT>0)  & (PROBNNghost<0.3)", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
-Lc_star3pi.CombinationCut = "ACUTDOCA(1.0*mm,'')"
+########## Tree for Lb0 -> (X+ -> Lc*+ pi0) mu- ###################################################
+X_pi0 = CombineParticles('X_pi0')
+X_pi0.Inputs = ['Phys/Lc_star/Particles']
+for pc in pi0containers:
+    X_pi0.Inputs.append('Phys/%s/Particles' % pc)
+X_pi0.DecayDescriptors = ['[Sigma_c+ -> Lambda_c(2595)+ pi0]cc']
+X_pi0.MotherCut = LbMassCut + " & (VFASPF(VCHI2/VDOF) < 7.0) "
+X_pi0.DaughtersCuts = { "pi0": "(PT>0)", 'Lambda_c(2595)+': '(PT>0)' }
+X_pi0.CombinationCut = "ACUTDOCA(1.0*mm,'')"
 
-Lb_SL_Sc3pi = CombineParticles("Lb_SL_Sc3pi")
-Lb_SL_Sc3pi.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/Lc_star3pi/Particles"]
-Lb_SL_Sc3pi.DecayDescriptors = [ "[Lambda_b0 -> Sigma_c+ mu-]cc" ]
-Lb_SL_Sc3pi.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
+Lb_X_pi0 = CombineParticles("Lb_X_pi0")
+Lb_X_pi0.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/X_pi0/Particles"]
+Lb_X_pi0.DecayDescriptors = [ "[Lambda_b0 -> Sigma_c+ mu-]cc" ]
+Lb_X_pi0.MotherCut = LbMassCut + " & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
 
-decay0_3pi = "[Lambda_b0 -> ^(Sigma_c+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^pi+ ^pi+ ^pi-) ^mu- ]CC"
+decay_pi0 = "[Lambda_b0 -> ^(Sigma_c+ -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^pi+ ^pi-) ^pi0) ^mu- ]CC"
 
-Lb_SL_Sc3pi_OfflineTree = DecayTreeTuple("Lb_SL_Sc3pi_OfflineTree")
-Lb_SL_Sc3pi_OfflineTree.Inputs = ["Phys/Lb_SL_Sc3pi/Particles"]
-Lb_SL_Sc3pi_OfflineTree.Decay = decay0_3pi
-Lb_SL_Sc3pi_OfflineTree.Branches = {
-    "Lambda_b0": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi+ pi-) mu- ]CC",
-    "Lambda_c_STAR": "[Lambda_b0 -> ^(Sigma_c+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi+ pi-) mu- ]CC",
-    "mu": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi+ pi-) ^mu- ]CC",
-    "Lambda_c": "[Lambda_b0 -> (Sigma_c+ -> ^(Lambda_c+ -> K- p+ pi+ ) pi+ pi+ pi-) mu- ]CC",
-    "pi_plus1": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c+ -> K- p+ pi+ ) ^pi+ pi+ pi-) mu- ]CC",
-    "pi_plus2": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ ^pi+ pi-) mu- ]CC",
-    "pi_minus": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi+ ^pi-) mu- ]CC",
-    "p": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c+ -> K- ^p+ pi+ ) pi+ pi+ pi-) mu- ]CC",
-    "pi_Lc": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c+ -> K- p+ ^pi+ ) pi+ pi+ pi-) mu- ]CC",
-    "K_Lc": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c+ -> ^K- p+ pi+ ) pi+ pi+ pi-) mu- ]CC"
+Lb_X_pi0_OfflineTree = DecayTreeTuple("Lb_X_pi0_OfflineTree")
+Lb_X_pi0_OfflineTree.Inputs = ["Phys/Lb_X_pi0/Particles"]
+Lb_X_pi0_OfflineTree.Decay = decay_pi0
+Lb_X_pi0_OfflineTree.Branches = {
+    "Lambda_b0": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) pi0) mu- ]CC",
+    "X": "[Lambda_b0 -> ^(Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) pi0) mu- ]CC",
+    "Lambda_c_STAR": "[Lambda_b0 -> (Sigma_c+ -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) pi0) mu- ]CC",
+    "pi0": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) ^pi0) mu- ]CC",
+    "pi_plus": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) ^pi+ pi-) pi0) mu- ]CC",
+    "pi_minus": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ ^pi-) pi0) mu- ]CC",
+    "pi_plus_Lc": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ ^pi+ ) pi+ pi-) pi0) mu- ]CC",
+    "p": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- ^p+ pi+ ) pi+ pi-) pi0) mu- ]CC",
+    "K": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) pi+ pi-) pi0) mu- ]CC",
+    "mu": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) pi0) ^mu- ]CC",
 }
 
+########## Tree for Lb0 -> (X+ -> Lc*+ pi+ pi-) mu- ###############################################
+X_2pi = CombineParticles('X_2pi')
+X_2pi.Inputs = ['Phys/Lc_star/Particles', 'Phys/StdAllLoosePions/Particles']
+X_2pi.DecayDescriptors = ['[Sigma_c+ -> Lambda_c(2595)+ pi+ pi-]cc']
+X_2pi.MotherCut = LbMassCut + " & (VFASPF(VCHI2/VDOF) < 7.0) "
+X_2pi.DaughtersCuts = { "pi+": "(PT>0) & (PROBNNghost<0.4)", 'pi-': "(PT>0) & (PROBNNghost<0.4)", 'Lambda_c(2595)+': '(PT>0)' }
+X_2pi.CombinationCut = "ACUTDOCA(1.0*mm,'')"
 
-# ########## Tree for Lc* -> Lc gamma ##################################################################
-# Lc_starG = CombineParticles("Lc_starG")
-# Lc_starG.Inputs = ["Phys/Lc2PKPiforB2DMuNuX/Particles"]
-# for g in ['StdAllLooseGammaDD', 'StdAllLooseGammaLL', 'StdAllTightGammaDD', 'StdAllTightGammaLL',
-#           'StdAllTightSymGammaDD', 'StdAllTightSymGammaLL', 'StdDiElectronGamma']:
-#     Lc_starG.Inputs.append('Phys/%s/Particles' % g)
-# Lc_starG.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ gamma]cc"]
-# Lc_starG.MotherCut = "(M<3500) & (M>2300) & (VFASPF(VCHI2/VDOF) < 7.0) "
-# Lc_starG.DaughtersCuts = { "gamma": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
+Lb_X_2pi = CombineParticles("Lb_X_2pi")
+Lb_X_2pi.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/X_2pi/Particles"]
+Lb_X_2pi.DecayDescriptors = [ "[Lambda_b0 -> Sigma_c+ mu-]cc" ]
+Lb_X_2pi.MotherCut = LbMassCut + " & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
 
-# Lb_SL_LcG = CombineParticles("Lb_SL_LcG")
-# Lb_SL_LcG.Inputs = [ "Phys/MuforB2DMuNuX/Particles", "Phys/Lc_starG/Particles"]
-# Lb_SL_LcG.DecayDescriptors = [ "[Lambda_b0 -> Lambda_c(2595)+ mu-]cc" ]
-# Lb_SL_LcG.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
+decay_2pi = "[Lambda_b0 -> ^(Sigma_c+ -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^pi+ ^pi-) ^pi+ ^pi-) ^mu- ]CC"
 
-# decay0_G = "[Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^gamma) ^mu- ]CC"
+Lb_X_2pi_OfflineTree = DecayTreeTuple("Lb_X_2pi_OfflineTree")
+Lb_X_2pi_OfflineTree.Inputs = ["Phys/Lb_X_2pi/Particles"]
+Lb_X_2pi_OfflineTree.Decay = decay_2pi
+Lb_X_2pi_OfflineTree.Branches = {
+    "Lambda_b0": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) pi+ pi-) mu- ]CC",
+    "X": "[Lambda_b0 -> ^(Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) pi+ pi-) mu- ]CC",
+    "Lambda_c_STAR": "[Lambda_b0 -> (Sigma_c+ -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) pi+ pi-) mu- ]CC",
+    "pi_plus_X": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) ^pi+ pi-) mu- ]CC",
+    "pi_minus_X": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) pi+ ^pi-) mu- ]CC",
+    "pi_plus": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) ^pi+ pi-) pi+ pi-) mu- ]CC",
+    "pi_minus": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ ^pi-) pi+ pi-) mu- ]CC",
+    "pi_plus_Lc": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ ^pi+ ) pi+ pi-) pi+ pi-) mu- ]CC",
+    "p": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- ^p+ pi+ ) pi+ pi-) pi+ pi-) mu- ]CC",
+    "K": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) pi+ pi-) pi+ pi-) mu- ]CC",
+    "mu": "[Lambda_b0 -> (Sigma_c+ -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) pi+ pi-) ^mu- ]CC",
+}
 
-# Lb_SL_LcG_OfflineTree = DecayTreeTuple("Lb_SL_LcG_OfflineTree")
-# Lb_SL_LcG_OfflineTree.Inputs = ["Phys/Lb_SL_LcG/Particles"]
-# Lb_SL_LcG_OfflineTree.Decay = decay0_G
-# Lb_SL_LcG_OfflineTree.Branches = {
-#     "Lambda_b0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) gamma) mu- ]CC",
-#     "Lambda_c_STAR": "[Lambda_b0 -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) gamma) mu- ]CC",
-#     "mu": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) gamma) ^mu- ]CC",
-#     "Lambda_c": "[Lambda_b0 -> (Lambda_c(2595)+ -> ^(Lambda_c+ -> K- p+ pi+ ) gamma) mu- ]CC",
-#     "gamma": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) ^gamma) mu- ]CC",
-#     "p": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- ^p+ pi+ ) gamma) mu- ]CC",
-#     "pi_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ ^pi+ ) gamma) mu- ]CC",
-#     "K_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) gamma) mu- ]CC"
-# }
 
 
 ############################## Configuring all DecayTreeTuples #######################################
@@ -356,14 +363,12 @@
 from Configurables import TupleToolGeometry
 from Configurables import TupleToolTrackIsolation
 
-Lc_star_list = [Lc_star3pi, Lc_starpi0pipi, Lc_star4pi, Lc_star]
-Lb_list = [Lb_SL_Sc3pi, Lb_SL_Lcpi0pipi, Lb_SL_Lc4pi, Lb_SL]
+combine_list = [Lc_star, X_pi0, X_2pi, Lc_starpi0pipi, Lc_star4pi]
+combine_list += [Lb_X_pi0, Lb_X_2pi, Lb_SL_Lcpi0pipi, Lb_SL_Lc4pi, Lb_SL]
 
-DTTlist = [# Lb_SL_LcG_OfflineTree,
-           Lb_SL_Sc3pi_OfflineTree,
-           # Lb_SL_Sc1pi_OfflineTree,
+DTTlist = [Lb_X_pi0_OfflineTree,
+           Lb_X_2pi_OfflineTree,
            Lb_SL_Lcpi0pipi_OfflineTree,
-           # Lb_SL_Lc1pi0_OfflineTree,
            Lb_SL_OfflineTree,
            Lb_SL_Lc4pi_OfflineTree]
 
@@ -431,10 +436,8 @@
                                 HltDecReportsDecoder()]
 # Add particles
 DaVinci().UserAlgorithms.append(omega_2pi)
-for lc in Lc_star_list:
-    DaVinci().UserAlgorithms.append(lc)
-for lb in Lb_list:
-    DaVinci().UserAlgorithms.append(lb)
+for part in combine_list:
+    DaVinci().UserAlgorithms.append(part)
 for DTT in DTTlist:
     DaVinci().UserAlgorithms.append(DTT)
 DaVinci().TupleFile = "Lb_Lcstarmunu.root"