From 51b4a15701e3c48b2ca4220ca867ec5e1dc5e7b7 Mon Sep 17 00:00:00 2001
From: Raphael van Laak <raphael.van.laak@cern.ch>
Date: Wed, 15 Mar 2023 16:44:08 +0100
Subject: [PATCH] This are decay files of (semi)leptonic decay of HNLs produced
 from Bu decay

---
 ...o2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec | 52 ++++++++++++++++++
 ...o2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec | 52 ++++++++++++++++++
 ...2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec | 52 ++++++++++++++++++
 ...2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec | 52 ++++++++++++++++++
 ...MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec | 52 ++++++++++++++++++
 ...MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec | 52 ++++++++++++++++++
 ...o2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec | 55 +++++++++++++++++++
 ...o2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec | 53 ++++++++++++++++++
 ...2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec | 55 +++++++++++++++++++
 ...2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec | 53 ++++++++++++++++++
 ...MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec | 55 +++++++++++++++++++
 ...MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec | 53 ++++++++++++++++++
 12 files changed, 636 insertions(+)
 create mode 100644 dkfiles/Bu_MuMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuXMajoranaNeutrino2EENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuXMajoranaNeutrino2EENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuXMajoranaNeutrino2MuENu,m=2000MeV,t=100ps,SS,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuXMajoranaNeutrino2MuENu,m=4000MeV,t=100ps,SS,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuXMajoranaNeutrino2MuMuNu,m=2000MeV,t=100ps,SS,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuXMajoranaNeutrino2MuMuNu,m=4000MeV,t=100ps,SS,DecProdCut.dec

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 000000000..0ed28de75
--- /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 000000000..e56216481
--- /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 000000000..f30e7a3ba
--- /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 000000000..4dfcc44b2
--- /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 000000000..e7d90fd29
--- /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 000000000..c7c79d53c
--- /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 000000000..63b50fb39
--- /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 000000000..9fa274e1a
--- /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 000000000..b6f33efad
--- /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 000000000..b1f0d174e
--- /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 000000000..65cfd4b51
--- /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 000000000..c9830c470
--- /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
+#
-- 
GitLab