Skip to content
Snippets Groups Projects

Draft: Port New B-->DD decfiles for the RD+ analysis to Sim10

Open Draft: Port New B-->DD decfiles for the RD+ analysis to Sim10
DDfix into Sim10
Open Michal Kreps requested to merge DDfix into Sim10
Compare
and
3 files
+ 59
− 66
Compare changes
  • Side-by-side
  • Inline
Files
3
dkfiles/Bd_DD,DD=cocktail,D+muRDplusCut.dec → dkfiles/Bd_DD,DD=cocktail,D+muRDplusCutFixed.dec
+ 27
− 31
# EventType: 11995204
# EventType: 11694605
#
# 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
# NickName: Bd_DD,DD=cocktail,D+muRDplusCutFixed
# Cuts: LoKi::GenCutTool/TightCut
#
# InsertPythonCode:
@@ -12,17 +12,18 @@
# gen = Generation(genAlgName)
# 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.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" ,
# "pipiKP = GCHILD(GP,1) + GCHILD(GP,2) + GCHILD(GP,3)" ,
# "pipiKPT = GCHILD(GPT,1) + GCHILD(GPT,2) + GCHILD(GPT,3)" ,
# "CS = LoKi.GenChild.Selector",
# 'pipiKP = GCHILD(GP, CS("[(Beauty) ==> ((Charm) ==> mu+ nu_mu {X} {X} {X} {X} ) (D- ==> ^K+ pi- pi-) {X} {X} {X} {X} {X} {X} {X} ]CC")) + GCHILD(GP, CS("[(Beauty) ==> ((Charm) ==> mu+ nu_mu {X} {X} {X} {X} ) (D- ==> K+ ^pi- pi-) {X} {X} {X} {X} {X} {X} {X} ]CC")) + GCHILD(GP, CS("[(Beauty) ==> ((Charm) ==> mu+ nu_mu {X} {X} {X} {X} ) (D- ==> K+ pi- ^pi-) {X} {X} {X} {X} {X} {X} {X} ]CC"))',
# 'pipiKPT = GCHILD(GPT, CS("[(Beauty) ==> ((Charm) ==> mu+ nu_mu {X} {X} {X} {X} ) (D- ==> ^K+ pi- pi-) {X} {X} {X} {X} {X} {X} {X} ]CC")) + GCHILD(GPT, CS("[(Beauty) ==> ((Charm) ==> mu+ nu_mu {X} {X} {X} {X} ) (D- ==> K+ ^pi- pi-) {X} {X} {X} {X} {X} {X} {X} ]CC")) + GCHILD(GPT, CS("[(Beauty) ==> ((Charm) ==> mu+ nu_mu {X} {X} {X} {X} ) (D- ==> K+ pi- ^pi-) {X} {X} {X} {X} {X} {X} {X} ]CC"))'
#]
#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) ",
# '[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
@@ -37,7 +38,7 @@
# Tested: Yes
# Responsible: Patrick Owen
# Email: patrick.haworth.owen@gmail.com
# Date: 20160802
# Date: 20230712
#
Alias MyD0 D0
Alias MyAntiD0 anti-D0
@@ -138,8 +139,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,48 +150,44 @@ 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.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.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
#
@@ -235,7 +232,6 @@ Enddecay
CDecay MyD'_1+
#
Decay MyD*+
0.67700 MyD0 pi+ VSS;
0.30700 MyD+ pi0 VSS;
0.01600 MyD+ gamma VSP_PWAVE;
Enddecay