From 835b60ccd422212dbd2d1e3730fd20c5f4f81db5 Mon Sep 17 00:00:00 2001
From: lfantini <lisa.fantini@cern.ch>
Date: Fri, 6 May 2022 10:57:22 +0200
Subject: [PATCH] Bu/Bc-> lepton MajoranaNeutrino(->lepton Pi) 3GeV 0/20ps
 SS/OS

---
 ...rino2MuPiSS,m=3000MeV,t=0ps,DecProdCut.dec | 61 +++++++++++++++++++
 ...ino2MuPiSS,m=3000MeV,t=20ps,DecProdCut.dec | 61 +++++++++++++++++++
 ...trino2EPiSS,m=3000MeV,t=0ps,DecProdCut.dec | 61 +++++++++++++++++++
 ...rino2EPiSS,m=3000MeV,t=20ps,DecProdCut.dec | 61 +++++++++++++++++++
 ...rino2MuPiOS,m=3000MeV,t=0ps,DecProdCut.dec | 61 +++++++++++++++++++
 ...ino2MuPiOS,m=3000MeV,t=20ps,DecProdCut.dec | 61 +++++++++++++++++++
 ...rino2MuPiSS,m=3000MeV,t=0ps,DecProdCut.dec | 61 +++++++++++++++++++
 ...ino2MuPiSS,m=3000MeV,t=20ps,DecProdCut.dec | 61 +++++++++++++++++++
 ...trino2EPiOS,m=3000MeV,t=0ps,DecProdCut.dec | 61 +++++++++++++++++++
 ...rino2EPiOS,m=3000MeV,t=20ps,DecProdCut.dec | 61 +++++++++++++++++++
 ...trino2EPiSS,m=3000MeV,t=0ps,DecProdCut.dec | 61 +++++++++++++++++++
 ...rino2EPiSS,m=3000MeV,t=20ps,DecProdCut.dec | 61 +++++++++++++++++++
 12 files changed, 732 insertions(+)
 create mode 100644 dkfiles/Bc_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=0ps,DecProdCut.dec
 create mode 100644 dkfiles/Bc_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=20ps,DecProdCut.dec
 create mode 100644 dkfiles/Bc_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=0ps,DecProdCut.dec
 create mode 100644 dkfiles/Bc_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=20ps,DecProdCut.dec
 create mode 100644 dkfiles/Bu_EMajoranaNeutrino2MuPiOS,m=3000MeV,t=0ps,DecProdCut.dec
 create mode 100644 dkfiles/Bu_EMajoranaNeutrino2MuPiOS,m=3000MeV,t=20ps,DecProdCut.dec
 create mode 100644 dkfiles/Bu_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=0ps,DecProdCut.dec
 create mode 100644 dkfiles/Bu_EMajoranaNeutrino2MuPiSS,m=3000MeV,t=20ps,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuMajoranaNeutrino2EPiOS,m=3000MeV,t=0ps,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuMajoranaNeutrino2EPiOS,m=3000MeV,t=20ps,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=0ps,DecProdCut.dec
 create mode 100644 dkfiles/Bu_MuMajoranaNeutrino2EPiSS,m=3000MeV,t=20ps,DecProdCut.dec

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 000000000..2cd5e5fd7
--- /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 000000000..ae1d0d5a1
--- /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 000000000..8e404a02c
--- /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 000000000..ea3cc4394
--- /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 000000000..22e973533
--- /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 000000000..218e7f3f8
--- /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 000000000..698740606
--- /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 000000000..925a9acee
--- /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 000000000..eb57be873
--- /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 000000000..3591c0d19
--- /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 000000000..fb8cf893f
--- /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 000000000..ffbf8cb45
--- /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
+#
-- 
GitLab