Functions
def	getChains2Monitor (inputFlags, monMode)

def	getEtaRange (chain)
	Helpful functions. More...

def	getBinningFromThreshold (chain, varname)

def	getHTBinning (chain, binwidth)

def	addFlavourTagVariables (conf, network_prefix, flavs="cub")

def	getJetCopyAlg (injets, outjets)

def	getL1JetCopyAlg (injets, outjets)

def	TrigJetMonConfig (inputFlags)

def	basicJetMonAlgSpec (jetcoll, isOnline)

def	basicHIJetMonAlgSpec (jetcoll, isOnline)

def	jetMonitoringConfig (inputFlags, jetcoll, jetCollDict, monMode)

def	l1JetMonitoringConfig (inputFlags, jetColl, jetDict, monMode, chain='', matched=False)

def	jetChainMonitoringConfig (inputFlags, jetcoll, chain, onlyUsePassingJets=True)

def	jetEfficiencyMonitoringConfig (inputFlags, onlinejetcoll, offlinejetcoll, chain, refChain)

Variables
	logger

	copySuffix
	constants More...

	OfflineJetCollections
	Offline jet collections to monitor. More...

	L1JetCollections
	L1 jet collections and chains to monitor. More...

	match_smallRL1_OfflineJets_List

	match_largeRL1_OfflineJets_List

	match_HIL1_OfflineJets_List

	items

	d_items

	errmsg

	l1Coll2MatcherKey

	Chain2L1JetCollDict

	JetCollections
	HLT jet collections and chains to monitor. More...

	ExtraOfflineHists

	ExtraSmallROnlineHists

	title

	ExtraLargeROnlineHists

	ExtraOnlineNJetHists

	OfflineScaleMomenta

	OnlineScaleMomenta

	parser

	action

	dest

	default

	args

	RunTruth

	GenOfflineR10PF

	PrintDetailedConfig

	inputFile

	flags

	Files

	isMC

	HISTFileName

	monMode

	Chains2Monitor

	cfg

	outputlist

	comp

	key

	EMPFlowCSSK

	label

	AntiKt10EMPFlowCSSK

	ptmin

	modifiers

	AntiKt10EMPFlowCSSKSoftDrop

	substrmods

	ZCut

	Beta

	ItemList

	helper

	scalestring

	name

	alg

	JetContainerName1

	hltColl

	JetContainerName2

	JetCalibScale

	sequenceName

	offjetColl

	L1JetContainerName1

	l1jetColl

	matchjetcoll

	MatchL1

	l1jetconf

	l1chainconf

	offlineMonitorConf

	monitorConf

	jetcoll

	chainMonitorConfT

	chainMonitorConfF

	effMonitorConf

	withDetails

Function Documentation

◆ addFlavourTagVariables()

def TrigJetMonitorAlgorithm.addFlavourTagVariables	(	conf,
		network_prefix,
		flavs = `"cub"`
	)

Definition at line 308 of file TrigJetMonitorAlgorithm.py.

 def addFlavourTagVariables(conf, network_prefix, flavs="cub"):
     cutname='pt20'
     fillerTools = []
     for f in flavs:
       xvar = f"{network_prefix}_p{f}"
       varname = f"ftag_p{f}"
       fillerTools += [HistoSpec(varname, xvar=xvar, bins=(70, -0.2, 1.2), title=f"{varname};{varname};;Entries")]
     fastDipsSelectSpec = SelectSpec(f"{network_prefix}_{cutname}", '20<pt:GeV&|eta|<3.2', path='NoTriggerSelection/'+cutname, FillerTools=fillerTools)
     conf.appendHistos(fastDipsSelectSpec)
  

◆ basicHIJetMonAlgSpec()

def TrigJetMonitorAlgorithm.basicHIJetMonAlgSpec	(	jetcoll,
		isOnline
	)

Definition at line 644 of file TrigJetMonitorAlgorithm.py.

 def basicHIJetMonAlgSpec(jetcoll,isOnline):
   # we use a specialized dictionnary (JetMonAlgSpec) which will be translated into the final C++ tool
   path = 'NoTriggerSelection' if isOnline else 'standardHistos/'
  
   TopLevelDir  = 'HLT/JetMon/'
   TopLevelDir += 'Online/' if isOnline else 'Offline/'
  
   jetcollFolder = jetcoll
   jetcollFolder=jetcoll.replace(f"_{copySuffix}","")
   Conf = JetMonAlgSpec(jetcoll+"Mon",JetContainerName = jetcoll, defaultPath = path, topLevelDir=TopLevelDir, bottomLevelDir=jetcollFolder, failureOnMissingContainer=False)
  
   # Now start filling the histo spec list
   Conf.appendHistos(
     
     #See knownHistos in JetStandardHistoSpecs.py for the list of standard specification.
     #or we can directly add our custom histo specification in the form of a HistoSpec:
     #the basic call is : HistoSpec( variable, histobins, title='histotile;xtitle,ytitle')
     
     #Say we want a 2nd 'pt' plot but with a different binning than in the standard spec.
     #WARNING : we can not re-use the same spec name in a given JetMonitoringAlg !!!
     #so we give a new name AND we specify the actual variable with the argument 'xvar'
     #(the ':GeV' means the variable is to be set at GeV scale)
     #HistoSpec( 'lowpt',  (100,0,150) , title='p_{T};p_{T} [GeV];', xvar='pt:GeV'),            
     #An equivalent solution would have been to clone the existing spec like in :
     #knownHistos.pt.clone('lowpt',bins= (100,0,200) ),
  
     #2D histos are usually refered to by concatenating vars with a ';' as in 'varx;vary' 
     #if the 'vax;vary' alias doesn't exist in knownHistos but 'varx' and 'vary'
     #do exist, then a spec fot 'vax;vary' will be automatically generated.
  
     #Jet multiplicity histograms can be added by using an EventHistoSpec
     #Their specifications (pT cut, ET cut, eta cuts) must be defined in the knownEventVar dictionary within JetStandardHistoSpecs.py
     #The following line is an example for a jet multiplicity histogram with ET>40 GeV, 1.0<|eta|<2.0, and binning of (10,0,10):
     #EventHistoSpec('njetsEt40Eta1_2', (10,0,10), title='NJetsEt40Eta1_2;NJetsEt40Eta1_2;Entries' ),
  
     #To select on multiple variables simultaneously, simply combine the selection strings via &
     #Example below to select on ET > 100 GeV and |eta| > 3.2:
     #SelectSpec( 'ETeta', '100<et:GeV&|eta|<3.2', path, FillerTools = ["pt","et","m","eta"] )
  
     #TProfile2D : just use 3 variables. For now the sytem will automatically
     #interpret it as a TProfile2D (the 3rd variable being profiled)
     #"phi;eta;e",  --> Average Energy vs pt and eta
      
     #another possible selections : only sub-leading jets and highJVF
     #SelectSpec( 'subleading',
     #           '', # no selection on variables
     #           SelectedIndex=1, # force 2nd (sub-leading) jet (we would set 0 for leading jets)
     #           path='standardHistos', # force the path where the histos are saved in the final ROOT file
     #           FillerTools = [
     #               "pt",
     #               "m",
     #           ] ),
     #SelectSpec( 'highJVF',
     #           '0.3<JVF[0]', # JVF is a vector<float> for each jets. Here we cut on the 0th entry of this vector
     #           FillerTools = [
     #               "pt",
     #           ] ),
  
     #1D histos
     "pt",  
     "m",
     "eta",
     "phi",
     "et",
     
     #2D histos
     "eta;phi", # phi vs eta
     HistoSpec( 'pt:GeV;m:GeV',  (100,0,400, 100,0,400) , title='p_{T} vs mass;p_{T} [GeV];m [GeV];Entries'),
     HistoSpec( 'eta;et:GeV',  (100,-5,5, 100,0,400) , title='#eta vs e_{T};#eta;E_{T} [GeV];Entries'),
     HistoSpec( 'phi;et:GeV',  (60,-math.pi,math.pi, 100,0,400) , title='#phi vs E_{T};#phi;E_{T} [GeV];Entries'),
  
     #Event selection
     SelectSpec( 'central', '|eta|<3.2', path, FillerTools = ["pt","et","m"] ),
     SelectSpec( 'forward', '3.2<|eta|', path, FillerTools = ["pt","et","m"] ),
     SelectSpec( 'pt60', '60<pt:GeV', path, FillerTools = ["m","phi","eta","eta;phi"]),
     
   )
  
   return Conf
  

◆ basicJetMonAlgSpec()

def TrigJetMonitorAlgorithm.basicJetMonAlgSpec	(	jetcoll,
		isOnline
	)

Definition at line 555 of file TrigJetMonitorAlgorithm.py.

 def basicJetMonAlgSpec(jetcoll,isOnline):
   # we use a specialized dictionnary (JetMonAlgSpec) which will be translated into the final C++ tool
   path = 'NoTriggerSelection' if isOnline else 'standardHistos/'
   minNjetBin = 1 if isOnline else 0
  
   TopLevelDir  = 'HLT/JetMon/'
   TopLevelDir += 'Online/' if isOnline else 'Offline/'
  
   jetcollFolder = jetcoll
   jetcollFolder=jetcoll.replace(f"_{copySuffix}","")
   Conf = JetMonAlgSpec(jetcoll+"Mon",JetContainerName = jetcoll, defaultPath = path, topLevelDir=TopLevelDir, bottomLevelDir=jetcollFolder, failureOnMissingContainer=False)
  
   # Now start filling the histo spec list
   Conf.appendHistos(
     
     #See knownHistos in JetStandardHistoSpecs.py for the list of standard specification.
     #or we can directly add our custom histo specification in the form of a HistoSpec:
     #the basic call is : HistoSpec( variable, histobins, title='histotile;xtitle,ytitle')
     
     #Say we want a 2nd 'pt' plot but with a different binning than in the standard spec.
     #WARNING : we can not re-use the same spec name in a given JetMonitoringAlg !!!
     #so we give a new name AND we specify the actual variable with the argument 'xvar'
     #(the ':GeV' means the variable is to be set at GeV scale)
     #HistoSpec( 'lowpt',  (100,0,150) , title='p_{T};p_{T} [GeV];', xvar='pt:GeV'),            
     #An equivalent solution would have been to clone the existing spec like in :
     #knownHistos.pt.clone('lowpt',bins= (100,0,200) ),
  
     #2D histos are usually refered to by concatenating vars with a ';' as in 'varx;vary' 
     #if the 'vax;vary' alias doesn't exist in knownHistos but 'varx' and 'vary'
     #do exist, then a spec fot 'vax;vary' will be automatically generated.
     
  
     #Jet multiplicity histograms can be added by using an EventHistoSpec
     #Their specifications (pT cut, ET cut, eta cuts) must be defined in the knownEventVar dictionary within JetStandardHistoSpecs.py
     #The following line is an example for a jet multiplicity histogram with ET>40 GeV, 1.0<|eta|<2.0, and binning of (10,0,10):
     #EventHistoSpec('njetsEt40Eta1_2', (10,0,10), title='NJetsEt40Eta1_2;NJetsEt40Eta1_2;Entries' ),
  
     #To select on multiple variables simultaneously, simply combine the selection strings via &
     #Example below to select on ET > 100 GeV and |eta| > 3.2:
     #SelectSpec( 'ETeta', '100<et:GeV&|eta|<3.2', path, FillerTools = ["pt","et","m","eta"] )
  
     #TProfile2D : just use 3 variables. For now the sytem will automatically
     #interpret it as a TProfile2D (the 3rd variable being profiled)
     #"phi;eta;e",  --> Average Energy vs pt and eta
      
     #another possible selections : only sub-leading jets and highJVF
     #SelectSpec( 'subleading',
     #           '', # no selection on variables
     #           SelectedIndex=1, # force 2nd (sub-leading) jet (we would set 0 for leading jets)
     #           path='standardHistos', # force the path where the histos are saved in the final ROOT file
     #           FillerTools = [
     #               "pt",
     #               "m",
     #           ] ),
     #SelectSpec( 'highJVF',
     #           '0.3<JVF[0]', # JVF is a vector<float> for each jets. Here we cut on the 0th entry of this vector
     #           FillerTools = [
     #               "pt",
     #           ] ),
  
     #1D histos
     "pt",  
     "m",
     "eta",
     "phi",
     "e",
     "et",
     
     #2D histos
     "pt;m",    # mass vs pt
     "eta;phi", # phi vs eta
     "eta;e",   # energy vs eta
     "phi;e",   # energy vs phi
  
     #Event selection
     SelectSpec( 'central', '|eta|<3.2', path, FillerTools = ["pt","et","m"] ),
     SelectSpec( 'forward', '3.2<|eta|', path, FillerTools = ["pt","et","m"] ),
     SelectSpec( 'lowmu', 'avgMu<30', path, isEventVariable=True, FillerTools = ["pt","et","m","phi","eta"]),
     SelectSpec( 'highmu', '30<avgMu', path, isEventVariable=True, FillerTools = ["pt","et","m","phi","eta"]),
     
     EventHistoSpec('njets', (25,minNjetBin,25), title='NJets;NJets;Entries' ),
     EventHistoSpec('njetsPt20', (25,minNjetBin,25), title='NJetsPt20;NJetsPt20;Entries' ),
     EventHistoSpec('njetsPt50', (25,minNjetBin,25), title='NJetsPt50;NJetsPt50;Entries' ),
     
   )
  
   return Conf
  
 # Basic selection of histograms common for heavy ion online and offline jets

◆ getBinningFromThreshold()

def TrigJetMonitorAlgorithm.getBinningFromThreshold	(	chain,
		varname
	)

Definition at line 263 of file TrigJetMonitorAlgorithm.py.

 def getBinningFromThreshold(chain,varname):
   #default binning if nothing below applies
   xbins, xmin, xmax = 160,0.,800000.
   #pt and et binning based on threshold
   if varname == "pt" or varname == "et":
     if 'noalg' in chain:
         if 'jJ500' in chain or 'J400' in chain: return 160,xmin,800000 
         else: return 100,xmin,500000 # good enough for L1 jJ40 & jJ100
     else:
         #threshold = int(chain.split("_")[1].split('j')[1])
         threshold = int(re.search(r'\d+',chain.split("_")[1].split('j')[1]).group())
     if threshold < 50:
       return 40, 0., 100000.
     if threshold < 120:
       return 36, 20000., 200000.
  
     xbins = 40
     xmin = 50000.+100000.*(int(threshold/100)-1) #example: threshold = 330 -> 250 to 450; threshold = 420 -> 350 to 550
     if threshold % 100 == 0: #gives enough low bins if threshold is an exact divider of 100 GeV such as 3j200
       xmin = 1000.*(threshold - 100.)
     xmax = xmin + 200000.
     if "a10" in chain: # efficiency curve broader for large-R jets
       xmin = xmin - 50000. 
       xmax = xmax + 50000. 
       if "pf" in chain:
         xmax = xmax + 50000. # needed to include efficiency plateau for large-R PFlow chains
         if "smc" in chain:
           xmax = xmax + 50000. # efficiency plateau even higher for a10 pdf smc chains due to imperfect calibration
   #mass binning for large-R smc chains
   elif varname == "m":
     xbins, xmin, xmax = 35, 30000., 100000.
   return xbins, xmin, xmax
  

◆ getChains2Monitor()

def TrigJetMonitorAlgorithm.getChains2Monitor	(	inputFlags,
		monMode
	)

Definition at line 175 of file TrigJetMonitorAlgorithm.py.

 def getChains2Monitor(inputFlags, monMode):
  
   Chains2Monitor = dict()
   from TrigConfigSvc.TriggerConfigAccess import getHLTMonitoringAccess
   monAccess = getHLTMonitoringAccess(inputFlags)
   
   # set HLT jet collection, reference chain and offline jet collection
   # for turn-on curves 
   ListOfMonChains = monAccess.monitoredChains(signatures="jetMon", monLevels = ["shifter","t0"])
  
   default_dict = {"HLTColl": "NONE", "RefChain": "NONE" ,"OfflineColl": "NONE"}
   Chains2Monitor[monMode] = dict((chain, default_dict.copy()) for chain in ListOfMonChains)  
  
   if monMode == 'HI':
     for chainName in Chains2Monitor['HI']:
       if '_ion_' in chainName:
           Chains2Monitor['HI'][chainName]["HLTColl"] = "HLT_AntiKt4HIJets"
           Chains2Monitor['HI'][chainName]["OfflineColl"] = "AntiKt4HIJets"
       else: 
           Chains2Monitor['HI'][chainName]["HLTColl"] = "HLT_AntiKt4EMPFlowJets_jes_ftf"
           if not inputFlags.Common.doExpressProcessing:
             Chains2Monitor['HI'][chainName]["OfflineColl"] = "AntiKt4EMPFlowJets" #we do not have EMPFlowJets in HI express stream
           else:
             Chains2Monitor['HI'][chainName]["OfflineColl"] = "AntiKt4HIJets"
     # only HLT_noalg get efficiency curves by default, so...
     # these are additional hard-coded chains for efficiency monitoring
     if Chains2Monitor['HI'].get('HLT_j60_ion_L1jJ40'): Chains2Monitor['HI']['HLT_j60_ion_L1jJ40'].update({"RefChain": "HLT_noalg_L1jTE50", "OfflineColl": "AntiKt4HIJets"})
     if Chains2Monitor['HI'].get('HLT_j75_ion_L1jJ60'): Chains2Monitor['HI']['HLT_j75_ion_L1jJ60'].update({"RefChain": "HLT_noalg_L1jTE50", "OfflineColl": "AntiKt4HIJets"})
     if Chains2Monitor['HI'].get('HLT_j85_ion_L1jJ60'): Chains2Monitor['HI']['HLT_j85_ion_L1jJ60'].update({"RefChain": "HLT_noalg_L1jTE50", "OfflineColl": "AntiKt4HIJets"})
     if Chains2Monitor['HI'].get('HLT_j150_ion_L1jJ90'): Chains2Monitor['HI']['HLT_j150_ion_L1jJ90'].update({"RefChain": "HLT_noalg_L1jTE50", "OfflineColl": "AntiKt4HIJets"})
     if Chains2Monitor['HI'].get('HLT_j200_ion_L1jJ90'): Chains2Monitor['HI']['HLT_j200_ion_L1jJ90'].update({"RefChain": "HLT_noalg_L1jTE50", "OfflineColl": "AntiKt4HIJets"})
  
     if Chains2Monitor['HI'].get('HLT_j50f_ion_L1jJ40p30ETA49'): Chains2Monitor['HI']['HLT_j50f_ion_L1jJ40p30ETA49'].update({"RefChain": "HLT_noalg_L1jTE50", "OfflineColl": "AntiKt4HIJets"})
     if Chains2Monitor['HI'].get('HLT_j60f_ion_L1jJ40p30ETA49'): Chains2Monitor['HI']['HLT_j60f_ion_L1jJ40p30ETA49'].update({"RefChain": "HLT_noalg_L1jTE50", "OfflineColl": "AntiKt4HIJets"})
  
   elif monMode == "pp":  
     # logic to define HLTColl, RefChain, OfflineColl
     for chainName in Chains2Monitor['pp']:
       Chains2Monitor['pp'][chainName]["HLTColl"] = "HLT_AntiKt4EMTopoJets_subjesIS"
       if '_pf_' in chainName and 'a10' not in chainName: 
         Chains2Monitor['pp'][chainName]["HLTColl"] = "HLT_AntiKt4EMPFlowJets_subresjesgscIS_ftf"
         Chains2Monitor['pp'][chainName]["OfflineColl"] = "AntiKt4EMPFlowJets"
       elif 'a10' in chainName:
         if 'a10t' in chainName: Chains2Monitor['pp'][chainName]["HLTColl"] = "HLT_AntiKt10LCTopoTrimmedPtFrac4SmallR20Jets_jes"
         elif 'sd_cssk_pf' in chainName: Chains2Monitor['pp'][chainName]["HLTColl"] = "HLT_AntiKt10EMPFlowCSSKSoftDropBeta100Zcut10Jets_jes_ftf"
         elif 'a10r' in chainName: Chains2Monitor['pp'][chainName]["HLTColl"] = "HLT_AntiKt10EMTopoRCJets_subjesIS"
         else: Chains2Monitor['pp'][chainName]["HLTColl"] = "HLT_AntiKt10LCTopoJets_subjes"
       elif '_noalg_' in chainName: 
         Chains2Monitor['pp'][chainName]["RefChain"] = "HLT_j45_pf_ftf_preselj20_L1jJ40" # temporarily modify to using small-R jet in turn-on for both small and large-R jets to fix tier0 jet mon crash ATR-25800!!
         Chains2Monitor['pp'][chainName]["OfflineColl"] = "AntiKt4EMPFlowJets"
         if 'gLJ' in chainName: Chains2Monitor['pp'][chainName]["HLTColl"] = "HLT_AntiKt10EMPFlowCSSKSoftDropBeta100Zcut10Jets_jes_ftf"
       else: continue
  
  
     if Chains2Monitor['pp'].get('HLT_j420_L1jJ160'): Chains2Monitor['pp']['HLT_j420_L1jJ160'].update({"RefChain": "HLT_j85_L1jJ50", "OfflineColl": "AntiKt4EMPFlowJets"})
     if Chains2Monitor['pp'].get('HLT_3j200_L1jJ160'): Chains2Monitor['pp']['HLT_3j200_L1jJ160'].update({"RefChain": "HLT_j85_L1jJ50", "OfflineColl": "AntiKt4EMPFlowJets"})    
     if Chains2Monitor['pp'].get('HLT_4j120_L13jJ90'): Chains2Monitor['pp']['HLT_4j120_L13jJ90'].update({"RefChain": "HLT_j85_L1jJ50", "OfflineColl": "AntiKt4EMPFlowJets"}) 
     if Chains2Monitor['pp'].get('HLT_5j80_pf_ftf_presel5j50_L14jJ40'): Chains2Monitor['pp']['HLT_5j80_pf_ftf_presel5j50_L14jJ40'].update({"RefChain": "HLT_j45_pf_ftf_preselj20_L1jJ40", "OfflineColl": "AntiKt4EMPFlowJets"})    
     if Chains2Monitor['pp'].get('HLT_j400_pf_ftf_L1jJ160'): Chains2Monitor['pp']['HLT_j400_pf_ftf_L1jJ160'].update({"RefChain": "HLT_j85_pf_ftf_preselj50_L1jJ50", "OfflineColl": "AntiKt4EMPFlowJets"})    
     if Chains2Monitor['pp'].get('HLT_j400_pf_ftf_preselj225_L1jJ160'): Chains2Monitor['pp']['HLT_j400_pf_ftf_preselj225_L1jJ160'].update({"RefChain": "HLT_j85_pf_ftf_preselj50_L1jJ50", "OfflineColl": "AntiKt4EMPFlowJets"})
  
   else: 
     errmsg = 'Returned empty Chains2Monitor due to invalid monMode'
     raise RuntimeError(errmsg)
   return Chains2Monitor
  

◆ getEtaRange()

def TrigJetMonitorAlgorithm.getEtaRange ( chain )

Helpful functions.

Definition at line 245 of file TrigJetMonitorAlgorithm.py.

 def getEtaRange(chain):
   etaMin,etaMax = 0,2.5 # central jets by default
  
   if 'eta490_j' in chain: #workaround for the UPC (ultra-peripheral) trigger chains
     etaMin,etaMax = 0,4.9
   elif 'eta' in chain:
     etaParts    = chain.split('eta')
     etaMinTemp  = etaParts[0].split('_')
     etaMin      = etaMinTemp[len(etaMinTemp)-1]
     etaMin      = int(etaMin)/10
     etaMax      = etaParts[1].split('_')[0]
     etaMax      = int(etaMax)/10
  
   if 'f_ion' in chain: #workaround for the HLT forward triggers
     etaMin,etaMax = 3.2,4.9
   
   return etaMin,etaMax
  

◆ getHTBinning()

def TrigJetMonitorAlgorithm.getHTBinning	(	chain,
		binwidth
	)

Definition at line 296 of file TrigJetMonitorAlgorithm.py.

 def getHTBinning(chain,binwidth):
   parts = chain.split('HT')
   threshold = parts[1].split('_')[0]
   if 'XX' in threshold:
     threshold = threshold.split('XX')[0]
   xmin = int(0.9 * int(threshold))  # xmin to make the threshold visible
   xmax = xmin + 500
   xbins = int((xmax-xmin)/binwidth)-1
   return xbins, xmin, xmax
  
 # Add fast flavour-tag monitoring.
 # Adds a 20 GeV jet pT cut to avoid FPE WARNINGS from jets below min jet pT for RoI track association

◆ getJetCopyAlg()

def TrigJetMonitorAlgorithm.getJetCopyAlg	(	injets,
		outjets
	)

  Schedules JetCopier tool to make a shallow copy of
  the original offline/HLT jet container, for the JetMatcherAlg to decorate.
  This prevents our jet monitoring from decorating
  the original jet containers, which may end up being
  persistified in AOD/ESD (ATLASRECTS-7168,ATR-27980,ATR-26076)

Definition at line 369 of file TrigJetMonitorAlgorithm.py.

 def getJetCopyAlg(injets,outjets):
     '''
       Schedules JetCopier tool to make a shallow copy of
       the original offline/HLT jet container, for the JetMatcherAlg to decorate.
       This prevents our jet monitoring from decorating
       the original jet containers, which may end up being
       persistified in AOD/ESD (ATLASRECTS-7168,ATR-27980,ATR-26076)
     '''
     jcopy = CompFactory.JetCopier(
         "copier",
         InputJets = injets,
         DecorDeps=[],
         ShallowCopy=True,
         ShallowIO=True)
  
     jprovider = CompFactory.JetRecAlg(
         "jetalg_copy_"+outjets,
         Provider = jcopy,
         Modifiers = [],
         OutputContainer = outjets,
         MonTool = None)
  
     return jprovider
  

◆ getL1JetCopyAlg()

def TrigJetMonitorAlgorithm.getL1JetCopyAlg	(	injets,
		outjets
	)

  Schedules L1JetCopyAlgorithm to make a shallow copy of
  the original L1 jet container, for the JetMatcherAlg to decorate.
  This prevents our jet monitoring from decorating
  the original jet containers, which may end up being
  persistified in AOD/ESD (ATLASRECTS-7168,ATR-27980,ATR-26076).
  The L1JetCopyAlgorithm is a templated class (e.g. L1JetCopyAlgorithm<JTM_JetRoIContainer>).
  The python class name is what is generated by Athena during build time.
  The template types are defined in JTMContainers.h.

Definition at line 393 of file TrigJetMonitorAlgorithm.py.

 def getL1JetCopyAlg(injets,outjets):
     '''
       Schedules L1JetCopyAlgorithm to make a shallow copy of
       the original L1 jet container, for the JetMatcherAlg to decorate.
       This prevents our jet monitoring from decorating
       the original jet containers, which may end up being
       persistified in AOD/ESD (ATLASRECTS-7168,ATR-27980,ATR-26076).
       The L1JetCopyAlgorithm is a templated class (e.g. L1JetCopyAlgorithm<JTM_JetRoIContainer>).
       The python class name is what is generated by Athena during build time.
       The template types are defined in JTMContainers.h.
     '''
     jcopy_alg = None
     jcopy_alg_name = "l1jetcopy_alg_"+injets  
     if injets == "L1_jFexSRJetRoI":
         jcopy_alg = CompFactory.L1JetCopyAlgorithm_JTM_jFexSRJetRoIContainer_(jcopy_alg_name)    
     elif injets in ["L1_gFexSRJetRoI", "L1_gFexLRJetRoI"]:
         jcopy_alg = CompFactory.L1JetCopyAlgorithm_JTM_gFexJetRoIContainer_(jcopy_alg_name)    
     else:
         raise ValueError(f"L1 jet container {injets} not recognised")
     jcopy_alg.JetInContainerName = injets
     jcopy_alg.JetOutContainerName = outjets
  
     return jcopy_alg
  

◆ jetChainMonitoringConfig()

def TrigJetMonitorAlgorithm.jetChainMonitoringConfig	(	inputFlags,
		jetcoll,
		chain,
		onlyUsePassingJets = `True`
	)

Function to configures some algorithms in the monitoring system.

Definition at line 904 of file TrigJetMonitorAlgorithm.py.

 def jetChainMonitoringConfig(inputFlags,jetcoll,chain,onlyUsePassingJets=True):
    '''Function to configures some algorithms in the monitoring system.'''
  
    jetcollFolder = jetcoll
    chainFolder = chain
  
    #if not athenaMT:
    #  onlyUsePassingJets = False #does not work for legacy samples yet
    jetMonAlgSpecName = chain+"TrigMon"
    if not onlyUsePassingJets:
      chainFolder = chainFolder + "/ExpertHistos"
      jetMonAlgSpecName = jetMonAlgSpecName + "_ExpertHistos"
  
    # Define helper functions to automatize ET & eta selection strings for NJet histograms of chains
    def getThreshold(parts):
      return parts[1].split('_')[0]
  
    def getEtaRangeString(chain):
      etaMin, etaMax = 0, 32
      if 'eta490_j' in chain: #workaround for the upc chains
        etaMin, etaMax = 0, 49
      elif 'eta' in chain:
        etaParts    = chain.split('eta')
        etaMinTemp  = etaParts[0].split('_')
        etaMin      = etaMinTemp[len(etaMinTemp)-1]
        etaMax      = etaParts[1].split('_')[0]
        if int(etaMin) > 0 : etaMin = str(int(int(etaMin)/10))
        if int(etaMax) > 0 : etaMax = str(int(int(etaMax)/10))
      if 'f_ion' in chain:
        etaMin, etaMax = 32, 49
      return 'Eta{}_{}'.format(etaMin,etaMax)
  
    def getNjetHistName(chain):
      NjetHistName = 'NONE'
      parts         = chain.split('j')
      # check if it is a multi-threshold multijet chain or a single-threshold multijet chain
      multiplicity  = parts[0].split('_')[1] # for single-threshold multijet chains
      if (chain.count('_j')-chain.count('_jes')) > 1  or multiplicity != '':
        NjetHistName = 'njetsEt{}{}'.format(getThreshold(parts),getEtaRangeString(chain))
      return NjetHistName
  
  
    trigConf = JetMonAlgSpec( # the usual JetMonAlgSpec 
      jetMonAlgSpecName,
      JetContainerName = jetcoll,
      TriggerChain = chain,
      defaultPath = chainFolder,
      topLevelDir="HLT/JetMon/Online/",
      bottomLevelDir=jetcollFolder,
      failureOnMissingContainer=True,
      onlyPassingJets=onlyUsePassingJets,
      isExpressStreamJob=inputFlags.Common.doExpressProcessing,
    )
    
    trigConf.appendHistos(
            "pt",
            "m",
            "eta",
            "et",
            "phi",
    )
    for hist in ExtraOnlineNJetHists: trigConf.appendHistos(EventHistoSpec(hist, (20,0,25), title=hist+';'+hist+';Entries'))
    # Add NjetEt and NjetPt histograms for simple scenarios
    if 'ht' not in chain and 'HT' not in chain and 'dijet' not in chain and 'DIJET' not in chain and 'fbdj' not in chain and 'noalg' not in chain:
      NjetHistName = getNjetHistName(chain)
      from JetMonitoring.JetStandardHistoSpecs import knownEventVar
      if knownEventVar.get(NjetHistName,None) is not None and NjetHistName not in ExtraOnlineNJetHists: #avoids duplication warnings for some chains
        trigConf.appendHistos(
          EventHistoSpec(NjetHistName, (25,0,25), title=NjetHistName+';'+NjetHistName+';Entries' ),
        )
      NjetHistName = NjetHistName.replace('Et','Pt')
      if knownEventVar.get(NjetHistName,None) is not None and NjetHistName not in ExtraOnlineNJetHists:
        trigConf.appendHistos(
          EventHistoSpec(NjetHistName, (25,0,25), title=NjetHistName+';'+NjetHistName+';Entries' ),
        )
    if 'ftf' in chain and 'a10' not in chain: # track-based JVT variables for FTF chains
      trigConf.appendHistos("Jvt")
      trigConf.appendHistos("JVFCorr")
      trigConf.appendHistos("JvtRpt")
  
    if 'ht' in chain or 'HT' in chain:
      def defineHistoForHTChain(conf, parentAlg, monhelper , path):
          # create a monitoring group with the histo path starting from the parentAlg
          group = monhelper.addGroup(parentAlg, conf.Group, conf.topLevelDir+jetcollFolder+'/')
          # define the histograms
          xbins, xmin, xmax = getHTBinning(chain,25) # bin width in GeV
          group.defineHistogram("jetHT;HT",title="Jet HT;H_{T} [GeV];Entries", type="TH1F", path=chainFolder, xbins=xbins , xmin=xmin, xmax=xmax ,)
      trigConf.appendHistos(ToolSpec('JetHistoHTFiller','JetHistoHTFiller_'+chain,MinPt=30.,MaxEta=3.2,FailureOnMissingContainer=False,
                                   defineHistoFunc=defineHistoForHTChain,Group='jetHT_'+jetcoll))
  
    return trigConf
  

◆ jetEfficiencyMonitoringConfig()

def TrigJetMonitorAlgorithm.jetEfficiencyMonitoringConfig	(	inputFlags,
		onlinejetcoll,
		offlinejetcoll,
		chain,
		refChain
	)

Function to configures some algorithms in the monitoring system.

Definition at line 996 of file TrigJetMonitorAlgorithm.py.

 def jetEfficiencyMonitoringConfig(inputFlags,onlinejetcoll,offlinejetcoll,chain,refChain):
    '''Function to configures some algorithms in the monitoring system.'''
  
    jetcollFolder = onlinejetcoll
    chainFolder = chain
  
    # We schedule a new JetAlg which will be acting only when a TriggerChain fired (using the TriggerChain from the base classes).
    # We'll plot 1 histo build by a dedicated JetHistoTriggEfficiency tool.
    # So we'll have to explicitely give a specification via the generic dicionnary 'ToolSpec'
    # This implies defining a little function which declares to the monitoring framework which variables to histogram and how.
    #  this is done here.
    def defineHistoForJetTrigg(conf, parentAlg, monhelper , path):
        # create a monitoring group with the histo path starting from the parentAlg
        group = monhelper.addGroup(parentAlg, conf.Group, conf.topLevelDir+jetcollFolder+'/')
        # define the histogram, give them individual names so they don't overwrite each other
        append = "offlineCut_"+conf.name.split("_")[-1] if "offlineCut" in conf.name else "noOfflineCut"
        histname = "trigEff_vs_"+conf.Var.Name+"_"+append
        xbins, xmin, xmax = getBinningFromThreshold(chain,conf.Var.Name)
        group.defineHistogram('trigPassed,jetVar;'+histname, title=histname, type="TEfficiency",
                              path=chainFolder,
                              xbins=xbins , xmin=xmin, xmax=xmax,)
        
    # Get jet index and eta selection for offline jets
    validchain = chain.replace('noalg','j0') 
    parts        = validchain.split('j')
    multiplicity = parts[0].split('_')[1]
    if multiplicity != '': index = int(multiplicity) - 1 # single-threhold multijet chains
    else: index = 0 # single-jet chain
    etaMin,etaMax = getEtaRange(chain)
  
    from JetMonitoring.JetMonitoringConfig import retrieveVarToolConf
    trigConf = JetMonAlgSpec( # the usual JetMonAlgSpec 
        chain+"TrigEffMon",
        JetContainerName          = offlinejetcoll,
        TriggerChain              = refChain, # reference chain
        defaultPath               = chainFolder,
        topLevelDir               = "HLT/JetMon/Online/",
        bottomLevelDir            = jetcollFolder,
        failureOnMissingContainer = True,
        onlyPassingJets           = False,
        )
    trigConf.appendHistos(
        SelectSpec( 'eff', '{}<|eta|<{}'.format(etaMin,etaMax), chainFolder, SelectedIndex=index, FillerTools = [
            # we pass directly the ToolSpec
            ToolSpec('JetHistoTriggEfficiency', chain,
                     # below we pass the Properties of this JetHistoTriggEfficiency tool :
                     Group='jetTrigGroup_'+chain,
                     Var=retrieveVarToolConf("pt"), # In this context we can not just pass a str alias to describe a histo variable
                                                    # so we use retrieveVarToolConf("pt") which returns a full specification for the "pt" histo variable.
                     ProbeTrigChain=chain,defineHistoFunc=defineHistoForJetTrigg),
        ] ),
    )
  
    if 'smc' in chain:
      trigConf.appendHistos(
              SelectSpec( 'm50', '50<m:GeV&{}<|eta|<{}'.format(etaMin,etaMax), chainFolder, SelectedIndex=index, FillerTools = [
                ToolSpec('JetHistoTriggEfficiency', chain+'_offlineCut_m50',
                  Group='jetTrigGroup_'+chain+'_m50',
                  Var=retrieveVarToolConf("pt"), # In this context we can not just pass a str alias to describe a histo variable
                  ProbeTrigChain=chain,defineHistoFunc=defineHistoForJetTrigg
                ),
              ] ),
              SelectSpec( 'et500', '500<et:GeV&{}<|eta|<{}'.format(etaMin,etaMax), chainFolder, SelectedIndex=index, FillerTools = [
                ToolSpec('JetHistoTriggEfficiency', chain+'_offlineCut_et500',
                  Group='jetTrigGroup_'+chain+'_et500',
                  Var=retrieveVarToolConf("m"), # In this context we can not just pass a str alias to describe a histo variable
                  SortJets=True,
                  ProbeTrigChain=chain,defineHistoFunc=defineHistoForJetTrigg
                ),
              ] ),
      )
  
    return trigConf
  

◆ jetMonitoringConfig()

def TrigJetMonitorAlgorithm.jetMonitoringConfig	(	inputFlags,
		jetcoll,
		jetCollDict,
		monMode
	)

Function to configures some algorithms in the monitoring system.

Definition at line 724 of file TrigJetMonitorAlgorithm.py.

 def jetMonitoringConfig(inputFlags,jetcoll,jetCollDict,monMode):
    '''Function to configures some algorithms in the monitoring system.'''
  
    isOnline  = True if 'HLT' in jetcoll else False
  
    if monMode == 'HI': # Heavy ion jet monitoring histos
       conf      = basicHIJetMonAlgSpec(jetcoll,isOnline)
    else:
       conf      = basicJetMonAlgSpec(jetcoll,isOnline)
  
    jetCollMonDetails = jetCollDict[monMode][jetcoll]
  
    # Declare a configuration dictionnary for a JetContainer
    if isOnline:
      if 'AntiKt4' in jetcoll or 'a4tcem' in jetcoll:
        if monMode == 'pp': #Use extra histos for pp only
         for hist in ExtraSmallROnlineHists: conf.appendHistos(hist)
        if 'ftf' in jetcoll: # dedicated histograms for FTF chains
          conf.appendHistos("Jvt")
          conf.appendHistos("JVFCorr")
          conf.appendHistos("JvtRpt")
          conf.appendHistos("SumPtTrkPt500[0]")
          conf.appendHistos("NumTrkPt1000[0]")
          conf.appendHistos("TrackWidthPt1000[0]")
          if 'PF' in jetcoll: # dedicated histograms for online PFlow jets
            conf.appendHistos("SumPtChargedPFOPt500[0]")
            conf.appendHistos("fCharged")
            if "subresjesgscIS" in jetcoll:
                addFlavourTagVariables(conf,"fastDIPS20211215")
                addFlavourTagVariables(conf,"GN120230331")
                addFlavourTagVariables(conf,"tlaGN220240122")
        if 'fastftag' in jetcoll:
            addFlavourTagVariables(conf,"fastDips")
            addFlavourTagVariables(conf, "fastGN120230327")
            addFlavourTagVariables(conf,"fastGN220240122")
            addFlavourTagVariables(conf,"fastGNTau20240216", flavs=["tau", "u"])
        if 'EMTopo' in jetcoll: #dedicated histograms for online EMTopo jets
            conf.appendHistos("Timing")
      else:
        for hist in ExtraLargeROnlineHists: conf.appendHistos(hist)
      # Add matched jets plots
      if jetCollMonDetails['MatchTo'] != 'NONE':
        def defineHistoForHLTJetMatch(conf, parentAlg, monhelper , path):
            # create a monitoring group with the histo path starting from the parentAlg
            group = monhelper.addGroup(parentAlg, conf.Group, conf.topLevelDir+'/'+conf.bottomLevelDir+'/NoTriggerSelection/')
            # define the histograms
            for histname in [ 'ptdiff', 'energydiff', 'massdiff' ]: #defines which variable difference will be plotted
              group.defineHistogram(histname,title=histname, type="TH1F",
                                    path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                    xbins=100 , xmin=-100000., xmax=100000. ,)
              
            for histname in [ 'ptresp', 'energyresp', 'massresp' ]:
              group.defineHistogram(histname,title=histname, type="TH1F",
                                    path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                    xbins=100 , xmin=-2., xmax=2. ,)
              
            group.defineHistogram('ptresp,ptref;ptresp_vs_ptRef',title='ptresponse vs ptRef', type="TH2F",
                                  path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                  xbins=10 , xmin=-2., xmax=2., ybins=10, ymin=0., ymax=500000.,)
            
            group.defineHistogram('ptresp,etaref;ptresp_vs_etaRef',title='ptresponse vs etaRef', type="TH2F",
                                  path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                  xbins=10 , xmin=-2., xmax=2., ybins=10, ymin=-5., ymax=5.,)
  
            group.defineHistogram('ptref,ptresp;ptRef_vs_ptresp',title='ptRef vs ptresponse', type="TH2F",
                                  path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                  xbins=100 , xmin=0., xmax=400000., ybins=80, ymin=-2., ymax=2.,)
                                 
            group.defineHistogram('etaref,ptresp;etaRef_vs_ptresp',title='etaRef vs ptresponse', type="TH2F",
                                  path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                  xbins=100 , xmin=-5., xmax=5., ybins=80, ymin=-2., ymax=2.,)
            
        matchedJetColl   = jetCollMonDetails['MatchTo']
  
        # we can get specific calibration scales by adding e.g. '_EtaJESScale' to the strings
        jetmatchKey      = '{}.matched_{}'.format(jetcoll,matchedJetColl)
        jetptdiffKey     = '{}.ptdiff_{}'.format(jetcoll,matchedJetColl)
        jetenergydiffKey = '{}.energydiff_{}'.format(jetcoll,matchedJetColl)
        jetmassdiffKey   = '{}.massdiff_{}'.format(jetcoll,matchedJetColl)
        jetptrespKey     = '{}.ptresp_{}'.format(jetcoll,matchedJetColl)
        jetenergyrespKey = '{}.energyresp_{}'.format(jetcoll,matchedJetColl)
        jetmassrespKey   = '{}.massresp_{}'.format(jetcoll,matchedJetColl)
        jetptrefKey      = '{}.ptRef_{}'.format(jetcoll,matchedJetColl)
        jetetarefKey     = '{}.etaRef_{}'.format(jetcoll,matchedJetColl)
        name = 'jetMatched_{}_{}'.format(jetcoll,matchedJetColl)
        conf.appendHistos(ToolSpec('JetHistoMatchedFiller', name,
                                   JetMatchedKey=jetmatchKey, JetPtDiffKey=jetptdiffKey,
                                   JetEnergyDiffKey=jetenergydiffKey,
                                   JetMassDiffKey=jetmassdiffKey, JetPtRespKey=jetptrespKey,
                                   JetEnergyRespKey=jetenergyrespKey, JetMassRespKey=jetmassrespKey,
                                   JetPtRefKey=jetptrefKey,JetEtaRefKey=jetetarefKey,
                                   defineHistoFunc=defineHistoForHLTJetMatch,Group='matchedJets_'+jetcoll)
        )
    else: # offline
      for hist in ExtraOfflineHists: conf.appendHistos(hist)
      if 'AntiKt4' in jetcoll and monMode=="pp":
          conf.appendHistos(SelectSpec('LooseBadFailedJets', 'LooseBad',
                                       InverseJetSel=True,
                                       FillerTools = ["pt",
                                                      "phi",
                                                      "eta"])) #cleaning variables not applicable for large-R collections
      
          if 'PF' in jetcoll: # dedicated histograms for offline PFlow jets
              conf.appendHistos("SumPtChargedPFOPt500[0]")
              conf.appendHistos("fCharged")
          elif 'EMTopo' in jetcoll:
              conf.appendHistos("Timing")
      if jetCollMonDetails['MatchTo'] != 'NONE':
        def defineHistoForOfflineJetMatch(conf, parentAlg, monhelper , path):
          # create a monitoring group with the histo path starting from the parentAlg
          group = monhelper.addGroup(parentAlg, conf.Group, conf.topLevelDir+'/'+conf.bottomLevelDir+'/standardHistos/')
          # define the histograms
          for histname in [ 'ptdiff', 'energydiff', 'massdiff' ]: #defines which variable difference will be plotted
            group.defineHistogram(histname,title=histname, type="TH1F",
                                  path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                  xbins=100 , xmin=-100000., xmax=100000. ,)
            
          for histname in [ 'ptresp', 'energyresp', 'massresp' ]:
            group.defineHistogram(histname,title=histname, type="TH1F",
                                  path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                  xbins=100 , xmin=-2., xmax=2. ,)
            
          group.defineHistogram('ptresp,ptref;ptresp_vs_ptRef',title='ptresp vs ptRef', type="TH2F",
                                path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                xbins=10 , xmin=-2., xmax=2., ybins=10, ymin=0., ymax=500000.,)
          
          group.defineHistogram('ptresp,etaref;ptresp_vs_etaRef',title='ptresp vs etaRef', type="TH2F",
                                path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                xbins=10 , xmin=-2., xmax=2., ybins=10, ymin=-5., ymax=5.,)
  
          group.defineHistogram('ptref,ptresp;ptRef_vs_ptresp',title='ptRef vs ptresponse', type="TH2F",
                                  path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                  xbins=100 , xmin=0., xmax=400000., ybins=80, ymin=-2., ymax=2.,)
                                 
          group.defineHistogram('etaref,ptresp;etaRef_vs_ptresp',title='etaRef vs ptresponse', type="TH2F",
                                  path='MatchedJets_{}'.format(jetCollMonDetails['MatchTo']),
                                  xbins=100 , xmin=-5., xmax=5., ybins=80, ymin=-2., ymax=2.,)
          
        matchedJetColl   = jetCollMonDetails['MatchTo']
        jetmatchKey      = '{}.matched_{}'.format(jetcoll,matchedJetColl)
        jetptdiffKey     = '{}.ptdiff_{}'.format(jetcoll,matchedJetColl)
        jetenergydiffKey = '{}.energydiff_{}'.format(jetcoll,matchedJetColl)
        jetmassdiffKey   = '{}.massdiff_{}'.format(jetcoll,matchedJetColl)
        jetptrespKey     = '{}.ptresp_{}'.format(jetcoll,matchedJetColl)
        jetenergyrespKey = '{}.energyresp_{}'.format(jetcoll,matchedJetColl)
        jetmassrespKey   = '{}.massresp_{}'.format(jetcoll,matchedJetColl)
        jetptrefKey      = '{}.ptRef_{}'.format(jetcoll,matchedJetColl)
        jetetarefKey     = '{}.etaRef_{}'.format(jetcoll,matchedJetColl)
        name = 'jetMatched_{}_{}'.format(jetcoll,matchedJetColl)
        conf.appendHistos(ToolSpec('JetHistoMatchedFiller',name,
                                   JetMatchedKey=jetmatchKey, JetPtDiffKey=jetptdiffKey,
                                   JetEnergyDiffKey=jetenergydiffKey,
                                   JetMassDiffKey=jetmassdiffKey, JetPtRespKey=jetptrespKey,
                                   JetEnergyRespKey=jetenergyrespKey,
                                   JetMassRespKey=jetmassrespKey,
                                   JetPtRefKey=jetptrefKey, JetEtaRefKey=jetetarefKey,
                                   defineHistoFunc=defineHistoForOfflineJetMatch,Group='matchedJets_'+jetcoll)
        )
  
    return conf
  

◆ l1JetMonitoringConfig()

def TrigJetMonitorAlgorithm.l1JetMonitoringConfig	(	inputFlags,
		jetColl,
		jetDict,
		monMode,
		chain = `''`,
		matched = `False`
	)

Definition at line 885 of file TrigJetMonitorAlgorithm.py.

 def l1JetMonitoringConfig(inputFlags,jetColl, jetDict, monMode,chain='',matched=False):
  
   from TrigJetMonitoring.L1JetMonitoringConfig import L1JetMonAlg
   name = jetColl if chain=='' else jetColl+'_'+chain
  
   jetCollKey = jetColl
   jetColl = jetColl.replace(f"_{copySuffix}","")
  
   if not jetDict[monMode][jetCollKey]['MatchTo']:
     conf = L1JetMonAlg(name,jetColl,jetCollKey,chain)
   else:
     assert  len(jetDict[monMode][jetCollKey]['MatchTo']) == 2
     
     conf = L1JetMonAlg(name,jetColl,jetCollKey,chain,
                        matched,jetDict[monMode][jetCollKey]['MatchTo'][0],
                        jetDict[monMode][jetCollKey]['MatchTo'][1])
  
   return conf
  

◆ TrigJetMonConfig()

def TrigJetMonitorAlgorithm.TrigJetMonConfig ( inputFlags )

Definition at line 417 of file TrigJetMonitorAlgorithm.py.

 def TrigJetMonConfig(inputFlags):
  
   from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
   cfg = ComponentAccumulator()
  
   monMode = 'pp'
   if inputFlags.Reco.EnableHI: monMode = 'HI'
  
   Chains2Monitor = getChains2Monitor(inputFlags, monMode)
  
   # Protections
   # Add missing jet collections to JetCollections dict
   # (this can happen if a given chain uses a jet collection that is not listed in JetCollections) 
   # TODO: make more general
   for chain,chaindict in Chains2Monitor[monMode].items():
     if chaindict['HLTColl'] not in JetCollections[case]: # chain will not be monitored unless HLT collection is present in JetCollections
       JetCollections[case][chaindict['HLTColl']] = {'MatchTo': 'NONE'}
  
   # Match HLT jets to offline jets
   CopiedJetCollections = copy.deepcopy(JetCollections)
   for hltColl,collDict in JetCollections[monMode].items():
     if collDict['MatchTo'] != 'NONE':
       copiedhltColl = f'{hltColl}_{copySuffix}'
       CopiedJetCollections[monMode][copiedhltColl] = CopiedJetCollections[monMode].pop(hltColl)
       jetcopyalg = getJetCopyAlg(hltColl,copiedhltColl)
       jetcopyalg.ExtraInputs.add(('xAOD::TrigCompositeContainer',
                                   'StoreGateSvc+%s' % getRun3NavigationContainerFromInput(inputFlags)))
       cfg.addEventAlgo(jetcopyalg)
       for jetcalibscale in OnlineScaleMomenta:
         scalestring = "_"+jetcalibscale if jetcalibscale != "" else ""
         name = 'Matching_{}{}_{}'.format(hltColl,scalestring,collDict['MatchTo'])
         alg = CompFactory.JetMatcherAlg(name,
                                         JetContainerName1=copiedhltColl,
                                         JetContainerName2=collDict['MatchTo'],
                                         JetCalibScale=jetcalibscale)
         
         alg.ExtraInputs.add(('xAOD::TrigCompositeContainer',
                              'StoreGateSvc+%s' % getRun3NavigationContainerFromInput(inputFlags)))
         cfg.addEventAlgo(alg)
  
   # Match offline to offline jets
   CopiedOfflineJetCollections = copy.deepcopy(OfflineJetCollections)
   for offjetColl,collDict in OfflineJetCollections[monMode].items():
     if collDict['MatchTo'] != 'NONE':
       copiedjetcoll = f'{offjetColl}_{copySuffix}'
       CopiedOfflineJetCollections[monMode][copiedjetcoll] = CopiedOfflineJetCollections[monMode].pop(offjetColl)
       jetcopyalg = getJetCopyAlg(offjetColl,copiedjetcoll)
       cfg.addEventAlgo(jetcopyalg)
       for jetcalibscale in OfflineScaleMomenta:
         scalestring = "_"+jetcalibscale if jetcalibscale != "" else ""
         name = 'Matching_{}{}_{}'.format(offjetColl,scalestring,collDict['MatchTo'])
         alg = CompFactory.JetMatcherAlg(name,
                                         JetContainerName1=copiedjetcoll,
                                         JetContainerName2=collDict['MatchTo'],
                                         JetCalibScale=jetcalibscale)
         
         alg.ExtraInputs.add(('xAOD::TrigCompositeContainer',
                              'StoreGateSvc+%s' % getRun3NavigationContainerFromInput(inputFlags)))
         cfg.addEventAlgo(alg)
  
   # Make copy of every L1 jet collection
   # Then match L1 to offline as well as HLT jets
   CopiedL1JetCollections = copy.deepcopy(L1JetCollections)
   for l1jetColl,collDict in L1JetCollections[monMode].items():
     copiedl1jetColl = f'{l1jetColl}_{copySuffix}'
     CopiedL1JetCollections[monMode][copiedl1jetColl] = CopiedL1JetCollections[monMode].pop(l1jetColl)
     l1jetcopyalg = getL1JetCopyAlg(l1jetColl,copiedl1jetColl)
     l1jetcopyalg.ExtraInputs.add(('xAOD::TrigCompositeContainer',
                                   'StoreGateSvc+%s' % getRun3NavigationContainerFromInput(inputFlags)))
     cfg.addEventAlgo(l1jetcopyalg)
     for matchjetcoll in collDict['MatchTo']:
  
       kwds = {'name': 'Matching_{}_{}'.format(l1jetColl,matchjetcoll),
               l1Coll2MatcherKey[l1jetColl]: copiedl1jetColl,
               'JetContainerName2': matchjetcoll,
               'MatchL1': True
              }
               
       alg = CompFactory.JetMatcherAlg(**kwds)
       alg.ExtraInputs.add(('xAOD::TrigCompositeContainer',
                            'StoreGateSvc+%s' % getRun3NavigationContainerFromInput(inputFlags)))
       cfg.addEventAlgo(alg)
  
   # The following class will make a sequence, configure algorithms, and link
   # them to GenericMonitoringTools
   from AthenaMonitoring import AthMonitorCfgHelper
   helper = AthMonitorCfgHelper(inputFlags,'TrigJetMonitorAlgorithm')
    # Configure filter tools                                                                                                                                                                     
   from AthenaMonitoring.EventFlagFilterToolConfig import EventFlagFilterToolCfg
   from AthenaMonitoring.BadLBFilterToolConfig import LArBadLBFilterToolCfg
   # Loop over L1 jet collections
   for jetcoll in CopiedL1JetCollections[monMode]:
     l1jetconf = l1JetMonitoringConfig(inputFlags,jetcoll,CopiedL1JetCollections,monMode,'',True)
     alg=l1jetconf.toAlg(helper)
     alg.FilterTools = [ EventFlagFilterToolCfg(inputFlags),helper.resobj.popToolsAndMerge(LArBadLBFilterToolCfg(inputFlags))]
  
   # Loop over L1 jet chains
   for chain,jetcolls in Chain2L1JetCollDict[monMode].items():
     for jetcoll in jetcolls:
       l1chainconf = l1JetMonitoringConfig(inputFlags,jetcoll,L1JetCollections,monMode,chain)
       alg=l1chainconf.toAlg(helper)
       alg.FilterTools = [ EventFlagFilterToolCfg(inputFlags),helper.resobj.popToolsAndMerge(LArBadLBFilterToolCfg(inputFlags))]
  
   # Loop over offline jet collections
   for jetcoll in CopiedOfflineJetCollections[monMode]:
     offlineMonitorConf = jetMonitoringConfig(inputFlags,jetcoll,CopiedOfflineJetCollections,monMode)
     alg=offlineMonitorConf.toAlg(helper)
     alg.FilterTools = [ EventFlagFilterToolCfg(inputFlags),helper.resobj.popToolsAndMerge(LArBadLBFilterToolCfg(inputFlags))]
  
   # Loop over HLT jet collections
   for jetcoll in CopiedJetCollections[monMode]:
     monitorConf = jetMonitoringConfig(inputFlags,jetcoll,CopiedJetCollections,monMode)
     # then we turn the full specification into properly configured algorithm and tools.
     # we use the method 'toAlg()' defined for the specialized dictionnary 'JetMonAlgSpec'
     monitorConf.toAlg(helper)
  
   # Loop over HLT jet chains
   for chain,chainDict in Chains2Monitor[monMode].items():
     jetcoll = chainDict['HLTColl']
     # kinematic plots
     # only use passing jets
     chainMonitorConfT = jetChainMonitoringConfig(inputFlags,jetcoll,chain,True)
     alg=chainMonitorConfT.toAlg(helper)
     alg.FilterTools = [ EventFlagFilterToolCfg(inputFlags),helper.resobj.popToolsAndMerge(LArBadLBFilterToolCfg(inputFlags))]
     # all jets
     chainMonitorConfF = jetChainMonitoringConfig(inputFlags,jetcoll,chain,False)
     alg=chainMonitorConfF.toAlg(helper)
     alg.FilterTools = [ EventFlagFilterToolCfg(inputFlags),helper.resobj.popToolsAndMerge(LArBadLBFilterToolCfg(inputFlags))]
     # efficiency plots
     if chainDict['RefChain'] != 'NONE' and chainDict['OfflineColl'] != 'NONE':
       effMonitorConf = jetEfficiencyMonitoringConfig(inputFlags,jetcoll,chainDict['OfflineColl'],chain,chainDict['RefChain'])
       effMonitorConf.toAlg(helper)
  
   cfg.merge(helper.result())
   return cfg
  
  
 # Basic selection of histograms common for online and offline jets

Variable Documentation

◆ action

TrigJetMonitorAlgorithm.action

Definition at line 1076 of file TrigJetMonitorAlgorithm.py.

◆ alg

TrigJetMonitorAlgorithm.alg

Definition at line 1172 of file TrigJetMonitorAlgorithm.py.

◆ AntiKt10EMPFlowCSSK

TrigJetMonitorAlgorithm.AntiKt10EMPFlowCSSK

Definition at line 1139 of file TrigJetMonitorAlgorithm.py.

◆ AntiKt10EMPFlowCSSKSoftDrop

TrigJetMonitorAlgorithm.AntiKt10EMPFlowCSSKSoftDrop

Definition at line 1143 of file TrigJetMonitorAlgorithm.py.

◆ args

TrigJetMonitorAlgorithm.args

Definition at line 1080 of file TrigJetMonitorAlgorithm.py.

◆ Beta

TrigJetMonitorAlgorithm.Beta

Definition at line 1143 of file TrigJetMonitorAlgorithm.py.

◆ cfg

TrigJetMonitorAlgorithm.cfg

Definition at line 1120 of file TrigJetMonitorAlgorithm.py.

◆ Chain2L1JetCollDict

TrigJetMonitorAlgorithm.Chain2L1JetCollDict

Definition at line 111 of file TrigJetMonitorAlgorithm.py.

◆ chainMonitorConfF

TrigJetMonitorAlgorithm.chainMonitorConfF

Definition at line 1225 of file TrigJetMonitorAlgorithm.py.

◆ chainMonitorConfT

TrigJetMonitorAlgorithm.chainMonitorConfT

Definition at line 1222 of file TrigJetMonitorAlgorithm.py.

◆ Chains2Monitor

TrigJetMonitorAlgorithm.Chains2Monitor

Definition at line 1107 of file TrigJetMonitorAlgorithm.py.

◆ comp

TrigJetMonitorAlgorithm.comp

Definition at line 1129 of file TrigJetMonitorAlgorithm.py.

◆ copySuffix

TrigJetMonitorAlgorithm.copySuffix

constants

Definition at line 25 of file TrigJetMonitorAlgorithm.py.

◆ d_items

TrigJetMonitorAlgorithm.d_items

Definition at line 77 of file TrigJetMonitorAlgorithm.py.

◆ default

TrigJetMonitorAlgorithm.default

Definition at line 1076 of file TrigJetMonitorAlgorithm.py.

◆ dest

TrigJetMonitorAlgorithm.dest

Definition at line 1076 of file TrigJetMonitorAlgorithm.py.

◆ effMonitorConf

TrigJetMonitorAlgorithm.effMonitorConf

Definition at line 1229 of file TrigJetMonitorAlgorithm.py.

◆ EMPFlowCSSK

TrigJetMonitorAlgorithm.EMPFlowCSSK

Definition at line 1138 of file TrigJetMonitorAlgorithm.py.

◆ errmsg

TrigJetMonitorAlgorithm.errmsg

Definition at line 82 of file TrigJetMonitorAlgorithm.py.

◆ ExtraLargeROnlineHists

TrigJetMonitorAlgorithm.ExtraLargeROnlineHists

Definition at line 349 of file TrigJetMonitorAlgorithm.py.

◆ ExtraOfflineHists

TrigJetMonitorAlgorithm.ExtraOfflineHists

Definition at line 325 of file TrigJetMonitorAlgorithm.py.

◆ ExtraOnlineNJetHists

TrigJetMonitorAlgorithm.ExtraOnlineNJetHists

Definition at line 352 of file TrigJetMonitorAlgorithm.py.

◆ ExtraSmallROnlineHists

TrigJetMonitorAlgorithm.ExtraSmallROnlineHists

Definition at line 337 of file TrigJetMonitorAlgorithm.py.

◆ Files

TrigJetMonitorAlgorithm.Files

Definition at line 1099 of file TrigJetMonitorAlgorithm.py.

◆ flags

TrigJetMonitorAlgorithm.flags

Definition at line 1098 of file TrigJetMonitorAlgorithm.py.

◆ GenOfflineR10PF

TrigJetMonitorAlgorithm.GenOfflineR10PF

Definition at line 1082 of file TrigJetMonitorAlgorithm.py.

◆ helper

TrigJetMonitorAlgorithm.helper

Definition at line 1163 of file TrigJetMonitorAlgorithm.py.

◆ HISTFileName

TrigJetMonitorAlgorithm.HISTFileName

Definition at line 1101 of file TrigJetMonitorAlgorithm.py.

◆ hltColl

TrigJetMonitorAlgorithm.hltColl

Definition at line 1172 of file TrigJetMonitorAlgorithm.py.

◆ inputFile

TrigJetMonitorAlgorithm.inputFile

Definition at line 1086 of file TrigJetMonitorAlgorithm.py.

◆ isMC

TrigJetMonitorAlgorithm.isMC

Definition at line 1100 of file TrigJetMonitorAlgorithm.py.

◆ ItemList

TrigJetMonitorAlgorithm.ItemList

Definition at line 1156 of file TrigJetMonitorAlgorithm.py.

◆ items

TrigJetMonitorAlgorithm.items

Definition at line 71 of file TrigJetMonitorAlgorithm.py.

◆ JetCalibScale

TrigJetMonitorAlgorithm.JetCalibScale

Definition at line 1172 of file TrigJetMonitorAlgorithm.py.

◆ jetcoll

TrigJetMonitorAlgorithm.jetcoll

Definition at line 1219 of file TrigJetMonitorAlgorithm.py.

◆ JetCollections

TrigJetMonitorAlgorithm.JetCollections

HLT jet collections and chains to monitor.

Definition at line 152 of file TrigJetMonitorAlgorithm.py.

◆ JetContainerName1

TrigJetMonitorAlgorithm.JetContainerName1

Definition at line 1172 of file TrigJetMonitorAlgorithm.py.

◆ JetContainerName2

TrigJetMonitorAlgorithm.JetContainerName2

Definition at line 1172 of file TrigJetMonitorAlgorithm.py.

◆ key

TrigJetMonitorAlgorithm.key

Definition at line 1132 of file TrigJetMonitorAlgorithm.py.

◆ l1chainconf

TrigJetMonitorAlgorithm.l1chainconf

Definition at line 1202 of file TrigJetMonitorAlgorithm.py.

◆ l1Coll2MatcherKey

TrigJetMonitorAlgorithm.l1Coll2MatcherKey

Definition at line 93 of file TrigJetMonitorAlgorithm.py.

◆ l1jetColl

TrigJetMonitorAlgorithm.l1jetColl

Definition at line 1191 of file TrigJetMonitorAlgorithm.py.

◆ L1JetCollections

TrigJetMonitorAlgorithm.L1JetCollections

L1 jet collections and chains to monitor.

Definition at line 51 of file TrigJetMonitorAlgorithm.py.

◆ l1jetconf

TrigJetMonitorAlgorithm.l1jetconf

Definition at line 1197 of file TrigJetMonitorAlgorithm.py.

◆ L1JetContainerName1

TrigJetMonitorAlgorithm.L1JetContainerName1

Definition at line 1191 of file TrigJetMonitorAlgorithm.py.

◆ label

TrigJetMonitorAlgorithm.label

Definition at line 1138 of file TrigJetMonitorAlgorithm.py.

◆ logger

TrigJetMonitorAlgorithm.logger

Definition at line 12 of file TrigJetMonitorAlgorithm.py.

◆ match_HIL1_OfflineJets_List

TrigJetMonitorAlgorithm.match_HIL1_OfflineJets_List

Definition at line 56 of file TrigJetMonitorAlgorithm.py.

◆ match_largeRL1_OfflineJets_List

TrigJetMonitorAlgorithm.match_largeRL1_OfflineJets_List

Definition at line 55 of file TrigJetMonitorAlgorithm.py.

◆ match_smallRL1_OfflineJets_List

TrigJetMonitorAlgorithm.match_smallRL1_OfflineJets_List

Definition at line 53 of file TrigJetMonitorAlgorithm.py.

◆ matchjetcoll

TrigJetMonitorAlgorithm.matchjetcoll

Definition at line 1191 of file TrigJetMonitorAlgorithm.py.

◆ MatchL1

TrigJetMonitorAlgorithm.MatchL1

Definition at line 1191 of file TrigJetMonitorAlgorithm.py.

◆ modifiers

TrigJetMonitorAlgorithm.modifiers

Definition at line 1140 of file TrigJetMonitorAlgorithm.py.

◆ monitorConf

TrigJetMonitorAlgorithm.monitorConf

Definition at line 1212 of file TrigJetMonitorAlgorithm.py.

◆ monMode

TrigJetMonitorAlgorithm.monMode

Definition at line 1104 of file TrigJetMonitorAlgorithm.py.

◆ name

TrigJetMonitorAlgorithm.name

Definition at line 1171 of file TrigJetMonitorAlgorithm.py.

◆ offjetColl

TrigJetMonitorAlgorithm.offjetColl

Definition at line 1182 of file TrigJetMonitorAlgorithm.py.

◆ OfflineJetCollections

TrigJetMonitorAlgorithm.OfflineJetCollections

Offline jet collections to monitor.

Definition at line 31 of file TrigJetMonitorAlgorithm.py.

◆ offlineMonitorConf

TrigJetMonitorAlgorithm.offlineMonitorConf

Definition at line 1207 of file TrigJetMonitorAlgorithm.py.

◆ OfflineScaleMomenta

TrigJetMonitorAlgorithm.OfflineScaleMomenta

Definition at line 357 of file TrigJetMonitorAlgorithm.py.

◆ OnlineScaleMomenta

TrigJetMonitorAlgorithm.OnlineScaleMomenta

Definition at line 358 of file TrigJetMonitorAlgorithm.py.

◆ outputlist

TrigJetMonitorAlgorithm.outputlist

Definition at line 1123 of file TrigJetMonitorAlgorithm.py.

◆ parser

TrigJetMonitorAlgorithm.parser

Definition at line 1075 of file TrigJetMonitorAlgorithm.py.

◆ PrintDetailedConfig

TrigJetMonitorAlgorithm.PrintDetailedConfig

Definition at line 1083 of file TrigJetMonitorAlgorithm.py.

◆ ptmin

TrigJetMonitorAlgorithm.ptmin

Definition at line 1139 of file TrigJetMonitorAlgorithm.py.

◆ RunTruth

TrigJetMonitorAlgorithm.RunTruth

Definition at line 1081 of file TrigJetMonitorAlgorithm.py.

◆ scalestring

TrigJetMonitorAlgorithm.scalestring

Definition at line 1170 of file TrigJetMonitorAlgorithm.py.

◆ sequenceName

TrigJetMonitorAlgorithm.sequenceName

Definition at line 1174 of file TrigJetMonitorAlgorithm.py.

◆ substrmods

TrigJetMonitorAlgorithm.substrmods

Definition at line 1143 of file TrigJetMonitorAlgorithm.py.

◆ title

TrigJetMonitorAlgorithm.title

Definition at line 338 of file TrigJetMonitorAlgorithm.py.

◆ withDetails

TrigJetMonitorAlgorithm.withDetails

Definition at line 1237 of file TrigJetMonitorAlgorithm.py.

◆ ZCut

TrigJetMonitorAlgorithm.ZCut

Definition at line 1143 of file TrigJetMonitorAlgorithm.py.

Functions

Variables

Function Documentation

◆ addFlavourTagVariables()

◆ basicHIJetMonAlgSpec()

◆ basicJetMonAlgSpec()

◆ getBinningFromThreshold()

◆ getChains2Monitor()

◆ getEtaRange()

◆ getHTBinning()

◆ getJetCopyAlg()

◆ getL1JetCopyAlg()

◆ jetChainMonitoringConfig()

◆ jetEfficiencyMonitoringConfig()

◆ jetMonitoringConfig()

◆ l1JetMonitoringConfig()

◆ TrigJetMonConfig()

Variable Documentation

◆ action

◆ alg

◆ AntiKt10EMPFlowCSSK

◆ AntiKt10EMPFlowCSSKSoftDrop

◆ args

◆ Beta

◆ cfg

◆ Chain2L1JetCollDict

◆ chainMonitorConfF

◆ chainMonitorConfT

◆ Chains2Monitor

◆ comp

◆ copySuffix

◆ d_items

◆ default

◆ dest

◆ effMonitorConf

◆ EMPFlowCSSK

◆ errmsg

◆ ExtraLargeROnlineHists

◆ ExtraOfflineHists

◆ ExtraOnlineNJetHists

◆ ExtraSmallROnlineHists

◆ Files

◆ flags

◆ GenOfflineR10PF

◆ helper

◆ HISTFileName

◆ hltColl

◆ inputFile

◆ isMC

◆ ItemList

◆ items

◆ JetCalibScale

◆ jetcoll

◆ JetCollections

◆ JetContainerName1

◆ JetContainerName2

◆ key

◆ l1chainconf

◆ l1Coll2MatcherKey

◆ l1jetColl

◆ L1JetCollections

◆ l1jetconf

◆ L1JetContainerName1

◆ label

◆ logger

◆ match_HIL1_OfflineJets_List

◆ match_largeRL1_OfflineJets_List

◆ match_smallRL1_OfflineJets_List

◆ matchjetcoll

◆ MatchL1

◆ modifiers

◆ monitorConf

◆ monMode

◆ name

◆ offjetColl

◆ OfflineJetCollections

◆ offlineMonitorConf

◆ OfflineScaleMomenta

◆ OnlineScaleMomenta

◆ outputlist