HadrontherapyLet Class Reference

#include <HadrontherapyLet.hh>

Collaboration diagram for HadrontherapyLet:

Public Member Functions
	~HadrontherapyLet ()
void	Initialize ()
void	Clear ()
void	Fill (G4int i, G4int j, G4int k, G4double DE, G4double DX)
void	FillEnergySpectrum (G4int trackID, G4ParticleDefinition *particleDef, G4double DE, G4double DX, G4int i, G4int j, G4int k)
void	LetOutput ()
void	StoreLetAscii ()
void	StoreLetRoot ()

Static Public Member Functions
static HadrontherapyLet *	GetInstance (HadrontherapyDetectorConstruction *)
static HadrontherapyLet *	GetInstance ()

Static Public Attributes
static G4bool	doCalculation = false

Private Member Functions
	HadrontherapyLet (HadrontherapyDetectorConstruction *)

Private Attributes
HadrontherapyPrimaryGeneratorAction *	pPGA
G4Material *	detectorMat
G4double	density
G4String	filename
std::ofstream	ofs
std::ofstream	stopFile
HadrontherapyMatrix *	matrix
G4int	nVoxels
G4int	numberOfVoxelAlongX
G4int	numberOfVoxelAlongY
G4int	numberOfVoxelAlongZ
G4double	primaryEnergy
G4double	energyLimit
G4double	binWidth
G4int	nBins
G4double	nT
G4double	dT
G4double	nD
G4double	dD
G4double	nSecondaryT
G4double	nSecondaryD
G4double	dSecondaryT
G4double	dSecondaryD
G4double	nPrimaryT
G4double	nPrimaryD
G4double	dPrimaryT
G4double	dPrimaryD
G4double *	secondaryLetT
G4double *	secondaryLetD
G4double *	totalLetT
G4double *	DtotalLetD
G4double *	totalLetD
G4String	nome_file
std::vector< ionLet >	ionLetStore

Static Private Attributes
static HadrontherapyLet *	instance = NULL

Detailed Description

Definition at line 61 of file HadrontherapyLet.hh.

Constructor & Destructor Documentation

HadrontherapyLet()

HadrontherapyLet::HadrontherapyLet ( HadrontherapyDetectorConstruction *  pDet )

private

Definition at line 54 of file HadrontherapyLet.cc.

  :filename("Let.out"),matrix(0) // Default output filename
{
    
    matrix = HadrontherapyMatrix::GetInstance();
    
    if (!matrix) 
      G4Exception("HadrontherapyLet::HadrontherapyLet",
          "Hadrontherapy0005", FatalException, 
          "HadrontherapyMatrix not found. Firstly create an instance of it.");
    
    nVoxels = matrix -> GetNvoxel();
    
    numberOfVoxelAlongX = matrix -> GetNumberOfVoxelAlongX();
    numberOfVoxelAlongY = matrix -> GetNumberOfVoxelAlongY();
    numberOfVoxelAlongZ = matrix -> GetNumberOfVoxelAlongZ();
    
    G4RunManager *runManager = G4RunManager::GetRunManager();
    pPGA = (HadrontherapyPrimaryGeneratorAction*)runManager -> GetUserPrimaryGeneratorAction();
    // Pointer to the detector material
    detectorMat = pDet -> GetDetectorLogicalVolume() -> GetMaterial();
    density = detectorMat -> GetDensity();
    // Instances for Total LET
    totalLetD =      new G4double[nVoxels];
    DtotalLetD =     new G4double[nVoxels];
    
}

Here is the call graph for this function:

Here is the caller graph for this function:

~HadrontherapyLet()

HadrontherapyLet::~HadrontherapyLet

(

)

Definition at line 82 of file HadrontherapyLet.cc.

{
    Clear();
    delete [] totalLetD;
    delete [] DtotalLetD;
}

Here is the call graph for this function:

Member Function Documentation

Clear()

void HadrontherapyLet::Clear

(

)

Clear all stored data

Definition at line 98 of file HadrontherapyLet.cc.

{
    for (size_t i=0; i < ionLetStore.size(); i++)
    {
        delete [] ionLetStore[i].letDN;
        delete [] ionLetStore[i].letDD;
    }
    ionLetStore.clear();
}

Here is the caller graph for this function:

Fill()

void HadrontherapyLet::Fill	(	G4int	i,
		G4int	j,
		G4int	k,
		G4double	DE,
		G4double	DX
	)

FillEnergySpectrum()

void HadrontherapyLet::FillEnergySpectrum	(	G4int	trackID,
		G4ParticleDefinition *	particleDef,
		G4double	DE,
		G4double	DX,
		G4int	i,
		G4int	j,
		G4int	k
	)

Definition at line 107 of file HadrontherapyLet.cc.

{
    if (DE <= 0. || DX <=0.) return;
    if (!doCalculation) return;
    G4int Z = particleDef -> GetAtomicNumber();
    
    
    G4int PDGencoding = particleDef -> GetPDGEncoding();
    PDGencoding -= PDGencoding%10;
    
    G4int voxel = matrix -> Index(i,j,k);
    // Total LET calculation...
    totalLetD[voxel]  += DE*(DE/DX);
    DtotalLetD[voxel] += DE;
    // Single ion LET
    if (Z>=1)
    {
        // Search for already allocated data...
        size_t l;
        for (l=0; l < ionLetStore.size(); l++)
        {
            if (ionLetStore[l].PDGencoding == PDGencoding)
                if ( ((trackID ==1) && (ionLetStore[l].isPrimary)) || ((trackID !=1) && (!ionLetStore[l].isPrimary)))
                    break;
        }
        
        if (l == ionLetStore.size()) // Just another type of ion/particle for our store...
        {
            
            G4int A = particleDef -> GetAtomicMass();
            
            G4String fullName = particleDef -> GetParticleName();
            G4String name = fullName.substr (0, fullName.find("[") ); // cut excitation energy [x.y]
            
            ionLet ion =
            {
                (trackID == 1) ? true:false, // is it the primary particle?
                PDGencoding,
                fullName,
                name,
                Z,
                A,
                new G4double[nVoxels], // Let Dose Numerator
                new G4double[nVoxels]  // Let Dose Denominator
            };
            
            // Initialize let
            for(G4int v=0; v < nVoxels; v++) ion.letDN[v] = ion.letDD[v] = 0.;
            ionLetStore.push_back(ion);
            //G4cout << "Allocated LET data for " << ion.name << G4endl;
            
        }
        ionLetStore[l].letDN[voxel] += DE*(DE/DX);
        ionLetStore[l].letDD[voxel] += DE;
    }
}

Here is the call graph for this function:

GetInstance() [1/2]

HadrontherapyLet * HadrontherapyLet::GetInstance ( HadrontherapyDetectorConstruction *  pDet )

static

Definition at line 42 of file HadrontherapyLet.cc.

{
    if (instance) delete instance;
    instance = new HadrontherapyLet(pDet);
    return instance;
}

Here is the call graph for this function:

GetInstance() [2/2]

HadrontherapyLet * HadrontherapyLet::GetInstance ( )

static

Definition at line 49 of file HadrontherapyLet.cc.

{
    return instance;
}

Here is the caller graph for this function:

Initialize()

void HadrontherapyLet::Initialize

(

)

Definition at line 90 of file HadrontherapyLet.cc.

{
    for(G4int v=0; v < nVoxels; v++) totalLetD[v] = DtotalLetD[v] = 0.;
    Clear();
}

Here is the call graph for this function:

LetOutput()

void HadrontherapyLet::LetOutput

(

)

Definition at line 172 of file HadrontherapyLet.cc.

{
    for(G4int v=0; v < nVoxels; v++) if (DtotalLetD[v]>0.) totalLetD[v] = totalLetD[v]/DtotalLetD[v];
    // Sort ions by A and then by Z ...
    std::sort(ionLetStore.begin(), ionLetStore.end());
    // Compute Let Track and Let Dose for any single ion
    
    for(G4int v=0; v < nVoxels; v++)
    {
        for (size_t ion=0; ion < ionLetStore.size(); ion++)
        {
            if (ionLetStore[ion].letDD[v] >0.) ionLetStore[ion].letDN[v] = ionLetStore[ion].letDN[v] / ionLetStore[ion].letDD[v];
            
        }// end loop over ions
    }
    
}// end loop over voxels

StoreLetAscii()

void HadrontherapyLet::StoreLetAscii

(

)

Definition at line 190 of file HadrontherapyLet.cc.

{
#define width 15L
    if(ionLetStore.size())
    {
        ofs.open(filename, std::ios::out);
        if (ofs.is_open())
        {
            
            // Write the voxels index and the list of particles/ions
            ofs << std::setprecision(6) << std::left <<
            "i\tj\tk\t";
            ofs <<  std::setw(width) << "LDT";
            for (size_t l=0; l < ionLetStore.size(); l++)
            {
                G4String a = (ionLetStore[l].isPrimary) ? "_1":"";
                ofs << std::setw(width) << ionLetStore[l].name  + a ;
            }
            ofs << G4endl;
            
            // Write data
            for(G4int i = 0; i < numberOfVoxelAlongX; i++)
                for(G4int j = 0; j < numberOfVoxelAlongY; j++)
                    for(G4int k = 0; k < numberOfVoxelAlongZ; k++)
                    {
                        G4int v = matrix -> Index(i, j, k);
                        // row write
                        for (size_t l=0; l < ionLetStore.size(); l++)
                        {
                            // Write only not identically null data lines
                            if(ionLetStore[l].letDN)
                            {
                                ofs << G4endl;
                                ofs << i << '\t' << j << '\t' << k << '\t';
                                
                                ofs << std::setw(width) << totalLetD[v]/(keV/um);
                                for (size_t ll=0; ll < ionLetStore.size(); ll++)
                                {
                                    ofs << std::setw(width) << ionLetStore[ll].letDN[v]/(keV/um) ;
                                }
                                break;
                            }
                        }
                    }
            ofs.close();
            G4cout << "Let is being written to " << filename << G4endl;
        }
        
    }
}

StoreLetRoot()

void HadrontherapyLet::StoreLetRoot

(

)

Definition at line 241 of file HadrontherapyLet.cc.

{
#ifdef G4ANALYSIS_USE_ROOT
    
    HadrontherapyAnalysisManager* analysis = HadrontherapyAnalysisManager::GetInstance();
    
    for(G4int i = 0; i < numberOfVoxelAlongX; i++)
        for(G4int j = 0; j < numberOfVoxelAlongY; j++) 
            for(G4int k = 0; k < numberOfVoxelAlongZ; k++) 
            {
                G4int v = matrix -> Index(i, j, k);
                for (size_t ion=0; ion < ionLetStore.size(); ion++)
                {
                    
                    analysis -> FillLetFragmentTuple( i, j, k, ionLetStore[ion].A, ionLetStore[ion].Z, ionLetStore[ion].letDN[v]);
                    
                    
                }
            }
    
#endif
}

Here is the call graph for this function:

Member Data Documentation

binWidth

G4double HadrontherapyLet::binWidth

private

Definition at line 98 of file HadrontherapyLet.hh.

dD

G4double HadrontherapyLet::dD

private

Definition at line 100 of file HadrontherapyLet.hh.

density

G4double HadrontherapyLet::density

private

Definition at line 91 of file HadrontherapyLet.hh.

detectorMat

G4Material* HadrontherapyLet::detectorMat

private

Definition at line 90 of file HadrontherapyLet.hh.

doCalculation

G4bool HadrontherapyLet::doCalculation = false

static

Definition at line 70 of file HadrontherapyLet.hh.

dPrimaryD

G4double HadrontherapyLet::dPrimaryD

private

Definition at line 102 of file HadrontherapyLet.hh.

dPrimaryT

G4double HadrontherapyLet::dPrimaryT

private

Definition at line 102 of file HadrontherapyLet.hh.

dSecondaryD

G4double HadrontherapyLet::dSecondaryD

private

Definition at line 101 of file HadrontherapyLet.hh.

dSecondaryT

G4double HadrontherapyLet::dSecondaryT

private

Definition at line 101 of file HadrontherapyLet.hh.

dT

G4double HadrontherapyLet::dT

private

Definition at line 100 of file HadrontherapyLet.hh.

DtotalLetD

G4double * HadrontherapyLet::DtotalLetD

private

Definition at line 104 of file HadrontherapyLet.hh.

energyLimit

G4double HadrontherapyLet::energyLimit

private

Definition at line 98 of file HadrontherapyLet.hh.

filename

G4String HadrontherapyLet::filename

private

Definition at line 92 of file HadrontherapyLet.hh.

instance

HadrontherapyLet * HadrontherapyLet::instance = NULL

staticprivate

Definition at line 87 of file HadrontherapyLet.hh.

ionLetStore

std::vector<ionLet> HadrontherapyLet::ionLetStore

private

Definition at line 107 of file HadrontherapyLet.hh.

matrix

HadrontherapyMatrix* HadrontherapyLet::matrix

private

Definition at line 96 of file HadrontherapyLet.hh.

nBins

G4int HadrontherapyLet::nBins

private

Definition at line 99 of file HadrontherapyLet.hh.

nD

G4double HadrontherapyLet::nD

private

Definition at line 100 of file HadrontherapyLet.hh.

nome_file

G4String HadrontherapyLet::nome_file

private

Definition at line 105 of file HadrontherapyLet.hh.

nPrimaryD

G4double HadrontherapyLet::nPrimaryD

private

Definition at line 102 of file HadrontherapyLet.hh.

nPrimaryT

G4double HadrontherapyLet::nPrimaryT

private

Definition at line 102 of file HadrontherapyLet.hh.

nSecondaryD

G4double HadrontherapyLet::nSecondaryD

private

Definition at line 101 of file HadrontherapyLet.hh.

nSecondaryT

G4double HadrontherapyLet::nSecondaryT

private

Definition at line 101 of file HadrontherapyLet.hh.

nT

G4double HadrontherapyLet::nT

private

Definition at line 100 of file HadrontherapyLet.hh.

numberOfVoxelAlongX

G4int HadrontherapyLet::numberOfVoxelAlongX

private

Definition at line 97 of file HadrontherapyLet.hh.

numberOfVoxelAlongY

G4int HadrontherapyLet::numberOfVoxelAlongY

private

Definition at line 97 of file HadrontherapyLet.hh.

numberOfVoxelAlongZ

G4int HadrontherapyLet::numberOfVoxelAlongZ

private

Definition at line 97 of file HadrontherapyLet.hh.

nVoxels

G4int HadrontherapyLet::nVoxels

private

Definition at line 97 of file HadrontherapyLet.hh.

ofs

std::ofstream HadrontherapyLet::ofs

private

Definition at line 94 of file HadrontherapyLet.hh.

pPGA

HadrontherapyPrimaryGeneratorAction* HadrontherapyLet::pPGA

private

Definition at line 88 of file HadrontherapyLet.hh.

primaryEnergy

G4double HadrontherapyLet::primaryEnergy

private

Definition at line 98 of file HadrontherapyLet.hh.

secondaryLetD

G4double * HadrontherapyLet::secondaryLetD

private

Definition at line 104 of file HadrontherapyLet.hh.

secondaryLetT

G4double* HadrontherapyLet::secondaryLetT

private

Definition at line 104 of file HadrontherapyLet.hh.

stopFile

std::ofstream HadrontherapyLet::stopFile

private

Definition at line 95 of file HadrontherapyLet.hh.

totalLetD

G4double * HadrontherapyLet::totalLetD

private

Definition at line 104 of file HadrontherapyLet.hh.

totalLetT

G4double * HadrontherapyLet::totalLetT

private

Definition at line 104 of file HadrontherapyLet.hh.

---

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

• HadrontherapyLet.hh
• HadrontherapyLet.cc