BrachyUserScoreWriter Class Reference

#include <BrachyUserScoreWriter.hh>

Inheritance diagram for BrachyUserScoreWriter:

Collaboration diagram for BrachyUserScoreWriter:

Public Member Functions
	BrachyUserScoreWriter ()
virtual	~BrachyUserScoreWriter ()
void	DumpQuantityToFile (const G4String &psName, const G4String &fileName, const G4String &option)
Public Member Functions inherited from G4VScoreWriter
	G4VScoreWriter ()
virtual	~G4VScoreWriter ()
virtual void	DumpAllQuantitiesToFile (const G4String &fileName, const G4String &option)
void	SetScoringMesh (G4VScoringMesh *sm)
void	SetVerboseLevel (G4int vl)

Additional Inherited Members
Protected Member Functions inherited from G4VScoreWriter
G4int	GetIndex (G4int x, G4int y, G4int z) const
Protected Attributes inherited from G4VScoreWriter
G4int	fNMeshSegments [3]
G4VScoringMesh *	fScoringMesh
G4int	verboseLevel

Detailed Description

Definition at line 44 of file BrachyUserScoreWriter.hh.

Constructor & Destructor Documentation

BrachyUserScoreWriter::BrachyUserScoreWriter

(

)

Definition at line 54 of file BrachyUserScoreWriter.cc.

                                            :
G4VScoreWriter() 
{
}

BrachyUserScoreWriter::~BrachyUserScoreWriter ( )

virtual

Definition at line 59 of file BrachyUserScoreWriter.cc.

{;}

Member Function Documentation

void BrachyUserScoreWriter::DumpQuantityToFile ( const G4String &  psName, const G4String &  fileName, const G4String &  option  )

virtual

Reimplemented from G4VScoreWriter.

Definition at line 62 of file BrachyUserScoreWriter.cc.

{
if(verboseLevel > 0) 
  {G4cout << "BrachyUserScorer-defined DumpQuantityToFile() method is invoked."
  << G4endl; 
  }

// change the option string into lowercase to the case-insensitive.
G4String opt = option;
std::transform(opt.begin(), opt.end(), opt.begin(), (int (*)(int))(tolower));

// confirm the option
if(opt.size() == 0) opt = "csv";

// open the file
std::ofstream ofile(fileName);
  
if(!ofile) 
{
   G4cerr << "ERROR : DumpToFile : File open error -> " << fileName << G4endl;
   return;
}
  ofile << "# mesh name: " << fScoringMesh->GetWorldName() << G4endl;

// retrieve the map
MeshScoreMap fSMap = fScoringMesh -> GetScoreMap();
  
MeshScoreMap::const_iterator msMapItr = fSMap.find(psName);
  
if(msMapItr == fSMap.end()) 
  {
   G4cerr << "ERROR : DumpToFile : Unknown quantity, \""<< psName 
   << "\"." << G4endl;
   return;
  }

std::map<G4int, G4StatDouble*> * score = msMapItr -> second-> GetMap();
  
ofile << "# primitive scorer name: " << msMapItr -> first << G4endl;
//
// Write quantity in the ASCII output file and in brachytherapy.root
//
ofile << std::setprecision(16); // for double value with 8 bytes
  
for(int x = 0; x < fNMeshSegments[0]; x++) {
   for(int y = 0; y < fNMeshSegments[1]; y++) {
     for(int z = 0; z < fNMeshSegments[2]; z++){
        G4int numberOfVoxel_x = fNMeshSegments[0];
        G4int numberOfVoxel_y = fNMeshSegments[1];
        G4int numberOfVoxel_z =fNMeshSegments[2];
        // If the voxel width is changed in the macro file, 
        // the voxel width variable must be updated
        G4double voxelWidth = 0.25 *mm;
        //
        G4double xx = ( - numberOfVoxel_x + 1+ 2*x )* voxelWidth/2;
        G4double yy = ( - numberOfVoxel_y + 1+ 2*y )* voxelWidth/2;
        G4double zz = ( - numberOfVoxel_z + 1+ 2*z )* voxelWidth/2;
        G4int idx = GetIndex(x, y, z);
        std::map<G4int, G4StatDouble*>::iterator value = score -> find(idx);

       if (value != score -> end()) 
        {
         // Print in the ASCII output file the information
 
         ofile << xx << "  " << yy << "  " << zz <<"  " 
               <<(value->second->sum_wx())/keV << G4endl;

#ifdef ANALYSIS_USE          
        // Save the same information in the output analysis file
       BrachyAnalysisManager* analysis = BrachyAnalysisManager::GetInstance();
   
       if(zz> -0.125 *mm && zz < 0.125*mm) analysis -> FillH2WithEnergyDeposition(xx,yy, (value->second->sum_wx())/keV);
#endif
}}}} 

ofile << std::setprecision(6);

// Close the output ASCII file
ofile.close();
}

Here is the call graph for this function:

---

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

• source/geant4.10.03.p03/examples/advanced/brachytherapy/include/BrachyUserScoreWriter.hh
• source/geant4.10.03.p03/examples/advanced/brachytherapy/src/BrachyUserScoreWriter.cc