B3Run Class Reference

#include <B3Run.hh>

Inheritance diagram for B3Run:

Public Member Functions
	B3Run ()
virtual	~B3Run ()
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

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 B3Run.hh.

Constructor & Destructor Documentation

B3Run::B3Run

(

)

Definition at line 43 of file B3Run.cc.

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

B3Run::~B3Run ( )

virtual

Definition at line 55 of file B3Run.cc.

{ }

Member Function Documentation

G4int B3Run::GetNbGoodEvents ( ) const

inline

Definition at line 52 of file B3Run.hh.

Referenced by B3RunAction::EndOfRunAction().

{ return fGoodEvents; }

G4double B3Run::GetSumDose ( ) const

inline

Definition at line 53 of file B3Run.hh.

Referenced by B3RunAction::EndOfRunAction().

{ return fSumDose; }    

void B3Run::Merge ( const G4Run *  aRun )

virtual

Reimplemented from G4Run.

Definition at line 119 of file B3Run.cc.

References G4Run::Merge().

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

  G4Run::Merge(aRun); 
} 

void B3Run::RecordEvent ( const G4Event *  event )

virtual

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

G4int copyNb = (itr->first);

Reimplemented from G4Run.

Definition at line 60 of file B3Run.cc.

References G4cout, G4endl, G4SDManager::GetCollectionID(), G4HCofThisEvent::GetHC(), G4THitsMap< T >::GetMap(), G4SDManager::GetSDMpointer(), python.hepunit::keV, and G4Run::RecordEvent().

{
  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 << "\n---> 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++;
    ///G4int copyNb  = (itr->first);
    ///G4cout << "\n  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;
  
  G4Run::RecordEvent(event);      
}  

---

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

• B3Run.hh
• B3Run.cc