##
## Make a HeavyFlavorHadronFilter and steer it to select C Hadrons
## with any pT within |eta| < 4
##
if not hasattr(filtSeq, "HeavyFlavorCHadronFilter"):
from GeneratorFilters.GeneratorFiltersConf import HeavyFlavorHadronFilter
HeavyFlavorCHadronFilter = HeavyFlavorHadronFilter(name="HeavyFlavorCHadronFilter")
## Default cut params
HeavyFlavorCHadronFilter.RequestBottom=False
HeavyFlavorCHadronFilter.RequestCharm=True
HeavyFlavorCHadronFilter.Request_cQuark=False
HeavyFlavorCHadronFilter.Request_bQuark=False
HeavyFlavorCHadronFilter.RequestSpecificPDGID=False
HeavyFlavorCHadronFilter.RequireTruthJet=False
HeavyFlavorCHadronFilter.CharmPtMin=0*GeV
HeavyFlavorCHadronFilter.CharmEtaMax=4.0
## Filter on presence of Dstar hadrond
if not hasattr(filtSeq, "HeavyFlavorHadronFilter"):
from GeneratorFilters.GeneratorFiltersConf import HeavyFlavorHadronFilter
filtSeq += HeavyFlavorHadronFilter()
pass
HeavyFlavorHadronFilter = filtSeq.HeavyFlavorHadronFilter
HeavyFlavorHadronFilter.RequestCharm = FALSE
HeavyFlavorHadronFilter.RequestBottom = FALSE
HeavyFlavorHadronFilter.Request_cQuark = FALSE
HeavyFlavorHadronFilter.Request_bQuark = FALSE
HeavyFlavorHadronFilter.RequestSpecificPDGID = TRUE
HeavyFlavorHadronFilter.PDGID = -413
HeavyFlavorHadronFilter.PDGAntiParticleToo = FALSE
HeavyFlavorHadronFilter.RequireTruthJet = FALSE
HeavyFlavorHadronFilter.PDGEtaMax = 3.0
HeavyFlavorHadronFilter.PDGPtMin = 4.0 * GeV
HeavyFlavorHadronFilter.RequireTruthJet = False
##
## Make a HeavyFlavorHadronFilter and steer it to select B Hadrons
## with any pT within |eta| < 4
##
if not hasattr(filtSeq, "HeavyFlavorBHadronFilter"):
from GeneratorFilters.GeneratorFiltersConf import HeavyFlavorHadronFilter
HeavyFlavorBHadronFilter = HeavyFlavorHadronFilter(
name="HeavyFlavorBHadronFilter")
## Default cut params
HeavyFlavorBHadronFilter.RequestBottom = True
HeavyFlavorBHadronFilter.RequestCharm = False
HeavyFlavorBHadronFilter.Request_cQuark = False
HeavyFlavorBHadronFilter.Request_bQuark = False
HeavyFlavorBHadronFilter.RequestSpecificPDGID = False
HeavyFlavorBHadronFilter.RequireTruthJet = False
HeavyFlavorBHadronFilter.BottomPtMin = 0 * GeV
HeavyFlavorBHadronFilter.BottomEtaMax = 4.0
##
## 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"
####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
##
## 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"
