In This Package:

HepEvt2HepMC.cc File Reference

#include "HepEvt2HepMC.h"
#include "GaudiKernel/IssueSeverity.h"
#include "HepMC/GenEvent.h"
#include "HepMC/GenParticle.h"
#include "HepMC/GenVertex.h"
#include "CLHEP/Units/SystemOfUnits.h"
#include <string>
#include <sstream>
#include <cstdio>

Include dependency graph for HepEvt2HepMC.cc:

Go to the source code of this file.

Functions
static FILE *	open_file (std::string filename)
static void	close_file (FILE *fp, std::string filename)
static bool	is_data (const char *line)
static bool	get_line (FILE *fp, std::string &line)
static bool	unpack (const std::string &line, HepMC::GenEvent &event)

---

Function Documentation

static FILE* open_file

(

std::string

filename

)

[static]

Definition at line 42 of file HepEvt2HepMC.cc.

{
    size_t siz = filename.size();
    if (filename[siz-1] == '|') {
        return popen(filename.substr(0,siz-1).c_str(), "r");
    }
    else {
        return fopen(filename.c_str(), "r");
    }
    return 0;
}

static void close_file	(	FILE *	fp,
		std::string	filename
	)			[static]

Definition at line 55 of file HepEvt2HepMC.cc.

{
    size_t siz = filename.size();
    if (filename[siz-1] == '|') {
        pclose(fp);
    }
    else {
        fclose(fp);
    }
}

static bool is_data

(

const char *

line

)

[static]

Definition at line 66 of file HepEvt2HepMC.cc.

{
    int numchkindx = 0;      // C/C++ generator output
    if (line[0] == ' ') numchkindx = 1; // Fortran output
    return line[numchkindx] >= '0' && line[numchkindx] <= '9';
}

static bool get_line	(	FILE *	fp,
		std::string &	line
	)			[static]

Definition at line 73 of file HepEvt2HepMC.cc.

{
    char buf[1024] = {'\0'};
    if (!fgets(buf,1024,fp)) {
        return false; // EOF or error
    }
    if (buf[1022]) {
        return false;          // Too large line, treat as error
    }
                                           // Non data line, keep fishing
    if (!is_data(buf)) return get_line(fp,line);

    line = buf;
    return true;
}

static bool unpack	(	const std::string &	line,
		HepMC::GenEvent &	event
	)			[static]

Definition at line 89 of file HepEvt2HepMC.cc.

{
    int status=0;     // status code
    int pdgid=0;      // HEP PDG code
    int oldest_daughter=0;    // first daughter
    int baby_daughter=0;    // last daughter
    double px=0.0; // px in GeV
    double py=0.0; // py in GeV
    double pz=0.0; // pz in GeV
    double req_novalue=1e30; // larger than the diameter of univers in mm
    double mass=req_novalue; // mass in GeV
    double vertexX=req_novalue; // x vertex in mm requested on this line
    double vertexY=req_novalue; // y vertex in mm "
    double vertexZ=req_novalue; // z vertex in mm "
    double vertex_dT=req_novalue; // vertex _delta_ time (in ns)
    double polx=0.0;  // x polarization
    double poly=0.0;  // y polarization
    double polz=0.0;  // z polarization
    double energy_unit= CLHEP::GeV;
    double position_unit= CLHEP::millimeter;
    double time_unit= CLHEP::nanosecond; /* used to be mm/c_light */

    std::istringstream is;
    is.str(line);
    is >> status >> pdgid >> oldest_daughter >> baby_daughter
       >> px >> py >> pz >> mass
       >> vertex_dT >> vertexX >> vertexY >> vertexZ // note order, T first
       >> polx >> poly >> polz;
        
    
    HepMC::GenParticle* particle = new HepMC::GenParticle
        (HepMC::FourVector(px*energy_unit,
                          py*energy_unit,
                          pz*energy_unit,
                          sqrt(px*px+py*py+pz*pz+mass*mass)*energy_unit),
         pdgid,status);
    HepMC::GenVertex* vertex = new HepMC::GenVertex
        (HepMC::FourVector(vertexX*position_unit,
                          vertexY*position_unit,
                          vertexZ*position_unit,
                          vertex_dT*time_unit));
    event.add_vertex(vertex);
    if (status == 1) {
        vertex->add_particle_out(particle);
        event.set_signal_process_vertex(vertex);
    }
    else {
        vertex->add_particle_in(particle);
    }
    return true;
}

---