B4dEventAction Class Reference

#include <B4dEventAction.hh>

Inheritance diagram for B4dEventAction:

Collaboration diagram for B4dEventAction:

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

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
G4int	fGapEdepHCID
G4int	fAbsoTrackLengthHCID
G4int	fGapTrackLengthHCID

Additional Inherited Members
Protected Attributes inherited from G4UserEventAction
G4EventManager *	fpEventManager

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 45 of file B4dEventAction.hh.

Constructor & Destructor Documentation

B4dEventAction()

B4dEventAction::B4dEventAction

(

)

Definition at line 45 of file B4dEventAction.cc.

 : G4UserEventAction(),
   fAbsoEdepHCID(-1),
   fGapEdepHCID(-1),
   fAbsoTrackLengthHCID(-1),
   fGapTrackLengthHCID(-1)
{}

~B4dEventAction()

B4dEventAction::~B4dEventAction ( )

virtual

Definition at line 55 of file B4dEventAction.cc.

{}

Member Function Documentation

BeginOfEventAction()

void B4dEventAction::BeginOfEventAction ( const G4Event *  event )

virtual

Reimplemented from G4UserEventAction.

Definition at line 113 of file B4dEventAction.cc.

{}

EndOfEventAction()

void B4dEventAction::EndOfEventAction ( const G4Event *  event )

virtual

Reimplemented from G4UserEventAction.

Definition at line 118 of file B4dEventAction.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
  //
  G4double absoEdep = GetSum(GetHitsCollection(fAbsoEdepHCID, event));
  G4double gapEdep = GetSum(GetHitsCollection(fGapEdepHCID, event));

  G4double absoTrackLength 
    = GetSum(GetHitsCollection(fAbsoTrackLengthHCID, event));
  G4double gapTrackLength 
    = GetSum(GetHitsCollection(fGapTrackLengthHCID, event));

  // get analysis manager
  G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();

  // fill histograms
  //  
  analysisManager->FillH1(1, absoEdep);
  analysisManager->FillH1(2, gapEdep);
  analysisManager->FillH1(3, absoTrackLength);
  analysisManager->FillH1(4, 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)
  //
  G4int eventID = event->GetEventID();
  G4int printModulo = G4RunManager::GetRunManager()->GetPrintProgress();
  if ( ( printModulo > 0 ) && ( eventID % printModulo == 0 ) ) {
    G4cout << "---> End of event: " << eventID << G4endl;     
    PrintEventStatistics(absoEdep, absoTrackLength, gapEdep, gapTrackLength);
  }
}  

Here is the call graph for this function:

GetHitsCollection()

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

private

Definition at line 61 of file B4dEventAction.cc.

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

  return hitsCollection;
}    

Here is the call graph for this function:

Here is the caller graph for this function:

GetSum()

G4double B4dEventAction::GetSum ( G4THitsMap< G4double > *  hitsMap ) const

private

Definition at line 80 of file B4dEventAction.cc.

{
  G4double sumValue = 0;
  std::map<G4int, G4double*>::iterator it;
  for ( it = hitsMap->GetMap()->begin(); it != hitsMap->GetMap()->end(); it++) {
    sumValue += *(it->second);
  }
  return sumValue;  
}  

Here is the call graph for this function:

Here is the caller graph for this function:

PrintEventStatistics()

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

private

Definition at line 92 of file B4dEventAction.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;
}

Here is the caller graph for this function:

Member Data Documentation

fAbsoEdepHCID

G4int B4dEventAction::fAbsoEdepHCID

private

Definition at line 63 of file B4dEventAction.hh.

fAbsoTrackLengthHCID

G4int B4dEventAction::fAbsoTrackLengthHCID

private

Definition at line 65 of file B4dEventAction.hh.

fGapEdepHCID

G4int B4dEventAction::fGapEdepHCID

private

Definition at line 64 of file B4dEventAction.hh.

fGapTrackLengthHCID

G4int B4dEventAction::fGapTrackLengthHCID

private

Definition at line 66 of file B4dEventAction.hh.

---

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

• B4dEventAction.hh
• B4dEventAction.cc