Skip to content
Snippets Groups Projects

Draft: Add B2Dh, D2Kspipi with BELLE2018 model

Open Camille Normand requested to merge canorman-B2Dh-BELLE2018 into master
1 unresolved thread
Files
2
+ 66
0
# 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
#
Loading