From 9ec8d10c388a62a59c28e2e1c34ced5ced50ae5d Mon Sep 17 00:00:00 2001
From: Patrick Haworth Owen <patrick.haworth.owen@gmail.com>
Date: Thu, 13 Jul 2023 15:41:56 +0200
Subject: [PATCH 1/8] Added fixed DD decfiles for RD+

---
 .../Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec  | 306 ++++++++++++++++++
 .../Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec  | 305 +++++++++++++++++
 2 files changed, 611 insertions(+)
 create mode 100644 dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
 create mode 100644 dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec

diff --git a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
new file mode 100644
index 000000000..c75b3e426
--- /dev/null
+++ b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -0,0 +1,306 @@
+# EventType: 11995205
+#
+# Descriptor: [[B0] ==> (D*- -> (D~0 -> K+ mu- nu_mu~)) (D_s*+ -> gamma (D_s+ -> K+ K- pi+))]cc
+#
+# NickName: Bd_DD,DD=cocktail,D+muRDplusCutFixed
+# Cuts: LoKi::GenCutTool/TightCut
+#
+# InsertPythonCode:
+#
+# from Configurables import LoKi__GenCutTool
+# from Gauss.Configuration import *
+# gen = Generation()
+# gen.SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' )
+# tightCut = gen.SignalRepeatedHadronization.TightCut
+# tightCut.Decay = "[(Beauty) ==> ((Charm) ==> ^mu+ nu_mu ... ) ^(D- ==> ^K+ ^pi- ^pi-) {X} {X} {X} {X} {X} {X} {X} ]CC"
+#tightCut.Preambulo += [
+#  "from LoKiCore.functions import in_range"  ,
+#  "from GaudiKernel.SystemOfUnits import GeV, MeV"  ,
+#  "pipiKP     = GCHILD(GP,1) + GCHILD(GP,2) + GCHILD(GP,3)" ,
+#  "pipiKPT     = GCHILD(GPT,1) + GCHILD(GPT,2) + GCHILD(GPT,3)" ,
+#]
+#tightCut.Cuts      =    {
+# '[pi+]cc'   : "( GPT > 0 * MeV )" ,
+# '[K-]cc'   : "( GPT > 0 * MeV )" ,
+# '[mu+]cc'  : "(GP > 0* MeV) ",
+# '[D+]cc'   : "( pipiKP > 0 * MeV)"
+#    }
+# EndInsertPythonCode
+#
+# Documentation: Sum of B0 decaying to two D mesons, one D decaying semileptonically. Filtered for a final state with D+ mu.
+# 		 Background for R(D+) analysis.
+#                
+# EndDocumentation
+#
+# CPUTime: <1min
+# PhysicsWG: B2SL
+# Tested: Yes
+# Responsible: Patrick Owen
+# Email: patrick.haworth.owen@gmail.com
+# Date:   20160802
+#
+Alias        MyD0		D0
+Alias        MyAntiD0    	anti-D0
+ChargeConj   MyD0	 	MyAntiD0
+#
+Alias        MyD+		D+
+Alias	     MyD-		D-		
+ChargeConj   MyD+		MyD-
+#
+Alias        MySLD+		D+
+Alias	     MySLD-		D-		
+ChargeConj   MySLD+		MySLD-
+#
+Alias        MyD_s+		D_s+
+Alias	     MyD_s-		D_s-		
+ChargeConj   MyD_s+		MyD_s-		
+#
+Alias	     MyD_1+		D_1+
+Alias	     MyD_1-		D_1-
+ChargeConj   MyD_1+		MyD_1-
+#
+Alias	     MyD_s1+		D_s1+
+Alias	     MyD_s1-		D_s1-		
+ChargeConj   MyD_s1+		MyD_s1-
+#
+Alias        MyD_10		D_10
+Alias	     MyAntiD_10		anti-D_10
+ChargeConj   MyD_10		MyAntiD_10
+#
+Alias	     MyD(2S)0		D(2S)0
+Alias	     MyAntiD(2S)0	anti-D(2S)0	
+ChargeConj   MyD(2S)0		MyAntiD(2S)0
+#
+Alias	     MyD(2S)+		D(2S)+
+Alias	     MyD(2S)-		D(2S)-		
+ChargeConj   MyD(2S)+		MyD(2S)-
+#
+Alias	     MyD*0		D*0
+Alias	     MyAntiD*0		anti-D*0
+ChargeConj   MyD*0	 	MyAntiD*0
+#
+Alias	     MyD_2*0		D_2*0
+Alias	     MyAntiD_2*0 	anti-D_2*0	
+ChargeConj   MyD_2*0		MyAntiD_2*0
+#
+Alias	     MyD_0*0		D_0*0
+Alias 	     MyAntiD_0*0	anti-D_0*0		
+ChargeConj   MyD_0*0		MyAntiD_0*0
+#
+Alias	     MyD*(2S)0		D*(2S)0
+Alias	     MyAntiD*(2S)0	anti-D*(2S)0	
+ChargeConj   MyD*(2S)0		MyAntiD*(2S)0
+#
+Alias	     MyD*+		D*+
+Alias	     MyD*-		D*-		
+ChargeConj   MyD*+		MyD*-
+#
+Alias	     MySLD*+		D*+
+Alias	     MySLD*-		D*-		
+ChargeConj   MySLD*+		MySLD*-
+#
+Alias	     MyD_2*+		D_2*+
+Alias	     MyD_2*-		D_2*-		
+ChargeConj   MyD_2*+		MyD_2*-
+#
+Alias        MyD_s*+		D_s*+
+Alias	     MyD_s*-		D_s*-		
+ChargeConj   MyD_s*+		MyD_s*-
+#
+Alias        MyD_s2*+		D_s2*+
+Alias	     MyD_s2*-		D_s2*-		
+ChargeConj   MyD_s2*+		MyD_s2*-
+#
+Alias	     MyD_0*+		D_0*+
+Alias	     MyD_0*-		D_0*-		
+ChargeConj   MyD_0*+		MyD_0*-
+#
+Alias	     MyD_s0*+		D_s0*+
+Alias	     MyD_s0*-		D_s0*-
+ChargeConj   MyD_s0*+		MyD_s0*-
+#
+Alias	     MyD*(2S)+		D*(2S)+
+Alias	     MyD*(2S)-		D*(2S)-		
+ChargeConj   MyD*(2S)+		MyD*(2S)-
+#
+Alias	     MyD'_10		D'_10
+Alias	     MyAntiD'_10	anti-D'_10
+ChargeConj   MyD'_10		MyAntiD'_10
+#
+Alias	     MyD'_s1+		D'_s1+
+Alias	     MyD'_s1-		D'_s1-		
+ChargeConj   MyD'_s1+		MyD'_s1-
+#
+Alias	     MyD'_1+		D'_1+
+Alias	     MyD'_1-		D'_1-		
+ChargeConj   MyD'_1+		MyD'_1-
+#
+Decay B0sig
+ 0.00720      MyD-	MyD_s+			PHSP;
+ 0.00740      MyD_s*+	MyD-			SVS;
+ 0.00800      MyD*-	MyD_s+			SVS;
+ 0.01770      MyD*-	MyD_s*+			SVV_HELAMP 0.4904 0.0 0.7204 0.0 0.4904 0.0;
+ 0.00060      MyD'_1-	MyD_s+			SVS;
+ 0.00120      MyD'_1-	MyD_s*+			SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
+ 0.00120      MyD_1-	MyD_s+			SVS;
+ 0.00240      MyD_1-	MyD_s*+			SVV_HELAMP  0.48 0.0 0.734 0.0 0.48 0.0;
+ 0.00420      MyD_2*-	MyD_s+			STS;
+ 0.00400      MyD_2*-	MyD_s*+			PHSP;
+
+ 0.00160	      MyD-	MySLD+	K0		PHSP;
+ 0.00110	      MySLD*-	MyD+	K0		PHSP;
+ 0.00110	      MyD*-	MySLD+	K0		PHSP;
+ 0.00190	      MySLD-	MyD*+	K0		PHSP;
+ 0.00190	      MyD-	MySLD*+	K0		PHSP;
+ 0.00920	      MyD*-	MySLD*+   K0		PHSP;
+ 0.00150	      MyD-	MyD0	K+		PHSP;
+ 0.00630	      MyD-	MyD*0	K+		PHSP;
+ 0.00230	      MyD*-	MyD0	K+		PHSP;
+ 0.01120	      MyD*-	MyD*0	K+		PHSP;
+
+ 0.00250	      MyD-	MyD0	K*+		PHSP;
+ 0.00250	      MyD*-	MyD0	K*+		PHSP;
+ 0.00250	      MyD-	MyD*0	K*+		PHSP;
+ 0.00500	      MyD*-	MyD*0	K*+		PHSP;
+
+ 0.00250	      MyD-	MySLD+	K*0		PHSP;
+ 0.00250	      MySLD-	MyD*+	K*0		PHSP;
+ 0.00250	      MyD*-	MySLD+	K*0		PHSP;
+ 0.00500	      MySLD*-	MyD*+	K*0		PHSP;
+
+ 0.00180       	MyD_s+	MyD-	pi0		PHSP;
+ 0.00180	      MyD_s*+ 	MyD-	   pi0		PHSP;
+ 0.00300      	      MyD_s+	MyD-	   pi-   pi+	PHSP;
+ 0.00220      	      MyD_s+	MyD-	   pi0	 pi0	PHSP;
+ 0.00300      	      MyD_s*+	MyD-	   pi-	 pi+	PHSP;
+ 0.00220      	      MyD_s*+	MyD-	   pi0	 pi0	PHSP;
+
+ 0.000211	      MyD-	MySLD+			PHSP;
+ 0.000200	      MySLD*-	MyD+			SVS;
+ 0.000600	      MyD*+	MySLD-			SVS;
+ 0.000800	      MySLD*-	MyD*+			SVV_HELAMP 0.4904 0.0 0.7204 0.0 0.4904 0.0;
+	      
+ 0.00075              MyD-      MyD_s0*+        PHSP;
+ 0.00090              MyD*-     MyD_s0*+        SVS;
+ 0.00310              MyD_s1+	MyD-		SVS;
+ 0.01200              MyD*-     MyD_s1+         SVV_HELAMP 0.4904 0. 0.7204 0. 0.4904 0.0;
+	      
+ 0.00045              MyD-      MyD'_s1+        PHSP;
+ 0.00094              MyD*-     MyD'_s1+        PHSP;
+Enddecay
+CDecay anti-B0sig
+#
+SetLineshapePW MyD_1+ MyD*+ pi0 2
+SetLineshapePW MyD_1- MyD*- pi0 2
+SetLineshapePW MyD_10 MyD*+ pi- 2
+SetLineshapePW MyAntiD_10 MyD*- pi+ 2
+#
+SetLineshapePW MyD_2*+ MyD*+ pi0 2
+SetLineshapePW MyD_2*- MyD*- pi0 2
+SetLineshapePW MyD_2*0 MyD*+ pi- 2
+SetLineshapePW MyAntiD_2*0 MyD*- pi+ 2 
+#
+Decay MyD0
+ 0.0333          K-      mu+     nu_mu   PHOTOS ISGW2;
+ 0.0192         K*-     mu+     nu_mu  PHOTOS ISGW2;
+ 0.000076       K_1-    mu+     nu_mu  PHOTOS ISGW2;
+ 0.00110        K_2*-   mu+     nu_mu  PHOTOS ISGW2;
+ 0.00238        pi-     mu+     nu_mu  PHOTOS PHSP;
+ 0.00040        K-      pi0     mu+    nu_mu  PHOTOS PHSP;
+Enddecay
+CDecay MyAntiD0
+#
+Decay MyD*0
+ 0.619          MyD0    pi0             VSS;
+ 0.381          MyD0    gamma           VSP_PWAVE;
+Enddecay
+CDecay MyAntiD*0
+#
+Decay MyD_2*+
+ 0.10300        MyD*+   pi0             TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
+ 0.20900        MyD*0   pi+             TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
+ 0.22900        MyD+    pi0             TSS;
+ 0.45900        MyD0    pi+             TSS;
+Enddecay
+CDecay MyD_2*-
+#
+Decay MyD'_1+
+ 0.33333        MyD*+   pi0             VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
+ 0.66667        MyD*0   pi+             VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
+Enddecay
+CDecay MyD'_1+
+#
+Decay MyD*+
+ 0.67700        MyD0    pi+             VSS;
+ 0.30700        MyD+    pi0             VSS;
+ 0.01600        MyD+    gamma           VSP_PWAVE;
+Enddecay
+CDecay MyD*-
+#
+Decay MySLD*+
+ 0.67700        MyD0    pi+             VSS;
+ 0.30700        MySLD+    pi0             VSS;
+ 0.01600        MySLD+    gamma           VSP_PWAVE;
+Enddecay
+CDecay MySLD*-
+#
+Decay MyD_s*+
+ 0.93500        MyD_s+  gamma           VSP_PWAVE;
+ 0.05800        MyD_s+  pi0             VSS;
+Enddecay
+CDecay MyD_s*-
+#
+Decay MyD_s+
+ 0.0259 phi      mu+     nu_mu                      PHOTOS  ISGW2;
+ 0.0267 eta      mu+     nu_mu                      PHOTOS  ISGW2;
+ 0.0099 eta'     mu+     nu_mu                      PHOTOS  ISGW2;
+ 0.0037 K0       mu+     nu_mu                      PHOTOS  ISGW2;  
+ 0.0018 K*0      mu+     nu_mu                      PHOTOS  ISGW2;
+ 0.0020 f_0      mu+     nu_mu                      PHOTOS  ISGW2; 
+ 0.0059 mu+      nu_mu                              PHOTOS   SLN; 
+Enddecay
+CDecay MyD_s-
+#
+Decay MySLD+
+  0.05280       anti-K*0        mu+     nu_mu   PHOTOS ISGW2;
+  0.09200       anti-K0         mu+     nu_mu   PHOTOS ISGW2;
+  0.00277       anti-K_10       mu+     nu_mu   PHOTOS ISGW2;
+  0.00293       anti-K_2*0      mu+     nu_mu   PHOTOS ISGW2;
+  0.00292       K-   pi+        mu+     nu_mu   PHOTOS PHSP;
+  0.00120       anti-K0   pi0   mu+     nu_mu   PHOTOS PHSP;
+Enddecay
+CDecay MySLD-
+#
+Decay MyD+
+  0.09460       K-        pi+   pi+             PHOTOS D_DALITZ;
+Enddecay
+CDecay MyD-
+#
+Decay MyD_s1+
+ 0.48000        MyD_s*+  pi0             PARTWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
+ 0.18000        MyD_s+   gamma           VSP_PWAVE;
+ 0.04300        MyD_s+   pi+    pi-      PHSP;
+ 0.0037         MyD_s0*+ gamma           PHSP;
+Enddecay
+CDecay MyD_s1-
+#
+Decay MyD'_s1+
+ 0.42500        MyD*+     K0            VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
+ 0.42500        MyD*0     K+            VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
+ 0.02800        MyD+      pi-    K+     PHSP;
+ 0.03400        MyD+      K0            PHSP;
+Enddecay
+CDecay MyD'_s1-
+#
+Decay MyD_1+
+0.3333    MyD*+		pi0              VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
+0.6667    MyD*0 	pi+              VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
+Enddecay
+CDecay MyD_1-
+#
+Decay MyD_s0*+
+ 1.00000        MyD_s+  pi0     PHSP;
+Enddecay
+CDecay MyD_s0*-
+#
+End
diff --git a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
new file mode 100644
index 000000000..8c0bee7f6
--- /dev/null
+++ b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -0,0 +1,305 @@
+# EventType: 12995605
+#
+# Descriptor: [[B+] ==> (D*(2007)~0 -> (D~0 -> K+ mu- nu_mu~) pi0) (D_s*+ -> (D_s+ -> K+ K- pi+) gamma)]cc
+#
+# NickName: Bu_DD,DD=cocktail,D+muRDplusCutFixed
+# Cuts: LoKi::GenCutTool/TightCut
+#
+# InsertPythonCode:
+#
+# from Configurables import LoKi__GenCutTool
+# from Gauss.Configuration import *
+# gen = Generation()
+# gen.SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' )
+# tightCut = gen.SignalRepeatedHadronization.TightCut
+# tightCut.Decay = "[(Beauty) ==> ^((Charm) -> ^mu+ nu_mu ... ) ^(D- ==> ^K+ ^pi- ^pi- ) {X} {X} {X} {X} {X} {X} {X} ]CC"
+# tightCut.Preambulo += [
+# "from GaudiKernel.SystemOfUnits import mrad" ,
+# "FilterD = GNINTREE (GCHARM, HepMC.parents)",
+# "FromD   = 1 == FilterD",
+# "BCut = (GTHETA < 4000.0*mrad)"
+#
+#]
+#tightCut.Cuts    =    {
+# '[mu+]cc'     : "FromD",
+# '[B+]cc'	: "BCut"
+#}
+#
+# EndInsertPythonCode
+#
+# Documentation: Sum of B+ decaying to two D mesons, one D decaying semileptonically. Filtered for a final state with D+mu.
+# 		 Background for B2DmunuX analyses.
+#
+# EndDocumentation
+#
+# CPUTime: <1min
+# PhysicsWG: B2SL
+# Tested: Yes
+# Responsible: Patrick Owen
+# Email: patrick.haworth.owen@gmail.com
+# Date:   20160803
+#
+Alias        MyD0		D0
+Alias        MyAntiD0    	anti-D0
+ChargeConj   MyD0	 	MyAntiD0
+#
+Alias        MyD+		D+
+Alias	     MyD-		D-		
+ChargeConj   MyD+		MyD-
+#
+Alias        MySLD+		D+
+Alias	     MySLD-		D-		
+ChargeConj   MySLD+		MySLD-
+#
+Alias        MyD_s+		D_s+
+Alias	     MyD_s-		D_s-		
+ChargeConj   MyD_s+		MyD_s-		
+#
+Alias	     MyD_1+		D_1+
+Alias	     MyD_1-		D_1-
+ChargeConj   MyD_1+		MyD_1-
+#
+Alias	     MyD_s1+		D_s1+
+Alias	     MyD_s1-		D_s1-		
+ChargeConj   MyD_s1+		MyD_s1-
+#
+Alias        MyD_10		D_10
+Alias	     MyAntiD_10		anti-D_10
+ChargeConj   MyD_10		MyAntiD_10
+#
+Alias	     MyD(2S)0		D(2S)0
+Alias	     MyAntiD(2S)0	anti-D(2S)0	
+ChargeConj   MyD(2S)0		MyAntiD(2S)0
+#
+Alias	     MyD(2S)+		D(2S)+
+Alias	     MyD(2S)-		D(2S)-		
+ChargeConj   MyD(2S)+		MyD(2S)-
+#
+Alias	     MyD*0		D*0
+Alias	     MyAntiD*0		anti-D*0
+ChargeConj   MyD*0	 	MyAntiD*0
+#
+Alias	     MyD_2*0		D_2*0
+Alias	     MyAntiD_2*0 	anti-D_2*0	
+ChargeConj   MyD_2*0		MyAntiD_2*0
+#
+Alias	     MyD_0*0		D_0*0
+Alias 	     MyAntiD_0*0	anti-D_0*0		
+ChargeConj   MyD_0*0		MyAntiD_0*0
+#
+Alias	     MyD*(2S)0		D*(2S)0
+Alias	     MyAntiD*(2S)0	anti-D*(2S)0	
+ChargeConj   MyD*(2S)0		MyAntiD*(2S)0
+#
+Alias	     MyD*+		D*+
+Alias	     MyD*-		D*-		
+ChargeConj   MyD*+		MyD*-
+#
+Alias	     MySLD*+		D*+
+Alias	     MySLD*-		D*-		
+ChargeConj   MySLD*+		MySLD*-
+#
+Alias	     MyD_2*+		D_2*+
+Alias	     MyD_2*-		D_2*-		
+ChargeConj   MyD_2*+		MyD_2*-
+#
+Alias        MyD_s*+		D_s*+
+Alias	     MyD_s*-		D_s*-		
+ChargeConj   MyD_s*+		MyD_s*-
+#
+Alias        MyD_s2*+		D_s2*+
+Alias	     MyD_s2*-		D_s2*-		
+ChargeConj   MyD_s2*+		MyD_s2*-
+#
+Alias	     MyD_0*+		D_0*+
+Alias	     MyD_0*-		D_0*-		
+ChargeConj   MyD_0*+		MyD_0*-
+#
+Alias	     MyD_s0*+		D_s0*+
+Alias	     MyD_s0*-		D_s0*-
+ChargeConj   MyD_s0*+		MyD_s0*-
+#
+Alias	     MyD*(2S)+		D*(2S)+
+Alias	     MyD*(2S)-		D*(2S)-		
+ChargeConj   MyD*(2S)+		MyD*(2S)-
+#
+Alias	     MyD'_10		D'_10
+Alias	     MyAntiD'_10	anti-D'_10
+ChargeConj   MyD'_10		MyAntiD'_10
+#
+Alias	     MyD'_s1+		D'_s1+
+Alias	     MyD'_s1-		D'_s1-		
+ChargeConj   MyD'_s1+		MyD'_s1-
+#
+Alias	     MyD'_1+		D'_1+
+Alias	     MyD'_1-		D'_1-		
+ChargeConj   MyD'_1+		MyD'_1-
+#
+Decay B+sig
+ 0.00120      MyAntiD'_10	MyD_s+			SVS;
+ 0.00240      MyAntiD'_10	MyD_s*+			SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
+ 0.00240      MyAntiD_10	MyD_s+			SVS;
+ 0.00480      MyAntiD_10	MyD_s*+			SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
+ 0.00840      MyAntiD_2*0	MyD_s+			STS;
+ 0.00800      MyAntiD_2*0	MyD_s*+			PHSP;
+
+ 0.00160      MyAntiD0		MyD+	K0	PHSP;
+ 0.00380      MyAntiD0		MyD*+	K0	PHSP;
+ 0.00210      MyAntiD*0		MyD+	K0	PHSP;
+ 0.00920      MyAntiD*0		MyD*+	K0	PHSP;
+
+ 0.00250      MyAntiD0		MyD+	K*0		PHSP;
+ 0.00250      MyAntiD0		MyD*+	K*0		PHSP;
+ 0.00250      MyAntiD*0		MyD+	K*0		PHSP;
+ 0.00500      MyAntiD*0		MyD*+	K*0		PHSP;
+
+ 0.00360      MyD_s+	MyD-		pi+		PHSP;
+ 0.00370      MyD_s*+ 	MyD- 	    	pi+		PHSP;
+
+ 0.00330      MyD_s+	MyD-	   pi0   pi+	PHSP;
+ 0.00330      MyD_s*+	MyD-	   pi0	 pi+	PHSP;
+
+ 0.000400      MyAntiD0		MyD+			PHSP;
+ 0.000600      MyAntiD*0	MyD+			SVS;
+ 0.000400      MyD*+		MyAntiD0		SVS;
+ 0.000800      MyAntiD*0	MyD*+			SVV_HELAMP 0.283 0.0 0.932 0.0 0.228 0.0;
+
+
+ 0.00160	MySLD-	MyD+	K+	PHSP;
+ 0.00380	MyD-	MySLD*+	K+	PHSP;
+ 0.00210	MySLD*-	MyD+	K+	PHSP;
+ 0.00920	MyD*-	MySLD*+	K+	PHSP;
+
+ 0.00250	MySLD-	MyD+	K*+	PHSP;
+ 0.00250	MyD-	MySLD*+	K*+	PHSP;
+ 0.00250	MySLD*-	MyD+	K*+	PHSP;
+ 0.00500	MyD*-	MySLD*+	K*+	PHSP;      
+Enddecay
+CDecay B-sig
+#
+SetLineshapePW MyD_1+ MyD*+ pi0 2
+SetLineshapePW MyD_1- MyD*- pi0 2
+SetLineshapePW MyD_10 MyD*+ pi- 2
+SetLineshapePW MyAntiD_10 MyD*- pi+ 2
+#
+SetLineshapePW MyD_2*+ MyD*+ pi0 2
+SetLineshapePW MyD_2*- MyD*- pi0 2
+SetLineshapePW MyD_2*0 MyD*+ pi- 2
+SetLineshapePW MyAntiD_2*0 MyD*- pi+ 2
+#
+Decay MyD0
+ 0.0333          K-      mu+     nu_mu		PHOTOS ISGW2;
+ 0.0192         K*-     mu+     nu_mu  		PHOTOS ISGW2;
+ 0.000076       K_1-    mu+     nu_mu  		PHOTOS ISGW2;
+ 0.00110        K_2*-   mu+     nu_mu  		PHOTOS ISGW2;
+ 0.00238        pi-     mu+     nu_mu  		PHOTOS PHSP;
+ 0.00040        K-      pi0     mu+    nu_mu  	PHOTOS PHSP;
+Enddecay
+CDecay MyAntiD0
+#
+Decay MyD*0
+ 0.619          MyD0    pi0             VSS;
+ 0.381          MyD0    gamma           VSP_PWAVE;
+Enddecay
+CDecay MyAntiD*0
+#
+Decay MyD_2*+
+ 0.10300        MyD*+   pi0             TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
+ 0.20900        MyD*0   pi+             TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
+ 0.22900        MyD+    pi0             TSS;
+ 0.45900        MyD0    pi+             TSS;
+Enddecay
+CDecay MyD_2*-
+#
+Decay MyD'_1+
+ 0.33333        MyD*+   pi0             VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
+ 0.66667        MyD*0   pi+             VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
+Enddecay
+CDecay MyD'_1+
+#
+Decay MyD*+
+ 0.67700        MyD0    pi+             VSS;
+ 0.30700        MyD+    pi0             VSS;
+ 0.01600        MyD+    gamma           VSP_PWAVE;
+Enddecay
+CDecay MyD*-
+#
+Decay MySLD*+
+ 0.67700        MyD0    pi+             VSS;
+ 0.30700        MySLD+    pi0             VSS;
+ 0.01600        MySLD+    gamma           VSP_PWAVE;
+Enddecay
+CDecay MySLD*-
+#
+Decay MyD_s*+
+ 0.93500        MyD_s+  gamma           VSP_PWAVE;
+ 0.05800        MyD_s+  pi0             VSS;
+Enddecay
+CDecay MyD_s*-
+#
+Decay MyD_s+
+ 0.0259 phi      mu+     nu_mu                      PHOTOS  ISGW2;
+ 0.0267 eta      mu+     nu_mu                      PHOTOS  ISGW2;
+ 0.0099 eta'     mu+     nu_mu                      PHOTOS  ISGW2;
+ 0.0037 K0       mu+     nu_mu                      PHOTOS  ISGW2;  
+ 0.0018 K*0      mu+     nu_mu                      PHOTOS  ISGW2;
+ 0.0020 f_0      mu+     nu_mu                      PHOTOS  ISGW2; 
+ 0.0059 mu+      nu_mu                              PHOTOS   SLN; 
+Enddecay
+CDecay MyD_s-
+#
+Decay MySLD+
+  0.05280       anti-K*0        mu+     nu_mu   PHOTOS ISGW2;
+  0.09200       anti-K0         mu+     nu_mu   PHOTOS ISGW2;
+  0.00277       anti-K_10       mu+     nu_mu   PHOTOS ISGW2;
+  0.00293       anti-K_2*0      mu+     nu_mu   PHOTOS ISGW2;
+  0.00292       K-   pi+        mu+     nu_mu   PHOTOS PHSP;
+  0.00120       anti-K0   pi0   mu+     nu_mu   PHOTOS PHSP;
+Enddecay
+CDecay MySLD-
+#
+Decay MyD+
+  0.09460       K-        pi+   pi+             PHOTOS D_DALITZ;
+Enddecay
+CDecay MyD-
+#
+Decay MyD_s1+
+ 0.48000        MyD_s*+  pi0             PARTWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
+ 0.18000        MyD_s+   gamma           VSP_PWAVE;
+ 0.04300        MyD_s+   pi+    pi-      PHSP;
+ 0.0037         MyD_s0*+ gamma           PHSP;
+Enddecay
+CDecay MyD_s1-
+#
+Decay MyD'_s1+
+ 0.42500        MyD*+     K0            VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
+ 0.42500        MyD*0     K+            VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
+ 0.02800        MyD+      pi-    K+     PHSP;
+ 0.03400        MyD+      K0            PHSP;
+Enddecay
+CDecay MyD'_s1-
+#
+Decay MyD_1+
+0.3333    MyD*+ pi0                        VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
+0.6667    MyD*0 pi+                        VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
+Enddecay
+CDecay MyD_1-
+#
+Decay MyD_s0*+
+ 1.00000        MyD_s+  pi0     PHSP;
+Enddecay
+CDecay MyD_s0*-
+#
+Decay MyD'_10
+ 0.66667	MyD*+	pi-	PHOTOS VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
+Enddecay
+CDecay MyAntiD'_10
+#
+Decay MyD_2*0
+  0.40157      MyD+    pi-     PHOTOS TSS;
+  0.26509      MyD*+   pi-     PHOTOS TVS_PWAVE  0.0 0.0 1.0 0.0 0.0 0.0;
+Enddecay
+CDecay MyAntiD_2*0
+#
+End
+
-- 
GitLab


From bc8d98e4ee79efb8efb06bcddba52eeba4d48562 Mon Sep 17 00:00:00 2001
From: Patrick Haworth Owen <patrick.haworth.owen@gmail.com>
Date: Thu, 13 Jul 2023 16:03:04 +0200
Subject: [PATCH 2/8] Solve Loki cuts issues with fixed decfiles

---
 .../Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec  | 10 ++++----
 .../Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec  | 24 +++++++++----------
 2 files changed, 17 insertions(+), 17 deletions(-)

diff --git a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
index c75b3e426..4b60bb460 100644
--- a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
+++ b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -12,7 +12,7 @@
 # gen = Generation()
 # gen.SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' )
 # tightCut = gen.SignalRepeatedHadronization.TightCut
-# tightCut.Decay = "[(Beauty) ==> ((Charm) ==> ^mu+ nu_mu ... ) ^(D- ==> ^K+ ^pi- ^pi-) {X} {X} {X} {X} {X} {X} {X} ]CC"
+# tightCut.Decay = "[(Beauty) -> ((Charm) ==> ^mu+ nu_mu ... ) ^(D- ==> ^K+ ^pi- ^pi-) {X} {X} {X} {X} {X} {X} {X} ]CC"
 #tightCut.Preambulo += [
 #  "from LoKiCore.functions import in_range"  ,
 #  "from GaudiKernel.SystemOfUnits import GeV, MeV"  ,
@@ -20,10 +20,10 @@
 #  "pipiKPT     = GCHILD(GPT,1) + GCHILD(GPT,2) + GCHILD(GPT,3)" ,
 #]
 #tightCut.Cuts      =    {
-# '[pi+]cc'   : "( GPT > 0 * MeV )" ,
-# '[K-]cc'   : "( GPT > 0 * MeV )" ,
-# '[mu+]cc'  : "(GP > 0* MeV) ",
-# '[D+]cc'   : "( pipiKP > 0 * MeV)"
+# '[pi+]cc'   : "in_range( 0.010 , GTHETA , 0.400 )& ( GPT > 150 * MeV )" ,
+# '[K-]cc'   : "in_range( 0.010 , GTHETA , 0.400 ) & ( GPT > 150 * MeV )" ,
+# '[mu+]cc'  : "in_range( 0.010 , GTHETA , 0.400 ) & (GP > 2500* MeV)",
+# '[D+]cc'   : "( pipiKP > 15000 * MeV ) & (pipiKPT > 2300 * MeV)"
 #    }
 # EndInsertPythonCode
 #
diff --git a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
index 8c0bee7f6..ba8cb4e06 100644
--- a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
+++ b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -12,19 +12,19 @@
 # gen = Generation()
 # gen.SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' )
 # tightCut = gen.SignalRepeatedHadronization.TightCut
-# tightCut.Decay = "[(Beauty) ==> ^((Charm) -> ^mu+ nu_mu ... ) ^(D- ==> ^K+ ^pi- ^pi- ) {X} {X} {X} {X} {X} {X} {X} ]CC"
-# tightCut.Preambulo += [
-# "from GaudiKernel.SystemOfUnits import mrad" ,
-# "FilterD = GNINTREE (GCHARM, HepMC.parents)",
-# "FromD   = 1 == FilterD",
-# "BCut = (GTHETA < 4000.0*mrad)"
-#
+# tightCut.Decay = "[(Beauty) ==> ((Charm) -> ^mu+ nu_mu ... ) ^(D- ==> ^K+ ^pi- ^pi- ) {X} {X} {X} {X} {X} {X} {X} ]CC"
+#tightCut.Preambulo += [
+#  "from LoKiCore.functions import in_range"  ,
+#  "from GaudiKernel.SystemOfUnits import GeV, MeV"  ,
+#  "pipiKP     = GCHILD(GP,1) + GCHILD(GP,2) + GCHILD(GP,3)" ,
+#  "pipiKPT     = GCHILD(GPT,1) + GCHILD(GPT,2) + GCHILD(GPT,3)" ,
 #]
-#tightCut.Cuts    =    {
-# '[mu+]cc'     : "FromD",
-# '[B+]cc'	: "BCut"
-#}
-#
+#tightCut.Cuts      =    {
+# '[pi+]cc'   : "in_range( 0.010 , GTHETA , 0.400 )& ( GPT > 150 * MeV )" ,
+# '[K-]cc'   : "in_range( 0.010 , GTHETA , 0.400 ) & ( GPT > 150 * MeV )" ,
+# '[mu+]cc'  : "in_range( 0.010 , GTHETA , 0.400 ) & (GP > 2500* MeV)",
+# '[D+]cc'   : "( pipiKP > 15000 * MeV ) & (pipiKPT > 2300 * MeV)"
+#    }
 # EndInsertPythonCode
 #
 # Documentation: Sum of B+ decaying to two D mesons, one D decaying semileptonically. Filtered for a final state with D+mu.
-- 
GitLab


From c9b38511a7a0fb966beac543072758dbfddced1e Mon Sep 17 00:00:00 2001
From: Patrick Haworth Owen <patrick.haworth.owen@gmail.com>
Date: Thu, 13 Jul 2023 17:16:44 +0200
Subject: [PATCH 3/8] fixed photos/inclusive clash

---
 dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
index 4b60bb460..1084de34b 100644
--- a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
+++ b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -12,7 +12,7 @@
 # gen = Generation()
 # gen.SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' )
 # tightCut = gen.SignalRepeatedHadronization.TightCut
-# tightCut.Decay = "[(Beauty) -> ((Charm) ==> ^mu+ nu_mu ... ) ^(D- ==> ^K+ ^pi- ^pi-) {X} {X} {X} {X} {X} {X} {X} ]CC"
+# tightCut.Decay = "[(Beauty) -> ((Charm) -> ^mu+ nu_mu ... ) ^(D- ==> ^K+ ^pi- ^pi-) {X} {X} {X} {X} {X} {X} {X} ]CC"
 #tightCut.Preambulo += [
 #  "from LoKiCore.functions import in_range"  ,
 #  "from GaudiKernel.SystemOfUnits import GeV, MeV"  ,
-- 
GitLab


From 19a6f1ebba532ea49c541c331e53ee4fcd27a243 Mon Sep 17 00:00:00 2001
From: Patrick Haworth Owen <patrick.haworth.owen@gmail.com>
Date: Thu, 13 Jul 2023 17:40:07 +0200
Subject: [PATCH 4/8] update dates

---
 dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec | 2 +-
 dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
index 1084de34b..86c71da5f 100644
--- a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
+++ b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -37,7 +37,7 @@
 # Tested: Yes
 # Responsible: Patrick Owen
 # Email: patrick.haworth.owen@gmail.com
-# Date:   20160802
+# Date:   20230712
 #
 Alias        MyD0		D0
 Alias        MyAntiD0    	anti-D0
diff --git a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
index ba8cb4e06..62ce97414 100644
--- a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
+++ b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -37,7 +37,7 @@
 # Tested: Yes
 # Responsible: Patrick Owen
 # Email: patrick.haworth.owen@gmail.com
-# Date:   20160803
+# Date:   20230712
 #
 Alias        MyD0		D0
 Alias        MyAntiD0    	anti-D0
-- 
GitLab


From e3cfa4ffc416ae748a487cad350b43795acd8f38 Mon Sep 17 00:00:00 2001
From: Patrick Haworth Owen <patrick.haworth.owen@gmail.com>
Date: Fri, 14 Jul 2023 10:03:43 +0200
Subject: [PATCH 5/8] remove defuct D0 decays to save timing

---
 .../Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec  | 37 +++++++++----------
 .../Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec  | 18 ++++-----
 2 files changed, 26 insertions(+), 29 deletions(-)

diff --git a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
index 86c71da5f..7ff6aa766 100644
--- a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
+++ b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -32,7 +32,7 @@
 #                
 # EndDocumentation
 #
-# CPUTime: <1min
+# CPUTime: <2min
 # PhysicsWG: B2SL
 # Tested: Yes
 # Responsible: Patrick Owen
@@ -138,8 +138,8 @@ ChargeConj   MyD'_1+		MyD'_1-
 Decay B0sig
  0.00720      MyD-	MyD_s+			PHSP;
  0.00740      MyD_s*+	MyD-			SVS;
- 0.00800      MyD*-	MyD_s+			SVS;
- 0.01770      MyD*-	MyD_s*+			SVV_HELAMP 0.4904 0.0 0.7204 0.0 0.4904 0.0;
+ 0.00250      MyD*-	MyD_s+			SVS;
+ 0.00570      MyD*-	MyD_s*+			SVV_HELAMP 0.4904 0.0 0.7204 0.0 0.4904 0.0;
  0.00060      MyD'_1-	MyD_s+			SVS;
  0.00120      MyD'_1-	MyD_s*+			SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
  0.00120      MyD_1-	MyD_s+			SVS;
@@ -149,24 +149,24 @@ Decay B0sig
 
  0.00160	      MyD-	MySLD+	K0		PHSP;
  0.00110	      MySLD*-	MyD+	K0		PHSP;
- 0.00110	      MyD*-	MySLD+	K0		PHSP;
- 0.00190	      MySLD-	MyD*+	K0		PHSP;
+ 0.00040	      MyD*-	MySLD+	K0		PHSP;
+ 0.00070	      MySLD-	MyD*+	K0		PHSP;
  0.00190	      MyD-	MySLD*+	K0		PHSP;
- 0.00920	      MyD*-	MySLD*+   K0		PHSP;
+ 0.00320	      MyD*-	MySLD*+   K0		PHSP;
  0.00150	      MyD-	MyD0	K+		PHSP;
  0.00630	      MyD-	MyD*0	K+		PHSP;
- 0.00230	      MyD*-	MyD0	K+		PHSP;
- 0.01120	      MyD*-	MyD*0	K+		PHSP;
+ 0.00070	      MyD*-	MyD0	K+		PHSP;
+ 0.00420	      MyD*-	MyD*0	K+		PHSP;
 
  0.00250	      MyD-	MyD0	K*+		PHSP;
- 0.00250	      MyD*-	MyD0	K*+		PHSP;
+ 0.00080	      MyD*-	MyD0	K*+		PHSP;
  0.00250	      MyD-	MyD*0	K*+		PHSP;
- 0.00500	      MyD*-	MyD*0	K*+		PHSP;
+ 0.00150	      MyD*-	MyD*0	K*+		PHSP;
 
  0.00250	      MyD-	MySLD+	K*0		PHSP;
- 0.00250	      MySLD-	MyD*+	K*0		PHSP;
- 0.00250	      MyD*-	MySLD+	K*0		PHSP;
- 0.00500	      MySLD*-	MyD*+	K*0		PHSP;
+ 0.00080	      MySLD-	MyD*+	K*0		PHSP;
+ 0.00080	      MyD*-	MySLD+	K*0		PHSP;
+ 0.00150	      MySLD*-	MyD*+	K*0		PHSP;
 
  0.00180       	MyD_s+	MyD-	pi0		PHSP;
  0.00180	      MyD_s*+ 	MyD-	   pi0		PHSP;
@@ -177,16 +177,16 @@ Decay B0sig
 
  0.000211	      MyD-	MySLD+			PHSP;
  0.000200	      MySLD*-	MyD+			SVS;
- 0.000600	      MyD*+	MySLD-			SVS;
- 0.000800	      MySLD*-	MyD*+			SVV_HELAMP 0.4904 0.0 0.7204 0.0 0.4904 0.0;
+ 0.000200	      MyD*+	MySLD-			SVS;
+ 0.000240	      MySLD*-	MyD*+			SVV_HELAMP 0.4904 0.0 0.7204 0.0 0.4904 0.0;
 	      
  0.00075              MyD-      MyD_s0*+        PHSP;
- 0.00090              MyD*-     MyD_s0*+        SVS;
+ 0.00030              MyD*-     MyD_s0*+        SVS;
  0.00310              MyD_s1+	MyD-		SVS;
- 0.01200              MyD*-     MyD_s1+         SVV_HELAMP 0.4904 0. 0.7204 0. 0.4904 0.0;
+ 0.00400              MyD*-     MyD_s1+         SVV_HELAMP 0.4904 0. 0.7204 0. 0.4904 0.0;
 	      
  0.00045              MyD-      MyD'_s1+        PHSP;
- 0.00094              MyD*-     MyD'_s1+        PHSP;
+ 0.00034              MyD*-     MyD'_s1+        PHSP;
 Enddecay
 CDecay anti-B0sig
 #
@@ -231,7 +231,6 @@ Enddecay
 CDecay MyD'_1+
 #
 Decay MyD*+
- 0.67700        MyD0    pi+             VSS;
  0.30700        MyD+    pi0             VSS;
  0.01600        MyD+    gamma           VSP_PWAVE;
 Enddecay
diff --git a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
index 62ce97414..3aa6489f1 100644
--- a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
+++ b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -144,15 +144,14 @@ Decay B+sig
  0.00800      MyAntiD_2*0	MyD_s*+			PHSP;
 
  0.00160      MyAntiD0		MyD+	K0	PHSP;
- 0.00380      MyAntiD0		MyD*+	K0	PHSP;
+ 0.00140      MyAntiD0		MyD*+	K0	PHSP;
  0.00210      MyAntiD*0		MyD+	K0	PHSP;
- 0.00920      MyAntiD*0		MyD*+	K0	PHSP;
+ 0.00320      MyAntiD*0		MyD*+	K0	PHSP;
 
  0.00250      MyAntiD0		MyD+	K*0		PHSP;
- 0.00250      MyAntiD0		MyD*+	K*0		PHSP;
+ 0.00080      MyAntiD0		MyD*+	K*0		PHSP;
  0.00250      MyAntiD*0		MyD+	K*0		PHSP;
- 0.00500      MyAntiD*0		MyD*+	K*0		PHSP;
-
+ 0.00150      MyAntiD*0		MyD*+	K*0		PHSP;
  0.00360      MyD_s+	MyD-		pi+		PHSP;
  0.00370      MyD_s*+ 	MyD- 	    	pi+		PHSP;
 
@@ -161,19 +160,19 @@ Decay B+sig
 
  0.000400      MyAntiD0		MyD+			PHSP;
  0.000600      MyAntiD*0	MyD+			SVS;
- 0.000400      MyD*+		MyAntiD0		SVS;
- 0.000800      MyAntiD*0	MyD*+			SVV_HELAMP 0.283 0.0 0.932 0.0 0.228 0.0;
+ 0.000150      MyD*+		MyAntiD0		SVS;
+ 0.000250      MyAntiD*0	MyD*+			SVV_HELAMP 0.283 0.0 0.932 0.0 0.228 0.0;
 
 
  0.00160	MySLD-	MyD+	K+	PHSP;
  0.00380	MyD-	MySLD*+	K+	PHSP;
  0.00210	MySLD*-	MyD+	K+	PHSP;
- 0.00920	MyD*-	MySLD*+	K+	PHSP;
+ 0.00320	MyD*-	MySLD*+	K+	PHSP;
 
  0.00250	MySLD-	MyD+	K*+	PHSP;
  0.00250	MyD-	MySLD*+	K*+	PHSP;
  0.00250	MySLD*-	MyD+	K*+	PHSP;
- 0.00500	MyD*-	MySLD*+	K*+	PHSP;      
+ 0.00150	MyD*-	MySLD*+	K*+	PHSP;      
 Enddecay
 CDecay B-sig
 #
@@ -218,7 +217,6 @@ Enddecay
 CDecay MyD'_1+
 #
 Decay MyD*+
- 0.67700        MyD0    pi+             VSS;
  0.30700        MyD+    pi0             VSS;
  0.01600        MyD+    gamma           VSP_PWAVE;
 Enddecay
-- 
GitLab


From c4bc15e819fd287c49b31cb6361a8797be1b4bfa Mon Sep 17 00:00:00 2001
From: Patrick Haworth Owen <patrick.haworth.owen@gmail.com>
Date: Mon, 17 Jul 2023 11:26:54 +0200
Subject: [PATCH 6/8] set old DD decfiles to be buggy

---
 dkfiles/Bd_DD,DD=cocktail,D+muRDplusCut.dec | 310 --------------------
 dkfiles/Bu_DD,DD=cocktail,D+muRDplusCut.dec | 309 -------------------
 doc/table_obsolete.sql                      |   2 +
 3 files changed, 2 insertions(+), 619 deletions(-)
 delete mode 100644 dkfiles/Bd_DD,DD=cocktail,D+muRDplusCut.dec
 delete mode 100644 dkfiles/Bu_DD,DD=cocktail,D+muRDplusCut.dec

diff --git a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCut.dec b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCut.dec
deleted file mode 100644
index c4276dcde..000000000
--- a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCut.dec
+++ /dev/null
@@ -1,310 +0,0 @@
-# EventType: 11995204
-#
-# Descriptor: [[B0] ==> (D*- -> (D~0 -> K+ mu- nu_mu~)) (D_s*+ -> gamma (D_s+ -> K+ K- pi+))]cc
-#
-# NickName: Bd_DD,DD=cocktail,D+muRDplusCut
-# Cuts: LoKi::GenCutTool/TightCut
-#
-# InsertPythonCode:
-#
-# from Configurables import LoKi__GenCutTool
-# from Gauss.Configuration import *
-# gen = Generation()
-# gen.SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' )
-# tightCut = gen.SignalRepeatedHadronization.TightCut
-# tightCut.Decay = "[B0 ==> ((Charm) -> ^mu+ nu_mu ... ) ^(D- -> ^K+ ^pi- ^pi-) {X} {X} {X} {X} {X} {X} {X} ]CC"
-#tightCut.Preambulo += [
-#  "from LoKiCore.functions import in_range"  ,
-#  "from GaudiKernel.SystemOfUnits import GeV, MeV"  ,
-#  "pipiKP     = GCHILD(GP,1) + GCHILD(GP,2) + GCHILD(GP,3)" ,
-#  "pipiKPT     = GCHILD(GPT,1) + GCHILD(GPT,2) + GCHILD(GPT,3)" ,
-#]
-#tightCut.Cuts      =    {
-# '[pi+]cc'   : " in_range( 0.010 , GTHETA , 0.400 )& ( GPT > 150 * MeV )" ,
-# '[K-]cc'   : " in_range( 0.010 , GTHETA , 0.400 ) & ( GPT > 150 * MeV )" ,
-# '[mu+]cc'  : " in_range( 0.010 , GTHETA , 0.400 ) & (GP > 2500* MeV) ",
-# '[D+]cc'   : "( pipiKP > 15000 * MeV ) & (pipiKPT > 2300 * MeV)"
-#    }
-# EndInsertPythonCode
-#
-# Documentation: Sum of B0 decaying to two D mesons, one D decaying semileptonically. Filtered for a final state with D+ mu.
-# 		 Background for R(D+) analysis.
-#                
-# EndDocumentation
-#
-# CPUTime: <1min
-# PhysicsWG: B2SL
-# Tested: Yes
-# Responsible: Patrick Owen
-# Email: patrick.haworth.owen@gmail.com
-# Date:   20160802
-#
-Alias        MyD0		D0
-Alias        MyAntiD0    	anti-D0
-ChargeConj   MyD0	 	MyAntiD0
-#
-Alias        MyD+		D+
-Alias	     MyD-		D-		
-ChargeConj   MyD+		MyD-
-#
-Alias        MySLD+		D+
-Alias	     MySLD-		D-		
-ChargeConj   MySLD+		MySLD-
-#
-Alias        MyD_s+		D_s+
-Alias	     MyD_s-		D_s-		
-ChargeConj   MyD_s+		MyD_s-		
-#
-Alias	     MyD_1+		D_1+
-Alias	     MyD_1-		D_1-
-ChargeConj   MyD_1+		MyD_1-
-#
-Alias	     MyD_s1+		D_s1+
-Alias	     MyD_s1-		D_s1-		
-ChargeConj   MyD_s1+		MyD_s1-
-#
-Alias        MyD_10		D_10
-Alias	     MyAntiD_10		anti-D_10
-ChargeConj   MyD_10		MyAntiD_10
-#
-Alias	     MyD(2S)0		D(2S)0
-Alias	     MyAntiD(2S)0	anti-D(2S)0	
-ChargeConj   MyD(2S)0		MyAntiD(2S)0
-#
-Alias	     MyD(2S)+		D(2S)+
-Alias	     MyD(2S)-		D(2S)-		
-ChargeConj   MyD(2S)+		MyD(2S)-
-#
-Alias	     MyD*0		D*0
-Alias	     MyAntiD*0		anti-D*0
-ChargeConj   MyD*0	 	MyAntiD*0
-#
-Alias	     MyD_2*0		D_2*0
-Alias	     MyAntiD_2*0 	anti-D_2*0	
-ChargeConj   MyD_2*0		MyAntiD_2*0
-#
-Alias	     MyD_0*0		D_0*0
-Alias 	     MyAntiD_0*0	anti-D_0*0		
-ChargeConj   MyD_0*0		MyAntiD_0*0
-#
-Alias	     MyD*(2S)0		D*(2S)0
-Alias	     MyAntiD*(2S)0	anti-D*(2S)0	
-ChargeConj   MyD*(2S)0		MyAntiD*(2S)0
-#
-Alias	     MyD*+		D*+
-Alias	     MyD*-		D*-		
-ChargeConj   MyD*+		MyD*-
-#
-Alias	     MySLD*+		D*+
-Alias	     MySLD*-		D*-		
-ChargeConj   MySLD*+		MySLD*-
-#
-Alias	     MyD_2*+		D_2*+
-Alias	     MyD_2*-		D_2*-		
-ChargeConj   MyD_2*+		MyD_2*-
-#
-Alias        MyD_s*+		D_s*+
-Alias	     MyD_s*-		D_s*-		
-ChargeConj   MyD_s*+		MyD_s*-
-#
-Alias        MyD_s2*+		D_s2*+
-Alias	     MyD_s2*-		D_s2*-		
-ChargeConj   MyD_s2*+		MyD_s2*-
-#
-Alias	     MyD_0*+		D_0*+
-Alias	     MyD_0*-		D_0*-		
-ChargeConj   MyD_0*+		MyD_0*-
-#
-Alias	     MyD_s0*+		D_s0*+
-Alias	     MyD_s0*-		D_s0*-
-ChargeConj   MyD_s0*+		MyD_s0*-
-#
-Alias	     MyD*(2S)+		D*(2S)+
-Alias	     MyD*(2S)-		D*(2S)-		
-ChargeConj   MyD*(2S)+		MyD*(2S)-
-#
-Alias	     MyD'_10		D'_10
-Alias	     MyAntiD'_10	anti-D'_10
-ChargeConj   MyD'_10		MyAntiD'_10
-#
-Alias	     MyD'_s1+		D'_s1+
-Alias	     MyD'_s1-		D'_s1-		
-ChargeConj   MyD'_s1+		MyD'_s1-
-#
-Alias	     MyD'_1+		D'_1+
-Alias	     MyD'_1-		D'_1-		
-ChargeConj   MyD'_1+		MyD'_1-
-#
-Decay B0sig
- 0.00720      MyD-	MyD_s+			PHSP;
- 0.00740      MyD_s*+	MyD-			SVS;
- 0.00800      MyD*-	MyD_s+			SVS;
- 0.01770      MyD*-	MyD_s*+			SVV_HELAMP 0.4904 0.0 0.7204 0.0 0.4904 0.0;
- 0.00060      MyD'_1-	MyD_s+			SVS;
- 0.00120      MyD'_1-	MyD_s*+			SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
- 0.00120      MyD_1-	MyD_s+			SVS;
- 0.00240      MyD_1-	MyD_s*+			SVV_HELAMP  0.48 0.0 0.734 0.0 0.48 0.0;
- 0.00420      MyD_2*-	MyD_s+			STS;
- 0.00400      MyD_2*-	MyD_s*+			PHSP;
-
- 0.00160	      MyD-	MySLD+	K0		PHSP;
- 0.00110	      MySLD*-	MyD+	K0		PHSP;
- 0.00110	      MyD*-	MySLD+	K0		PHSP;
- 0.00190	      MySLD-	MyD*+	K0		PHSP;
- 0.00190	      MyD-	MySLD*+	K0		PHSP;
- 0.00920	      MyD*-	MySLD*+   K0		PHSP;
- 0.00150	      MyD-	MyD0	K+		PHSP;
- 0.00630	      MyD-	MyD*0	K+		PHSP;
- 0.00230	      MyD*-	MyD0	K+		PHSP;
- 0.01120	      MyD*-	MyD*0	K+		PHSP;
-
- 0.00250	      MyD-	MyD0	K*+		PHSP;
- 0.00250	      MyD*-	MyD0	K*+		PHSP;
- 0.00250	      MyD-	MyD*0	K*+		PHSP;
- 0.00500	      MyD*-	MyD*0	K*+		PHSP;
-
- 0.00250	      MyD-	MySLD+	K*0		PHSP;
- 0.00250	      MySLD-	MyD*+	K*0		PHSP;
- 0.00250	      MyD*-	MySLD+	K*0		PHSP;
- 0.00500	      MySLD*-	MyD*+	K*0		PHSP;
-
- 0.00180       	MyD_s+	MyD-	pi0		PHSP;
- 0.00370	      MyD_s+	MyAntiD0   pi-		PHSP;
- 0.00180	      MyD_s*+ 	MyD-	   pi0		PHSP;
- 0.00370   	      MyD_s*+	MyAntiD0   pi-		PHSP;
- 0.00300      	      MyD_s+	MyD-	   pi-   pi+	PHSP;
- 0.00220      	      MyD_s+	MyD-	   pi0	 pi0	PHSP;
- 0.00220      	      MyD_s+	MyAntiD0   pi-	 pi0	PHSP;
- 0.00300      	      MyD_s*+	MyD-	   pi-	 pi+	PHSP;
- 0.00220      	      MyD_s*+	MyD-	   pi0	 pi0	PHSP;
- 0.00220      	      MyD_s*+	MyAntiD0   pi-	 pi0	PHSP;
-
- 0.000211	      MyD-	MySLD+			PHSP;
- 0.000200	      MySLD*-	MyD+			SVS;
- 0.000600	      MyD*+	MySLD-			SVS;
- 0.000800	      MySLD*-	MyD*+			SVV_HELAMP 0.4904 0.0 0.7204 0.0 0.4904 0.0;
-	      
- 0.00075              MyD-      MyD_s0*+        PHSP;
- 0.00090              MyD*-     MyD_s0*+        SVS;
- 0.00310              MyD_s1+	MyD-		SVS;
- 0.01200              MyD*-     MyD_s1+         SVV_HELAMP 0.4904 0. 0.7204 0. 0.4904 0.0;
-	      
- 0.00045              MyD-      MyD'_s1+        PHSP;
- 0.00094              MyD*-     MyD'_s1+        PHSP;
-Enddecay
-CDecay anti-B0sig
-#
-SetLineshapePW MyD_1+ MyD*+ pi0 2
-SetLineshapePW MyD_1- MyD*- pi0 2
-SetLineshapePW MyD_10 MyD*+ pi- 2
-SetLineshapePW MyAntiD_10 MyD*- pi+ 2
-#
-SetLineshapePW MyD_2*+ MyD*+ pi0 2
-SetLineshapePW MyD_2*- MyD*- pi0 2
-SetLineshapePW MyD_2*0 MyD*+ pi- 2
-SetLineshapePW MyAntiD_2*0 MyD*- pi+ 2 
-#
-Decay MyD0
- 0.0333          K-      mu+     nu_mu   PHOTOS ISGW2;
- 0.0192         K*-     mu+     nu_mu  PHOTOS ISGW2;
- 0.000076       K_1-    mu+     nu_mu  PHOTOS ISGW2;
- 0.00110        K_2*-   mu+     nu_mu  PHOTOS ISGW2;
- 0.00238        pi-     mu+     nu_mu  PHOTOS PHSP;
- 0.00040        K-      pi0     mu+    nu_mu  PHOTOS PHSP;
-Enddecay
-CDecay MyAntiD0
-#
-Decay MyD*0
- 0.619          MyD0    pi0             VSS;
- 0.381          MyD0    gamma           VSP_PWAVE;
-Enddecay
-CDecay MyAntiD*0
-#
-Decay MyD_2*+
- 0.10300        MyD*+   pi0             TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
- 0.20900        MyD*0   pi+             TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
- 0.22900        MyD+    pi0             TSS;
- 0.45900        MyD0    pi+             TSS;
-Enddecay
-CDecay MyD_2*-
-#
-Decay MyD'_1+
- 0.33333        MyD*+   pi0             VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
- 0.66667        MyD*0   pi+             VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
-Enddecay
-CDecay MyD'_1+
-#
-Decay MyD*+
- 0.67700        MyD0    pi+             VSS;
- 0.30700        MyD+    pi0             VSS;
- 0.01600        MyD+    gamma           VSP_PWAVE;
-Enddecay
-CDecay MyD*-
-#
-Decay MySLD*+
- 0.67700        MyD0    pi+             VSS;
- 0.30700        MySLD+    pi0             VSS;
- 0.01600        MySLD+    gamma           VSP_PWAVE;
-Enddecay
-CDecay MySLD*-
-#
-Decay MyD_s*+
- 0.93500        MyD_s+  gamma           VSP_PWAVE;
- 0.05800        MyD_s+  pi0             VSS;
-Enddecay
-CDecay MyD_s*-
-#
-Decay MyD_s+
- 0.0259 phi      mu+     nu_mu                      PHOTOS  ISGW2;
- 0.0267 eta      mu+     nu_mu                      PHOTOS  ISGW2;
- 0.0099 eta'     mu+     nu_mu                      PHOTOS  ISGW2;
- 0.0037 K0       mu+     nu_mu                      PHOTOS  ISGW2;  
- 0.0018 K*0      mu+     nu_mu                      PHOTOS  ISGW2;
- 0.0020 f_0      mu+     nu_mu                      PHOTOS  ISGW2; 
- 0.0059 mu+      nu_mu                              PHOTOS   SLN; 
-Enddecay
-CDecay MyD_s-
-#
-Decay MySLD+
-  0.05280       anti-K*0        mu+     nu_mu   PHOTOS ISGW2;
-  0.09200       anti-K0         mu+     nu_mu   PHOTOS ISGW2;
-  0.00277       anti-K_10       mu+     nu_mu   PHOTOS ISGW2;
-  0.00293       anti-K_2*0      mu+     nu_mu   PHOTOS ISGW2;
-  0.00292       K-   pi+        mu+     nu_mu   PHOTOS PHSP;
-  0.00120       anti-K0   pi0   mu+     nu_mu   PHOTOS PHSP;
-Enddecay
-CDecay MySLD-
-#
-Decay MyD+
-  0.09460       K-        pi+   pi+             PHOTOS D_DALITZ;
-Enddecay
-CDecay MyD-
-#
-Decay MyD_s1+
- 0.48000        MyD_s*+  pi0             PARTWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
- 0.18000        MyD_s+   gamma           VSP_PWAVE;
- 0.04300        MyD_s+   pi+    pi-      PHSP;
- 0.0037         MyD_s0*+ gamma           PHSP;
-Enddecay
-CDecay MyD_s1-
-#
-Decay MyD'_s1+
- 0.42500        MyD*+     K0            VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
- 0.42500        MyD*0     K+            VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
- 0.02800        MyD+      pi-    K+     PHSP;
- 0.03400        MyD+      K0            PHSP;
-Enddecay
-CDecay MyD'_s1-
-#
-Decay MyD_1+
-0.3333    MyD*+		pi0              VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
-0.6667    MyD*0 	pi+              VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
-Enddecay
-CDecay MyD_1-
-#
-Decay MyD_s0*+
- 1.00000        MyD_s+  pi0     PHSP;
-Enddecay
-CDecay MyD_s0*-
-#
-End
diff --git a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCut.dec b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCut.dec
deleted file mode 100644
index 10cc752f0..000000000
--- a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCut.dec
+++ /dev/null
@@ -1,309 +0,0 @@
-# EventType: 12995604
-#
-# Descriptor: [[B+] ==> (D*(2007)~0 -> (D~0 -> K+ mu- nu_mu~) pi0) (D_s*+ -> (D_s+ -> K+ K- pi+) gamma)]cc
-#
-# NickName: Bu_DD,DD=cocktail,D+muRDplusCut
-# Cuts: LoKi::GenCutTool/TightCut
-#
-# InsertPythonCode:
-#
-# from Configurables import LoKi__GenCutTool
-# from Gauss.Configuration import *
-# gen = Generation()
-# gen.SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' )
-# tightCut = gen.SignalRepeatedHadronization.TightCut
-# tightCut.Decay = "[B+ ==> ^((Charm) -> ^mu+ nu_mu ... ) ^(D- -> ^K+ ^pi- ^pi- ) {X} {X} {X} {X} {X} {X} {X} ]CC"
-# tightCut.Preambulo += [
-# "from GaudiKernel.SystemOfUnits import mrad" ,
-# "FilterD = GNINTREE (GCHARM, HepMC.parents)",
-# "FromD   = 1 == FilterD",
-# "BCut = (GTHETA < 400.0*mrad)"
-#
-#]
-#tightCut.Cuts    =    {
-# '[mu+]cc'     : "FromD",
-# '[B+]cc'	: "BCut"
-#}
-#
-# EndInsertPythonCode
-#
-# Documentation: Sum of B+ decaying to two D mesons, one D decaying semileptonically. Filtered for a final state with D+mu.
-# 		 Background for B2DmunuX analyses.
-#
-# EndDocumentation
-#
-# CPUTime: <1min
-# PhysicsWG: B2SL
-# Tested: Yes
-# Responsible: Patrick Owen
-# Email: patrick.haworth.owen@gmail.com
-# Date:   20160803
-#
-Alias        MyD0		D0
-Alias        MyAntiD0    	anti-D0
-ChargeConj   MyD0	 	MyAntiD0
-#
-Alias        MyD+		D+
-Alias	     MyD-		D-		
-ChargeConj   MyD+		MyD-
-#
-Alias        MySLD+		D+
-Alias	     MySLD-		D-		
-ChargeConj   MySLD+		MySLD-
-#
-Alias        MyD_s+		D_s+
-Alias	     MyD_s-		D_s-		
-ChargeConj   MyD_s+		MyD_s-		
-#
-Alias	     MyD_1+		D_1+
-Alias	     MyD_1-		D_1-
-ChargeConj   MyD_1+		MyD_1-
-#
-Alias	     MyD_s1+		D_s1+
-Alias	     MyD_s1-		D_s1-		
-ChargeConj   MyD_s1+		MyD_s1-
-#
-Alias        MyD_10		D_10
-Alias	     MyAntiD_10		anti-D_10
-ChargeConj   MyD_10		MyAntiD_10
-#
-Alias	     MyD(2S)0		D(2S)0
-Alias	     MyAntiD(2S)0	anti-D(2S)0	
-ChargeConj   MyD(2S)0		MyAntiD(2S)0
-#
-Alias	     MyD(2S)+		D(2S)+
-Alias	     MyD(2S)-		D(2S)-		
-ChargeConj   MyD(2S)+		MyD(2S)-
-#
-Alias	     MyD*0		D*0
-Alias	     MyAntiD*0		anti-D*0
-ChargeConj   MyD*0	 	MyAntiD*0
-#
-Alias	     MyD_2*0		D_2*0
-Alias	     MyAntiD_2*0 	anti-D_2*0	
-ChargeConj   MyD_2*0		MyAntiD_2*0
-#
-Alias	     MyD_0*0		D_0*0
-Alias 	     MyAntiD_0*0	anti-D_0*0		
-ChargeConj   MyD_0*0		MyAntiD_0*0
-#
-Alias	     MyD*(2S)0		D*(2S)0
-Alias	     MyAntiD*(2S)0	anti-D*(2S)0	
-ChargeConj   MyD*(2S)0		MyAntiD*(2S)0
-#
-Alias	     MyD*+		D*+
-Alias	     MyD*-		D*-		
-ChargeConj   MyD*+		MyD*-
-#
-Alias	     MySLD*+		D*+
-Alias	     MySLD*-		D*-		
-ChargeConj   MySLD*+		MySLD*-
-#
-Alias	     MyD_2*+		D_2*+
-Alias	     MyD_2*-		D_2*-		
-ChargeConj   MyD_2*+		MyD_2*-
-#
-Alias        MyD_s*+		D_s*+
-Alias	     MyD_s*-		D_s*-		
-ChargeConj   MyD_s*+		MyD_s*-
-#
-Alias        MyD_s2*+		D_s2*+
-Alias	     MyD_s2*-		D_s2*-		
-ChargeConj   MyD_s2*+		MyD_s2*-
-#
-Alias	     MyD_0*+		D_0*+
-Alias	     MyD_0*-		D_0*-		
-ChargeConj   MyD_0*+		MyD_0*-
-#
-Alias	     MyD_s0*+		D_s0*+
-Alias	     MyD_s0*-		D_s0*-
-ChargeConj   MyD_s0*+		MyD_s0*-
-#
-Alias	     MyD*(2S)+		D*(2S)+
-Alias	     MyD*(2S)-		D*(2S)-		
-ChargeConj   MyD*(2S)+		MyD*(2S)-
-#
-Alias	     MyD'_10		D'_10
-Alias	     MyAntiD'_10	anti-D'_10
-ChargeConj   MyD'_10		MyAntiD'_10
-#
-Alias	     MyD'_s1+		D'_s1+
-Alias	     MyD'_s1-		D'_s1-		
-ChargeConj   MyD'_s1+		MyD'_s1-
-#
-Alias	     MyD'_1+		D'_1+
-Alias	     MyD'_1-		D'_1-		
-ChargeConj   MyD'_1+		MyD'_1-
-#
-Decay B+sig
- 0.00900      MyAntiD0		MyD_s+			PHSP;
- 0.00760      MyD_s*+		MyAntiD0		SVS;
- 0.00820      MyAntiD*0		MyD_s+			SVS;
- 0.01710      MyAntiD*0		MyD_s*+			SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
- 0.00120      MyAntiD'_10	MyD_s+			SVS;
- 0.00240      MyAntiD'_10	MyD_s*+			SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
- 0.00240      MyAntiD_10	MyD_s+			SVS;
- 0.00480      MyAntiD_10	MyD_s*+			SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
- 0.00840      MyAntiD_2*0	MyD_s+			STS;
- 0.00800      MyAntiD_2*0	MyD_s*+			PHSP;
-
- 0.00160      MyAntiD0		MyD+	K0	PHSP;
- 0.00380      MyAntiD0		MyD*+	K0	PHSP;
- 0.00210      MyAntiD*0		MyD+	K0	PHSP;
- 0.00920      MyAntiD*0		MyD*+	K0	PHSP;
-
- 0.00250      MyAntiD0		MyD+	K*0		PHSP;
- 0.00250      MyAntiD0		MyD*+	K*0		PHSP;
- 0.00250      MyAntiD*0		MyD+	K*0		PHSP;
- 0.00500      MyAntiD*0		MyD*+	K*0		PHSP;
-
- 0.00360      MyD_s+	MyD-		pi+		PHSP;
- 0.00370      MyD_s*+ 	MyD- 	    	pi+		PHSP;
-
- 0.00330      MyD_s+	MyD-	   pi0   pi+	PHSP;
- 0.00330      MyD_s*+	MyD-	   pi0	 pi+	PHSP;
-
- 0.000400      MyAntiD0		MyD+			PHSP;
- 0.000600      MyAntiD*0	MyD+			SVS;
- 0.000400      MyD*+		MyAntiD0		SVS;
- 0.000800      MyAntiD*0	MyD*+			SVV_HELAMP 0.283 0.0 0.932 0.0 0.228 0.0;
-
-
- 0.00020	MySLD-	MyD+	K+	PHSP;
- 0.00060	MyD-	MySLD*+	K+	PHSP;
- 0.00060	MySLD*-	MyD+	K+	PHSP;
- 0.00130	MyD*-	MySLD*+	K+	PHSP;
-
- 0.00050	MySLD-	MyD+	K*+	PHSP;
- 0.00050	MyD-	MySLD*+	K*+	PHSP;
- 0.00050	MySLD*-	MyD+	K*+	PHSP;
- 0.00100	MyD*-	MySLD*+	K*+	PHSP;      
-Enddecay
-CDecay B-sig
-#
-SetLineshapePW MyD_1+ MyD*+ pi0 2
-SetLineshapePW MyD_1- MyD*- pi0 2
-SetLineshapePW MyD_10 MyD*+ pi- 2
-SetLineshapePW MyAntiD_10 MyD*- pi+ 2
-#
-SetLineshapePW MyD_2*+ MyD*+ pi0 2
-SetLineshapePW MyD_2*- MyD*- pi0 2
-SetLineshapePW MyD_2*0 MyD*+ pi- 2
-SetLineshapePW MyAntiD_2*0 MyD*- pi+ 2
-#
-Decay MyD0
- 0.0333          K-      mu+     nu_mu		PHOTOS ISGW2;
- 0.0192         K*-     mu+     nu_mu  		PHOTOS ISGW2;
- 0.000076       K_1-    mu+     nu_mu  		PHOTOS ISGW2;
- 0.00110        K_2*-   mu+     nu_mu  		PHOTOS ISGW2;
- 0.00238        pi-     mu+     nu_mu  		PHOTOS PHSP;
- 0.00040        K-      pi0     mu+    nu_mu  	PHOTOS PHSP;
-Enddecay
-CDecay MyAntiD0
-#
-Decay MyD*0
- 0.619          MyD0    pi0             VSS;
- 0.381          MyD0    gamma           VSP_PWAVE;
-Enddecay
-CDecay MyAntiD*0
-#
-Decay MyD_2*+
- 0.10300        MyD*+   pi0             TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
- 0.20900        MyD*0   pi+             TVS_PWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
- 0.22900        MyD+    pi0             TSS;
- 0.45900        MyD0    pi+             TSS;
-Enddecay
-CDecay MyD_2*-
-#
-Decay MyD'_1+
- 0.33333        MyD*+   pi0             VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
- 0.66667        MyD*0   pi+             VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
-Enddecay
-CDecay MyD'_1+
-#
-Decay MyD*+
- 0.67700        MyD0    pi+             VSS;
- 0.30700        MyD+    pi0             VSS;
- 0.01600        MyD+    gamma           VSP_PWAVE;
-Enddecay
-CDecay MyD*-
-#
-Decay MySLD*+
- 0.67700        MyD0    pi+             VSS;
- 0.30700        MySLD+    pi0             VSS;
- 0.01600        MySLD+    gamma           VSP_PWAVE;
-Enddecay
-CDecay MySLD*-
-#
-Decay MyD_s*+
- 0.93500        MyD_s+  gamma           VSP_PWAVE;
- 0.05800        MyD_s+  pi0             VSS;
-Enddecay
-CDecay MyD_s*-
-#
-Decay MyD_s+
- 0.0259 phi      mu+     nu_mu                      PHOTOS  ISGW2;
- 0.0267 eta      mu+     nu_mu                      PHOTOS  ISGW2;
- 0.0099 eta'     mu+     nu_mu                      PHOTOS  ISGW2;
- 0.0037 K0       mu+     nu_mu                      PHOTOS  ISGW2;  
- 0.0018 K*0      mu+     nu_mu                      PHOTOS  ISGW2;
- 0.0020 f_0      mu+     nu_mu                      PHOTOS  ISGW2; 
- 0.0059 mu+      nu_mu                              PHOTOS   SLN; 
-Enddecay
-CDecay MyD_s-
-#
-Decay MySLD+
-  0.05280       anti-K*0        mu+     nu_mu   PHOTOS ISGW2;
-  0.09200       anti-K0         mu+     nu_mu   PHOTOS ISGW2;
-  0.00277       anti-K_10       mu+     nu_mu   PHOTOS ISGW2;
-  0.00293       anti-K_2*0      mu+     nu_mu   PHOTOS ISGW2;
-  0.00292       K-   pi+        mu+     nu_mu   PHOTOS PHSP;
-  0.00120       anti-K0   pi0   mu+     nu_mu   PHOTOS PHSP;
-Enddecay
-CDecay MySLD-
-#
-Decay MyD+
-  0.09460       K-        pi+   pi+             PHOTOS D_DALITZ;
-Enddecay
-CDecay MyD-
-#
-Decay MyD_s1+
- 0.48000        MyD_s*+  pi0             PARTWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
- 0.18000        MyD_s+   gamma           VSP_PWAVE;
- 0.04300        MyD_s+   pi+    pi-      PHSP;
- 0.0037         MyD_s0*+ gamma           PHSP;
-Enddecay
-CDecay MyD_s1-
-#
-Decay MyD'_s1+
- 0.42500        MyD*+     K0            VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
- 0.42500        MyD*0     K+            VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
- 0.02800        MyD+      pi-    K+     PHSP;
- 0.03400        MyD+      K0            PHSP;
-Enddecay
-CDecay MyD'_s1-
-#
-Decay MyD_1+
-0.3333    MyD*+ pi0                        VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
-0.6667    MyD*0 pi+                        VVS_PWAVE  0.0 0.0 0.0 0.0 1.0 0.0;
-Enddecay
-CDecay MyD_1-
-#
-Decay MyD_s0*+
- 1.00000        MyD_s+  pi0     PHSP;
-Enddecay
-CDecay MyD_s0*-
-#
-Decay MyD'_10
- 0.66667	MyD*+	pi-	PHOTOS VVS_PWAVE 1.0 0.0 0.0 0.0 0.0 0.0;
-Enddecay
-CDecay MyAntiD'_10
-#
-Decay MyD_2*0
-  0.40157      MyD+    pi-     PHOTOS TSS;
-  0.26509      MyD*+   pi-     PHOTOS TVS_PWAVE  0.0 0.0 1.0 0.0 0.0 0.0;
-Enddecay
-CDecay MyAntiD_2*0
-#
-End
-
diff --git a/doc/table_obsolete.sql b/doc/table_obsolete.sql
index 63d48e338..630e64a5d 100644
--- a/doc/table_obsolete.sql
+++ b/doc/table_obsolete.sql
@@ -405,3 +405,5 @@ EVTTYPEID = 39122249, DESCRIPTION = "pi0_gammaA,ee,displaced,mA=20MeV,badlifetim
 EVTTYPEID = 39122940, DESCRIPTION = "pi0_gammaA,ee,displaced,mA=50MeV,badlifetime", PRIMARY = "[eta -> gamma (A' -> e+ e-)]"
 EVTTYPEID = 39122246, DESCRIPTION = "pi0_gammaA,ee,displaced,mA=5MeV,badlifetime", PRIMARY = "[eta -> gamma (A' -> e+ e-)]"
 EVTTYPEID = 12513011, DESCRIPTION = "Bu_phimunu=TightCut,BToVlnuBall,buggy", PRIMARY = "[B+ -> (phi(1020) -> K+ K-) mu+ nu_mu]cc"
+EVTTYPEID = 11995204, DESCRIPTION = "Bd_DD,DD=cocktail,D+muRDplusCut", PRIMARY = "[B0 -> (D- -> K+ pi- pi) D_s+]cc"
+EVTTYPEID = 12995604, DESCRIPTION = "Bu_DD,DD=cocktail,D+muRDplusCut", PRIMARY = "[B+ -> (D- -> K+ pi- pi) D_s+ K+]cc"
-- 
GitLab


From 456d943e96d184e5c841d60d0fbd85a62262a078 Mon Sep 17 00:00:00 2001
From: Patrick Haworth Owen <patrick.haworth.owen@gmail.com>
Date: Wed, 19 Jul 2023 13:55:37 +0200
Subject: [PATCH 7/8] fix decaynumbers based on review

---
 dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec | 2 +-
 dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
index 7ff6aa766..1c7f5cccf 100644
--- a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
+++ b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -1,4 +1,4 @@
-# EventType: 11995205
+# EventType: 11694605
 #
 # Descriptor: [[B0] ==> (D*- -> (D~0 -> K+ mu- nu_mu~)) (D_s*+ -> gamma (D_s+ -> K+ K- pi+))]cc
 #
diff --git a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
index 3aa6489f1..be39df0b9 100644
--- a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
+++ b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -1,4 +1,4 @@
-# EventType: 12995605
+# EventType: 12695405
 #
 # Descriptor: [[B+] ==> (D*(2007)~0 -> (D~0 -> K+ mu- nu_mu~) pi0) (D_s*+ -> (D_s+ -> K+ K- pi+) gamma)]cc
 #
-- 
GitLab


From 05c598878aa8ca1e3e605763a97b1745c9e6e74b Mon Sep 17 00:00:00 2001
From: Patrick Haworth Owen <patrick.haworth.owen@gmail.com>
Date: Thu, 3 Aug 2023 10:37:57 +0200
Subject: [PATCH 8/8] fixes for sim10 compatibility

---
 dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec | 4 ++--
 dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec | 2 +-
 2 files changed, 3 insertions(+), 3 deletions(-)

diff --git a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
index 1c7f5cccf..2e6dc3d46 100644
--- a/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
+++ b/dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -12,7 +12,7 @@
 # gen = Generation()
 # gen.SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' )
 # tightCut = gen.SignalRepeatedHadronization.TightCut
-# tightCut.Decay = "[(Beauty) -> ((Charm) -> ^mu+ nu_mu ... ) ^(D- ==> ^K+ ^pi- ^pi-) {X} {X} {X} {X} {X} {X} {X} ]CC"
+# tightCut.Decay = "[(Beauty) ==> ((Charm) ==> ^mu+ nu_mu {X} {X} {X} {X} ) ^(D- ==> ^K+ ^pi- ^pi-) {X} {X} {X} {X} {X} {X} {X} ]CC"
 #tightCut.Preambulo += [
 #  "from LoKiCore.functions import in_range"  ,
 #  "from GaudiKernel.SystemOfUnits import GeV, MeV"  ,
@@ -32,7 +32,7 @@
 #                
 # EndDocumentation
 #
-# CPUTime: <2min
+# CPUTime: <1min
 # PhysicsWG: B2SL
 # Tested: Yes
 # Responsible: Patrick Owen
diff --git a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
index be39df0b9..4755ce7fc 100644
--- a/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
+++ b/dkfiles/Bu_DD,DD=cocktail,D+muRDplusCutFixed.dec
@@ -12,7 +12,7 @@
 # gen = Generation()
 # gen.SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' )
 # tightCut = gen.SignalRepeatedHadronization.TightCut
-# tightCut.Decay = "[(Beauty) ==> ((Charm) -> ^mu+ nu_mu ... ) ^(D- ==> ^K+ ^pi- ^pi- ) {X} {X} {X} {X} {X} {X} {X} ]CC"
+# tightCut.Decay = "[(Beauty) ==> ((Charm) ==> ^mu+ nu_mu {X} {X} {X} {X} ) ^(D- ==> ^K+ ^pi- ^pi- ) {X} {X} {X} {X} {X} {X} {X} ]CC"
 #tightCut.Preambulo += [
 #  "from LoKiCore.functions import in_range"  ,
 #  "from GaudiKernel.SystemOfUnits import GeV, MeV"  ,
-- 
GitLab