B4c::EventAction Class Reference

Event action class. More...

#include <Doxymodules_basic.h>

Inheritance diagram for B4c::EventAction:

Public Member Functions
	EventAction ()=default
	~EventAction () override=default
void	BeginOfEventAction (const G4Event *event) override
void	EndOfEventAction (const G4Event *event) override

Private Member Functions
CalorHitsCollection *	GetHitsCollection (G4int hcID, const G4Event *event) const
void	PrintEventStatistics (G4double absoEdep, G4double absoTrackLength, G4double gapEdep, G4double gapTrackLength) const

Private Attributes
G4int	fAbsHCID = -1
G4int	fGapHCID = -1

Detailed Description

Event action class.

In EndOfEventAction(), it prints the accumulated quantities of the energy deposit and track lengths of charged particles in Absober and Gap layers stored in the hits collections.

Definition at line 160 of file Doxymodules_basic.h.

Constructor & Destructor Documentation

EventAction()

B4c::EventAction::EventAction ( )

default

~EventAction()

B4c::EventAction::~EventAction ( )

overridedefault

Member Function Documentation

BeginOfEventAction()

void B4c::EventAction::BeginOfEventAction ( const G4Event *  event )

override

Definition at line 89 of file EventAction.cc.

{}

EndOfEventAction()

void B4c::EventAction::EndOfEventAction ( const G4Event *  event )

override

Definition at line 94 of file EventAction.cc.

{
  // Get hits collections IDs (only once)
  if ( fAbsHCID == -1 ) {
    fAbsHCID
      = G4SDManager::GetSDMpointer()->GetCollectionID("AbsorberHitsCollection");
    fGapHCID
      = G4SDManager::GetSDMpointer()->GetCollectionID("GapHitsCollection");
  }
 
  // Get hits collections
  auto absoHC = GetHitsCollection(fAbsHCID, event);
  auto gapHC = GetHitsCollection(fGapHCID, event);
 
  // Get hit with total values
  auto absoHit = (*absoHC)[absoHC->entries()-1];
  auto gapHit = (*gapHC)[gapHC->entries()-1];
 
  // Print per event (modulo n)
  //
  auto eventID = event->GetEventID();
  auto printModulo = G4RunManager::GetRunManager()->GetPrintProgress();
  if ( ( printModulo > 0 ) && ( eventID % printModulo == 0 ) ) {
    PrintEventStatistics(
      absoHit->GetEdep(), absoHit->GetTrackLength(),
      gapHit->GetEdep(), gapHit->GetTrackLength());
    G4cout << "--> End of event: " << eventID << "\n" << G4endl;      
  }
 
  // Fill histograms, ntuple
  //
 
  // get analysis manager
  auto analysisManager = G4AnalysisManager::Instance();
 
  // fill histograms
  analysisManager->FillH1(0, absoHit->GetEdep());
  analysisManager->FillH1(1, gapHit->GetEdep());
  analysisManager->FillH1(2, absoHit->GetTrackLength());
  analysisManager->FillH1(3, gapHit->GetTrackLength());
 
  // fill ntuple
  analysisManager->FillNtupleDColumn(0, absoHit->GetEdep());
  analysisManager->FillNtupleDColumn(1, gapHit->GetEdep());
  analysisManager->FillNtupleDColumn(2, absoHit->GetTrackLength());
  analysisManager->FillNtupleDColumn(3, gapHit->GetTrackLength());
  analysisManager->AddNtupleRow();
}

GetHitsCollection()

CalorHitsCollection * B4c::EventAction::GetHitsCollection ( G4int  hcID, const G4Event *  event  ) const

private

Definition at line 50 of file EventAction.cc.

{
  auto hitsCollection
    = static_cast<CalorHitsCollection*>(
        event->GetHCofThisEvent()->GetHC(hcID));
 
  if ( ! hitsCollection ) {
    G4ExceptionDescription msg;
    msg << "Cannot access hitsCollection ID " << hcID;
    G4Exception("EventAction::GetHitsCollection()",
      "MyCode0003", FatalException, msg);
  }
 
  return hitsCollection;
}

PrintEventStatistics()

void B4c::EventAction::PrintEventStatistics ( G4double  absoEdep, G4double  absoTrackLength, G4double  gapEdep, G4double  gapTrackLength  ) const

private

Definition at line 69 of file EventAction.cc.

{
  // print event statistics
  G4cout
     << "   Absorber: total energy: "
     << std::setw(7) << G4BestUnit(absoEdep, "Energy")
     << "       total track length: "
     << std::setw(7) << G4BestUnit(absoTrackLength, "Length")
     << G4endl
     << "        Gap: total energy: "
     << std::setw(7) << G4BestUnit(gapEdep, "Energy")
     << "       total track length: "
     << std::setw(7) << G4BestUnit(gapTrackLength, "Length")
     << G4endl;
}

Member Data Documentation

fAbsHCID

G4int B4c::EventAction::fAbsHCID = -1

private

Definition at line 65 of file EventAction.hh.

fGapHCID

G4int B4c::EventAction::fGapHCID = -1

private

Definition at line 66 of file EventAction.hh.

---

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

• .doxygen/Doxymodules_basic.h
• basic/B4/B4c/include/EventAction.hh
• basic/B4/B4c/src/EventAction.cc