from D3PDMakerConfig.D3PDMakerFlags import D3PDMakerFlags
from AthenaCommon.AlgSequence import AlgSequence, AthSequencer
# create electron and muon selectors
sequencer = AthSequencer("HSG5ZBB_Sequencer", StopOverride=False)
from D2PDMaker.D2PDMakerConf import D2PDJetSelector
boostedJetCollectionName = 'HSG5ZBB_boosted_LooseBoostedJetLinkCollection'
sequencer += D2PDJetSelector("HSG5ZBB_BoostedJetFilter",
inputCollection='CamKt12LCTopoJets',
outputLinkCollection=boostedJetCollectionName,
minNumberPassed=1,
ptMin=150.0 * Units.GeV)
HSG5D3PD_Stream.RequireAlgs.append("HSG5ZBB_BoostedJetFilter")
from AssociationComps.AssociationCompsConf import DeltaRAssociationTool
ToolSvc += DeltaRAssociationTool(
"HSG5ZBB_subjetBoostedJetDeltaRAssociationTool",
OutputLevel=INFO,
inputAssociateToCollection=boostedJetCollectionName,
deltaRMax=1.3,
writeUserData=False)
if not 'IS_SIMULATION' in inputFileSummary['evt_type']:
sequencer += D2PDJetSelector(
"HSG5ZBB_btaggedSubjetFilter",
OutputLevel=DEBUG,
inputCollection='CamKt12LCTopoSplitFiltSubjetsminSplitR0Jets',
outputLinkCollection='HSG5ZBB_btaggedSubjetLinkCollection',
associationToolList=[
ToolSvc.HSG5ZBB_subjetBoostedJetDeltaRAssociationTool
include.block( "PhysicsD3PDMaker/JetSelector.py" )
from D2PDMaker.D2PDMakerConf import D2PDJetSelector
from AthenaCommon.AlgSequence import AlgSequence
preseq = AlgSequence (D3PDMakerFlags.PreD3PDAlgSeqName())
preseq += D2PDJetSelector( "SMWZ_HighPtJetFilter",
inputCollection = 'AntiKt4TopoEMJets',
outputCollection = 'HighPtJets',
ptMin = 20.0*Units.GeV )
preseq += D2PDJetSelector( "SMWZ_HighPtJetEM6Filter",
inputCollection = 'AntiKt6TopoEMJets',
outputCollection = 'HighPtJetsEM6',
ptMin = 20.0*Units.GeV )
preseq += D2PDJetSelector( "SMWZ_HighPtJetLC4Filter",
inputCollection = 'AntiKt4LCTopoJets',
outputCollection = 'HighPtJetsLC4',
ptMin = 20.0*Units.GeV )
preseq += D2PDJetSelector( "SMWZ_HighPtJetLC6Filter",
inputCollection = 'AntiKt4LCTopoJets',
outputCollection = 'HighPtJetsLC6',
ptMin = 20.0*Units.GeV )
jetAssociationToolList += [getattr(ToolSvc, filterName)]
jetOutputAssociationContainerList += [
jetCollName + "_Match_" + offlineName
]
pass
from D2PDMaker.D2PDMakerConf import D2PDJetSelector
# Create an instance of this filter and add it to the AlgSequence
# Also, configure the filter to your needs. This are just some examples
filterName = jetCollName + "SelectorInHighMassDiJetStream"
topSequence += D2PDJetSelector(
filterName,
OutputLevel=INFO,
inputCollection=jetCollName,
outputLinkCollection=jetCollName +
'LinkCollectionInHighMassDiJetStream',
associationToolList=jetAssociationToolList,
outputAssociationContainerList=jetOutputAssociationContainerList,
minNumberPassed=2,
absEtaMax=2.8,
ptMin=30.0 * Units.GeV)
algsToBookkeep.append(filterName)
filterName = "HighPt" + jetCollName + "SelectorInHighMassDiJetStream"
topSequence += D2PDJetSelector(filterName,
OutputLevel=INFO,
inputCollection=jetCollName,
outputLinkCollection='HighP' + jetCollName +
'LinkCollectionInHighMassDiJetStream',
minNumberPassed=1,
absEtaMax=2.8,
from D3PDMakerConfig.D3PDMakerFlags import D3PDMakerFlags
from AthenaCommon.AlgSequence import AlgSequence,AthSequencer
sequencer = AthSequencer("HSG5ZHMET_Sequencer",
StopOverride=False)
# boosted jet selector
from D2PDMaker.D2PDMakerConf import D2PDJetSelector
sequencer += D2PDJetSelector( "HSG5ZHMET_BoostedJetFilter",
inputCollection = 'CamKt12LCTopoJets',
outputLinkCollection = 'HSG5ZHMET_boosted_LooseBoostedJetLinkCollection',
minNumberPassed = 1,
ptMin = 20.0*Units.GeV)
HSG5D3PD_Stream.RequireAlgs.append("HSG5ZHMET_BoostedJetFilter")
# MET selector
from D2PDMaker.D2PDMakerConf import D2PDMissingEtSelector
sequencer += D2PDMissingEtSelector( "HSG5ZHMET_MissingEtFilter",
inputCollection = 'MET_RefFinal',
outputCollection = 'HSG5ZHMET_NeutrinoRefFinal',
minNumberPassed = 1,
missingEtMin = 40.0*Units.GeV)
#missingEtSignificanceMin = 5.0)
HSG5D3PD_Stream.RequireAlgs.append("HSG5ZHMET_MissingEtFilter")
if False:
# (for private production only) insert in beginning of PreD3PDSequencer
mainSequencer = AlgSequence(D3PDMakerFlags.PreD3PDAlgSeqName(),
#photonAuthor = egammaParameters.AuthorPhoton,
#requested by Jovan
etMin = 20.0*Units.GeV,
clusterEMCaloLayer = 2,
clusterEMCaloLayerEtaMin = -2.5,
clusterEMCaloLayerEtaMax = 2.5
)
# Jet filters
from D2PDMaker.D2PDMakerConf import D2PDJetSelector
JetSelectorSUSY = D2PDJetSelector( "JetSelectorSUSY",
inputCollection = D2PDFlags.InputSUSYJetCollection() ,
outputLinkCollection = 'AntiKt4H1TopoJets_SelectedJetLinks',
minNumberPassed = 1,
useJetSignalState = 0,
jetScaledEtMin = 80.0*Units.GeV,
jetSignalStateAbsEtaMax = 3.0
)
DiJetSelectorSUSY = D2PDJetSelector( "DiJetSelectorSUSY",
inputCollection = D2PDFlags.InputSUSYJetCollection() ,
outputLinkCollection = 'AntiKt4H1TopoJets_SelectedDijetLinks',
minNumberPassed = 2,
useJetSignalState = 0,
jetScaledEtMin = 30.0*Units.GeV,
jetSignalStateAbsEtaMax = 3.0
)
#FourJet not used yet, but likely to be soon when dijet filter tightens
FourJetSelectorSUSY = D2PDJetSelector( "FourJetSelectorSUSY",
from UserDataUtils.UserDataUtilsConf import UserDataCalcJetQualityTool
ToolSvc += UserDataCalcJetQualityTool( "UserDataCalcJetQualityToolInWenuStream",
OutputLevel = INFO,
userDataSvc = "UserDataInWenuStream",
userDataPrefix = "jet_",
userDataJetQualityName = "isGoodUglyBad"
)
# Now, use the JetSelector and use the UserDataCalcJetQualityTool with it
from D2PDMaker.D2PDMakerConf import D2PDJetSelector
topSequence += D2PDJetSelector( "JetMetCleaningSelectorInWenuStream",
OutputLevel = INFO,
userDataSvc = "UserDataInWenuStream",
inputCollection = 'AntiKt4TopoJets',
outputLinkCollection = 'WenuStreamEMScalePt10BadJetCollection',
userDataCalcToolList = [ ToolSvc.UserDataCalcJetQualityToolInWenuStream ],
minNumberPassed = 1,
useJetSignalState = PyAthena.P4SignalState.JETEMSCALE,
jetSignalStatePtMin = 10.0*Units.GeV,
userDataNameList = [ "jet_isGoodUglyBad" ],
userDataMinCutList = [ 3.5 ] # good=1, ugly=2, bad=4 (combinations add up)
)
algsToBookkeep.append( "JetMetCleaningSelectorInWenuStream" )
#--------------------------------------------------------------------
# Set up the cumulative cut-flow for the standard pre-selection.
# This will use the Logical Filter Combiners to add individual cuts.
#--------------------------------------------------------------------
cmdString_GRL_PrimVtx_Trigger_JetMetClean = cmdString_GRL_PrimVtx_Trigger + " and not JetMetCleaningSelectorInWenuStream"
topSequence += LogicalFilterCombiner( "StreamWenuFilterCombiner_GRL_PrimVtx_Trigger_JetMetClean",
cmdstring = cmdString_GRL_PrimVtx_Trigger_JetMetClean )
algsToBookkeep.append( "StreamWenuFilterCombiner_GRL_PrimVtx_Trigger_JetMetClean" )
trigger="EF_g([2-9][0-9][0-9]|1[2-9][0-9]+).*") # 120 or higher
HSG5D3PD_Stream.RequireAlgs.append("HSG5GAMZ_triggerFilter")
else:
print "no trigger selection"
zbb_or_zllSequencer = AthSequencer("HSG5GAMZ_zbb_or_zllSequencer",
StopOverride=True)
zbbSequencer = AthSequencer("HSG5GAMZ_zbbSequencer", StopOverride=False)
boostedJetCollectionName = 'HSG5GAMZ_boosted_LooseBoostedJetLinkCollection'
from D2PDMaker.D2PDMakerConf import D2PDJetSelector
zbbSequencer += D2PDJetSelector("HSG5GAMZ_BoostedJetFilter",
OutputLevel=outputLevel,
inputCollection='CamKt12LCTopoJets',
outputLinkCollection=boostedJetCollectionName,
minNumberPassed=1,
ptMin=80.0 * Units.GeV)
from AssociationComps.AssociationCompsConf import DeltaRAssociationTool
ToolSvc += DeltaRAssociationTool(
"HSG5GAMZ_subjetBoostedJetDeltaRAssociationTool",
OutputLevel=outputLevel,
inputAssociateToCollection=boostedJetCollectionName,
deltaRMax=1.3,
writeUserData=False)
zbbSequencer += D2PDJetSelector(
"HSG5GAMZ_btaggedSubjetFilter",
OutputLevel=outputLevel,
inputCollection='CamKt12LCTopoSplitFiltSubjetsminSplitR0Jets',
# jet selector
from AssociationComps.AssociationCompsConf import DeltaRAssociationTool
ToolSvc += DeltaRAssociationTool(
"HSG5WHU_emDeltaRAssociationTool",
OutputLevel=INFO,
inputAssociateToCollection='HSG5WHU_LooseElectronLinkCollection',
deltaRMax=0.3,
writeUserData=False)
from D2PDMaker.D2PDMakerConf import D2PDJetSelector
sequencer += D2PDJetSelector(
"HSG5WHU_JetFilter",
inputCollection='AntiKt4TopoEMJets',
outputLinkCollection='HSG5WHU_JetLinkCollection',
associationToolList=[ToolSvc.HSG5WHU_emDeltaRAssociationTool],
outputAssociationContainerList=["HSG5WHU_jetsMatchedToElectrons"],
numberOfAssociationsMaxCutList=[0],
minNumberPassed=1,
ptMin=20.0 * Units.GeV,
etaMax=2.5)
HSG5D3PD_Stream.RequireAlgs.append("HSG5WHU_JetFilter")
if False:
# (for private production ony) insert in beginning of PreD3PDSequencer
mainSequencer = AlgSequence(D3PDMakerFlags.PreD3PDAlgSeqName(),
StopOverride=False)
if not hasattr(topSequence, D3PDMakerFlags.PreD3PDAlgSeqName()):
topSequence += mainSequencer