xAODJetROIRetriever.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#include "TrigJiveXML/xAODJetROIRetriever.h"
 
#include <string>
 
#include "CLHEP/Units/SystemOfUnits.h"
 
#include "xAODTrigger/JetRoIContainer.h"
 
namespace JiveXML {
 
  //--------------------------------------------------------------------------
 
  xAODJetROIRetriever::xAODJetROIRetriever(const std::string& type, const std::string& name, const IInterface* parent):
    AthAlgTool(type, name, parent), m_typeName("JetROI")
  {
 
    declareInterface<IDataRetriever>(this);
  }
 
  //--------------------------------------------------------------------------
 
  StatusCode xAODJetROIRetriever::retrieve(ToolHandle<IFormatTool> &FormatTool) {
 
  DataVect phi;
  DataVect eta; 
  DataVect energy; 
  DataVect roiWord; 
  DataVect thrPattern; 
 
 
  const xAOD::JetRoIContainer* jetROIs = 0; 
 
 
    // L1JetObject -not- available
    m_sgKey = "LVL1JetRoIs"; 
    if ( evtStore()->retrieve(jetROIs,m_sgKey).isFailure() ) {
      if (msgLvl(MSG::DEBUG)) msg(MSG::DEBUG) <<  "No LVL1JetROIs found in SG " << endmsg;
      return StatusCode::SUCCESS;
    } 
    if (msgLvl(MSG::DEBUG)) msg(MSG::DEBUG) <<  "Found LVL1JetROIs in SG ! " << endmsg;
 
    xAOD::JetRoIContainer::const_iterator itJET  = jetROIs->begin();
    xAOD::JetRoIContainer::const_iterator itJETe = jetROIs->end();
 
 
    int counter = 0;
    for (; itJET != itJETe; ++itJET)
        {
        phi.push_back(DataType( (*itJET)->phi()) );
    eta.push_back(DataType( (*itJET)->eta()) );
 
       if (msgLvl(MSG::DEBUG)) msg(MSG::DEBUG) << "xAOD JetROI #" << counter++ 
          << ", eta: " << (*itJET)->eta() << ", phi: " << (*itJET)->phi() << endmsg;
 
// Placeholders ! No direct access to those in Run-2
// Reference:
//    Event/xAOD/xAODTrigger/trunk/Root/JetRoI_v2.cxx
// jpt 9Jan15:
        energy.push_back(DataType( 0. ));
        roiWord.push_back(DataType( 0. ));
    thrPattern.push_back(DataType( 0. ));
 
// Run-1: info from Alan Watson: pT is getET8x8();
//        energy.push_back(DataType(  itJET->pt()/1000. ));
//      energy.push_back(DataType( (*itJET)->getET4x4()/1000. )); 
//      roiWord.push_back(DataType( (*itJET)->getRoI()));
//  thrPattern.push_back(DataType( (*itJET)->getThrPattern()));
 
      }
 
    DataMap myDataMap;
    const auto nPhi = phi.size();
    myDataMap["phi"] = std::move(phi);
    myDataMap["eta"] = std::move(eta);
    myDataMap["energy"] = std::move(energy);
    myDataMap["roiWord"] = std::move(roiWord);
    myDataMap["thrPattern"] = std::move(thrPattern);
 
    if (msgLvl(MSG::DEBUG)) msg(MSG::DEBUG) << dataTypeName() << ": "<< nPhi
                        << " from: " << m_sgKey << endmsg;
 
    //forward data to formating tool
    return FormatTool->AddToEvent(dataTypeName(), m_sgKey, &myDataMap);
  }
}