#!/usr/bin/env python3
##########################################################################
# basf2 (Belle II Analysis Software Framework) #
# Author: The Belle II Collaboration #
# #
# See git log for contributors and copyright holders. #
# This file is licensed under LGPL-3.0, see LICENSE.md. #
##########################################################################
"""
Skim list building functions for :math:`B\\to X_u +h` analyses
"""
import modularAnalysis as ma
from skim import BaseSkim, fancy_skim_header
from stdV0s import stdKshorts
from skim.standardlists.charmless import (
loadStdVeryLooseTracks,
loadStdPi0ForBToCharmless,
loadStdVeryLooseKstar0,
loadStdVeryLooseRho0,
loadStdVeryLooseKstarPlus,
loadStdVeryLooseKstarPlusPi0,
loadStdVeryLooseRhoPlus)
from skim.standardlists.lightmesons import (
loadStdSkimHighEffTracks,
loadStdSkimHighEffEtaPrime,
loadStdSkimHighEffEta)
from stdPi0s import stdPi0s
from stdPhotons import stdPhotons
__liaison__ = "Benedikt Wach <benedikt.wach@desy.de>"
_VALIDATION_SAMPLE = "mdst16.root"
[docs]
@fancy_skim_header
class BtoPi0Pi0(BaseSkim):
"""
Reconstructed decay modes:
* :math:`B^{0}\\to \\pi^0 \\pi^0`
Cuts applied:
* ``5.20 < Mbc < 5.29``
* ``abs(deltaE) < 0.5``
"""
__authors__ = ["Fernando Abudinen", "Riccardo Manfredi", "Sebastiano Raiz", "Francis Pham", "Benedikt Wach"]
__description__ = "Skim list definitions for neutral B to two neutral pions."
__contact__ = __liaison__
__category__ = "physics, hadronic B to charmless"
ApplyHLTHadronCut = False
NoisyModules = ["ParticleLoader", "RootOutput"]
[docs]
def load_standard_lists(self, path):
loadStdPi0ForBToCharmless(path=path)
[docs]
def build_lists(self, path):
Bcuts = '5.20 < Mbc < 5.29 and abs(deltaE) < 0.5'
BsigList = []
ma.reconstructDecay('B0:Pi0Pi0 -> pi0:charmlessFit pi0:charmlessFit', Bcuts, path=path)
BsigList.append('B0:Pi0Pi0')
return BsigList
[docs]
@fancy_skim_header
class BtoPi0Eta(BaseSkim):
"""
Reconstructed decay modes:
* :math:`B^{0}\\to \\pi^0 \\eta`
* :math:`B^{0}\\to \\pi^0 \\eta^{'}`
Cuts applied:
* ``5.20 < Mbc < 5.29``
* ``abs(deltaE) < 0.5``
"""
__authors__ = ["Longke Li", "Yinghui Guan"]
__description__ = "Skim list definitions for neutral B to pi0 eta and B to pi0 eta'."
__contact__ = __liaison__
__category__ = "physics, hadronic B to charmless"
ApplyHLTHadronCut = False
NoisyModules = ["ParticleLoader", "RootOutput"]
[docs]
def load_standard_lists(self, path):
loadStdPi0ForBToCharmless(path=path)
stdPhotons('tight', path=path)
stdPi0s("eff40_May2020", path=path)
loadStdSkimHighEffTracks('pi', path=path)
loadStdSkimHighEffEta(path=path)
loadStdSkimHighEffEtaPrime(path=path)
[docs]
def build_lists(self, path):
Bcuts = '5.20 < Mbc < 5.29 and abs(deltaE) < 0.5'
BsigList = []
ma.reconstructDecay('B0:Pi0Eta -> pi0:charmlessFit eta:SkimHighEff', Bcuts, path=path)
ma.reconstructDecay('B0:Pi0Etap -> pi0:charmlessFit eta\':SkimHighEff', Bcuts, path=path)
BsigList.append('B0:Pi0Eta')
BsigList.append('B0:Pi0Etap')
return BsigList
[docs]
@fancy_skim_header
class BtoHadTracks(BaseSkim):
"""
Reconstructed decay modes:
* :math:`B^{0}\\to \\pi^+ \\pi^-`
* :math:`B^{0}\\to K^+ \\pi^-`
* :math:`B^{0}\\to K^+ K^-`
* :math:`B^{+}\\to \\pi^+ \\pi^+ \\pi^-`
* :math:`B^{+}\\to K_{\\rm S}^0 \\pi^+`
* :math:`B^{+}\\to K_{\\rm S}^0 K^+`
* :math:`B^{+}\\to K^+ \\pi^+ \\pi^-`
* :math:`B^{+}\\to K^+ K^+ \\pi^-`
* :math:`B^{+}\\to K^+ K^- \\pi^+`
* :math:`B^{+}\\to K^+ K^+ K^-`
* :math:`B^{0}\\to \\pi^+ \\pi^- \\pi^+ \\pi^-`
* :math:`B^{0}\\to K^+ \\pi^- \\pi^+ \\pi^-`
* :math:`B^{0}\\to K^+ K^- \\pi^+ \\pi^-`
* :math:`B^{0}\\to K^+ \\pi^- K^+ \\pi^-`
* :math:`B^{0}\\to K^+ K^- K^+ \\pi^-`
* :math:`B^{0}\\to K^+ K^- K^+ K^-`
* :math:`B^{0}\\to K_{\\rm S}^0 \\pi^+ \\pi^-`
* :math:`B^{0}\\to K_{\\rm S}^0 K^+ \\pi^-`
* :math:`B^{0}\\to K_{\\rm S}^0 K^+ K^-`
* :math:`B^{0}\\to K_{\\rm S}^0 K_{\\rm S}^0`
* :math:`B^{0}\\to K^{*0} K^- K^+`
* :math:`B^{0}\\to K^{*0} \\pi^- \\pi^+`
* :math:`B^{0}\\to K^{*0} K^+ \\pi^-`
* :math:`B^{0}\\to \\rho^0 \\rho^0`
* :math:`B^{+}\\to K^{*+} K^- K^+`
* :math:`B^{+}\\to K^{*+} \\pi^- \\pi^+`
* :math:`B^{+}\\to K^{*+} K^+ \\pi^-`
Cuts applied:
* ``nTracks >= 2 (event level)``
* ``5.20 < Mbc < 5.29``
* ``abs(deltaE) < 0.3``
"""
__authors__ = ["Fernando Abudinen", "Riccardo Manfredi", "Sebastiano Raiz", "Benedikt Wach"]
__description__ = "Skim list definitions for all B to charmless modes with tracks only in the final state."
__contact__ = __liaison__
__category__ = "physics, hadronic B to charmless"
ApplyHLTHadronCut = True
NoisyModules = ["ParticleLoader", "RootOutput"]
[docs]
def load_standard_lists(self, path):
loadStdVeryLooseTracks('K', path=path)
loadStdVeryLooseTracks('pi', path=path)
stdKshorts(path=path)
loadStdVeryLooseKstar0(path=path)
loadStdVeryLooseRho0(path=path)
loadStdVeryLooseKstarPlus(path=path)
[docs]
def build_lists(self, path):
Bcuts = "5.2 < Mbc < 5.29 and abs(deltaE) < 0.3"
channels = {'B0:Charmless_b2pipi': 'pi+:SkimVeryLoose pi-:SkimVeryLoose', # 1
'B0:Charmless_b2Kpi': 'K+:SkimVeryLoose pi-:SkimVeryLoose', # 2
'B0:Charmless_b2KK': 'K+:SkimVeryLoose K-:SkimVeryLoose', # 3
'B+:Charmless_b2pipipi': 'pi+:SkimVeryLoose pi+:SkimVeryLoose pi-:SkimVeryLoose', # 4
'B+:Charmless_b2Kspi': 'K_S0:merged pi+:SkimVeryLoose', # 5
'B+:Charmless_b2KsK': 'K_S0:merged K+:SkimVeryLoose', # 6
'B+:Charmless_b2Kpipi': 'K+:SkimVeryLoose pi+:SkimVeryLoose pi-:SkimVeryLoose ', # 7
'B+:Charmless_b2KKpi_0': 'K+:SkimVeryLoose K+:SkimVeryLoose pi-:SkimVeryLoose', # 8
'B+:Charmless_b2KKpi_1': 'K+:SkimVeryLoose K-:SkimVeryLoose pi+:SkimVeryLoose', # 9
'B+:Charmless_b2KKK': 'K+:SkimVeryLoose K+:SkimVeryLoose K-:SkimVeryLoose', # 10
'B0:Charmless_b2pipipipi': 'pi+:SkimVeryLoose pi-:SkimVeryLoose pi+:SkimVeryLoose pi-:SkimVeryLoose', # 11
'B0:Charmless_b2Kpipipi': 'K+:SkimVeryLoose pi-:SkimVeryLoose pi+:SkimVeryLoose pi-:SkimVeryLoose', # 12
'B0:Charmless_b2KKpipi': 'K+:SkimVeryLoose K-:SkimVeryLoose pi+:SkimVeryLoose pi-:SkimVeryLoose', # 13
'B0:Charmless_b2KpiKpi': 'K+:SkimVeryLoose pi-:SkimVeryLoose K+:SkimVeryLoose pi-:SkimVeryLoose', # 14
'B0:Charmless_b2KKKpi': 'K+:SkimVeryLoose K-:SkimVeryLoose K+:SkimVeryLoose pi-:SkimVeryLoose', # 15
'B0:Charmless_b2KKKK': 'K+:SkimVeryLoose K-:SkimVeryLoose K+:SkimVeryLoose K-:SkimVeryLoose', # 16
'B0:Charmless_b2Kspipi': 'K_S0:merged pi+:SkimVeryLoose pi-:SkimVeryLoose', # 17
'B0:Charmless_b2KsKpi': 'K_S0:merged K+:SkimVeryLoose pi-:SkimVeryLoose', # 18
'B0:Charmless_b2KsKK': 'K_S0:merged K+:SkimVeryLoose K-:SkimVeryLoose', # 19
'B0:Charmless_b2KsKs': 'K_S0:merged K_S0:merged', # 20
'B0:Charmless_b2K*KK': 'K*0:veryLoose K+:SkimVeryLoose K-:SkimVeryLoose', # 21
'B0:Charmless_b2K*pipi': 'K*0:veryLoose pi+:SkimVeryLoose pi-:SkimVeryLoose', # 22
'B0:Charmless_b2K*Kpi': 'K*0:veryLoose K+:SkimVeryLoose pi-:SkimVeryLoose', # 23
'B0:Charmless_b2rhorho': 'rho0:veryLoose rho0:veryLoose', # 24
'B+:Charmless_b2K*pKK': 'K*+:veryLoose K+:SkimVeryLoose K-:SkimVeryLoose', # 25
'B+:Charmless_b2K*ppipi': 'K*+:veryLoose pi+:SkimVeryLoose pi-:SkimVeryLoose', # 26
'B+:Charmless_b2K*pKpi': 'K*+:veryLoose K+:SkimVeryLoose pi-:SkimVeryLoose', # 27
}
BsigList = []
for chID, channel in enumerate(channels.keys()):
ma.reconstructDecay(decayString=channel + ' -> ' + channels[channel],
cut=Bcuts, dmID=chID, path=path)
BsigList.append(channel)
path = self.skim_event_cuts("nTracks >= 2", path=path)
return BsigList
[docs]
@fancy_skim_header
class BtoHad1Pi0(BaseSkim):
"""
Reconstructed decay modes:
* :math:`B^+\\to \\pi^+ \\pi^0`
* :math:`B^+\\to K^+ \\pi^0`
* :math:`B^0\\to \\pi^+ \\pi^- \\pi^0`
* :math:`B^0\\to K^+ \\pi^- \\pi^0`
* :math:`B^0\\to K^- \\pi^+ \\pi^0`
* :math:`B^0\\to K^+ K^- \\pi^0`
* :math:`B^0\\to K_{\\rm S}^0 \\pi^0`
* :math:`B^+\\to K_{\\rm S}^{0} \\pi^+ \\pi^0`
* :math:`B^+\\to K_{\\rm S}^{0} K^+ \\pi^0`
* :math:`B^+\\to \\rho^{+} \\rho^{0}`
Cuts applied:
* ``5.20 < Mbc < 5.29``
* ``abs(deltaE) < 0.5``
"""
__authors__ = ["Fernando Abudinen", "Riccardo Manfredi", "Sebastiano Raiz", "Benedikt Wach"]
__description__ = "Skim list for B to charmless modes including 1 pi0 in the final state."
__contact__ = __liaison__
__category__ = "physics, hadronic B to charmless"
ApplyHLTHadronCut = True
NoisyModules = ["ParticleLoader", "RootOutput"]
[docs]
def load_standard_lists(self, path):
loadStdVeryLooseTracks('K', path=path)
loadStdVeryLooseTracks('pi', path=path)
loadStdPi0ForBToCharmless(path=path)
stdKshorts(path=path)
loadStdVeryLooseRho0(path=path)
loadStdVeryLooseRhoPlus(path=path)
[docs]
def build_lists(self, path):
Bcuts = '5.20 < Mbc < 5.29 and abs(deltaE) < 0.5'
BsigList = []
channels = {
'B+:Charmless_b2pipi0': 'pi+:SkimVeryLoose pi0:charmlessFit', # 1
'B+:Charmless_b2Kpi0': 'K+:SkimVeryLoose pi0:charmlessFit', # 2
'B0:Charmless_b2pipipi0': 'pi+:SkimVeryLoose pi-:SkimVeryLoose pi0:charmlessFit', # 1
'B0:Charmless_b2Kpipi0_0': 'K+:SkimVeryLoose pi-:SkimVeryLoose pi0:charmlessFit', # 2
'B0:Charmless_b2KKpi0': 'K+:SkimVeryLoose K-:SkimVeryLoose pi0:charmlessFit', # 4
'B0:Charmless_b2Kspi0': 'K_S0:merged pi0:charmlessFit', # 5
'B+:Charmless_b2Kspipi0': 'K_S0:merged pi+:SkimVeryLoose pi0:charmlessFit', # 6
'B+:Charmless_b2KsKpi0': 'K_S0:merged K+:SkimVeryLoose pi0:charmlessFit', # 7
'B+:Charmless_b2rhorho0': 'rho+:veryLoose rho0:veryLoose', # 8
}
for chID, channel in enumerate(channels.keys()):
ma.reconstructDecay(decayString=channel + ' -> ' + channels[channel],
cut=Bcuts, dmID=chID, path=path)
BsigList.append(channel)
return BsigList
[docs]
@fancy_skim_header
class BtoHad3Tracks1Pi0(BaseSkim):
"""
Reconstructed decay modes:
* :math:`B^+\\to K^{*+} K^+ K^-`
* :math:`B^+\\to K^{*+} K^+ \\pi^-`
* :math:`B^+\\to K^{*+} \\pi^+ \\pi^-`
Cuts applied:
* ``5.20 < Mbc < 5.29``
* ``abs(deltaE) < 0.5``
"""
__authors__ = ["Fernando Abudinen", "Riccardo Manfredi", "Sebastiano Raiz", "Benedikt Wach"]
__description__ = "Skim list for all charged B to charmless modes with 3 tracks and 1 Pi0."
__contact__ = __liaison__
__category__ = "physics, hadronic B to charmless"
ApplyHLTHadronCut = True
NoisyModules = ["ParticleLoader", "RootOutput"]
[docs]
def load_standard_lists(self, path):
loadStdVeryLooseTracks('K', path=path)
loadStdVeryLooseTracks('pi', path=path)
loadStdPi0ForBToCharmless(path=path)
loadStdVeryLooseKstarPlusPi0(path=path)
[docs]
def build_lists(self, path):
Bcuts = '5.20 < Mbc < 5.29 and abs(deltaE) < 0.5'
BsigList = []
channels = {
'B+:Charmless_b2K*KK_2': 'K*+:veryLoosePi0 K-:SkimVeryLoose K+:SkimVeryLoose', # 1
'B+:Charmless_b2K*pipi_2': 'K*+:veryLoosePi0 pi-:SkimVeryLoose pi+:SkimVeryLoose', # 2
'B+:Charmless_b2K*Kpi_2': 'K*+:veryLoosePi0 K+:SkimVeryLoose pi-:SkimVeryLoose', # 3
}
for chID, channel in enumerate(channels.keys()):
ma.reconstructDecay(decayString=channel + ' -> ' + channels[channel],
cut=Bcuts, dmID=chID, path=path)
ma.rankByHighest(particleList=channel, variable="cos(mdstIndex)", numBest=3,
outputVariable="cosMdstIndex_rank", path=path)
BsigList.append(channel)
path = self.skim_event_cuts("nTracks >= 2", path=path)
return BsigList
[docs]
@fancy_skim_header
class BtoRhopRhom(BaseSkim):
"""
Reconstructed decay modes:
* :math:`B^+ \\to \\rho^{+} \\rho^{-}`
Cuts applied:
* ``5.20 < Mbc < 5.29``
* ``abs(deltaE) < 0.5``
"""
__authors__ = ["Yu Nakazawa", "Benedikt Wach"]
__description__ = "Skim list for B to rho+ rho-."
__contact__ = __liaison__
__category__ = "physics, hadronic B to charmless"
ApplyHLTHadronCut = True
NoisyModules = ["ParticleLoader", "RootOutput"]
[docs]
def load_standard_lists(self, path):
loadStdVeryLooseTracks('pi', path=path)
loadStdPi0ForBToCharmless(path=path)
loadStdVeryLooseRhoPlus(path=path)
[docs]
def build_lists(self, path):
Bcuts = '5.20 < Mbc < 5.29 and abs(deltaE) < 0.5'
BsigList = []
ma.reconstructDecay('B0:Charmless_b2rr -> rho+:veryLoose rho-:veryLoose', Bcuts, path=path)
BsigList.append('B0:Charmless_b2rr')
path = self.skim_event_cuts("nTracks >= 2", path=path)
return BsigList
[docs]
@fancy_skim_header
class BtoEtapKstp(BaseSkim):
"""
Reconstructed decay mode:
* :math:`B^{+}\\to \\eta^{'} K^{*+}`
Cuts applied:
* ``5.20 < Mbc < 5.29``
* ``abs(deltaE) < 0.5``
* ``0.4 < M_{eta} < 0.6``
* ``0.8 < M_{eta'} < 1.1``
* ``0.7 < M_{K*} < 1.6``
Note:
This skim uses `skim.standardlists.lightmesons.loadStdSkimHighEffEtaPrime`
and `skim.standardlists.charmless.loadStdVeryLooseKstarPlus` where
:math:`\\eta^{'}` and :math:`K^{*+}` are defined, respectively.
:math:`\\eta` channel is defined in `skim.standardlists.lightmesons.loadStdSkimHighEffEta`.
"""
__authors__ = ["Chia-Ling Hsu"]
__description__ = "Skim list for B+ to eta' K*+."
__contact__ = __liaison__
__category__ = "physics, hadronic B to charmless"
ApplyHLTHadronCut = True
NoisyModules = ["ParticleLoader", "RootOutput"]
[docs]
def load_standard_lists(self, path):
stdKshorts(path=path)
stdPhotons('tight', path=path)
stdPi0s("eff40_May2020", path=path)
loadStdVeryLooseTracks('pi', path=path)
loadStdSkimHighEffTracks('pi', path=path)
loadStdSkimHighEffEta(path=path)
loadStdSkimHighEffEtaPrime(path=path)
loadStdVeryLooseKstarPlus(path=path)
[docs]
def build_lists(self, path):
Bcuts = '5.20 < Mbc < 5.29 and abs(deltaE) < 0.5'
BsigList = []
ma.reconstructDecay("B+:Charmless_b2etapkst -> eta':SkimHighEff K*+:veryLoose", Bcuts, path=path)
BsigList.append('B+:Charmless_b2etapkst')
return BsigList
[docs]
@fancy_skim_header
class BtoEtapRhop(BaseSkim):
"""
Reconstructed decay mode:
* :math:`B^{+}\\to \\eta^{'} (\\to \\pi^{+} \\pi^{-}
\\eta (\\to \\gamma \\gamma) ) \\rho^{+} (\\pi^{+} \\pi^{0})`
* :math:`B^{+}\\to \\eta^{'} (\\to \\rho^{0} (\\to \\pi^{+}
\\pi^{-}) \\gamma ) \\rho^{+} (\\pi^{+} \\pi^{0})`
Cuts applied:
* ``5.20 < Mbc < 5.29``
* ``abs(deltaE) < 0.5``
* ``0.8< M_{etap} < 1.1``
* ``0.4< M_{eta} < 0.6``
* ``0.47 < M_{rhop} < 1.15``
* ``0.47 < M_{rho0} < 1.15``
"""
__authors__ = ["Ruei-Si Hsiao"]
__description__ = "Skim list for B+ to etap rho+."
__contact__ = __liaison__
__category__ = "physics, hadronic B to charmless"
ApplyHLTHadronCut = True
NoisyModules = ["ParticleLoader", "RootOutput"]
[docs]
def load_standard_lists(self, path):
stdPhotons('tight', path=path)
stdPi0s("eff40_May2020", path=path)
loadStdSkimHighEffTracks('pi', path=path)
[docs]
def build_lists(self, path):
Bcuts = '5.20 < Mbc < 5.29 and abs(deltaE) < 0.5'
BsigList = []
ma.reconstructDecay("eta:b2etaprhop_gamgam -> gamma:tight gamma:tight", '0.4 < M < 0.6', path=path)
ma.reconstructDecay('rho0:b2etaprhop_pipi -> pi+:SkimHighEff pi-:SkimHighEff',
'0.47 < M < 1.15', path=path)
ma.reconstructDecay("eta\':b2etaprhop_loose1 -> pi+:SkimHighEff pi-:SkimHighEff eta:b2etaprhop_gamgam",
'0.8 < M < 1.1', 1, path=path)
ma.reconstructDecay("eta\':b2etaprhop_loose2 -> rho0:b2etaprhop_pipi gamma:tight",
'0.8 < M < 1.1', 2, path=path)
ma.copyLists("eta\':b2etaprhop_loose", ["eta\':b2etaprhop_loose1", "eta\':b2etaprhop_loose2"],
path=path)
ma.reconstructDecay('rho+:b2etaprhop_loose -> pi+:SkimHighEff pi0:eff40_May2020', '0.47 < M < 1.15',
path=path)
ma.reconstructDecay("B+:Charmless_b2etaprhop -> eta\':b2etaprhop_loose rho+:b2etaprhop_loose",
Bcuts, path=path)
BsigList.append("B+:Charmless_b2etaprhop")
return BsigList
