simcore::hepmc::HepMCReader Class Reference

Reads HepMC event data into a HepMCEvent object. More...

#include <HepMCReader.h>

Public Member Functions
	HepMCReader (std::string &fileName)
	Class constructor.
virtual	~HepMCReader ()=default
	Class destructor.
std::unique_ptr< HepMCEvent >	readNextEvent ()
	Read the next event.

Private Attributes
std::shared_ptr< HepMC3::Reader >	reader_
	The HepMC3 reader.
int	event_counter_ {0}
	Event counter for tracking progress.

Detailed Description

Reads HepMC event data into a HepMCEvent object.

Definition at line 29 of file HepMCReader.h.

Constructor & Destructor Documentation

HepMCReader()

simcore::hepmc::HepMCReader::HepMCReader

(

std::string &

fileName

)

Class constructor.

Parameters

fileName

The input file name.

Definition at line 9 of file HepMCReader.cxx.

                                            {
  ldmx_log(info) << "Opening HepMC file " << filename;
 
  // Try to detect the file format by reading the first non-empty line
  std::ifstream test_file(filename);
  std::string first_line;
  bool found_line = false;
 
  if (test_file.is_open()) {
    // Skip empty lines and find the first non-empty line
    while (std::getline(test_file, first_line)) {
      // Trim whitespace
      first_line.erase(0, first_line.find_first_not_of(" \t\r\n"));
      first_line.erase(first_line.find_last_not_of(" \t\r\n") + 1);
 
      if (!first_line.empty()) {
        found_line = true;
        break;
      }
    }
    test_file.close();
  }
 
  if (found_line) {
    ldmx_log(debug) << "First non-empty line: " << first_line;
 
    // Check if it's HepMC2 format
    if (first_line.find("HepMC::") != std::string::npos ||
        first_line.find("IO_GenEvent") != std::string::npos) {
      ldmx_log(info) << "Detected HepMC2 format, using ReaderAsciiHepMC2";
      reader_ = std::make_shared<HepMC3::ReaderAsciiHepMC2>(filename);
    } else {
      ldmx_log(info) << "Assuming HepMC3 format, using ReaderAscii";
      reader_ = std::make_shared<HepMC3::ReaderAscii>(filename);
    }
  } else {
    ldmx_log(warn)
        << "Could not find non-empty line, defaulting to HepMC3 format";
    reader_ = std::make_shared<HepMC3::ReaderAscii>(filename);
  }
 
  if (reader_->failed()) {
    EXCEPTION_RAISE("BadFile", "Failed to open HepMC file: " + filename);
  }
}

References reader_.

Member Function Documentation

readNextEvent()

std::unique_ptr< HepMCEvent > simcore::hepmc::HepMCReader::readNextEvent

(

)

Read the next event.

Returns

The next HepMC event.

Definition at line 55 of file HepMCReader.cxx.

                                                     {
  // Create a new HepMC3 GenEvent
  auto gen_event = std::make_shared<HepMC3::GenEvent>();
 
  // Read the next event from the file
  if (!reader_->read_event(*gen_event)) {
    ldmx_log(warn) << "No next event was found by the HepMC reader.";
    return nullptr;
  }
 
  // Check if the event is empty
  if (gen_event->particles().empty()) {
    ldmx_log(warn) << "Empty event found by HepMC reader.";
    return nullptr;
  }
 
  // Increment event counter and log progress every 100 events
  event_counter_++;
  if (event_counter_ % 100 == 0) {
    ldmx_log(debug) << "HepMCReader: Read " << event_counter_ << " events";
  }
 
  // Create and return the wrapped event
  return std::make_unique<HepMCEvent>(gen_event);
}

References event_counter_, and reader_.

Referenced by simcore::generators::HepMCPrimaryGenerator::GeneratePrimaryVertex().

Member Data Documentation

event_counter_

int simcore::hepmc::HepMCReader::event_counter_ {0}

mutableprivate

Event counter for tracking progress.

Definition at line 57 of file HepMCReader.h.

{0};

Referenced by readNextEvent().

reader_

std::shared_ptr<HepMC3::Reader> simcore::hepmc::HepMCReader::reader_

private

The HepMC3 reader.

Definition at line 52 of file HepMCReader.h.

Referenced by HepMCReader(), and readNextEvent().

---

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

• SimCore/include/SimCore/HepMC/HepMCReader.h
• SimCore/src/SimCore/HepMC/HepMCReader.cxx