#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""(Semi-)Leptonic Working Group Skim list building functions for semi-leptonic analyses.
"""
__authors__ = [
"Racha Cheaib",
"Hannah Wakeling",
"Phil Grace"
]
import modularAnalysis as ma
from skim.standardlists.charm import (loadKForBtoHadrons, loadPiForBtoHadrons,
loadStdD0, loadStdDplus, loadStdDstar0,
loadStdDstarPlus, loadPiSkimHighEff,
loadKSkimHighEff, loadSlowPi,
loadSkimHighEffD0_Kpi, loadSkimHighEffDstarPlus_D0pi_Kpi,
loadSkimHighEffD0_Kpipipi, loadSkimHighEffDstarPlus_D0pi_Kpipipi,
loadStdD0_eff20_Kpipi0, loadStdDstarPlus_D0pi_Kpipi0_eff20)
from skim.standardlists.lightmesons import loadStdPi0ForBToHadrons
from skimExpertFunctions import BaseSkim, fancy_skim_header
from stdCharged import stdE, stdK, stdMu, stdPi
from stdPhotons import stdPhotons
from stdPi0s import stdPi0s
from stdV0s import stdKshorts
from variables import variables as vm
__liaison__ = "Shanette De La Motte <shanette.delamotte@adelaide.edu.au>"
[docs]@fancy_skim_header
class PRsemileptonicUntagged(BaseSkim):
"""
Reconstructed decay modes:
* :math:`B^0 \\to \\pi^- e^+`
* :math:`B^0 \\to \\pi^- \\mu^+`
Event-level cuts:
* :math:`\\text{foxWolframR2} > 0.5` constructed using tracks with
:math:`p_T>0.1\\,\\text{GeV}` and clusters with :math:`E>0.1\\,\\text{GeV}`.
* :math:`n_{\\text{tracks}} > 4`
Cuts on electrons:
* :math:`\\text{electronID} > 0.5`
* :math:`p > 1.5\\,\\text{GeV}` in CMS frame
Cuts on muons:
* :math:`\\text{muonID} > 0.5`
* :math:`p > 1.5\\,\\text{GeV}` in CMS frame
Cuts on pions:
* :math:`\\text{pionID}>0.5`
* :math:`\\text{muonID}<0.2`
* :math:`0.060\\,\\text{GeV}<p<0.220\\,\\text{GeV}` in CMS frame
Cuts on partially reconstructed :math:`B` mesons:
* :math:`\\cos\\theta_{\\ell,\\,\\pi}<0` in CMS frame.
"""
__authors__ = ["Lucien Cremaldi", "Racha Cheaib", "Romulus Godang"]
__description__ = "Skim for partial reconstruction analysis in leptonic group."
__contact__ = __liaison__
__category__ = "physics, semileptonic"
[docs] def load_standard_lists(self, path):
stdE("all", path=path)
stdMu("all", path=path)
stdPi("all", path=path)
[docs] def build_lists(self, path):
ma.fillParticleList(decayString="pi+:PRSL_eventshape",
cut="pt> 0.1", path=path)
ma.fillParticleList(decayString="gamma:PRSL_eventshape",
cut="E > 0.1", path=path)
ma.buildEventShape(inputListNames=["pi+:PRSL_eventshape", "gamma:PRSL_eventshape"],
allMoments=False,
foxWolfram=True,
harmonicMoments=False,
cleoCones=False,
thrust=False,
collisionAxis=False,
jets=False,
sphericity=False,
checkForDuplicates=False,
path=path)
path = self.skim_event_cuts("foxWolframR2<0.5 and nTracks>4", path=path)
ma.cutAndCopyList("e+:PRSemileptonic_1", "e+:all",
"useCMSFrame(p) > 1.50 and electronID > 0.5", path=path)
ma.cutAndCopyList("mu+:PRSemileptonic_1", "mu+:all",
"useCMSFrame(p) > 1.50 and muonID > 0.5", path=path)
ma.cutAndCopyList("pi-:PRSemileptonic_1", "pi-:all",
"pionID>0.5 and muonID<0.2 and 0.060<useCMSFrame(p)<0.220", path=path)
ma.cutAndCopyList("e+:PRSemileptonic_2", "e+:all",
"0.600 < useCMSFrame(p) <= 1.50 and electronID > 0.5", path=path)
ma.cutAndCopyList("mu+:PRSemileptonic_2", "mu+:all",
"0.350 < useCMSFrame(p) <= 1.50 and muonID > 0.5", path=path)
ma.cutAndCopyList("pi-:PRSemileptonic_2", "pi-:all",
"pionID>0.5 and muonID<0.2 and 0.060<useCMSFrame(p)<0.160", path=path)
ma.reconstructDecay("B0:PRSemileptonic_1 -> pi-:PRSemileptonic_1 e+:PRSemileptonic_1",
"useCMSFrame(cos(daughterAngle(0,1)))<0.00", 1, path=path)
ma.reconstructDecay("B0:PRSemileptonic_2 -> pi-:PRSemileptonic_1 mu+:PRSemileptonic_1",
"useCMSFrame(cos(daughterAngle(0,1)))<0.00", 2, path=path)
ma.reconstructDecay("B0:PRSemileptonic_3 -> pi-:PRSemileptonic_2 e+:PRSemileptonic_2",
"useCMSFrame(cos(daughterAngle(0,1)))<1.00", 3, path=path)
ma.reconstructDecay("B0:PRSemileptonic_4 -> pi-:PRSemileptonic_2 mu+:PRSemileptonic_2",
"useCMSFrame(cos(daughterAngle(0,1)))<1.00", 4, path=path)
self.SkimLists = ["B0:PRSemileptonic_1", "B0:PRSemileptonic_2"]
[docs] def validation_histograms(self, path):
# NOTE: the validation package is not part of the light releases, so this import
# must be made here rather than at the top of the file.
from validation_tools.metadata import create_validation_histograms
ma.cutAndCopyLists("B0:PRSemileptonic_semileptonic",
["B0:PRSemileptonic_1", "B0:PRSemileptonic_2"], "", path=path)
ma.buildRestOfEvent("B0:PRSemileptonic_semileptonic", path=path)
ma.appendROEMask("B0:PRSemileptonic_semileptonic", "basic",
"pt>0.05 and -2<dr<2 and -4.0<dz<4.0",
"E>0.05",
path=path)
ma.buildContinuumSuppression("B0:PRSemileptonic_semileptonic", "basic", path=path)
vm.addAlias("d0_p", "daughter(0, p)")
vm.addAlias("d1_p", "daughter(1, p)")
vm.addAlias("MissM2", "weMissM2(basic,0)")
histogramFilename = f"{self}_Validation.root"
email = "Phil Grace <philip.grace@adelaide.edu.au>"
create_validation_histograms(
rootfile=histogramFilename,
particlelist="B0:PRSemileptonic",
variables_1d=[
("Mbc", 100, 4.0, 5.3, "Mbc", email, "", ""),
("d0_p", 100, 0, 5.2, "Signal-side pion momentum", email, "", ""),
("d1_p", 100, 0, 5.2, "Signal-side lepton momentum", email, "", ""),
("MissM2", 100, -5, 5, "Missing mass squared", email, "", "")
],
variables_2d=[("deltaE", 100, -5, 5, "Mbc", 100, 4.0, 5.3, "Mbc vs deltaE", email, "", "")],
path=path)
[docs]@fancy_skim_header
class SLUntagged(BaseSkim):
"""
Cuts applied:
* :math:`p_{\\ell} > 0.35\\,\\text{GeV}`
* :math:`5.24 < M_{\\text{bc}} < 5.29`
* :math:`|\\Delta E | < 0.5`
* :math:`n_{\\text{tracks}} > 4`
Reconstructed decays:
* :math:`B^+ \\to \\overline{D}^{0} e^+`
* :math:`B^+ \\to \\overline{D}^{0} \\mu^+`
* :math:`B^+ \\to \\overline{D}^{*0} e^+`
* :math:`B^+ \\to \\overline{D}^{*0} \\mu^+`
* :math:`B^0 \\to D^{-} e^+`
* :math:`B^0 \\to D^{-} \\mu^+`
* :math:`B^0 \\to D^{*-} e^+`
* :math:`B^0 \\to D^{*-} \\mu^+`
"""
__authors__ = ["Phillip Urquijo", "Racha Cheaib"]
__description__ = (
"Skim for semileptonic decays, :math:`B` decays "
"(:math:`B \\to D \\ell\\nu`, where :math:`\\ell=e,\\mu`)"
)
__contact__ = __liaison__
__category__ = "physics, semileptonic"
[docs] def load_standard_lists(self, path):
stdE("all", path=path)
stdK("all", path=path)
stdMu("all", path=path)
stdPi("all", path=path)
stdPi("loose", path=path)
stdPhotons("loose", path=path)
stdPi0s("eff40_Jan2020", path=path)
stdKshorts(path=path)
loadStdPi0ForBToHadrons(path=path)
loadPiForBtoHadrons(path=path)
loadKForBtoHadrons(path=path)
loadStdD0(path=path)
loadStdDstar0(path=path)
loadStdDplus(path=path)
loadStdDstarPlus(path=path)
[docs] def build_lists(self, path):
ma.cutAndCopyList("e+:SLUntagged", "e+:all", "p>0.35", True, path=path)
ma.cutAndCopyList("mu+:SLUntagged", "mu+:all", "p>0.35", True, path=path)
Bcuts = "5.24 < Mbc < 5.29 and abs(deltaE) < 0.5"
BplusChannels = ["anti-D0:all e+:SLUntagged",
"anti-D0:all mu+:SLUntagged",
"anti-D*0:all e+:SLUntagged",
"anti-D*0:all mu+:SLUntagged"
]
B0Channels = ["D-:all e+:SLUntagged",
"D-:all mu+:SLUntagged",
"D*-:all e+:SLUntagged",
"D*-:all mu+:SLUntagged"
]
bplusList = []
for chID, channel in enumerate(BplusChannels):
ma.reconstructDecay(f"B+:SLUntagged_{chID} -> {channel}", Bcuts, chID, path=path)
ma.applyCuts(f"B+:SLUntagged_{chID}", "nTracks>4", path=path)
bplusList.append(f"B+:SLUntagged_{chID}")
b0List = []
for chID, channel in enumerate(B0Channels):
ma.reconstructDecay(f"B0:SLUntagged_{chID} -> {channel}", Bcuts, chID, path=path)
ma.applyCuts(f"B0:SLUntagged_{chID}", "nTracks>4", path=path)
b0List.append(f"B0:SLUntagged_{chID}")
self.SkimLists = b0List + bplusList
[docs] def validation_histograms(self, path):
# NOTE: the validation package is not part of the light releases, so this import
# must be made here rather than at the top of the file.
from validation_tools.metadata import create_validation_histograms
ma.cutAndCopyLists("B+:SLUntagged",
["B+:SLUntagged_0", "B+:SLUntagged_1", "B+:SLUntagged_2", "B+:SLUntagged_3"],
"", path=path)
ma.buildRestOfEvent("B+:SLUntagged", path=path)
ma.appendROEMask("B+:SLUntagged", "basic",
"pt>0.05 and -2<dr<2 and -4.0<dz<4.0",
"E>0.05",
path=path)
ma.buildContinuumSuppression("B+:SLUntagged", "basic", path=path)
vm.addAlias("d1_p", "daughter(1,p)")
vm.addAlias("MissM2", "weMissM2(basic,0)")
histogramFilename = f"{self}_Validation.root"
myEmail = "Phil Grace <philip.grace@adelaide.edu.au>"
create_validation_histograms(
rootfile=histogramFilename,
particlelist="B+:SLUntagged",
variables_1d=[
("cosThetaBetweenParticleAndNominalB", 100, -6.0, 4.0, "cosThetaBY", myEmail, "", ""),
("Mbc", 100, 4.0, 5.3, "Mbc", myEmail, "", ""),
("d1_p", 100, 0, 5.2, "Signal-side lepton momentum", myEmail, "", ""),
("MissM2", 100, -5, 5, "Missing mass squared", myEmail, "", "")
],
variables_2d=[("deltaE", 100, -5, 5, "Mbc", 100, 4.0, 5.3, "Mbc vs deltaE", myEmail, "", "")],
path=path)
[docs]@fancy_skim_header
class B0toDstarl_Kpi_Kpipi0_Kpipipi(BaseSkim):
"""
Cuts applied:
* ``SkimHighEff tracks thetaInCDCAcceptance AND abs(dr) < 2 AND abs(dz) < 5 AND PID>=0.01``
* ``slowPi tracks thetaInCDCAcceptance AND abs(dr) < 2 AND abs(dz) < 5 AND useCMSFrame(p) < 0.4``
* :math:`2.5 > p_{\\ell} > 1.1\\,\\text{GeV}`
* ``lepton with abs(d0) < 0.5 AND abs(z0) < 2 AND thetaInCDCAcceptance AND ID >= 0.95 AND 1.1 < useCMSFrame(p) < 2.5``
* ``1.8 < M_D0 < 2.0``
* ``DM_Dstar_D < 0.16``
Reconstructed decays:
* :math:`B^{0}\\to D^{*-} (D^{0} \\to K^+ \\pi^-) e^+`,
* :math:`B^{0}\\to D^{*-} (D^{0} \\to K^+ \\pi^- \\pi^0) e^+`,
* :math:`B^{0}\\to D^{*-} (D^{0} \\to K^+ \\pi^- \\pi^- \\pi^+) e^+`,
* :math:`B^{0}\\to D^{*-} (D^{0} \\to K^+ \\pi^-) mu^+`,
* :math:`B^{0}\\to D^{*-} (D^{0} \\to K^+ \\pi^- \\pi^0) mu^+`,
* :math:`B^{0}\\to D^{*-} (D^{0} \\to K^+ \\pi^- \\pi^- \\pi^+) mu^+`,
Note:
This skim uses `skim.standardlists.charm.loadSkimHighEffD0_Kpi`,
`skim.standardlists.charm.loadSkimHighEffD0_Kpipipi` and
`skim.standardlists.charm.loadStdD0_eff20_Kpipi0`, where :math:`D^0`
channel is defined.
`skim.standardlists.charm.loadSkimHighEffDstarPlus_D0pi_Kpi`,
`skim.standardlists.charm.loadSkimHighEffDstarPlus_D0pi_Kpipipi`,
`skim.standardlists.charm.loadStdDstarPlus_D0pi_Kpipi0_eff20`,where the
:math:`D^{*-}` channel is defined.
The pion and kaon lists used to define :math:`D^0` and :math:`D^{*-}` are:
`skim.standardlists.charm.loadPiSkimHighEff`, `skim.standardlists.charm.loadKSkimHighEff` and
`skim.standardlists.charm.loadSlowPi`
"""
__authors__ = ["Bae Hanwook, Chiara La Licata"]
__description__ = ""
__contact__ = __liaison__
__category__ = "physics, semileptonic"
ApplyHLTHadronCut = True
produce_on_tau_samples = False # retention is very close to zero on taupair
[docs] def load_standard_lists(self, path):
stdE("all", path=path)
stdMu("all", path=path)
stdPi("all", path=path)
stdK("all", path=path)
stdPi0s("eff20_Jan2020", path=path)
loadPiSkimHighEff(path=path)
loadKSkimHighEff(path=path)
loadSlowPi(path=path)
loadSkimHighEffD0_Kpi(path=path)
loadSkimHighEffDstarPlus_D0pi_Kpi(path=path)
loadSkimHighEffD0_Kpipipi(path=path)
loadSkimHighEffDstarPlus_D0pi_Kpipipi(path=path)
loadStdD0_eff20_Kpipi0(path=path)
loadStdDstarPlus_D0pi_Kpipi0_eff20(path=path)
[docs] def build_lists(self, path):
ma.cutAndCopyList(
'e+:sig',
'e+:all',
'abs(d0) < 0.5 and abs(z0) < 2 and thetaInCDCAcceptance and electronID >= 0.95 and 1.1 < useCMSFrame(p) < 2.5 ',
path=path)
ma.cutAndCopyList(
'mu+:sig',
'mu+:all',
'abs(d0) < 0.5 and abs(z0) < 2 and thetaInCDCAcceptance and muonID >= 0.95 and 1.1 < useCMSFrame(p) < 2.5',
path=path)
B0_channels = ["D*-:D0_Kpi_skimhigheff e+:sig",
"D*-:D0_Kpipi0_eff20 e+:sig",
"D*-:D0_Kpipipi_skimhigheff e+:sig",
"D*-:D0_Kpi_skimhigheff mu+:sig",
"D*-:D0_Kpipi0_eff20 mu+:sig",
"D*-:D0_Kpipipi_skimhigheff mu+:sig"]
B0_list = []
for chID, channel in enumerate(B0_channels):
ma.reconstructDecay("B0:Dstl_kpi_kpipi0_kpipipi" + str(chID) + " -> " + channel, "", chID, path=path)
B0_list.append("B0:Dstl_kpi_kpipi0_kpipipi" + str(chID))
self.SkimLists = B0_list