ATLAS Offline Software
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("doVeryLoose", False)
10  kwargs.setdefault("doLoose", False)
11  kwargs.setdefault("doMedium", False)
12  kwargs.setdefault("doTight", False)
13 
14  acc = ComponentAccumulator()
15 
16  # tau selection relies on RNN electron veto, we must decorate the fixed eveto WPs before applying tau selection
17  acc.merge(AddTauIDDecorationCfg(flags, TauContainerName="TauJets"))
18 
19  from DerivationFrameworkTools.DerivationFrameworkToolsConfig import AsgSelectionToolWrapperCfg
20  from TauAnalysisTools.TauAnalysisToolsConfig import TauSelectionToolCfg
21 
22  TauAugmentationTools = []
23 
24  if kwargs["doVeryLoose"]:
25  TauSelectorVeryLoose = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
26  name = 'TauSelectorVeryLoose',
27  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_veryloose.conf'))
28  acc.addPublicTool(TauSelectorVeryLoose)
29 
30  TauVeryLooseWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
31  name = "TauVeryLooseWrapper",
32  AsgSelectionTool = TauSelectorVeryLoose,
33  StoreGateEntryName = "DFTauVeryLoose",
34  ContainerName = "TauJets"))
35  TauAugmentationTools.append(TauVeryLooseWrapper)
36 
37  if kwargs["doLoose"]:
38  TauSelectorLoose = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
39  name = 'TauSelectorLoose',
40  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_loose.conf'))
41  acc.addPublicTool(TauSelectorLoose)
42 
43  TauLooseWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
44  name = "TauLooseWrapper",
45  AsgSelectionTool = TauSelectorLoose,
46  StoreGateEntryName = "DFTauLoose",
47  ContainerName = "TauJets"))
48  TauAugmentationTools.append(TauLooseWrapper)
49 
50  if kwargs["doMedium"]:
51  TauSelectorMedium = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
52  name = 'TauSelectorMedium',
53  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_medium.conf'))
54  acc.addPublicTool(TauSelectorMedium)
55 
56  TauMediumWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
57  name = "TauMediumWrapper",
58  AsgSelectionTool = TauSelectorMedium,
59  StoreGateEntryName = "DFTauMedium",
60  ContainerName = "TauJets"))
61  TauAugmentationTools.append(TauMediumWrapper)
62 
63  if kwargs["doTight"]:
64  TauSelectorTight = acc.popToolsAndMerge(TauSelectionToolCfg(flags,
65  name = 'TauSelectorTight',
66  ConfigPath = 'TauAnalysisAlgorithms/tau_selection_tight.conf'))
67  acc.addPublicTool(TauSelectorTight)
68 
69  TauTightWrapper = acc.getPrimaryAndMerge(AsgSelectionToolWrapperCfg(flags,
70  name = "TauTightWrapper",
71  AsgSelectionTool = TauSelectorTight,
72  StoreGateEntryName = "DFTauTight",
73  ContainerName = "TauJets"))
74  TauAugmentationTools.append(TauTightWrapper)
75 
76  if TauAugmentationTools:
77  CommonAugmentation = CompFactory.DerivationFramework.CommonAugmentation
78  acc.addEventAlgo(CommonAugmentation(f"{prefix}_TauAugmentationKernel", AugmentationTools = TauAugmentationTools))
79 
80  return acc
81 
82 
83 # Low pT di-taus
84 def AddDiTauLowPtCfg(flags, **kwargs):
85  """Configure the low-pt di-tau building"""
86 
87  acc = ComponentAccumulator()
88 
89  from JetRecConfig.JetRecConfig import JetRecCfg
90  from JetRecConfig.StandardLargeRJets import AntiKt10LCTopo
91  acc.merge(JetRecCfg(flags,AntiKt10LCTopo))
92 
93  from DiTauRec.DiTauBuilderConfig import DiTauBuilderLowPtCfg
94  acc.merge(DiTauBuilderLowPtCfg(flags, name="DiTauLowPtBuilder"))
95 
96  return acc
97 
98 
99 def AddTauIDDecorationCfg(flags, **kwargs):
100  """Decorate tau ID scores and working points"""
101 
102  kwargs.setdefault("evetoFix", True)
103  kwargs.setdefault("GNNTauID", True)
104  kwargs.setdefault("TauContainerName", "TauJets")
105  kwargs.setdefault("prefix", kwargs['TauContainerName'])
106 
107  acc = ComponentAccumulator()
108 
109  import tauRec.TauToolHolder as tauTools
110  tools = []
111 
112  if kwargs['evetoFix']:
113  tools.append( acc.popToolsAndMerge(tauTools.TauWPDecoratorEleRNNFixCfg(flags)) )
114 
115  if kwargs['GNNTauID']:
116  # vertex-corrected clusters must be rebuilt for tau ID
117  tools.append( acc.popToolsAndMerge(tauTools.TauVertexedClusterDecoratorCfg(flags)) )
118  # Add in GNTau!
119  tools.append( acc.popToolsAndMerge(tauTools.TauGNNEvaluatorCfg(flags)) )
120  tools.append( acc.popToolsAndMerge(tauTools.TauWPDecoratorGNNCfg(flags)) )
121 
122  if tools:
123  for tool in tools:
124  acc.addPublicTool(tool)
125 
126  TauIDDecoratorWrapper = CompFactory.DerivationFramework.TauIDDecoratorWrapper
127  TauIDDecoratorKernel = CompFactory.DerivationFramework.CommonAugmentation
128 
129  prefix = kwargs['prefix']
130  tauIDDecoratorWrapper = TauIDDecoratorWrapper(name = f"{prefix}_TauIDDecoratorWrapper",
131  TauContainerName = kwargs['TauContainerName'],
132  TauIDTools = tools)
133 
134  acc.addPublicTool(tauIDDecoratorWrapper)
135  acc.addEventAlgo(TauIDDecoratorKernel(name = f"{prefix}_TauIDDecorKernel",
136  AugmentationTools = [tauIDDecoratorWrapper]))
137 
138  return acc
139 
140 
141 # TauJets_MuonRM steering
142 def AddMuonRemovalTauAODReRecoAlgCfg(flags, **kwargs):
143  """Configure the MuonRM AOD tau building"""
144 
145  acc = ComponentAccumulator()
146 
147  # get tools from holder
148  import tauRec.TauToolHolder as tauTools
149  tools_mod = []
150  tools_mod.append( acc.popToolsAndMerge(tauTools.TauAODMuonRemovalCfg(flags)) )
151  tools_after = []
152  tools_after.append( acc.popToolsAndMerge(tauTools.TauVertexedClusterDecoratorCfg(flags)) )
153  tools_after.append( acc.popToolsAndMerge(tauTools.TauTrackRNNClassifierCfg(flags)) )
154  tools_after.append( acc.popToolsAndMerge(tauTools.EnergyCalibrationLCCfg(flags)) )
155  tools_after.append( acc.popToolsAndMerge(tauTools.TauCommonCalcVarsCfg(flags)) )
156  tools_after.append( acc.popToolsAndMerge(tauTools.TauSubstructureCfg(flags)) )
157  tools_after.append( acc.popToolsAndMerge(tauTools.Pi0ClusterCreatorCfg(flags)) )
158  tools_after.append( acc.popToolsAndMerge(tauTools.Pi0ClusterScalerCfg(flags)) )
159  tools_after.append( acc.popToolsAndMerge(tauTools.Pi0ScoreCalculatorCfg(flags)) )
160  tools_after.append( acc.popToolsAndMerge(tauTools.Pi0SelectorCfg(flags)) )
161  tools_after.append( acc.popToolsAndMerge(tauTools.TauVertexVariablesCfg(flags)) )
162  import PanTauAlgs.JobOptions_Main_PanTau_New as pantau
163  tools_after.append( acc.popToolsAndMerge(pantau.PanTauCfg(flags)) )
164  tools_after.append( acc.popToolsAndMerge(tauTools.TauCombinedTESCfg(flags)) )
165  tools_after.append( acc.popToolsAndMerge(tauTools.MvaTESVariableDecoratorCfg(flags)) )
166  tools_after[-1].EventShapeKey = ''
167  tools_after.append( acc.popToolsAndMerge(tauTools.MvaTESEvaluatorCfg(flags)) )
168  tools_after.append( acc.popToolsAndMerge(tauTools.TauIDVarCalculatorCfg(flags)) )
169  tools_after.append( acc.popToolsAndMerge(tauTools.TauJetRNNEvaluatorCfg(flags)) )
170  tools_after.append( acc.popToolsAndMerge(tauTools.TauWPDecoratorJetRNNCfg(flags)) )
171  tools_after.append( acc.popToolsAndMerge(tauTools.TauEleRNNEvaluatorCfg(flags)) )
172  tools_after.append( acc.popToolsAndMerge(tauTools.TauWPDecoratorEleRNNCfg(flags)) )
173  tools_after.append( acc.popToolsAndMerge(tauTools.TauDecayModeNNClassifierCfg(flags)) )
174  TauAODRunnerAlg=CompFactory.getComp("TauAODRunnerAlg")
175  for tool in tools_mod:
176  tool.inAOD = True
177  for tool in tools_after:
178  tool.inAOD = True
179  myTauAODRunnerAlg = TauAODRunnerAlg(
180  name = "MuonRemovalTauAODReRecoAlg",
181  Key_tauOutputContainer = "TauJets_MuonRM",
182  Key_pi0OutputContainer = "TauFinalPi0s_MuonRM",
183  Key_neutralPFOOutputContainer = "TauNeutralParticleFlowObjects_MuonRM",
184  Key_chargedPFOOutputContainer = "TauChargedParticleFlowObjects_MuonRM",
185  Key_hadronicPFOOutputContainer = "TauHadronicParticleFlowObjects_MuonRM",
186  Key_tauTrackOutputContainer = "TauTracks_MuonRM",
187  Key_vertexOutputContainer = "TauSecondaryVertices_MuonRM",
188  modificationTools = tools_mod,
189  officialTools = tools_after
190  )
191  acc.addEventAlgo(myTauAODRunnerAlg)
192  return acc
193 
194 
195 def TauThinningCfg(flags, name, **kwargs):
196  """configure tau thinning"""
197 
198  acc = ComponentAccumulator()
199  TauThinningTool = CompFactory.DerivationFramework.TauThinningTool
200  acc.addPublicTool(TauThinningTool(name, **kwargs), primary=True)
201  return acc
python.JetAnalysisCommon.ComponentAccumulator
ComponentAccumulator
Definition: JetAnalysisCommon.py:302
TauAnalysisToolsConfig.TauSelectionToolCfg
def TauSelectionToolCfg(ConfigFlags, name, **kwargs)
Definition: TauAnalysisToolsConfig.py:6
python.TauCommonConfig.AddMuonRemovalTauAODReRecoAlgCfg
def AddMuonRemovalTauAODReRecoAlgCfg(flags, **kwargs)
Definition: TauCommonConfig.py:142
DiTauBuilderConfig.DiTauBuilderLowPtCfg
def DiTauBuilderLowPtCfg(flags, name="DiTauLowPtBuilder", **kwargs)
Definition: DiTauBuilderConfig.py:35
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:37
python.DerivationFrameworkToolsConfig.AsgSelectionToolWrapperCfg
def AsgSelectionToolWrapperCfg(ConfigFlags, name, **kwargs)
Definition: DerivationFrameworkToolsConfig.py:11
python.TauCommonConfig.AddTauIDDecorationCfg
def AddTauIDDecorationCfg(flags, **kwargs)
Definition: TauCommonConfig.py:99
TauAODRunnerAlg
The implementation of the TauAODRunnerAlg, which is meant to run at AOD level.
Definition: TauAODRunnerAlg.h:40
python.TauCommonConfig.TauThinningCfg
def TauThinningCfg(flags, name, **kwargs)
Definition: TauCommonConfig.py:195
python.TauCommonConfig.AddDiTauLowPtCfg
def AddDiTauLowPtCfg(flags, **kwargs)
Definition: TauCommonConfig.py:84