GtHepMCDumper.cc

Go to the documentation of this file.

#include "GtHepMCDumper.h"

#include "Event/GenHeader.h"

#include "HepMC/GenEvent.h"
#include "HepMC/GenVertex.h"
#include "HepMC/GenParticle.h"
using namespace HepMC;

GtHepMCDumper::GtHepMCDumper(const std::string& name, ISvcLocator* pSvcLocator)
    : GaudiAlgorithm(name,pSvcLocator)
{
    declareProperty("Location",m_tesLoc = DayaBay::GenHeaderLocation::Default,
                    "Location in TES of HepMCEvents");
    declareProperty("StandardDumper",m_standardDumper=false,
                    "Use standard HepMC dumper (def=false)");
}
GtHepMCDumper::~GtHepMCDumper()
{
}
StatusCode GtHepMCDumper::initialize()
{
    return this->GaudiAlgorithm::initialize();
}

#include "CLHEP/Units/SystemOfUnits.h"

StatusCode GtHepMCDumper::execute()
{
    DayaBay::GenHeader* header = 0;
    try {
        header = get<DayaBay::GenHeader>(m_tesLoc);
    }
    catch (const GaudiException& eg) {
        fatal() << "Failed to get GenHeader from TES at " << m_tesLoc << endreq;
        return StatusCode::FAILURE;
    }

    const TimeStamp& reftime = header->timeStamp();

    info () << "Generator: " << header->generatorName() 
            << " @ reference time: " << reftime
            << " made:" << endreq;
    if (m_standardDumper) {
        header->event()->print();
        return StatusCode::SUCCESS;
    }

    GenEvent* event = header->event();
    if (!event) {
        info() << "\nEvent #(unknown), the GenEvent must be pruned"
               << endreq;
        return StatusCode::SUCCESS;
    }
    info() << "\nEvent #" << event->event_number() 
           << " signal_process_id=" << event->signal_process_id() << ", " 
           << event->vertices_size() << " vertices\n";
    GenEvent::vertex_iterator vit, vdone = event->vertices_end();
    for (vit=event->vertices_begin(); vit != vdone; ++vit) {
        GenVertex* vtx = *vit;
        info() << "\t#" << vtx->id() << " bc="<< vtx->barcode() << " @ x=(" 
               << vtx->position().x()/CLHEP::cm << ","
               << vtx->position().y()/CLHEP::cm << ","
               << vtx->position().z()/CLHEP::cm << ","
               << vtx->position().t()/CLHEP::second << ") [cm,cm,cm,second] (abstime="
               << vtx->position().t()/CLHEP::second + (double)reftime << ")\n"
               << "\t " << vtx->particles_in_size() << " IN particles:\n";
        GenVertex::particles_in_const_iterator piit, pidone = vtx->particles_in_const_end();
        for (piit=vtx->particles_in_const_begin(); piit != pidone; ++piit) {
            const GenParticle* part = *piit;
            info() << "\t\tPDGID=" << part->pdg_id() << " P=("
                   << part->momentum().x() << ","
                   << part->momentum().y() << ","
                   << part->momentum().z() << ","
                   << part->momentum().t() << ")\n";
        }
        info() << "\t " << vtx->particles_out_size() << " OUT particles:\n";
        GenVertex::particles_out_const_iterator poit, podone = vtx->particles_out_const_end();
        for (poit=vtx->particles_out_const_begin(); poit != podone; ++poit) {
            const GenParticle* part = *poit;
            info() << "\t\tPDGID=" << part->pdg_id() << " P=("
                   << part->momentum().x() << ","
                   << part->momentum().y() << ","
                   << part->momentum().z() << ","
                   << part->momentum().t() << ")\n";
        }

    } // loop over vertices
    info() << endreq;

    return StatusCode::SUCCESS;
}
StatusCode GtHepMCDumper::finalize()
{
    return this->GaudiAlgorithm::finalize();
}