Loading [MathJax]/extensions/tex2jax.js
ATLAS Offline Software
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
TauCommonConfig.py
Go to the documentation of this file.
1 # Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
2 
3 from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
4 from AthenaConfiguration.ComponentFactory import CompFactory
5 
6 def AddTauAugmentationCfg(flags, **kwargs):
7 
8  prefix = kwargs["prefix"]
9  kwargs.setdefault("TauContainerName", "TauJets")
10  kwargs.setdefault("doRNNVeryLoose", False)
11  kwargs.setdefault("doRNNLoose", False)
12  kwargs.setdefault("doRNNMedium", False)
13  kwargs.setdefault("doRNNTight", False)
14  kwargs.setdefault("doGNTauVeryLoose", False)
15  kwargs.setdefault("doGNTauLoose", False)
16  kwargs.setdefault("doGNTauMedium", False)
17  kwargs.setdefault("doGNTauTight", False)
18 
19  acc = ComponentAccumulator()
20 
21  # tau selection relies on RNN electron veto, we must decorate the fixed eveto WPs before applying tau selection
22  acc.merge(AddTauIDDecorationCfg(flags, TauContainerName=kwargs["TauContainerName"]))
23 
24  from DerivationFrameworkTools.DerivationFrameworkToolsConfig import AsgSelectionToolWrapperCfg
25  from TauAnalysisTools.TauAnalysisToolsConfig import TauSelectionToolCfg
26 
27  TauAugmentationTools = []
28 
29  # RNN TauID WPs
30  if kwargs["doRNNVeryLoose"]:
31  TauSelectorRNNVeryLoose = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
32  name = 'TauSelectorRNNVeryLoose',
33  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_veryloose_noeleid.conf'))
34  acc.addPublicTool(TauSelectorRNNVeryLoose)
35 
36  TauRNNVeryLooseWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
37  name = "TauRNNVeryLooseWrapper",
38  AsgSelectionTool = TauSelectorRNNVeryLoose,
39  StoreGateEntryName = "DFTauRNNVeryLoose",
40  ContainerName = kwargs["TauContainerName"]))
41  TauAugmentationTools.append(TauRNNVeryLooseWrapper)
42 
43  if kwargs["doRNNLoose"]:
44  TauSelectorRNNLoose = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
45  name = 'TauSelectorRNNLoose',
46  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_loose_noeleid.conf'))
47  acc.addPublicTool(TauSelectorRNNLoose)
48 
49  TauRNNLooseWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
50  name = "TauRNNLooseWrapper",
51  AsgSelectionTool = TauSelectorRNNLoose,
52  StoreGateEntryName = "DFTauRNNLoose",
53  ContainerName = kwargs["TauContainerName"]))
54  TauAugmentationTools.append(TauRNNLooseWrapper)
55 
56  if kwargs["doRNNMedium"]:
57  TauSelectorRNNMedium = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
58  name = 'TauSelectorRNNMedium',
59  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_medium_noeleid.conf'))
60  acc.addPublicTool(TauSelectorRNNMedium)
61 
62  TauRNNMediumWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
63  name = "TauRNNMediumWrapper",
64  AsgSelectionTool = TauSelectorRNNMedium,
65  StoreGateEntryName = "DFTauRNNMedium",
66  ContainerName = kwargs["TauContainerName"]))
67  TauAugmentationTools.append(TauRNNMediumWrapper)
68 
69  if kwargs["doRNNTight"]:
70  TauSelectorRNNTight = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
71  name = 'TauSelectorRNNTight',
72  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_tight_noeleid.conf'))
73  acc.addPublicTool(TauSelectorRNNTight)
74 
75  TauRNNTightWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
76  name = "TauRNNTightWrapper",
77  AsgSelectionTool = TauSelectorRNNTight,
78  StoreGateEntryName = "DFTauRNNTight",
79  ContainerName = kwargs["TauContainerName"]))
80  TauAugmentationTools.append(TauRNNTightWrapper)
81 
82  # GNTau TauID WPs
83  if kwargs["doGNTauVeryLoose"]:
84  TauSelectorGNTauVeryLoose = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
85  name = 'TauSelectorGNTauVeryLoose',
86  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_gntau_veryloose_noeleid.conf'))
87  acc.addPublicTool(TauSelectorGNTauVeryLoose)
88 
89  TauGNTauVeryLooseWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
90  name = "TauGNTauVeryLooseWrapper",
91  AsgSelectionTool = TauSelectorGNTauVeryLoose,
92  StoreGateEntryName = "DFTauGNTauVeryLoose",
93  ContainerName = kwargs["TauContainerName"]))
94  TauAugmentationTools.append(TauGNTauVeryLooseWrapper)
95 
96 
97  if kwargs["doGNTauLoose"]:
98  TauSelectorGNTauLoose = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
99  name = 'TauSelectorGNTauLoose',
100  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_gntau_loose_noeleid.conf'))
101  acc.addPublicTool(TauSelectorGNTauLoose)
102 
103  TauGNTauLooseWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
104  name = "TauGNTauLooseWrapper",
105  AsgSelectionTool = TauSelectorGNTauLoose,
106  StoreGateEntryName = "DFTauGNTauLoose",
107  ContainerName = kwargs["TauContainerName"]))
108  TauAugmentationTools.append(TauGNTauLooseWrapper)
109 
110  if kwargs["doGNTauMedium"]:
111  TauSelectorGNTauMedium = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
112  name = 'TauSelectorGNTauMedium',
113  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_gntau_medium_noeleid.conf'))
114  acc.addPublicTool(TauSelectorGNTauMedium)
115 
116  TauGNTauMediumWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
117  name = "TauGNTauMediumWrapper",
118  AsgSelectionTool = TauSelectorGNTauMedium,
119  StoreGateEntryName = "DFTauGNTauMedium",
120  ContainerName = kwargs["TauContainerName"]))
121  TauAugmentationTools.append(TauGNTauMediumWrapper)
122 
123  if kwargs["doGNTauTight"]:
124  TauSelectorGNTauTight = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
125  name = 'TauSelectorGNTauTight',
126  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_gntau_tight_noeleid.conf'))
127  acc.addPublicTool(TauSelectorGNTauTight)
128 
129  TauGNTauTightWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
130  name = "TauGNTauTightWrapper",
131  AsgSelectionTool = TauSelectorGNTauTight,
132  StoreGateEntryName = "DFTauGNTauTight",
133  ContainerName = kwargs["TauContainerName"]))
134  TauAugmentationTools.append(TauGNTauTightWrapper)
135 
136 
137  if TauAugmentationTools:
138  CommonAugmentation = CompFactory.DerivationFramework.CommonAugmentation
139  acc.addEventAlgo(CommonAugmentation(f"{prefix}_TauAugmentationKernel", AugmentationTools = TauAugmentationTools))
140 
141  return acc
142 
143 
144 # Low pT di-taus
145 def AddDiTauLowPtCfg(flags, **kwargs):
146  """Configure the low-pt di-tau building"""
147 
148  acc = ComponentAccumulator()
149 
150  from JetRecConfig.JetRecConfig import JetRecCfg
151  from JetRecConfig.StandardLargeRJets import AntiKt10LCTopo
152  acc.merge(JetRecCfg(flags,AntiKt10LCTopo))
153 
154  from DiTauRec.DiTauBuilderConfig import DiTauBuilderCfg
155  acc.merge(DiTauBuilderCfg(flags, name="DiTauLowPtBuilder", doLowPt=True))
156 
157  return acc
158 
159 
160 def AddTauIDDecorationCfg(flags, **kwargs):
161  """Decorate tau ID scores and working points"""
162 
163  kwargs.setdefault("evetoFix", True)
164  kwargs.setdefault("GNNTauID", True)
165  kwargs.setdefault("TauContainerName", "TauJets")
166  kwargs.setdefault("prefix", kwargs['TauContainerName'])
167 
168  acc = ComponentAccumulator()
169 
170  import tauRec.TauToolHolder as tauTools
171  tools = []
172 
173  if kwargs['evetoFix']:
174  tools.append( acc.popToolsAndMerge(tauTools.TauWPDecoratorEleRNNFixCfg(flags)) )
175 
176  if kwargs['GNNTauID']:
177  # vertex-corrected clusters must be rebuilt for tau ID
178  tools.append( acc.popToolsAndMerge(tauTools.TauVertexedClusterDecoratorCfg(flags)) )
179  # Add in GNTau!
180  tools.append( acc.popToolsAndMerge(tauTools.TauGNNEvaluatorCfg(flags,0)) )
181  tools.append( acc.popToolsAndMerge(tauTools.TauGNNEvaluatorCfg(flags,1)) )
182  tools.append( acc.popToolsAndMerge(tauTools.TauWPDecoratorGNNCfg(flags,0)) )
183  tools.append( acc.popToolsAndMerge(tauTools.TauWPDecoratorGNNCfg(flags,1)) )
184 
185  if tools:
186  for tool in tools:
187  acc.addPublicTool(tool)
188 
189  TauIDDecoratorWrapper = CompFactory.DerivationFramework.TauIDDecoratorWrapper
190  TauIDDecoratorKernel = CompFactory.DerivationFramework.CommonAugmentation
191 
192  prefix = kwargs['prefix']
193  tauIDDecoratorWrapper = TauIDDecoratorWrapper(name = f"{prefix}_TauIDDecoratorWrapper",
194  TauContainerName = kwargs['TauContainerName'],
195  TauIDTools = tools)
196 
197  acc.addPublicTool(tauIDDecoratorWrapper)
198  acc.addEventAlgo(TauIDDecoratorKernel(name = f"{prefix}_TauIDDecorKernel",
199  AugmentationTools = [tauIDDecoratorWrapper]))
200 
201  return acc
202 
203 
204 # TauJets_MuonRM steering
205 def AddMuonRemovalTauAODReRecoAlgCfg(flags, **kwargs):
206  """Configure the MuonRM AOD tau building"""
207 
208  acc = ComponentAccumulator()
209 
210  # get tools from holder
211  import tauRec.TauToolHolder as tauTools
212  tools_mod = []
213  tools_mod.append( acc.popToolsAndMerge(tauTools.TauAODMuonRemovalCfg(flags)) )
214  tools_after = []
215  tools_after.append( acc.popToolsAndMerge(tauTools.TauVertexedClusterDecoratorCfg(flags)) )
216  tools_after.append( acc.popToolsAndMerge(tauTools.TauTrackRNNClassifierCfg(flags)) )
217  tools_after.append( acc.popToolsAndMerge(tauTools.EnergyCalibrationLCCfg(flags)) )
218  tools_after.append( acc.popToolsAndMerge(tauTools.TauCommonCalcVarsCfg(flags)) )
219  tools_after.append( acc.popToolsAndMerge(tauTools.TauSubstructureCfg(flags)) )
220  tools_after.append( acc.popToolsAndMerge(tauTools.Pi0ClusterCreatorCfg(flags)) )
221  tools_after.append( acc.popToolsAndMerge(tauTools.Pi0ClusterScalerCfg(flags)) )
222  tools_after.append( acc.popToolsAndMerge(tauTools.Pi0ScoreCalculatorCfg(flags)) )
223  tools_after.append( acc.popToolsAndMerge(tauTools.Pi0SelectorCfg(flags)) )
224  tools_after.append( acc.popToolsAndMerge(tauTools.TauVertexVariablesCfg(flags)) )
225  import PanTauAlgs.JobOptions_Main_PanTau as pantau
226  tools_after.append( acc.popToolsAndMerge(pantau.PanTauCfg(flags)) )
227  tools_after.append( acc.popToolsAndMerge(tauTools.TauCombinedTESCfg(flags)) )
228  tools_after.append( acc.popToolsAndMerge(tauTools.MvaTESVariableDecoratorCfg(flags)) )
229  tools_after[-1].EventShapeKey = ''
230  tools_after.append( acc.popToolsAndMerge(tauTools.MvaTESEvaluatorCfg(flags)) )
231  tools_after.append( acc.popToolsAndMerge(tauTools.TauIDVarCalculatorCfg(flags)) )
232  tools_after.append( acc.popToolsAndMerge(tauTools.TauJetRNNEvaluatorCfg(flags)) )
233  tools_after.append( acc.popToolsAndMerge(tauTools.TauWPDecoratorJetRNNCfg(flags)) )
234  tools_after.append( acc.popToolsAndMerge(tauTools.TauEleRNNEvaluatorCfg(flags)) )
235  tools_after.append( acc.popToolsAndMerge(tauTools.TauWPDecoratorEleRNNCfg(flags)) )
236  tools_after.append( acc.popToolsAndMerge(tauTools.TauDecayModeNNClassifierCfg(flags)) )
237  TauAODRunnerAlg=CompFactory.getComp("TauAODRunnerAlg")
238  for tool in tools_mod:
239  tool.inAOD = True
240  for tool in tools_after:
241  tool.inAOD = True
242  myTauAODRunnerAlg = TauAODRunnerAlg(
243  name = "MuonRemovalTauAODReRecoAlg",
244  Key_tauOutputContainer = "TauJets_MuonRM",
245  Key_pi0OutputContainer = "TauFinalPi0s_MuonRM",
246  Key_neutralPFOOutputContainer = "TauNeutralParticleFlowObjects_MuonRM",
247  Key_chargedPFOOutputContainer = "TauChargedParticleFlowObjects_MuonRM",
248  Key_hadronicPFOOutputContainer = "TauHadronicParticleFlowObjects_MuonRM",
249  Key_tauTrackOutputContainer = "TauTracks_MuonRM",
250  Key_vertexOutputContainer = "TauSecondaryVertices_MuonRM",
251  modificationTools = tools_mod,
252  officialTools = tools_after
253  )
254  acc.addEventAlgo(myTauAODRunnerAlg)
255  return acc
256 
257 
258 def TauThinningCfg(flags, name, **kwargs):
259  """configure tau thinning"""
260 
261  acc = ComponentAccumulator()
262  TauThinningTool = CompFactory.DerivationFramework.TauThinningTool
263  acc.addPublicTool(TauThinningTool(name, **kwargs), primary=True)
264  return acc
python.JetAnalysisCommon.ComponentAccumulator
ComponentAccumulator
Definition: JetAnalysisCommon.py:302
python.TauCommonConfig.AddMuonRemovalTauAODReRecoAlgCfg
def AddMuonRemovalTauAODReRecoAlgCfg(flags, **kwargs)
Definition: TauCommonConfig.py:205
TauAnalysisToolsConfig.TauSelectionToolCfg
def TauSelectionToolCfg(flags, name, **kwargs)
Definition: TauAnalysisToolsConfig.py:6
DiTauBuilderConfig.DiTauBuilderCfg
def DiTauBuilderCfg(flags, name="DiTauBuilder", doLowPt=False, **kwargs)
Definition: DiTauBuilderConfig.py:6
python.TauCommonConfig.AddTauAugmentationCfg
def AddTauAugmentationCfg(flags, **kwargs)
Definition: TauCommonConfig.py:6
python.JetRecConfig.JetRecCfg
def JetRecCfg(flags, jetdef, returnConfiguredDef=False)
Top level functions returning ComponentAccumulator out of JetDefinition.
Definition: JetRecConfig.py:36
python.DerivationFrameworkToolsConfig.AsgSelectionToolWrapperCfg
def AsgSelectionToolWrapperCfg(ConfigFlags, name, **kwargs)
Definition: DerivationFrameworkToolsConfig.py:11
python.TauCommonConfig.AddTauIDDecorationCfg
def AddTauIDDecorationCfg(flags, **kwargs)
Definition: TauCommonConfig.py:160
TauAODRunnerAlg
The implementation of the TauAODRunnerAlg, which is meant to run at AOD level.
Definition: TauAODRunnerAlg.h:44
python.TauCommonConfig.TauThinningCfg
def TauThinningCfg(flags, name, **kwargs)
Definition: TauCommonConfig.py:258
python.TauCommonConfig.AddDiTauLowPtCfg
def AddDiTauLowPtCfg(flags, **kwargs)
Definition: TauCommonConfig.py:145
Generated on Mon May 12 2025 21:18:33 for ATLAS Offline Software by   doxygen 1.8.18