diff --git a/dkfiles/Bu_D0K,KSpipi=TightCut,BELLE2018.dec b/dkfiles/Bu_D0K,KSpipi=TightCut,BELLE2018.dec
new file mode 100644
index 0000000000000000000000000000000000000000..46953b01ee0d425cc63444d464509d4f88e9ba21
--- /dev/null
+++ b/dkfiles/Bu_D0K,KSpipi=TightCut,BELLE2018.dec
@@ -0,0 +1,66 @@
+# EventType: 12165163
+#
+# Descriptor: [B+ -> K+ (anti-D0 -> (K_S0 -> pi+ pi-) pi+ pi-)]cc
+#
+# Documentation: B+ -> K+ D0b with D0b going to KsPiPi following the BELLE2018 model https://arxiv.org/abs/1804.06153, Ks -> pi+pi-
+# EndDocumentation
+#
+# NickName: Bu_D0K,KSpipi=TightCut,BELLE2018
+#
+# Cuts: LoKi::GenCutTool/TightCut
+#
+# InsertPythonCode:
+# #
+# from Configurables import LoKi__GenCutTool
+# from Gauss.Configuration import *
+# Generation().SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' )
+# tightCut = Generation().SignalRepeatedHadronization.TightCut
+# tightCut.Decay     = '^[B+ => ^(D~0 => ^(KS0 => ^pi+ ^pi-) ^pi+ ^pi-) ^K+]CC'
+# tightCut.Preambulo += [
+#     'GVZ = LoKi.GenVertices.PositionZ() ' ,
+#     'from GaudiKernel.SystemOfUnits import millimeter',
+#     'inAcc        = (in_range (0.005, GTHETA, 0.400))',
+#     'goodB        = (GP > 25000 * MeV) & (GPT > 1500 * MeV) & (GTIME > 0.05 * millimeter)',
+#     'goodD        = (GP > 10000 * MeV) & (GPT > 500 * MeV)',
+#     'goodKS       = (GP > 4000 * MeV) & (GPT > 250 * MeV)',
+#     'goodDDaugPi  = (GNINTREE (("pi+" == GABSID) & (GP > 1000 * MeV) & inAcc, 1) > 1.5)',
+#     'goodKsDaugPi = (GNINTREE (("pi+" == GABSID) & (GP > 1750 * MeV) & inAcc, 1) > 1.5)',
+#     'goodBachPi   = (GNINTREE (("K+" == GABSID) & (GP > 4000 * MeV) & (GPT > 400 * MeV) & inAcc, 1) > 0.5)'
+# ]
+# tightCut.Cuts      =    {
+#     '[B+]cc'         : 'goodB  & goodBachPi',
+#     '[D0]cc'         : 'goodD  & goodDDaugPi',
+#     '[KS0]cc'        : 'goodKS & goodKsDaugPi',
+#     '[pi+]cc'        : 'inAcc'
+#     }
+# EndInsertPythonCode
+#
+# CPUTime: < 1 min
+# PhysicsWG: B2OC
+# Tested: Yes
+# Responsible: Camille Normand
+# Email: camille.normand@cern.ch
+# Date: 20240219
+#
+Alias      MyD0      D0
+Alias      Myanti-D0 anti-D0
+Alias      MyKs      K_S0
+ChargeConj MyD0      Myanti-D0
+ChargeConj MyKs      MyKs
+#
+Decay B+sig
+  1.000    Myanti-D0  K+  BTODDALITZCPK 1.22 2.27 0.10;
+Enddecay
+CDecay B-sig
+#
+Decay MyD0
+  1.000    MyKs       pi-         pi+          D0TOKSPIPI;
+Enddecay
+CDecay Myanti-D0
+#
+Decay MyKs
+  1.000    pi+        pi-         PHSP;
+Enddecay
+#
+End
+#
diff --git a/dkfiles/Bu_D0pi,KSpipi=TightCut,BELLE2018.dec b/dkfiles/Bu_D0pi,KSpipi=TightCut,BELLE2018.dec
new file mode 100644
index 0000000000000000000000000000000000000000..e25d788f793eb86b5eb0b164a62ee29fa1a697bf
--- /dev/null
+++ b/dkfiles/Bu_D0pi,KSpipi=TightCut,BELLE2018.dec
@@ -0,0 +1,66 @@
+# EventType: 12165164
+#
+# Descriptor: [B+ -> pi+ (anti-D0 -> (K_S0 -> pi+ pi-) pi+ pi-)]cc
+#
+# Documentation: B+ -> pi+ D0b with D0b going to KsPiPi following the BELLE2018 model https://arxiv.org/abs/1804.06153, Ks -> pi+pi-
+# EndDocumentation
+#
+# NickName: Bu_D0pi,KSpipi=TightCut,BELLE2018
+#
+# Cuts: LoKi::GenCutTool/TightCut
+#
+# InsertPythonCode:
+# #
+# from Configurables import LoKi__GenCutTool
+# from Gauss.Configuration import *
+# Generation().SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' )
+# tightCut = Generation().SignalRepeatedHadronization.TightCut
+# tightCut.Decay     = '^[B+ => ^(D~0 => ^(KS0 => ^pi+ ^pi-) ^pi+ ^pi-) ^pi+]CC'
+# tightCut.Preambulo += [
+#     'GVZ = LoKi.GenVertices.PositionZ() ' ,
+#     'from GaudiKernel.SystemOfUnits import millimeter',
+#     'inAcc        = (in_range (0.005, GTHETA, 0.400))',
+#     'goodB        = (GP > 25000 * MeV) & (GPT > 1500 * MeV) & (GTIME > 0.05 * millimeter)',
+#     'goodD        = (GP > 10000 * MeV) & (GPT > 500 * MeV)',
+#     'goodKS       = (GP > 4000 * MeV) & (GPT > 250 * MeV)',
+#     'goodDDaugPi  = (GNINTREE (("pi+" == GABSID) & (GP > 1000 * MeV) & inAcc, 1) > 1.5)',
+#     'goodKsDaugPi = (GNINTREE (("pi+" == GABSID) & (GP > 1750 * MeV) & inAcc, 1) > 1.5)',
+#     'goodBachPi   = (GNINTREE (("pi+" == GABSID) & (GP > 4000 * MeV) & (GPT > 400 * MeV) & inAcc, 1) > 0.5)'
+# ]
+# tightCut.Cuts      =    {
+#     '[B+]cc'         : 'goodB  & goodBachPi',
+#     '[D0]cc'         : 'goodD  & goodDDaugPi',
+#     '[KS0]cc'        : 'goodKS & goodKsDaugPi',
+#     '[pi+]cc'        : 'inAcc'
+#     }
+# EndInsertPythonCode
+#
+# CPUTime: < 1 min
+# PhysicsWG: B2OC
+# Tested: Yes
+# Responsible: Camille Normand
+# Email: camille.normand@cern.ch
+# Date: 20240219
+#
+Alias      MyD0      D0
+Alias      Myanti-D0 anti-D0
+Alias      MyKs      K_S0
+ChargeConj MyD0      Myanti-D0
+ChargeConj MyKs      MyKs
+#
+Decay B+sig
+  1.000    Myanti-D0  pi+  BTODDALITZCPK 1.22 -1.206 0.005;
+Enddecay
+CDecay B-sig
+#
+Decay MyD0
+  1.000    MyKs       pi-         pi+          D0TOKSPIPI;
+Enddecay
+CDecay Myanti-D0
+#
+Decay MyKs
+  1.000    pi+        pi-         PHSP;
+Enddecay
+#
+End
+#