AthenaHepMCtoG4EventAction Class Reference

#include <AthenaHepMCtoG4EventAction.h>

Inheritance diagram for AthenaHepMCtoG4EventAction:

Public Member Functions
	AthenaHepMCtoG4EventAction ()
	AthenaHepMCtoG4EventAction (int)
	~AthenaHepMCtoG4EventAction ()
void	GeneratePrimaries (G4Event *anEvent)
void	SetHepMCEvent (McEventCol::iterator)

Detailed Description

Definition at line 28 of file AthenaHepMCtoG4EventAction.h.

Constructor & Destructor Documentation

AthenaHepMCtoG4EventAction() [1/2]

AthenaHepMCtoG4EventAction::AthenaHepMCtoG4EventAction

(

)

Definition at line 37 of file AthenaHepMCtoG4EventAction.cpp.

                                                      :
m_logLevel( int(MSG::INFO) ) {
}

AthenaHepMCtoG4EventAction() [2/2]

AthenaHepMCtoG4EventAction::AthenaHepMCtoG4EventAction

(

int

level

)

Definition at line 33 of file AthenaHepMCtoG4EventAction.cpp.

                                                                :
  m_logLevel(level) {
}

~AthenaHepMCtoG4EventAction()

AthenaHepMCtoG4EventAction::~AthenaHepMCtoG4EventAction

(

)

Definition at line 41 of file AthenaHepMCtoG4EventAction.cpp.

                                                        {
 
}

Member Function Documentation

GeneratePrimaries()

void AthenaHepMCtoG4EventAction::GeneratePrimaries

(

G4Event *

anEvent

)

Definition at line 152 of file AthenaHepMCtoG4EventAction.cpp.

{
  if ( m_logLevel <= int(MSG::DEBUG) ) {
    std::cout << "AthenaHepMCtoG4EventAction: - Generating Primaries" 
          << std::endl;
  }
 
  
  p_evt = (const HepMC::GenEvent *)(*p_evtCollItr);
 
  
  if (p_evt == 0) { 
 
    if ( m_logLevel <= int(MSG::WARNING) ) {
      std::cout << "AthenaHepMCtoG4EventAction: - HepMC ptr == null" 
        << std::endl;
    }
 
    return; 
 
  }
  
  G4PrimaryVertex* vertex = new G4PrimaryVertex(0.,0.,0.,0.);
  
  int n=0;
  for ( HepMC::GenEvent::vertex_const_iterator v = p_evt->vertices_begin();
        v != p_evt->vertices_end(); ++v ) {
    bool IsGood=false;
    for (HepMC::GenVertex::particle_iterator it=
       (*v)->particles_begin(HepMC::children);
     it!=(*v)->particles_end(HepMC::children);
     it++) {
      
      if (!(*it)->end_vertex() && (*it)->status()==1 ) { 
    IsGood=true;
    break;
      }
    }
 
    if (IsGood) {
//        std::cout << "found a good vertex, barcode: " << (*v)->barcode() 
//          << std::endl;
      
      for (HepMC::GenVertex::particle_iterator it=
         (*v)->particles_begin(HepMC::children);
       it!=(*v)->particles_end(HepMC::children);
       it++) {
    
        if ( (*it)->status() != 1) continue;
 
    int pdgcode=(*it)->pdg_id();
//      GenVertex *dcy = (*it)->end_vertex();
//      std::cout << "found a good particle, barcode: " << (*it)->barcode()
//            << "  pdg: " << pdgcode << "  status: " << (*it)->status()
//            << "  end: ";
//      if (dcy != 0) {
//        std::cout << dcy->barcode();
//      } else {
//        std::cout << "nil";
//      }
//      std::cout << std::endl;
 
    const G4ThreeVector& p=(*it)->momentum();
    G4PrimaryParticle *part;
    if (pdgcode!=999) {
      part=new G4PrimaryParticle (pdgcode,p.x()*GeV,p.y()*GeV,p.z()*GeV);
    } else {
      G4ParticleDefinition *pdef = G4Geantino::GeantinoDefinition();
      part=new G4PrimaryParticle (pdef,p.x()*GeV,p.y()*GeV,p.z()*GeV);
    }
 
      vertex->SetPosition((*v)->position().x(),
                          (*v)->position().y(),
                          (*v)->position().z());
      vertex->SetT0((*v)->position().t());
 
    vertex->SetPrimary(part);
 
    n++;
      }
    }
  }       
  anEvent->AddPrimaryVertex(vertex);
 
  if ( m_logLevel <= int(MSG::DEBUG) ) {
    std::cout << "AthenaHepMCtoG4EventAction: - done. G4Event has " 
          << n << "  particles. "
          << std::endl;
  }
 
  p_evtCollItr++;
  
}

SetHepMCEvent()

void AthenaHepMCtoG4EventAction::SetHepMCEvent

(

McEventCol::iterator

p

)

Definition at line 48 of file AthenaHepMCtoG4EventAction.cpp.

                                                                   { 
  
  if ( m_logLevel <= int(MSG::DEBUG) ) {
    std::cout << "AthenaHepMCtoG4EventAction: - Setting HepMC Event ptr" 
          << std::endl;
  }
 
  //  p_evt = p;
  p_evtCollItr = p;
 
}

---

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

• AthenaHepMCtoG4EventAction.h
• AthenaHepMCtoG4EventAction.cpp