148def EGammaCookieCutClusterToolCfg(flags, name = 'EGCookieCutTool', **kwargs):
149 acc = ComponentAccumulator()
150
151 from LArGeoAlgsNV.LArGMConfig import LArGMCfg
152 acc.merge(LArGMCfg(flags))
153 from TileGeoModel.TileGMConfig import TileGMCfg
154 acc.merge(TileGMCfg(flags))
155
156 clusterKey = kwargs.setdefault("ClusterContainerName", "ForwardElectronCookieCutClusters")
157 kwargs.setdefault("ClusterContainerLinksName", clusterKey + "_links")
158
159 if kwargs["ClusterContainerLinksName"] != (clusterKey + "_links"):
160 raise AttributeError("ClusterContainerLinksName ({}) is not syncrhonised with ClusterContainerName ({})".format(kwargs["ClusterContainerLinksName"], clusterKey))
161 kwargs.setdefault('StoreCookedMoments', False)
162 kwargs.setdefault('StoreInputMoments', False)
163 kwargs.setdefault("SGKey_electrons", flags.Egamma.Keys.Output.ForwardElectrons)
164 momentNames = [
165 "CENTER_X", "CENTER_Y", "CENTER_Z",
166 "SECOND_LAMBDA", "LATERAL", "LONGITUDINAL", "ENG_FRAC_MAX",
167 "SECOND_R", "CENTER_LAMBDA", "SECOND_ENG_DENS", "SIGNIFICANCE" ]
168 import ROOT
169 moments = [
170 ROOT.xAOD.CaloCluster.CENTER_X, ROOT.xAOD.CaloCluster.CENTER_Y, ROOT.xAOD.CaloCluster.CENTER_Z,
171 ROOT.xAOD.CaloCluster.SECOND_LAMBDA, ROOT.xAOD.CaloCluster.LATERAL, ROOT.xAOD.CaloCluster.LONGITUDINAL, ROOT.xAOD.CaloCluster.ENG_FRAC_MAX,
172 ROOT.xAOD.CaloCluster.SECOND_R, ROOT.xAOD.CaloCluster.CENTER_LAMBDA, ROOT.xAOD.CaloCluster.SECOND_ENG_DENS, ROOT.xAOD.CaloCluster.SIGNIFICANCE ]
173 kwargs.setdefault("MomentNames", momentNames)
174 kwargs.setdefault("Moments", moments)
175 cookedMoments = [ "cookiecut" + moment for moment in momentNames] if kwargs['StoreCookedMoments'] else []
176 originalMoments = [ "original" + moment for moment in momentNames] if kwargs['StoreInputMoments'] else []
177 electronDecorations = [i for sublist in zip(cookedMoments, originalMoments) for i in sublist]
178 electronDecorations += ["cookiecutClusterLink"]
179 kwargs.setdefault("SGKey_electrons_decorations", electronDecorations)
180
181 from CaloTools.CaloNoiseCondAlgConfig import CaloNoiseCondAlgCfg
182 acc.merge(CaloNoiseCondAlgCfg(flags,"totalNoise"))
183 from CaloRec.CaloTopoClusterConfig import getTopoMoments
184 momentsMaker = acc.popToolsAndMerge(getTopoMoments(flags))
185 kwargs.setdefault("ClusterMomentMaker",[momentsMaker])
186
187 acc.setPrivateTools(
188 CompFactory.DerivationFramework.EGammaCookieCutClusterTool(
189 name, **kwargs))
190 return acc
191