python.HLT.Jet.JetMenuSequencesConfig Namespace Reference

Classes
class	JetHypoAlgType
	— Menu Sequence helpers — More...

Functions
	getCaloInputMaker ()
	— Input Makers -—
	getTrackingInputMaker (AthConfigFlags flags, str trkopt)
	getHitZInputMaker (AthConfigFlags flags)
	jetSelectionCfg (flags, jetDefStr, jetsIn, hypoType=JetHypoAlgType.STANDARD)
	selName (recoSequenceName, hypoType=JetHypoAlgType.STANDARD)
	hypoToolGenerator (hypoType)
	jetCaloPreselMenuSequenceGenCfg (flags, **jetDefDict)
	— Menu Sequence getters —
	jetCaloRecoMenuSequenceGenCfg (flags, clusterCalib)
	jetCaloHypoMenuSequenceGenCfg (flags, isPerf, **jetDefDict)
	jetHICaloHypoMenuSequenceGenCfg (flags, isPerf, **jetRecoDict)
	jetFSTrackingHypoMenuSequenceGenCfg (flags, isPerf, **jetDefDict)
	jetRoITrackJetTagSelCfg (flags, preselJetDef, isPresel=True)
	jetRoITrackJetTagHypoMenuSequenceGenCfg (flags, jetDef, isPresel=True)
	jetRoIHitZRecoCfg (flags, preselJetDef)

Variables
	log = logging.getLogger(__name__)

Function Documentation

getCaloInputMaker()

getCaloInputMaker

(

)

— Input Makers -—

Definition at line 31 of file JetMenuSequencesConfig.py.

def getCaloInputMaker():
    from TrigT2CaloCommon.CaloDef import clusterFSInputMaker
    InputMakerAlg = clusterFSInputMaker()
    return InputMakerAlg
 
# For later steps, where calo reco should not be run
# The same instance of an algorithm cannot be run in different steps
# Used for chains that use tracking

getHitZInputMaker()

getHitZInputMaker

(

AthConfigFlags

flags

)

Definition at line 88 of file JetMenuSequencesConfig.py.

def getHitZInputMaker(flags : AthConfigFlags):
    return CompFactory.EventViewCreatorAlgorithm(
        "IMJetRoIHitZ",
        mergeUsingFeature = False,
        RoITool = CompFactory.ViewCreatorJetSuperROITool(
            'ViewCreatorJetSuperRoIHitZ',
            RoisWriteHandleKey  = recordable( 'HLT_Roi_JetSuperHitZ' ),
            RoIEtaWidth = flags.Trigger.InDetTracking.jetSuper.etaHalfWidth,
            RoIPhiWidth = flags.Trigger.InDetTracking.jetSuper.phiHalfWidth,
            RoIZWidth   = flags.Trigger.InDetTracking.jetSuper.zedHalfWidth,
        ),
        Views = "JetSuperRoIHitZViews",
        InViewRoIs = "InViewRoIs",
        RequireParentView = False,
        ViewFallThrough = True,
    )
 

getTrackingInputMaker()

getTrackingInputMaker	(	AthConfigFlags	flags,
		str	trkopt )

Definition at line 39 of file JetMenuSequencesConfig.py.

def getTrackingInputMaker(flags : AthConfigFlags, trkopt : str):
    if trkopt=="ftf":
        log.debug( "jet FS tracking: useDynamicRoiZWidth: %s", flags.Trigger.InDetTracking.fullScan.useDynamicRoiZWidth )
        
        roiUpdater = None
        if flags.Trigger.InDetTracking.fullScan.useDynamicRoiZWidth:
            roiUpdater = CompFactory.RoiUpdaterTool( useBeamSpot=True )
 
            log.info( roiUpdater )
 
            InputMakerAlg = CompFactory.InputMakerForRoI(
                "IM_Jet_TrackingStep",
                mergeUsingFeature = False,
                RoITool = CompFactory.ViewCreatorFSROITool(
                    name="RoiTool_FS",
                    RoiUpdater=roiUpdater,
                    RoisWriteHandleKey=recordable( flags.Trigger.InDetTracking.fullScan.roi )
                ),
                RoIs = trkFSRoI,
            )
        else: 
            InputMakerAlg = CompFactory.InputMakerForRoI(
                "IM_Jet_TrackingStep",
                mergeUsingFeature = False,
                RoITool = CompFactory.ViewCreatorInitialROITool(),
                RoIs = trkFSRoI,
            )
 
    elif trkopt=="roiftf":
        InputMakerAlg = CompFactory.EventViewCreatorAlgorithm(
            "IMJetRoIFTF",
            mergeUsingFeature = False,
            RoITool = CompFactory.ViewCreatorJetSuperROITool(
                'ViewCreatorJetSuperRoI',
                RoisWriteHandleKey  = recordable( flags.Trigger.InDetTracking.jetSuper.roi ),
                RoIEtaWidth = flags.Trigger.InDetTracking.jetSuper.etaHalfWidth,
                RoIPhiWidth = flags.Trigger.InDetTracking.jetSuper.phiHalfWidth,
                RoIZWidth   = flags.Trigger.InDetTracking.jetSuper.zedHalfWidth,
            ),
            Views = "JetSuperRoIViews",
            InViewRoIs = "InViewRoIs",
            RequireParentView = False,
            ViewFallThrough = True,
        )
    else:
        raise RuntimeError(f"Unrecognised trkopt '{trkopt}' provided, choices are ['ftf','roiftf']")
    return InputMakerAlg
 
# The same super-RoI as the RoI FTF step, in views of its own

hypoToolGenerator()

hypoToolGenerator

(

hypoType

)

returns function (that in turn returns hypo tool) for menu sequence

Definition at line 167 of file JetMenuSequencesConfig.py.

def hypoToolGenerator(hypoType):
    """returns function (that in turn returns hypo tool) for menu sequence"""
    def trigStreamerHypoTool(flags, chainDict):
        return CompFactory.TrigStreamerHypoTool(chainDict["chainName"])
    return {
        JetHypoAlgType.STANDARD:    trigJetHypoToolFromDict,
        JetHypoAlgType.PASSTHROUGH: trigStreamerHypoTool,
        JetHypoAlgType.CALOPRESEL:  caloPreselJetHypoToolFromDict,
        JetHypoAlgType.ROIPRESEL:   roiPreselJetHypoToolFromDict,
    }[hypoType]
    
 
 

jetCaloHypoMenuSequenceGenCfg()

jetCaloHypoMenuSequenceGenCfg	(		flags,
			isPerf,
		**	jetDefDict )

Definition at line 231 of file JetMenuSequencesConfig.py.

def jetCaloHypoMenuSequenceGenCfg(flags, isPerf, **jetDefDict):
    jetsOut, jetDef = jetDefDict['final']
    jetDefStr = jetDefToString(jetDef)
    reco = InEventRecoCA(f"jetSeqCaloHypo_{jetDefStr}{'_perf' if isPerf else ''}_RecoSequence", inputMaker=getCaloInputMaker())
 
    if 'LC' in jetDef.inputdef.label:
        reco.mergeReco(jetmetTopoClusteringCfg_LC(flags, RoIs=''))
    else:
        reco.mergeReco(jetmetTopoClusteringCfg(flags, RoIs=''))
 
    from .JetRecoSequencesConfig import JetRecoCfg
    jetreco = JetRecoCfg(flags, **jetDefDict)
    reco.mergeReco(jetreco)
    log.debug("Generating jet calo hypo menu sequence for reco %s",jetDef.fullname())
 
    hypoType = JetHypoAlgType.PASSTHROUGH if isPerf else JetHypoAlgType.STANDARD
    selAcc = SelectionCA(selName(reco.name, hypoType=hypoType))
    selAcc.mergeReco(reco)
    selAcc.mergeHypo(jetSelectionCfg(flags, jetDefStr=jetDefStr, jetsIn=jetsOut, hypoType=hypoType))
 
    return MenuSequence(flags, selAcc, HypoToolGen=hypoToolGenerator(hypoType))
 
 
# A full hypo selecting only on heavy ion calo jets (step 1)
# Passing isPerf = True disables the hypo
# We set RoIs='' for same reason as described for jetCaloPreselMenuSequence
@AccumulatorCache

jetCaloPreselMenuSequenceGenCfg()

jetCaloPreselMenuSequenceGenCfg	(		flags,
		**	jetDefDict )

— Menu Sequence getters —

Definition at line 188 of file JetMenuSequencesConfig.py.

def jetCaloPreselMenuSequenceGenCfg(flags, **jetDefDict):
    jetsOut, jetDef = jetDefDict['final']
    jetDefStr = jetDefToString(jetDef)
    reco = InEventRecoCA(f"jetSeqCaloPresel_{jetDefStr}_RecoSequence", inputMaker=getCaloInputMaker())
 
    if 'LC' in jetDef.inputdef.label:
        reco.mergeReco(jetmetTopoClusteringCfg_LC(flags, RoIs=''))
    else:
        reco.mergeReco(jetmetTopoClusteringCfg(flags, RoIs=''))
    
    from .JetRecoSequencesConfig import JetRecoCfg
    jetreco = JetRecoCfg(flags, **jetDefDict)
    reco.mergeReco(jetreco)
 
    log.debug("Generating jet preselection menu sequence for reco %s",jetDef.fullname())
    selAcc = SelectionCA(selName(reco.name, hypoType=JetHypoAlgType.CALOPRESEL))
    selAcc.mergeReco(reco)
    selAcc.mergeHypo(jetSelectionCfg(flags, jetDefStr=jetDefStr, jetsIn=jetsOut, hypoType=JetHypoAlgType.CALOPRESEL))
 
    return MenuSequence(flags, selAcc, HypoToolGen=hypoToolGenerator(hypoType=JetHypoAlgType.CALOPRESEL))
 
# A null preselection, which will only run the cluster making (step 1)
# We set RoIs='' for same reason as described for jetCaloPreselMenuSequence
@AccumulatorCache

jetCaloRecoMenuSequenceGenCfg()

jetCaloRecoMenuSequenceGenCfg	(		flags,
			clusterCalib )

Definition at line 212 of file JetMenuSequencesConfig.py.

def jetCaloRecoMenuSequenceGenCfg(flags, clusterCalib):
    reco = InEventRecoCA(f"jetSeqCaloReco_{clusterCalib}_RecoSequence", inputMaker=getCaloInputMaker())
 
    if clusterCalib=='lcw':
        reco.mergeReco(jetmetTopoClusteringCfg_LC(flags, RoIs=''))
    else:
        reco.mergeReco(jetmetTopoClusteringCfg(flags, RoIs=''))
 
    selAcc = SelectionCA(selName(reco.name, hypoType=JetHypoAlgType.PASSTHROUGH))
    selAcc.mergeReco(reco)
    selAcc.mergeHypo(jetSelectionCfg(flags, jetDefStr="caloreco", jetsIn=None, hypoType=JetHypoAlgType.PASSTHROUGH))
 
    return MenuSequence(flags, selAcc, HypoToolGen=hypoToolGenerator(hypoType=JetHypoAlgType.PASSTHROUGH))
 
 
# A full hypo selecting only on calo jets (step 1)
# Passing isPerf = True disables the hypo
# We set RoIs='' for same reason as described for jetCaloPreselMenuSequence
@AccumulatorCache

jetFSTrackingHypoMenuSequenceGenCfg()

jetFSTrackingHypoMenuSequenceGenCfg	(		flags,
			isPerf,
		**	jetDefDict )

Definition at line 278 of file JetMenuSequencesConfig.py.

def jetFSTrackingHypoMenuSequenceGenCfg(flags, isPerf, **jetDefDict):
    jetsOut, jetDef = jetDefDict['final']
    jetDefStr = jetDefToString(jetDef)
    trkopt = jetDef.context
    reco = InEventRecoCA(f"jetFSTrackingHypo_{jetDefStr}{'_perf' if isPerf else ''}_RecoSequence", inputMaker=getTrackingInputMaker(flags,trkopt))
 
    assert trkopt != "notrk"
    from .JetTrackingConfig import JetFSTrackingCfg
    trk_acc = JetFSTrackingCfg(flags, trkopt, trkFSRoI)
    reco.mergeReco(trk_acc)
 
    from .JetRecoSequencesConfig import JetRecoCfg
    jetreco = JetRecoCfg(flags, **jetDefDict)
    reco.mergeReco(jetreco)
    log.debug("Generating jet tracking hypo menu sequence for reco %s",jetDef.fullname())
 
    if 'PFlow' in jetDef.basename and jetDef.basename.startswith('AntiKt4') and 'sub' in jetCalibFromJetDef(jetDef):
        pvKey = flags.Trigger.InDetTracking.fullScan.vertex_jet
        trig_evt_info_key = recordable("HLT_TCEventInfo_jet")
 
        # Can encapsulate in another CA if necessary but only this instance
        # currently needed
        reco.mergeReco(
            TrigEventInfoRecorderAlgCfg(
                flags,
                name="TrigEventInfoRecorderAlg_jet",
                decoratePFlowInfo=True,
                decorateEMTopoInfo=False,
                trigEventInfoKey=trig_evt_info_key, primaryVertexInputName=pvKey,
                RhoKey_EMTopo='HLT_Kt4EMTopoEventShape', RhoKey_PFlow='HLT_Kt4EMPFlowEventShape'
            )
        )
 
    hypoType = JetHypoAlgType.PASSTHROUGH if isPerf else JetHypoAlgType.STANDARD
    selAcc = SelectionCA(selName(reco.name, hypoType=hypoType))
    selAcc.mergeReco(reco)
    selAcc.mergeHypo(jetSelectionCfg(flags, jetDefStr=jetDefStr, jetsIn=jetsOut, hypoType=hypoType))
 
    return MenuSequence(flags, selAcc, HypoToolGen=hypoToolGenerator(hypoType))
 
 
# A full hypo selecting on jets with RoI track reco (step 2)
# Needs to be preceded by a presel sequence, and be provided
# with the input jets from which to define RoIs
# Presel jets to be reused, which makes ghost association impossible
# Substitute DR association decorator
@AccumulatorCache

jetHICaloHypoMenuSequenceGenCfg()

jetHICaloHypoMenuSequenceGenCfg	(		flags,
			isPerf,
		**	jetRecoDict )

Definition at line 258 of file JetMenuSequencesConfig.py.

def jetHICaloHypoMenuSequenceGenCfg(flags, isPerf, **jetRecoDict):
    reco = InEventRecoCA(f"jetSeqHICaloHypo_{jetRecoDict['jetDefStr']}{'_perf' if isPerf else ''}_RecoSequence", inputMaker=getCaloInputMaker())
 
    reco.mergeReco( HICaloTowerCfg(flags) )
 
    from .JetHIConfig import jetHIRecoSequenceCA
    jetreco, jetsOut, jetDef = jetHIRecoSequenceCA(flags, clustersKey="HLT_HICaloClustersFS",towerKey = fs_towers, **jetRecoDict)
    reco.mergeReco(jetreco)        
    log.debug("Generating jet HI calo hypo menu sequence for reco %s",jetDef.fullname())
    hypoType = JetHypoAlgType.PASSTHROUGH if isPerf else JetHypoAlgType.STANDARD
    selAcc = SelectionCA(selName(reco.name, hypoType=hypoType))
    selAcc.mergeReco(reco)
    selAcc.mergeHypo(jetSelectionCfg(flags, jetDefStr=jetRecoDict['jetDefStr'], jetsIn=jetsOut, hypoType=hypoType))
 
    return MenuSequence(flags, selAcc, HypoToolGen=hypoToolGenerator(hypoType))
 
 
# A full hypo selecting on jets with FS track reco (step 2)
# To combine either with a presel or a passthrough sequence
@AccumulatorCache

jetRoIHitZRecoCfg()

jetRoIHitZRecoCfg	(		flags,
			preselJetDef )

Definition at line 384 of file JetMenuSequencesConfig.py.

def jetRoIHitZRecoCfg(flags, preselJetDef):
    jetDefStr = jetDefToString(preselJetDef)+'_roihitz'
    reco = InEventRecoCA(
        f"jetRoIHitZHypo_{jetDefStr}_RecoSequence",
        inputMaker=getHitZInputMaker(flags)
    )
 
    from .JetRecoSequencesConfig import (
        HitZJetCopyAlgCfg, formatFilteredJetsName, JET_DEFAULT_VIEW_PT_MIN_GEV
    )
    from .JetRecoCommon import getDecorList
    hitzjet_acc, hitzJetDef = HitZJetCopyAlgCfg(flags,preselJetDef)
    hitzJetName = recordable(hitzJetDef.fullname())
    reco.mergeReco(hitzjet_acc)
 
    seqname = "JetRoIHitZ_RecoSequence"
    reco.addSequence(parOR(seqname),primary=True)
    from ..Bjet.BjetHitZConfig import hitZRoICfg
    reco.merge(hitZRoICfg(flags, hitzJetName, RoIs=reco.inputMaker().InViewRoIs), seqname)
    reco.inputMaker().ViewNodeName = seqname
 
    # Not JetViewAlgCfg: its cache keys on the JetDefinition, whose hash ignores
    # the suffix, so it would hand back the RoI FTF step's _fastftag view alg
    filtered_jetsIn = formatFilteredJetsName(hitzJetName,JET_DEFAULT_VIEW_PT_MIN_GEV)
    reco.addEventAlgo(CompFactory.JetViewAlg(
        "jetview_"+filtered_jetsIn,
        InputContainer=hitzJetName,
        OutputContainer=filtered_jetsIn,
        PtMin=JET_DEFAULT_VIEW_PT_MIN_GEV*1e3,
        DecorDeps=getDecorList(hitzJetDef),
    ))
    return reco, filtered_jetsIn

jetRoITrackJetTagHypoMenuSequenceGenCfg()

jetRoITrackJetTagHypoMenuSequenceGenCfg	(		flags,
			jetDef,
			isPresel = True )

Definition at line 377 of file JetMenuSequencesConfig.py.

def jetRoITrackJetTagHypoMenuSequenceGenCfg(flags, jetDef, isPresel=True):
    selAcc, hypoType = jetRoITrackJetTagSelCfg(flags, jetDef, isPresel)
    return MenuSequence(flags, selAcc, HypoToolGen=hypoToolGenerator(hypoType))
 
 
 
# HitZ on the super-RoI hits, before any tracking: returns the reco and the decorated jets

jetRoITrackJetTagSelCfg()

jetRoITrackJetTagSelCfg	(		flags,
			preselJetDef,
			isPresel = True )

Definition at line 325 of file JetMenuSequencesConfig.py.

def jetRoITrackJetTagSelCfg(flags, preselJetDef, isPresel=True):
    # Seems odd, but we have to combine event and view execution here
    # where InViewRecoCA will do all in view
    jetDefStr = jetDefToString(preselJetDef)+'_roiftf'
    trkopt = 'roiftf'
    reco = InEventRecoCA(
        f"jetRoITrackJetTagHypo_{jetDefStr}_RecoSequence",
        inputMaker=getTrackingInputMaker(flags,trkopt)
    )
 
    # Decorate EventInfo with online beamspot info in global context -- avoid multiple executes
    from ..Bjet.BjetFlavourTaggingConfig import OnlineBeamspotAugmenterCfg
    reco.mergeReco(OnlineBeamspotAugmenterCfg(flags))
 
    # Add to top-level serial sequence to ensure it is ready for in-view reco
    from .JetRecoSequencesConfig import (
        FastFtaggedJetCopyAlgCfg, JetRoITrackJetTagSequenceCfg, JetViewAlgCfg, formatFilteredJetsName, JET_DEFAULT_VIEW_PT_MIN_GEV
    )
    ftagjet_acc, ftaggedJetDef = FastFtaggedJetCopyAlgCfg(flags,preselJetDef)
    ftaggedJetName = recordable(ftaggedJetDef.fullname())
    reco.mergeReco(ftagjet_acc)
 
    track_acc = JetRoITrackJetTagSequenceCfg(
        flags,
        ftaggedJetName,
        trkopt,
        RoIs=reco.inputMaker().InViewRoIs)
    # Explicitly add the sequence here that is to run in the super-RoI view
    seqname = f"JetRoITrackJetTag_{trkopt}_RecoSequence"
    reco.addSequence(parOR(seqname),primary=True)
    verifier = CompFactory.AthViews.ViewDataVerifier("roiftf_ftag_ViewDataVerifier")
    verifier.DataObjects = [
        ('SG::AuxElement',f'EventInfo.{dec}')
        for dec in ['onlineBeamPosSigmaXY','onlineBeamPosX','onlineBeamPosY','onlineBeamPosZ']
    ]
    reco.addEventAlgo(verifier,seqname)
    reco.merge(track_acc,seqname)
    reco.inputMaker().ViewNodeName = seqname
 
    # Run the JetViewAlg sequence to filter out low pT jets
    # Have to run it outside of JetRoITrackJetTagSequence (which runs in EventView), so that hypo recognises the filtered jets.
    jetview_Acc = JetViewAlgCfg(flags,jetDef=ftaggedJetDef)
    filtered_jetsIn = formatFilteredJetsName(ftaggedJetName,JET_DEFAULT_VIEW_PT_MIN_GEV)
    reco.merge(jetview_Acc)
 
    # Needs track-to-jet association here, maybe with dR decorator
    hypoType = JetHypoAlgType.ROIPRESEL if isPresel else JetHypoAlgType.STANDARD
    selAcc = SelectionCA(selName(reco.name, hypoType=hypoType))
    selAcc.mergeReco(reco)
    selAcc.mergeHypo(jetSelectionCfg(flags, jetDefStr=jetDefStr, jetsIn=filtered_jetsIn, hypoType=hypoType))
    return selAcc, hypoType
 

jetSelectionCfg()

jetSelectionCfg	(		flags,
			jetDefStr,
			jetsIn,
			hypoType = JetHypoAlgType.STANDARD )

constructs CA with hypo alg given arguments 

Definition at line 134 of file JetMenuSequencesConfig.py.

def jetSelectionCfg(flags, jetDefStr, jetsIn, hypoType=JetHypoAlgType.STANDARD):
    """constructs CA with hypo alg given arguments """
    if hypoType==JetHypoAlgType.PASSTHROUGH:
        hyponame = f"TrigStreamerHypoAlg_{jetDefStr}_passthrough"
        hypo = CompFactory.TrigStreamerHypoAlg(hyponame)
    else:
        assert jetsIn is not None
        if hypoType==JetHypoAlgType.CALOPRESEL:
            hyponame = f"TrigJetHypoAlg_{jetDefStr}_calopresel"
            hypo = CompFactory.TrigJetHypoAlg(hyponame, Jets=jetsIn, DoPresel=True)
        elif hypoType==JetHypoAlgType.ROIPRESEL:
            hyponame = f"TrigJetHypoAlg_{jetDefStr}_roipresel"
            hypo = CompFactory.TrigJetHypoAlg(hyponame, Jets=jetsIn, DoPresel=True)
        else:
            hyponame = f"TrigJetHypoAlg_{jetDefStr}"
            hypo = CompFactory.TrigJetHypoAlg(hyponame, Jets=jetsIn)
    ca = ComponentAccumulator()
    ca.addEventAlgo(hypo)
    return ca
 

selName()

selName	(		recoSequenceName,
			hypoType = JetHypoAlgType.STANDARD )

Construct selection (the name passed to SelectionCA) given reco sequence and hypo type

Definition at line 154 of file JetMenuSequencesConfig.py.

def selName(recoSequenceName, hypoType=JetHypoAlgType.STANDARD):
    """Construct selection (the name passed to SelectionCA) given reco sequence and hypo type"""
    selname = recoSequenceName.replace('RecoSequence','MenuSequence')
    if hypoType==JetHypoAlgType.PASSTHROUGH:
        selname += "_passthrough"
    else:
        if hypoType==JetHypoAlgType.CALOPRESEL:
            selname += "_calopresel"
        elif hypoType==JetHypoAlgType.ROIPRESEL:
            selname += "_roipresel"
    return selname
 
 

Variable Documentation

log

python.HLT.Jet.JetMenuSequencesConfig.log = logging.getLogger(__name__)

Definition at line 24 of file JetMenuSequencesConfig.py.