B3bRun Class Reference

#include <B3bRun.hh>

Inheritance diagram for B3bRun:

Collaboration diagram for B3bRun:

Public Member Functions
	B3bRun ()
virtual	~B3bRun ()
virtual void	RecordEvent (const G4Event *)
virtual void	Merge (const G4Run *)
G4int	GetNbGoodEvents () const
G4double	GetSumDose () const
Public Member Functions inherited from G4Run
	G4Run ()
virtual	~G4Run ()
G4int	GetRunID () const
G4int	GetNumberOfEvent () const
G4int	GetNumberOfEventToBeProcessed () const
const G4HCtable *	GetHCtable () const
const G4DCtable *	GetDCtable () const
const G4String &	GetRandomNumberStatus () const
void	SetRunID (G4int id)
void	SetNumberOfEventToBeProcessed (G4int n_ev)
void	SetHCtable (G4HCtable *HCtbl)
void	SetDCtable (G4DCtable *DCtbl)
void	SetRandomNumberStatus (G4String &st)
void	StoreEvent (G4Event *evt)
const std::vector< const G4Event * > *	GetEventVector () const

Private Attributes
G4int	fCollID_cryst
G4int	fCollID_patient
G4int	fPrintModulo
G4int	fGoodEvents
G4double	fSumDose

Additional Inherited Members
Protected Attributes inherited from G4Run
G4int	runID
G4int	numberOfEvent
G4int	numberOfEventToBeProcessed
G4HCtable *	HCtable
G4DCtable *	DCtable
G4String	randomNumberStatus
std::vector< const G4Event * > *	eventVector

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 42 of file B3bRun.hh.

Constructor & Destructor Documentation

B3bRun()

B3bRun::B3bRun

(

)

Definition at line 43 of file B3bRun.cc.

 : G4Run(), 
   fCollID_cryst(-1),
   fCollID_patient(-1),
   fPrintModulo(10000),
   fGoodEvents(0),
   fSumDose(0.)
{ }

~B3bRun()

B3bRun::~B3bRun ( )

virtual

Definition at line 54 of file B3bRun.cc.

{ }

Member Function Documentation

GetNbGoodEvents()

G4int B3bRun::GetNbGoodEvents ( ) const

inline

Definition at line 52 of file B3bRun.hh.

{ return fGoodEvents; }

Here is the caller graph for this function:

GetSumDose()

G4double B3bRun::GetSumDose ( ) const

inline

Definition at line 53 of file B3bRun.hh.

{ return fSumDose; }    

Here is the caller graph for this function:

Merge()

void B3bRun::Merge ( const G4Run *  aRun )

virtual

Reimplemented from G4Run.

Definition at line 118 of file B3bRun.cc.

{
  const B3bRun* localRun = static_cast<const B3bRun*>(aRun);
  fGoodEvents += localRun->fGoodEvents;
  fSumDose    += localRun->fSumDose;

  G4Run::Merge(aRun); 
} 

Here is the call graph for this function:

RecordEvent()

void B3bRun::RecordEvent ( const G4Event *  event )

virtual

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

G4int copyNb = (itr->first);

Reimplemented from G4Run.

Definition at line 59 of file B3bRun.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;
   
  G4THitsMap<G4double>* 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++;
  }  
  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++) {
    dose = *(itr->second);
  }
  fSumDose += dose;
  
  G4Run::RecordEvent(event);      
}  

Here is the call graph for this function:

Member Data Documentation

fCollID_cryst

G4int B3bRun::fCollID_cryst

private

Definition at line 56 of file B3bRun.hh.

fCollID_patient

G4int B3bRun::fCollID_patient

private

Definition at line 57 of file B3bRun.hh.

fGoodEvents

G4int B3bRun::fGoodEvents

private

Definition at line 59 of file B3bRun.hh.

fPrintModulo

G4int B3bRun::fPrintModulo

private

Definition at line 58 of file B3bRun.hh.

fSumDose

G4double B3bRun::fSumDose

private

Definition at line 60 of file B3bRun.hh.

---

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

• B3bRun.hh
• B3bRun.cc