A base class for primary generation via HepMC object. More...

#include <Doxymodules_eventgenerator.h>

Inheritance diagram for HepMCG4Interface:

Public Member Functions
	HepMCG4Interface ()

virtual	~HepMCG4Interface ()

HepMC::GenEvent *	GetHepMCGenEvent () const

virtual void	GeneratePrimaryVertex (G4Event *anEvent)

Protected Member Functions
virtual G4bool	CheckVertexInsideWorld (const G4ThreeVector &pos) const

void	HepMC2G4 (const HepMC::GenEvent hepmcevt, G4Event g4event)

virtual HepMC::GenEvent *	GenerateHepMCEvent ()

Protected Attributes
HepMC::GenEvent *	hepmcEvent

Detailed Description

A base class for primary generation via HepMC object.

This class is derived from G4VPrimaryGenerator.

Definition at line 58 of file Doxymodules_eventgenerator.h.

Constructor & Destructor Documentation

◆ HepMCG4Interface()

HepMCG4Interface::HepMCG4Interface ( )

Definition at line 43 of file HepMCG4Interface.cc.

  : hepmcEvent(0)
{
}

◆ ~HepMCG4Interface()

HepMCG4Interface::~HepMCG4Interface ( )

virtual

Definition at line 49 of file HepMCG4Interface.cc.

{
  delete hepmcEvent;
}

Member Function Documentation

◆ CheckVertexInsideWorld()

G4bool HepMCG4Interface::CheckVertexInsideWorld ( const G4ThreeVector & pos ) const

protectedvirtual

Definition at line 55 of file HepMCG4Interface.cc.

{
  G4Navigator* navigator= G4TransportationManager::GetTransportationManager()
                                                 -> GetNavigatorForTracking();
 
  G4VPhysicalVolume* world= navigator-> GetWorldVolume();
  G4VSolid* solid= world-> GetLogicalVolume()-> GetSolid();
  EInside qinside= solid-> Inside(pos);
 
  if( qinside != kInside) return false;
  else return true;
}

◆ HepMC2G4()

void HepMCG4Interface::HepMC2G4	(	const HepMC::GenEvent *	hepmcevt,
		G4Event *	g4event
	)

protected

Definition at line 70 of file HepMCG4Interface.cc.

{
  for(HepMC::GenEvent::vertex_const_iterator vitr= hepmcevt->vertices_begin();
      vitr != hepmcevt->vertices_end(); ++vitr ) { // loop for vertex ...
 
    // real vertex?
    G4bool qvtx=false;
    for (HepMC::GenVertex::particle_iterator
           pitr= (*vitr)->particles_begin(HepMC::children);
         pitr != (*vitr)->particles_end(HepMC::children); ++pitr) {
 
      if (!(*pitr)->end_vertex() && (*pitr)->status()==1) {
        qvtx=true;
        break;
      }
    }
    if (!qvtx) continue;
 
    // check world boundary
    HepMC::FourVector pos= (*vitr)-> position();
    G4LorentzVector xvtx(pos.x(), pos.y(), pos.z(), pos.t());
    if (! CheckVertexInsideWorld(xvtx.vect()*mm)) continue;
 
    // create G4PrimaryVertex and associated G4PrimaryParticles
    G4PrimaryVertex* g4vtx=
      new G4PrimaryVertex(xvtx.x()*mm, xvtx.y()*mm, xvtx.z()*mm,
                          xvtx.t()*mm/c_light);
 
    for (HepMC::GenVertex::particle_iterator
           vpitr= (*vitr)->particles_begin(HepMC::children);
         vpitr != (*vitr)->particles_end(HepMC::children); ++vpitr) {
 
      if( (*vpitr)->status() != 1 ) continue;
 
      G4int pdgcode= (*vpitr)-> pdg_id();
      pos= (*vpitr)-> momentum();
      G4LorentzVector p(pos.px(), pos.py(), pos.pz(), pos.e());
      G4PrimaryParticle* g4prim=
        new G4PrimaryParticle(pdgcode, p.x()*GeV, p.y()*GeV, p.z()*GeV);
 
      g4vtx-> SetPrimary(g4prim);
    }
    g4event-> AddPrimaryVertex(g4vtx);
  }
}

◆ GenerateHepMCEvent()

HepMC::GenEvent * HepMCG4Interface::GenerateHepMCEvent ( )

protectedvirtual

Reimplemented in HepMCG4AsciiReader, and HepMCG4PythiaInterface.

Definition at line 118 of file HepMCG4Interface.cc.

{
  HepMC::GenEvent* aevent= new HepMC::GenEvent();
  return aevent;
}

◆ GetHepMCGenEvent()

HepMC::GenEvent * HepMCG4Interface::GetHepMCGenEvent ( ) const

inline

Definition at line 73 of file HepMCG4Interface.hh.

{
  return hepmcEvent;
}

◆ GeneratePrimaryVertex()

void HepMCG4Interface::GeneratePrimaryVertex ( G4Event * anEvent )

virtual

Definition at line 125 of file HepMCG4Interface.cc.

{
  // delete previous event object
  delete hepmcEvent;
 
  // generate next event
  hepmcEvent= GenerateHepMCEvent();
  if(! hepmcEvent) {
    G4cout << "HepMCInterface: no generated particles. run terminated..."
           << G4endl;
    G4RunManager::GetRunManager()-> AbortRun();
    return;
  }
  HepMC2G4(hepmcEvent, anEvent);
}

Member Data Documentation

◆ hepmcEvent

HepMC::GenEvent* HepMCG4Interface::hepmcEvent