HadrontherapyLet Class Reference

#include <HadrontherapyLet.hh>

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

Detailed Description

Definition at line 59 of file HadrontherapyLet.hh.

Constructor & Destructor Documentation

HadrontherapyLet::~HadrontherapyLet

(

)

Definition at line 76 of file HadrontherapyLet.cc.

References Clear().

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

Member Function Documentation

void HadrontherapyLet::Clear

(

)

Clear all stored data

Definition at line 92 of file HadrontherapyLet.cc.

Referenced by Initialize(), and ~HadrontherapyLet().

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

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

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

Definition at line 101 of file HadrontherapyLet.cc.

References doCalculation, ionLet::letDD, ionLet::letDN, and test::v.

{
    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;
    }
}

HadrontherapyLet * HadrontherapyLet::GetInstance ( HadrontherapyDetectorConstruction *  pDet )

static

Definition at line 40 of file HadrontherapyLet.cc.

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

HadrontherapyLet * HadrontherapyLet::GetInstance ( )

static

Definition at line 47 of file HadrontherapyLet.cc.

Referenced by main(), HadrontherapyDetectorSD::ProcessHits(), HadrontherapyAnalysisFileMessenger::SetNewValue(), and HadrontherapyDetectorConstruction::UpdateGeometry().

{
    return instance;
}

void HadrontherapyLet::Initialize

(

)

Definition at line 84 of file HadrontherapyLet.cc.

References Clear(), and test::v.

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

void HadrontherapyLet::LetOutput

(

)

Definition at line 166 of file HadrontherapyLet.cc.

References sort(), and test::v.

{
    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

void HadrontherapyLet::StoreLetAscii

(

)

Definition at line 184 of file HadrontherapyLet.cc.

References test::a, G4cout, G4endl, python.hepunit::keV, left, test::v, and width.

{
#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;
        }

    }
}

void HadrontherapyLet::StoreLetRoot

(

)

Definition at line 235 of file HadrontherapyLet.cc.

References HadrontherapyAnalysisManager::GetInstance(), and test::v.

{
#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
}

Field Documentation

G4bool HadrontherapyLet::doCalculation = false

static

Definition at line 68 of file HadrontherapyLet.hh.

Referenced by FillEnergySpectrum().

---

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

• HadrontherapyLet.hh
• HadrontherapyLet.cc