## ## 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"