ATLAS Offline Software
TrigEgammaPrecisionPhotonCaloIsoHypoAlg.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
3 */
4 
9 #include "AthViews/ViewHelper.h"
10 
11 namespace TCU = TrigCompositeUtils;
12 
14  ISvcLocator* pSvcLocator ) :
15  ::HypoBase( name, pSvcLocator ) {}
16 
17 
19 {
20  ATH_MSG_DEBUG ( "Initializing " << name() << "..." );
21  ATH_CHECK( m_hypoTools.retrieve() );
22  ATH_CHECK( m_photonsKey.initialize() );
23  ATH_CHECK( m_IsophotonsKey.initialize() );
24  renounce( m_photonsKey );// photons are made in views, so they are not in the EvtStore: hide them
25  renounce( m_IsophotonsKey );// photons are made in views, so they are not in the EvtStore: hide them
26 
27  return StatusCode::SUCCESS;
28 }
29 
30 
31 StatusCode TrigEgammaPrecisionPhotonCaloIsoHypoAlg::execute( const EventContext& context ) const
32 {
33  ATH_MSG_DEBUG ( "Executing " << name() << "..." );
34  auto previousDecisionsHandle = SG::makeHandle( decisionInput(), context );
35  ATH_CHECK( previousDecisionsHandle.isValid() );
36  ATH_MSG_DEBUG( "Running with "<< previousDecisionsHandle->size() <<" previous decisions");
37 
38 
39  // create handles for isolated photons
41  //make the output photon container
42  ATH_CHECK(h_isoPhotons.record(std::make_unique<xAOD::PhotonContainer>(),
43  std::make_unique<xAOD::PhotonAuxContainer>()));
44 
45  // new decisions
46 
47  // new output decisions
49  auto decisions = outputHandle.ptr();
50 
51  // input for decision
52  std::vector<ITrigEgammaPrecisionPhotonCaloIsoHypoTool::PhotonInfo> toolInput;
53 
54  // loop over previous decisions
55  size_t counter=0;
56  for ( auto previousDecision: *previousDecisionsHandle ) {
57 
58  //get updated RoI
59  auto roiELInfo = TCU::findLink<TrigRoiDescriptorCollection>( previousDecision, TCU::roiString() );
60 
61  ATH_CHECK( roiELInfo.isValid() );
62  const TrigRoiDescriptor* roi = *(roiELInfo.link);
63 
64  // not using View so commenting out following lines
65  const auto viewEL = previousDecision->objectLink<ViewContainer>( TCU::viewString() );
66  ATH_CHECK( viewEL.isValid() );
67  auto photonHandle = ViewHelper::makeHandle( *viewEL, m_photonsKey, context);
68  ATH_CHECK( photonHandle.isValid() );
69  ATH_MSG_DEBUG ( "Photon handle size: " << photonHandle->size() << "..." );
70  // Loop over the photonHandles
71  size_t validphotons=0;
72  for (size_t cl=0; cl< photonHandle->size(); cl++){
73 
74  {
75  auto ph = ViewHelper::makeLink( *viewEL, photonHandle, cl );
76  ATH_MSG_DEBUG ( "Checking ph.isValid()...");
77  if( !ph.isValid() ) {
78  ATH_MSG_DEBUG ( "PhotonHandle in position " << cl << " -> invalid ElemntLink!. Skipping...");
79  }
80 
81  ATH_CHECK(ph.isValid());
82 
83  ATH_MSG_DEBUG ( "PhotonHandle in position " << cl << " processing...");
84  ATH_MSG_DEBUG ( "Copying xAOD photon to new isolated container " << m_IsophotonsKey << " to run isolation tool decoration later...");
85  // Copy the photon in non-iso container to the To-Be-Isolation-Evaluated IsoPhoton container
86  xAOD::Photon *isoPhoton = new xAOD::Photon();
87  h_isoPhotons->push_back(isoPhoton);
88  const xAOD::Photon *originalPhoton = *ph;
89  *isoPhoton = *originalPhoton; // Deep copy
90 
92  d->setObjectLink( TCU::featureString(), ElementLink<xAOD::PhotonContainer>(*h_isoPhotons, isoPhoton->index(), context) );
94 
95  ITrigEgammaPrecisionPhotonCaloIsoHypoTool::PhotonInfo info(d, roi, isoPhoton, previousDecision);
96 
97 
98  toolInput.push_back(info);
99  validphotons++;
100  }
101  }
102 
103  ATH_MSG_DEBUG( "Photons with valid links: " << validphotons );
104 
105  ATH_MSG_DEBUG( "roi, photon, previous decision to new decision " << counter << " for roi " );
106  counter++;
107  }
108 
109  ATH_MSG_DEBUG( "Found "<<toolInput.size()<<" inputs to tools");
110 
111  for ( auto& tool: m_hypoTools ) {
112  ATH_CHECK( tool->decide( toolInput ) );
113  }
114 
115  ATH_CHECK( hypoBaseOutputProcessing(outputHandle) );
116  return StatusCode::SUCCESS;
117 }
TrigCompositeUtils::featureString
const std::string & featureString()
Definition: TrigCompositeUtils.h:420
TrigEgammaPrecisionPhotonCaloIsoHypoAlg::m_photonsKey
SG::ReadHandleKey< xAOD::PhotonContainer > m_photonsKey
Definition: TrigEgammaPrecisionPhotonCaloIsoHypoAlg.h:32
TrigEgammaPrecisionPhotonCaloIsoHypoAlg::TrigEgammaPrecisionPhotonCaloIsoHypoAlg
TrigEgammaPrecisionPhotonCaloIsoHypoAlg(const std::string &name, ISvcLocator *pSvcLocator)
Definition: TrigEgammaPrecisionPhotonCaloIsoHypoAlg.cxx: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
hist_file_dump.d
d
Definition: hist_file_dump.py:142
TrigEgammaPrecisionPhotonCaloIsoHypoAlg::m_IsophotonsKey
SG::WriteHandleKey< xAOD::PhotonContainer > m_IsophotonsKey
Definition: TrigEgammaPrecisionPhotonCaloIsoHypoAlg.h:33
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
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:18
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
TrigCompositeUtils::roiString
const std::string & roiString()
Definition: TrigCompositeUtils.h:418
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:22
SG::makeHandle
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
Definition: ReadCondHandle.h:274
TrigCompositeUtils.h
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:35
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
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
TrigEgammaPrecisionPhotonCaloIsoHypoAlg::initialize
virtual StatusCode initialize() override
Definition: TrigEgammaPrecisionPhotonCaloIsoHypoAlg.cxx:18
SG::AuxElement::index
size_t index() const
Return the index of this element within its container.
xAOD::decisions
decisions
Definition: TrigComposite_v1.cxx:101
TrigEgammaPrecisionPhotonCaloIsoHypoAlg::m_hypoTools
ToolHandleArray< ITrigEgammaPrecisionPhotonCaloIsoHypoTool > m_hypoTools
Definition: TrigEgammaPrecisionPhotonCaloIsoHypoAlg.h:31
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
DataVector::push_back
value_type push_back(value_type pElem)
Add an element to the end of the collection.
TrigCompositeUtils::viewString
const std::string & viewString()
Definition: TrigCompositeUtils.h:419
TrigEgammaPrecisionPhotonCaloIsoHypoAlg::execute
virtual StatusCode execute(const EventContext &context) const override
Definition: TrigEgammaPrecisionPhotonCaloIsoHypoAlg.cxx:31
TrigEgammaPrecisionPhotonCaloIsoHypoAlg.h
xAOD::Photon
Photon_v1 Photon
Definition of the current "egamma version".
Definition: Event/xAOD/xAODEgamma/xAODEgamma/Photon.h:17
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
SG::WriteHandle::record
StatusCode record(std::unique_ptr< T > data)
Record a const object to the store.
xAOD::Photon_v1
Definition: Photon_v1.h:37
TrigCompositeUtils::hypoAlgNodeName
const std::string & hypoAlgNodeName()
Definition: TrigCompositeUtils.h:426
HLTIdentifier.h
TrigCompositeUtils
Definition: Event/xAOD/xAODTrigger/xAODTrigger/TrigComposite.h:19
ITrigEgammaPrecisionPhotonCaloIsoHypoTool::PhotonInfo
Definition: ITrigEgammaPrecisionPhotonCaloIsoHypoTool.h:25
TrigRoiDescriptor.h
test_pyathena.counter
counter
Definition: test_pyathena.py:15
dq_make_web_display.cl
cl
print [x.__class__ for x in toList(dqregion.getSubRegions()) ]
Definition: dq_make_web_display.py:25
python.ParticleTypeUtil.info
def info
Definition: ParticleTypeUtil.py:87
mapkey::key
key
Definition: TElectronEfficiencyCorrectionTool.cxx:37
ViewContainer
Definition: View.h:158