Skip to content
Snippets Groups Projects
Commit 3ede0aa8 authored by Hanae Tilquin's avatar Hanae Tilquin Committed by Rosen Matev
Browse files

Added XTauTau electronic lines

parent 372847db
No related branches found
No related tags found
1 merge request!1451Added XTauTau electronic lines
Expand all Hide whitespace changes
Inline Side-by-side
Showing
with 6786 additions and 40 deletions
Hlt/Hlt2Conf/python/Hlt2Conf/lines/rd/__init__.py
+ 10
− 0
View file @ 3ede0aa8
......@@ -30,6 +30,9 @@ from . import qqbar_to_ll
from . import b_to_ll_LFV
from . import b_to_kstarmumu
from . import hlt2_btosmutau_exclusive
from . import btosetau_exclusive_hlt2
from . import btosmutau_tau_to_e_exclusive_hlt2
from . import btosetau_tau_to_e_exclusive_hlt2
from . import hlt2_btostautau_exclusive
from . import b_to_xgamma_exclusive_hlt2
from . import b_to_xgamma_exclusive_spruce
......@@ -43,6 +46,8 @@ from . import baryonic_radiative
from . import bnv_lines_hlt2
from . import lfv_lines_hlt2
from . import b_to_tautau_hlt2
from . import btostautau_mue_exclusive_hlt2
from . import btostautau_ee_exclusive_hlt2
# provide "all_lines" for correct registration by the overall HLT2 lines module
all_lines = {}
......@@ -61,6 +66,9 @@ all_lines.update(qqbar_to_ll.all_lines)
all_lines.update(b_to_ll_LFV.hlt2_lines)
all_lines.update(b_to_kstarmumu.all_lines)
all_lines.update(hlt2_btosmutau_exclusive.all_lines)
all_lines.update(btosetau_exclusive_hlt2.all_lines)
all_lines.update(btosmutau_tau_to_e_exclusive_hlt2.all_lines)
all_lines.update(btosetau_tau_to_e_exclusive_hlt2.all_lines)
all_lines.update(hlt2_btostautau_exclusive.all_lines)
all_lines.update(b_to_xgamma_exclusive_hlt2.all_lines)
all_lines.update(b_to_xll_hlt2.all_lines)
......@@ -72,6 +80,8 @@ all_lines.update(baryonic_radiative.all_lines)
all_lines.update(bnv_lines_hlt2.all_lines)
all_lines.update(lfv_lines_hlt2.all_lines)
all_lines.update(b_to_tautau_hlt2.hlt2_lines)
all_lines.update(btostautau_mue_exclusive_hlt2.all_lines)
all_lines.update(btostautau_ee_exclusive_hlt2.all_lines)
sprucing_lines = {}
sprucing_lines.update(b_to_ll_spruce.sprucing_lines)
......
Hlt/Hlt2Conf/python/Hlt2Conf/lines/rd/btosetau_exclusive_hlt2.py 0 → 100644
+ 1378
− 0
View file @ 3ede0aa8
This diff is collapsed.
Hlt/Hlt2Conf/python/Hlt2Conf/lines/rd/btosetau_tau_to_e_exclusive_hlt2.py 0 → 100644
+ 1388
− 0
View file @ 3ede0aa8
This diff is collapsed.
Hlt/Hlt2Conf/python/Hlt2Conf/lines/rd/btosmutau_tau_to_e_exclusive_hlt2.py 0 → 100644
+ 1398
− 0
View file @ 3ede0aa8
This diff is collapsed.
Hlt/Hlt2Conf/python/Hlt2Conf/lines/rd/btostautau_ee_exclusive_hlt2.py 0 → 100644
+ 895
− 0
View file @ 3ede0aa8
This diff is collapsed.
Hlt/Hlt2Conf/python/Hlt2Conf/lines/rd/btostautau_mue_exclusive_hlt2.py 0 → 100644
+ 1412
− 0
View file @ 3ede0aa8
This diff is collapsed.
Hlt/Hlt2Conf/python/Hlt2Conf/lines/rd/builders/btosetau_exclusive.py 0 → 100644
+ 220
− 0
View file @ 3ede0aa8
###############################################################################
# (c) Copyright 2020-2021 CERN for the benefit of the LHCb Collaboration #
# #
# This software is distributed under the terms of the GNU General Public #
# Licence version 3 (GPL Version 3), copied verbatim in the file "COPYING". #
# #
# In applying this licence, CERN does not waive the privileges and immunities #
# granted to it by virtue of its status as an Intergovernmental Organization #
# or submit itself to any jurisdiction. #
###############################################################################
from GaudiKernel.SystemOfUnits import MeV, mm
import Functors as F
from Hlt2Conf.algorithms_thor import ParticleCombiner
def make_phie(kaons1,
kaons2,
electrons,
pvs,
comb_m_max=3600 * MeV,
vchi2pdof_max=8,
docachi2_max=8.0,
doca_max=0.8 * mm,
dira_min=0.5,
phi_docachi2_max=4.0,
phi_doca_max=0.15 * mm,
decay_descriptor="[B0 -> K+ K- e+]cc",
name="rd_dikaon_electron_for_btosetau_{hash}"):
"""Builder for X -> K K e decays, used in XETau lines"""
combination_code = F.require_all(F.MASS < comb_m_max,
F.MAXSDOCACHI2CUT(docachi2_max),
F.MAXSDOCACUT(doca_max))
two_body_combination_code = F.require_all(
F.MASS < comb_m_max, F.MAXSDOCACHI2CUT(phi_docachi2_max),
F.MAXSDOCACUT(phi_doca_max))
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min)
phie = ParticleCombiner([kaons1, kaons2, electrons],
DecayDescriptor=decay_descriptor,
Combination12Cut=two_body_combination_code,
CombinationCut=combination_code,
CompositeCut=vertex_code,
name=name)
return phie
def make_ke(kaons,
electrons,
pvs,
comb_m_max=3550 * MeV,
vchi2pdof_max=4,
docachi2_max=4.0,
doca_max=0.6 * mm,
dira_min=0.5,
decay_descriptor="[B0 -> e+ K-]cc",
name="rd_kaon_electron_for_btosetau_{hash}"):
"""Builder for X -> K e decays, used in XETau lines"""
combination_code = F.require_all(F.MASS < comb_m_max,
F.MAXSDOCACHI2CUT(docachi2_max),
F.MAXSDOCACUT(doca_max))
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min)
ke = ParticleCombiner([electrons, kaons],
DecayDescriptor=decay_descriptor,
CombinationCut=combination_code,
CompositeCut=vertex_code,
name=name)
return ke
def make_kste(kaons,
pions,
electrons,
pvs,
comb_m_max=3550 * MeV,
vchi2pdof_max=6,
docachi2_max=6.0,
doca_max=0.6 * mm,
dira_min=0.5,
kst_docachi2_max=3.0,
kst_doca_max=0.12 * mm,
decay_descriptor="[B0 -> K+ pi- e+]cc",
name="rd_kpi_electron_for_btosetau_{hash}"):
"""Builder for X -> K pi e decays, used in XETau lines"""
combination_code = F.require_all(F.MASS < comb_m_max,
F.MAXSDOCACHI2CUT(docachi2_max),
F.MAXSDOCACUT(doca_max))
two_body_combination_code = F.require_all(
F.MASS < comb_m_max, F.MAXSDOCACHI2CUT(kst_docachi2_max),
F.MAXSDOCACUT(kst_doca_max))
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min)
kste = ParticleCombiner([kaons, pions, electrons],
DecayDescriptor=decay_descriptor,
Combination12Cut=two_body_combination_code,
CombinationCut=combination_code,
CompositeCut=vertex_code,
name=name)
return kste
def make_pke(kaons,
protons,
electrons,
pvs,
comb_m_max=4000 * MeV,
vchi2pdof_max=8,
docachi2_max=8.0,
doca_max=0.8 * mm,
dira_min=0.5,
pk_docachi2_max=4.0,
pk_doca_max=0.15 * mm,
decay_descriptor="[B0 -> p+ K- e+]cc",
name="rd_pk_electron_for_btosetau_{hash}"):
"""Builder for X -> p K e decays, used in XETau lines"""
combination_code = F.require_all(F.MASS < comb_m_max,
F.MAXSDOCACHI2CUT(docachi2_max),
F.MAXSDOCACUT(doca_max))
two_body_combination_code = F.require_all(
F.MASS < comb_m_max, F.MAXSDOCACHI2CUT(pk_docachi2_max),
F.MAXSDOCACUT(pk_doca_max))
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min)
pke = ParticleCombiner([protons, kaons, electrons],
DecayDescriptor=decay_descriptor,
Combination12Cut=two_body_combination_code,
CombinationCut=combination_code,
CompositeCut=vertex_code,
name=name)
return pke
def make_bs(kkes,
electrons,
pvs,
comb_m_max=7000 * MeV,
vchi2pdof_max=75,
dira_min=0.999,
bpvfdchi2_min=50,
decay_descriptor="[B_s0 -> B0 e-]cc",
name="rd_make_bs_to_kktaue_{hash}"):
"""Builder for Bs -> X e tau(-> e nu nu) decays, where X is K K"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min,
F.BPVFDCHI2(pvs) > bpvfdchi2_min)
bs = ParticleCombiner([kkes, electrons],
DecayDescriptor=decay_descriptor,
CombinationCut=combination_code,
CompositeCut=vertex_code,
name=name)
return bs
def make_bu(kes,
electrons,
pvs,
comb_m_max=6750 * MeV,
vchi2pdof_max=25,
dira_min=0.999,
bpvfdchi2_min=50,
decay_descriptor="[B_s0 -> B0 e-]cc",
name="rd_make_bu_to_ketau_{hash}"):
"""Builder for Bs -> X e tau(-> e nu nu) decays, where X is K"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min,
F.BPVFDCHI2(pvs) > bpvfdchi2_min)
bu = ParticleCombiner([kes, electrons],
DecayDescriptor=decay_descriptor,
CombinationCut=combination_code,
CompositeCut=vertex_code,
name=name)
return bu
def make_lb(pkes,
electrons,
pvs,
comb_m_max=7250 * MeV,
vchi2pdof_max=50,
dira_min=0.999,
bpvfdchi2_min=50,
decay_descriptor="[Lambda_b0 -> B0 e-]cc",
name="rd_make_lb_to_pketau_{hash}"):
"""Builder for Lb -> X e tau(-> e nu nu) decays, where X is p K"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min,
F.BPVFDCHI2(pvs) > bpvfdchi2_min)
lb = ParticleCombiner([pkes, electrons],
DecayDescriptor=decay_descriptor,
CombinationCut=combination_code,
CompositeCut=vertex_code,
name=name)
return lb
def make_bd(kstes,
electrons,
pvs,
comb_m_max=6750 * MeV,
vchi2pdof_max=50,
dira_min=0.999,
bpvfdchi2_min=50,
decay_descriptor="[B_s0 -> B0 e+]cc",
name="rd_make_bd_to_kpietau_{hash}"):
"""Builder for Bd -> X e tau(-> e nu nu) decays, where X is K pi"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min,
F.BPVFDCHI2(pvs) > bpvfdchi2_min)
bd = ParticleCombiner([kstes, electrons],
DecayDescriptor=decay_descriptor,
CombinationCut=combination_code,
CompositeCut=vertex_code,
name=name)
return bd
Hlt/Hlt2Conf/python/Hlt2Conf/lines/rd/builders/btosmutau_exclusive.py
+ 17
− 15
View file @ 3ede0aa8
......@@ -13,6 +13,7 @@ from GaudiKernel.SystemOfUnits import MeV, mm
import Functors as F
from Hlt2Conf.algorithms_thor import ParticleCombiner
from PyConf import configurable
def make_phimu(kaons1,
......@@ -26,7 +27,7 @@ def make_phimu(kaons1,
phi_docachi2_max=4.0,
phi_doca_max=0.15 * mm,
decay_descriptor="[B0 -> K+ K- mu+]cc",
name='rd_dikaon_muon_for_btosmutau_{hash}'):
name="rd_dikaon_muon_for_btosmutau_{hash}"):
"""Builder for X -> K K mu decays, used in XMuTau lines"""
combination_code = F.require_all(F.MASS < comb_m_max,
F.MAXSDOCACHI2CUT(docachi2_max),
......@@ -52,7 +53,7 @@ def make_kmu(kaons,
docachi2_max=4.0,
doca_max=0.6 * mm,
decay_descriptor="[B0 -> mu+ K-]cc",
name='rd_kaon_muon_for_btosmutau_{hash}'):
name="rd_kaon_muon_for_btosmutau_{hash}"):
"""Builder for X -> K mu decays, used in XMuTau lines"""
combination_code = F.require_all(F.MASS < comb_m_max,
F.MAXSDOCACHI2CUT(docachi2_max),
......@@ -77,7 +78,7 @@ def make_kstmu(kaons,
kst_docachi2_max=3.0,
kst_doca_max=0.12 * mm,
decay_descriptor="[B0 -> K+ pi- mu+]cc",
name='rd_kpi_muon_for_btosmutau_{hash}'):
name="rd_kpi_muon_for_btosmutau_{hash}"):
"""Builder for X -> K pi mu decays, used in XMuTau lines"""
combination_code = F.require_all(F.MASS < comb_m_max,
F.MAXSDOCACHI2CUT(docachi2_max),
......@@ -106,7 +107,7 @@ def make_pkmu(kaons,
pk_docachi2_max=4.0,
pk_doca_max=0.15 * mm,
decay_descriptor="[B0 -> p+ K- mu+]cc",
name='rd_pk_muon_for_btosmutau_{hash}'):
name="rd_pk_muon_for_btosmutau_{hash}"):
"""Builder for X -> p K mu decays, used in XMuTau lines"""
combination_code = F.require_all(F.MASS < comb_m_max,
F.MAXSDOCACHI2CUT(docachi2_max),
......@@ -127,12 +128,12 @@ def make_pkmu(kaons,
def make_bs(kkmus,
muons,
pvs,
comb_m_max=7250 * MeV,
vchi2pdof_max=100,
comb_m_max=7000 * MeV,
vchi2pdof_max=75,
dira_min=0.999,
bpvfdchi2_min=16,
decay_descriptor="[B_s0 -> B0 mu-]cc",
name='rd_make_bs_to_kktaumu_{hash}'):
name="rd_make_bs_to_kktaumu_{hash}"):
"""Builder for Bs -> X mu tau(-> mu nu nu) decays, where X is K K"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
......@@ -149,12 +150,12 @@ def make_bs(kkmus,
def make_bu(kmus,
muons,
pvs,
comb_m_max=7000 * MeV,
vchi2pdof_max=75,
comb_m_max=6750 * MeV,
vchi2pdof_max=50,
dira_min=0.999,
bpvfdchi2_min=50,
decay_descriptor="[B_s0 -> B0 mu-]cc",
name='rd_make_bu_to_kmutau_{hash}'):
name="rd_make_bu_to_kmutau_{hash}"):
"""Builder for Bs -> X mu tau(-> mu nu nu) decays, where X is K"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
......@@ -171,12 +172,12 @@ def make_bu(kmus,
def make_lb(pkmus,
muons,
pvs,
comb_m_max=7500 * MeV,
vchi2pdof_max=100,
comb_m_max=7250 * MeV,
vchi2pdof_max=75,
dira_min=0.999,
bpvfdchi2_min=16,
decay_descriptor="[Lambda_b0 -> B0 mu-]cc",
name='rd_make_lb_to_pkmutau_{hash}'):
name="rd_make_lb_to_pkmutau_{hash}"):
"""Builder for Lb -> X mu tau(-> mu nu nu) decays, where X is p K"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
......@@ -190,15 +191,16 @@ def make_lb(pkmus,
return lb
@configurable
def make_bd(kstmus,
muons,
pvs,
comb_m_max=6800 * MeV,
comb_m_max=6750 * MeV,
vchi2pdof_max=50,
dira_min=0.999,
bpvfdchi2_min=25,
decay_descriptor="[B_s0 -> B0 mu+]cc",
name='rd_make_bd_to_kpimutau_{hash}'):
name="rd_make_bd_to_kpimutau_{hash}"):
"""Builder for Bd -> X mu tau(-> mu nu nu) decays, where X is K pi"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
......
Hlt/Hlt2Conf/python/Hlt2Conf/lines/rd/builders/btostautau_exclusive.py
+ 68
− 25
View file @ 3ede0aa8
......@@ -14,8 +14,10 @@ import Functors as F
from Functors.math import in_range
from Hlt2Conf.algorithms_thor import ParticleCombiner
from PyConf import configurable
@configurable
def make_phi(kaons1,
kaons2,
pvs,
......@@ -25,7 +27,7 @@ def make_phi(kaons1,
bpvfdchi2_min=16,
dira_min=0.96,
decay_descriptor="[phi(1020) -> K+ K-]cc",
name="rd_dikaons_for_btostautau"):
name="rd_dikaons_for_btostautau_{hash}"):
"""Builder for X -> K K decays"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max, F.PT > comb_pt_min,
......@@ -38,6 +40,7 @@ def make_phi(kaons1,
name=name)
@configurable
def make_kst(kaons,
pions,
pvs,
......@@ -48,7 +51,7 @@ def make_kst(kaons,
bpvfdchi2_min=75,
dira_min=0.96,
decay_descriptor="[K*(892)0 -> K+ pi-]cc",
name="rd_kpi_for_btostautau"):
name="rd_kpi_for_btostautau_{hash}"):
"""Builder for X -> K pi decays"""
combination_code = in_range(comb_m_min, F.MASS, comb_m_max)
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max, F.PT > comb_pt_min,
......@@ -61,6 +64,7 @@ def make_kst(kaons,
name=name)
@configurable
def make_lst(protons,
kaons,
pvs,
......@@ -70,7 +74,7 @@ def make_lst(protons,
bpvfdchi2_min=16,
dira_min=0.96,
decay_descriptor="[Lambda(1520)0 -> p+ K-]cc",
name="rd_pk_for_btostautau"):
name="rd_pk_for_btostautau_{hash}"):
"""Builder for X -> p K decays"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max, F.PT > comb_pt_min,
......@@ -90,7 +94,7 @@ def make_dimuon(muons1,
vchi2pdof_max=25,
dira_min=0,
decay_descriptor="[D0 -> mu+ mu-]cc",
name="rd_dimuons_for_btostautau"):
name="rd_dimuons_for_btostautau_{hash}"):
"""Builder for detached dimuons"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
......@@ -102,36 +106,75 @@ def make_dimuon(muons1,
name=name)
def make_dielectron(electrons1,
electrons2,
pvs,
comb_m_max=4500 * MeV,
vchi2pdof_max=16,
dira_min=0,
decay_descriptor="[D0 -> e+ e-]cc",
name="rd_dielectrons_for_btostautau_{hash}"):
"""Builder for detached dielectrons"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min)
return ParticleCombiner([electrons1, electrons2],
DecayDescriptor=decay_descriptor,
CombinationCut=combination_code,
CompositeCut=vertex_code,
name=name)
def make_mue(muons,
electrons,
pvs,
comb_m_max=4500 * MeV,
vchi2pdof_max=20,
dira_min=0,
decay_descriptor="[D0 -> mu+ e-]cc",
name="rd_mue_for_btostautau_{hash}"):
"""Builder for detached mue"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min)
return ParticleCombiner([muons, electrons],
DecayDescriptor=decay_descriptor,
CombinationCut=combination_code,
CompositeCut=vertex_code,
name=name)
def make_bs(phis,
dimuons,
dileptons,
pvs,
comb_m_max=7500 * MeV,
vchi2pdof_max=150,
dira_min=0.999,
name="rd_make_bs_to_kktautau"):
"""Builder for Bs -> K K tau (-> mu nu nu) tau (-> mu nu nu) decays"""
name="rd_make_bs_to_kktautau_{hash}"):
"""Builder for Bs -> K K tau (-> l nu nu) tau (-> l nu nu) decays"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min)
return ParticleCombiner([dimuons, phis],
return ParticleCombiner([dileptons, phis],
DecayDescriptor="[B_s0 -> D0 phi(1020)]cc",
CombinationCut=combination_code,
CompositeCut=vertex_code,
name=name)
@configurable
def make_bu(kaons,
dimuons,
dileptons,
pvs,
comb_m_max=7000 * MeV,
vchi2pdof_max=75,
dira_min=0.999,
name="rd_make_bu_to_ktautau"):
"""Builder for Bu -> K tau (-> mu nu nu) tau (-> mu nu nu) decays"""
name="rd_make_bu_to_ktautau_{hash}"):
"""Builder for Bu -> K tau (-> l nu nu) tau (-> l nu nu) decays"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min)
return ParticleCombiner([dimuons, kaons],
return ParticleCombiner([dileptons, kaons],
DecayDescriptor="[B+ -> D0 K+]cc",
CombinationCut=combination_code,
CompositeCut=vertex_code,
......@@ -139,18 +182,18 @@ def make_bu(kaons,
def make_lb(pks,
dimuons,
dileptons,
pvs,
comb_m_max=7750 * MeV,
vchi2pdof_max=150,
dira_min=0.999,
name="rd_make_lb_to_pktautau"):
"""Builder for Lb -> p K tau (-> mu nu nu) tau (-> mu nu nu) decays"""
name="rd_make_lb_to_pktautau_{hash}"):
"""Builder for Lb -> p K tau (-> l nu nu) tau (-> l nu nu) decays"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min)
return ParticleCombiner(
[dimuons, pks],
[dileptons, pks],
DecayDescriptor="[Lambda_b0 -> D0 Lambda(1520)0]cc",
CombinationCut=combination_code,
CompositeCut=vertex_code,
......@@ -158,17 +201,17 @@ def make_lb(pks,
def make_bd(ksts,
dimuons,
dileptons,
pvs,
comb_m_max=7250 * MeV,
vchi2pdof_max=100,
dira_min=0.999,
name="rd_make_bd_to_kpitautau"):
"""Builder for Bd -> K pi tau (-> mu nu nu) tau (-> mu nu nu) decays"""
name="rd_make_bd_to_kpitautau_{hash}"):
"""Builder for Bd -> K pi tau (-> l nu nu) tau (-> l nu nu) decays"""
combination_code = F.MASS < comb_m_max
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min)
return ParticleCombiner([dimuons, ksts],
return ParticleCombiner([dileptons, ksts],
DecayDescriptor="[B0 -> D0 K*(892)0]cc",
CombinationCut=combination_code,
CompositeCut=vertex_code,
......@@ -176,7 +219,7 @@ def make_bd(ksts,
def make_bs_to_kstkst(ksts,
dimuons,
dileptons,
pvs,
comb_m_max=7000 * MeV,
kstkst_m_max=3500 * MeV,
......@@ -184,14 +227,14 @@ def make_bs_to_kstkst(ksts,
dira_min=0.995,
docachi2_max=25.0,
decay_descriptor="[B_s0 -> K*(892)0 K*(892)~0 D0]cc",
name="rd_make_bs_to_kstksttautau"):
"""Builder for Bs -> Kst Kst tau (-> mu nu nu) tau (-> mu nu nu) decays"""
name="rd_make_bs_to_kstksttautau_{hash}"):
"""Builder for Bs -> Kst Kst tau (-> l nu nu) tau (-> l nu nu) decays"""
combination_code = F.MASS < comb_m_max
two_body_combination_code = F.require_all(F.MASS < kstkst_m_max,
F.MAXDOCACHI2CUT(docachi2_max))
F.MAXSDOCACHI2CUT(docachi2_max))
vertex_code = F.require_all(F.CHI2DOF < vchi2pdof_max,
F.BPVDIRA(pvs) > dira_min)
return ParticleCombiner([ksts, ksts, dimuons],
return ParticleCombiner([ksts, ksts, dileptons],
DecayDescriptor=decay_descriptor,
CombinationCut=combination_code,
Combination12Cut=two_body_combination_code,
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment