B3b::Run Class Reference

Run class. More...

#include <Doxymodules_basic.h>

Inheritance diagram for B3b::Run:

Public Member Functions
	Run ()=default
	~Run () override=default
void	RecordEvent (const G4Event *) override
void	Merge (const G4Run *) override
G4int	GetNbGoodEvents () const
G4double	GetSumDose () const
G4StatAnalysis	GetStatDose () const

Private Attributes
G4int	fCollID_cryst = -1
G4int	fCollID_patient = -1
G4int	fPrintModulo = 10000
G4int	fGoodEvents = 0
G4double	fSumDose = 0.
G4StatAnalysis	fStatDose

Detailed Description

Run class.

In RecordEvent() there is collected information event per event from Hits Collections, and accumulated statistic for the run

Definition at line 105 of file Doxymodules_basic.h.

Constructor & Destructor Documentation

Run()

B3b::Run::Run ( )

default

~Run()

B3b::Run::~Run ( )

overridedefault

Member Function Documentation

RecordEvent()

void B3b::Run::RecordEvent ( const G4Event *  event )

override

G4int copyNb = (itr->first); G4cout << G4endl << " cryst" << copyNb << ": " << edep/keV << " keV ";

G4int copyNb = (itr->first);

Definition at line 45 of file Run.cc.

{
  if ( fCollID_cryst < 0 ) {
   fCollID_cryst
     = G4SDManager::GetSDMpointer()->GetCollectionID("crystal/edep");
   //G4cout << " fCollID_cryst: " << fCollID_cryst << G4endl;
  }
 
  if ( fCollID_patient < 0 ) {
   fCollID_patient
     = G4SDManager::GetSDMpointer()->GetCollectionID("patient/dose");
   //G4cout << " fCollID_patient: " << fCollID_patient << G4endl;
  }
 
  G4int evtNb = event->GetEventID();
 
  if (evtNb%fPrintModulo == 0) {
    G4cout << G4endl << "---> end of event: " << evtNb << G4endl;
  }
 
  //Hits collections
  //
  G4HCofThisEvent* HCE = event->GetHCofThisEvent();
  if(!HCE) return;
 
  //Energy in crystals : identify 'good events'
  //
  const G4double eThreshold = 500*keV;
  G4int nbOfFired = 0;
 
  auto evtMap = static_cast<G4THitsMap<G4double>*>(HCE->GetHC(fCollID_cryst));
 
  std::map<G4int,G4double*>::iterator itr;
  for (itr = evtMap->GetMap()->begin(); itr != evtMap->GetMap()->end(); itr++) {
    G4double edep = *(itr->second);
    if (edep > eThreshold) nbOfFired++;
    ///G4int copyNb  = (itr->first);
    ///G4cout << G4endl << "  cryst" << copyNb << ": " << edep/keV << " keV ";
  }
  if (nbOfFired == 2) fGoodEvents++;
 
  //Dose deposit in patient
  //
  G4double dose = 0.;
 
  evtMap = static_cast<G4THitsMap<G4double>*>(HCE->GetHC(fCollID_patient));
 
  for (itr = evtMap->GetMap()->begin(); itr != evtMap->GetMap()->end(); itr++) {
    ///G4int copyNb  = (itr->first);
    dose = *(itr->second);
  }
  fSumDose += dose;
  fStatDose += dose;
 
  G4Run::RecordEvent(event);
}

Merge()

void B3b::Run::Merge ( const G4Run *  aRun )

override

Definition at line 104 of file Run.cc.

{
  const Run* localRun = static_cast<const Run*>(aRun);
  fGoodEvents += localRun->fGoodEvents;
  fSumDose    += localRun->fSumDose;
  fStatDose   += localRun->fStatDose;
  G4Run::Merge(aRun);
}

GetNbGoodEvents()

G4int B3b::Run::GetNbGoodEvents ( ) const

inline

Definition at line 55 of file Run.hh.

{ return fGoodEvents; }

GetSumDose()

G4double B3b::Run::GetSumDose ( ) const

inline

Definition at line 56 of file Run.hh.

{ return fSumDose; }

GetStatDose()

G4StatAnalysis B3b::Run::GetStatDose ( ) const

inline

Definition at line 57 of file Run.hh.

{ return fStatDose; }

Member Data Documentation

fCollID_cryst

G4int B3b::Run::fCollID_cryst = -1

private

Definition at line 60 of file Run.hh.

fCollID_patient

G4int B3b::Run::fCollID_patient = -1

private

Definition at line 61 of file Run.hh.

fPrintModulo

G4int B3b::Run::fPrintModulo = 10000

private

Definition at line 62 of file Run.hh.

fGoodEvents

G4int B3b::Run::fGoodEvents = 0

private

Definition at line 63 of file Run.hh.

fSumDose

G4double B3b::Run::fSumDose = 0.

private

Definition at line 64 of file Run.hh.

fStatDose

G4StatAnalysis B3b::Run::fStatDose

private

Definition at line 65 of file Run.hh.

---

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

• .doxygen/Doxymodules_basic.h
• basic/B3/B3b/include/Run.hh
• basic/B3/B3b/src/Run.cc