def minbiasSpacePointAthSequence(ConfigFlags):
InputMakerAlg = clusterFSInputMaker()
(recoSequence, sequenceOut) = MinBiasRecoSequences()
minbiasSpacePointAthSequence = seqAND("minbiasSpacePointAthSequence",
[InputMakerAlg, recoSequence])
return (minbiasSpacePointAthSequence, InputMakerAlg, sequenceOut)
def inputMaker(self):
""" The input maker for this configuration """
if hasattr(self, "_inputMaker"):
return self._inputMaker
from TrigT2CaloCommon.CaloDef import clusterFSInputMaker
self._inputMaker = clusterFSInputMaker()
return self._inputMaker
def jetAthSequence(dummyFlags, **jetRecoDict):
from TrigT2CaloCommon.CaloDef import clusterFSInputMaker
InputMakerAlg = conf2toConfigurable(clusterFSInputMaker())
# Does nothing now, but may need to use this to toggle Insitu calib step
# in which case it should be deduced from input flags
dataSource = "data" # or mc
(recoSequence,
sequenceOut) = RecoFragmentsPool.retrieve(jetRecoSequence,
None,
dataSource=dataSource,
**jetRecoDict)
jetDefString = jetRecoDictToString(jetRecoDict)
JetAthSequence = seqAND("jetAthSequence_" + jetDefString,
[InputMakerAlg, recoSequence])
return (JetAthSequence, InputMakerAlg, sequenceOut)
def generateChains(flags, chainDict):
stepName = getChainStepName('Jet', 1)
stepReco, stepView = createStepView(stepName)
acc = ComponentAccumulator()
acc.addSequence(stepView)
# All this should be some common FS cell module?
from TrigT2CaloCommon.TrigCaloDataAccessConfig import trigCaloDataAccessSvcCfg
acc.merge(trigCaloDataAccessSvcCfg(flags))
cdaSvc = acc.getService(
"TrigCaloDataAccessSvc") # should be made primary component
acc.printConfig()
from TrigT2CaloCommon.CaloDef import clusterFSInputMaker
inEventReco = InEventReco("JetReco", inputMaker=clusterFSInputMaker())
cellsname = "CaloCellsFS"
clustersname = "HLT_CaloTopoClustersFS"
cellmakerCfg = HLTCaloCellMakerCfg(cellsname, cdaSvc)
inEventReco.mergeReco(cellmakerCfg)
from CaloRec.CaloTopoClusterConfig import CaloTopoClusterCfg
inEventReco.mergeReco(
CaloTopoClusterCfg(flags,
cellsname=cellsname,
clustersname=clustersname,
doLCCalib=False,
sequenceName=inEventReco.recoSeq.name))
#sequencing of actual jet reconstruction
from JetRecConfig import JetRecConfig
from JetRecConfig.JetDefinition import JetConstit, JetDefinition, xAODType
#hardcoded jet collection for now
clustermods = ["ECPSFrac", "ClusterMoments"]
trigMinPt = 7e3
HLT_EMTopo = JetConstit(xAODType.CaloCluster, ["EM"])
HLT_EMTopo.rawname = clustersname
HLT_EMTopo.inputname = clustersname
HLT_AntiKt4EMTopo_subjesIS = JetDefinition("AntiKt",
0.4,
HLT_EMTopo,
ptmin=trigMinPt,
ptminfilter=trigMinPt)
HLT_AntiKt4EMTopo_subjesIS.modifiers = [
"Calib:TrigRun2:data:JetArea_EtaJES_GSC_Insitu:HLT_Kt4EMTopoEventShape",
"Sort"
] + clustermods
jetprefix = "HLT_"
jetsuffix = "_subjesIS"
evsprefix = "HLT_"
# May need a switch to disable automatic modifier prerequisite generation
jetRecoComps = JetRecConfig.JetRecCfg(HLT_AntiKt4EMTopo_subjesIS, flags,
jetprefix, jetsuffix, evsprefix)
inEventReco.mergeReco(jetRecoComps)
acc.merge(inEventReco, stepReco.getName())
#hypo
from TrigHLTJetHypo.TrigJetHypoToolConfig import trigJetHypoToolFromDict
hypo = CompFactory.TrigJetHypoAlgMT("TrigJetHypoAlgMT_a4tcem_subjesIS")
jetsfullname = jetprefix + HLT_AntiKt4EMTopo_subjesIS.basename + "Jets" + jetsuffix
hypo.Jets = jetsfullname
acc.addEventAlgo(hypo)
jetSequence = CAMenuSequence(Sequence=inEventReco.sequence(),
Maker=inEventReco.inputMaker(),
Hypo=hypo,
HypoToolGen=trigJetHypoToolFromDict,
CA=acc)
jetStep = ChainStep(name=stepName,
Sequences=[jetSequence],
chainDicts=[chainDict])
l1Thresholds = []
for part in chainDict['chainParts']:
l1Thresholds.append(part['L1threshold'])
log.debug('dictionary is: %s\n', pprint.pformat(chainDict))
acc.printConfig()
chain = Chain(chainDict['chainName'],
L1Thresholds=l1Thresholds,
ChainSteps=[jetStep])
return chain