RE01EventAction Class Reference

#include <RE01EventAction.hh>

Inheritance diagram for RE01EventAction:

Public Member Functions
	RE01EventAction ()
virtual	~RE01EventAction ()
virtual void	BeginOfEventAction (const G4Event *)
virtual void	EndOfEventAction (const G4Event *)
Public Member Functions inherited from G4UserEventAction
	G4UserEventAction ()
virtual	~G4UserEventAction ()
void	SetEventManager (G4EventManager *value)

Additional Inherited Members
Protected Attributes inherited from G4UserEventAction
G4EventManager *	fpEventManager

Detailed Description

Definition at line 40 of file RE01EventAction.hh.

Constructor & Destructor Documentation

RE01EventAction::RE01EventAction

(

)

Definition at line 51 of file RE01EventAction.cc.

  :G4UserEventAction(),
   fTrackerCollID(-1),fCalorimeterCollID(-1),fMuonCollID(-1)
{;}

RE01EventAction::~RE01EventAction ( )

virtual

Definition at line 57 of file RE01EventAction.cc.

{;}

Member Function Documentation

void RE01EventAction::BeginOfEventAction ( const G4Event *  )

virtual

Reimplemented from G4UserEventAction.

Definition at line 61 of file RE01EventAction.cc.

References G4SDManager::GetCollectionID(), and G4SDManager::GetSDMpointer().

{
  G4SDManager * SDman = G4SDManager::GetSDMpointer();
  if(fTrackerCollID<0||fCalorimeterCollID<0||fMuonCollID<0)
  {
    G4String colNam;
    fTrackerCollID = SDman->GetCollectionID(colNam="trackerCollection");
    fCalorimeterCollID = SDman->GetCollectionID(colNam="calCollection");
  }
}

void RE01EventAction::EndOfEventAction ( const G4Event *  evt )

virtual

Reimplemented from G4UserEventAction.

Definition at line 73 of file RE01EventAction.cc.

References G4TrajectoryContainer::entries(), G4THitsCollection< T >::entries(), G4cout, G4endl, G4Event::GetEventID(), G4HCofThisEvent::GetHC(), G4Event::GetHCofThisEvent(), G4Event::GetNumberOfPrimaryVertex(), G4PrimaryVertex::GetPrimary(), G4Event::GetPrimaryVertex(), G4PrimaryVertex::GetT0(), G4Event::GetTrajectoryContainer(), G4PrimaryVertex::GetX0(), G4PrimaryVertex::GetY0(), G4PrimaryVertex::GetZ0(), python.hepunit::GeV, ns, G4InuclParticleNames::pp, and RE01Trajectory::ShowTrajectory().

{
  G4cout << ">>> Summary of Event " << evt->GetEventID() << G4endl;
  
  if(fTrackerCollID<0||fCalorimeterCollID<0) return;

  G4HCofThisEvent * HCE = evt->GetHCofThisEvent();
  RE01TrackerHitsCollection* THC = 0;
  RE01CalorimeterHitsCollection* CHC = 0;
  if(HCE)
  {
    THC = (RE01TrackerHitsCollection*)(HCE->GetHC(fTrackerCollID));
    CHC = (RE01CalorimeterHitsCollection*)(HCE->GetHC(fCalorimeterCollID));
  }

  if(THC)
  {
    int n_hit = THC->entries();
    G4cout << G4endl;
    G4cout << "Tracker hits " <<
      "--------------------------------------------------------------" 
           << G4endl;
    G4cout << n_hit << " hits are stored in RE01TrackerHitsCollection." 
           << G4endl;
    G4cout << "List of hits in tracker" << G4endl;
    for(int i=0;i<n_hit;i++)
    { (*THC)[i]->Print(); }
  }
  if(CHC)
  {
    int n_hit = CHC->entries();
    G4cout << G4endl;
    G4cout << "Calorimeter hits "<<
      "--------------------------------------------------------------" 
           << G4endl;
    G4cout << n_hit << " hits are stored in RE01CalorimeterHitsCollection." 
           << G4endl;
    G4double totE = 0;
    for(int i=0;i<n_hit;i++)
    { totE += (*CHC)[i]->GetEdep(); }
    G4cout << "     Total energy deposition in calorimeter : "
         << totE / GeV << " (GeV)" << G4endl;
  }

  // get number of stored trajectories
  G4TrajectoryContainer* trajectoryContainer = evt->GetTrajectoryContainer();
  G4int n_trajectories = 0;
  if (trajectoryContainer) n_trajectories = trajectoryContainer->entries();
  // extract the trajectories and print them out
  G4cout << G4endl;
  G4cout << "Trajectories in tracker "<<
    "--------------------------------------------------------------" 
         << G4endl;
  for(G4int i=0; i<n_trajectories; i++) 
  {
    RE01Trajectory* trj = 
      (RE01Trajectory*)((*(evt->GetTrajectoryContainer()))[i]);
    trj->ShowTrajectory();
  }
    
  G4cout << G4endl;
  G4cout << "Primary particles "<<
    "--------------------------------------------------------------" 
         << G4endl;
  G4int n_vertex = evt->GetNumberOfPrimaryVertex();
  for(G4int iv=0;iv<n_vertex;iv++)
  {
    G4PrimaryVertex* pv = evt->GetPrimaryVertex(iv);
    G4cout << G4endl;
    G4cout << "Primary vertex "
           << G4ThreeVector(pv->GetX0(),pv->GetY0(),pv->GetZ0())
           << "   at t = " << (pv->GetT0())/ns << " [ns]" << G4endl;
    G4PrimaryParticle* pp = pv->GetPrimary();
    while(pp)
    {
      PrintPrimary(pp,0);
      pp = pp->GetNext();
    }
  }
}

---

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

• RE01EventAction.hh
• RE01EventAction.cc