ElectronRun Class Reference

#include <ElectronRun.hh>

Inheritance diagram for ElectronRun:

Collaboration diagram for ElectronRun:

Public Member Functions
	ElectronRun ()
virtual	~ElectronRun ()
virtual void	RecordEvent (const G4Event *)
virtual void	Merge (const G4Run *)
void	DumpData (G4String &) 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 Member Functions
void	Print (const std::vector< G4String > &title, const std::map< G4int, std::vector< G4double > > &out, G4String &) const

Private Attributes
std::map< G4int, G4THitsMap< G4double > *>	fMap

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

Definition at line 43 of file ElectronRun.hh.

Constructor & Destructor Documentation

ElectronRun()

ElectronRun::ElectronRun

(

)

Definition at line 40 of file ElectronRun.cc.

: G4Run(), fMap()
{
    fMap[0] = new G4THitsMap<G4double>("MyDetector", "cell flux");
    fMap[1] = new G4THitsMap<G4double>("MyDetector", "e cell flux");
    fMap[2] = new G4THitsMap<G4double>("MyDetector", "population");
    fMap[3] = new G4THitsMap<G4double>("MyDetector", "e population");
}

~ElectronRun()

ElectronRun::~ElectronRun ( )

virtual

Definition at line 51 of file ElectronRun.cc.

{
    // Important to clean up the map
    std::map<G4int, G4THitsMap<G4double>* >::iterator iter = fMap.begin();
    
    while (iter != fMap.end()) {
        delete iter->second;
        iter++;
    }
}

Member Function Documentation

DumpData()

void ElectronRun::DumpData

(

G4String &

outputFileSpec

)

const

Definition at line 95 of file ElectronRun.cc.

{
    // Titles
    std::vector<G4String> title;
    title.push_back("Radius");
    
    // Output map - energy binning on x axis, theta on y
    std::map< G4int, std::vector<G4double> > output;
    
    G4int nThetaBins = 233;
    
    // Energy bins depends on the number of scorers
    G4int nEnergyBins = fMap.size();
    
    G4int i(0), j(0);
    
    // Initialise current to 0 in all bins
    for (i=0; i<nThetaBins; i++) {
        for (j=0; j<nEnergyBins; j++) {
            output[i].push_back(0);
        }
    }
    
    i=0;
    
    // Fill the output map
    std::map< G4int, G4THitsMap<G4double>* >::const_iterator iter = fMap.begin();
    
    while (iter != fMap.end()) {
        G4THitsMap<G4double>* hitMap = iter->second;
        
        title.push_back(hitMap->GetName());
        
        std::map<G4int,G4double*>* myMap = hitMap->GetMap();
        
        for (j=0; j<nThetaBins; j++) {
            G4double* current = (*myMap)[j];
            if (0 != current) output[j][i] = (*current);
        }
        
        i++;
        iter++;
    }
    
    Print(title, output, outputFileSpec);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Merge()

void ElectronRun::Merge ( const G4Run *  aRun )

virtual

Reimplemented from G4Run.

Definition at line 198 of file ElectronRun.cc.

{
    // This method is called at the end of the run for each worker thread.
    // It accumulates the worker's results into global results.
    const ElectronRun* localRun = static_cast<const ElectronRun*>(aRun);
    const std::map< G4int, G4THitsMap<G4double>* >& localMap = localRun->fMap;
    std::map< G4int, G4THitsMap<G4double>* >::const_iterator iter = localMap.begin();
    for ( ; iter != localMap.end() ; ++iter)
        (*(fMap[iter->first])) += (*(iter->second));
    
    // This call lets Geant4 maintain overall summary information.
    G4Run::Merge(aRun);
}

Here is the call graph for this function:

Print()

void ElectronRun::Print ( const std::vector< G4String > &  title, const std::map< G4int, std::vector< G4double > > &  out, G4String &  outputFileSpec  ) const

private

Definition at line 144 of file ElectronRun.cc.

{
    // Print to G4cout and an output file
    std::ofstream outFile(outputFileSpec);
    
    // Print title vector
    std::vector<G4String>::const_iterator titleIter = title.begin();
    
    while (titleIter != title.end()) {
        G4cout << std::setw(8)<<*titleIter<<" ";
        titleIter++;
    }
    
    G4cout<<G4endl;
    
    // Print current data
    std::map< G4int, std::vector<G4double> >::const_iterator iter = myMap.begin();
    
    while (iter != myMap.end()) {
        G4cout << std::setw(8)<<std::setprecision(3)<< iter->first<<" ";
        
        std::vector<G4double>::const_iterator energyBinIter = iter->second.begin();
        
        // The built-in scorers do not automatically account for the area of the
        // cylinder replica rings. We must account for this now by multiplying our result
        // by the ratio of the area of the full cylinder end over the area of the actual
        // scoring ring.
        // In this ratio, PI divides out, as does the width of the scoring rings.
        // Left with just the number of rings squared divided by the ring index plus
        // 1 squared minus ring index squared.
        G4int ringNum = iter->first;
        G4double areaCorrection = 233.*233. /
        ( (ringNum+1)*(ringNum+1) - ringNum*ringNum );
        G4int counter = 0;
        
        while (energyBinIter != iter->second.end()) {
            G4double value = *energyBinIter;
            if (counter < 2) value = value*areaCorrection;
            G4cout << std::setw(10)<<std::setprecision(5)<< value*mm*mm<<" ";
            outFile << value*mm*mm;
            if (counter < 3) outFile <<",";
            counter++;
            energyBinIter++;
        }
        outFile<<G4endl;
        G4cout<<G4endl;
        iter++;
    }
}

Here is the caller graph for this function:

RecordEvent()

void ElectronRun::RecordEvent ( const G4Event *  anEvent )

virtual

Reimplemented from G4Run.

Definition at line 64 of file ElectronRun.cc.

{
    // Get the hits collection
    G4HCofThisEvent* eventHitCollection = anEvent->GetHCofThisEvent();
    
    if (!eventHitCollection) return;
    
    // Update our private fMap
    std::map< G4int, G4THitsMap<G4double>* >::iterator iter = fMap.begin();
    
    while (iter != fMap.end()) {
        G4int id = iter->first;
        
        // Get the hit collection corresponding to "id"
        G4THitsMap<G4double>* eventHitsMap
        = dynamic_cast< G4THitsMap<G4double>* >(eventHitCollection->GetHC(id));
        
        // Expect this to exist
        assert (0 != eventHitsMap);
        
        // Accumulate event data into our G4THitsMap<G4double> map
        *(iter->second) += *eventHitsMap;
        
        iter++;
    }
    
    G4Run::RecordEvent(anEvent);
}

Here is the call graph for this function:

Member Data Documentation

fMap

std::map<G4int, G4THitsMap<G4double>* > ElectronRun::fMap

private

Definition at line 58 of file ElectronRun.hh.

---

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

• ElectronRun.hh
• ElectronRun.cc