ATLAS Offline Software
TrigEgammaFastCaloHypoAlg.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 #include <algorithm>
10 #include "AthViews/ViewHelper.h"
13 
14 
15 namespace TCU = TrigCompositeUtils;
16 
18  ISvcLocator* pSvcLocator ) :
19  ::HypoBase( name, pSvcLocator ),
20  m_lumiBlockMuTool("LumiBlockMuTool/LumiBlockMuTool")
21  {
22  declareProperty("LumiBlockMuTool", m_lumiBlockMuTool, "Luminosity Tool" );
23  }
24 
25 
27 
28  ATH_CHECK( m_hypoTools.retrieve() );
29  ATH_CHECK( m_ringerNNTools.retrieve());
30  ATH_CHECK( m_clustersKey.initialize() );
31  ATH_CHECK( m_ringsKey.initialize(SG::AllowEmpty));
32 
33  renounce( m_clustersKey );// clusters are made in views, so they are not in the EvtStore: hide them
35 
36  if ( not m_monTool.name().empty() )
37  CHECK( m_monTool.retrieve() );
38 
39 
40  ATH_MSG_DEBUG( "Initialization of FastCaloHypoAlg completed successfully");
41  return StatusCode::SUCCESS;
42 }
43 
44 
45 
46 StatusCode TrigEgammaFastCaloHypoAlg::execute( const EventContext& context ) const
47 {
48  ATH_MSG_DEBUG ( "Executing " << name() << "..." );
49 
50  auto timer = Monitored::Timer("TIME_exec");
51  auto timer_predict = Monitored::Timer("TIME_NN_exec");
52  auto monitoring = Monitored::Group( m_monTool, timer, timer_predict);
53 
54  auto previousDecisionsHandle = SG::makeHandle( decisionInput(), context );
55  ATH_CHECK( previousDecisionsHandle.isValid() );
56  ATH_MSG_DEBUG( "Running with "<< previousDecisionsHandle->size() <<" previous decisions");
57 
58  // new decisions
59  // new output decisions
61  auto decisions = outputHandle.ptr();
62 
63  // input for decision
64  std::vector<ITrigEgammaFastCaloHypoTool::FastClusterInfo> toolInput;
65 
66  // loop over previous decisions
67  size_t counter=0;
68  for ( const auto previousDecision: *previousDecisionsHandle )
69  {
70  //get RoI
71  auto roiELInfo = TCU::findLink<TrigRoiDescriptorCollection>( previousDecision, TCU::initialRoIString() );
72 
73  ATH_CHECK( roiELInfo.isValid() );
74  const TrigRoiDescriptor* roi = *(roiELInfo.link);
75 
76  // get View
77  const auto viewEL = previousDecision->objectLink<ViewContainer>( TCU::viewString() );
78  ATH_CHECK( viewEL.isValid() );
79 
80  // get cluster
81  auto clusterHandle = ViewHelper::makeHandle( *viewEL, m_clustersKey, context);
82  ATH_CHECK( clusterHandle.isValid() );
83  ATH_MSG_DEBUG ( "Cluster handle size: " << clusterHandle->size() << "..." );
84  auto emCluster = clusterHandle.cptr()->at(0);
85  ATH_MSG_DEBUG("Event cluster with eta = "<< emCluster->eta());
86 
87  // get rings
88  auto ringsHandle = ViewHelper::makeHandle( *viewEL, m_ringsKey, context);
89  ATH_CHECK( ringsHandle.isValid());
90  ATH_MSG_DEBUG ( "Rings handle size: " << ringsHandle->size() << "..." );
91 
92  auto ringsCluster = ringsHandle.cptr()->at(0);
93 
94  // create new decision (node)
96 
97 
98  ITrigEgammaFastCaloHypoTool::FastClusterInfo info(d, roi, emCluster, ringsCluster , previousDecision);
99 
100  float avgmu = m_lumiBlockMuTool->averageInteractionsPerCrossing(context);
101  info.valueDecorator["avgmu"] = avgmu;
102 
103  // Decorate the info object with NN ringer decision. if rings size is zero, this object is a dummy and
104  // should not compute any decision. the hypo tool will not found any pid name and will reprove this.
105  if(ringsCluster && emCluster && !ringsCluster->rings().empty() )
106  {
107  int idx=0;
108  for (auto& pidname : m_pidNames ){
109  timer_predict.start();
110  float nnOutput = m_ringerNNTools[idx]->predict(ringsCluster);
111  timer_predict.stop();
112  info.pidDecorator[pidname] = (bool)m_ringerNNTools[idx]->accept(ringsCluster, nnOutput, avgmu);
113  info.valueDecorator[pidname+"NNOutput"] = nnOutput;
114  idx++;
115  }
116  }
117 
118  toolInput.push_back( info );
119 
120 
121  // link the cluster
122  {
123  auto clus = ViewHelper::makeLink( *viewEL, clusterHandle, 0 );
124  ATH_CHECK( clus.isValid() );
125  d->setObjectLink( TCU::featureString(), clus );
126  }
127 
128  // link the rings
129  {
130  auto rings = ViewHelper::makeLink( *viewEL, ringsHandle, 0 );
131  ATH_CHECK( rings.isValid() );
132  d->setObjectLink( "ringer", rings );
133  }
134 
135  // link the RoI
136  {
137  d->setObjectLink( TCU::roiString(), roiELInfo.link );
138  }
139 
140  TCU::linkToPrevious( d, previousDecision, context );
141 
142  ATH_MSG_DEBUG( "Added view, roi, cluster, previous decision to new decision " << counter << " for view " << (*viewEL)->name() );
143  counter++;
144  }
145 
146  ATH_MSG_DEBUG( "Found "<<toolInput.size()<<" inputs to tools");
147 
148 
149  for ( auto& tool: m_hypoTools ) {
150  ATH_CHECK( tool->decide( toolInput ) );
151  }
152 
153  ATH_CHECK( hypoBaseOutputProcessing(outputHandle) );
154 
155  return StatusCode::SUCCESS;
156 }
ITrigEgammaFastCaloHypoTool::FastClusterInfo
Definition: ITrigEgammaFastCaloHypoTool.h:28
TrigEgammaFastCaloHypoAlg::execute
virtual StatusCode execute(const EventContext &context) const override
Definition: TrigEgammaFastCaloHypoAlg.cxx:46
TrigDefs::Group
Group
Properties of a chain group.
Definition: GroupProperties.h:13
TrigCompositeUtils::newDecisionIn
Decision * newDecisionIn(DecisionContainer *dc, const std::string &name)
Helper method to create a Decision object, place it in the container and return a pointer to it.
Definition: TrigCompositeUtilsRoot.cxx:43
TrigEgammaFastCaloHypoAlg::TrigEgammaFastCaloHypoAlg
TrigEgammaFastCaloHypoAlg(const std::string &name, ISvcLocator *pSvcLocator)
Definition: TrigEgammaFastCaloHypoAlg.cxx:17
TrigEgammaFastCaloHypoAlg::m_monTool
ToolHandle< GenericMonitoringTool > m_monTool
Definition: TrigEgammaFastCaloHypoAlg.h:36
hist_file_dump.d
d
Definition: hist_file_dump.py:142
ViewHelper::makeLink
ElementLink< T > makeLink(const SG::View *view, const SG::ReadHandle< T > &handle, size_t index)
Create EL to a collection in view.
Definition: ViewHelper.h:297
TrigCompositeUtils::hypoAlgNodeName
const std::string & hypoAlgNodeName()
Definition: TrigCompositeUtilsRoot.cxx:919
CutsMETMaker::accept
StatusCode accept(const xAOD::Muon *mu)
Definition: CutsMETMaker.cxx:18
ViewHelper::makeHandle
SG::ReadHandle< T > makeHandle(const SG::View *view, const SG::ReadHandleKey< T > &rhKey, const EventContext &context)
navigate from the TrigComposite to nearest view and fetch object from it
Definition: ViewHelper.h:265
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounce
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce(T &h)
Definition: AthCommonDataStore.h:380
HypoBase::decisionInput
const SG::ReadHandleKey< TrigCompositeUtils::DecisionContainer > & decisionInput() const
methods for derived classes to access handles of the base class input other read/write handles may be...
Definition: HypoBase.cxx:16
TrigRoiDescriptor
nope - should be used for standalone also, perhaps need to protect the class def bits #ifndef XAOD_AN...
Definition: TrigRoiDescriptor.h:56
ViewHelper.h
TrigCompositeUtils::createAndStore
SG::WriteHandle< DecisionContainer > createAndStore(const SG::WriteHandleKey< DecisionContainer > &key, const EventContext &ctx)
Creates and right away records the DecisionContainer with the key.
Definition: TrigCompositeUtilsRoot.cxx:27
TrigEgammaFastCaloHypoAlg::m_lumiBlockMuTool
ToolHandle< ILumiBlockMuTool > m_lumiBlockMuTool
Definition: TrigEgammaFastCaloHypoAlg.h:35
HypoBase::decisionOutput
const SG::WriteHandleKey< TrigCompositeUtils::DecisionContainer > & decisionOutput() const
methods for derived classes to access handles of the base class output other read/write handles may b...
Definition: HypoBase.cxx:20
python.utils.AtlRunQueryTimer.timer
def timer(name, disabled=False)
Definition: AtlRunQueryTimer.py:85
SG::makeHandle
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
Definition: ReadCondHandle.h:274
TrigCompositeUtils.h
TrigEgammaFastCaloHypoAlg::m_ringsKey
SG::ReadHandleKey< xAOD::TrigRingerRingsContainer > m_ringsKey
Definition: TrigEgammaFastCaloHypoAlg.h:32
TrigCompositeUtils::initialRoIString
const std::string & initialRoIString()
Definition: TrigCompositeUtilsRoot.cxx:883
TrigEgammaFastCaloHypoAlg::initialize
virtual StatusCode initialize() override
Definition: TrigEgammaFastCaloHypoAlg.cxx:26
HypoBase::hypoBaseOutputProcessing
StatusCode hypoBaseOutputProcessing(SG::WriteHandle< TrigCompositeUtils::DecisionContainer > &outputHandle, MSG::Level lvl=MSG::DEBUG) const
Base class function to be called once slice specific code has finished. Handles debug printing and va...
Definition: HypoBase.cxx:33
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
SG::WriteHandle::ptr
pointer_type ptr()
Dereference the pointer.
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition: AthCommonDataStore.h:145
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
CHECK
#define CHECK(...)
Evaluate an expression and check for errors.
Definition: Control/AthenaKernel/AthenaKernel/errorcheck.h:422
Monitored.h
Header file to be included by clients of the Monitored infrastructure.
xAOD::decisions
decisions
Definition: TrigComposite_v1.cxx:101
TrigCompositeContainer.h
TrigEgammaFastCaloHypoAlg::m_clustersKey
SG::ReadHandleKey< xAOD::TrigEMClusterContainer > m_clustersKey
Definition: TrigEgammaFastCaloHypoAlg.h:31
TrigEgammaFastCaloHypoAlg.h
TrigCompositeUtils::featureString
const std::string & featureString()
Definition: TrigCompositeUtilsRoot.cxx:899
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
HypoBase
Hypothesis algorithms take the output of reco algorithms and the decision from the preceeding InputMa...
Definition: HypoBase.h:13
AtlCoolConsole.tool
tool
Definition: AtlCoolConsole.py:452
SG::WriteHandle
Definition: StoreGate/StoreGate/WriteHandle.h:73
TrigCompositeUtils::linkToPrevious
void linkToPrevious(Decision *d, const std::string &previousCollectionKey, size_t previousIndex)
Links to the previous object, location of previous 'seed' decision supplied by hand.
Definition: TrigCompositeUtilsRoot.cxx:139
HLTIdentifier.h
TrigEgammaFastCaloHypoAlg::m_pidNames
Gaudi::Property< std::vector< std::string > > m_pidNames
Definition: TrigEgammaFastCaloHypoAlg.h:37
LArNewCalib_DelayDump_OFC_Cali.idx
idx
Definition: LArNewCalib_DelayDump_OFC_Cali.py:69
TrigCompositeUtils
Definition: Event/xAOD/xAODTrigger/xAODTrigger/TrigComposite.h:19
TrigCompositeUtils::roiString
const std::string & roiString()
Definition: TrigCompositeUtilsRoot.cxx:891
test_pyathena.counter
counter
Definition: test_pyathena.py:15
xAOD::bool
setBGCode setTAP setLVL2ErrorBits bool
Definition: TrigDecision_v1.cxx:60
SG::AllowEmpty
@ AllowEmpty
Definition: StoreGate/StoreGate/VarHandleKey.h:30
TrigCompositeUtils::viewString
const std::string & viewString()
Definition: TrigCompositeUtilsRoot.cxx:895
python.ParticleTypeUtil.info
def info
Definition: ParticleTypeUtil.py:87
TrigRoiDescriptorCollection.h
Monitored::Timer
A monitored timer.
Definition: MonitoredTimer.h:32
TrigEgammaFastCaloHypoAlg::m_hypoTools
ToolHandleArray< ITrigEgammaFastCaloHypoTool > m_hypoTools
Definition: TrigEgammaFastCaloHypoAlg.h:33
TrigEgammaFastCaloHypoAlg::m_ringerNNTools
ToolHandleArray< Ringer::IAsgRingerSelectorTool > m_ringerNNTools
Definition: TrigEgammaFastCaloHypoAlg.h:34
ViewContainer
Definition: View.h:158