Source code for skim.WGs.semileptonic

#!/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.                  #
##########################################################################

"""(Semi-)Leptonic Working Group Skim list building functions for semi-leptonic analyses.
"""

import basf2 as b2
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 skim 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
# from ROOT import Belle2

__liaison__ = "Cameron Harris <cameron.harris@adelaide.edu.au>, Tommy Martinov <tommy.martinov@desy.de>"
_VALIDATION_SAMPLE = "mdst16.root"


[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" ApplyHLTHadronCut = False validation_sample = _VALIDATION_SAMPLE
[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.cutAndCopyList("pi+:PRSL_eventshape", "pi+:all", 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) return ["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.copyLists('B0:all', self.SkimLists, path=path) ma.buildRestOfEvent('B0:all', path=path) ma.appendROEMask('B0:all', 'basic', 'pt>0.05 and dr<2 and -4.0<dz<4.0', 'E>0.05', 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' create_validation_histograms( rootfile=histogramFilename, particlelist='B0:all', variables_1d=[ ('Mbc', 100, 4.0, 5.3, 'Mbc', __liaison__, '', ''), ('d0_p', 100, 0, 0.3, 'Signal-side pion momentum', __liaison__, '', ''), ('d1_p', 100, 1, 4, 'Signal-side lepton momentum', __liaison__, '', ''), ('MissM2', 100, -5, 5, 'Missing mass squared', __liaison__, '', '') ], variables_2d=[('deltaE', 100, -5, 5, 'Mbc', 100, 4.0, 5.3, 'Mbc vs deltaE', __liaison__, '', '')], 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" ApplyHLTHadronCut = False validation_sample = _VALIDATION_SAMPLE
[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_May2020", 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): path = self.skim_event_cuts("nTracks > 4", path=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) bplusList.append(f"B+:SLUntagged_{chID}") b0List = [] for chID, channel in enumerate(B0Channels): ma.reconstructDecay(f"B0:SLUntagged_{chID} -> {channel}", Bcuts, chID, path=path) b0List.append(f"B0:SLUntagged_{chID}") return 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.copyLists('B+:all', [lst for lst in self.SkimLists if "B+" in lst], path=path) ma.buildRestOfEvent('B+:all', path=path) ma.appendROEMask('B+:all', 'basic', 'pt>0.05 and -2<dr<2 and -4.0<dz<4.0', 'E>0.05', path=path) vm.addAlias('d1_p', 'daughter(1,p)') vm.addAlias('MissM2', 'weMissM2(basic,0)') histogramFilename = f'{self}_Validation.root' create_validation_histograms( rootfile=histogramFilename, particlelist='B+:all', variables_1d=[ ('cosThetaBetweenParticleAndNominalB', 100, -3.0, 3.0, 'cosThetaBY', __liaison__, '', ''), ('Mbc', 100, 5.24, 5.3, 'Mbc', __liaison__, '', ''), ('d1_p', 100, 1, 4, 'Signal-side lepton momentum', __liaison__, '', ''), ('MissM2', 100, -5, 5, 'Missing mass squared', __liaison__, '', '') ], variables_2d=[('deltaE', 100, -1, 1, 'Mbc', 100, 5.2, 5.3, 'Mbc vs deltaE', __liaison__, '', '')], path=path)
[docs] @fancy_skim_header class B0toDstarl_Kpi_Kpipi0_Kpipipi(BaseSkim): """ Cuts applied: * ``SkimHighEff tracks thetaInCDCAcceptance AND dr < 2 AND abs(dz) < 5 AND PID>=0.01`` * ``slowPi tracks thetaInCDCAcceptance AND dr < 2 AND abs(dz) < 5 AND useCMSFrame(p) < 0.4`` * :math:`2.5 > p_{\\ell} > 1.1\\,\\text{GeV}` * ``lepton with abs(dr) < 0.5 AND abs(dz) < 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_May2020", 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(dr) < 0.5 and abs(dz) < 2 and thetaInCDCAcceptance and electronID >= 0.95 and 1.1 < useCMSFrame(p) < 2.5 ', path=path) ma.cutAndCopyList( 'mu+:sig', 'mu+:all', 'abs(dr) < 0.5 and abs(dz) < 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)) return B0_list
# new SLskim
[docs] @fancy_skim_header class BtoDl_and_ROE_e_or_mu_or_lowmult(BaseSkim): """ Reconstructed decay modes: * :math:`B^- \\to D^0 e^-`, :math:`B^- \\to D^0 \\mu^-` * anti-:math:`B^0 \\to D^+ e^-`, anti-:math:`B^0 \\to D^+ \\mu^-` * anti-:math:`B^0 \\to D^{*+}(D^0 \\pi^+) e^-`, anti-:math:`B^0 \\to D^{*+}(D^0 \\pi^+) \\mu^-` Additional ROE cuts: * The logical OR of the following: * identified :math:`e^{\\pm}` with :math:`p(CM) < 3.0` GeV * identified :math:`\\mu^{\\pm}` with :math:`p(CM) < 3.0` GeV * identified :math:`\\gamma` with :math:`E(CM) > 1.4` GeV * Two or fewer charged tracks * :math:`E_{ECL} < 2.0` GeV Cuts on electrons: * :math:`\\text{electronID} > 0.3` * :math:`p_t > 0.3\\,\\text{GeV}` in lab frame, :math:`p > 0.5\\,\\text{GeV}` in lab frame * :math:`dr < 0.5`, :math:`|dz| < 2` * :math:`\\text{thetaInCDCAcceptance}` Cuts on muons: * :math:`\\text{muonID} > 0.9` * :math:`p_t > 0.4\\,\\text{GeV}` in lab frame, :math:`p > 0.7\\,\\text{GeV}` in lab frame * :math:`dr < 0.5`, :math:`|dz| < 2` Charged mask for ROE: * :math:`p_t>0.05\\,\\text{GeV}` * :math:`dr < 5\\,\\text{cm}`, :math:`|dz| < 10\\,\\text{cm}` ECL cluster mask for ROE: * :math:`\\text{clusterNHits}>1.5`, :math:`\\theta` in CDC acceptance * :math:`E>0.080,\\text{GeV}` """ __authors__ = ["Bob Kowalewski"] __description__ = "Skim for semileptonic tags with an ROE electron, muon, photon or a low-multiplicity signal decay" __contact__ = __liaison__ __category__ = "physics, semileptonic" ApplyHLTHadronCut = True TestSampleProcess = "charged" validation_sample = _VALIDATION_SAMPLE
[docs] def load_standard_lists(self, path): stdKshorts(path=path)
[docs] def build_lists(self, path): eIDCut = "electronID > 0.3" muIDCut = "muonID > 0.9" ePCut = "p > 0.5" ePtCut = "pt > 0.3" muPCut = "p > 0.7" muPtCut = "pt > 0.4" lepTrkCuts = "dr < 0.5 and abs(dz) < 2" hadTrkCuts = "dr < 2.0 and abs(dz) < 4" gammaCuts = "[clusterNHits>1.5] and thetaInCDCAcceptance" gammaECuts = "[[clusterReg==1 and E>0.025] or [clusterReg==2 and E>0.025] or [clusterReg==3 and E>0.040]]" cleanMask = ('cleanMask', 'pt>0.05 and dr < 5 and abs(dz) < 10', f'{gammaCuts} and E>0.080 and minC2TDist>20.0 and abs(clusterTiming)<200') gammaSignalECut = "1.4" BSLRecoCut = "InvM > 2.5 and cosBY > -3.0 and cosBY < 1.5" BSLSkimCut = f"e_ROE_pCM < 3.0 or mu_ROE_pCM < 3.0 or gamma_ROE_ECM > {gammaSignalECut}" +\ " or nROE_Ch < 2.5 or E_extra_ROE < 1.6" vm.addAlias('pCM', 'useCMSFrame(p)') vm.addAlias('ECM', 'useCMSFrame(E)') vm.addAlias('cosBY', 'cosThetaBetweenParticleAndNominalB') # electrons and muons ma.fillParticleList("e-:sig", f"{lepTrkCuts} and thetaInCDCAcceptance and {ePtCut} and {ePCut}", path=path) ma.applyCuts('e-:sig', eIDCut, path=path) ma.fillParticleList("mu-:sig", f"{muIDCut} and {lepTrkCuts} and {muPtCut} and {muPCut}", path=path) total_leptons = "formula( nParticlesInList(e-:sig) + nParticlesInList(mu-:sig) )" # don't waste time if there are no leptons path = self.skim_event_cuts(f"{total_leptons} > 0", path=path) # charged K and pi candidates ma.fillParticleList("K-:BtoDl", f"binaryPID(321,211) > 0.1 and {hadTrkCuts} and thetaInCDCAcceptance and pt>0.1", path=path) ma.fillParticleList( "pi+:BtoDl", f"binaryPID(211,321) > 0.1 and {hadTrkCuts} and thetaInCDCAcceptance and pt>0.1", path=path) ma.fillParticleList("pi+:slow", f"{hadTrkCuts} and thetaInCDCAcceptance and p>0.05", path=path) # gamma and pi0 candidates ma.fillParticleList("gamma:eff50", f'{gammaCuts} and {gammaECuts}', path=path) ma.reconstructDecay("pi0:eff50 -> gamma:eff50 gamma:eff50", '0.114<InvM<0.150', path=path) ma.cutAndCopyList("gamma:e90", "gamma:eff50", "E>0.090", path=path) ma.reconstructDecay("pi0:e90 -> gamma:e90 gamma:e90", '0.118<InvM<0.148', path=path) ma.cutAndCopyList('K_S0:good', 'K_S0:merged', cut='significanceOfDistance>2.5', path=path) # D reconstruction - only clean modes ma.reconstructDecay("D0:Kpi_BtoDl -> K-:BtoDl pi+:BtoDl", cut="abs(dM) < 0.025", dmID=1, path=path) ma.rankByLowest("D0:Kpi_BtoDl", variable="abs(dM)", numBest=15, path=path) ma.reconstructDecay("D0:KK -> K-:BtoDl K+:BtoDl", cut="abs(dM) < 0.020", dmID=2, path=path) ma.rankByLowest("D0:KK", variable="abs(dM)", numBest=15, path=path) ma.reconstructDecay("D0:K3pi -> K-:BtoDl pi+:BtoDl pi-:BtoDl pi+:BtoDl", cut="abs(dM) < 0.020", dmID=3, path=path) ma.rankByLowest("D0:K3pi", variable="abs(dM)", numBest=15, path=path) ma.reconstructDecay("D0:Kpi0pi -> K-:BtoDl pi0:e90 pi+:BtoDl", cut="abs(dM)-0.005 < 0.050", dmID=4, path=path) ma.rankByLowest("D0:Kpi0pi", variable="abs(dM)", numBest=15, path=path) ma.reconstructDecay("D0:piKspi -> pi+:BtoDl K_S0:good pi-:BtoDl", cut="abs(dM) < 0.020", dmID=5, path=path) ma.rankByLowest("D0:piKspi", variable="abs(dM)", numBest=15, path=path) ma.reconstructDecay("D0:Kspi0_BtoDl -> K_S0:good pi0:e90", cut="abs(dM)-0.005 < 0.050", dmID=6, path=path) ma.rankByLowest("D0:Kspi0_BtoDl", variable="abs(dM)", numBest=15, path=path) ma.copyLists("D0:sig", ['D0:Kpi_BtoDl', 'D0:KK', 'D0:K3pi', 'D0:Kpi0pi', 'D0:piKspi', 'D0:Kspi0_BtoDl'], path=path) ma.reconstructDecay("D+:Kpipi_BtoDl -> K-:BtoDl pi+:BtoDl pi+:BtoDl", cut="abs(dM) < 0.020", dmID=1, path=path) ma.rankByLowest("D+:Kpipi_BtoDl", variable="abs(dM)", numBest=15, path=path) ma.reconstructDecay("D+:Kspi_BtoDl -> K_S0:good pi+:BtoDl", cut="abs(dM) < 0.020", dmID=2, path=path) ma.rankByLowest("D+:Kspi_BtoDl", variable="abs(dM)", numBest=15, path=path) ma.reconstructDecay("D+:KKpi -> K-:BtoDl K+:BtoDl pi+:BtoDl", cut="abs(dM) < 0.015", dmID=3, path=path) ma.rankByLowest("D+:KKpi", variable="abs(dM)", numBest=15, path=path) ma.copyLists("D+:sig", ['D+:Kpipi_BtoDl', 'D+:Kspi_BtoDl', 'D+:KKpi'], path=path) # No explicit reconstruction of D*0 to D0 pi0/gamma or D*+ to D0 pi+ (included in feeddown) ma.reconstructDecay("D*+:D0 -> pi+:slow D0:sig", cut="abs(dQ) < 0.0035", dmID=11, path=path) ma.rankByLowest("D*+:D0", variable="abs(dQ)", numBest=15, path=path) ma.copyLists("D*+:sig", ['D*+:D0'], path=path) # B reconstruction. Use wide cut on cosThetaBY. ma.reconstructDecay('anti-B0:Dpe -> e-:sig D+:sig ?nu ', BSLRecoCut, dmID=200, path=path) ma.reconstructDecay('anti-B0:Dpmu -> mu-:sig D+:sig ?nu ', BSLRecoCut, dmID=201, path=path) ma.copyLists('anti-B0:Dpl', ['anti-B0:Dpe', 'anti-B0:Dpmu'], path=path) ma.rankByHighest('anti-B0:Dpl', variable='InvM', numBest=9, path=path) ma.reconstructDecay('anti-B0:Dspe -> e-:sig D*+:sig ?nu ', BSLRecoCut, dmID=100, path=path) ma.reconstructDecay('anti-B0:Dspmu -> mu-:sig D*+:sig ?nu ', BSLRecoCut, dmID=101, path=path) ma.copyLists('anti-B0:Dspl', ['anti-B0:Dspe', 'anti-B0:Dspmu'], path=path) ma.rankByHighest('anti-B0:Dspl', variable='InvM', numBest=9, path=path) ma.reconstructDecay('B-:D0e -> e-:sig D0:sig ?nu ', BSLRecoCut, dmID=400, path=path) ma.reconstructDecay('B-:D0mu -> mu-:sig D0:sig ?nu ', BSLRecoCut, dmID=401, path=path) ma.copyLists('B-:D0l', ['B-:D0e', 'B-:D0mu'], path=path) ma.rankByHighest('B-:D0l', variable='InvM', numBest=9, path=path) # main SL B lists ma.copyLists('B-:SL', ['B-:D0l'], path=path) ma.copyLists('anti-B0:SL', ['anti-B0:Dspl', 'anti-B0:Dpl'], path=path) # Build ROE and look for e-, mu-, charged multiplicity and E_ECL for Btype in ['B-:SL', 'anti-B0:SL']: roe_path = b2.create_path() deadEndPath = b2.create_path() # Execute the filter module: ma.buildRestOfEvent(Btype, fillWithMostLikely=True, path=path) ma.appendROEMasks(Btype, [cleanMask], path=path) ma.signalSideParticleFilter(Btype, '', roe_path, deadEndPath) ma.fillParticleList( "e-:roeB", f"isInRestOfEvent == 1 and {eIDCut} and {lepTrkCuts} and thetaInCDCAcceptance and {ePtCut} and {ePCut}", path=roe_path) ma.fillParticleList( "mu-:roeB", f"isInRestOfEvent == 1 and {muIDCut} and {lepTrkCuts} and {muPtCut} and {muPCut}", path=roe_path) ma.fillParticleList("gamma:roeB", f'{gammaCuts} and ECM > {gammaSignalECut}', path=roe_path) ma.rankByHighest('e-:roeB', 'pCM', numBest=1, path=roe_path) ma.rankByHighest('mu-:roeB', 'pCM', numBest=1, path=roe_path) ma.rankByHighest('gamma:roeB', 'ECM', numBest=1, path=roe_path) ma.variableToSignalSideExtraInfo('e-:roeB', {'pCM': 'e_ROEB_pCM'}, path=roe_path) ma.variableToSignalSideExtraInfo('mu-:roeB', {'pCM': 'mu_ROEB_pCM'}, path=roe_path) ma.variableToSignalSideExtraInfo('gamma:roeB', {'ECM': 'gamma_ROEB_ECM'}, path=roe_path) path.for_each('RestOfEvent', 'RestOfEvents', path=roe_path) vm.addAlias('e_ROE_pCM', 'extraInfo(e_ROEB_pCM)') vm.addAlias('mu_ROE_pCM', 'extraInfo(mu_ROEB_pCM)') vm.addAlias('gamma_ROE_ECM', 'extraInfo(gamma_ROEB_ECM)') vm.addAlias('nROE_Ch', 'nROE_Charged(cleanMask)') vm.addAlias('E_extra_ROE', 'useCMSFrame(roeEextra(cleanMask))') # Only keep events whose ROE has an extra e or mu, or Ntrk<3, or E_ECL<2.0 GeV ma.applyCuts('B-:SL', f'passesCut({BSLSkimCut})', path=path) ma.applyCuts('anti-B0:SL', f'passesCut({BSLSkimCut})', path=path) b_sl_list = ['B-:SL', 'anti-B0:SL'] return b_sl_list
[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 vm.addAlias('d0_p', 'daughter(0, p)') vm.addAlias('d1_p', 'daughter(1, p)') histogramFilename = f'{self}_Validation.root' create_validation_histograms( rootfile=histogramFilename, particlelist='B-:SL', variables_1d=[ ('cosBY', 90, -3.0, 1.5, 'cosThetaBY', __liaison__, 'Cosine of angle between the reconstructed B and the nominal B', '', 'cos#theta_{BY}', 'Candidates'), ('InvM', 100, 2.5, 5.3, 'M(D(*)+l)', __liaison__, 'Invariant mass distribution of D(*)l system', '', 'm(D(*)l) [GeV]', 'Candidates / 28 MeV'), ('useCMSFrame(d0_p)', 100, 0, 3.0, 'Tag lepton momentum in CMS', __liaison__, 'Momentum of tag-side lepton in CMS', '', 'p(l) [GeV]', 'Candidates / 30 MeV'), ('useCMSFrame(d1_p)', 100, 0, 3.0, 'Tag D(*) momentum in CMS', __liaison__, 'Momentum of tag-side D(*) meson in CMS', '', 'p(D(*)) [GeV]', 'Candidates / 30 MeV'), ('e_ROE_pCM', 100, 0, 3.0, 'ROE e momentum in CMS', __liaison__, 'Momentum of highest-momentum ROE electron in CMS', '', 'p(e) [GeV]', 'Candidates / 30 MeV'), ('mu_ROE_pCM', 100, 0, 3.0, 'ROE mu momentum in CMS', __liaison__, 'Momentum of highest-momentum ROE muon in CMS', '', 'p($\\mu$) [GeV]', 'Candidates / 30 MeV'), ('gamma_ROE_ECM', 100, 0, 5.0, 'ROE gamma energy in CMS', __liaison__, 'Energy of most-energetic ROE photon', '', 'E($\\gamma$) [GeV]', 'Candidates / 50 MeV'), ('nROE_Ch', 20, 0, 20.0, 'N(trk) in ROE', __liaison__, 'Number of tracks in ROE', '', 'tracks', 'Candidates'), ('E_extra_ROE', 100, 0, 2.5, 'E_ECL (CMS) in ROE', __liaison__, 'Extra energy in ROE', '', 'E_{extra} [GeV]', 'Candidates / 25 MeV'), ], path=path)