JetMonitoringStandard Namespace Reference

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

Variables
list	commonHistoSpecs
list	jvfHistosSpec
list	topoHistosSpec
list	pflowHistosSpec

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()

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

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

Definition at line 175 of file JetMonitoringStandard.py.

def jetMonAlgConfig(  jetName, inputFlags, truthJetName='', trigger=''):
    """returns a specification of a JetMonitoringAlg (in the form of a JetMonAlgSpec dictionnary).
    """
    
    # we use a specialized dictionnary (JetMonAlgSpec) which will be translated into the final C++ tool
    jetAlgConfig = JetMonAlgSpec(
        jetName+"MonAlg",
        JetContainerName = jetName,
        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()

JetMonitoringStandard.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 217 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

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.