G4HEPEvtInterface Class Reference

#include <G4HEPEvtInterface.hh>

Inheritance diagram for G4HEPEvtInterface:

Public Member Functions
	G4HEPEvtInterface (const char *evfile, G4int vl=0)
	~G4HEPEvtInterface ()
void	GeneratePrimaryVertex (G4Event *evt)
Public Member Functions inherited from G4VPrimaryGenerator
	G4VPrimaryGenerator ()
virtual	~G4VPrimaryGenerator ()
G4ThreeVector	GetParticlePosition ()
G4double	GetParticleTime ()
void	SetParticlePosition (G4ThreeVector aPosition)
void	SetParticleTime (G4double aTime)

Additional Inherited Members
Static Public Member Functions inherited from G4VPrimaryGenerator
static G4bool	CheckVertexInsideWorld (const G4ThreeVector &pos)
Protected Attributes inherited from G4VPrimaryGenerator
G4ThreeVector	particle_position
G4double	particle_time

Detailed Description

Definition at line 78 of file G4HEPEvtInterface.hh.

Constructor & Destructor Documentation

G4HEPEvtInterface::G4HEPEvtInterface	(	const char *	evfile,
		G4int	vl = 0
	)

Definition at line 43 of file G4HEPEvtInterface.cc.

References FatalException, G4cout, G4endl, G4Exception(), G4VPrimaryGenerator::particle_position, and G4VPrimaryGenerator::particle_time.

:vLevel(vl)
{
  inputFile.open((char*)evfile);
  if (inputFile.is_open()) {
    fileName = evfile;
    if(vl>0) G4cout << "G4HEPEvtInterface - " << fileName << " is open." << G4endl;
  }
  else {
    G4Exception("G4HEPEvtInterface::G4HEPEvtInterface","Event0201",FatalException,
    "G4HEPEvtInterface:: cannot open file.");
  }
  G4ThreeVector zero;
  particle_position = zero;
  particle_time = 0.0;

}

G4HEPEvtInterface::~G4HEPEvtInterface

(

)

Definition at line 77 of file G4HEPEvtInterface.cc.

{;}

Member Function Documentation

void G4HEPEvtInterface::GeneratePrimaryVertex ( G4Event *  evt )

virtual

Implements G4VPrimaryGenerator.

Definition at line 80 of file G4HEPEvtInterface.cc.

References G4Event::AddPrimaryVertex(), FatalException, G4cout, G4endl, G4Exception(), python.hepunit::GeV, JustWarning, G4VPrimaryGenerator::particle_position, G4VPrimaryGenerator::particle_time, G4PrimaryParticle::SetDaughter(), G4PrimaryParticle::SetMass(), G4PrimaryParticle::SetMomentum(), and G4PrimaryVertex::SetPrimary().

{
  G4int NHEP = 0;  // number of entries
  if (inputFile.is_open()) {
    inputFile >> NHEP;
  }
  else {
    G4Exception("G4HEPEvtInterface::G4HEPEvtInterface","Event0201",FatalException,
    "G4HEPEvtInterface:: cannot open file.");
  }
  if( inputFile.eof() ) 
  {
    G4Exception("G4HEPEvtInterface::GeneratePrimaryVertex","Event0202",
    JustWarning,"End-Of-File : HEPEvt input file");
    return;
  }

  if(vLevel > 0)
  {
    G4cout << "G4HEPEvtInterface - reading " << NHEP << " HEPEvt particles from "
           << fileName << "." << G4endl;
  }
  for( G4int IHEP=0; IHEP<NHEP; IHEP++ )
  {
    G4int ISTHEP;   // status code
    G4int IDHEP;    // PDG code
    G4int JDAHEP1;  // first daughter
    G4int JDAHEP2;  // last daughter
    G4double PHEP1; // px in GeV
    G4double PHEP2; // py in GeV
    G4double PHEP3; // pz in GeV
    G4double PHEP5; // mass in GeV

    inputFile >> ISTHEP >> IDHEP >> JDAHEP1 >> JDAHEP2
       >> PHEP1 >> PHEP2 >> PHEP3 >> PHEP5;
    if( inputFile.eof() ) 
    {
      G4Exception("G4HEPEvtInterface::GeneratePrimaryVertex","Event0202",
        FatalException,"Unexpected End-Of-File : HEPEvt input file");
    }
    if(vLevel > 1)
    {
      G4cout << " " << ISTHEP << " " << IDHEP << " " << JDAHEP1 << " " << JDAHEP2
             << " " << PHEP1 << " " << PHEP2 << " " << PHEP3 << " " << PHEP5
             << G4endl;
    }

    // create G4PrimaryParticle object
    G4PrimaryParticle* particle 
      = new G4PrimaryParticle( IDHEP );
    particle->SetMass( PHEP5*GeV );
    particle->SetMomentum(PHEP1*GeV, PHEP2*GeV, PHEP3*GeV );

    // create G4HEPEvtParticle object
    G4HEPEvtParticle* hepParticle
      = new G4HEPEvtParticle( particle, ISTHEP, JDAHEP1, JDAHEP2 );

    // Store
    HPlist.push_back( hepParticle );
  }

  // check if there is at least one particle
  if( HPlist.size() == 0 ) return; 

  // make connection between daughter particles decayed from 
  // the same mother
  for( size_t i=0; i<HPlist.size(); i++ )
  {
    if( HPlist[i]->GetJDAHEP1() > 0 ) //  it has daughters
    {
      G4int jda1 = HPlist[i]->GetJDAHEP1()-1; // FORTRAN index starts from 1
      G4int jda2 = HPlist[i]->GetJDAHEP2()-1; // but C++ starts from 0.
      G4PrimaryParticle* mother = HPlist[i]->GetTheParticle();
      for( G4int j=jda1; j<=jda2; j++ )
      {
        G4PrimaryParticle* daughter = HPlist[j]->GetTheParticle();
        if(HPlist[j]->GetISTHEP()>0)
        {
          mother->SetDaughter( daughter );
          HPlist[j]->Done();
        }
      }
    }
  }

  // create G4PrimaryVertex object
  G4PrimaryVertex* vertex = new G4PrimaryVertex(particle_position,particle_time);

  // put initial particles to the vertex
  for( size_t ii=0; ii<HPlist.size(); ii++ )
  {
    if( HPlist[ii]->GetISTHEP() > 0 ) // ISTHEP of daughters had been 
                                       // set to negative
    {
      G4PrimaryParticle* initialParticle = HPlist[ii]->GetTheParticle();
      vertex->SetPrimary( initialParticle );
    }
  }

  // clear G4HEPEvtParticles
  //HPlist.clearAndDestroy();
  for(size_t iii=0;iii<HPlist.size();iii++)
  { delete HPlist[iii]; }
  HPlist.clear();

  // Put the vertex to G4Event object
  evt->AddPrimaryVertex( vertex );
}

---

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

• G4HEPEvtInterface.hh
• G4HEPEvtInterface.cc