////////////////////////////////////////////////////////////
// $Id: IG4PrimaryParticle.cxx,v 1.7 2008/09/04 22:34:25 mcgrew Exp $
//
#include <TClass.h>
#include <TStreamerInfo.h>
#include <TDatabasePDG.h>
#include <TParticlePDG.h>

#include "IG4PrimaryParticle.hxx"
#include "ICOMETLog.hxx"

ClassImp(COMET::IG4PrimaryParticle);

COMET::IG4PrimaryParticle::IG4PrimaryParticle(void)
    : fTrackId(-1), fParentId(-1), fPDGCode(0), 
      fMomentum(0,0,0,0), IDatum("primary_particle") {}

COMET::IG4PrimaryParticle::~IG4PrimaryParticle() {}

const TParticlePDG* COMET::IG4PrimaryParticle::GetParticle(void) const {
    TParticlePDG* particle = TDatabasePDG::Instance()->GetParticle(fPDGCode);
    if (!particle) {
        COMETLog("IG4PrimaryParticle:: Warning -- Undefined PDG code: "
                 << fPDGCode << " [RETURNING NULL]");
    }
    return particle;
}

void COMET::IG4PrimaryParticle::ls(Option_t *opt) const {
    COMET::ls_header(this,opt);
    std::cout << " Id: " << fTrackId
              << " Particle: " << fPDGCode << " (" << GetTitle() << ")"
              << " E: " << fMomentum.E()
              << std::endl;
}