DerivationFramework::HICentralityDecorationTool Class Reference

#include <HICentralityDecorationTool.h>

Inheritance diagram for DerivationFramework::HICentralityDecorationTool:

Collaboration diagram for DerivationFramework::HICentralityDecorationTool:

Public Member Functions
virtual StatusCode	initialize () override final
virtual StatusCode	addBranches (const EventContext &ctx) const override

Private Attributes
SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfoKey { this, "EventInfoKey", "EventInfo", "" }
Gaudi::Property< std::string >	m_centralityDefinitionFile {this, "centralityDefinitionFile", "HIEventUtils/HeavyIonAnalysis2015_centrality_cuts_Gv32_proposed.txt", "File for centrality definitions"}
std::vector< float >	m_centralityPercentiles
std::vector< float >	m_fCalValues

Detailed Description

Definition at line 19 of file HICentralityDecorationTool.h.

Member Function Documentation

addBranches()

StatusCode DerivationFramework::HICentralityDecorationTool::addBranches ( const EventContext & ctx ) const

overridevirtual

Definition at line 56 of file HICentralityDecorationTool.cxx.

  {
    // Load event EventInfo
    SG::ReadHandle<xAOD::EventInfo> eventInfo (m_eventInfoKey, ctx);
 
    // Set up the decorators for centrality
    const static SG::AuxElement::Decorator< float > ecCentralityMin("CentralityMin") ;
    const static SG::AuxElement::Decorator< float > ecCentralityMax("CentralityMax") ;
 
    const static SG::AuxElement::ConstAccessor<float>  acc_FCalEtA("FCalEtA");
    const static SG::AuxElement::ConstAccessor<float>  acc_FCalEtC("FCalEtC");
 
    // Calculate total FCal ET
    float total_fcal_et = (acc_FCalEtA(*eventInfo) + acc_FCalEtC(*eventInfo)) / 1.e6;
 
    float centralityMin = 0.0;
    float centralityMax = 100.0;
    bool foundRange = false;
    for (size_t i = 0; i < m_fCalValues.size(); ++i) {
      if (total_fcal_et < m_fCalValues[i]) {
        centralityMin = m_centralityPercentiles[i];
        foundRange = true;
        break;
      }
      centralityMax = m_centralityPercentiles[i];
    }
    if (!foundRange) {
      // Top possible range
      centralityMin = 0.;
    }
 
    // Decorate eventInfo with centrality values
    ecCentralityMin(*eventInfo) = centralityMin;
    ecCentralityMax(*eventInfo) = centralityMax;
 
    return StatusCode::SUCCESS;
  }

initialize()

StatusCode DerivationFramework::HICentralityDecorationTool::initialize ( )

finaloverridevirtual

Definition at line 14 of file HICentralityDecorationTool.cxx.

                                                    {
 
    ATH_CHECK( m_eventInfoKey.initialize() );
    // Resolve the path to the centrality definition file
    std::string resolvedPath = PathResolver::find_file(m_centralityDefinitionFile, "CALIBPATH");
 
    // Debug message to print the resolved path or indicate failure
    if (resolvedPath.empty()) {
      ATH_MSG_ERROR("Could not find centrality definition file: " << m_centralityDefinitionFile);
      return StatusCode::FAILURE;
    }
 
    std::ifstream infile(resolvedPath);
    if (!infile.is_open()) {
      ATH_MSG_ERROR("Could not open centrality definition file: " << resolvedPath);
      return StatusCode::FAILURE;
    }
 
    std::string line;
    for (int i = 0; i < 3; ++i) {
      std::getline(infile, line); // Skip header lines
    }
 
    while (std::getline(infile, line)) {
      std::istringstream iss(line);
      std::string centileStr;
      float fCal, centile;
 
      if (iss >> centileStr >> std::skipws >> fCal) {
        centileStr.pop_back(); // Remove '%' character
        centile = std::stof(centileStr);
        m_centralityPercentiles.push_back(centile);
        m_fCalValues.push_back(fCal);
      } else {
        ATH_MSG_WARNING("Could not parse line: " << line);
      }
    }
 
    infile.close();
    return StatusCode::SUCCESS;
  }

Member Data Documentation

m_centralityDefinitionFile

Gaudi::Property<std::string> DerivationFramework::HICentralityDecorationTool::m_centralityDefinitionFile {this, "centralityDefinitionFile", "HIEventUtils/HeavyIonAnalysis2015_centrality_cuts_Gv32_proposed.txt", "File for centrality definitions"}

private

Definition at line 32 of file HICentralityDecorationTool.h.

{this, "centralityDefinitionFile", "HIEventUtils/HeavyIonAnalysis2015_centrality_cuts_Gv32_proposed.txt", "File for centrality definitions"};

m_centralityPercentiles

std::vector<float> DerivationFramework::HICentralityDecorationTool::m_centralityPercentiles

private

Definition at line 34 of file HICentralityDecorationTool.h.

m_eventInfoKey

SG::ReadHandleKey<xAOD::EventInfo> DerivationFramework::HICentralityDecorationTool::m_eventInfoKey { this, "EventInfoKey", "EventInfo", "" }

private

Definition at line 31 of file HICentralityDecorationTool.h.

{ this, "EventInfoKey", "EventInfo", "" };

m_fCalValues

std::vector<float> DerivationFramework::HICentralityDecorationTool::m_fCalValues

private

Definition at line 35 of file HICentralityDecorationTool.h.

---

The documentation for this class was generated from the following files:

• HICentralityDecorationTool.h
• HICentralityDecorationTool.cxx