diff --git a/dkfiles/Bc_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=0ps,DecProdCut.dec b/dkfiles/Bc_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=0ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..2cd5e5fd7571dd14365c4ffe22f24d205e68fa3e
--- /dev/null
+++ b/dkfiles/Bc_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=0ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bc+ -> E+ Majorana neutrino (-> PI- Mu+) SS
+#
+# EventType: 14113024
+#
+# Descriptor: [B_c+ -> e+ (Higgs0 -> pi- mu+) ]cc
+#
+# NickName: Bc_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=0ps,DecProdCut
+#
+# Production: BcVegPy
+#
+# Cuts: BcDaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     0.0000e+00     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B_c+ -> E+ Majorana neutrino (-> Pi-Mu+)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         mu+      pi-       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         mu-      pi+       PHSP;
+Enddecay
+#
+#
+Decay B_c+sig
+  1.000         e+    MyHNL         PHSP;
+Enddecay
+CDecay B_c-sig
+#
+#
+End
+#
diff --git a/dkfiles/Bc_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=20ps,DecProdCut.dec b/dkfiles/Bc_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=20ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..ae1d0d5a18fcbe5721d086540319b37def78bb22
--- /dev/null
+++ b/dkfiles/Bc_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=20ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bc+ -> E+ Majorana neutrino (-> PI- Mu+) SS
+#
+# EventType: 14113025
+#
+# Descriptor: [B_c+ -> e+ (Higgs0 -> pi- mu+) ]cc
+#
+# NickName: Bc_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=20ps,DecProdCut
+#
+# Production: BcVegPy
+#
+# Cuts: BcDaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     2.0000e-11     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B_c+ -> E+ Majorana neutrino (-> Pi-Mu+)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         mu+      pi-       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         mu-      pi+       PHSP;
+Enddecay
+#
+#
+Decay B_c+sig
+  1.000         e+    MyHNL         PHSP;
+Enddecay
+CDecay B_c-sig
+#
+#
+End
+#
diff --git a/dkfiles/Bc_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=0ps,DecProdCut.dec b/dkfiles/Bc_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=0ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..8e404a02cfea604c13af91d80ee04a0718a09c4b
--- /dev/null
+++ b/dkfiles/Bc_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=0ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bc+ -> Mu+ Majorana neutrino (-> PI- E+) SS
+#
+# EventType: 14113022
+#
+# Descriptor: [B_c+ -> mu+ (Higgs0 -> pi- e+) ]cc
+#
+# NickName: Bc_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=0ps,DecProdCut
+#
+# Production: BcVegPy
+#
+# Cuts: BcDaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     0.0000e+00     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B_c+ -> Mu+ Majorana neutrino (-> Pi-E+)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         e+      pi-       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         e-      pi+       PHSP;
+Enddecay
+#
+#
+Decay B_c+sig
+  1.000         mu+    MyHNL         PHSP;
+Enddecay
+CDecay B_c-sig
+#
+#
+End
+#
diff --git a/dkfiles/Bc_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=20ps,DecProdCut.dec b/dkfiles/Bc_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=20ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..ea3cc4394335574641c10a3fe14d37dcbe9d2e95
--- /dev/null
+++ b/dkfiles/Bc_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=20ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bc+ -> Mu+ Majorana neutrino (-> PI- E+) SS
+#
+# EventType: 14113023
+#
+# Descriptor: [B_c+ -> mu+ (Higgs0 -> pi- e+) ]cc
+#
+# NickName: Bc_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=20ps,DecProdCut
+#
+# Production: BcVegPy
+#
+# Cuts: BcDaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     2.0000e-11     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B_c+ -> Mu+ Majorana neutrino (-> Pi-E+)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         e+      pi-       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         e-      pi+       PHSP;
+Enddecay
+#
+#
+Decay B_c+sig
+  1.000         mu+    MyHNL         PHSP;
+Enddecay
+CDecay B_c-sig
+#
+#
+End
+#
diff --git a/dkfiles/Bu_EMajoranaNeutrino2MuPiOS,m=3000MeV,t=0ps,DecProdCut.dec b/dkfiles/Bu_EMajoranaNeutrino2MuPiOS,m=3000MeV,t=0ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..22e973533f2ef2786d8f6d47035e03e73a8ea81f
--- /dev/null
+++ b/dkfiles/Bu_EMajoranaNeutrino2MuPiOS,m=3000MeV,t=0ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bu+ -> E+ Majorana neutrino (-> PI+ Mu-) OS
+#
+# EventType: 12113080
+#
+# Descriptor: [B+ -> e+ (Higgs0 -> pi+ mu-) ]cc
+#
+# NickName: Bu_EMajoranaNeutrino2MuPiOS,m=3000MeV,t=0ps,DecProdCut
+#
+# 
+#
+# Cuts: DaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     0.0000e+00     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B+ -> E+ Majorana neutrino (-> Pi+Mu-)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         mu-      pi+       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         mu+      pi-       PHSP;
+Enddecay
+#
+#
+Decay B+sig
+  1.000         e+    MyHNL         PHSP;
+Enddecay
+CDecay B-sig
+#
+#
+End
+#
diff --git a/dkfiles/Bu_EMajoranaNeutrino2MuPiOS,m=3000MeV,t=20ps,DecProdCut.dec b/dkfiles/Bu_EMajoranaNeutrino2MuPiOS,m=3000MeV,t=20ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..218e7f3f8674e7ae1a9fc87a18cac2d499ea75bf
--- /dev/null
+++ b/dkfiles/Bu_EMajoranaNeutrino2MuPiOS,m=3000MeV,t=20ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bu+ -> E+ Majorana neutrino (-> PI+ Mu-) OS
+#
+# EventType: 12113000
+#
+# Descriptor: [B+ -> e+ (Higgs0 -> pi+ mu-) ]cc
+#
+# NickName: Bu_EMajoranaNeutrino2MuPiOS,m=3000MeV,t=20ps,DecProdCut
+#
+# 
+#
+# Cuts: DaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     2.0000e-11     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B+ -> E+ Majorana neutrino (-> Pi+Mu-)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         mu-      pi+       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         mu+      pi-       PHSP;
+Enddecay
+#
+#
+Decay B+sig
+  1.000         e+    MyHNL         PHSP;
+Enddecay
+CDecay B-sig
+#
+#
+End
+#
diff --git a/dkfiles/Bu_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=0ps,DecProdCut.dec b/dkfiles/Bu_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=0ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..698740606d965be7d15890b5caccd33f5a15bd32
--- /dev/null
+++ b/dkfiles/Bu_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=0ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bu+ -> E+ Majorana neutrino (-> PI- Mu+) SS
+#
+# EventType: 12113078
+#
+# Descriptor: [B+ -> e+ (Higgs0 -> pi- mu+) ]cc
+#
+# NickName: Bu_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=0ps,DecProdCut
+#
+# 
+#
+# Cuts: DaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     0.0000e+00     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B+ -> E+ Majorana neutrino (-> Pi-Mu+)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         mu+      pi-       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         mu-      pi+       PHSP;
+Enddecay
+#
+#
+Decay B+sig
+  1.000         e+    MyHNL         PHSP;
+Enddecay
+CDecay B-sig
+#
+#
+End
+#
diff --git a/dkfiles/Bu_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=20ps,DecProdCut.dec b/dkfiles/Bu_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=20ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..925a9acee41d2aa46159a0469e64926fb7091528
--- /dev/null
+++ b/dkfiles/Bu_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=20ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bu+ -> E+ Majorana neutrino (-> PI- Mu+) SS
+#
+# EventType: 12113079
+#
+# Descriptor: [B+ -> e+ (Higgs0 -> pi- mu+) ]cc
+#
+# NickName: Bu_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=20ps,DecProdCut
+#
+# 
+#
+# Cuts: DaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     2.0000e-11     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B+ -> E+ Majorana neutrino (-> Pi-Mu+)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         mu+      pi-       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         mu-      pi+       PHSP;
+Enddecay
+#
+#
+Decay B+sig
+  1.000         e+    MyHNL         PHSP;
+Enddecay
+CDecay B-sig
+#
+#
+End
+#
diff --git a/dkfiles/Bu_MuMajoranaNeutrino2EPiOS,m=3000MeV,t=0ps,DecProdCut.dec b/dkfiles/Bu_MuMajoranaNeutrino2EPiOS,m=3000MeV,t=0ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..eb57be8737c51f396923b9d2246493dbe43e7cab
--- /dev/null
+++ b/dkfiles/Bu_MuMajoranaNeutrino2EPiOS,m=3000MeV,t=0ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bu+ -> Mu+ Majorana neutrino (-> PI+ E-) OS
+#
+# EventType: 12113076
+#
+# Descriptor: [B+ -> mu+ (Higgs0 -> pi+ e-) ]cc
+#
+# NickName: Bu_MuMajoranaNeutrino2EPiOS,m=3000MeV,t=0ps,DecProdCut
+#
+# 
+#
+# Cuts: DaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     0.0000e+00     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B+ -> Mu+ Majorana neutrino (-> Pi+E-)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         e-      pi+       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         e+      pi-       PHSP;
+Enddecay
+#
+#
+Decay B+sig
+  1.000         mu+    MyHNL         PHSP;
+Enddecay
+CDecay B-sig
+#
+#
+End
+#
diff --git a/dkfiles/Bu_MuMajoranaNeutrino2EPiOS,m=3000MeV,t=20ps,DecProdCut.dec b/dkfiles/Bu_MuMajoranaNeutrino2EPiOS,m=3000MeV,t=20ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..3591c0d197d3d75a8936e125f5d66a516e1f56d5
--- /dev/null
+++ b/dkfiles/Bu_MuMajoranaNeutrino2EPiOS,m=3000MeV,t=20ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bu+ -> Mu+ Majorana neutrino (-> PI+ E-) OS
+#
+# EventType: 12113077
+#
+# Descriptor: [B+ -> mu+ (Higgs0 -> pi+ e-) ]cc
+#
+# NickName: Bu_MuMajoranaNeutrino2EPiOS,m=3000MeV,t=20ps,DecProdCut
+#
+# 
+#
+# Cuts: DaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     2.0000e-11     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B+ -> Mu+ Majorana neutrino (-> Pi+E-)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         e-      pi+       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         e+      pi-       PHSP;
+Enddecay
+#
+#
+Decay B+sig
+  1.000         mu+    MyHNL         PHSP;
+Enddecay
+CDecay B-sig
+#
+#
+End
+#
diff --git a/dkfiles/Bu_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=0ps,DecProdCut.dec b/dkfiles/Bu_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=0ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..fb8cf893f5180a26429295ebc9af633a0ce6537b
--- /dev/null
+++ b/dkfiles/Bu_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=0ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bu+ -> Mu+ Majorana neutrino (-> PI- E+) SS
+#
+# EventType: 12113074
+#
+# Descriptor: [B+ -> mu+ (Higgs0 -> pi- e+) ]cc
+#
+# NickName: Bu_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=0ps,DecProdCut
+#
+# 
+#
+# Cuts: DaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     0.0000e+00     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B+ -> Mu+ Majorana neutrino (-> Pi-E+)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         e+      pi-       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         e-      pi+       PHSP;
+Enddecay
+#
+#
+Decay B+sig
+  1.000         mu+    MyHNL         PHSP;
+Enddecay
+CDecay B-sig
+#
+#
+End
+#
diff --git a/dkfiles/Bu_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=20ps,DecProdCut.dec b/dkfiles/Bu_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=20ps,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..ffbf8cb45ae4e7534eefeb252042237130fa8feb
--- /dev/null
+++ b/dkfiles/Bu_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=20ps,DecProdCut.dec
@@ -0,0 +1,61 @@
+#
+# This is the decay file for the decay Bu+ -> Mu+ Majorana neutrino (-> PI- E+) SS
+#
+# EventType: 12113075
+#
+# Descriptor: [B+ -> mu+ (Higgs0 -> pi- e+) ]cc
+#
+# NickName: Bu_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=20ps,DecProdCut
+#
+# 
+#
+# Cuts: DaughtersInLHCb
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import Gauss, PrintMCTree, PrintMCDecayTreeTool, HistogramPersistencySvc, NTupleSvc, DumpHepMCDecay, DumpHepMCTree, GaussMonitor__CheckLifeTimeHepMC, GaussMonitor__CheckLifeTimeMC, GiGa, GiGaPhysListModular, GiGaHiggsParticles, GenerationToSimulation, PythiaProduction
+#
+#
+# ParticlePropertySvc().Particles = [ "H_10     87     25     0.0     3.00     2.0000e-11     Higgs0     25     0.000000e+000" ]
+# ApplicationMgr().ExtSvc    += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"] # H_10, H_20, H_30
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]#
+#
+#
+# EndInsertPythonCode
+#
+#
+#
+# Documentation:  Decay B+ -> Mu+ Majorana neutrino (-> Pi-E+)
+# EndDocumentation
+#
+# PhysicsWG: RD
+# Tested: Yes
+# Responsible: Lisa Fantini, Valeriia Lukashenko
+# Email: lisa.fantini@cern.ch, valeriia.lukashenko@cern.ch
+# Date: 20220506
+# CPUTime: <1min
+#
+Alias      MyHNL  Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+#
+Decay MyHNL
+  1.000         e+      pi-       PHSP;
+Enddecay
+Decay Myanti-HNL
+  1.000         e-      pi+       PHSP;
+Enddecay
+#
+#
+Decay B+sig
+  1.000         mu+    MyHNL         PHSP;
+Enddecay
+CDecay B-sig
+#
+#
+End
+#