ATLAS Offline Software
Loading...
Searching...
No Matches
python.PHYS Namespace Reference

Functions

 PHYSKernelCfg (flags, name='PHYSKernel', **kwargs)
 PHYSCoreCfg (flags, name_tag='PHYS', StreamName='StreamDAOD_PHYS', TriggerListsHelper=None, addExtraVariables=None)
 PHYSCfg (flags)

Variables

 logPHYS = logging.getLogger('PHYS')

Function Documentation

◆ PHYSCfg()

python.PHYS.PHYSCfg ( flags)

Definition at line 265 of file PHYS.py.

265def PHYSCfg(flags):
266
267 logPHYS.info('****************** STARTING PHYS *****************')
268
269 stream_name = 'StreamDAOD_PHYS'
270 acc = ComponentAccumulator()
271
272 # Get the lists of triggers needed for trigger matching.
273 # This is needed at this scope (for the slimming) and further down in the config chain
274 # for actually configuring the matching, so we create it here and pass it down
275 # TODO: this should ideally be called higher up to avoid it being run multiple times in a train
276 from DerivationFrameworkPhys.TriggerListsHelper import TriggerListsHelper
277 PHYSTriggerListsHelper = TriggerListsHelper(flags)
278
279 # Common augmentations
280 acc.merge(PHYSKernelCfg(
281 flags,
282 name="PHYSKernel",
283 StreamName = stream_name,
284 TriggerListsHelper = PHYSTriggerListsHelper
285 ))
286 # PHYS content
287 acc.merge(PHYSCoreCfg(
288 flags,
289 "PHYS",
290 StreamName = stream_name,
291 TriggerListsHelper = PHYSTriggerListsHelper
292 ))
293
294 return acc

◆ PHYSCoreCfg()

python.PHYS.PHYSCoreCfg ( flags,
name_tag = 'PHYS',
StreamName = 'StreamDAOD_PHYS',
TriggerListsHelper = None,
addExtraVariables = None )

Definition at line 63 of file PHYS.py.

63def PHYSCoreCfg(flags, name_tag='PHYS', StreamName='StreamDAOD_PHYS', TriggerListsHelper=None, addExtraVariables=None):
64
65 if TriggerListsHelper is None:
66 from DerivationFrameworkPhys.TriggerListsHelper import TriggerListsHelper
67 TriggerListsHelper = TriggerListsHelper(flags)
68
69 acc = ComponentAccumulator()
70
71
72 from DerivationFrameworkHiggs.HiggsPhysContent import HiggsAugmentationAlgsCfg
73 acc.merge(HiggsAugmentationAlgsCfg(flags))
74
75
77 from IsolationSelection.IsolationSelectionConfig import IsoCloseByAlgsCfg
78 acc.merge(IsoCloseByAlgsCfg(flags, isPhysLite = False, stream_name = StreamName))
79
80
81 from LeptonTaggers.LeptonTaggersConfig import DecoratePLITAlgsCfg
82 acc.merge(DecoratePLITAlgsCfg(flags))
83
84 #===================================================
85 # HEAVY FLAVOR CLASSIFICATION FOR ttbar+jets EVENTS
86 #===================================================
87 from DerivationFrameworkMCTruth.HFClassificationCommonConfig import HFClassificationCommonCfg
88 acc.merge(HFClassificationCommonCfg(flags))
89
90 # ============================
91 # Define contents of the format
92 # =============================
93 from OutputStreamAthenaPool.OutputStreamConfig import OutputStreamCfg
94 from xAODMetaDataCnv.InfileMetaDataConfig import SetupMetaDataForStreamCfg
95 from DerivationFrameworkCore.SlimmingHelper import SlimmingHelper
96
97 PHYSSlimmingHelper = SlimmingHelper(name_tag+"SlimmingHelper", NamesAndTypes = flags.Input.TypedCollections, flags = flags)
98 PHYSSlimmingHelper.SmartCollections = ["EventInfo",
99 "Electrons",
100 "Photons",
101 "Muons",
102 "PrimaryVertices",
103 "InDetTrackParticles",
104 "AntiKt4EMTopoJets",
105 "AntiKt4EMPFlowJets",
106 "MET_Baseline_AntiKt4EMTopo",
107 "MET_Baseline_AntiKt4EMPFlow",
108 "TauJets",
109 "TauJets_MuonRM",
110 "DiTauJets",
111 "DiTauJetsLowPt",
112 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets",
113 "AntiKtVR30Rmax4Rmin02PV0TrackJets",
114 ]
115 if flags.Tau.TauEleRM_isAvailable:
116 PHYSSlimmingHelper.SmartCollections.append("TauJets_EleRM")
117
118 excludedVertexAuxData = "-vxTrackAtVertex.-MvfFitInfo.-isInitialized.-VTAV"
119 StaticContent = []
120 StaticContent += ["xAOD::VertexContainer#SoftBVrtClusterTool_Tight_Vertices"]
121 StaticContent += ["xAOD::VertexAuxContainer#SoftBVrtClusterTool_Tight_VerticesAux." + excludedVertexAuxData]
122 StaticContent += ["xAOD::VertexContainer#SoftBVrtClusterTool_Medium_Vertices"]
123 StaticContent += ["xAOD::VertexAuxContainer#SoftBVrtClusterTool_Medium_VerticesAux." + excludedVertexAuxData]
124 StaticContent += ["xAOD::VertexContainer#SoftBVrtClusterTool_Loose_Vertices"]
125 StaticContent += ["xAOD::VertexAuxContainer#SoftBVrtClusterTool_Loose_VerticesAux." + excludedVertexAuxData]
126 StaticContent += ["xAOD::VertexContainer#NVSI_SecVrt_Tight"]
127 StaticContent += ["xAOD::VertexAuxContainer#NVSI_SecVrt_TightAux."+excludedVertexAuxData]
128
129 PHYSSlimmingHelper.StaticContent = StaticContent
130
131 # Extra content
132 PHYSSlimmingHelper.ExtraVariables += ["AntiKt4EMTopoJets.ConeExclBHadronsFinal.ConeExclCHadronsFinal.GhostBHadronsFinal.GhostCHadronsFinal.GhostBHadronsFinalCount.GhostBHadronsFinalPt.GhostCHadronsFinalCount.GhostCHadronsFinalPt.IsoFixedCone5PtPUsub",
133 "AntiKt4EMPFlowJets.ConeExclBHadronsFinal.ConeExclCHadronsFinal.GhostBHadronsFinal.GhostCHadronsFinal.GhostBHadronsFinalCount.GhostBHadronsFinalPt.GhostCHadronsFinalCount.GhostCHadronsFinalPt.isJvtHS.isJvtPU.IsoFixedCone5PtPUsub",
134 "TruthPrimaryVertices.t.x.y.z",
135 "InDetTrackParticles.TTVA_AMVFVertices.TTVA_AMVFWeights.eProbabilityHT.numberOfTRTHits.numberOfTRTOutliers",
136 "EventInfo.GenFiltHT.GenFiltMET.GenFiltHTinclNu.GenFiltPTZ.GenFiltFatJ.HF_Classification.HF_SimpleClassification.HF_ClassificationC5J20.HF_ClassificationC5J25.HF_ClassificationC15J20.HF_ClassificationC15J25",
137 "TauJets.dRmax.etOverPtLeadTrk",
138 "TauJets_MuonRM.dRmax.etOverPtLeadTrk",
139 "HLT_xAOD__TrigMissingETContainer_TrigEFMissingET.ex.ey",
140 "HLT_xAOD__TrigMissingETContainer_TrigEFMissingET_mht.ex.ey",
141 "HLT_AnomDet_ComboHypo.adScore"]
142
143 if addExtraVariables:
144 PHYSSlimmingHelper.ExtraVariables += addExtraVariables
145
146 # HSGNN Score
147 from TrkConfig.VertexFindingFlags import VertexSortingSetup
148 if flags.Tracking.PriVertex.sortingSetup is VertexSortingSetup.GNNSorting:
149 PHYSSlimmingHelper.ExtraVariables += ["PrimaryVertices.gnnScore"]
150
151 if flags.Tau.TauEleRM_isAvailable:
152 PHYSSlimmingHelper.ExtraVariables += ["TauJets_EleRM.dRmax.etOverPtLeadTrk"]
153
154 # IFF extra content
155 from LeptonTaggers.LeptonTaggersConfig import GetExtraPLITVariablesForDxAOD
156 PHYSSlimmingHelper.ExtraVariables += GetExtraPLITVariablesForDxAOD()
157
158 # Truth extra content
159 if flags.Input.isMC:
160
161 from DerivationFrameworkMCTruth.MCTruthCommonConfig import addTruth3ContentToSlimmerTool
162 addTruth3ContentToSlimmerTool(PHYSSlimmingHelper)
163 PHYSSlimmingHelper.AllVariables += ['TruthLHEParticles','InTimeAntiKt4TruthJets','OutOfTimeAntiKt4TruthJets']
164 PHYSSlimmingHelper.ExtraVariables += ["Electrons.TruthLink",
165 "Muons.TruthLink",
166 "Photons.TruthLink",
167 "AntiKt4TruthDressedWZJets.IsoFixedCone5Pt.HFHadronOriginID",
168 "TruthHFWithDecayParticles.prodVtxLink.prodVtxLink.prodVtxLink.decayVtxLink.decayVtxLink.decayVtxLink.parentLinks.childLinks.m.px.py.pz.e.pdgId.Classification.uid.classifierParticleOrigin.classifierParticleType.classifierParticleOutCome.status",
169 "TruthHFWithDecayVertices.incomingParticleLinks.outgoingParticleLinks.uid.status.x.y.z.t",
170 "TruthCharm.prodVtxLink.prodVtxLink.prodVtxLink.decayVtxLink.decayVtxLink.decayVtxLink.parentLinks.childLinks.m.px.py.pz.e.pdgId.Classification.uid.classifierParticleOrigin.classifierParticleType.classifierParticleOutCome.status.barcode.polarizationPhi.polarizationTheta",
171 "TruthPileupParticles.prodVtxLink.prodVtxLink.prodVtxLink.decayVtxLink.decayVtxLink.decayVtxLink.m.px.py.pz.e.pdgId.Classification.uid.classifierParticleOrigin.classifierParticleType.classifierParticleOutCome.status.barcode.PVz.pileupEventNumber.parentHadronID"]
172
173 from DerivationFrameworkMCTruth.MCTruthCommonConfig import AddTauAndDownstreamParticlesCfg
174 acc.merge(AddTauAndDownstreamParticlesCfg(flags))
175 PHYSSlimmingHelper.ExtraVariables += ["TruthTausWithDecayParticles.prodVtxLink.prodVtxLink.prodVtxLink.decayVtxLink.decayVtxLink.decayVtxLink.m.px.py.pz.e.pdgId.Classification.uid.classifierParticleOrigin.classifierParticleType.classifierParticleOutCome.status",
176 "TruthTausWithDecayVertices.incomingParticleLinks.outgoingParticleLinks.uid.status.x.y.z.t"]
177
178 #Truth jet matching variables - largeR
179 PHYSSlimmingHelper.ExtraVariables += ["AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10TruthLabel_R22v1_MatchedTruthJetMass",
180 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10TruthLabel_R22v1_MatchedTruthJetPt",
181 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10TruthLabel_R22v1_MatchedTruthJetEta",
182 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10TruthLabel_R22v1_MatchedTruthJetPhi",
183 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10TruthLabel_R22v1_MatchedTruthJetDR",
184 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10TruthLabel_R22v1_MatchedTruthJetGF",
185 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10TruthLabel_R22v1_MatchedTruthGroomedJetMass",
186 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10TruthLabel_R22v1_MatchedTruthGroomedJetPt",
187 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10WZTruthLabel_R22v1_MatchedTruthJetMass",
188 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10WZTruthLabel_R22v1_MatchedTruthJetPt",
189 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10WZTruthLabel_R22v1_MatchedTruthJetEta",
190 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10WZTruthLabel_R22v1_MatchedTruthJetPhi",
191 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10WZTruthLabel_MatchedTruthJetDR",
192 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10WZTruthLabel_MatchedTruthJetGF",
193 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10WZTruthLabel_R22v1_MatchedTruthGroomedJetMass",
194 "AntiKt10UFOCSSKSoftDropBeta100Zcut10Jets.R10WZTruthLabel_R22v1_MatchedTruthGroomedJetPt"]
195
196 #Truth jet matching variables - smallR
197 PHYSSlimmingHelper.ExtraVariables += ["AntiKt4EMPFlowJets.R4TruthLabel_MatchedTruthJetMass",
198 "AntiKt4EMPFlowJets.R4TruthLabel_MatchedTruthJetPt",
199 "AntiKt4EMPFlowJets.R4TruthLabel_MatchedTruthJetEta",
200 "AntiKt4EMPFlowJets.R4TruthLabel_MatchedTruthJetPhi",
201 "AntiKt4EMPFlowJets.R4TruthLabel_MatchedTruthJetDR",
202 "AntiKt4EMPFlowJets.R4TruthLabel_MatchedTruthJetGF"]
203
204
205 from DerivationFrameworkHiggs.HiggsPhysContent import setupHiggsSlimmingVariables
206 setupHiggsSlimmingVariables(flags, PHYSSlimmingHelper)
207
208
209 PHYSSlimmingHelper.AllVariables += [ 'AFPSiHitContainer', 'AFPToFHitContainer' ]
210
211
212 PHYSSlimmingHelper.AppendToDictionary.update({'MET_Core_AntiKt4EMPFlowHR':'xAOD::MissingETContainer', 'MET_Core_AntiKt4EMPFlowHRAux':'xAOD::MissingETAuxContainer',
213 'METAssoc_AntiKt4EMPFlowHR':'xAOD::MissingETAssociationMap', 'METAssoc_AntiKt4EMPFlowHRAux':'xAOD::MissingETAuxAssociationMap'})
214
215 PHYSSlimmingHelper.AllVariables += ['METAssoc_AntiKt4EMPFlowHR']
216
217 PHYSSlimmingHelper.ExtraVariables += ['Muons.UEcorr_Pt','Electrons.UEcorr_Pt','MET_Core_AntiKt4EMPFlowHR.name.mpx.mpy.sumet.source']
218
219 # Trigger content
220 PHYSSlimmingHelper.IncludeTriggerNavigation = False
221 PHYSSlimmingHelper.IncludeJetTriggerContent = False
222 PHYSSlimmingHelper.IncludeMuonTriggerContent = False
223 PHYSSlimmingHelper.IncludeEGammaTriggerContent = False
224 PHYSSlimmingHelper.IncludeTauTriggerContent = False
225 PHYSSlimmingHelper.IncludeEtMissTriggerContent = False
226 PHYSSlimmingHelper.IncludeBJetTriggerContent = False
227 PHYSSlimmingHelper.IncludeBPhysTriggerContent = False
228 PHYSSlimmingHelper.IncludeMinBiasTriggerContent = False
229 # Compact b-jet trigger matching info
230 PHYSSlimmingHelper.IncludeBJetTriggerByYearContent = True
231
232 # Trigger matching
233 # Run 2
234 if flags.Trigger.EDMVersion == 2:
235 from DerivationFrameworkPhys.TriggerMatchingCommonConfig import AddRun2TriggerMatchingToSlimmingHelper
236 AddRun2TriggerMatchingToSlimmingHelper(SlimmingHelper = PHYSSlimmingHelper,
237 OutputContainerPrefix = "TrigMatch_",
238 TriggerList = TriggerListsHelper.Run2TriggerNamesTau)
239 AddRun2TriggerMatchingToSlimmingHelper(SlimmingHelper = PHYSSlimmingHelper,
240 OutputContainerPrefix = "TrigMatch_",
241 TriggerList = TriggerListsHelper.Run2TriggerNamesNoTau)
242 # Run 3, or Run 2 with navigation conversion
243 if flags.Trigger.EDMVersion == 3 or (flags.Trigger.EDMVersion == 2 and flags.Trigger.doEDMVersionConversion):
244 from TrigNavSlimmingMT.TrigNavSlimmingMTConfig import AddRun3TrigNavSlimmingCollectionsToSlimmingHelper
245 AddRun3TrigNavSlimmingCollectionsToSlimmingHelper(PHYSSlimmingHelper)
246
247 # L1 trigger objects
248 from Campaigns.Utils import getDataYear
249 if getDataYear(flags) >= 2024:
250 # Run 3 with Phase I jet RoIs.
251 from DerivationFrameworkPhys.TriggerMatchingCommonConfig import AddjFexRoIsToSlimmingHelper
252 AddjFexRoIsToSlimmingHelper(SlimmingHelper = PHYSSlimmingHelper)
253 elif getDataYear(flags) >= 2015:
254 # Run 2 and early Run 3, legacy L1 RoIs
255 from DerivationFrameworkPhys.TriggerMatchingCommonConfig import AddLegacyL1JetRoIsToSlimmingHelper
256 AddLegacyL1JetRoIsToSlimmingHelper(SlimmingHelper = PHYSSlimmingHelper)
257
258 # Output stream
259 PHYSItemList = PHYSSlimmingHelper.GetItemList()
260 acc.merge(OutputStreamCfg(flags, "DAOD_"+name_tag, ItemList=PHYSItemList, AcceptAlgs=[name_tag+"Kernel"]))
261 acc.merge(SetupMetaDataForStreamCfg(flags, "DAOD_"+name_tag, AcceptAlgs=[name_tag+"Kernel"], createMetadata=[MetadataCategory.CutFlowMetaData, MetadataCategory.TruthMetaData]))
262
263 return acc
264

◆ PHYSKernelCfg()

python.PHYS.PHYSKernelCfg ( flags,
name = 'PHYSKernel',
** kwargs )
Configure the derivation framework driving algorithm (kernel) for PHYS

Definition at line 17 of file PHYS.py.

17def PHYSKernelCfg(flags, name='PHYSKernel', **kwargs):
18 """Configure the derivation framework driving algorithm (kernel) for PHYS"""
19 acc = ComponentAccumulator()
20
21 from TrkConfig.VertexFindingFlags import VertexSortingSetup
22 if flags.Tracking.PriVertex.sortingSetup is VertexSortingSetup.GNNSorting:
23 from DerivationFrameworkPhys.GNNVertexConfig import GNNVertexCfg
24 acc.merge(GNNVertexCfg(flags))
25
26 # Common augmentations
27 from DerivationFrameworkPhys.PhysCommonConfig import PhysCommonAugmentationsCfg
28 acc.merge(PhysCommonAugmentationsCfg(
29 flags,
30 TriggerListsHelper = kwargs['TriggerListsHelper']
31 ))
32
33 # Thinning tools
34 # These are set up in PhysCommonThinningConfig. Only thing needed here the list of tools to schedule
35 nametag = name.replace('Kernel', '') #get the name to label the tools below such that other formats can use this KernelCfg
36 thinningToolsArgs = {
37 'TrackParticleThinningToolName' : nametag+"TrackParticleThinningTool",
38 'MuonTPThinningToolName' : nametag+"MuonTPThinningTool",
39 'TauJetThinningToolName' : nametag+"TauJetThinningTool",
40 'TauJets_MuonRMThinningToolName' : nametag+"TauJets_MuonRMThinningTool",
41 'DiTauThinningToolName' : nametag+"DiTauThinningTool",
42 'DiTauTPThinningToolName' : nametag+"DiTauTPThinningTool",
43 'DiTauLowPtThinningToolName' : nametag+"DiTauLowPtThinningTool",
44 'DiTauLowPtTPThinningToolName' : nametag+"DiTauLowPtTPThinningTool",
45 }
46 # for AOD produced before 24.0.17, the electron removal tau is not available
47 if flags.Tau.TauEleRM_isAvailable:
48 thinningToolsArgs['TauJets_EleRMThinningToolName'] = nametag+"TauJets_EleRMThinningTool"
49 # Configure the thinning tools
50 from DerivationFrameworkPhys.PhysCommonThinningConfig import PhysCommonThinningCfg
51 acc.merge(PhysCommonThinningCfg(flags, StreamName = kwargs['StreamName'], **thinningToolsArgs))
52 # Get them from the CA so they can be added to the kernel
53 thinningTools = []
54 for key in thinningToolsArgs:
55 thinningTools.append(acc.getPublicTool(thinningToolsArgs[key]))
56
57 # The kernel algorithm itself
58 DerivationKernel = CompFactory.DerivationFramework.DerivationKernel
59 acc.addEventAlgo(DerivationKernel(name, ThinningTools = thinningTools))
60 return acc
61
62

Variable Documentation

◆ logPHYS

python.PHYS.logPHYS = logging.getLogger('PHYS')

Definition at line 14 of file PHYS.py.