Functions
	jetMonAlgConfig (jetName, inputFlags, truthJetName='', trigger='')
	standardJetMonitoring (inputFlags)

Variables
list	commonHistoSpecs
list	jvfHistosSpec
list	topoHistosSpec
list	pflowHistosSpec
dict	calibToolConfigurations

Detailed Description

@file JetMonitoringExample.py
@author P-A. Delsart
@date 2019-03-12
@brief Main python configuration for the Run III Jet Monitoring

Function Documentation

◆ jetMonAlgConfig()

jetMonAlgConfig	(	jetName,
		inputFlags,
		truthJetName = '',
		trigger = '' )

returns a specification of a JetMonitoringAlg (in the form of a JetMonAlgSpec dictionnary).

Definition at line 179 of file JetMonitoringStandard.py.

def jetMonAlgConfig(  jetName, inputFlags, truthJetName='', trigger=''):
    """returns a specification of a JetMonitoringAlg (in the form of a JetMonAlgSpec dictionnary).
    """
    
    # Just a protection in case PhysVal is using this as well
    if not inputFlags.Input.isMC and jetName != "AntiKt4HIJets":
        from AthenaConfiguration.ComponentFactory import CompFactory
        configfile, calibArea, calibseq_def = calibToolConfigurations[jetName]
        toolname = "jetcalib_DQ_{0}".format(jetName)
        jct = CompFactory.JetCalibrationTool(toolname,
                                             JetCollection = jetName.replace('Jets',''),
                                             ConfigFile = configfile,
                                             CalibArea = calibArea,
                                             CalibSequence = calibseq_def,
                                             IsData = True)
 
        # we use a specialized dictionnary (JetMonAlgSpec) which will be translated into the final C++ tool
        jetAlgConfig = JetMonAlgSpec(
            jetName+"MonAlg",
            JetContainerName = jetName,
            applyLatestCalibration = True,
            JetCalibTool = jct,
            TriggerChain = trigger,
        )
    else:
        jetAlgConfig = JetMonAlgSpec(
            jetName+"MonAlg",
            JetContainerName = jetName,
            applyLatestCalibration = False,
            TriggerChain = trigger,
        )
 
 
    # the list of histos specifications
    histoSpecs = []
 
    # then add pre-defined lists as defined above :
    histoSpecs += commonHistoSpecs 
 
    if inputFlags.DQ.DataType is not DQDataType.Cosmics:
        histoSpecs += jvfHistosSpec
 
    if 'Topo' in jetName:
        histoSpecs += topoHistosSpec
    if 'PFlow' in jetName:
        histoSpecs += pflowHistosSpec
 
 
    if truthJetName != "" :
        # then add histos showing efficiency and pT responses vs True
        from JetMonitoring.JetStandardHistoSpecs import responseAndEffSpecMap
        if truthJetName not in responseAndEffSpecMap:
            print( "ERROR !! can't schedule a JetHistoResponseAndEff for truth container : ",truthJetName, ". No specification available" )
            return None
 
        histoSpecs +=[ responseAndEffSpecMap[truthJetName]  ]
 
    # finally all all histos specs to our JetMonitoringAlg specification :
    jetAlgConfig.appendHistos( * histoSpecs)
 
    return jetAlgConfig
    
 
 

◆ standardJetMonitoring()

standardJetMonitoring ( inputFlags )

Standard jet monitoring function to be inserted from top-level algs. 
returns an a component accumulator as given by AthMonitorCfgHelper.result()
Details of what goes into jet monitoring is implemented by dedicated functions such as jetMonAlgConfig().

Definition at line 243 of file JetMonitoringStandard.py.

def standardJetMonitoring(inputFlags):
    """Standard jet monitoring function to be inserted from top-level algs. 
    returns an a component accumulator as given by AthMonitorCfgHelper.result()
    Details of what goes into jet monitoring is implemented by dedicated functions such as jetMonAlgConfig().
    """
 
    from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
    rv = ComponentAccumulator()
 
    # do not run monitoring in RAWtoESD
    if inputFlags.DQ.Environment == 'tier0Raw':
        return rv
 
    from AthenaMonitoring import AthMonitorCfgHelper
    helper = AthMonitorCfgHelper(inputFlags,'JetMonitoring')
 
    # create a list of JetMonitoringAlg specifications
    jetAlgConfs = []
 
    if inputFlags.Reco.EnableHI:
        if inputFlags.Tracking.doUPC:
            jetAlgConfs.append(jetMonAlgConfig( "AntiKt4EMPFlowJets", inputFlags))
        else:
            jetAlgConfs.append(jetMonAlgConfig( "AntiKt4HIJets", inputFlags))
    else:
        # use the helper function defined above :
        #jetMonAlgConfig( "AntiKt4LCTopoJets", truthJetName="AntiKt4TruthJets"),     #How can we make sure truth jets are available ??
        jetAlgConfs.append(jetMonAlgConfig( "AntiKt4LCTopoJets", inputFlags))
        jetAlgConfs.append(jetMonAlgConfig( "AntiKt4EMTopoJets", inputFlags))
        jetAlgConfs.append(jetMonAlgConfig( "AntiKt4EMPFlowJets", inputFlags))
    
    # Configure filter tools 
    from AthenaMonitoring.EventFlagFilterToolConfig import EventFlagFilterToolCfg
#    from AthenaMonitoring.AtlasReadyFilterConfig import AtlasReadyFilterCfg
    from AthenaMonitoring.BadLBFilterToolConfig import LArBadLBFilterToolCfg
#    from AthenaMonitoring.FilledBunchFilterToolConfig import FilledBunchFilterToolCfg 
 
    # schedule each JetMonitoringAlg by invoking the toAlg() methods of the config specification
    for conf in jetAlgConfs:        
       alg = conf.toAlg(helper)
       alg.FilterTools = [ EventFlagFilterToolCfg(inputFlags),helper.resobj.popToolsAndMerge(LArBadLBFilterToolCfg(inputFlags))]    
 
    rv.merge(helper.result())  # the AthMonitorCfgHelper returns an accumulator to be used by the general configuration system.
    return rv

Variable Documentation

◆ calibToolConfigurations

dict JetMonitoringStandard.calibToolConfigurations

Initial value:

=  {
    "AntiKt4EMPFlowJets":("AntiKt4EMPFlow_MC23a_PreRecR22_Phase2_CalibConfig_ResPU_EtaJES_GSC_241208_InSitu.config","00-04-83","JetArea_Residual_EtaJES"),
    "AntiKt4EMTopoJets":("PreRec_R22_EMTopo_ResPU_EtaJES_October23_231024.config","00-04-82","JetArea_Residual_EtaJES"),
    "AntiKt4LCTopoJets":("JES_MC15cRecommendation_May2016_rel21.config","00-04-77","JetArea_Residual_EtaJES")
}

Definition at line 173 of file JetMonitoringStandard.py.

◆ commonHistoSpecs

list JetMonitoringStandard.commonHistoSpecs

Definition at line 16 of file JetMonitoringStandard.py.

◆ jvfHistosSpec

list JetMonitoringStandard.jvfHistosSpec

Initial value:

=  [
    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","m","eta","phi",
                ] ),
]

Definition at line 129 of file JetMonitoringStandard.py.

◆ pflowHistosSpec

list JetMonitoringStandard.pflowHistosSpec

Initial value:

=  [
    # histos common to all PFlow jets. These are defined in JetStandardHistoSpecs.py
    "SumPtChargedPFOPt500[0]",
    "SumPtTrkPt500[0]",
    "NumTrkPt500[0]",
    "NumTrkPt1000[0]",
    "fCharged",
]

Definition at line 164 of file JetMonitoringStandard.py.

◆ topoHistosSpec

list JetMonitoringStandard.topoHistosSpec

Initial value:

=  [
    # histos common to all topo jets. Energy By Sampling layer plots.
    HistoSpec("PreSamplerB",(100,0.,200.),title="Energy in PreSamplerB;E;Entries",xvar=VarSpec('EnergyPerSampling[0]')),
    HistoSpec("EMB1",(100,0.,200.),title="Energy in EMB1;E;Entries",xvar=VarSpec('EnergyPerSampling[1]')),
    HistoSpec("EMB2",(100,0.,200.),title="Energy in EMB2;E;Entries",xvar=VarSpec('EnergyPerSampling[2]')),
    HistoSpec("EMB3",(100,0.,200.),title="Energy in EMB3;E;Entries",xvar=VarSpec('EnergyPerSampling[3]')),
    HistoSpec("PreSamplerE",(100,0.,200.),title="Energy in PreSamplerE;E;Entries",xvar=VarSpec('EnergyPerSampling[4]')),
    HistoSpec("EME1",(100,0.,200.),title="Energy in EME2;E;Entries",xvar=VarSpec('EnergyPerSampling[5]')),
    HistoSpec("EME2",(100,0.,200.),title="Energy in EME2;E;Entries",xvar=VarSpec('EnergyPerSampling[6]')),
    HistoSpec("EME3",(100,0.,200.),title="Energy in EME3;E;Entries",xvar=VarSpec('EnergyPerSampling[7]')),
    HistoSpec("HEC0",(100,0.,200.),title="Energy in HEC0;E;Entries",xvar=VarSpec('EnergyPerSampling[8]')),
    HistoSpec("HEC1",(100,0.,200.),title="Energy in HEC1;E;Entries",xvar=VarSpec('EnergyPerSampling[9]')),
    HistoSpec("HEC2",(100,0.,200.),title="Energy in HEC2;E;Entries",xvar=VarSpec('EnergyPerSampling[10]')),
    HistoSpec("HEC3",(100,0.,200.),title="Energy in HEC3;E;Entries",xvar=VarSpec('EnergyPerSampling[11]')),
    HistoSpec("TileBar0",(100,0.,200.),title="Energy in TileBar0;E;Entries",xvar=VarSpec('EnergyPerSampling[12]')),
    HistoSpec("TileBar1",(100,0.,200.),title="Energy in TileBar1;E;Entries",xvar=VarSpec('EnergyPerSampling[13]')),
    HistoSpec("TileBar2",(100,0.,200.),title="Energy in TileBar1;E;Entries",xvar=VarSpec('EnergyPerSampling[14]')),
    HistoSpec("TileGap1",(100,0.,200.),title="Energy in TileGap1;E;Entries",xvar=VarSpec('EnergyPerSampling[15]')),
    HistoSpec("TileGap2",(100,0.,200.),title="Energy in TileGap2;E;Entries",xvar=VarSpec('EnergyPerSampling[16]')),
    HistoSpec("TileGap3",(100,0.,200.),title="Energy in TileGap3;E;Entries",xvar=VarSpec('EnergyPerSampling[17]')),
    HistoSpec("TileExt0",(100,0.,200.),title="Energy in TileExt0;E;Entries",xvar=VarSpec('EnergyPerSampling[18]')),
    HistoSpec("TileExt1",(100,0.,200.),title="Energy in TileExt1;E;Entries",xvar=VarSpec('EnergyPerSampling[19]')),
    HistoSpec("TileExt2",(100,0.,200.),title="Energy in TileExt2;E;Entries",xvar=VarSpec('EnergyPerSampling[20]')),
    HistoSpec("FCAL0",(100,0.,200.),title="Energy in FCAL0;E;Entries",xvar=VarSpec('EnergyPerSampling[21]')),
    HistoSpec("FCAL1",(100,0.,200.),title="Energy in FCAL1;E;Entries",xvar=VarSpec('EnergyPerSampling[22]')),
    HistoSpec("FCAL2",(100,0.,200.),title="Energy in FCAL2;E;Entries",xvar=VarSpec('EnergyPerSampling[23]')),
]

Definition at line 136 of file JetMonitoringStandard.py.

Functions

Variables

Detailed Description

Function Documentation

◆ jetMonAlgConfig()

◆ standardJetMonitoring()

Variable Documentation

◆ calibToolConfigurations

◆ commonHistoSpecs

◆ jvfHistosSpec

◆ pflowHistosSpec

◆ topoHistosSpec