#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.