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 35 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 37 of file AthenaHepMCtoG4EventAction.cpp.

{}

Member Function Documentation

GeneratePrimaries()

void AthenaHepMCtoG4EventAction::GeneratePrimaries

(

G4Event *

anEvent

)

Definition at line 143 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 42 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