Skip to content
Snippets Groups Projects
Commit fe7907bc authored by Jonathan Bossio's avatar Jonathan Bossio Committed by Frank Winklmeier
Browse files

Simplify conf of offline MT jet trigger monitoring

parent 022a08ba
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Trigger/TrigMonitoring/TrigJetMonitoring/python/TrigJetMonitorAlgorithm.py
+ 64
− 98
View file @ fe7907bc
......@@ -35,82 +35,54 @@ Chain2L1JetCollDict = { # set L1 jet collection name for L1 jet chains
# HLT jet collections and chains to monitor
############################################
Chain2JetCollDict = dict() # set HLT jet collection for AT and legacy master HLT jet chains
JetCollections = dict() # List of HLT jet collections for AT and legacy master
TurnOnCurves = dict() # List reference chains and offline jet collections to be used for producing turn-on curves
JetColls2Match = dict()
Chains2Monitor = dict() # set HLT jet collection, reference chain and offline jet collection for turn-on curves, for AT and legacy master HLT jet chains
JetCollections = dict() # List of HLT jet collections for AT and legacy master (stating which should be matched and to which offline jet collection
# AthenaMT
JetCollections['MT'] = [
'HLT_AntiKt4EMTopoJets_subjesIS', # default small-R EM
'HLT_AntiKt10JetRCJets_subjesIS', # a10r
'HLT_AntiKt10LCTopoJets_subjes', # a10
'HLT_AntiKt10LCTopoTrimmedPtFrac4SmallR20Jets_jes', # a10t
'HLT_AntiKt4EMPFlowJets_subjesIS_ftf', # pflow w/o calo+track GSC
'HLT_AntiKt4EMPFlowJets_subjesgscIS_ftf', # pflow w/ calo+track GSC
'HLT_AntiKt4EMPFlowJets_subresjesgscIS_ftf', # pflow w/ residual + calo+track GSC
'HLT_AntiKt4EMPFlowJets_nojcalib_ftf', # pflow nojcalib
'HLT_AntiKt4EMPFlowCSSKJets_nojcalib_ftf', # pflow cssk nojcalib
]
Chain2JetCollDict['MT'] = {
'HLT_j420_L1J100' : 'HLT_AntiKt4EMTopoJets_subjesIS',
'HLT_j260_320eta490_L1J75_31ETA49' : 'HLT_AntiKt4EMTopoJets_subjesIS',
'HLT_5j70_0eta240_L14J20' : 'HLT_AntiKt4EMTopoJets_subjesIS',
'HLT_3j200_L1J100' : 'HLT_AntiKt4EMTopoJets_subjesIS',
'HLT_j460_a10r_L1J100' : 'HLT_AntiKt10JetRCJets_subjesIS',
'HLT_j460_a10_lcw_subjes_L1J100' : 'HLT_AntiKt10LCTopoJets_subjes',
'HLT_j460_a10t_lcw_jes_L1J100' : 'HLT_AntiKt10LCTopoTrimmedPtFrac4SmallR20Jets_jes',
'HLT_2j330_a10t_lcw_jes_35smcINF_L1J100' : 'HLT_AntiKt10LCTopoTrimmedPtFrac4SmallR20Jets_jes',
'HLT_j45_pf_ftf_L1J15' : 'HLT_AntiKt4EMPFlowJets_subjesIS_ftf',
'HLT_j45_pf_subjesgscIS_ftf_L1J15' : 'HLT_AntiKt4EMPFlowJets_subjesgscIS_ftf',
'HLT_j85_pf_ftf_L1J20' : 'HLT_AntiKt4EMPFlowJets_subjesIS_ftf',
'HLT_j45_pf_nojcalib_ftf_L1J15' : 'HLT_AntiKt4EMPFlowJets_nojcalib_ftf',
'HLT_j45_csskpf_nojcalib_ftf_L1J15' : 'HLT_AntiKt4EMPFlowCSSKJets_nojcalib_ftf',
'HLT_10j40_pf_subresjesgscIS_ftf_L14J15' : 'HLT_AntiKt4EMPFlowJets_subresjesgscIS_ftf',
}
TurnOnCurves['MT'] = { # ref chain, offline jet coll
'HLT_j420_L1J100' : ['HLT_j85_L1J20','AntiKt4EMTopoJets'],
'HLT_3j200_L1J100' : ['HLT_j85_L1J20','AntiKt4EMTopoJets'],
'HLT_j460_a10r_L1J100' : ['HLT_j85_L1J20','AntiKt4EMTopoJets'],
'HLT_j460_a10_lcw_subjes_L1J100' : ['HLT_j85_L1J20','AntiKt4EMTopoJets'],
'HLT_j460_a10t_lcw_jes_L1J100' : ['HLT_j85_L1J20','AntiKt4EMTopoJets'],
'HLT_2j330_a10t_lcw_jes_35smcINF_L1J100' : ['HLT_j85_L1J20','AntiKt4EMTopoJets'],
'HLT_j85_pf_ftf_L1J20' : ['HLT_j45_pf_ftf_L1J15','AntiKt4EMPFlowJets'],
JetCollections['MT'] = {
'HLT_AntiKt4EMTopoJets_subjesIS' : { 'MatchTo' : 'AntiKt4EMPFlowJets' }, # default small-R EM
'HLT_AntiKt10JetRCJets_subjesIS' : { 'MatchTo' : 'NONE' }, # a10r
'HLT_AntiKt10LCTopoJets_subjes' : { 'MatchTo' : 'NONE' }, # a10
'HLT_AntiKt10LCTopoTrimmedPtFrac4SmallR20Jets_jes' : { 'MatchTo' : 'NONE' }, # a10t
'HLT_AntiKt4EMPFlowJets_subjesIS_ftf' : { 'MatchTo' : 'AntiKt4EMPFlowJets' }, # pflow w/o calo+track GSC
'HLT_AntiKt4EMPFlowJets_subjesgscIS_ftf' : { 'MatchTo' : 'AntiKt4EMPFlowJets' }, # pflow w/ calo+track GSC
'HLT_AntiKt4EMPFlowJets_subresjesgscIS_ftf' : { 'MatchTo' : 'AntiKt4EMPFlowJets' }, # pflow w/ residual + calo+track GSC
'HLT_AntiKt4EMPFlowJets_nojcalib_ftf' : { 'MatchTo' : 'NONE' }, # pflow nojcalib
'HLT_AntiKt4EMPFlowCSSKJets_nojcalib_ftf' : { 'MatchTo' : 'NONE' }, # pflow cssk nojcalib
}
JetColls2Match['MT'] = {
'HLT_AntiKt4EMTopoJets_subjesIS' : 'AntiKt4EMPFlowJets',
'HLT_AntiKt4EMPFlowJets_subjesIS_ftf' : 'AntiKt4EMPFlowJets',
'HLT_AntiKt4EMPFlowJets_subjesgscIS_ftf' : 'AntiKt4EMPFlowJets',
'HLT_AntiKt4EMPFlowJets_subresjesgscIS_ftf' : 'AntiKt4EMPFlowJets',
Chains2Monitor['MT'] = {
'HLT_j420_L1J100' : { 'HLTColl' : 'HLT_AntiKt4EMTopoJets_subjesIS', 'RefChain' : 'HLT_j85_L1J20', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_j260_320eta490_L1J75_31ETA49' : { 'HLTColl' : 'HLT_AntiKt4EMTopoJets_subjesIS', 'RefChain' : 'NONE', 'OfflineColl' : 'NONE' },
'HLT_5j70_0eta240_L14J20' : { 'HLTColl' : 'HLT_AntiKt4EMTopoJets_subjesIS', 'RefChain' : 'NONE', 'OfflineColl' : 'NONE' },
'HLT_3j200_L1J100' : { 'HLTColl' : 'HLT_AntiKt4EMTopoJets_subjesIS', 'RefChain' : 'HLT_j85_L1J20', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_j460_a10r_L1J100' : { 'HLTColl' : 'HLT_AntiKt10JetRCJets_subjesIS', 'RefChain' : 'HLT_j85_L1J20', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_j460_a10_lcw_subjes_L1J100' : { 'HLTColl' : 'HLT_AntiKt10LCTopoJets_subjes', 'RefChain' : 'HLT_j85_L1J20', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_j460_a10t_lcw_jes_L1J100' : { 'HLTColl' : 'HLT_AntiKt10LCTopoTrimmedPtFrac4SmallR20Jets_jes', 'RefChain' : 'HLT_j85_L1J20', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_2j330_a10t_lcw_jes_35smcINF_L1J100' : { 'HLTColl' : 'HLT_AntiKt10LCTopoTrimmedPtFrac4SmallR20Jets_jes', 'RefChain' : 'HLT_j85_L1J20', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_j45_pf_ftf_L1J15' : { 'HLTColl' : 'HLT_AntiKt4EMPFlowJets_subjesIS_ftf', 'RefChain' : 'NONE', 'OfflineColl' : 'NONE' },
'HLT_j45_pf_subjesgscIS_ftf_L1J15' : { 'HLTColl' : 'HLT_AntiKt4EMPFlowJets_subjesgscIS_ftf', 'RefChain' : 'NONE', 'OfflineColl' : 'NONE' },
'HLT_j85_pf_ftf_L1J20' : { 'HLTColl' : 'HLT_AntiKt4EMPFlowJets_subjesIS_ftf', 'RefChain' : 'HLT_j45_pf_ftf_L1J15', 'OfflineColl' : 'AntiKt4EMPFlowJets' },
'HLT_j45_pf_nojcalib_ftf_L1J15' : { 'HLTColl' : 'HLT_AntiKt4EMPFlowJets_nojcalib_ftf', 'RefChain' : 'NONE', 'OfflineColl' : 'NONE' },
'HLT_j45_csskpf_nojcalib_ftf_L1J15' : { 'HLTColl' : 'HLT_AntiKt4EMPFlowCSSKJets_nojcalib_ftf', 'RefChain' : 'NONE', 'OfflineColl' : 'NONE' },
'HLT_10j40_pf_subresjesgscIS_ftf_L14J15' : { 'HLTColl' : 'HLT_AntiKt4EMPFlowJets_subresjesgscIS_ftf', 'RefChain' : 'NONE', 'OfflineColl' : 'NONE' },
}
# Legacy
JetCollections['Legacy'] = [
'HLT_xAOD__JetContainer_a4tcemsubjesISFS', # default small-R
'HLT_xAOD__JetContainer_a10r_tcemsubjesISFS', # a10r
'HLT_xAOD__JetContainer_a10tclcwsubjesFS', # a10
'HLT_xAOD__JetContainer_a10ttclcwjesFS', # a10t
]
Chain2JetCollDict['Legacy'] = {
'HLT_j420' : 'HLT_xAOD__JetContainer_a4tcemsubjesISFS',
'HLT_j260_320eta490' : 'HLT_xAOD__JetContainer_a4tcemsubjesISFS',
'HLT_j460_a10r_L1J100' : 'HLT_xAOD__JetContainer_a10r_tcemsubjesISFS',
'HLT_j460_a10_lcw_subjes_L1J100' : 'HLT_xAOD__JetContainer_a10tclcwsubjesFS',
'HLT_j460_a10t_lcw_jes_L1J100' : 'HLT_xAOD__JetContainer_a10ttclcwjesFS',
'HLT_2j330_a10t_lcw_jes_35smcINF_L1J100' : 'HLT_xAOD__JetContainer_a10ttclcwjesFS',
'HLT_5j70_0eta240' : 'HLT_xAOD__JetContainer_a4tcemsubjesISFS',
'HLT_3j200' : 'HLT_xAOD__JetContainer_a4tcemsubjesISFS',
}
TurnOnCurves['Legacy'] = { # ref chain, offline jet coll
'HLT_j420' : ['HLT_j175','AntiKt4EMTopoJets'],
'HLT_j260_320eta490' : ['HLT_j45_320eta490','AntiKt4EMTopoJets'],
'HLT_j460_a10r_L1J100' : ['HLT_j175','AntiKt4EMTopoJets'],
'HLT_j460_a10_lcw_subjes_L1J100' : ['HLT_j175','AntiKt4EMTopoJets'],
'HLT_j460_a10t_lcw_jes_L1J100' : ['HLT_j175','AntiKt4EMTopoJets'],
'HLT_2j330_a10t_lcw_jes_35smcINF_L1J100' : ['HLT_j175','AntiKt4EMTopoJets'],
'HLT_3j200' : ['HLT_j175','AntiKt4EMTopoJets'],
JetCollections['Legacy'] = {
'HLT_xAOD__JetContainer_a4tcemsubjesISFS' : { 'MatchTo' : 'NONE' }, # default small-R
'HLT_xAOD__JetContainer_a10r_tcemsubjesISFS' : { 'MatchTo' : 'NONE' }, # a10r
'HLT_xAOD__JetContainer_a10tclcwsubjesFS' : { 'MatchTo' : 'NONE' }, # a10
'HLT_xAOD__JetContainer_a10ttclcwjesFS' : { 'MatchTo' : 'NONE' }, # a10t
}
JetColls2Match['Legacy'] = {
Chains2Monitor['Legacy'] = {
'HLT_j420' : { 'HLTColl' : 'HLT_xAOD__JetContainer_a4tcemsubjesISFS', 'RefChain' : 'HLT_j175', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_j260_320eta490' : { 'HLTColl' : 'HLT_xAOD__JetContainer_a4tcemsubjesISFS', 'RefChain' : 'HLT_j45_320eta490', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_j460_a10r_L1J100' : { 'HLTColl' : 'HLT_xAOD__JetContainer_a10r_tcemsubjesISFS', 'RefChain' : 'HLT_j175', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_j460_a10_lcw_subjes_L1J100' : { 'HLTColl' : 'HLT_xAOD__JetContainer_a10tclcwsubjesFS', 'RefChain' : 'HLT_j175', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_j460_a10t_lcw_jes_L1J100' : { 'HLTColl' : 'HLT_xAOD__JetContainer_a10ttclcwjesFS', 'RefChain' : 'HLT_j175', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_2j330_a10t_lcw_jes_35smcINF_L1J100' : { 'HLTColl' : 'HLT_xAOD__JetContainer_a10ttclcwjesFS', 'RefChain' : 'HLT_j175', 'OfflineColl' : 'AntiKt4EMTopoJets' },
'HLT_5j70_0eta240' : { 'HLTColl' : 'HLT_xAOD__JetContainer_a4tcemsubjesISFS', 'RefChain' : 'NONE', 'OfflineColl' : 'NONE' },
'HLT_3j200' : { 'HLTColl' : 'HLT_xAOD__JetContainer_a4tcemsubjesISFS', 'RefChain' : 'HLT_j175', 'OfflineColl' : 'AntiKt4EMTopoJets' },
}
#########################################################
......@@ -187,11 +159,12 @@ def TrigJetMonConfig(inputFlags):
from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
cfg = ComponentAccumulator()
# Match HLT to offline jets
for j1,j2 in JetColls2Match[InputType].items():
name = 'Matching_{}_{}'.format(j1,j2)
alg = CompFactory.JetMatcherAlg(name, JetContainerName1=j1,JetContainerName2=j2)
cfg.addEventAlgo(alg)
# Match HLT jets to offline jets
for hltColl,collDict in JetCollections[InputType].items():
if collDict['MatchTo'] != 'NONE':
name = 'Matching_{}_{}'.format(hltColl,collDict['MatchTo'])
alg = CompFactory.JetMatcherAlg(name, JetContainerName1=hltColl,JetContainerName2=collDict['MatchTo'])
cfg.addEventAlgo(alg)
# The following class will make a sequence, configure algorithms, and link
# them to GenericMonitoringTools
......@@ -221,7 +194,8 @@ def TrigJetMonConfig(inputFlags):
monitorConf.toAlg(helper)
# Loop over HLT jet chains
for chain,jetcoll in Chain2JetCollDict[InputType].items():
for chain,chainDict in Chains2Monitor[InputType].items():
jetcoll = chainDict['HLTColl']
# kinematic plots
if AthenaMT:
chainMonitorConfT = jetChainMonitoringConfig(inputFlags,jetcoll,chain,AthenaMT,True)
......@@ -229,13 +203,8 @@ def TrigJetMonConfig(inputFlags):
chainMonitorConfF = jetChainMonitoringConfig(inputFlags,jetcoll,chain,AthenaMT,False)
chainMonitorConfF.toAlg(helper)
# efficiency plots
refChain = 'NONE'
offlineJetColl = 'NONE'
if chain in TurnOnCurves[InputType]:
refChain = TurnOnCurves[InputType][chain][0]
offlineJetColl = TurnOnCurves[InputType][chain][1]
if offlineJetColl != 'NONE' and refChain != 'NONE':
effMonitorConf = jetEfficiencyMonitoringConfig(inputFlags,jetcoll,offlineJetColl,chain,refChain,AthenaMT)
if chainDict['RefChain'] != 'NONE' and chainDict['OfflineColl'] != 'NONE':
effMonitorConf = jetEfficiencyMonitoringConfig(inputFlags,jetcoll,chainDict['OfflineColl'],chain,chainDict['RefChain'],AthenaMT)
effMonitorConf.toAlg(helper)
cfg.merge(helper.result())
......@@ -334,7 +303,7 @@ def jetMonitoringConfig(inputFlags,jetcoll,athenaMT):
print('bottomlevel = '+str(conf.bottomLevelDir))
group = monhelper.addGroup(parentAlg, conf.Group, conf.topLevelDir+'/'+conf.bottomLevelDir+'/NoTriggerSelection/')
# define the histogram
group.defineHistogram('ptdiff',title='title', type="TH1F", path='MatchedJets_{}'.format(JetColls2Match[InputType][jetcoll]), xbins=100 , xmin=-100000, xmax=100000. ,)
group.defineHistogram('ptdiff',title='title', type="TH1F", path='MatchedJets_{}'.format(JetCollections[InputType][jetcoll]['MatchTo']), xbins=100 , xmin=-100000, xmax=100000. ,)
# Declare a configuration dictionnary for a JetContainer
if isOnline:
......@@ -352,8 +321,8 @@ def jetMonitoringConfig(inputFlags,jetcoll,athenaMT):
else:
for hist in ExtraLargeROnlineHists: conf.appendHistos(hist)
# Add matched jets plots
if jetcoll in JetColls2Match[InputType]:
matchedJetColl = JetColls2Match[InputType][jetcoll]
if JetCollections[InputType][jetcoll]['MatchTo'] != 'NONE':
matchedJetColl = JetCollections[InputType][jetcoll]['MatchTo']
name = 'jetMatched_{}_{}'.format(jetcoll,matchedJetColl)
jetmatchKey = '{}.matched_{}'.format(jetcoll,matchedJetColl)
jetptdiffKey = '{}.ptdiff_{}'.format(jetcoll,matchedJetColl)
......@@ -598,11 +567,12 @@ if __name__=='__main__':
helper = AthMonitorCfgHelper(ConfigFlags,'TrigJetMonitorAlgorithm')
cfg.merge(helper.result()) # merge it to add the sequence needed to add matchers
# Match HLT to offline jets
for j1,j2 in JetColls2Match[InputType].items():
name = 'Matching_{}_{}'.format(j1,j2)
alg = CompFactory.JetMatcherAlg(name, JetContainerName1=j1,JetContainerName2=j2)
cfg.addEventAlgo(alg,sequenceName='AthMonSeq_TrigJetMonitorAlgorithm') # Add matchers to monitoring alg sequence
# Match HLT jets to offline jets
for hltColl,collDict in JetCollections[InputType].items():
if collDict['MatchTo'] != 'NONE':
name = 'Matching_{}_{}'.format(hltColl,collDict['MatchTo'])
alg = CompFactory.JetMatcherAlg(name, JetContainerName1=hltColl,JetContainerName2=collDict['MatchTo'])
cfg.addEventAlgo(alg,sequenceName='AthMonSeq_TrigJetMonitorAlgorithm') # Add matchers to monitoring alg sequence
# Loop over L1 jet collectoins
for jetcoll in L1JetCollections:
......@@ -627,7 +597,8 @@ if __name__=='__main__':
monitorConf.toAlg(helper)
# Loop over HLT jet chains
for chain,jetcoll in Chain2JetCollDict[InputType].items():
for chain,chainDict in Chains2Monitor[InputType].items():
jetcoll = chainDict['HLTColl']
# kinematic plots
if AthenaMT:
chainMonitorConfT = jetChainMonitoringConfig(ConfigFlags,jetcoll,chain,AthenaMT,True)
......@@ -635,13 +606,8 @@ if __name__=='__main__':
chainMonitorConfF = jetChainMonitoringConfig(ConfigFlags,jetcoll,chain,AthenaMT,False)
chainMonitorConfF.toAlg(helper)
# efficiency plots
refChain = 'NONE'
offlineJetColl = 'NONE'
if chain in TurnOnCurves[InputType]:
refChain = TurnOnCurves[InputType][chain][0]
offlineJetColl = TurnOnCurves[InputType][chain][1]
if offlineJetColl != 'NONE' and refChain != 'NONE':
effMonitorConf = jetEfficiencyMonitoringConfig(ConfigFlags,jetcoll,offlineJetColl,chain,refChain,AthenaMT)
if chainDict['RefChain'] != 'NONE' and chainDict['OfflineColl'] != 'NONE':
effMonitorConf = jetEfficiencyMonitoringConfig(ConfigFlags,jetcoll,chainDict['OfflineColl'],chain,chainDict['RefChain'],AthenaMT)
effMonitorConf.toAlg(helper)
cfg.merge(helper.result())
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment