RE01EventAction Class Reference

#include <Doxymodules_runAndEvent.h>

Inheritance diagram for RE01EventAction:

Public Member Functions
	RE01EventAction ()
virtual	~RE01EventAction ()
virtual void	BeginOfEventAction (const G4Event *)
virtual void	EndOfEventAction (const G4Event *)

Private Member Functions
void	PrintPrimary (G4PrimaryParticle *pp, G4int ind)

Private Attributes
G4int	fTrackerCollID
G4int	fCalorimeterCollID

Detailed Description

Definition at line 25 of file Doxymodules_runAndEvent.h.

Constructor & Destructor Documentation

RE01EventAction()

RE01EventAction::RE01EventAction

(

)

Definition at line 50 of file RE01EventAction.cc.

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

~RE01EventAction()

RE01EventAction::~RE01EventAction ( )

virtual

Definition at line 56 of file RE01EventAction.cc.

{;}

Member Function Documentation

BeginOfEventAction()

void RE01EventAction::BeginOfEventAction ( const G4Event *  )

virtual

Definition at line 60 of file RE01EventAction.cc.

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

EndOfEventAction()

void RE01EventAction::EndOfEventAction ( const G4Event *  evt )

virtual

Definition at line 72 of file RE01EventAction.cc.

{
  G4cout << ">>> Summary of Event " << evt->GetEventID() << G4endl;
  if(evt->GetNumberOfPrimaryVertex()==0)
  {
    G4cout << "Event is empty." << G4endl;
    return;
  }
  
  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;
    if(fpEventManager->GetVerboseLevel()>0)
    {
      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;
  if(fpEventManager->GetVerboseLevel()>0)
  {
    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;
    if(fpEventManager->GetVerboseLevel()>0)
    {
      G4PrimaryParticle* pp = pv->GetPrimary();
      while(pp)
      {
        PrintPrimary(pp,0);
        pp = pp->GetNext();
      }
    }
  }
}

PrintPrimary()

void RE01EventAction::PrintPrimary ( G4PrimaryParticle *  pp, G4int  ind  )

private

Definition at line 168 of file RE01EventAction.cc.

{
  for(G4int ii=0;ii<=ind;ii++)
  { G4cout << "  "; }
  G4cout << "==PDGcode " << pp->GetPDGcode() << " ";
  if(pp->GetG4code()!=0)
  { G4cout << "(" << pp->GetG4code()->GetParticleName() << ")"; }
  else
  { G4cout << "is not defined in G4"; }
  G4cout << " " << pp->GetMomentum()/GeV << " [GeV] ";
  if(pp->GetTrackID()<0)
  { G4cout << G4endl; }
  else
  { G4cout << ">>> G4Track ID " << pp->GetTrackID() << G4endl; }
 
  G4PrimaryParticle* daughter = pp->GetDaughter();
  while(daughter)
  {
    PrintPrimary(daughter,ind+1);
    daughter = daughter->GetNext();
  }
}

Member Data Documentation

fTrackerCollID

G4int RE01EventAction::fTrackerCollID

private

Definition at line 50 of file RE01EventAction.hh.

fCalorimeterCollID

G4int RE01EventAction::fCalorimeterCollID

private

Definition at line 51 of file RE01EventAction.hh.

---

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

• .doxygen/Doxymodules_runAndEvent.h
• extended/runAndEvent/RE01/include/RE01EventAction.hh
• extended/runAndEvent/RE01/src/RE01EventAction.cc