How to make track jets from AOD
From ATLAS-TRIUMF
- Event selector
include( "AthenaPoolCnvSvc/ReadAthenaPool_jobOptions.py" ) EventSelector = Service( "EventSelector" )
- EventSelector.BackNavigation = True
- Particle Properties
include( "PartPropSvc/PartPropSvc.py" )
- the POOL converters
include( "ParticleBuilderOptions/ESD_PoolCnv_jobOptions.py" ) include( "ParticleBuilderOptions/AOD_PoolCnv_jobOptions.py" ) include( "ParticleBuilderOptions/McAOD_PoolCnv_jobOptions.py" ) include( "EventAthenaPool/EventAthenaPool_joboptions.py" )
- The AOD input file
EventSelector.InputCollections = [ "/home/yoshio/carnage/users/yoshio/vbfhww/FakeMet/AODtoNTPL/misal1_csc11.005015.J6_pythia_jetjet.KAOD.v12000301_tid003096._00033.1.root" ]
- The libraries to loaded at run time
theApp.Dlls+=[ "AnalysisTools" ] theApp.Dlls+=[ "UserAnalysisUtils" ] theApp.Dlls+=[ "UserAnalysis" ]
- Athena-Aware NTuple making Tools
CBNTAthenaAware = True include("CBNT_Athena/CBNT_AthenaAware_jobOptions.py" ) include("CBNT_Athena/CBNT_EventInfo_jobOptions.py" )
- include ( "JetRec/ConeTowerJet_jobOptions.py" )
include ( "ConeTowerJet_jobOptions.py" ) ConeTowerJets.InputCollectionNames = [ "TrackParticleCandidate" ] ConeTowerJets.JetCollectionName = "TrackJets"
- include ( "Cone4TopoJet_jobOptions.py" )
- Cone4TopoJetAlg.InputCollectionNames = [ "TrackParticleCandidate" ]
- Cone4TopoJetAlg.JetCollectionName = "TrackJets"
- from(PhysicsAnalysis/JetMissingEtID/JetMissingEtAlgs/share/ParticleJetBuilder_jobOptions.py)
import EventKernel.ParticleDataType
theApp.Dlls += ["JetMissingEtUtils"] theApp.Dlls += ["JetMissingEtAlgs"]
ParticleJetList = [ "TrackJets" ]
for key in ParticleJetList:
algname=key.replace("Jets","ParticleJet")+"Builder"
theApp.TopAlg += ["ParticleJetBuilder/"+algname]
thisAlg = Algorithm( algname )
thisAlg.JetCollection=key
thisAlg.ParticleJetContainer = key.replace("Jet","ParticleJet")
thisAlg.ParticleJetBuilderTool.dataType=EventKernel.ParticleDataType.Full
- thisAlg.ParticleJetBuilderTool.CellCalibratorName="H1WeightToolCSC12"
- list of the algorithms to be executed at run time
CBNT_AthenaAware=Algorithm( "CBNT_AthenaAware" ) CBNT_AthenaAware.Members+=[ "AnalysisSkeleton" ]
- The properties of the AnalysisSkeleton Algorithm
AnalysisSkeleton = Algorithm( "AnalysisSkeleton" ) AnalysisSkeleton.ElectronContainer ="ElectronCollection"
- AnalysisSkeleton.ElectronTruthContainer ="SpclMC"
AnalysisSkeleton.MuonContainer ="MuidMuonCollection"
- AnalysisSkeleton.MuonTruthContainer ="SpclMC"
AnalysisSkeleton.ParticleJetContainer ="ConeTowerParticleJets" AnalysisSkeleton.ParticleJetTruthContainer ="ConeTruthParticleJets" AnalysisSkeleton.MissingEtContainer ="MET_Final" AnalysisSkeleton.MissingEtMuonContainer ="MET_Muon" AnalysisSkeleton.MissingEtBaseContainer ="MET_Base" AnalysisSkeleton.MissingEtTruthContainer ="MET_Truth" AnalysisSkeleton.MissingEtMuonTruthContainer="MET_Truth" AnalysisSkeleton.TrackParticleContainer ="TrackParticleCandidate" AnalysisSkeleton.TrackParticleTruthContainer="TrackParticleTruthCollection" AnalysisSkeleton.TrackParticleJetContainer ="TrackParticleJets" AnalysisSkeleton.OutputLevel =INFO
- setup TTree registration Service
- save ROOT histograms and Tuple
THistSvc=Service( "THistSvc" ) THistSvc.Output=[ "AANT DATAFILE='test.root' OPT='RECREATE'" ] theApp.TopAlg+=[ "AANTupleStream" ] AANTupleStream=Algorithm( "AANTupleStream" ) AANTupleStream.ExtraRefNames=[ "StreamESD","Stream1" ] AANTupleStream.OutputName='test.root' AANTupleStream.WriteInputDataHeader=True AANTupleStream.OutputLevel=WARNING
- Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
MessageSvc = Service( "MessageSvc" ) MessageSvc.OutputLevel=ERROR
- Number of Events to process
theApp.EvtMax=500
- For interactive analysis
- include ("PyAnalysisCore/InitPyAnalysisCore.py")
- Detail time measurement and auditors
- Use auditors
theApp.AuditAlgorithms=True theApp.Dlls+=[ "GaudiAud" ]
theAuditorSvc=AuditorSvc() theAuditorSvc.Auditors+=[ "ChronoAuditor" ]
AthenaPoolCnvSvc=Service( "AthenaPoolCnvSvc" ) AthenaPoolCnvSvc.UseDetailChronoStat=TRUE
- StoreGateSvc=Service( "StoreGateSvc" )
- StoreGateSvc.Dump=True
- MessageSvc.OutputLevel=DEBUG
- print "HERE"
- print theApp.TopAlg
theApp.TopAlg = [ 'JetAlgorithm/ConeTowerJets',
'ParticleJetBuilder/TrackParticleJetBuilder',
'CBNT_AthenaAware',
'AANTupleStream' ]
- theApp.TopAlg = [ 'JetAlgorithm/Cone4TopoJetAlg',
- 'ParticleJetBuilder/TrackParticleJetBuilder',
- 'CBNT_AthenaAware',
- 'AANTupleStream' ]

