#include <HICentralityDecorationTool.h>

Inheritance diagram for DerivationFramework::HICentralityDecorationTool:

Collaboration diagram for DerivationFramework::HICentralityDecorationTool:

Public Member Functions
	HICentralityDecorationTool (const std::string &type, const std::string &name, const IInterface *parent)

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.

Constructor & Destructor Documentation

◆ HICentralityDecorationTool()

DerivationFramework::HICentralityDecorationTool::HICentralityDecorationTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 13 of file HICentralityDecorationTool.cxx.

14 : base_class(type, name, parent) {

15 }

Member Function Documentation

◆ addBranches()

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

overridevirtual

Definition at line 60 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 ( )

finaloverride

Definition at line 18 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 31 of file HICentralityDecorationTool.h.