diff --git a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_1pi.py b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_1pi.py
index a239110..d70ca2a 100644
--- a/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_1pi.py
+++ b/ganga/Lb2Lcmunu/Lb2Lcmunu_S20_data_NEWJET_1pi.py
@@ -136,7 +136,7 @@
 Lc_star4pi = CombineParticles("Lc_star4pi")
 Lc_star4pi.Inputs = ["Phys/StdAllLoosePions/Particles", "Phys/Lc2PKPiforB2DMuNuX/Particles"]
 Lc_star4pi.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ pi+ pi- pi+ pi-]cc"]
-Lc_star4pi.MotherCut = "(M<3000) & (M>2300) & (VFASPF(VCHI2/VDOF) < 7.0) "
+Lc_star4pi.MotherCut = "(M<3500) & (M>2300) & (VFASPF(VCHI2/VDOF) < 15.0) "
 Lc_star4pi.DaughtersCuts = { "pi+": "PT>0", "pi-": "PT>0", "Lambda_c+": "ADMASS('Lambda_c+') < 50.0 *MeV" }
 
 Lb_SL_Lc4pi = CombineParticles("Lb_SL_Lc4pi")
@@ -197,7 +197,7 @@
 Lc_starpi0pipi = CombineParticles("Lc_starpi0pipi")
 Lc_starpi0pipi.Inputs = ["Phys/StdLooseMergedPi0/Particles", "Phys/Lc2PKPiforB2DMuNuX/Particles", "Phys/StdAllLoosePions/Particles"]
 Lc_starpi0pipi.DecayDescriptors = ["[Lambda_c(2595)+ -> Lambda_c+ pi0 pi+ pi-]cc"]
-Lc_starpi0pipi.MotherCut = "(M<3000) & (M>2300) & (VFASPF(VCHI2/VDOF) < 7.0) "
+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" }
 
 Lb_SL_Lcpi0pipi = CombineParticles("Lb_SL_Lcpi0pipi")
@@ -322,13 +322,18 @@
 from Configurables import TupleToolGeometry
 from Configurables import TupleToolTrackIsolation
 
-for DTT in [Lb_SL_LcG_OfflineTree,
-            Lb_SL_Sc3pi_OfflineTree,
-            Lb_SL_Sc1pi_OfflineTree,
-            Lb_SL_Lcpi0pipi_OfflineTree,
-            Lb_SL_Lc1pi0_OfflineTree,
-            Lb_SL_OfflineTree,
-            Lb_SL_Lc4pi_OfflineTree]:
+Lc_star_list = [Lc_starG, Lc_star3pi, Lc_star1pi, Lc_starpi0pipi, Lc_star1pi0, Lc_star4pi, Lc_star]
+Lb_list = [Lb_SL_LcG, Lb_SL_Sc3pi, Lb_SL_Sc1pi, Lb_SL_Lcpi0pipi, Lb_SL_Lc1pi0, Lb_SL_Lc4pi, Lb_SL]
+
+DTTlist = [Lb_SL_LcG_OfflineTree,
+           Lb_SL_Sc3pi_OfflineTree,
+           Lb_SL_Sc1pi_OfflineTree,
+           Lb_SL_Lcpi0pipi_OfflineTree,
+           Lb_SL_Lc1pi0_OfflineTree,
+           Lb_SL_OfflineTree,
+           Lb_SL_Lc4pi_OfflineTree]
+
+for DTT in DTTlist:
 
     DTT.ToolList = [
         "TupleToolKinematic",
@@ -378,11 +383,6 @@
     DTT.TrackIsolation.Verbose = True
 
 ##################################################
-#DaVinci().Input = ([
-#'root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2012/ALLSTREAMS.DST/00025296/0000/00025296_00000002_1.allstreams.dst'
-#'root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2012/ALLSTREAMS.DST/00025296/0000/00025296_00000003_1.allstreams.dst',
-#'root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2012/ALLSTREAMS.DST/00025296/0000/00025296_00000004_1.allstreams.dst'
-#])
 from Configurables import DaVinci
 #DaVinci().EventPreFilters = fltrs.filters('Filters')
 DaVinci().EvtMax = -1                          # Number of events
@@ -394,9 +394,13 @@
 DaVinci().appendToMainSequence( [ sc.sequence() ] )
 DaVinci().UserAlgorithms = [    HltSelReportsDecoder(),
                                 HltVertexReportsDecoder(),
-                                HltDecReportsDecoder(),
-                                Lc_star, Lb_SL,
-                                Lb_SL_OfflineTree]
+                                HltDecReportsDecoder()]
+for lc in Lc_star_list:
+    DaVinci().UserAlgorithms.append(lc)
+for lb in Lb_list:
+    DaVinci().UserAlgorithms.append(lb)
+for DTT in DTTlist:
+    DaVinci().UserAlgorithms.append(DTT)
 DaVinci().TupleFile = "Lb_Lcstarmunu.root"
 
 #DaVinci().CondDBtag = "Sim08-20130503-1-vc-md100"
diff --git a/ganga/Lb2Lcmunu/testfile.py b/ganga/Lb2Lcmunu/testfile.py
new file mode 100644
index 0000000..c707696
--- /dev/null
+++ b/ganga/Lb2Lcmunu/testfile.py
@@ -0,0 +1,13 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+# @Author: Elena Graverini
+# @Date:   2016-01-14 12:04:32
+# @Last Modified by:   Elena Graverini
+# @Last Modified time: 2016-01-14 12:05:23
+from Configurables import DaVinci
+DaVinci().Input = ([
+    "DATAFILE='PFN:00021317_00001395_1.semileptonic.dst' TYP='POOL_ROOTTREE' OPT='READ'"
+    #'root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2012/ALLSTREAMS.DST/00025296/0000/00025296_00000002_1.allstreams.dst'
+    #'root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2012/ALLSTREAMS.DST/00025296/0000/00025296_00000003_1.allstreams.dst',
+    #'root://eoslhcb.cern.ch//eos/lhcb/grid/prod/lhcb/MC/2012/ALLSTREAMS.DST/00025296/0000/00025296_00000004_1.allstreams.dst'
+])