B4d::EventAction Class Reference

Event action class. More...

#include <Doxymodules_basic.h>

Inheritance diagram for B4d::EventAction:

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

Private Member Functions
G4THitsMap< G4double > *	GetHitsCollection (G4int hcID, const G4Event *event) const
G4double	GetSum (G4THitsMap< G4double > *hitsMap) const
void	PrintEventStatistics (G4double absoEdep, G4double absoTrackLength, G4double gapEdep, G4double gapTrackLength) const

Private Attributes
G4int	fAbsoEdepHCID = -1
G4int	fGapEdepHCID = -1
G4int	fAbsoTrackLengthHCID = -1
G4int	fGapTrackLengthHCID = -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 174 of file Doxymodules_basic.h.

Constructor & Destructor Documentation

EventAction()

B4d::EventAction::EventAction ( )

default

~EventAction()

B4d::EventAction::~EventAction ( )

overridedefault

Member Function Documentation

BeginOfEventAction()

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

override

Definition at line 100 of file EventAction.cc.

{}

EndOfEventAction()

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

override

Definition at line 105 of file EventAction.cc.

{
   // Get hist collections IDs
  if ( fAbsoEdepHCID == -1 ) {
    fAbsoEdepHCID
      = G4SDManager::GetSDMpointer()->GetCollectionID("Absorber/Edep");
    fGapEdepHCID
      = G4SDManager::GetSDMpointer()->GetCollectionID("Gap/Edep");
    fAbsoTrackLengthHCID
      = G4SDManager::GetSDMpointer()->GetCollectionID("Absorber/TrackLength");
    fGapTrackLengthHCID
      = G4SDManager::GetSDMpointer()->GetCollectionID("Gap/TrackLength");
  }
 
  // Get sum values from hits collections
  //
  auto absoEdep = GetSum(GetHitsCollection(fAbsoEdepHCID, event));
  auto gapEdep = GetSum(GetHitsCollection(fGapEdepHCID, event));
 
  auto absoTrackLength
    = GetSum(GetHitsCollection(fAbsoTrackLengthHCID, event));
  auto gapTrackLength
    = GetSum(GetHitsCollection(fGapTrackLengthHCID, event));
 
  // get analysis manager
  auto analysisManager = G4AnalysisManager::Instance();
 
  // fill histograms
  //
  analysisManager->FillH1(0, absoEdep);
  analysisManager->FillH1(1, gapEdep);
  analysisManager->FillH1(2, absoTrackLength);
  analysisManager->FillH1(3, gapTrackLength);
 
  // fill ntuple
  //
  analysisManager->FillNtupleDColumn(0, absoEdep);
  analysisManager->FillNtupleDColumn(1, gapEdep);
  analysisManager->FillNtupleDColumn(2, absoTrackLength);
  analysisManager->FillNtupleDColumn(3, gapTrackLength);
  analysisManager->AddNtupleRow();
 
  //print per event (modulo n)
  //
  auto eventID = event->GetEventID();
  auto printModulo = G4RunManager::GetRunManager()->GetPrintProgress();
  if ( ( printModulo > 0 ) && ( eventID % printModulo == 0 ) ) {
    PrintEventStatistics(absoEdep, absoTrackLength, gapEdep, gapTrackLength);
    G4cout << "--> End of event: " << eventID << "\n" << G4endl;  
  }
}

GetHitsCollection()

G4THitsMap< G4double > * B4d::EventAction::GetHitsCollection ( G4int  hcID, const G4Event *  event  ) const

private

Definition at line 48 of file EventAction.cc.

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

GetSum()

G4double B4d::EventAction::GetSum ( G4THitsMap< G4double > *  hitsMap ) const

private

Definition at line 67 of file EventAction.cc.

{
  G4double sumValue = 0.;
  for ( auto it : *hitsMap->GetMap() ) {
    // hitsMap->GetMap() returns the map of std::map<G4int, G4double*>
    sumValue += *(it.second);
  }
  return sumValue;
}

PrintEventStatistics()

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

private

Definition at line 79 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

fAbsoEdepHCID

G4int B4d::EventAction::fAbsoEdepHCID = -1

private

Definition at line 65 of file EventAction.hh.

fGapEdepHCID

G4int B4d::EventAction::fGapEdepHCID = -1

private

Definition at line 66 of file EventAction.hh.

fAbsoTrackLengthHCID

G4int B4d::EventAction::fAbsoTrackLengthHCID = -1

private

Definition at line 67 of file EventAction.hh.

fGapTrackLengthHCID

G4int B4d::EventAction::fGapTrackLengthHCID = -1

private

Definition at line 68 of file EventAction.hh.

---

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

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