diff --git a/Lb2pKll/__init__.py b/Lb2pKll/__init__.py
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/Lb2pKll/dv_pKll.py b/Lb2pKll/dv_pKll.py
new file mode 100755
index 0000000000000000000000000000000000000000..953f931e9d4cc97e17571dffee85e3b3409332c4
--- /dev/null
+++ b/Lb2pKll/dv_pKll.py
@@ -0,0 +1,570 @@
+# References:
+# https://twiki.cern.ch/twiki/pub/LHCbPhysics/Lst2pkll/LHCb-ANA-2019-007_v3.1.pdf
+# https://gitlab.cern.ch/LHCb-RD/ewp-rpk-fullr1r2/-/blob/master/DaVinci/LeptonU-options-NOPID.py
+# https://gitlab.cern.ch/LHCb-RD/ewp-rks-rkstplus/-/blob/master/Produce_Ntuples/Data/Bu2LLK_DaVinci.py
+# https://gitlab.cern.ch/LHCb-RD/ewp-bd2ksteeangular-legacy/-/blob/master/samples_and_selection/DaVinci_scripts/options_Bu2LLK_EE.py
+
+import sys
+sys.path.append("./")
+#from steering_files.MC_2018_MD_ReStrip0_Filt1_mumu import * # BD localtest
+#from steering_files.MC_2017_MU_ReStrip0_Filt1_ee import * # BD localtest
+#from steering_files.MC_2012_MD_ReStrip0_Filt1_mumu import * # BD localtest
+#from steering_files.MC_2012_MD_ReStrip1_Filt0_mumu import * # BD localtest
+#from steering_files.Data_2012 import * # BD localtest
+#from steering_files.Data_2017 import * # BD localtest
+
+import os
+isRun2      = (os.getenv("ISRUN2", 'False') == 'True')     
+isMC        = (os.getenv("ISMC", 'False') == 'True') 
+year        = os.getenv("YEAR")   
+inputType   = os.getenv("INPUTTYPE")  
+isReDecay   = (os.getenv("ISREDECAY", 'False') == 'True')  
+isfiltered  = (os.getenv("ISFILTERED", 'False') == 'True')  
+magpos      = os.getenv("MAGPOS")    
+restrip     = (os.getenv("RESTRIP", 'False') == 'True') 
+doMuon      = (os.getenv("DOMUON", 'False') == 'True')  
+doElectron  = (os.getenv("DOELECTRON", 'False') == 'True') 
+
+
+from Lb2pKll.helpers.helper import *
+#from helpers.helper import *  # BD localtest
+from LHCbKernel.Configuration import *
+from Configurables import DaVinci, TupleToolDecayTreeFitter, DecayTreeTuple, TupleToolDecay, TupleToolTrigger, TupleToolRecoStats, TupleToolTISTOS, TupleToolGeometry, TupleToolTrackInfo, TupleToolMCTruth, TupleToolL0Data, SubstitutePID
+if isReDecay: from Configurables import MCTupleToolRedecay
+
+from Configurables import DstConf, CaloDstUnPackConf, LoKi__Hybrid__TupleTool
+from Configurables import EventNodeKiller, ProcStatusCheck
+from Configurables import StrippingReport, TimingAuditor, SequencerTimerTool
+
+from StrippingConf.Configuration import StrippingConf, StrippingStream
+from StrippingSettings.Utils import strippingConfiguration
+from StrippingArchive.Utils import buildStreams
+from StrippingArchive import strippingArchive
+
+from Gaudi.Configuration import *
+from DecayTreeTuple.Configuration import *
+
+from Configurables import CombineParticles
+from Configurables import TupleToolTrigger, SubstitutePID
+from Configurables import FilterDesktop
+from PhysSelPython.Wrappers import Selection, SelectionSequence, DataOnDemand, MergedSelection
+
+##################################################################
+# Initialize some settings independent of LHCb tools
+evtmax         =  -1 # change this for localtest. has to be -1 for AP
+###################################################################
+
+# RpK
+linename_e_Run1   = "Bu2LLK_eeLine2"
+linename_mu_Run1  = "Bu2LLK_mmLine"
+# new lines
+linename_e_Run2   = "Bu2LLK_ee_Lb2pkeeLooseLine"
+linename_mu_Run2  = "Bu2LLK_mm_Lb2pkmmLooseLine"
+
+#######################################################################
+
+# applied on the B head
+trig_L0TIS        = ['L0GlobalDecision','L0HadronDecision',"L0MuonDecision", "L0ElectronDecision"] # TIS-ed on
+trig_Run1_B_mu    = trig_L0TIS + ['Hlt1TrackAllL0Decision', 'Hlt1TrackMuonDecision', 'Hlt1TrackDiMuonDecision', 'Hlt2DiMuonDetachedDecision', 'Hlt2Topo2BodyBBDTDecision', 'Hlt2Topo3BodyBBDTDecision', 'Hlt2Topo4BodyBBDTDecision', 'Hlt2TopoMu2BodyBBDTDecision', 'Hlt2TopoMu3BodyBBDTDecision', 'Hlt2TopoMu4BodyBBDTDecision']
+trig_Run2_B_mu    = trig_L0TIS + ['Hlt1TrackMuonDecision','Hlt1TrackMVADecision', 'Hlt1TwoTrackMVADecision', 'Hlt1TrackMuonMVADecision', 'Hlt1DiMuonHighMassDecision', 'Hlt2DiMuonDetachedDecision', 'Hlt2Topo2BodyDecision','Hlt2Topo3BodyDecision','Hlt2Topo4BodyDecision','Hlt2TopoMu2BodyDecision','Hlt2TopoMu3BodyDecision','Hlt2TopoMu4BodyDecision','Hlt2TopoMuMu2BodyDecision','Hlt2TopoMuMu3BodyDecision','Hlt2TopoMuMu4BodyDecision']
+trig_Run1_B_e     = trig_L0TIS + ['Hlt1TrackAllL0Decision', 'Hlt1TrackAllL0TightDecision','Hlt2Topo2BodyBBDTDecision', 'Hlt2Topo3BodyBBDTDecision', 'Hlt2Topo4BodyBBDTDecision', 'Hlt2TopoE2BodyBBDTDecision', 'Hlt2TopoE3BodyBBDTDecision', 'Hlt2TopoE4BodyBBDTDecision']
+trig_Run2_B_e     = trig_L0TIS + ['Hlt1TrackMVADecision', 'Hlt1TwoTrackMVADecision', 'Hlt2Topo2BodyDecision', 'Hlt2Topo3BodyDecision', 'Hlt2Topo4BodyDecision', 'Hlt2TopoE2BodyDecision', 'Hlt2TopoE3BodyDecision', 'Hlt2TopoE4BodyDecision', 'Hlt2TopoEE2BodyDecision', 'Hlt2TopoEE3BodyDecision', 'Hlt2TopoEE4BodyDecision']
+
+# some Hlt1/2 muon lines for the dimuon
+trig_Run1_Jpsi_mu = ["L0MuonDecision","L0DiMuonDecision", "Hlt1DiMuonLowMassDecision","Hlt1DiMuonHighMassDecision","Hlt1TrackMuonDecision","Hlt2SingleMuonDecision","Hlt2DiMuonDetachedDecision"]
+trig_Run2_Jpsi_mu = ["L0MuonDecision","L0DiMuonDecision", "Hlt1TrackMuonMVADecision","Hlt1TrackMuonDecision","Hlt1SingleMuonHighPTDecision","Hlt1DiMuonLowMassDecision","Hlt1DiMuonHighMassDecision","Hlt2DiMuonDetachedDecision","Hlt2SingleMuonDecision"]
+trig_e            = ["L0ElectronDecision","L0ElectronHiDecision", "L0PhotonDecision", "L0PhotonHiDecision", "Hlt1SingleElectronNoIPDecision"]
+
+# applied on the pK
+trig_Lst          = ['L0HadronDecision']
+
+
+def addNtupleToDaVinci(name, isRun2, isMC, isDST):
+    trig_B = trig_Run1_B_mu
+    trig_Jpsi = trig_Run1_Jpsi_mu
+    line   = linename_mu_Run1
+    if "ee" in name and isRun2 == 0: 
+       trig_B = trig_Run1_B_e 
+       line   = linename_e_Run1
+    elif "mm" in name and isRun2 == 1:
+       trig_B = trig_Run2_B_mu 
+       trig_Jpsi = trig_Run2_Jpsi_mu
+       line   = linename_mu_Run2
+    elif "ee" in name and isRun2 == 1:
+       trig_B = trig_Run2_B_e 
+       line   = linename_e_Run2
+    theTuple = DecayTreeTuple(name)
+    if "ee" in name:
+       trig_Jpsi = trig_e
+       theTuple.Decay = "[Lambda_b0 -> ^(J/psi(1S) -> ^e+ ^e-) ^(Lambda(1520)0 -> ^p+ ^K-)]CC"
+       theTuple.addBranches ({
+         "Lb"  : "[Lambda_b0->  (J/psi(1S) -> e+ e-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "L1"  : "[Lambda_b0->  (J/psi(1S) -> ^e+ e-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "L2"  : "[Lambda_b0->  (J/psi(1S) -> e+ ^e-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "Jpsi": "[Lambda_b0->  ^(J/psi(1S) -> e+ e-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "p"   : "[Lambda_b0->  (J/psi(1S) -> e+ e-)  ( Lambda(1520)0 -> ^p+  K-) ]CC",
+         "k"   : "[Lambda_b0->  (J/psi(1S) -> e+ e-)  ( Lambda(1520)0 -> p+  ^K-) ]CC",
+         "Lst" : "[Lambda_b0->  (J/psi(1S) -> e+ e-)  ^( Lambda(1520)0 -> p+  K-)]CC"
+       })
+    else:
+        theTuple.Decay = "[Lambda_b0 -> ^(J/psi(1S) -> ^mu+ ^mu-) ^(Lambda(1520)0 -> ^p+ ^K-)]CC" 
+        theTuple.addBranches ({
+         "Lb"  : "[Lambda_b0->  (J/psi(1S) -> mu+ mu-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "L1"  : "[Lambda_b0->  (J/psi(1S) -> ^mu+ mu-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "L2"  : "[Lambda_b0->  (J/psi(1S) -> mu+ ^mu-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "Jpsi": "[Lambda_b0->  ^(J/psi(1S) -> mu+ mu-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "p"   : "[Lambda_b0->  (J/psi(1S) -> mu+ mu-)  ( Lambda(1520)0 -> ^p+  K-) ]CC",
+         "k"   : "[Lambda_b0->  (J/psi(1S) -> mu+ mu-)  ( Lambda(1520)0 -> p+  ^K-) ]CC",
+         "Lst" : "[Lambda_b0->  (J/psi(1S) -> mu+ mu-)  ^( Lambda(1520)0 -> p+  K-)]CC"
+        })
+    if isMC and restrip: 
+        if year == "2011" : 
+              strip = "stripping21r1p2"
+        elif year == "2012" : 
+              strip = "stripping21r0p2"
+        elif year == "2016" :
+              strip = "stripping28r2p1"
+        elif year == "2017" :
+              strip = "stripping29r2p2"
+        elif year == "2018" :
+              strip = "stripping34r0p2"
+        # The event node killer should be run before
+        event_node_killer = EventNodeKiller("StripKiller")
+        event_node_killer.Nodes = ["/Event/AllStreams", "/Event/Strip"]
+
+        DaVinci().appendToMainSequence( [ event_node_killer ] )
+        #print(year, strip)
+        streams = buildStreams(stripping=strippingConfiguration(strip), archive=strippingArchive(strip))
+
+        custom_stream = StrippingStream("CustomStream")
+        custom_line = "Stripping" + line
+        print("line: ", line, " custom_line ", custom_line)
+
+        for stream in streams:
+           for line2 in stream.lines:
+              if line2.name() == custom_line:
+                custom_stream.appendLines([line2])
+                print("line2 name during strip lines loop:", line2.name())
+
+        # Create the actual Stripping configurable
+        filterBadEvents = ProcStatusCheck()
+        sc = StrippingConf(Streams=[custom_stream], MaxCandidates=2000, AcceptBadEvents=False, BadEventSelection=filterBadEvents)
+        # The output is placed directly into Phys, so we only need to
+
+        sr = StrippingReport(Selections = sc.selections())
+        DaVinci().appendToMainSequence( [ sc.sequence() ] )
+
+    relinfo_prefix = ""
+    if isMC:
+       if restrip:
+          theTuple.Inputs = ["/Event/Phys/{}/Particles".format(line)]
+          relinfo_prefix = "/Event/Phys/{}/".format(line)
+       elif isfiltered: # filtering script
+          #stream = "Bu2KLL_NoPID.Strip" 
+          #line_noPID = line.replace("Bu2LLK", "Bu2LLKNoPID")
+          #theTuple.Inputs = ['/Event/{}/Phys/{}/Particles'.format(stream,line_noPID)]
+          #stream = "Leptonic.Strip" 
+          stream = "Bu2LLK_B2LLXBDT_B2XMuMu.Strip" 
+          theTuple.Inputs = ['/Event/{}/Phys/{}/Particles'.format(stream,line)]
+          relinfo_prefix = '/Event/{}/Phys/{}/'.format(stream,line)
+          #DaVinci().RootInTES = "/Event/%s/"%stream #---> only for MDST!!! 
+          #if year == "2012" : # Run1 TM issue
+          #  from Configurables import ToolSvc,Calo2MCTool
+          #  ToolSvc().addTool(Calo2MCTool,name="Calo2MCTool")
+          #  ToolSvc().Calo2MCTool.Hypo2Cluster=True
+       else: #flagged MC
+          theTuple.Inputs = ['/Event/AllStreams/Phys/{}/Particles'.format(line)]
+          relinfo_prefix = '/Event/AllStreams/Phys/{}/'.format(line)
+    else:
+       # local tests fail with this from RpK
+       #DaVinci().RootInTES = "/Event/Leptonic/"  
+       #theInput = "Phys/" + line + "/Particles"     
+       theInput = "/Event/Leptonic/Phys/" + line + "/Particles"     
+       theTuple.Inputs = [ theInput ]
+       relinfo_prefix = '/Event/Leptonic/Phys/{}/'.format(line)
+    print("relinfo prefix is %s"%relinfo_prefix)
+
+    theTuple.ReFitPVs = True 
+    theTuple.ToolList = ["TupleToolRecoStats","TupleToolEventInfo","TupleToolGeometry","TupleToolTrackInfo","TupleToolTrackPosition",]
+    if "mm" in name: theTuple.ToolList += ["TupleToolMuonPid"]
+    # Tools for particle combinations
+    theTuple.Lb.ToolList   = ["TupleToolKinematic"]
+    theTuple.Jpsi.ToolList = ["TupleToolKinematic"]
+    theTuple.Lst.ToolList  = ["TupleToolKinematic"]
+
+    # Tools for the final tracks
+    theTuple.L1.ToolList    =  ["TupleToolKinematic","TupleToolANNPID"] 
+    tt_L1_pid = theTuple.L1.addTupleTool('TupleToolPid')    
+    tt_L1_pid.Verbose = True
+
+    theTuple.L2.ToolList    =  ["TupleToolKinematic","TupleToolANNPID"] 
+    tt_L2_pid = theTuple.L2.addTupleTool('TupleToolPid')    
+    tt_L2_pid.Verbose = True
+
+    theTuple.p.ToolList    =  ["TupleToolKinematic","TupleToolANNPID"] 
+    tt_p_pid = theTuple.p.addTupleTool('TupleToolPid')    
+    tt_p_pid.Verbose = True
+
+    theTuple.k.ToolList    =  ["TupleToolKinematic","TupleToolANNPID"] 
+    tt_k_pid = theTuple.k.addTupleTool('TupleToolPid')    
+    tt_k_pid.Verbose = True
+   
+    theTuple.addTool(TupleToolGeometry,name = 'TupleToolGeometry') 
+    theTuple.TupleToolGeometry.RefitPVs = True
+    if isMC: 
+         theTuple.ToolList += [ "TupleToolMCBackgroundInfo"] 
+         MCTruth=theTuple.addTupleTool("TupleToolMCTruth")
+         MCTruth.addTupleTool("MCTupleToolHierarchy")
+         if isReDecay: theTuple.ToolList += ["MCTupleToolRedecay"]
+
+    LoKiVariables = LoKi__Hybrid__TupleTool('LoKiVariables')
+    LoKiVariables.Variables = {
+                                "ETA"        : "ETA",
+                                "LOKI_PHI"    : "PHI",
+                                "inMuon"      : "switch(INMUON, 1, 0)"
+                              }
+
+    theTuple.addTool(LoKiVariables , name = "LoKiVariables" )
+    theTuple.ToolList   += [ "LoKi::Hybrid::TupleTool/LoKiVariables"]
+
+    # TupleToolTISTOS
+    theTuple.Lb.addTool(TupleToolTISTOS)
+    theTuple.Lb.TupleToolTISTOS.Verbose = True     
+    theTuple.Lb.TupleToolTISTOS.VerboseL0 = True
+    theTuple.Lb.TupleToolTISTOS.VerboseHlt1 = True
+    theTuple.Lb.TupleToolTISTOS.VerboseHlt2 = True
+    theTuple.Lb.TupleToolTISTOS.FillL0 = True
+    theTuple.Lb.TupleToolTISTOS.FillHlt1 = True
+    theTuple.Lb.TupleToolTISTOS.FillHlt2 = True
+    theTuple.Lb.TupleToolTISTOS.OutputLevel = INFO
+    theTuple.Lb.TupleToolTISTOS.TriggerList = trig_B     
+    theTuple.Lb.ToolList +=["TupleToolTISTOS"]     
+
+    theTuple.Lst.addTool(TupleToolTISTOS)
+    theTuple.Lst.TupleToolTISTOS.Verbose = True     
+    theTuple.Lst.TupleToolTISTOS.TriggerList = trig_Lst     
+    theTuple.Lst.ToolList += ["TupleToolTISTOS"]     
+
+    theTuple.Jpsi.addTool(TupleToolTISTOS)
+    theTuple.Jpsi.TupleToolTISTOS.Verbose = True     
+    theTuple.Jpsi.TupleToolTISTOS.TriggerList = trig_Jpsi     
+    theTuple.Jpsi.ToolList +=["TupleToolTISTOS"]     
+
+    #if year == "2012" and isMC: # HLT1TrackAllL0DecisionTOS12MC correction
+    theTuple.L1.addTool(TupleToolTISTOS)
+    theTuple.L1.TupleToolTISTOS.Verbose = True     
+    theTuple.L1.TupleToolTISTOS.TriggerList = trig_Jpsi #["Hlt1TrackAllL0Decision"]    
+    theTuple.L1.ToolList +=["TupleToolTISTOS"]     
+
+    theTuple.L2.addTool(TupleToolTISTOS)
+    theTuple.L2.TupleToolTISTOS.Verbose = True     
+    theTuple.L2.TupleToolTISTOS.TriggerList = trig_Jpsi #["Hlt1TrackAllL0Decision"]    
+    theTuple.L2.ToolList +=["TupleToolTISTOS"]     
+
+    theTuple.p.addTool(TupleToolTISTOS)
+    theTuple.p.TupleToolTISTOS.Verbose = True     
+    theTuple.p.TupleToolTISTOS.TriggerList = trig_Lst #["Hlt1TrackAllL0Decision"]    
+    theTuple.p.ToolList +=["TupleToolTISTOS"]     
+
+    theTuple.k.addTool(TupleToolTISTOS)
+    theTuple.k.TupleToolTISTOS.Verbose = True     
+    theTuple.k.TupleToolTISTOS.TriggerList = trig_Lst #["Hlt1TrackAllL0Decision"]    
+    theTuple.k.ToolList +=["TupleToolTISTOS"]     
+
+    theTuple.ToolList +=["TupleToolL0Data"]     
+
+    # Brem stuff ---------------- 
+    if "ee" in name:
+
+         theTuple.addTupleTool('TupleToolHOP')
+         theTuple.TupleToolHOP.Verbose=True
+
+         theTuple.addTupleTool('TupleToolBremInfo')
+         theTuple.TupleToolBremInfo.Verbose=True
+
+         theTuple.addTupleTool('TupleToolPhotonInfo')
+         theTuple.TupleToolPhotonInfo.Verbose=True
+
+         # track variables: including track momentum before brem recovery and some calo-track matching info
+         LoKiCaloVars = LoKi__Hybrid__TupleTool("LoKiCaloVars")
+         LoKiCaloVars.Preambulo += ['from LoKiTracks.decorators import *']
+
+         LoKiCaloVars.Variables = {
+             "TRACK_PX" : "TRFUN(TrPX)",
+             "TRACK_PY" : "TRFUN(TrPY)",
+             "TRACK_PZ" : "TRFUN(TrPZ)",
+             "TrP"      : "TRFUN(TrP)",
+             "TrQ"      : "TRFUN(TrQ)",
+
+             "LoKi_InAccBrem" : "PPINFO( LHCb.ProtoParticle.InAccBrem, -1 )",
+
+             "CaloTrMatch"       : "PPINFO( LHCb.ProtoParticle.CaloTrMatch,       -999 )",
+             "CaloElectronMatch" : "PPINFO( LHCb.ProtoParticle.CaloElectronMatch, -999 )",
+             "CaloBremMatch"     : "PPINFO( LHCb.ProtoParticle.CaloBremMatch,     -999 )"
+         }
+         
+         theTuple.L1.addTool(LoKiCaloVars , name="LoKiCaloVars" )
+         theTuple.L1.ToolList   += [ "LoKi::Hybrid::TupleTool/LoKiCaloVars"]
+         theTuple.L2.addTool(LoKiCaloVars , name="LoKiCaloVars" )
+         theTuple.L2.ToolList   += [ "LoKi::Hybrid::TupleTool/LoKiCaloVars"]
+       
+         theTuple.L1.addTupleTool( 'TupleToolL0Calo', name = "L1L0ECalo" )
+         theTuple.L1.ToolList += [ "TupleToolL0Calo/L1L0ECalo" ]
+         theTuple.L1.L1L0ECalo.WhichCalo = "ECAL"
+         theTuple.L2.addTupleTool( 'TupleToolL0Calo', name = "L2L0ECalo" )
+         theTuple.L2.ToolList += [ "TupleToolL0Calo/L2L0ECalo" ]
+         theTuple.L2.L2L0ECalo.WhichCalo = "ECAL"
+
+         theTuple.L1.addTupleTool( 'TupleToolL0Calo', name = "L1L0HCalo" )
+         theTuple.L1.ToolList += [ "TupleToolL0Calo/L1L0HCalo" ]
+         theTuple.L1.L1L0HCalo.WhichCalo = "HCAL"
+         theTuple.L2.addTupleTool( 'TupleToolL0Calo', name = "L2L0HCalo" )
+         theTuple.L2.ToolList += [ "TupleToolL0Calo/L2L0HCalo" ]
+         theTuple.L2.L2L0HCalo.WhichCalo = "HCAL"
+
+         theTuple.k.addTupleTool( 'TupleToolL0Calo', name = "KaonL0Calo" )
+         theTuple.k.ToolList += [ "TupleToolL0Calo/KaonL0Calo" ]
+         theTuple.k.KaonL0Calo.WhichCalo = "HCAL"
+         theTuple.p.addTupleTool( 'TupleToolL0Calo', name = "ProtonL0Calo" )
+         theTuple.p.ToolList += [ "TupleToolL0Calo/ProtonL0Calo" ]
+         theTuple.p.ProtonL0Calo.WhichCalo = "HCAL"
+
+    #if not isDST: # isolations 
+    LoKi_Isols =  LoKi__Hybrid__TupleTool("LoKi_Isols")
+    LoKi_Isols.Variables = fillIsoVars(relinfo_prefix,isRun2)    
+    theTuple.Lb.addTool(LoKi_Isols , name = "LoKi_Isols" )
+    theTuple.Lb.ToolList   += [ "LoKi::Hybrid::TupleTool/LoKi_Isols"]
+
+    ##### mass hypothesis substitutions, to study misID backgrounds #####
+    theTuple.Lb.addTupleTool( 'TupleToolSubMass' )
+    theTuple.Lb.ToolList += [ "TupleToolSubMass" ]
+  
+    ## common for mm and ee modes:  
+    theTuple.Lb.TupleToolSubMass.Substitution += ["K+ => pi+"]+["K+ => p+"]
+    theTuple.Lb.TupleToolSubMass.Substitution += ["p+ => K+"]+["p+ => pi+"]
+
+    theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["K-/p+ => p+/K-"]+["K-/p+ => pi+/K-"]+["K-/p+ => p+/pi-"]+["K-/p+ => pi+/pi-"]+["K-/p+ => mu+/mu-"]+["K-/p+ => e+/e-"]
+
+    ##involving specific leptons:
+    if "ee" in name: #electrons
+         theTuple.Lb.TupleToolSubMass.Substitution += ["K+ => e+"]+["p+ => e+"]
+         theTuple.Lb.TupleToolSubMass.Substitution += ["e+ => K+"]+["e+ => pi+"]+["e+ => p+"]+["e+ => mu+"]
+   
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["e+/e- => pi+/pi-"]+["e+/e- => K+/K-"]+["e+/e- => K+/pi-"]+["e+/e- => pi+/K-"]+["e+/e- => mu+/mu-"]
+    
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["K+/e+ => e+/K+"]+["K+/e+ => pi+/pi+"]+["K+/e+ => p+/p+"]
+  
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["p+/e+ => e+/p+"]+["p+/e+ => pi+/pi+"]+["p+/e+ => K+/K+"]+["p+/e+ => K+/pi+"]
+  
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["K+/e- => e+/K-"]+["K+/e- => pi+/pi-"]+["K+/e- => p+/p~-"]+["K+/e- => pi+/K-"]+["K+/e- => pi+/p~-"]+["K+/e- => p+/pi-"]+["K+/e- => p+/K-"]+["K+/e- => mu+/mu-"]
+    
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["p+/e- => e+/p~-"]+["p+/e- => K+/K-"]+["p+/e- => pi+/pi-"]+["p+/e- => pi+/p~-"]+["p+/e- => K+/pi-"]+["p+/e- => K+/p~-"]+["p+/e- => pi+/K-"]+["p+/e- => mu+/mu-"]
+    else: #muons
+         theTuple.Lb.TupleToolSubMass.Substitution += ["K+ => mu+"]+["p+ => mu+"]
+         theTuple.Lb.TupleToolSubMass.Substitution += ["mu+ => K+"]+["mu+ => pi+"]+["mu+ => p+"]+["mu+ => e+"]
+
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["mu+/mu- => pi+/pi-"]+["mu+/mu- => K+/K-"]+["mu+/mu- => K+/pi-"]+["mu+/mu- => pi+/K-"]+["mu+/mu- => e+/e-"]
+    
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["K+/mu+ => mu+/K+"]+["K+/mu+ => pi+/pi+"]+["K+/mu+ => p+/p+"]
+  
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["p+/mu+ => mu+/p+"]+["p+/mu+ => pi+/pi+"]+["p+/mu+ => K+/K+"]+["p+/mu+ => K+/pi+"]
+ 
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["K+/mu- => mu+/K-"]+["K+/mu- => pi+/pi-"]+["K+/mu- => p+/p~-"]+["K+/mu- => pi+/K-"]+["K+/mu- => pi+/p~-"]+["K+/mu- => p+/pi-"]+["K+/mu- => p+/K-"]+["K+/mu- => e+/e-"] 
+    
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["p+/mu- => mu+/p~-"]+["p+/mu- => pi+/pi-"]+["p+/mu- => pi+/p~-"]+["p+/mu- => K+/pi-"]+["p+/mu- => K+/p~-"]+["p+/mu- => pi+/K-"]+["p+/mu- => e+/e-"]
+
+
+    # the signal modes
+    # --- no PV constraint
+    #theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_NoPV"]
+    #dtfSigB_NoPV = TupleToolDecayTreeFitter("SigB_NoPV",
+    #                                        Verbose=True ,
+    #                                        UpdateDaughters=True
+    #                                       )
+    #theTuple.Lb.addTool(dtfSigB_NoPV)
+
+
+    #theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_noPV_Jpsi"]
+    #dtfSigB_noPV_Jpsi = TupleToolDecayTreeFitter("SigB_noPV_Jpsi",
+    #                                        Verbose=True ,
+    #                                        UpdateDaughters=True,
+    #                                        daughtersToConstrain = [ "J/psi(1S)" ] 
+    #                                       )
+    #theTuple.Lb.addTool(dtfSigB_noPV_Jpsi)
+
+    #theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_noPV_psi2S"]
+    #dtfSigB_noPV_psi2S = TupleToolDecayTreeFitter("SigB_noPV_psi2S",
+    #                                        Verbose=True ,
+    #                                        UpdateDaughters=True,
+    #                                        Substitutions = { "[ [Lambda_b0]cc -> ^J/psi(1S) Lambda(1520)0 ]CC" : "psi(2S)" },
+    #                                        daughtersToConstrain = [ "psi(2S)" ] 
+    #                                       )
+    #theTuple.Lb.addTool(dtfSigB_noPV_psi2S)
+
+    # --- with PV constraint (nominal)
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_PV"]
+    dtfSigB_PV = TupleToolDecayTreeFitter("SigB_PV",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            UpdateDaughters=True
+                                           )
+    theTuple.Lb.addTool(dtfSigB_PV)
+
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_PV_Lb_mc"]
+    dtfSigB_PV_Lb_mc = TupleToolDecayTreeFitter("SigB_PV_Lb_mc",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            UpdateDaughters=True,
+                                            daughtersToConstrain = [ "Lambda_b0" ] 
+                                           )
+    theTuple.Lb.addTool(dtfSigB_PV_Lb_mc)
+
+
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_PV_Jpsi"]
+    dtfSigB_PV_Jpsi = TupleToolDecayTreeFitter("SigB_PV_Jpsi",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            UpdateDaughters=True,
+                                            daughtersToConstrain = [ "J/psi(1S)" ] 
+                                           )
+    theTuple.Lb.addTool(dtfSigB_PV_Jpsi)
+
+
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_PV_Jpsi_Lb_mc"]
+    dtfSigB_PV_Jpsi_Lb_mc = TupleToolDecayTreeFitter("SigB_PV_Jpsi_Lb_mc",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            UpdateDaughters=True,
+                                            daughtersToConstrain = [ "J/psi(1S)", "Lambda_b0" ] 
+                                           )
+    theTuple.Lb.addTool(dtfSigB_PV_Jpsi_Lb_mc)
+
+
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_PV_psi2S"]
+    dtfSigB_PV_psi2S = TupleToolDecayTreeFitter("SigB_PV_psi2S",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            UpdateDaughters=True,
+                                            Substitutions = { "Lambda_b0 -> ^J/psi(1S) Lambda(1520)0" : "psi(2S)",
+                                                              "Lambda_b~0 -> ^J/psi(1S) Lambda(1520)~0" : "psi(2S)" 
+                                                            },
+                                            daughtersToConstrain = [ "psi(2S)" ] 
+                                           )
+    theTuple.Lb.addTool(dtfSigB_PV_psi2S)
+
+
+    """
+    # the B0 -> K pi l+ l- bkgd 
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SubBd_PV"]
+    dtfSubBd_PV = TupleToolDecayTreeFitter("SubBd_PV",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            Substitutions = { "Lambda_b0 -> J/psi(1S)  (Lambda(1520)0 -> ^p+ X-)"  : "pi+",
+                                                              "Lambda_b~0 -> J/psi(1S)    (Lambda(1520)~0 -> ^p~- X+)" : "pi-" 
+                                                            }
+                                           )
+    theTuple.Lb.addTool(dtfSubBd_PV)
+
+    # the B0s -> K K l+ l- bkgd 
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SubBs_PV"]
+    dtfSubBs_PV = TupleToolDecayTreeFitter("SubBs_PV",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            Substitutions = { "Lambda_b0 -> J/psi(1S)   (Lambda(1520)0 -> ^p+ X-)"  : "K+",
+                                                              "Lambda_b~0 -> J/psi(1S)  (Lambda(1520)~0 -> ^p~- X+)" : "K-" 
+                                                            }
+                                           )
+    theTuple.Lb.addTool(dtfSubBs_PV)
+    """
+
+
+    return theTuple
+
+
+###################################################
+# automatically_configure
+# Use the CondDB to get the right database tags for data
+# from Configurables import CondDB
+# if not isMC: CondDB(LatestGlobalTagByDataType = year )
+###################################################
+
+dv = DaVinci(
+     #InputType     = inputType, #automatically_configure # BD localtest
+     #DataType      =  year  , # automatically_configure # BD localtest
+     EvtMax        =  evtmax  , # -1 for all files
+     TupleFile     = 'DTT.root',
+     PrintFreq     = 100000,
+     Simulation    = isMC,
+     Lumi          = False,
+)
+
+#MessageSvc().OutputLevel = DEBUG
+
+"""
+# Matt Needham(?)'s momentum scale calibration for data
+from Configurables import TrackScaleState
+if not isMC:
+   scaler            = TrackScaleState('scaler')
+   #scaler.Input = "/Event/Leptonic/Rec/Track/Best"
+   dv.UserAlgorithms =  [ scaler ]
+"""
+
+# smear proto particles and their associated tracks for MC
+# https://gitlab.cern.ch/LHCb-RD/ewp-bd2ksteeangular-legacy/-/blob/master/samples_and_selection/DaVinci_scripts/options_Bu2LLK_EE.py#L121
+#from Configurables import TrackSmearState
+# removed 10th march 2022
+#if isMC:
+#   smear = TrackSmearState("TrackSmearState")
+#   dv.UserAlgorithms =  [ smear ]
+
+# scaler/smear must be executed before the analysis algorithms
+isDST = True
+if inputType == "MDST": isDST = False
+
+if doMuon:
+   #print(isRun2, isMC, isDST) # BD localtest
+   Lbtuple_mm = addNtupleToDaVinci("Lb2pkmm", isRun2, isMC, isDST)
+   dv.UserAlgorithms += [Lbtuple_mm]
+
+if doElectron: 
+   #print("doElectron: ",doElectron) # BD localtest
+   Lbtuple_ee = addNtupleToDaVinci("Lb2pkee", isRun2, isMC, isDST)
+   dv.UserAlgorithms += [Lbtuple_ee]
+
+import glob
+#evt+std://LHCb/Collision12/90000000/Beam4000GeV-VeloClosed-MagUp/Real Data/Reco14/Stripping21r0p2/LEPTONIC.MDST
+#import glob
+#flist    = glob.glob("/eos/lhcb/grid/prod//lhcb/LHCb/Collision12/LEPTONIC.MDST/00095380/0000/0*.leptonic.mdst")
+#dv.Input = flist[:1]
+
+# evt+std://LHCb/Collision17/90000000/Beam6500GeV-VeloClosed-MagDown/Real Data/Reco17/Stripping29r2p2/LEPTONIC.MDST 
+#dv.Input = ["/eos/lhcb/grid/prod//lhcb/LHCb/Collision17/LEPTONIC.MDST/00141012/0000/00141012_00000060_1.leptonic.mdst"]
+#flist    = glob.glob("/eos/lhcb/grid/prod//lhcb/LHCb/Collision17/LEPTONIC.MDST/00141012/0000/0*.leptonic.mdst")
+#dv.Input = flist[:1]
+
+#MC_15144001_Lb_JpsipKmmphspDecProdCut_ReStrip
+# evt+std://MC/2018/15144001/Beam6500GeV-2018-MagDown-Nu1.6-25ns-Pythia8/Sim09f/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34NoPrescalingFlagged/ALLSTREAMS.DST
+#dv.Input = ["/eos/lhcb/grid/prod//lhcb/MC/2018/ALLSTREAMS.DST/00086923/0000/00086923_00000406_7.AllStreams.dst"]
+
+# MC_15114401_Lb_pKstmmKpi0phspDecProdCut
+# evt+std://MC/2012/15114401/Beam4000GeV-2012-MagUp-Nu2.5-Pythia8/Sim09k/Trig0x409f0045/Reco14c/Stripping21r0p2NoPrescalingFlagged/ALLSTREAMS.DST
+#dv.Input = ["/eos/lhcb/grid/prod//lhcb/MC/2012/ALLSTREAMS.DST/00149395/0000/00149395_00000110_5.AllStreams.dst"]
+
+#MC_15144059_Lb_JpsipKmmXLL2_ReStrip
+# evt+std://MC/2012/15144059/Beam4000GeV-2012-MagDown-Nu2.5-Pythia8/Sim09h/Trig0x409f0045/Reco14c/Stripping21NoPrescalingFlagged/ALLSTREAMS.DST
+#dv.Input = ["/eos/lhcb/grid/prod//lhcb//MC/2012/ALLSTREAMS.DST/00095494/0000/00095494_00000001_5.AllStreams.dst"]
+
+#evt+std://MC/2012/15154011/Beam4000GeV-2012-MagUp-Nu2.5-Pythia8/Sim09l-ReDecay01/Trig0x409f0045/Reco14c/Stripping21r0p2Filtered/BU2LLK_B2LLXBDT_B2XMUMU.STRIP.MDST
+# Lb2pKee
+#flist    = glob.glob("/eos/lhcb/grid/prod//lhcb/MC/2017/BU2LLK_B2LLXBDT_B2XMUMU.STRIP.MDST/00158235/0000/00158235_00000001_1.bu2llk_b2llxbdt_b2xmumu.strip.mdst")
+# Lb2psi2SpK,ee
+#flist    = glob.glob("/eos/lhcb/grid/prod//lhcb/MC/2012/BU2LLK_B2LLXBDT_B2XMUMU.STRIP.MDST/00157640/0000/00157640_000000*_1.bu2llk_b2llxbdt_b2xmumu.strip.mdst")
+
+#15114012 (Lb_pKmumu=XLL,DecProdCut) Run1 dst
+#flist    = glob.glob("/eos/lhcb/grid/prod//lhcb/MC/2012/BU2LLK_B2LLXBDT_B2XMUMU.STRIP.DST/00166104/0000/00166104_00000*bu2llk_b2llxbdt_b2xmumu.strip.dst")
+#dv.Input = flist[:1]
+
+
+
+
+
diff --git a/Lb2pKll/dv_pKll_mdst.py b/Lb2pKll/dv_pKll_mdst.py
new file mode 100755
index 0000000000000000000000000000000000000000..1f84110e313c9e507bad75479b49982e90da902d
--- /dev/null
+++ b/Lb2pKll/dv_pKll_mdst.py
@@ -0,0 +1,570 @@
+# References:
+# https://twiki.cern.ch/twiki/pub/LHCbPhysics/Lst2pkll/LHCb-ANA-2019-007_v3.1.pdf
+# https://gitlab.cern.ch/LHCb-RD/ewp-rpk-fullr1r2/-/blob/master/DaVinci/LeptonU-options-NOPID.py
+# https://gitlab.cern.ch/LHCb-RD/ewp-rks-rkstplus/-/blob/master/Produce_Ntuples/Data/Bu2LLK_DaVinci.py
+# https://gitlab.cern.ch/LHCb-RD/ewp-bd2ksteeangular-legacy/-/blob/master/samples_and_selection/DaVinci_scripts/options_Bu2LLK_EE.py
+
+import sys
+sys.path.append("./")
+#from steering_files.MC_2018_MD_ReStrip0_Filt1_mumu import * # BD localtest
+#from steering_files.MC_2017_MU_ReStrip0_Filt1_ee import * # BD localtest
+#from steering_files.MC_2012_MD_ReStrip0_Filt1_mumu import * # BD localtest
+#from steering_files.MC_2012_MD_ReStrip1_Filt0_mumu import * # BD localtest
+#from steering_files.Data_2012 import * # BD localtest
+#from steering_files.Data_2017 import * # BD localtest
+
+import os
+isRun2      = (os.getenv("ISRUN2", 'False') == 'True')     
+isMC        = (os.getenv("ISMC", 'False') == 'True') 
+year        = os.getenv("YEAR")   
+inputType   = os.getenv("INPUTTYPE")  
+isReDecay   = (os.getenv("ISREDECAY", 'False') == 'True')  
+isfiltered  = (os.getenv("ISFILTERED", 'False') == 'True')  
+magpos      = os.getenv("MAGPOS")    
+restrip     = (os.getenv("RESTRIP", 'False') == 'True') 
+doMuon      = (os.getenv("DOMUON", 'False') == 'True')  
+doElectron  = (os.getenv("DOELECTRON", 'False') == 'True') 
+
+
+from Lb2pKll.helpers.helper import *
+#from helpers.helper import *  # BD localtest
+from LHCbKernel.Configuration import *
+from Configurables import DaVinci, TupleToolDecayTreeFitter, DecayTreeTuple, TupleToolDecay, TupleToolTrigger, TupleToolRecoStats, TupleToolTISTOS, TupleToolGeometry, TupleToolTrackInfo, TupleToolMCTruth, TupleToolL0Data, SubstitutePID
+if isReDecay: from Configurables import MCTupleToolRedecay
+
+from Configurables import DstConf, CaloDstUnPackConf, LoKi__Hybrid__TupleTool
+from Configurables import EventNodeKiller, ProcStatusCheck
+from Configurables import StrippingReport, TimingAuditor, SequencerTimerTool
+
+from StrippingConf.Configuration import StrippingConf, StrippingStream
+from StrippingSettings.Utils import strippingConfiguration
+from StrippingArchive.Utils import buildStreams
+from StrippingArchive import strippingArchive
+
+from Gaudi.Configuration import *
+from DecayTreeTuple.Configuration import *
+
+from Configurables import CombineParticles
+from Configurables import TupleToolTrigger, SubstitutePID
+from Configurables import FilterDesktop
+from PhysSelPython.Wrappers import Selection, SelectionSequence, DataOnDemand, MergedSelection
+
+##################################################################
+# Initialize some settings independent of LHCb tools
+evtmax         =  -1 # change this for localtest. has to be -1 for AP
+###################################################################
+
+# RpK
+linename_e_Run1   = "Bu2LLK_eeLine2"
+linename_mu_Run1  = "Bu2LLK_mmLine"
+# new lines
+linename_e_Run2   = "Bu2LLK_ee_Lb2pkeeLooseLine"
+linename_mu_Run2  = "Bu2LLK_mm_Lb2pkmmLooseLine"
+
+#######################################################################
+
+# applied on the B head
+trig_L0TIS        = ['L0GlobalDecision','L0HadronDecision',"L0MuonDecision", "L0ElectronDecision"] # TIS-ed on
+trig_Run1_B_mu    = trig_L0TIS + ['Hlt1TrackAllL0Decision', 'Hlt1TrackMuonDecision', 'Hlt1TrackDiMuonDecision', 'Hlt2DiMuonDetachedDecision', 'Hlt2Topo2BodyBBDTDecision', 'Hlt2Topo3BodyBBDTDecision', 'Hlt2Topo4BodyBBDTDecision', 'Hlt2TopoMu2BodyBBDTDecision', 'Hlt2TopoMu3BodyBBDTDecision', 'Hlt2TopoMu4BodyBBDTDecision']
+trig_Run2_B_mu    = trig_L0TIS + ['Hlt1TrackMuonDecision','Hlt1TrackMVADecision', 'Hlt1TwoTrackMVADecision', 'Hlt1TrackMuonMVADecision', 'Hlt1DiMuonHighMassDecision', 'Hlt2DiMuonDetachedDecision', 'Hlt2Topo2BodyDecision','Hlt2Topo3BodyDecision','Hlt2Topo4BodyDecision','Hlt2TopoMu2BodyDecision','Hlt2TopoMu3BodyDecision','Hlt2TopoMu4BodyDecision','Hlt2TopoMuMu2BodyDecision','Hlt2TopoMuMu3BodyDecision','Hlt2TopoMuMu4BodyDecision']
+trig_Run1_B_e     = trig_L0TIS + ['Hlt1TrackAllL0Decision', 'Hlt1TrackAllL0TightDecision','Hlt2Topo2BodyBBDTDecision', 'Hlt2Topo3BodyBBDTDecision', 'Hlt2Topo4BodyBBDTDecision', 'Hlt2TopoE2BodyBBDTDecision', 'Hlt2TopoE3BodyBBDTDecision', 'Hlt2TopoE4BodyBBDTDecision']
+trig_Run2_B_e     = trig_L0TIS + ['Hlt1TrackMVADecision', 'Hlt1TwoTrackMVADecision', 'Hlt2Topo2BodyDecision', 'Hlt2Topo3BodyDecision', 'Hlt2Topo4BodyDecision', 'Hlt2TopoE2BodyDecision', 'Hlt2TopoE3BodyDecision', 'Hlt2TopoE4BodyDecision', 'Hlt2TopoEE2BodyDecision', 'Hlt2TopoEE3BodyDecision', 'Hlt2TopoEE4BodyDecision']
+
+# some Hlt1/2 muon lines for the dimuon
+trig_Run1_Jpsi_mu = ["L0MuonDecision","L0DiMuonDecision", "Hlt1DiMuonLowMassDecision","Hlt1DiMuonHighMassDecision","Hlt1TrackMuonDecision","Hlt2SingleMuonDecision","Hlt2DiMuonDetachedDecision"]
+trig_Run2_Jpsi_mu = ["L0MuonDecision","L0DiMuonDecision", "Hlt1TrackMuonMVADecision","Hlt1TrackMuonDecision","Hlt1SingleMuonHighPTDecision","Hlt1DiMuonLowMassDecision","Hlt1DiMuonHighMassDecision","Hlt2DiMuonDetachedDecision","Hlt2SingleMuonDecision"]
+trig_e            = ["L0ElectronDecision","L0ElectronHiDecision", "L0PhotonDecision", "L0PhotonHiDecision", "Hlt1SingleElectronNoIPDecision"]
+
+# applied on the pK
+trig_Lst          = ['L0HadronDecision']
+
+
+def addNtupleToDaVinci(name, isRun2, isMC, isDST):
+    trig_B = trig_Run1_B_mu
+    trig_Jpsi = trig_Run1_Jpsi_mu
+    line   = linename_mu_Run1
+    if "ee" in name and isRun2 == 0: 
+       trig_B = trig_Run1_B_e 
+       line   = linename_e_Run1
+    elif "mm" in name and isRun2 == 1:
+       trig_B = trig_Run2_B_mu 
+       trig_Jpsi = trig_Run2_Jpsi_mu
+       line   = linename_mu_Run2
+    elif "ee" in name and isRun2 == 1:
+       trig_B = trig_Run2_B_e 
+       line   = linename_e_Run2
+    theTuple = DecayTreeTuple(name)
+    if "ee" in name:
+       trig_Jpsi = trig_e
+       theTuple.Decay = "[Lambda_b0 -> ^(J/psi(1S) -> ^e+ ^e-) ^(Lambda(1520)0 -> ^p+ ^K-)]CC"
+       theTuple.addBranches ({
+         "Lb"  : "[Lambda_b0->  (J/psi(1S) -> e+ e-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "L1"  : "[Lambda_b0->  (J/psi(1S) -> ^e+ e-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "L2"  : "[Lambda_b0->  (J/psi(1S) -> e+ ^e-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "Jpsi": "[Lambda_b0->  ^(J/psi(1S) -> e+ e-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "p"   : "[Lambda_b0->  (J/psi(1S) -> e+ e-)  ( Lambda(1520)0 -> ^p+  K-) ]CC",
+         "k"   : "[Lambda_b0->  (J/psi(1S) -> e+ e-)  ( Lambda(1520)0 -> p+  ^K-) ]CC",
+         "Lst" : "[Lambda_b0->  (J/psi(1S) -> e+ e-)  ^( Lambda(1520)0 -> p+  K-)]CC"
+       })
+    else:
+        theTuple.Decay = "[Lambda_b0 -> ^(J/psi(1S) -> ^mu+ ^mu-) ^(Lambda(1520)0 -> ^p+ ^K-)]CC" 
+        theTuple.addBranches ({
+         "Lb"  : "[Lambda_b0->  (J/psi(1S) -> mu+ mu-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "L1"  : "[Lambda_b0->  (J/psi(1S) -> ^mu+ mu-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "L2"  : "[Lambda_b0->  (J/psi(1S) -> mu+ ^mu-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "Jpsi": "[Lambda_b0->  ^(J/psi(1S) -> mu+ mu-)  ( Lambda(1520)0 -> p+  K-) ]CC",
+         "p"   : "[Lambda_b0->  (J/psi(1S) -> mu+ mu-)  ( Lambda(1520)0 -> ^p+  K-) ]CC",
+         "k"   : "[Lambda_b0->  (J/psi(1S) -> mu+ mu-)  ( Lambda(1520)0 -> p+  ^K-) ]CC",
+         "Lst" : "[Lambda_b0->  (J/psi(1S) -> mu+ mu-)  ^( Lambda(1520)0 -> p+  K-)]CC"
+        })
+    if isMC and restrip: 
+        if year == "2011" : 
+              strip = "stripping21r1p2"
+        elif year == "2012" : 
+              strip = "stripping21r0p2"
+        elif year == "2016" :
+              strip = "stripping28r2p1"
+        elif year == "2017" :
+              strip = "stripping29r2p2"
+        elif year == "2018" :
+              strip = "stripping34r0p2"
+        # The event node killer should be run before
+        event_node_killer = EventNodeKiller("StripKiller")
+        event_node_killer.Nodes = ["/Event/AllStreams", "/Event/Strip"]
+
+        DaVinci().appendToMainSequence( [ event_node_killer ] )
+        #print(year, strip)
+        streams = buildStreams(stripping=strippingConfiguration(strip), archive=strippingArchive(strip))
+
+        custom_stream = StrippingStream("CustomStream")
+        custom_line = "Stripping" + line
+        print("line: ", line, " custom_line ", custom_line)
+
+        for stream in streams:
+           for line2 in stream.lines:
+              if line2.name() == custom_line:
+                custom_stream.appendLines([line2])
+                print("line2 name during strip lines loop:", line2.name())
+
+        # Create the actual Stripping configurable
+        filterBadEvents = ProcStatusCheck()
+        sc = StrippingConf(Streams=[custom_stream], MaxCandidates=2000, AcceptBadEvents=False, BadEventSelection=filterBadEvents)
+        # The output is placed directly into Phys, so we only need to
+
+        sr = StrippingReport(Selections = sc.selections())
+        DaVinci().appendToMainSequence( [ sc.sequence() ] )
+
+    relinfo_prefix = ""
+    if isMC:
+       if restrip:
+          theTuple.Inputs = ["/Event/Phys/{}/Particles".format(line)]
+          relinfo_prefix = "/Event/Phys/{}/".format(line)
+       elif isfiltered: # filtering script
+          #stream = "Bu2KLL_NoPID.Strip" 
+          #line_noPID = line.replace("Bu2LLK", "Bu2LLKNoPID")
+          #theTuple.Inputs = ['/Event/{}/Phys/{}/Particles'.format(stream,line_noPID)]
+          #stream = "Leptonic.Strip" 
+          stream = "Bu2LLK_B2LLXBDT_B2XMuMu.Strip" 
+          theTuple.Inputs = ['/Event/{}/Phys/{}/Particles'.format(stream,line)]
+          relinfo_prefix = '/Event/{}/Phys/{}/'.format(stream,line)
+          DaVinci().RootInTES = "/Event/%s/"%stream #---> only for MDST!!! 
+          #if year == "2012" : # Run1 TM issue
+          #  from Configurables import ToolSvc,Calo2MCTool
+          #  ToolSvc().addTool(Calo2MCTool,name="Calo2MCTool")
+          #  ToolSvc().Calo2MCTool.Hypo2Cluster=True
+       else: #flagged MC
+          theTuple.Inputs = ['/Event/AllStreams/Phys/{}/Particles'.format(line)]
+          relinfo_prefix = '/Event/AllStreams/Phys/{}/'.format(line)
+    else:
+       # local tests fail with this from RpK
+       #DaVinci().RootInTES = "/Event/Leptonic/"  
+       #theInput = "Phys/" + line + "/Particles"     
+       theInput = "/Event/Leptonic/Phys/" + line + "/Particles"     
+       theTuple.Inputs = [ theInput ]
+       relinfo_prefix = '/Event/Leptonic/Phys/{}/'.format(line)
+    print("relinfo prefix is %s"%relinfo_prefix)
+
+    theTuple.ReFitPVs = True 
+    theTuple.ToolList = ["TupleToolRecoStats","TupleToolEventInfo","TupleToolGeometry","TupleToolTrackInfo","TupleToolTrackPosition",]
+    if "mm" in name: theTuple.ToolList += ["TupleToolMuonPid"]
+    # Tools for particle combinations
+    theTuple.Lb.ToolList   = ["TupleToolKinematic"]
+    theTuple.Jpsi.ToolList = ["TupleToolKinematic"]
+    theTuple.Lst.ToolList  = ["TupleToolKinematic"]
+
+    # Tools for the final tracks
+    theTuple.L1.ToolList    =  ["TupleToolKinematic","TupleToolANNPID"] 
+    tt_L1_pid = theTuple.L1.addTupleTool('TupleToolPid')    
+    tt_L1_pid.Verbose = True
+
+    theTuple.L2.ToolList    =  ["TupleToolKinematic","TupleToolANNPID"] 
+    tt_L2_pid = theTuple.L2.addTupleTool('TupleToolPid')    
+    tt_L2_pid.Verbose = True
+
+    theTuple.p.ToolList    =  ["TupleToolKinematic","TupleToolANNPID"] 
+    tt_p_pid = theTuple.p.addTupleTool('TupleToolPid')    
+    tt_p_pid.Verbose = True
+
+    theTuple.k.ToolList    =  ["TupleToolKinematic","TupleToolANNPID"] 
+    tt_k_pid = theTuple.k.addTupleTool('TupleToolPid')    
+    tt_k_pid.Verbose = True
+   
+    theTuple.addTool(TupleToolGeometry,name = 'TupleToolGeometry') 
+    theTuple.TupleToolGeometry.RefitPVs = True
+    if isMC: 
+         theTuple.ToolList += [ "TupleToolMCBackgroundInfo"] 
+         MCTruth=theTuple.addTupleTool("TupleToolMCTruth")
+         MCTruth.addTupleTool("MCTupleToolHierarchy")
+         if isReDecay: theTuple.ToolList += ["MCTupleToolRedecay"]
+
+    LoKiVariables = LoKi__Hybrid__TupleTool('LoKiVariables')
+    LoKiVariables.Variables = {
+                                "ETA"        : "ETA",
+                                "LOKI_PHI"    : "PHI",
+                                "inMuon"      : "switch(INMUON, 1, 0)"
+                              }
+
+    theTuple.addTool(LoKiVariables , name = "LoKiVariables" )
+    theTuple.ToolList   += [ "LoKi::Hybrid::TupleTool/LoKiVariables"]
+
+    # TupleToolTISTOS
+    theTuple.Lb.addTool(TupleToolTISTOS)
+    theTuple.Lb.TupleToolTISTOS.Verbose = True     
+    theTuple.Lb.TupleToolTISTOS.VerboseL0 = True
+    theTuple.Lb.TupleToolTISTOS.VerboseHlt1 = True
+    theTuple.Lb.TupleToolTISTOS.VerboseHlt2 = True
+    theTuple.Lb.TupleToolTISTOS.FillL0 = True
+    theTuple.Lb.TupleToolTISTOS.FillHlt1 = True
+    theTuple.Lb.TupleToolTISTOS.FillHlt2 = True
+    theTuple.Lb.TupleToolTISTOS.OutputLevel = INFO
+    theTuple.Lb.TupleToolTISTOS.TriggerList = trig_B     
+    theTuple.Lb.ToolList +=["TupleToolTISTOS"]     
+
+    theTuple.Lst.addTool(TupleToolTISTOS)
+    theTuple.Lst.TupleToolTISTOS.Verbose = True     
+    theTuple.Lst.TupleToolTISTOS.TriggerList = trig_Lst     
+    theTuple.Lst.ToolList += ["TupleToolTISTOS"]     
+
+    theTuple.Jpsi.addTool(TupleToolTISTOS)
+    theTuple.Jpsi.TupleToolTISTOS.Verbose = True     
+    theTuple.Jpsi.TupleToolTISTOS.TriggerList = trig_Jpsi     
+    theTuple.Jpsi.ToolList +=["TupleToolTISTOS"]     
+
+    #if year == "2012" and isMC: # HLT1TrackAllL0DecisionTOS12MC correction
+    theTuple.L1.addTool(TupleToolTISTOS)
+    theTuple.L1.TupleToolTISTOS.Verbose = True     
+    theTuple.L1.TupleToolTISTOS.TriggerList = trig_Jpsi #["Hlt1TrackAllL0Decision"]    
+    theTuple.L1.ToolList +=["TupleToolTISTOS"]     
+
+    theTuple.L2.addTool(TupleToolTISTOS)
+    theTuple.L2.TupleToolTISTOS.Verbose = True     
+    theTuple.L2.TupleToolTISTOS.TriggerList = trig_Jpsi #["Hlt1TrackAllL0Decision"]    
+    theTuple.L2.ToolList +=["TupleToolTISTOS"]     
+
+    theTuple.p.addTool(TupleToolTISTOS)
+    theTuple.p.TupleToolTISTOS.Verbose = True     
+    theTuple.p.TupleToolTISTOS.TriggerList = trig_Lst #["Hlt1TrackAllL0Decision"]    
+    theTuple.p.ToolList +=["TupleToolTISTOS"]     
+
+    theTuple.k.addTool(TupleToolTISTOS)
+    theTuple.k.TupleToolTISTOS.Verbose = True     
+    theTuple.k.TupleToolTISTOS.TriggerList = trig_Lst #["Hlt1TrackAllL0Decision"]    
+    theTuple.k.ToolList +=["TupleToolTISTOS"]     
+
+    theTuple.ToolList +=["TupleToolL0Data"]     
+
+    # Brem stuff ---------------- 
+    if "ee" in name:
+
+         theTuple.addTupleTool('TupleToolHOP')
+         theTuple.TupleToolHOP.Verbose=True
+
+         theTuple.addTupleTool('TupleToolBremInfo')
+         theTuple.TupleToolBremInfo.Verbose=True
+
+         theTuple.addTupleTool('TupleToolPhotonInfo')
+         theTuple.TupleToolPhotonInfo.Verbose=True
+
+         # track variables: including track momentum before brem recovery and some calo-track matching info
+         LoKiCaloVars = LoKi__Hybrid__TupleTool("LoKiCaloVars")
+         LoKiCaloVars.Preambulo += ['from LoKiTracks.decorators import *']
+
+         LoKiCaloVars.Variables = {
+             "TRACK_PX" : "TRFUN(TrPX)",
+             "TRACK_PY" : "TRFUN(TrPY)",
+             "TRACK_PZ" : "TRFUN(TrPZ)",
+             "TrP"      : "TRFUN(TrP)",
+             "TrQ"      : "TRFUN(TrQ)",
+
+             "LoKi_InAccBrem" : "PPINFO( LHCb.ProtoParticle.InAccBrem, -1 )",
+
+             "CaloTrMatch"       : "PPINFO( LHCb.ProtoParticle.CaloTrMatch,       -999 )",
+             "CaloElectronMatch" : "PPINFO( LHCb.ProtoParticle.CaloElectronMatch, -999 )",
+             "CaloBremMatch"     : "PPINFO( LHCb.ProtoParticle.CaloBremMatch,     -999 )"
+         }
+         
+         theTuple.L1.addTool(LoKiCaloVars , name="LoKiCaloVars" )
+         theTuple.L1.ToolList   += [ "LoKi::Hybrid::TupleTool/LoKiCaloVars"]
+         theTuple.L2.addTool(LoKiCaloVars , name="LoKiCaloVars" )
+         theTuple.L2.ToolList   += [ "LoKi::Hybrid::TupleTool/LoKiCaloVars"]
+       
+         theTuple.L1.addTupleTool( 'TupleToolL0Calo', name = "L1L0ECalo" )
+         theTuple.L1.ToolList += [ "TupleToolL0Calo/L1L0ECalo" ]
+         theTuple.L1.L1L0ECalo.WhichCalo = "ECAL"
+         theTuple.L2.addTupleTool( 'TupleToolL0Calo', name = "L2L0ECalo" )
+         theTuple.L2.ToolList += [ "TupleToolL0Calo/L2L0ECalo" ]
+         theTuple.L2.L2L0ECalo.WhichCalo = "ECAL"
+
+         theTuple.L1.addTupleTool( 'TupleToolL0Calo', name = "L1L0HCalo" )
+         theTuple.L1.ToolList += [ "TupleToolL0Calo/L1L0HCalo" ]
+         theTuple.L1.L1L0HCalo.WhichCalo = "HCAL"
+         theTuple.L2.addTupleTool( 'TupleToolL0Calo', name = "L2L0HCalo" )
+         theTuple.L2.ToolList += [ "TupleToolL0Calo/L2L0HCalo" ]
+         theTuple.L2.L2L0HCalo.WhichCalo = "HCAL"
+
+         theTuple.k.addTupleTool( 'TupleToolL0Calo', name = "KaonL0Calo" )
+         theTuple.k.ToolList += [ "TupleToolL0Calo/KaonL0Calo" ]
+         theTuple.k.KaonL0Calo.WhichCalo = "HCAL"
+         theTuple.p.addTupleTool( 'TupleToolL0Calo', name = "ProtonL0Calo" )
+         theTuple.p.ToolList += [ "TupleToolL0Calo/ProtonL0Calo" ]
+         theTuple.p.ProtonL0Calo.WhichCalo = "HCAL"
+
+    #if not isDST: # isolations 
+    LoKi_Isols =  LoKi__Hybrid__TupleTool("LoKi_Isols")
+    LoKi_Isols.Variables = fillIsoVars(relinfo_prefix,isRun2)    
+    theTuple.Lb.addTool(LoKi_Isols , name = "LoKi_Isols" )
+    theTuple.Lb.ToolList   += [ "LoKi::Hybrid::TupleTool/LoKi_Isols"]
+
+    ##### mass hypothesis substitutions, to study misID backgrounds #####
+    theTuple.Lb.addTupleTool( 'TupleToolSubMass' )
+    theTuple.Lb.ToolList += [ "TupleToolSubMass" ]
+  
+    ## common for mm and ee modes:  
+    theTuple.Lb.TupleToolSubMass.Substitution += ["K+ => pi+"]+["K+ => p+"]
+    theTuple.Lb.TupleToolSubMass.Substitution += ["p+ => K+"]+["p+ => pi+"]
+
+    theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["K-/p+ => p+/K-"]+["K-/p+ => pi+/K-"]+["K-/p+ => p+/pi-"]+["K-/p+ => pi+/pi-"]+["K-/p+ => mu+/mu-"]+["K-/p+ => e+/e-"]
+
+    ##involving specific leptons:
+    if "ee" in name: #electrons
+         theTuple.Lb.TupleToolSubMass.Substitution += ["K+ => e+"]+["p+ => e+"]
+         theTuple.Lb.TupleToolSubMass.Substitution += ["e+ => K+"]+["e+ => pi+"]+["e+ => p+"]+["e+ => mu+"]
+   
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["e+/e- => pi+/pi-"]+["e+/e- => K+/K-"]+["e+/e- => K+/pi-"]+["e+/e- => pi+/K-"]+["e+/e- => mu+/mu-"]
+    
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["K+/e+ => e+/K+"]+["K+/e+ => pi+/pi+"]+["K+/e+ => p+/p+"]
+  
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["p+/e+ => e+/p+"]+["p+/e+ => pi+/pi+"]+["p+/e+ => K+/K+"]+["p+/e+ => K+/pi+"]
+  
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["K+/e- => e+/K-"]+["K+/e- => pi+/pi-"]+["K+/e- => p+/p~-"]+["K+/e- => pi+/K-"]+["K+/e- => pi+/p~-"]+["K+/e- => p+/pi-"]+["K+/e- => p+/K-"]+["K+/e- => mu+/mu-"]
+    
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["p+/e- => e+/p~-"]+["p+/e- => K+/K-"]+["p+/e- => pi+/pi-"]+["p+/e- => pi+/p~-"]+["p+/e- => K+/pi-"]+["p+/e- => K+/p~-"]+["p+/e- => pi+/K-"]+["p+/e- => mu+/mu-"]
+    else: #muons
+         theTuple.Lb.TupleToolSubMass.Substitution += ["K+ => mu+"]+["p+ => mu+"]
+         theTuple.Lb.TupleToolSubMass.Substitution += ["mu+ => K+"]+["mu+ => pi+"]+["mu+ => p+"]+["mu+ => e+"]
+
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["mu+/mu- => pi+/pi-"]+["mu+/mu- => K+/K-"]+["mu+/mu- => K+/pi-"]+["mu+/mu- => pi+/K-"]+["mu+/mu- => e+/e-"]
+    
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["K+/mu+ => mu+/K+"]+["K+/mu+ => pi+/pi+"]+["K+/mu+ => p+/p+"]
+  
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["p+/mu+ => mu+/p+"]+["p+/mu+ => pi+/pi+"]+["p+/mu+ => K+/K+"]+["p+/mu+ => K+/pi+"]
+ 
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["K+/mu- => mu+/K-"]+["K+/mu- => pi+/pi-"]+["K+/mu- => p+/p~-"]+["K+/mu- => pi+/K-"]+["K+/mu- => pi+/p~-"]+["K+/mu- => p+/pi-"]+["K+/mu- => p+/K-"]+["K+/mu- => e+/e-"] 
+    
+         theTuple.Lb.TupleToolSubMass.DoubleSubstitution += ["p+/mu- => mu+/p~-"]+["p+/mu- => pi+/pi-"]+["p+/mu- => pi+/p~-"]+["p+/mu- => K+/pi-"]+["p+/mu- => K+/p~-"]+["p+/mu- => pi+/K-"]+["p+/mu- => e+/e-"]
+
+
+    # the signal modes
+    # --- no PV constraint
+    #theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_NoPV"]
+    #dtfSigB_NoPV = TupleToolDecayTreeFitter("SigB_NoPV",
+    #                                        Verbose=True ,
+    #                                        UpdateDaughters=True
+    #                                       )
+    #theTuple.Lb.addTool(dtfSigB_NoPV)
+
+
+    #theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_noPV_Jpsi"]
+    #dtfSigB_noPV_Jpsi = TupleToolDecayTreeFitter("SigB_noPV_Jpsi",
+    #                                        Verbose=True ,
+    #                                        UpdateDaughters=True,
+    #                                        daughtersToConstrain = [ "J/psi(1S)" ] 
+    #                                       )
+    #theTuple.Lb.addTool(dtfSigB_noPV_Jpsi)
+
+    #theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_noPV_psi2S"]
+    #dtfSigB_noPV_psi2S = TupleToolDecayTreeFitter("SigB_noPV_psi2S",
+    #                                        Verbose=True ,
+    #                                        UpdateDaughters=True,
+    #                                        Substitutions = { "[ [Lambda_b0]cc -> ^J/psi(1S) Lambda(1520)0 ]CC" : "psi(2S)" },
+    #                                        daughtersToConstrain = [ "psi(2S)" ] 
+    #                                       )
+    #theTuple.Lb.addTool(dtfSigB_noPV_psi2S)
+
+    # --- with PV constraint (nominal)
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_PV"]
+    dtfSigB_PV = TupleToolDecayTreeFitter("SigB_PV",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            UpdateDaughters=True
+                                           )
+    theTuple.Lb.addTool(dtfSigB_PV)
+
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_PV_Lb_mc"]
+    dtfSigB_PV_Lb_mc = TupleToolDecayTreeFitter("SigB_PV_Lb_mc",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            UpdateDaughters=True,
+                                            daughtersToConstrain = [ "Lambda_b0" ] 
+                                           )
+    theTuple.Lb.addTool(dtfSigB_PV_Lb_mc)
+
+
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_PV_Jpsi"]
+    dtfSigB_PV_Jpsi = TupleToolDecayTreeFitter("SigB_PV_Jpsi",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            UpdateDaughters=True,
+                                            daughtersToConstrain = [ "J/psi(1S)" ] 
+                                           )
+    theTuple.Lb.addTool(dtfSigB_PV_Jpsi)
+
+
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_PV_Jpsi_Lb_mc"]
+    dtfSigB_PV_Jpsi_Lb_mc = TupleToolDecayTreeFitter("SigB_PV_Jpsi_Lb_mc",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            UpdateDaughters=True,
+                                            daughtersToConstrain = [ "J/psi(1S)", "Lambda_b0" ] 
+                                           )
+    theTuple.Lb.addTool(dtfSigB_PV_Jpsi_Lb_mc)
+
+
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SigB_PV_psi2S"]
+    dtfSigB_PV_psi2S = TupleToolDecayTreeFitter("SigB_PV_psi2S",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            UpdateDaughters=True,
+                                            Substitutions = { "Lambda_b0 -> ^J/psi(1S) Lambda(1520)0" : "psi(2S)",
+                                                              "Lambda_b~0 -> ^J/psi(1S) Lambda(1520)~0" : "psi(2S)" 
+                                                            },
+                                            daughtersToConstrain = [ "psi(2S)" ] 
+                                           )
+    theTuple.Lb.addTool(dtfSigB_PV_psi2S)
+
+
+    """
+    # the B0 -> K pi l+ l- bkgd 
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SubBd_PV"]
+    dtfSubBd_PV = TupleToolDecayTreeFitter("SubBd_PV",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            Substitutions = { "Lambda_b0 -> J/psi(1S)  (Lambda(1520)0 -> ^p+ X-)"  : "pi+",
+                                                              "Lambda_b~0 -> J/psi(1S)    (Lambda(1520)~0 -> ^p~- X+)" : "pi-" 
+                                                            }
+                                           )
+    theTuple.Lb.addTool(dtfSubBd_PV)
+
+    # the B0s -> K K l+ l- bkgd 
+    theTuple.Lb.ToolList +=  ["TupleToolDecayTreeFitter/SubBs_PV"]
+    dtfSubBs_PV = TupleToolDecayTreeFitter("SubBs_PV",
+                                            Verbose=True ,
+                                            constrainToOriginVertex=True ,
+                                            Substitutions = { "Lambda_b0 -> J/psi(1S)   (Lambda(1520)0 -> ^p+ X-)"  : "K+",
+                                                              "Lambda_b~0 -> J/psi(1S)  (Lambda(1520)~0 -> ^p~- X+)" : "K-" 
+                                                            }
+                                           )
+    theTuple.Lb.addTool(dtfSubBs_PV)
+    """
+
+
+    return theTuple
+
+
+###################################################
+# automatically_configure
+# Use the CondDB to get the right database tags for data
+# from Configurables import CondDB
+# if not isMC: CondDB(LatestGlobalTagByDataType = year )
+###################################################
+
+dv = DaVinci(
+     #InputType     = inputType, #automatically_configure # BD localtest
+     #DataType      =  year  , # automatically_configure # BD localtest
+     EvtMax        =  evtmax  , # -1 for all files
+     TupleFile     = 'DTT.root',
+     PrintFreq     = 100000,
+     Simulation    = isMC,
+     Lumi          = False,
+)
+
+#MessageSvc().OutputLevel = DEBUG
+
+"""
+# Matt Needham(?)'s momentum scale calibration for data
+from Configurables import TrackScaleState
+if not isMC:
+   scaler            = TrackScaleState('scaler')
+   #scaler.Input = "/Event/Leptonic/Rec/Track/Best"
+   dv.UserAlgorithms =  [ scaler ]
+"""
+
+# smear proto particles and their associated tracks for MC
+# https://gitlab.cern.ch/LHCb-RD/ewp-bd2ksteeangular-legacy/-/blob/master/samples_and_selection/DaVinci_scripts/options_Bu2LLK_EE.py#L121
+#from Configurables import TrackSmearState
+# removed 10th march 2022
+#if isMC:
+#   smear = TrackSmearState("TrackSmearState")
+#   dv.UserAlgorithms =  [ smear ]
+
+# scaler/smear must be executed before the analysis algorithms
+isDST = True
+if inputType == "MDST": isDST = False
+
+if doMuon:
+   #print(isRun2, isMC, isDST) # BD localtest
+   Lbtuple_mm = addNtupleToDaVinci("Lb2pkmm", isRun2, isMC, isDST)
+   dv.UserAlgorithms += [Lbtuple_mm]
+
+if doElectron: 
+   #print("doElectron: ",doElectron) # BD localtest
+   Lbtuple_ee = addNtupleToDaVinci("Lb2pkee", isRun2, isMC, isDST)
+   dv.UserAlgorithms += [Lbtuple_ee]
+
+import glob
+#evt+std://LHCb/Collision12/90000000/Beam4000GeV-VeloClosed-MagUp/Real Data/Reco14/Stripping21r0p2/LEPTONIC.MDST
+#import glob
+#flist    = glob.glob("/eos/lhcb/grid/prod//lhcb/LHCb/Collision12/LEPTONIC.MDST/00095380/0000/0*.leptonic.mdst")
+#dv.Input = flist[:1]
+
+# evt+std://LHCb/Collision17/90000000/Beam6500GeV-VeloClosed-MagDown/Real Data/Reco17/Stripping29r2p2/LEPTONIC.MDST 
+#dv.Input = ["/eos/lhcb/grid/prod//lhcb/LHCb/Collision17/LEPTONIC.MDST/00141012/0000/00141012_00000060_1.leptonic.mdst"]
+#flist    = glob.glob("/eos/lhcb/grid/prod//lhcb/LHCb/Collision17/LEPTONIC.MDST/00141012/0000/0*.leptonic.mdst")
+#dv.Input = flist[:1]
+
+#MC_15144001_Lb_JpsipKmmphspDecProdCut_ReStrip
+# evt+std://MC/2018/15144001/Beam6500GeV-2018-MagDown-Nu1.6-25ns-Pythia8/Sim09f/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34NoPrescalingFlagged/ALLSTREAMS.DST
+#dv.Input = ["/eos/lhcb/grid/prod//lhcb/MC/2018/ALLSTREAMS.DST/00086923/0000/00086923_00000406_7.AllStreams.dst"]
+
+# MC_15114401_Lb_pKstmmKpi0phspDecProdCut
+# evt+std://MC/2012/15114401/Beam4000GeV-2012-MagUp-Nu2.5-Pythia8/Sim09k/Trig0x409f0045/Reco14c/Stripping21r0p2NoPrescalingFlagged/ALLSTREAMS.DST
+#dv.Input = ["/eos/lhcb/grid/prod//lhcb/MC/2012/ALLSTREAMS.DST/00149395/0000/00149395_00000110_5.AllStreams.dst"]
+
+#MC_15144059_Lb_JpsipKmmXLL2_ReStrip
+# evt+std://MC/2012/15144059/Beam4000GeV-2012-MagDown-Nu2.5-Pythia8/Sim09h/Trig0x409f0045/Reco14c/Stripping21NoPrescalingFlagged/ALLSTREAMS.DST
+#dv.Input = ["/eos/lhcb/grid/prod//lhcb//MC/2012/ALLSTREAMS.DST/00095494/0000/00095494_00000001_5.AllStreams.dst"]
+
+#evt+std://MC/2012/15154011/Beam4000GeV-2012-MagUp-Nu2.5-Pythia8/Sim09l-ReDecay01/Trig0x409f0045/Reco14c/Stripping21r0p2Filtered/BU2LLK_B2LLXBDT_B2XMUMU.STRIP.MDST
+# Lb2pKee
+#flist    = glob.glob("/eos/lhcb/grid/prod//lhcb/MC/2017/BU2LLK_B2LLXBDT_B2XMUMU.STRIP.MDST/00158235/0000/00158235_00000001_1.bu2llk_b2llxbdt_b2xmumu.strip.mdst")
+# Lb2psi2SpK,ee
+#flist    = glob.glob("/eos/lhcb/grid/prod//lhcb/MC/2012/BU2LLK_B2LLXBDT_B2XMUMU.STRIP.MDST/00157640/0000/00157640_000000*_1.bu2llk_b2llxbdt_b2xmumu.strip.mdst")
+
+#15114012 (Lb_pKmumu=XLL,DecProdCut) Run1 dst
+#flist    = glob.glob("/eos/lhcb/grid/prod//lhcb/MC/2012/BU2LLK_B2LLXBDT_B2XMUMU.STRIP.DST/00166104/0000/00166104_00000*bu2llk_b2llxbdt_b2xmumu.strip.dst")
+#dv.Input = flist[:1]
+
+
+
+
+
diff --git a/Lb2pKll/helpers/__init__.py b/Lb2pKll/helpers/__init__.py
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/Lb2pKll/helpers/helper.py b/Lb2pKll/helpers/helper.py
new file mode 100644
index 0000000000000000000000000000000000000000..95de5f90e549fdad8ba006c0fe33cd530603c72b
--- /dev/null
+++ b/Lb2pKll/helpers/helper.py
@@ -0,0 +1,177 @@
+def fillIsoVars(line, isRun2):
+    lineQ = "'%s"%line
+    print("RELINFO("+lineQ+"TrackIsoInfoL1', 'CONEPT', -999.)")
+    var_dict = {
+    # https://gitlab.cern.ch/LHCb-RD/ewp-rpk/-/blob/master/DaVinci/LeptonU-options-NOPID.py#L207
+                 #track
+ "L1_cpt_0.50": "RELINFO("+lineQ+"TrackIsoInfoL1', 'CONEPT', -999.)",
+ "L2_cpt_0.50": "RELINFO("+lineQ+"TrackIsoInfoL2', 'CONEPT', -999.)",
+ "Proton_cpt_0.50": "RELINFO("+lineQ+"TrackIsoInfoH1', 'CONEPT', -999.)",
+ "Kaon_cpt_0.50": "RELINFO("+lineQ+"TrackIsoInfoH2', 'CONEPT', -999.)",
+
+ "L1_cp_0.50": "RELINFO("+lineQ+"TrackIsoInfoL1', 'CONEP', -999.)",
+ "L2_cp_0.50": "RELINFO("+lineQ+"TrackIsoInfoL2', 'CONEP', -999.)",
+ "Proton_cp_0.50": "RELINFO("+lineQ+"TrackIsoInfoH1', 'CONEP', -999.)",
+ "Kaon_cp_0.50": "RELINFO("+lineQ+"TrackIsoInfoH2', 'CONEP', -999.)",
+
+ "L1_cmult_0.50": "RELINFO("+lineQ+"TrackIsoInfoL1', 'CONEMULT', -1.)",
+ "L2_cmult_0.50": "RELINFO("+lineQ+"TrackIsoInfoL2', 'CONEMULT', -1.)",
+ "Proton_cmult_0.50": "RELINFO("+lineQ+"TrackIsoInfoH1', 'CONEMULT', -1.)",
+ "Kaon_cmult_0.50": "RELINFO("+lineQ+"TrackIsoInfoH2', 'CONEMULT', -1.)",
+
+ "L1_ptasy_0.50": "RELINFO("+lineQ+"TrackIsoInfoL1', 'CONEPTASYM', -999.)",
+ "L2_ptasy_0.50": "RELINFO("+lineQ+"TrackIsoInfoL2', 'CONEPTASYM', -999.)",
+ "Proton_ptasy_0.50": "RELINFO("+lineQ+"TrackIsoInfoH1', 'CONEPTASYM', -999.)",
+ "Kaon_ptasy_0.50": "RELINFO("+lineQ+"TrackIsoInfoH2', 'CONEPTASYM', -999.)",
+
+ "L1_pasy_0.50": "RELINFO("+lineQ+"TrackIsoInfoL1', 'CONEPASYM', -999.)",
+ "L2_pasy_0.50": "RELINFO("+lineQ+"TrackIsoInfoL2', 'CONEPASYM', -999.)",
+ "Proton_pasy_0.50": "RELINFO("+lineQ+"TrackIsoInfoH1', 'CONEPASYM', -999.)",
+ "Kaon_pasy_0.50": "RELINFO("+lineQ+"TrackIsoInfoH2', 'CONEPASYM', -999.)",
+
+ "L1_deltaEta_0.50":  "RELINFO("+lineQ+"TrackIsoInfoL1', 'CONEDELTAETA', -999.)",
+ "L2_deltaEta_0.50":  "RELINFO("+lineQ+"TrackIsoInfoL2', 'CONEDELTAETA', -999.)",
+ "Proton_deltaEta_0.50":  "RELINFO("+lineQ+"TrackIsoInfoH1', 'CONEDELTAETA', -999.)",
+ "Kaon_deltaEta_0.50":  "RELINFO("+lineQ+"TrackIsoInfoH2', 'CONEDELTAETA', -999.)",
+
+ "L1_deltaPhi_0.50":  "RELINFO("+lineQ+"TrackIsoInfoL1', 'CONEDELTAPHI', -999.)",
+ "L2_deltaPhi_0.50":  "RELINFO("+lineQ+"TrackIsoInfoL2', 'CONEDELTAPHI', -999.)",
+ "Proton_deltaPhi_0.50":  "RELINFO("+lineQ+"TrackIsoInfoH1', 'CONEDELTAPHI', -999.)",
+ "Kaon_deltaPhi_0.50":  "RELINFO("+lineQ+"TrackIsoInfoH2', 'CONEDELTAPHI', -999.)",
+
+                 #cone
+ "L1_0.50_cc_mult": "RELINFO("+lineQ+"ConeIsoInfoL1', 'CC_MULT', -1.)",
+ "L2_0.50_cc_mult": "RELINFO("+lineQ+"ConeIsoInfoL2', 'CC_MULT', -1.)",
+ "Proton_0.50_cc_mult": "RELINFO("+lineQ+"ConeIsoInfoH1', 'CC_MULT', -1.)",
+ "Kaon_0.50_cc_mult": "RELINFO("+lineQ+"ConeIsoInfoH2', 'CC_MULT', -1.)",
+
+ "L1_0.50_cc_sPT": "RELINFO("+lineQ+"ConeIsoInfoL1', 'CC_SPT', -1.)",
+ "L2_0.50_cc_sPT": "RELINFO("+lineQ+"ConeIsoInfoL2', 'CC_SPT', -1.)",
+ "Proton_0.50_cc_sPT": "RELINFO("+lineQ+"ConeIsoInfoH1', 'CC_SPT', -1.)",
+ "Kaon_0.50_cc_sPT": "RELINFO("+lineQ+"ConeIsoInfoH2', 'CC_SPT', -1.)",
+
+ "L1_0.50_cc_vPT": "RELINFO("+lineQ+"ConeIsoInfoL1', 'CC_VPT', -1.)",
+ "L2_0.50_cc_vPT": "RELINFO("+lineQ+"ConeIsoInfoL2', 'CC_VPT', -1.)",
+ "Proton_0.50_cc_vPT": "RELINFO("+lineQ+"ConeIsoInfoH1', 'CC_VPT', -1.)",
+ "Kaon_0.50_cc_vPT": "RELINFO("+lineQ+"ConeIsoInfoH2', 'CC_VPT', -1.)",
+
+ "L1_0.50_cc_asy_P": "RELINFO("+lineQ+"ConeIsoInfoL1', 'CC_PASYM', -1.)",
+ "L2_0.50_cc_asy_P": "RELINFO("+lineQ+"ConeIsoInfoL2', 'CC_PASYM', -1.)",
+ "Proton_0.50_cc_asy_P": "RELINFO("+lineQ+"ConeIsoInfoH1', 'CC_PASYM', -1.)",
+ "Kaon_0.50_cc_asy_P": "RELINFO("+lineQ+"ConeIsoInfoH2', 'CC_PASYM', -1.)",
+
+ "L1_0.50_cc_asy_PT": "RELINFO("+lineQ+"ConeIsoInfoL1', 'CC_PTASYM', -1.)",
+ "L2_0.50_cc_asy_PT": "RELINFO("+lineQ+"ConeIsoInfoL2', 'CC_PTASYM', -1.)",
+ "Proton_0.50_cc_asy_PT": "RELINFO("+lineQ+"ConeIsoInfoH1', 'CC_PTASYM', -1.)",
+ "Kaon_0.50_cc_asy_PT": "RELINFO("+lineQ+"ConeIsoInfoH2', 'CC_PTASYM', -1.)",
+
+ "L1_0.50_cc_deltaEta": "RELINFO("+lineQ+"ConeIsoInfoL1', 'CC_DELTAETA', -1.)",
+ "L2_0.50_cc_deltaEta": "RELINFO("+lineQ+"ConeIsoInfoL2', 'CC_DELTAETA', -1.)",
+ "Proton_0.50_cc_deltaEta": "RELINFO("+lineQ+"ConeIsoInfoH1', 'CC_DELTAETA', -1.)",
+ "Kaon_0.50_cc_deltaEta": "RELINFO("+lineQ+"ConeIsoInfoH2', 'CC_DELTAETA', -1.)",
+
+ "L1_0.50_cc_deltaPhi": "RELINFO("+lineQ+"ConeIsoInfoL1', 'CC_DELTAPHI', -1.)",
+ "L2_0.50_cc_deltaPhi": "RELINFO("+lineQ+"ConeIsoInfoL2', 'CC_DELTAPHI', -1.)",
+ "Proton_0.50_cc_deltaPhi": "RELINFO("+lineQ+"ConeIsoInfoH1', 'CC_DELTAPHI', -1.)",
+ "Kaon_0.50_cc_deltaPhi": "RELINFO("+lineQ+"ConeIsoInfoH2', 'CC_DELTAPHI', -1.)",
+
+ "L1_0.50_cc_IT": "RELINFO("+lineQ+"ConeIsoInfoL1', 'CC_IT', -1.)",
+ "L2_0.50_cc_IT": "RELINFO("+lineQ+"ConeIsoInfoL2', 'CC_IT', -1.)",
+ "Proton_0.50_cc_IT": "RELINFO("+lineQ+"ConeIsoInfoH1', 'CC_IT', -1.)",
+ "Kaon_0.50_cc_IT": "RELINFO("+lineQ+"ConeIsoInfoH2', 'CC_IT', -1.)",
+
+ "L1_0.50_cc_maxPt_Q": "RELINFO("+lineQ+"ConeIsoInfoL1', 'CC_MAXPT_Q', -1.)",
+ "L2_0.50_cc_maxPt_Q": "RELINFO("+lineQ+"ConeIsoInfoL2', 'CC_MAXPT_Q', -1.)",
+ "Proton_0.50_cc_maxPt_Q": "RELINFO("+lineQ+"ConeIsoInfoH1', 'CC_MAXPT_Q', -1.)",
+ "Kaon_0.50_cc_maxPt_Q": "RELINFO("+lineQ+"ConeIsoInfoH2', 'CC_MAXPT_Q', -1.)",
+
+ "L1_0.50_cc_maxPt_PT": "RELINFO("+lineQ+"ConeIsoInfoL1', 'CC_MAXPT_PT', -1.)",
+ "L2_0.50_cc_maxPt_PT": "RELINFO("+lineQ+"ConeIsoInfoL2', 'CC_MAXPT_PT', -1.)",
+ "Proton_0.50_cc_maxPt_PT": "RELINFO("+lineQ+"ConeIsoInfoH1', 'CC_MAXPT_PT', -1.)",
+ "Kaon_0.50_cc_maxPt_PT": "RELINFO("+lineQ+"ConeIsoInfoH2', 'CC_MAXPT_PT', -1.)",
+
+ "L1_0.50_cc_maxPt_PE": "RELINFO("+lineQ+"ConeIsoInfoL1', 'CC_MAXPT_PE', -1.)",
+ "L2_0.50_cc_maxPt_PE": "RELINFO("+lineQ+"ConeIsoInfoL2', 'CC_MAXPT_PE', -1.)",
+ "Proton_0.50_cc_maxPt_PE": "RELINFO("+lineQ+"ConeIsoInfoH1', 'CC_MAXPT_PE', -1.)",
+ "Kaon_0.50_cc_maxPt_PE": "RELINFO("+lineQ+"ConeIsoInfoH2', 'CC_MAXPT_PE', -1.)",
+
+                 #cone
+ "L1_0.50_nc_mult": "RELINFO("+lineQ+"ConeIsoInfoL1', 'NC_MULT', -1.)",
+ "L2_0.50_nc_mult": "RELINFO("+lineQ+"ConeIsoInfoL2', 'NC_MULT', -1.)",
+ "Proton_0.50_nc_mult": "RELINFO("+lineQ+"ConeIsoInfoH1', 'NC_MULT', -1.)",
+ "Kaon_0.50_nc_mult": "RELINFO("+lineQ+"ConeIsoInfoH2', 'NC_MULT', -1.)",
+
+ "L1_0.50_nc_sPT": "RELINFO("+lineQ+"ConeIsoInfoL1', 'NC_SPT', -1.)",
+ "L2_0.50_nc_sPT": "RELINFO("+lineQ+"ConeIsoInfoL2', 'NC_SPT', -1.)",
+ "Proton_0.50_nc_sPT": "RELINFO("+lineQ+"ConeIsoInfoH1', 'NC_SPT', -1.)",
+ "Kaon_0.50_nc_sPT": "RELINFO("+lineQ+"ConeIsoInfoH2', 'NC_SPT', -1.)",
+
+
+ "L1_0.50_nc_vPT": "RELINFO("+lineQ+"ConeIsoInfoL1', 'NC_VPT', -1.)",
+ "L2_0.50_nc_vPT": "RELINFO("+lineQ+"ConeIsoInfoL2', 'NC_VPT', -1.)",
+ "Proton_0.50_nc_vPT": "RELINFO("+lineQ+"ConeIsoInfoH1', 'NC_VPT', -1.)",
+ "Kaon_0.50_nc_vPT": "RELINFO("+lineQ+"ConeIsoInfoH2', 'NC_VPT', -1.)",
+
+ "L1_0.50_nc_asy_P": "RELINFO("+lineQ+"ConeIsoInfoL1', 'NC_PASYM', -1.)",
+ "L2_0.50_nc_asy_P": "RELINFO("+lineQ+"ConeIsoInfoL2', 'NC_PASYM', -1.)",
+ "Proton_0.50_nc_asy_P": "RELINFO("+lineQ+"ConeIsoInfoH1', 'NC_PASYM', -1.)",
+ "Kaon_0.50_nc_asy_P": "RELINFO("+lineQ+"ConeIsoInfoH2', 'NC_PASYM', -1.)",
+
+ "L1_0.50_nc_asy_PT": "RELINFO("+lineQ+"ConeIsoInfoL1', 'NC_PTASYM', -1.)",
+ "L2_0.50_nc_asy_PT": "RELINFO("+lineQ+"ConeIsoInfoL2', 'NC_PTASYM', -1.)",
+ "Proton_0.50_nc_asy_PT": "RELINFO("+lineQ+"ConeIsoInfoH1', 'NC_PTASYM', -1.)",
+ "Kaon_0.50_nc_asy_PT": "RELINFO("+lineQ+"ConeIsoInfoH2', 'NC_PTASYM', -1.)",
+
+ "L1_0.50_nc_deltaEta": "RELINFO("+lineQ+"ConeIsoInfoL1', 'NC_DELTAETA', -1.)",
+ "L2_0.50_nc_deltaEta": "RELINFO("+lineQ+"ConeIsoInfoL2', 'NC_DELTAETA', -1.)",
+ "Proton_0.50_nc_deltaEta": "RELINFO("+lineQ+"ConeIsoInfoH1', 'NC_DELTAETA', -1.)",
+ "Kaon_0.50_nc_deltaEta": "RELINFO("+lineQ+"ConeIsoInfoH2', 'NC_DELTAETA', -1.)",
+
+ "L1_0.50_nc_deltaPhi": "RELINFO("+lineQ+"ConeIsoInfoL1', 'NC_DELTAPHI', -1.)",
+ "L2_0.50_nc_deltaPhi": "RELINFO("+lineQ+"ConeIsoInfoL2', 'NC_DELTAPHI', -1.)",
+ "Proton_0.50_nc_deltaPhi": "RELINFO("+lineQ+"ConeIsoInfoH1', 'NC_DELTAPHI', -1.)",
+ "Kaon_0.50_nc_deltaPhi": "RELINFO("+lineQ+"ConeIsoInfoH2', 'NC_DELTAPHI', -1.)",
+
+ "L1_0.50_nc_IT": "RELINFO("+lineQ+"ConeIsoInfoL1', 'NC_IT', -1.)",
+ "L2_0.50_nc_IT": "RELINFO("+lineQ+"ConeIsoInfoL2', 'NC_IT', -1.)",
+ "Proton_0.50_nc_IT": "RELINFO("+lineQ+"ConeIsoInfoH1', 'NC_IT', -1.)",
+ "Kaon_0.50_nc_IT": "RELINFO("+lineQ+"ConeIsoInfoH2', 'NC_IT', -1.)",
+
+ "L1_0.50_nc_maxPt_PT": "RELINFO("+lineQ+"ConeIsoInfoL1', 'NC_MAXPT_PT', -1.)",
+ "L2_0.50_nc_maxPt_PT": "RELINFO("+lineQ+"ConeIsoInfoL2', 'NC_MAXPT_PT', -1.)",
+ "Proton_0.50_nc_maxPt_PT": "RELINFO("+lineQ+"ConeIsoInfoH1', 'NC_MAXPT_PT', -1.)",
+ "Kaon_0.50_nc_maxPt_PT": "RELINFO("+lineQ+"ConeIsoInfoH2', 'NC_MAXPT_PT', -1.)",
+
+###vertex variable
+ "NumVtxWithinChi2WindowOneTrack":  "RELINFO("+lineQ+"VertexIsoInfo', 'VTXISONUMVTX', -999.)",
+ "SmallestDeltaChi2OneTrack":  "RELINFO("+lineQ+"VertexIsoInfo', 'VTXISODCHI2ONETRACK', -999.)",
+ "SmallestDeltaChi2MassOneTrack":  "RELINFO("+lineQ+"VertexIsoInfo', 'VTXISODCHI2MASSONETRACK', -999.)",
+ "SmallestDeltaChi2TwoTracks":  "RELINFO("+lineQ+"VertexIsoInfo', 'VTXISODCHI2TWOTRACK', -999.)",
+ "SmallestDeltaChi2MassTwoTracks":  "RELINFO("+lineQ+"VertexIsoInfo', 'VTXISODCHI2MASSTWOTRACK', -999.)",
+ "NumVtxWithinChi2WindowOneTrack":  "RELINFO("+lineQ+"VertexIsoInfo', 'VTXISONUMVTX', -999.)",
+               }
+
+    var_dict_R1 = {
+#  new Bs2MM vars for leptons
+ "TRKISOBDTLONG_12":  "RELINFO("+lineQ+"TrackIsolationBDT2', 'TRKISOBDTLONG', -999.)",
+ "TRKISOBDTVELO_12":  "RELINFO("+lineQ+"TrackIsolationBDT2', 'TRKISOBDTVELO', -999.)",
+    }
+
+    var_dict_R2 = {
+#  new Bs2MM vars for all combinations
+ "TRKISOBDTLONG_12":  "RELINFO("+lineQ+"TrackIsolationBDT2_12', 'TRKISOBDTLONG', -999.)",
+ "TRKISOBDTVELO_12":  "RELINFO("+lineQ+"TrackIsolationBDT2_12', 'TRKISOBDTVELO', -999.)",
+ "TRKISOBDTLONG_14":  "RELINFO("+lineQ+"TrackIsolationBDT2_14', 'TRKISOBDTLONG', -999.)",
+ "TRKISOBDTVELO_14":  "RELINFO("+lineQ+"TrackIsolationBDT2_14', 'TRKISOBDTVELO', -999.)",
+ "TRKISOBDTLONG_25":  "RELINFO("+lineQ+"TrackIsolationBDT2_25', 'TRKISOBDTLONG', -999.)",
+ "TRKISOBDTVELO_25":  "RELINFO("+lineQ+"TrackIsolationBDT2_25', 'TRKISOBDTVELO', -999.)",
+ "TRKISOBDTLONG_45":  "RELINFO("+lineQ+"TrackIsolationBDT2_45', 'TRKISOBDTLONG', -999.)",
+ "TRKISOBDTVELO_45":  "RELINFO("+lineQ+"TrackIsolationBDT2_45', 'TRKISOBDTVELO', -999.)",
+               }
+
+    if not isRun2: var_dict.update(var_dict_R1)
+    if isRun2: var_dict.update(var_dict_R2)
+    
+    return var_dict
+
+
diff --git a/Lb2pKll/info.yaml b/Lb2pKll/info.yaml
new file mode 100644
index 0000000000000000000000000000000000000000..7f80c8ccfea892eae3f32c7314368bfab7fb053c
--- /dev/null
+++ b/Lb2pKll/info.yaml
@@ -0,0 +1,184 @@
+defaults:
+    application: DaVinci/v46r5
+    output: RD_LBTOPKLL_DVNTUPLE.ROOT
+    wg: RD
+    input:
+        n_test_lfns: 6
+    automatically_configure: yes
+    inform:
+        - biplab.dey@cern.ch
+
+
+{%- set MC13154041BsJKKeeRStrp = [
+    (2016, 1, 0, "Sim09k", "S28r2", "v44r11p4", "MU", "/MC/2016/Beam6500GeV-2016-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/13154041/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09k", "S28r2", "v44r11p4", "MD", "/MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/13154041/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09k", "S29r2p1", "v42r11p2", "MU", "/MC/2017/Beam6500GeV-2017-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2p1NoPrescalingFlagged/13154041/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09k", "S29r2p1", "v42r11p2", "MD", "/MC/2017/Beam6500GeV-2017-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2p1NoPrescalingFlagged/13154041/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09k", "S34r0p1", "v44r11p1", "MU", "/MC/2018/Beam6500GeV-2018-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34r0p1NoPrescalingFlagged/13154041/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09k", "S34r0p1", "v44r11p1", "MD", "/MC/2018/Beam6500GeV-2018-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34r0p1NoPrescalingFlagged/13154041/ALLSTREAMS.DST"),
+]%}
+
+
+{%- set MC13144041BsJKKmmRStrp = [
+    (2012, 1, 0, "Sim09h", "S21r0p2", "v39r1p6", "MU", "/MC/2012/Beam4000GeV-2012-MagUp-Nu2.5-Pythia8/Sim09h/Trig0x409f0045/Reco14c/Stripping21r0p2NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+    (2012, 1, 0, "Sim09h", "S21r0p2", "v39r1p6", "MD", "/MC/2012/Beam4000GeV-2012-MagDown-Nu2.5-Pythia8/Sim09h/Trig0x409f0045/Reco14c/Stripping21r0p2NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+    (2011, 1, 0, "Sim09k", "S21r1p2", "v39r1p6", "MU", "/MC/2011/Beam3500GeV-2011-MagUp-Nu2-Pythia8/Sim09k/Trig0x40760037/Reco14c/Stripping21r1p2NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+    (2011, 1, 0, "Sim09k", "S21r1p2", "v39r1p6", "MD", "/MC/2011/Beam3500GeV-2011-MagDown-Nu2-Pythia8/Sim09k/Trig0x40760037/Reco14c/Stripping21r1p2NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+    (2012, 1, 0, "Sim09k", "S21r0p2", "v39r1p6", "MU", "/MC/2012/Beam4000GeV-2012-MagUp-Nu2.5-Pythia8/Sim09k/Trig0x409f0045/Reco14c/Stripping21r0p2NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+    (2012, 1, 0, "Sim09k", "S21r0p2", "v39r1p6", "MD", "/MC/2012/Beam4000GeV-2012-MagDown-Nu2.5-Pythia8/Sim09k/Trig0x409f0045/Reco14c/Stripping21r0p2NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09k", "S28r2", "v44r11p4", "MU", "/MC/2016/Beam6500GeV-2016-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09k", "S28r2", "v44r11p4", "MD", "/MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09k", "S29r2p1", "v42r11p2", "MU", "/MC/2017/Beam6500GeV-2017-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2p1NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09k", "S29r2p1", "v42r11p2", "MD", "/MC/2017/Beam6500GeV-2017-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2p1NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09k", "S34r0p1", "v44r11p1", "MU", "/MC/2018/Beam6500GeV-2018-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34r0p1NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09k", "S34r0p1", "v44r11p1", "MD", "/MC/2018/Beam6500GeV-2018-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34r0p1NoPrescalingFlagged/13144041/ALLSTREAMS.DST"),
+]%}
+
+
+{%- set MC11144042BdJKpimmXLLRStrp = [
+    (2017, 1, 0, "Sim09h", "S29r2", "v42r11p2", "MU", "/MC/2017/Beam6500GeV-2017-MagUp-Nu1.6-25ns-Pythia8/Sim09h/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2NoPrescalingFlagged/11144042/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09h", "S29r2", "v42r11p2", "MD", "/MC/2017/Beam6500GeV-2017-MagDown-Nu1.6-25ns-Pythia8/Sim09h/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2NoPrescalingFlagged/11144042/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09h", "S34", "v44r11p1", "MU", "/MC/2018/Beam6500GeV-2018-MagUp-Nu1.6-25ns-Pythia8/Sim09h/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34NoPrescalingFlagged/11144042/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09h", "S28r2", "v44r11p4", "MU", "/MC/2016/Beam6500GeV-2016-MagUp-Nu1.6-25ns-Pythia8/Sim09h/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/11144042/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09h", "S28r2", "v44r11p4", "MD", "/MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09h/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/11144042/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09h-ReDecay01", "S34", "v44r11p1", "MD", "/MC/2018/Beam6500GeV-2018-MagDown-Nu1.6-25ns-Pythia8/Sim09h-ReDecay01/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34NoPrescalingFlagged/11144042/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09h-ReDecay01", "S34", "v44r11p1", "MU", "/MC/2018/Beam6500GeV-2018-MagUp-Nu1.6-25ns-Pythia8/Sim09h-ReDecay01/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34NoPrescalingFlagged/11144042/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09h-ReDecay01", "S28r2", "v44r11p4", "MD", "/MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09h-ReDecay01/Trig0x6138160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/11144042/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09h-ReDecay01", "S28r2", "v44r11p4", "MU", "/MC/2016/Beam6500GeV-2016-MagUp-Nu1.6-25ns-Pythia8/Sim09h-ReDecay01/Trig0x6138160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/11144042/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09h-ReDecay01", "S29r2", "v42r11p2", "MD", "/MC/2017/Beam6500GeV-2017-MagDown-Nu1.6-25ns-Pythia8/Sim09h-ReDecay01/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2NoPrescalingFlagged/11144042/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09h-ReDecay01", "S29r2", "v42r11p2", "MU", "/MC/2017/Beam6500GeV-2017-MagUp-Nu1.6-25ns-Pythia8/Sim09h-ReDecay01/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2NoPrescalingFlagged/11144042/ALLSTREAMS.DST"),
+]%}
+
+
+{%- set MC11144050BdJKpimmRStrp = [
+    (2017, 1, 0, "Sim09h", "S29r2", "v42r11p2", "MU", "/MC/2017/Beam6500GeV-2017-MagUp-Nu1.6-25ns-Pythia8/Sim09h/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09h", "S29r2", "v42r11p2", "MD", "/MC/2017/Beam6500GeV-2017-MagDown-Nu1.6-25ns-Pythia8/Sim09h/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09h", "S28r1", "v44r11p4", "MU", "/MC/2016/Beam6500GeV-2016-MagUp-Nu1.6-25ns-Pythia8/Sim09h/Trig0x6139160F/Reco16/Turbo03a/Stripping28r1NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09h", "S28r1", "v44r11p4", "MD", "/MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09h/Trig0x6139160F/Reco16/Turbo03a/Stripping28r1NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2011, 1, 0, "Sim09k", "S21r1", "v39r1p6", "MU", "/MC/2011/Beam3500GeV-2011-MagUp-Nu2-Pythia8/Sim09k/Trig0x40760037/Reco14c/Stripping21r1NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2011, 1, 0, "Sim09k", "S21r1", "v39r1p6", "MD", "/MC/2011/Beam3500GeV-2011-MagDown-Nu2-Pythia8/Sim09k/Trig0x40760037/Reco14c/Stripping21r1NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2012, 1, 0, "Sim09k", "S21", "v39r1p6", "MU", "/MC/2012/Beam4000GeV-2012-MagUp-Nu2.5-Pythia8/Sim09k/Trig0x409f0045/Reco14c/Stripping21NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2012, 1, 0, "Sim09k", "S21", "v39r1p6", "MD", "/MC/2012/Beam4000GeV-2012-MagDown-Nu2.5-Pythia8/Sim09k/Trig0x409f0045/Reco14c/Stripping21NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09l", "S34", "v44r11p1", "MD", "/MC/2018/Beam6500GeV-2018-MagDown-Nu1.6-25ns-Pythia8/Sim09l/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09l", "S34", "v44r11p1", "MU", "/MC/2018/Beam6500GeV-2018-MagUp-Nu1.6-25ns-Pythia8/Sim09l/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09l", "S29r2", "v42r11p2", "MD", "/MC/2017/Beam6500GeV-2017-MagDown-Nu1.6-25ns-Pythia8/Sim09l/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09l", "S29r2", "v42r11p2", "MU", "/MC/2017/Beam6500GeV-2017-MagUp-Nu1.6-25ns-Pythia8/Sim09l/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09l", "S28r2", "v44r11p4", "MD", "/MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09l/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09l", "S28r2", "v44r11p4", "MU", "/MC/2016/Beam6500GeV-2016-MagUp-Nu1.6-25ns-Pythia8/Sim09l/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/11144050/ALLSTREAMS.DST"),
+]%}
+
+
+{%- set MC11154050BdJKpieeRStrp = [
+    (2018, 1, 0, "Sim09l", "S34r0p1", "v44r11p1", "MU", "/MC/2018/Beam6500GeV-2018-MagUp-Nu1.6-25ns-Pythia8/Sim09l/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34r0p1NoPrescalingFlagged/11154050/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09l", "S34r0p1", "v44r11p1", "MD", "/MC/2018/Beam6500GeV-2018-MagDown-Nu1.6-25ns-Pythia8/Sim09l/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34r0p1NoPrescalingFlagged/11154050/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09l", "S29r2p1", "v42r11p2", "MU", "/MC/2017/Beam6500GeV-2017-MagUp-Nu1.6-25ns-Pythia8/Sim09l/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2p1NoPrescalingFlagged/11154050/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09l", "S29r2p1", "v42r11p2", "MD", "/MC/2017/Beam6500GeV-2017-MagDown-Nu1.6-25ns-Pythia8/Sim09l/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2p1NoPrescalingFlagged/11154050/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09l", "S28r2", "v44r11p4", "MU", "/MC/2016/Beam6500GeV-2016-MagUp-Nu1.6-25ns-Pythia8/Sim09l/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/11154050/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09l", "S28r2", "v44r11p4", "MD", "/MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09l/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/11154050/ALLSTREAMS.DST"),
+]%}
+
+
+{%- set MC15144401Lb2JpKpi0KstmmRStrp = [
+    (2011, 1, 0, "Sim09k", "S21r1p2", "v39r1p6", "MU", "/MC/2011/Beam3500GeV-2011-MagUp-Nu2-Pythia8/Sim09k/Trig0x40760037/Reco14c/Stripping21r1p2NoPrescalingFlagged/15144401/ALLSTREAMS.DST"),
+    (2011, 1, 0, "Sim09k", "S21r1p2", "v39r1p6", "MD", "/MC/2011/Beam3500GeV-2011-MagDown-Nu2-Pythia8/Sim09k/Trig0x40760037/Reco14c/Stripping21r1p2NoPrescalingFlagged/15144401/ALLSTREAMS.DST"),
+    (2012, 1, 0, "Sim09k", "S21r0p2", "v39r1p6", "MU", "/MC/2012/Beam4000GeV-2012-MagUp-Nu2.5-Pythia8/Sim09k/Trig0x409f0045/Reco14c/Stripping21r0p2NoPrescalingFlagged/15144401/ALLSTREAMS.DST"),
+    (2012, 1, 0, "Sim09k", "S21r0p2", "v39r1p6", "MD", "/MC/2012/Beam4000GeV-2012-MagDown-Nu2.5-Pythia8/Sim09k/Trig0x409f0045/Reco14c/Stripping21r0p2NoPrescalingFlagged/15144401/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09k", "S28r2", "v44r11p4", "MU", "/MC/2016/Beam6500GeV-2016-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/15144401/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09k", "S28r2", "v44r11p4", "MD", "/MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/15144401/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09k", "S29r2p1", "v42r11p2", "MU", "/MC/2017/Beam6500GeV-2017-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2p1NoPrescalingFlagged/15144401/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09k", "S29r2p1", "v42r11p2", "MD", "/MC/2017/Beam6500GeV-2017-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2p1NoPrescalingFlagged/15144401/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09k", "S34r0p1", "v44r11p1", "MU", "/MC/2018/Beam6500GeV-2018-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34r0p1NoPrescalingFlagged/15144401/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09k", "S34r0p1", "v44r11p1", "MD", "/MC/2018/Beam6500GeV-2018-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34r0p1NoPrescalingFlagged/15144401/ALLSTREAMS.DST"),
+]%}
+
+
+{%- set MC15154401Lb2JpKpi0KsteeRStrp = [
+    (2016, 1, 0, "Sim09k", "S28r2", "v44r11p4", "MU", "/MC/2016/Beam6500GeV-2016-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/15154401/ALLSTREAMS.DST"),
+    (2016, 1, 0, "Sim09k", "S28r2", "v44r11p4", "MD", "/MC/2016/Beam6500GeV-2016-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x6139160F/Reco16/Turbo03a/Stripping28r2NoPrescalingFlagged/15154401/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09k", "S29r2p1", "v42r11p2", "MU", "/MC/2017/Beam6500GeV-2017-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2p1NoPrescalingFlagged/15154401/ALLSTREAMS.DST"),
+    (2017, 1, 0, "Sim09k", "S29r2p1", "v42r11p2", "MD", "/MC/2017/Beam6500GeV-2017-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x62661709/Reco17/Turbo04a-WithTurcal/Stripping29r2p1NoPrescalingFlagged/15154401/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09k", "S34r0p1", "v44r11p1", "MU", "/MC/2018/Beam6500GeV-2018-MagUp-Nu1.6-25ns-Pythia8/Sim09k/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34r0p1NoPrescalingFlagged/15154401/ALLSTREAMS.DST"),
+    (2018, 1, 0, "Sim09k", "S34r0p1", "v44r11p1", "MD", "/MC/2018/Beam6500GeV-2018-MagDown-Nu1.6-25ns-Pythia8/Sim09k/Trig0x617d18a4/Reco18/Turbo05-WithTurcal/Stripping34r0p1NoPrescalingFlagged/15154401/ALLSTREAMS.DST"),
+]%}
+
+
+
+{%- for year, restrip, filt, simver, strpver, dvver, polarity, bkk_path in MC13154041BsJKKeeRStrp %}
+MC13154041BsJKKeeRStrp{{year}}{{polarity}}{{simver}}{{strpver}}DV{{dvver}}:
+    application: DaVinci/{{ dvver }}
+    input:
+        bk_query: {{bkk_path}}
+    options:
+        - helpers/helper.py
+        - steering_files/MC_{{year}}_{{polarity}}_ReStrip{{restrip}}_Filt{{filt}}_ee.py
+        - dv_pKll.py
+{%- endfor %}
+
+{%- for year, restrip, filt, simver, strpver, dvver, polarity, bkk_path in MC13144041BsJKKmmRStrp %}
+MC13144041BsJKKmmRStrp{{year}}{{polarity}}{{simver}}{{strpver}}DV{{dvver}}:
+    application: DaVinci/{{ dvver }}
+    input:
+        bk_query: {{bkk_path}}
+    options:
+        - helpers/helper.py
+        - steering_files/MC_{{year}}_{{polarity}}_ReStrip{{restrip}}_Filt{{filt}}_mumu.py
+        - dv_pKll.py
+{%- endfor %}
+
+{%- for year, restrip, filt, simver, strpver, dvver, polarity, bkk_path in MC11144042BdJKpimmXLLRStrp %}
+MC11144042BdJKpimmXLLRStrp{{year}}{{polarity}}{{simver}}{{strpver}}DV{{dvver}}:
+    application: DaVinci/{{ dvver }}
+    input:
+        bk_query: {{bkk_path}}
+    options:
+        - helpers/helper.py
+        - steering_files/MC_{{year}}_{{polarity}}_ReStrip{{restrip}}_Filt{{filt}}_mumu.py
+        - dv_pKll.py
+{%- endfor %}
+
+{%- for year, restrip, filt, simver, strpver, dvver, polarity, bkk_path in MC11144050BdJKpimmRStrp %}
+MC11144050BdJKpimmRStrp{{year}}{{polarity}}{{simver}}{{strpver}}DV{{dvver}}:
+    application: DaVinci/{{ dvver }}
+    input:
+        bk_query: {{bkk_path}}
+    options:
+        - helpers/helper.py
+        - steering_files/MC_{{year}}_{{polarity}}_ReStrip{{restrip}}_Filt{{filt}}_mumu.py
+        - dv_pKll.py
+{%- endfor %}
+
+{%- for year, restrip, filt, simver, strpver, dvver, polarity, bkk_path in MC11154050BdJKpieeRStrp %}
+MC11154050BdJKpieeRStrp{{year}}{{polarity}}{{simver}}{{strpver}}DV{{dvver}}:
+    application: DaVinci/{{ dvver }}
+    input:
+        bk_query: {{bkk_path}}
+    options:
+        - helpers/helper.py
+        - steering_files/MC_{{year}}_{{polarity}}_ReStrip{{restrip}}_Filt{{filt}}_ee.py
+        - dv_pKll.py
+{%- endfor %}
+
+
+{%- for year, restrip, filt, simver, strpver, dvver, polarity, bkk_path in MC15144401Lb2JpKpi0KstmmRStrp %}
+MC15144401Lb2JpKpi0KstmmRStrp{{year}}{{polarity}}{{simver}}{{strpver}}DV{{dvver}}:
+    application: DaVinci/{{ dvver }}
+    input:
+        bk_query: {{bkk_path}}
+    options:
+        - helpers/helper.py
+        - steering_files/MC_{{year}}_{{polarity}}_ReStrip{{restrip}}_Filt{{filt}}_mumu.py
+        - dv_pKll.py
+{%- endfor %}
+
+{%- for year, restrip, filt, simver, strpver, dvver, polarity, bkk_path in MC15154401Lb2JpKpi0KsteeRStrp %}
+MC15154401Lb2JpKpi0KsteeRStrp{{year}}{{polarity}}{{simver}}{{strpver}}DV{{dvver}}:
+    application: DaVinci/{{ dvver }}
+    input:
+        bk_query: {{bkk_path}}
+    options:
+        - helpers/helper.py
+        - steering_files/MC_{{year}}_{{polarity}}_ReStrip{{restrip}}_Filt{{filt}}_ee.py
+        - dv_pKll.py
+{%- endfor %}
+
+
diff --git a/Lb2pKll/steering_files/Data_2011.py b/Lb2pKll/steering_files/Data_2011.py
new file mode 100644
index 0000000000000000000000000000000000000000..69cfd14366f651ec4ce7e3ee24c8417e349e1422
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2011.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/Data_2011_ee.py b/Lb2pKll/steering_files/Data_2011_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..5db1f654d80c92153b24318d7e98006cbd48e1f2
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2011_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/Data_2011_mumu.py b/Lb2pKll/steering_files/Data_2011_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..5efe76eab0bfb30d14eab0f19cab36743e0b46d9
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2011_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/Data_2012.py b/Lb2pKll/steering_files/Data_2012.py
new file mode 100644
index 0000000000000000000000000000000000000000..bddb8e8a3bcba2246224fb702121caab0c2594d2
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2012.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/Data_2012_ee.py b/Lb2pKll/steering_files/Data_2012_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..0be68a1abd5c3081d23605ffadc25b4ea1f0e8f9
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2012_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/Data_2012_mumu.py b/Lb2pKll/steering_files/Data_2012_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..dc1cb989886e49088e1c2e09ccb6aa5c990f7bd7
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2012_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/Data_2015.py b/Lb2pKll/steering_files/Data_2015.py
new file mode 100644
index 0000000000000000000000000000000000000000..ac46db2f31c13aab2a9f1f22c72b51e1b5098739
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2015.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2015'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/Data_2016.py b/Lb2pKll/steering_files/Data_2016.py
new file mode 100644
index 0000000000000000000000000000000000000000..68ab0d8db0a06ba61ee8cc4eeb16337f9cad42ce
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2016.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/Data_2016_ee.py b/Lb2pKll/steering_files/Data_2016_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..bddbeabc22e2dadc06ff4f04a4dcec4139027fb0
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2016_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/Data_2016_mumu.py b/Lb2pKll/steering_files/Data_2016_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..538634fbbbdd1d1e725a9926e2afe77edaee61d5
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2016_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/Data_2017.py b/Lb2pKll/steering_files/Data_2017.py
new file mode 100644
index 0000000000000000000000000000000000000000..9404a3b0a02fa86b51ec4e3eb575b15f5e34a68f
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2017.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/Data_2017_ee.py b/Lb2pKll/steering_files/Data_2017_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..8206f079771746182d978fd343f0db766b7f4807
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2017_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/Data_2017_mumu.py b/Lb2pKll/steering_files/Data_2017_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..a2269c9f0bae16412a7ae86bec36feddedad8954
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2017_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/Data_2018.py b/Lb2pKll/steering_files/Data_2018.py
new file mode 100644
index 0000000000000000000000000000000000000000..e2456f5a74b321fc8e64fc52dae00834b712aa48
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2018.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/Data_2018_ee.py b/Lb2pKll/steering_files/Data_2018_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..caf238f172607fd9d1f12e3af6223a162e68eba3
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2018_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/Data_2018_mumu.py b/Lb2pKll/steering_files/Data_2018_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..f3e8978b1a366fc71a5fc01bd6e5f6684bd0c83b
--- /dev/null
+++ b/Lb2pKll/steering_files/Data_2018_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "False"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt0_ee.py b/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..5097d6fb9ea2fd6e000747fe31995340ef2b0d04
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..eac909f3e533b6c6d308b4d0f4e3c8c6b41504dd
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt1_ee.py b/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..e4f593f8fb533cae4f8861c6edc9b179a6c1c701
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..e5bb6e43aa4bd4b2fcdb22bc63cbf056f5df0726
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MD_ReStrip0_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt0_ee.py b/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..e63a332e0f0c67c9038cb9dd52b3c1d0a7ce128a
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..7fea205f835330d49ba450aa6f98b8bd43528eb5
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt1_ee.py b/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..18d447c6315d95cb507e2fd581264e955ba42e12
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..c796240ad885a4f291b1a3254980108b826d5fdd
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MD_ReStrip1_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt0_ee.py b/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..9f500c961182af5130604b248316eb52d8d68401
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..d15bfdbeabb0b034cf2de1ea2ce34294e9550d35
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt1_ee.py b/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..f81a9131e26381681244db0bc9bfc87d63ceb546
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..4e15a5147c3972cf9574de0f2ba0f123585d8d30
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MU_ReStrip0_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt0_ee.py b/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..73fa7e07fe91042b9463241cdf799dd800937c5f
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..a46e1703c67aafb0d4c9cd9d069b62ab89dd99e5
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt1_ee.py b/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..406de8e1eaa6c7998d3e9288c2d2802b00456a1b
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..9ff5a4d89284fb05f29ed46f0c9417d3e7833da0
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2011_MU_ReStrip1_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2011'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt0_ee.py b/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..cf4e6985d1696cbf626b5b206bf0bf564c17bcaf
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..b6ef0fc2649f467288dffb862f51f44583d9d302
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt1_ee.py b/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..4827e0b86d1915c6cb2d37c6dc7d12183669f513
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..6d3f875074e1b681302a7a5d4681867baa682d28
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MD_ReStrip0_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt0_ee.py b/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..983e7493a7dac7d47c911665c0b6ff3ff4b49443
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..ad11d80b7de583f07374bdf79beb6f4323688e2f
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt1_ee.py b/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..1d79afd34dac080842db6f4f3aadc952ffe3c6f9
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..be75be54f28273c98b276026917a6f7bbbbbce26
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MD_ReStrip1_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt0_ee.py b/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..4412f156c78808dcf6e43c16e64ce5f08bad0924
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..edcde154e939c462203f1a8f13ed470129ae585c
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt1_ee.py b/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..6975f59ed431afe33e6e25bb3000a4a5736598bf
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..af88a6fe0956e30dd3a51c72a8ba2a0a4c0bd69a
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MU_ReStrip0_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt0_ee.py b/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..a977e4207713dbe9f48e95c712262c36f058bc47
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..1889d575584ac56493a40db5dd6a3e22501eb2b8
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt1_ee.py b/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..986b2cfbdce3615ddf2b4af3c834167b26ac09e7
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..b1f065a9f76d21ab7df3e84898001c4a81a55f0d
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2012_MU_ReStrip1_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "False"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2012'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt0_ee.py b/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..ec702ca6e450e3d411e7fa331f4308f685425f1d
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..5de1759852e4bd8a496fc23639fef6375226d150
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt1_ee.py b/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..bc8371d573b4c451c4404b4cc1b0865c2684c80d
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..4cd66fa7e1b35680a8eb7bfe00ecbb39bbfb22be
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MD_ReStrip0_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt0_ee.py b/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..2c052043f2c0c7052957097e7b1e396ff0e63cfb
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..5b9921b6244bedec5163fa4d33280d904a9acbb4
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt1_ee.py b/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..919b4de91e8beffe9ed6dd82772ee92ea97fa172
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..d52582fac78c21cf3b955c0186a3063451a46246
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MD_ReStrip1_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt0_ee.py b/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..c14a5c12b54f4fc185df753ba2a2297cf25e9026
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..f952bb5d3b22492dfe437af3164d1f65d9752488
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt1_ee.py b/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..4fd234c9bbb59040e97f6d3b8c2247b8b531d07c
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..4623632fbb829f677f5d21651080499aea02c5e9
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MU_ReStrip0_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt0_ee.py b/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..c331123042f8191516548acfa1eace12231564af
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..a29300b9dead0623ff588e0027881df1a1bfa135
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt1_ee.py b/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..b086927e2b9d183b22d6e20f7aaa60345f65cb64
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..c33de9a124883ce0fb9f312429557941acac86e9
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2016_MU_ReStrip1_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2016'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt0_ee.py b/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..0f2b0548e79d9fa7352f97727fa25536aa0e85af
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..a1f5584d7afaa6b673f2591327027c2c39ba40c8
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt1_ee.py b/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..5aa615f974d03702b75295464b31e771cdf43bce
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..d7ec9ad947784ce037d7b1669adbf1b17566a9e9
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MD_ReStrip0_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt0_ee.py b/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..46d569d4a69f19c65443f9d2616f00c75d3bcf50
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..aa06232da126b54120d7b9579639b34b9109b435
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt1_ee.py b/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..7ec28f43e85d15bade2ee7d8e11210ef8b84b359
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..b9adb7a587889dace9f77fca6425c70d875fde65
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MD_ReStrip1_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt0_ee.py b/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..eb4c722e1dc74c0772dbb77c4323265acadbd116
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..513e298946000e10b124f823a405fe8004e62db5
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt1_ee.py b/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..1d81757d3984b4ae64d55c161a2384c71fba463b
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..cb6adbec77aba9419890b60c40be5b3de35a05d7
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MU_ReStrip0_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt0_ee.py b/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..feccea999119cf155ef22436ae2eb5d765a99c56
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..77709c4fe91e0a42e9f8f44ba3810b338c1f6df9
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt1_ee.py b/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..d9dcb78bdd3933148495aac9828eee431440865a
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..f24fa4fe6cb7ff2041ec81044234c86847290505
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2017_MU_ReStrip1_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2017'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt0_ee.py b/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..b6a7b490a777b567b4e6420ae41f354e70bae701
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..7d0b9e6d9029d503bbc9fcda3a2aa4f19f5e6a99
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt1_ee.py b/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..cd4e87c06fd03f0d6ee91fcce203cf56ca225e30
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..180fb33ee16671f892b282e098f72df916d42af7
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MD_ReStrip0_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt0_ee.py b/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..e962b82523d8975e7273e79352140c75b6732d46
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..ef5de379acc03c218b6efee70c63670044920481
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt1_ee.py b/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..6b85580757192f0ca0ac63e2ac94fff17f76187d
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..07361f83a742b74d194105dede17ba0413365eae
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MD_ReStrip1_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MD'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt0_ee.py b/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..7b28c4a2887f295b6c039ad3fd7d25aba2e22200
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..ed2b94a6409ff94b79f4ad7d619b2ae7a3dda84f
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt1_ee.py b/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..99ebbe4849790a5e963842df096ccbbf4985061a
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..b908e3f877dd873f2b30847e431a6fac1ab00af0
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MU_ReStrip0_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "False"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt0_ee.py b/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt0_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..9f97b926f9f160324a03ac002626a3056883d049
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt0_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt0_mumu.py b/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt0_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..c1ab0a452ae42030b2542af12e7858f1c2276e0a
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt0_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'DST'
+os.environ['ISREDECAY']           = "False"
+os.environ['ISFILTERED']          = "False"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt1_ee.py b/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt1_ee.py
new file mode 100644
index 0000000000000000000000000000000000000000..29e67c98b96fcf57a931aed4cc2e19ccae584dab
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt1_ee.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "False"
+os.environ['DOELECTRON']          = "True"
diff --git a/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt1_mumu.py b/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt1_mumu.py
new file mode 100644
index 0000000000000000000000000000000000000000..bd4bd6b85f9e91becb40a71afccc6ed842a9f053
--- /dev/null
+++ b/Lb2pKll/steering_files/MC_2018_MU_ReStrip1_Filt1_mumu.py
@@ -0,0 +1,11 @@
+import os
+os.environ['ISRUN2']              = "True"
+os.environ['ISMC']                = "True"
+os.environ['YEAR']                = '2018'
+os.environ['INPUTTYPE']           = 'MDST'
+os.environ['ISREDECAY']           = "True"
+os.environ['ISFILTERED']          = "True"
+os.environ['MAGPOS']              = 'MU'
+os.environ['RESTRIP']             = "True"
+os.environ['DOMUON']              = "True"
+os.environ['DOELECTRON']          = "False"
diff --git a/Lb2pKll/steering_files/__init__.py b/Lb2pKll/steering_files/__init__.py
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391