def JetAlgFromTools(rtools, suffix="HI", persistify=True):
#insert exe tools at front of list, e.g. tracksel and tvassoc for HI etc.
HIJet_exe_tools = []
from JetRec.JetRecFlags import jetFlags
if jetFlags.useTruth() and not jetFlags.Enabled():
HIJet_exe_tools += HITruthParticleCopy()
#if jetFlags.useCells(): HIJet_exe_tools += [jtm.missingcells]
if HIJetFlags.UseHITracks():
HIJet_exe_tools += [jtm.tracksel_HI, jtm.gtracksel_HI, jtm.tvassoc_HI]
rtools = HIJet_exe_tools + rtools
from JetRec.JetRecConf import JetToolRunner
runner = JetToolRunner("jetrun" + suffix,
Tools=rtools,
Timer=jetFlags.timeJetToolRunner())
jtm.add(runner)
from JetRec.JetRecConf import JetAlgorithm
theAlg = JetAlgorithm("jetalg" + suffix)
theAlg.Tools = [runner]
from AthenaCommon.AlgSequence import AlgSequence
topsequence = AlgSequence()
topsequence += theAlg
from GaudiKernel.Constants import DEBUG
if jetFlags.debug > 0:
jtm.setOutputLevel(runner, DEBUG)
theAlg.OutputLevel = DEBUG
if jetFlags.debug > 3:
jtm.setOutputLevel(jtm.jetBuilderWithoutArea, DEBUG)
if persistify:
for t in rtools:
if hasattr(t, "OutputContainer"):
AddToOutputList(t.OutputContainer)
return theAlg
ParticleCaloExtensionTool=ParticleToCaloExtrapolationTool,
CaloClusterLocation="TimedCaloCalTopoClusters",
ClustersInConeTool=CaloClustersInCone,
ConeSize=0.1,
UseCovariance=True)
ParticleCaloCellAssociation.OutputLevel = DEBUG
ToolSvc += ParticleCaloCellAssociation
print ParticleCaloCellAssociation
# Access the algorithm sequence:
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
# Add the algorithm.
jetalg = JetAlgorithm("JetAlg_OriginAndTime")
jetalg.OutputLevel = INFO
jetalg.Tools += [PFSequence]
topSequence += jetalg
from TrackParticleAssociationAlgs.TrackParticleAssociationAlgsConf import TrackParticleClusterAssociationAlg
TrackParticleClusterAssociation = TrackParticleClusterAssociationAlg(
name="TrackParticleClusterAssociationInDet",
ParticleCaloClusterAssociationTool=ParticleCaloCellAssociation,
TrackParticleContainerName="InDetTrackParticles",
PtCut=400.,
OutputCollectionPostFix="Test",
CaloClusterLocation="TimedCaloCalTopoClusters")
# TrackParticleClusterAssociation.OutputLevel = DEBUG
topSequence += TrackParticleClusterAssociation
print TrackParticleClusterAssociation
