diff --git a/ganga/Lb2Lcmunu/Lb2Lcmunu_S24_data2015_NEWJET.py b/ganga/Lb2Lcmunu/Lb2Lcmunu_S24_data2015_NEWJET.py
index 82fa4dd..6de8562 100644
--- a/ganga/Lb2Lcmunu/Lb2Lcmunu_S24_data2015_NEWJET.py
+++ b/ganga/Lb2Lcmunu/Lb2Lcmunu_S24_data2015_NEWJET.py
@@ -127,6 +127,13 @@
 Lb_SL.DecayDescriptors = [ "[Lambda_b0 -> Lambda_c(2595)+ mu-]cc" ]
 Lb_SL.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
 
+############## WRONG SIGN COMBINATION #######################################
+
+Lb_SL_WS = CombineParticles("Lb_SL_WS")
+Lb_SL_WS.Inputs = [ "Phys/MuforB2DMuForTauMu/Particles", "Phys/Lc_star/Particles"]
+Lb_SL_WS.DecayDescriptors = [ "[Lambda_b0 -> Lambda_c(2595)+ mu+]cc" ]
+Lb_SL_WS.MotherCut = "(M<6500) & (VFASPF(VCHI2/VDOF) < 7.0) & (BPVDIRA>0.999)"
+
 
 '''
 ################################################################################
@@ -143,19 +150,14 @@
 
 
 ########################################################################
-#decay = "[ B*_s20 -> ^(B- -> ^(J/psi(1S) -> ^mu+ ^mu-) ^K-) ^[pi+]CC]CC"
 decay0 = "[Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^pi+ ^pi-) ^mu- ]CC"
-#decay0 = "[H_10 -> ^( B*_s20 -> ^K+ ^(B- -> ^(D0 -> ^K- ^pi+) ^mu-)) ^CELLjet ]CC"
-#decay0 = "[H_10 -> ^( Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ ->  K- p+ pi+  )   ) ^mu- ) ^CELLjet ]CC"
-#decay1 = "[H_10 -> ^(  B*_s20-> ^K- ^(B- -> ^(D0 -> ^K- ^pi+) ^mu-)) ^CELLjet ]CC"
+decay0_WS = "[Lambda_b0 -> ^(Lambda_c(2595)+ -> ^(Lambda_c+ -> ^K- ^p+ ^pi+ ) ^pi+ ^pi-) ^mu+ ]CC"
 ########################################################################
 Lb_SL_OfflineTree = DecayTreeTuple("Lb_SL_OfflineTree")
 Lb_SL_OfflineTree.Inputs = ["Phys/Lb_SL/Particles"]
 Lb_SL_OfflineTree.Decay = decay0
 Lb_SL_OfflineTree.Branches = {
-    #   "H_10": "[^H_10 -> (Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ) CELLjet]CC",
     "Lambda_b0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ]CC",
-    #   "CELLjet": "[H_10 -> (Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ) ^CELLjet]CC",
     "Lambda_c_STAR": "[Lambda_b0 -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ]CC",
     "mu": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) ^mu- ]CC",
     "Lambda_c": "[Lambda_b0 -> (Lambda_c(2595)+ -> ^(Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu- ]CC",
@@ -165,67 +167,75 @@
     "pi_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ ^pi+ ) pi+ pi-) mu- ]CC",
     "K_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) pi+ pi-) mu- ]CC"
 }
-Lb_SL_OfflineTree.ToolList = [
-    "TupleToolKinematic",
-    "TupleToolPid",
-    "TupleToolGeometry",
-    "TupleToolPrimaries",
-    "TupleToolTrackInfo",
-    "TupleToolEventInfo",
-    "TupleToolRecoStats",
-    "TupleToolAngles",
-    "TupleToolNeutrinoReco",
-    "TupleToolMCBackgroundInfo",
-    "TupleToolMCTruth",
-    "TupleToolJets",
-    'TupleToolTrigger',
-    'TupleToolTISTOS',
-    'TupleToolRICHPid',
-    'TupleToolDira',
-    'TupleToolEventInfo',
-    'TupleToolPropertime',
-    'TupleToolRecoStats',
-    'TupleToolTrackPosition'
-]
 
+Lb_SL_WS_OfflineTree = DecayTreeTuple("Lb_SL_WS_OfflineTree")
+Lb_SL_WS_OfflineTree.Inputs = ["Phys/Lb_SL_WS/Particles"]
+Lb_SL_WS_OfflineTree.Decay = decay0_WS
+Lb_SL_WS_OfflineTree.Branches = {
+    "Lambda_b0": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu+ ]CC",
+    "Lambda_c_STAR": "[Lambda_b0 -> ^(Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu+ ]CC",
+    "mu": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ pi-) ^mu+ ]CC",
+    "Lambda_c": "[Lambda_b0 -> (Lambda_c(2595)+ -> ^(Lambda_c+ -> K- p+ pi+ ) pi+ pi-) mu+ ]CC",
+    "pi_plus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) ^pi+ pi-) mu+ ]CC",
+    "pi_minus": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ pi+ ) pi+ ^pi-) mu+ ]CC",
+    "p": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- ^p+ pi+ ) pi+ pi-) mu+ ]CC",
+    "pi_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> K- p+ ^pi+ ) pi+ pi-) mu+ ]CC",
+    "K_Lc": "[Lambda_b0 -> (Lambda_c(2595)+ -> (Lambda_c+ -> ^K- p+ pi+ ) pi+ pi-) mu+ ]CC"
+}
 
-Lb_SL_OfflineTree.addTool( TupleToolTrigger, name='TupleToolTrigger' )
-Lb_SL_OfflineTree.TupleToolTrigger.Verbose = True
-Lb_SL_OfflineTree.TupleToolTrigger.TriggerList = trigger_list
-
-
-Lb_SL_OfflineTree.addTool( TupleToolTISTOS, name='TupleToolTISTOS' )
-Lb_SL_OfflineTree.TupleToolTISTOS.Verbose = True
-Lb_SL_OfflineTree.TupleToolTISTOS.TriggerList = trigger_list
-Lb_SL_OfflineTree.ToolList += [ "TupleToolTISTOS" ]
-
-
-########################################################################
 from Configurables import TupleToolNeutrinoReco
-Lb_SL_OfflineTree.addTool( TupleToolNeutrinoReco, name = "TupleToolNeutrinoReco" )
-Lb_SL_OfflineTree.TupleToolNeutrinoReco.Verbose = True
-Lb_SL_OfflineTree.TupleToolNeutrinoReco.MotherMass = 5619.5  # MeV
-Lb_SL_OfflineTree.TupleToolNeutrinoReco.ExtraName = 'mLb'
-
 from Configurables import TupleToolKinematic
-Lb_SL_OfflineTree.addTool( TupleToolKinematic, name = "TupleToolKinematic" )
-Lb_SL_OfflineTree.TupleToolKinematic.Verbose = True
-
 from Configurables import TupleToolGeometry
-Lb_SL_OfflineTree.addTool( TupleToolGeometry, name = "TupleToolGeometry" )
-Lb_SL_OfflineTree.TupleToolGeometry.Verbose = True
-
-Lb_SL_OfflineTree.addTool(TupleToolRecoStats, name="TupleToolRecoStats")
-Lb_SL_OfflineTree.TupleToolRecoStats.Verbose = True
-Lb_SL_OfflineTree.UseLabXSyntax = True
-Lb_SL_OfflineTree.RevertToPositiveID = False
-
 from Configurables import TupleToolTrackIsolation
-Lb_SL_OfflineTree.addTupleTool('TupleToolTrackIsolation/TrackIsolation')
-Lb_SL_OfflineTree.TrackIsolation.MinConeAngle = 0.5
-Lb_SL_OfflineTree.TrackIsolation.MaxConeAngle = 1.5
-Lb_SL_OfflineTree.TrackIsolation.StepSize = 0.1
-Lb_SL_OfflineTree.TrackIsolation.Verbose = True
+
+
+for DTT in [Lb_SL_OfflineTree, Lb_SL_WS_OfflineTree]:
+    DTT.ToolList = [
+        "TupleToolKinematic",
+        "TupleToolPid",
+        "TupleToolGeometry",
+        "TupleToolPrimaries",
+        "TupleToolTrackInfo",
+        "TupleToolEventInfo",
+        "TupleToolRecoStats",
+        "TupleToolAngles",
+        "TupleToolNeutrinoReco",
+        "TupleToolMCBackgroundInfo",
+        "TupleToolMCTruth",
+        "TupleToolJets",
+        'TupleToolTrigger',
+        'TupleToolTISTOS',
+        'TupleToolRICHPid',
+        'TupleToolDira',
+        'TupleToolEventInfo',
+        'TupleToolPropertime',
+        'TupleToolRecoStats',
+        'TupleToolTrackPosition'
+    ]
+    DTT.addTool( TupleToolTrigger, name='TupleToolTrigger' )
+    DTT.TupleToolTrigger.Verbose = True
+    DTT.TupleToolTrigger.TriggerList = trigger_list
+    DTT.addTool( TupleToolTISTOS, name='TupleToolTISTOS' )
+    DTT.TupleToolTISTOS.Verbose = True
+    DTT.TupleToolTISTOS.TriggerList = trigger_list
+    DTT.ToolList += [ "TupleToolTISTOS" ]
+    DTT.addTool( TupleToolNeutrinoReco, name = "TupleToolNeutrinoReco" )
+    DTT.TupleToolNeutrinoReco.Verbose = True
+    DTT.TupleToolNeutrinoReco.MotherMass = 5619.5  # MeV
+    DTT.TupleToolNeutrinoReco.ExtraName = 'mLb'
+    DTT.addTool( TupleToolKinematic, name = "TupleToolKinematic" )
+    DTT.TupleToolKinematic.Verbose = True
+    DTT.addTool( TupleToolGeometry, name = "TupleToolGeometry" )
+    DTT.TupleToolGeometry.Verbose = True
+    DTT.addTool(TupleToolRecoStats, name="TupleToolRecoStats")
+    DTT.TupleToolRecoStats.Verbose = True
+    DTT.UseLabXSyntax = True
+    DTT.RevertToPositiveID = False
+    DTT.addTupleTool('TupleToolTrackIsolation/TrackIsolation')
+    DTT.TrackIsolation.MinConeAngle = 0.5
+    DTT.TrackIsolation.MaxConeAngle = 1.5
+    DTT.TrackIsolation.StepSize = 0.1
+    DTT.TrackIsolation.Verbose = True
 
 
 ##################################################
@@ -246,8 +256,8 @@
 DaVinci().UserAlgorithms = [    HltSelReportsDecoder(),
                                 HltVertexReportsDecoder(),
                                 HltDecReportsDecoder(),
-                                Lc_star, Lb_SL,
-                                Lb_SL_OfflineTree]
+                                Lc_star, Lb_SL, Lb_SL_WS
+                                Lb_SL_OfflineTree, Lb_SL_WS_OfflineTree]
 DaVinci().TupleFile = "Lb_Lcstarmunu.root"
 
 #DaVinci().CondDBtag = "Sim08-20130503-1-vc-md100"