globalSelection = ""
##
# necessary updates for validation directories
###
athenaMonTool = JetMonitoringTool(
HistoTools=[
commonPhysValTool("AntiKt4LCTopoJets",
akt4refContainer,
globalSelection=globalSelection),
commonPhysValTool("AntiKt4EMTopoJets",
akt4refContainer,
globalSelection=globalSelection),
# the containers that are commented out are kept so as to make it possible to swicth them on in the future if needed
#commonPhysValTool( "AntiKt10LCTopoJets" ),
commonPhysValTool("AntiKt10LCTopoTrimmedPtFrac5SmallR20Jets"),
#commonPhysValTool( "AntiKt10PV0TracksJets" ),
#commonPhysValTool( "AntiKt10TruthJets" ),
commonPhysValTool("AntiKt10TruthTrimmedPtFrac5SmallR20Jets"),
commonPhysValTool("AntiKt4EMPFlowJets"),
commonPhysValTool("AntiKt4TruthJets"),
#commonPhysValTool( "AntiKt10TruthWZJets" ),
],
IntervalType=8) # 8 == HistoGroupBase::all
if (isMC == False):
athenaMonTool = JetMonitoringTool(
HistoTools=[
commonPhysValTool("AntiKt4LCTopoJets",
akt4refContainer,
]),
selectionAndHistos(jhm.alljets, [jhm.Timing]),
# Other types of selection with simple expression.
# Only jets satisfying the expression will fill the histos
# expresion must be in the form : '1.23<var<3.45'
# 'var' must be a jet attribute (float, int or vector<float>)
selectionAndHistos("-1<eta<1", [
"basickinematics",
]),
selectionAndHistos("eta<1", [
"basickinematics",
]),
selectionAndHistos("JVF[0]<0.8", ["basickinematics", "JVF[0]"]),
]
jetMonTool = JetMonitoringTool(HistoTools=[akt4Filler])
ToolSvc += jetMonTool
monMan.AthenaMonTools += [jetMonTool]
from AthenaCommon.AppMgr import ServiceMgr
from GaudiSvc.GaudiSvcConf import THistSvc
ServiceMgr += THistSvc()
svcMgr.THistSvc.Output += ["PhysVal DATAFILE='PhysVal.root' OPT='RECREATE'"]
from AthenaCommon.AppMgr import theApp
ServiceMgr.MessageSvc.OutputLevel = INFO
ServiceMgr.MessageSvc.defaultLimit = 1000000
PlotAverageE=True,
PlotNJet=True),
jhm.ptN,
jhm.Timing,
jhm.EMFrac,
jhm.HECFrac,
jhm.LArQuality,
]
return filler
athenaMonTool = JetMonitoringTool(
HistoTools=[
commonMonitoringTool(
"AntiKt4LCTopoJets"
), # if truth is present, we could add : , "AntiKt4TruthJets" ,
commonMonitoringTool("AntiKt4EMTopoJets"),
commonMonitoringTool("AntiKt10LCTopoJets"),
],
IntervalType=6,
) # 6 is 'Interval_t::run' interval
if monitorTracks:
athenaMonTool.HistoTools += [commonMonitoringTool("AntiKt3PV0TrackJets")]
#if jetFlags.useTracks:
# athenaMonTool.HistoTools += [ commonMonitoringTool( "AntiKt3PV0TrackJets" ) ]
#cbg
athenaMonTool.FilterTools += [monbadlb]
ToolSvc += athenaMonTool
athenaMonTool_LB = JetMonitoringTool(
if refcontainer:
# efficiency
filler.HistoTools += [
jhm.effresponse,
]
filler.HistoTools['effresponse'].RefContainer = refcontainer
#filler.OutputLevel =2
return filler
athenaMonTool = JetMonitoringTool(
HistoTools=[
commonMonitoringTool(
"AntiKt4HIJets"
), # if truth is present, we could add : , "AntiKt4TruthJets" ,
],
IntervalType=6,
) # 6 is 'Interval_t::run' interval
if monitorTracks:
athenaMonTool.HistoTools += [commonMonitoringTool("AntiKt4HITrackJets")]
#cbg
athenaMonTool.FilterTools += [monbadlb]
#ToolSvc += athenaMonTool
athenaMonTool_LB = JetMonitoringTool("JetMonitoring_LB",
HistoTools=[
"ptN", "Timing", "EMFrac", "HECFrac",
"LArQuality",