ATLAS Offline Software
Loading...
Searching...
No Matches
TrigTauRecToolsConfig.py
Go to the documentation of this file.
1# Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
2
3from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
4from AthenaConfiguration.ComponentFactory import CompFactory
5from AthenaCommon.SystemOfUnits import mm
6
7def trigTauVertexFinderCfg(flags, name=''):
8 '''Algorithm that overwrites numTrack() and charge() of tauJets in container'''
9 acc = ComponentAccumulator()
10
11 acc.setPrivateTools(CompFactory.TauVertexFinder(
12 name = name,
13 UseTJVA = False,
14 AssociatedTracks = 'GhostTrack',
15 InDetTrackSelectionToolForTJVA = '',
16 Key_vertexInputContainer = '',
17 TVATool = '',
18 ))
19
20 return acc
21
22def trigTauTrackFinderCfg(flags, name='', TrackParticlesContainer=''):
23 '''Tau track association'''
24 acc = ComponentAccumulator()
25
26 from TrkConfig.TrkVertexFitterUtilsConfig import AtlasTrackToVertexIPEstimatorCfg
27 AtlasTrackToVertexIPEstimator = acc.popToolsAndMerge(AtlasTrackToVertexIPEstimatorCfg(flags))
28
29 from TrackToVertex.TrackToVertexConfig import TrackToVertexCfg
30 TrackToVertexTool = acc.popToolsAndMerge(TrackToVertexCfg(flags))
31
32 from TrackToCalo.TrackToCaloConfig import ParticleCaloExtensionToolCfg
33 ParticleCaloExtensionTool = acc.popToolsAndMerge(ParticleCaloExtensionToolCfg(flags))
34
35 from InDetConfig.InDetTrackSelectorToolConfig import TrigTauInDetTrackSelectorToolCfg
36 TrigTauInDetTrackSelectorTool = acc.popToolsAndMerge(TrigTauInDetTrackSelectorToolCfg(flags))
37
38 acc.setPrivateTools(CompFactory.TauTrackFinder(
39 name = name,
40 MaxJetDrTau = 0.2,
41 MaxJetDrWide = 0.4,
42 TrackSelectorToolTau = TrigTauInDetTrackSelectorTool,
43 TrackToVertexTool = TrackToVertexTool,
44 Key_trackPartInputContainer = TrackParticlesContainer,
45 maxDeltaZ0wrtLeadTrk = 0.75*mm,
46 removeTracksOutsideZ0wrtLeadTrk = True,
47 ParticleCaloExtensionTool = ParticleCaloExtensionTool,
48 BypassExtrapolator = True,
49 tauParticleCache = "",
50 TrackToVertexIPEstimator = AtlasTrackToVertexIPEstimator,
51 ))
52
53 return acc
54
55def tauVertexVariablesCfg(flags, name=''):
56 '''Vertex variables calculation'''
57 acc = ComponentAccumulator()
58
59 from TrkConfig.TrkVertexFittersConfig import TauAdaptiveVertexFitterCfg
60 TauAdaptiveVertexFitter = acc.popToolsAndMerge(TauAdaptiveVertexFitterCfg(flags))
61
62 from TrkConfig.TrkVertexSeedFinderToolsConfig import CrossDistancesSeedFinderCfg
63 CrossDistancesSeedFinder = acc.popToolsAndMerge(CrossDistancesSeedFinderCfg(flags))
64
65 acc.setPrivateTools(CompFactory.TauVertexVariables(
66 name = name,
67 VertexFitter = TauAdaptiveVertexFitter,
68 SeedFinder = CrossDistancesSeedFinder,
69 ))
70
71 return acc
72
73def trigTauJetLVNNEvaluatorCfg(flags, tau_id='', use_taujet_rnnscore=True):
74 '''TauJet identification inference based on LVNN models, for RNNs and DeepSets'''
75 acc = ComponentAccumulator()
76
77 try: id_flags = getattr(flags.Trigger.Offline.Tau, tau_id)
78 except NameError: raise ValueError(f'Invalid LVNN (RNN, DeepSet) TauID configuration: {tau_id}')
79
80 # For legacy mediumRNN/tightRNN_tracktwoMVA/tracktwoLLP/trackLRT chains
81 output_variable = 'RNNJetScore' if use_taujet_rnnscore else f'{tau_id}_Score'
82
83 # This is to prevent clashes between the same ID running over the legacy chains, and the new decorator-based chains
84 sfx = '_RNNJetScore' if use_taujet_rnnscore else ''
85
86 acc.setPrivateTools(CompFactory.TauJetRNNEvaluator(
87 name = f'TrigTau_TauJetLVNNEvaluator_{tau_id}{sfx}',
88 useTRT = flags.Detector.EnableTRT,
89
90 # Network config:
91 NetworkFile0P = id_flags.NetworkConfig[0],
92 NetworkFile1P = id_flags.NetworkConfig[1],
93 NetworkFile3P = id_flags.NetworkConfig[2],
94 InputLayerScalar = 'scalar',
95 InputLayerTracks = 'tracks',
96 InputLayerClusters = 'clusters',
97 OutputLayer = 'rnnid_output',
98 OutputNode = 'sig_prob',
99
100 # Inputs:
101 MaxTracks = id_flags.MaxTracks,
102 MaxClusters = id_flags.MaxClusters,
103 MaxClusterDR = 1.0,
104 VertexCorrection = False,
105 TrackClassification = False,
106
107 # Decorated TauJet variable names:
108 OutputVarname = output_variable,
109 ))
110
111 return acc
112
113def trigTauJetONNXEvaluatorCfg(flags, tau_id=''):
114 '''TauJet identification inference based on ONNX models, for GNNs, transformers, etc...'''
115 acc = ComponentAccumulator()
116
117 try: id_flags = getattr(flags.Trigger.Offline.Tau, tau_id)
118 except NameError: raise ValueError(f'Invalid ONNX TauID configuration: {tau_id}')
119
120 if isinstance(id_flags.ONNXConfig, str):
121 network_config = {'NetworkFileInclusive': id_flags.ONNXConfig}
122 elif isinstance(id_flags.ONNXConfig, (list, tuple)) and len(id_flags.ONNXConfig) == 3:
123 network_config = {
124 'NetworkFile0P': id_flags.ONNXConfig[0],
125 'NetworkFile1P': id_flags.ONNXConfig[1],
126 'NetworkFile3P': id_flags.ONNXConfig[2],
127 }
128 else:
129 raise ValueError(f'Invalid {tau_id} ONNX network config file')
130
131 acc.setPrivateTools(CompFactory.TauGNNEvaluator(
132 name = f'TrigTau_TauJetONNXEvaluator_{tau_id}',
133 useTRT = flags.Detector.EnableTRT,
134
135 # Network config:
136 **network_config,
137 InputLayerScalar = 'tau_vars',
138 InputLayerTracks = 'track_vars',
139 InputLayerClusters = 'cluster_vars',
140 NodeNameTau = id_flags.NodeNameTau if hasattr(id_flags, 'NodeNameTau') else 'pTau',
141 NodeNameJet = id_flags.NodeNameJet if hasattr(id_flags, 'NodeNameJet') else 'pJet',
142
143
144 # Inputs:
145 MaxTracks = id_flags.MaxTracks,
146 MaxClusters = id_flags.MaxClusters,
147 MaxClusterDR = 1.0,
148 VertexCorrection = False,
149 TrackClassification = False,
150
151 # Decorated TauJet variable names:
152 OutputVarname = f'{tau_id}_Score',
153 OutputDiscriminant = id_flags.OutputDiscriminant,
154 OutputPTau = f'{tau_id}_ProbTau',
155 OutputPJet = f'{tau_id}_ProbJet',
156
157 # Tau prongness selection
158 MinProngTrackPt = id_flags.MinProngTrackPt if hasattr(id_flags, 'MinProngTrackPt') else 0,
159 ))
160
161 return acc
162
163def trigTauWPDecoratorRNNCfg(flags, tau_id: str, precision_seq_name: str):
164 '''
165 TauJet signal transformed score and ID WPs decorator tool,
166 for the legacy mediumRNN/tightRNN_tracktwoMVA/tracktwoLLP/trackLRT chains ONLY!
167 '''
168
169 acc = ComponentAccumulator()
170
171 try: id_flags = getattr(flags.Trigger.Offline.Tau, tau_id)
172 except NameError: raise ValueError(f'Invalid TauID configuration: {tau_id}')
173
174 # In this version we store the WPs as flags accessable through
175 # tau->isTau(xAOD::TauJetParameters::IsTauFlags::JetRNNSig...)
176 import PyUtils.RootUtils as ru
177 ROOT = ru.import_root()
178 import cppyy
179 cppyy.load_library('libxAODTau_cDict')
180 WPEnumVals = [
181 ROOT.xAOD.TauJetParameters.IsTauFlag.JetRNNSigVeryLoose,
182 ROOT.xAOD.TauJetParameters.IsTauFlag.JetRNNSigLoose,
183 ROOT.xAOD.TauJetParameters.IsTauFlag.JetRNNSigMedium,
184 ROOT.xAOD.TauJetParameters.IsTauFlag.JetRNNSigTight,
185 ]
186
187 acc.setPrivateTools(CompFactory.TauWPDecorator(
188 name=f'TrigTau_TauWPDecoratorRNN_{precision_seq_name}_{tau_id}',
189 flatteningFile0Prong=id_flags.ScoreFlatteningConfig[0],
190 flatteningFile1Prong=id_flags.ScoreFlatteningConfig[1],
191 flatteningFile3Prong=id_flags.ScoreFlatteningConfig[2],
192 CutEnumVals=WPEnumVals,
193 SigEff0P=id_flags.TargetEff[0],
194 SigEff1P=id_flags.TargetEff[1],
195 SigEff3P=id_flags.TargetEff[2],
196 ScoreName='RNNJetScore',
197 NewScoreName='RNNJetScoreSigTrans',
198 DefineWPs=True,
199 ))
200
201 return acc
202
203def trigTauWPDecoratorCfg(flags, tau_id: str, precision_seq_name: str, tauContainerName: str):
204 '''TauJet signal transformed score and ID WPs decorator tool'''
205 acc = ComponentAccumulator()
206
207 try: id_flags = getattr(flags.Trigger.Offline.Tau, tau_id)
208 except NameError: raise ValueError(f'Invalid TauID configuration: {tau_id}')
209
210 if isinstance(id_flags.ScoreFlatteningConfig, str):
211 cfg = {
212 'flatteningFile0Prong': id_flags.ScoreFlatteningConfig,
213 'flatteningFile1Prong': id_flags.ScoreFlatteningConfig,
214 'flatteningFile3Prong': id_flags.ScoreFlatteningConfig,
215 }
216 elif isinstance(id_flags.ScoreFlatteningConfig, (list, tuple)) and len(id_flags.ScoreFlatteningConfig) == 3:
217 cfg = {
218 'flatteningFile0Prong': id_flags.ScoreFlatteningConfig[0],
219 'flatteningFile1Prong': id_flags.ScoreFlatteningConfig[1],
220 'flatteningFile3Prong': id_flags.ScoreFlatteningConfig[2],
221 }
222 else:
223 raise ValueError(f'Invalid {tau_id} WP decorator flattening config')
224
225
226 acc.setPrivateTools(CompFactory.TauWPDecorator(
227 name=f'TrigTau_TauWPDecorator_{precision_seq_name}_{tau_id}',
228 **cfg,
229 TauContainerName=tauContainerName,
230 DecorWPNames=[f'{tau_id}_{wp}' for wp in id_flags.TargetWPs],
231 DecorWPCutEffs0P=[eff[0] for eff in id_flags.TargetWPs.values()],
232 DecorWPCutEffs1P=[eff[1] for eff in id_flags.TargetWPs.values()],
233 DecorWPCutEffs3P=[eff[2] for eff in id_flags.TargetWPs.values()],
234 ScoreName=f'{tau_id}_Score',
235 NewScoreName=f'{tau_id}_ScoreSigTrans',
236 DefineWPs=True,
237 ))
238
239 return acc
TrigTauRecToolsConfig.trigTauVertexFinderCfg
trigTauVertexFinderCfg(flags, name='')
Definition TrigTauRecToolsConfig.py:7
Generated on for ATLAS Offline Software by doxygen 1.14.0