diff --git a/dkfiles/Bu_MuMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..0ed28de7540afbe0d538e8e3bbaea4776ebb6c6f
--- /dev/null
+++ b/dkfiles/Bu_MuMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,52 @@
+# $Id: Bu_MuMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Anti-Nu_Mu E- E+) Mu+, with CC preserving same sign Muon
+#
+# EventType: 12513082
+#
+# Descriptor: [B+ -> (Higgs0 -> anti-nu_mu e- e+) mu+]cc
+#
+# NickName: Bu_MuMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 2.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Anti-Nu_Mu E- E+) Mu+, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+1	MyHNL	mu+	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 anti-nu_mu e- e+	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#
diff --git a/dkfiles/Bu_MuMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..e56216481caba70eaef1c55f675f1b0abfa778bd
--- /dev/null
+++ b/dkfiles/Bu_MuMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,52 @@
+# $Id: Bu_MuMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Anti-Nu_Mu E- E+) Mu+, with CC preserving same sign Muon
+#
+# EventType: 12513085
+#
+# Descriptor: [B+ -> (Higgs0 -> anti-nu_mu e- e+) mu+]cc
+#
+# NickName: Bu_MuMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 4.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Anti-Nu_Mu E- E+) Mu+, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+1	MyHNL	mu+	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 anti-nu_mu e- e+	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#
diff --git a/dkfiles/Bu_MuMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..f30e7a3ba214fe859a8fbcd17473a4dbb3628e09
--- /dev/null
+++ b/dkfiles/Bu_MuMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,52 @@
+# $Id: Bu_MuMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Mu+ E- Anti-Nu_E) Mu+, with CC preserving same sign Muon
+#
+# EventType: 12513081
+#
+# Descriptor: [B+ -> (Higgs0 -> mu+ e- anti-nu_e) mu+]cc
+#
+# NickName: Bu_MuMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 2.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Mu+ E- Anti-Nu_E) Mu+, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+1	MyHNL	mu+	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 mu+ e- anti-nu_e	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#
diff --git a/dkfiles/Bu_MuMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..4dfcc44b2257f15dab9ef28f8a1d1ad4a3968420
--- /dev/null
+++ b/dkfiles/Bu_MuMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,52 @@
+# $Id: Bu_MuMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Mu+ E- Anti-Nu_E) Mu+, with CC preserving same sign Muon
+#
+# EventType: 12513084
+#
+# Descriptor: [B+ -> (Higgs0 -> mu+ e- anti-nu_e) mu+]cc
+#
+# NickName: Bu_MuMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 4.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Mu+ E- Anti-Nu_E) Mu+, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+1	MyHNL	mu+	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 mu+ e- anti-nu_e	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#
diff --git a/dkfiles/Bu_MuMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..e7d90fd290a8b474bac79f2ac57b352313b52a5b
--- /dev/null
+++ b/dkfiles/Bu_MuMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,52 @@
+# $Id: Bu_MuMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Mu+ Mu- Anti-Nu_Mu) Mu+, with CC preserving same sign Muon
+#
+# EventType: 12513080
+#
+# Descriptor: [B+ -> (Higgs0 -> mu+ mu- anti-nu_mu) mu+]cc
+#
+# NickName: Bu_MuMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 2.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Mu+ Mu- Anti-Nu_Mu) Mu+, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+1	MyHNL	mu+	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 mu+ mu- anti-nu_mu	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#
diff --git a/dkfiles/Bu_MuMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..c7c79d53c8dc5f9e49c7fac6a90498b5f4f0670f
--- /dev/null
+++ b/dkfiles/Bu_MuMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,52 @@
+# $Id: Bu_MuMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Mu+ Mu- Anti-Nu_Mu) Mu+, with CC preserving same sign Muon
+#
+# EventType: 12513083
+#
+# Descriptor: [B+ -> (Higgs0 -> mu+ mu- anti-nu_mu) mu+]cc
+#
+# NickName: Bu_MuMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 4.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Mu+ Mu- Anti-Nu_Mu) Mu+, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+1	MyHNL	mu+	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 mu+ mu- anti-nu_mu	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#
diff --git a/dkfiles/Bu_MuXMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuXMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..63b50fb3989959c1a4257acf965b99646ab20fb4
--- /dev/null
+++ b/dkfiles/Bu_MuXMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,55 @@
+# $Id: Bu_MuXMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Anti-Nu_Mu E- E+) Mu+ X, with CC preserving same sign Muon
+#
+# EventType: 12773422
+#
+# Descriptor: [B+ -> (Higgs0 -> anti-nu_mu e- e+) mu+ X]cc
+#
+# NickName: Bu_MuXMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 2.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Anti-Nu_Mu E- E+) Mu+ X, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+0.0030	MyHNL	mu+	pi0	PHSP;
+0.0058	MyHNL	mu+	rho0	PHSP;
+0.3378	MyHNL	mu+	D0	PHSP;
+0.6534	MyHNL	mu+	D*0	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 anti-nu_mu e- e+	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#
diff --git a/dkfiles/Bu_MuXMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuXMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..9fa274e1a03000b9612072d3cfb646821cab5f2b
--- /dev/null
+++ b/dkfiles/Bu_MuXMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,53 @@
+# $Id: Bu_MuXMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Anti-Nu_Mu E- E+) Mu+ X, with CC preserving same sign Muon
+#
+# EventType: 12713402
+#
+# Descriptor: [B+ -> (Higgs0 -> anti-nu_mu e- e+) mu+ X]cc
+#
+# NickName: Bu_MuXMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 4.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Anti-Nu_Mu E- E+) Mu+ X, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+0.8748	MyHNL	mu+	pi0	PHSP;
+0.1252	MyHNL	mu+	rho0	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 anti-nu_mu e- e+	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#
diff --git a/dkfiles/Bu_MuXMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuXMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..b6f33efad144cba95a5cae00b1a211719e9e2eb6
--- /dev/null
+++ b/dkfiles/Bu_MuXMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,55 @@
+# $Id: Bu_MuXMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Mu+ E- Anti-Nu_E) Mu+ X, with CC preserving same sign Muon
+#
+# EventType: 12773421
+#
+# Descriptor: [B+ -> (Higgs0 -> mu+ e- anti-nu_e) mu+ X]cc
+#
+# NickName: Bu_MuXMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 2.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Mu+ E- Anti-Nu_E) Mu+ X, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+0.0030	MyHNL	mu+	pi0	PHSP;
+0.0058	MyHNL	mu+	rho0	PHSP;
+0.3378	MyHNL	mu+	D0	PHSP;
+0.6534	MyHNL	mu+	D*0	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 mu+ e- anti-nu_e	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#
diff --git a/dkfiles/Bu_MuXMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuXMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..b1f0d174e6d14ab0cbcfb3d2baa3d41e4ad703af
--- /dev/null
+++ b/dkfiles/Bu_MuXMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,53 @@
+# $Id: Bu_MuXMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Mu+ E- Anti-Nu_E) Mu+ X, with CC preserving same sign Muon
+#
+# EventType: 12713401
+#
+# Descriptor: [B+ -> (Higgs0 -> mu+ e- anti-nu_e) mu+ X]cc
+#
+# NickName: Bu_MuXMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 4.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Mu+ E- Anti-Nu_E) Mu+ X, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+0.8748	MyHNL	mu+	pi0	PHSP;
+0.1252	MyHNL	mu+	rho0	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 mu+ e- anti-nu_e	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#
diff --git a/dkfiles/Bu_MuXMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuXMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..65cfd4b519f2ea76d58c77efee8a3d9a79724dc3
--- /dev/null
+++ b/dkfiles/Bu_MuXMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,55 @@
+# $Id: Bu_MuXMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Mu+ Mu- Anti-Nu_Mu) Mu+ X, with CC preserving same sign Muon
+#
+# EventType: 12773420
+#
+# Descriptor: [B+ -> (Higgs0 -> mu+ mu- anti-nu_mu) mu+ X]cc
+#
+# NickName: Bu_MuXMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 2.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Mu+ Mu- Anti-Nu_Mu) Mu+ X, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+0.0030	MyHNL	mu+	pi0	PHSP;
+0.0058	MyHNL	mu+	rho0	PHSP;
+0.3378	MyHNL	mu+	D0	PHSP;
+0.6534	MyHNL	mu+	D*0	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 mu+ mu- anti-nu_mu	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#
diff --git a/dkfiles/Bu_MuXMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec b/dkfiles/Bu_MuXMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec
new file mode 100644
index 0000000000000000000000000000000000000000..c9830c47072d62e0ba4a77cd10673ae3779f9d84
--- /dev/null
+++ b/dkfiles/Bu_MuXMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec
@@ -0,0 +1,53 @@
+# $Id: Bu_MuXMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec   $
+#
+# This is the decay file for the decay B+ -> Majorana_Neutrino (-> Mu+ Mu- Anti-Nu_Mu) Mu+ X, with CC preserving same sign Muon
+#
+# EventType: 12713400
+#
+# Descriptor: [B+ -> (Higgs0 -> mu+ mu- anti-nu_mu) mu+ X]cc
+#
+# NickName: Bu_MuXMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut
+#
+# Cuts: DaughtersInLHCb
+#
+# InsertPythonCode:
+# from Gauss.Configuration import *
+# from Configurables import LHCb__ParticlePropertySvc as ParticlePropertySvc
+# from Configurables import GiGaPhysListModular, GiGaHiggsParticles
+#
+# ParticlePropertySvc().Particles = [ "H_10 87 25 0.0000 4.0000 1.0000e-10 Higgs0 25 0.0000" ]
+# ApplicationMgr().ExtSvc += [ ParticlePropertySvc() ]
+#
+# gigaHiggsPart = GiGaHiggsParticles()
+# gigaHiggsPart.Higgses = ["H_10"]
+# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHiggsPart ]
+# EndInsertPythonCode
+#
+#
+# Documentation: File for B+ -> Majorana_Neutrino (-> Mu+ Mu- Anti-Nu_Mu) Mu+ X, with CC preserving same sign Muon
+# EndDocumentation
+#
+# PhysicsWG: Exotica
+# Tested: Yes
+# Responsible: Raphael van Laak
+# Email: raphael.van.laak@cern.ch
+# Date:   20230307
+# CPUTime: <3min
+#
+Alias      MyHNL       Higgs0
+Alias      Myanti-HNL  Higgs0
+ChargeConj MyHNL Myanti-HNL
+#
+Decay B+sig
+0.8748	MyHNL	mu+	pi0	PHSP;
+0.1252	MyHNL	mu+	rho0	PHSP;
+Enddecay
+CDecay B-sig
+#
+Decay MyHNL
+1 mu+ mu- anti-nu_mu	PHSP;
+Enddecay
+CDecay Myanti-HNL
+#
+End
+#