####EvtInclusiveDecay.userDecayFile = "DstarP2D0PiP_D02KpiPlusAnti.DEC"
EvtInclusiveDecay.pdtFile = "inclusive.pdt"
EvtInclusiveDecay.decayFile = "inclusive.dec"
# Build jets for the filter
try:
from JetRec.JetGetters import *
a4alg = make_StandardJetGetter('AntiKt', 0.4, 'Truth').jetAlgorithmHandle()
except Exception, e:
pass
#--------------------------------------------------------------
# Filter on the B and charm hadrons
#--------------------------------------------------------------
from GeneratorFilters.GeneratorFiltersConf import HeavyFlavorHadronFilter
topAlg += HeavyFlavorHadronFilter()
topAlg.HeavyFlavorHadronFilter.OutputLevel = INFO
topAlg.HeavyFlavorHadronFilter.RequestCharm = TRUE
topAlg.HeavyFlavorHadronFilter.RequestBottom = TRUE
topAlg.HeavyFlavorHadronFilter.Request_cQuark = TRUE
topAlg.HeavyFlavorHadronFilter.Request_bQuark = TRUE
topAlg.HeavyFlavorHadronFilter.RequestSpecificPDGID = FALSE
topAlg.HeavyFlavorHadronFilter.RequireTruthJet = TRUE
topAlg.HeavyFlavorHadronFilter.TruthContainerName = "AntiKt4TruthJets"
topAlg.HeavyFlavorHadronFilter.JetPtMin = 8.0 * GeV
topAlg.HeavyFlavorHadronFilter.JetEtaMax = 5.0
topAlg.HeavyFlavorHadronFilter.BottomPtMin = 3.0 * GeV
topAlg.HeavyFlavorHadronFilter.CharmPtMin = 3.0 * GeV
topAlg.HeavyFlavorHadronFilter.BottomEtaMax = 2.7
topAlg.HeavyFlavorHadronFilter.CharmEtaMax = 2.7
topAlg.HeavyFlavorHadronFilter.DeltaRFromTruth = 0.3
## Example configuration for HeavyFlavorHadronFilter setting up defaults
include("MC15JobOptions/AntiKt4TruthJets_pileup.py")
include("MC15JobOptions/AntiKt6TruthJets_pileup.py")
include("MC15JobOptions/JetFilter_Fragment.py")
if not hasattr(filtSeq, "HeavyFlavorHadronFilter"):
from GeneratorFilters.GeneratorFiltersConf import HeavyFlavorHadronFilter
filtSeq += HeavyFlavorHadronFilter()
pass
# Setup defaults as example
"""
HeavyFlavorHadronFilter = filtSeq.HeavyFlavorHadronFilter
HeavyFlavorHadronFilter.CharmPtMin = 2.0 * Gaudi::Units::GeV
HeavyFlavorHadronFilter.BottomPtMin = 5.0 * Gaudi::Units::GeV
HeavyFlavorHadronFilter.CharmEtaMax = 3.0
HeavyFlavorHadronFilter.BottomEtaMax = 3.0
HeavyFlavorHadronFilter.cPtMin = 0.0 * Gaudi::Units::GeV
HeavyFlavorHadronFilter.bPtMin = 0.0 * Gaudi::Units::GeV
HeavyFlavorHadronFilter.cEtaMax = 5.0
HeavyFlavorHadronFilter.bEtaMax = 5.0
HeavyFlavorHadronFilter.JetPtMin = 20.0 * Gaudi::Units::GeV
HeavyFlavorHadronFilter.JetEtaMax = 2.5
HeavyFlavorHadronFilter.PDGPtMin = 2.0 * Gaudi::Units::GeV
HeavyFlavorHadronFilter.PDGEtaMax = 2.5
HeavyFlavorHadronFilter.PDGID = 0
HeavyFlavorHadronFilter.PDGAntiParticleToo = True
HeavyFlavorHadronFilter.RequestCharm = True
HeavyFlavorHadronFilter.RequestBottom = True
HeavyFlavorHadronFilter.Request_cQuark = True
##
## Make a HeavyFlavorHadronFilter and steer it to select at least
## one AntiKt jet with pT>15 GeV, |eta|<3 containing stable
## C Hadrons with pT>4 GeV matched within DR<0.5
##
if not hasattr(filtSeq, "HeavyFlavorCHadronPt4Eta4_JetFilter"):
from GeneratorFilters.GeneratorFiltersConf import HeavyFlavorHadronFilter
HeavyFlavorCHadronPt4Eta4_JetFilter = HeavyFlavorHadronFilter(
name="HeavyFlavorCHadronPt4Eta4_JetFilter")
# Anti-kt truth jets with pt>7 GeV and R=0.4
include("MC15JobOptions/AntiKt4TruthJets.py")
#
HeavyFlavorCHadronPt4Eta4_JetFilter.RequestBottom = False
HeavyFlavorCHadronPt4Eta4_JetFilter.RequestCharm = True
HeavyFlavorCHadronPt4Eta4_JetFilter.Request_cQuark = False
HeavyFlavorCHadronPt4Eta4_JetFilter.Request_bQuark = False
HeavyFlavorCHadronPt4Eta4_JetFilter.RequestSpecificPDGID = False
HeavyFlavorCHadronPt4Eta4_JetFilter.RequireTruthJet = True
HeavyFlavorCHadronPt4Eta4_JetFilter.CharmPtMin = 4 * GeV
HeavyFlavorCHadronPt4Eta4_JetFilter.CharmEtaMax = 4.0
#
HeavyFlavorCHadronPt4Eta4_JetFilter.DeltaRFromTruth = 0.5
HeavyFlavorCHadronPt4Eta4_JetFilter.JetPtMin = 15.0 * GeV
HeavyFlavorCHadronPt4Eta4_JetFilter.JetEtaMax = 3.0
HeavyFlavorCHadronPt4Eta4_JetFilter.TruthContainerName = "AntiKt4TruthJets"
