GammaKnifeController Class Reference

#include <GammaKnifeController.hh>

Collaboration diagram for GammaKnifeController:

Public Member Functions
	GammaKnifeController (GammaKnifeDetectorConstruction *)
	~GammaKnifeController ()
void	BeamOn (G4int)
void	ReadFile (std::string fileName)

Private Member Functions
void	StoreHits ()
void	PrepareHitsAccumulation ()
void	AccumulateAllHits ()
void	RotateForward (G4int)
void	RotateBack (G4int)

Private Attributes
std::vector< G4double >	phiAngles
std::vector< G4double >	thetaAngles
GammaKnifeDetectorConstruction *	detector
GammaKnifeMessenger *	messenger
MeshScoreMap *	scoreMaps

Detailed Description

Definition at line 44 of file GammaKnifeController.hh.

Constructor & Destructor Documentation

GammaKnifeController()

GammaKnifeController::GammaKnifeController

(

GammaKnifeDetectorConstruction *

det

)

Definition at line 35 of file GammaKnifeController.cc.

                                                                                : scoreMaps( 0 )
{
    detector = det;

    messenger = new GammaKnifeMessenger( this );
}

~GammaKnifeController()

GammaKnifeController::~GammaKnifeController

(

)

Definition at line 42 of file GammaKnifeController.cc.

{
    delete messenger;
}

Member Function Documentation

AccumulateAllHits()

void GammaKnifeController::AccumulateAllHits ( )

private

Definition at line 140 of file GammaKnifeController.cc.

{
    G4ScoringManager* scm = G4ScoringManager::GetScoringManagerIfExist();
    if (scm)
    {
        for (size_t i = 0; i < scm->GetNumberOfMesh(); i++)
        {
            G4VScoringMesh* mesh = scm->GetMesh(i);
            MeshScoreMap& storedScoreMap = scoreMaps[i];
            MeshScoreMap::iterator it = storedScoreMap.begin();
            for( ; it != storedScoreMap.end(); it++)
            {
                mesh->Accumulate( it->second );
            }
        }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

BeamOn()

void GammaKnifeController::BeamOn

(

G4int

n_event

)

Definition at line 47 of file GammaKnifeController.cc.

{
    PrepareHitsAccumulation();
    for (G4int i = 0; i < GAMMAKNIFE_SOURCES; i++)
    {
        RotateForward(i);
    G4RunManager::GetRunManager()->BeamOn(n_event);

        if (i != (GAMMAKNIFE_SOURCES - 1))        
      StoreHits();
        
        RotateBack(i);
    }
    AccumulateAllHits();
}

Here is the call graph for this function:

Here is the caller graph for this function:

PrepareHitsAccumulation()

void GammaKnifeController::PrepareHitsAccumulation ( )

private

Definition at line 92 of file GammaKnifeController.cc.

{
    G4ScoringManager* scm = G4ScoringManager::GetScoringManagerIfExist();
    if (scm)
    {
        size_t size = scm->GetNumberOfMesh();
        scoreMaps = new MeshScoreMap[size];
        for (size_t i = 0; i < size; i++)
        {
            G4VScoringMesh * mesh = scm->GetMesh(i);

            MeshScoreMap scoreMap = mesh->GetScoreMap();
            MeshScoreMap& storedScoreMap = scoreMaps[i];

            MeshScoreMap::iterator it = scoreMap.begin();
            for( ; it != scoreMap.end(); it++)
            {
                std::string hitMapName = it->first;
                G4THitsMap<G4double>* hitMapToStore = new G4THitsMap<G4double>("GammaKnifeController", hitMapName);
                storedScoreMap[ hitMapName ] = hitMapToStore;
            }
        }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

ReadFile()

void GammaKnifeController::ReadFile

(

std::string

fileName

)

Definition at line 158 of file GammaKnifeController.cc.

{
  //G4cout << "Enter ReadFile()...";
  const int SZ = 100;
  char buf[SZ];
  
  phiAngles.clear();    // If called for the second time
  thetaAngles.clear();  // we won't have 402 positions...
  
  std::ifstream ifs;
  ifs.open( fileName.c_str() );
  
  for (G4int i = 0; i < GAMMAKNIFE_SOURCES; i++)
    {
      G4double phi, theta;
      
      /* Skip the "Axx" at the beginning of the line */
      for (G4int c = 0; c < 4; c++) ifs.get();
      
      ifs >> phi >> theta;
      ifs.getline(buf, SZ); // Next line
      
      phiAngles.push_back( phi * degree );
      thetaAngles.push_back( theta * degree );
    }
  ifs.close();
  //G4cout << "... done " << G4endl;
}

Here is the caller graph for this function:

RotateBack()

void GammaKnifeController::RotateBack ( G4int  position )

private

Definition at line 78 of file GammaKnifeController.cc.

{
    G4ScoringManager* scm = G4ScoringManager::GetScoringManagerIfExist();
    if (scm)
    {
        for (size_t i = 0; i < scm->GetNumberOfMesh(); i++)
        {
            G4VScoringMesh * mesh = scm->GetMesh(i);
            mesh->RotateZ( - phiAngles[position] );
            mesh->RotateX( - thetaAngles[position] );
        }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

RotateForward()

void GammaKnifeController::RotateForward ( G4int  position )

private

Definition at line 63 of file GammaKnifeController.cc.

{
    // Rotate all scoring meshes to the right position
    G4ScoringManager* scm = G4ScoringManager::GetScoringManagerIfExist();
    if (scm)
    {
        for (size_t i = 0; i < scm->GetNumberOfMesh(); i++)
        {
            G4VScoringMesh * mesh = scm->GetMesh(i);
            mesh->RotateX( thetaAngles[position] );
            mesh->RotateZ( phiAngles[position] );
        }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

StoreHits()

void GammaKnifeController::StoreHits ( )

private

Definition at line 117 of file GammaKnifeController.cc.

{
    G4ScoringManager* scm = G4ScoringManager::GetScoringManagerIfExist();
    if (scm)
    {
        for (size_t i = 0; i < scm->GetNumberOfMesh(); i++)
        {
            G4VScoringMesh* mesh = scm->GetMesh(i);


            MeshScoreMap scoreMap = mesh->GetScoreMap();
            MeshScoreMap& storedScoreMap = scoreMaps[i];

            MeshScoreMap::iterator it = scoreMap.begin();
            for( ; it != scoreMap.end(); it++)
            {
                std::string hitMapName = it->first;
                *storedScoreMap[hitMapName] += *(it->second);
            }
        }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

detector

GammaKnifeDetectorConstruction* GammaKnifeController::detector

private

Definition at line 72 of file GammaKnifeController.hh.

messenger

GammaKnifeMessenger* GammaKnifeController::messenger

private

Definition at line 74 of file GammaKnifeController.hh.

phiAngles

std::vector<G4double> GammaKnifeController::phiAngles

private

Definition at line 68 of file GammaKnifeController.hh.

scoreMaps

MeshScoreMap* GammaKnifeController::scoreMaps

private

Definition at line 76 of file GammaKnifeController.hh.

thetaAngles

std::vector<G4double> GammaKnifeController::thetaAngles

private

Definition at line 70 of file GammaKnifeController.hh.

---

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

• GammaKnifeController.hh
• GammaKnifeController.cc