CML2SDWithVoxels Class Reference

#include <ML2SDWithVoxels.hh>

Inheritance diagram for CML2SDWithVoxels:

Collaboration diagram for CML2SDWithVoxels:

Public Member Functions
	CML2SDWithVoxels (G4String name, G4int saving_in_ROG_Voxels_every_events, G4int seed, G4String ROGOutFile, G4bool bSaveROG, G4ThreeVector centre, G4ThreeVector halfSize, G4int NumberOfVoxelsAlongX, G4int NumberOfVoxelsAlongY, G4int NumberOfVoxelsAlongZ)
	~CML2SDWithVoxels (void)
G4bool	ProcessHits (G4Step *aStep, G4TouchableHistory *ROHist)
void	Initialize (G4HCofThisEvent *)
void	EndOfEvent (G4HCofThisEvent *)
G4int	getTotalNumberOfEvents ()
void	setActive (G4bool act)
void	save ()
void	setRecycling (int recycling)
void	setVolumeNameIdLink (std::vector< SvolumeNameId > volNameIdLink)
void	resetVoxelsSingle ()
void	setFullOutFileDataSingle (G4String val)
Public Member Functions inherited from G4VSensitiveDetector
	G4VSensitiveDetector (G4String name)
	G4VSensitiveDetector (const G4VSensitiveDetector &right)
virtual	~G4VSensitiveDetector ()
G4VSensitiveDetector &	operator= (const G4VSensitiveDetector &right)
G4int	operator== (const G4VSensitiveDetector &right) const
G4int	operator!= (const G4VSensitiveDetector &right) const
virtual void	clear ()
virtual void	DrawAll ()
virtual void	PrintAll ()
G4bool	Hit (G4Step *aStep)
void	SetROgeometry (G4VReadOutGeometry *value)
void	SetFilter (G4VSDFilter *value)
G4int	GetNumberOfCollections () const
G4String	GetCollectionName (G4int id) const
void	SetVerboseLevel (G4int vl)
void	Activate (G4bool activeFlag)
G4bool	isActive () const
G4String	GetName () const
G4String	GetPathName () const
G4String	GetFullPathName () const
G4VReadOutGeometry *	GetROgeometry () const
G4VSDFilter *	GetFilter () const
virtual G4VSensitiveDetector *	Clone () const

Private Member Functions
void	saveData (G4String Filename, Svoxel ***voxels)
G4int	getIdFromVolumeName (G4String name)

Private Attributes
G4ThreeVector	halfSize
G4ThreeVector	centre
G4ThreeVector	pos
G4double	halfXVoxelDimensionX
G4double	halfXVoxelDimensionY
G4double	halfXVoxelDimensionZ
G4int	NumberOfVoxelsAlongX
G4int	NumberOfVoxelsAlongY
G4int	NumberOfVoxelsAlongZ
Svoxel ***	voxelsSum
Svoxel ***	voxelsSingle
G4String	fullOutFileData
G4String	fullOutFileDataSingle
G4int	nTotalEvents
G4int	nSingleTotalEvents
G4int	nParticle
G4int	nParticleValatile
G4int	saving_in_ROG_Voxels_every_events
G4int	nRecycling
G4bool	bActive
G4bool	bSaveROG
G4double	voxelMass
G4double	density
G4double	voxelVolume
std::vector< SvolumeNameId >	volumeNameIdLink

Additional Inherited Members
Protected Member Functions inherited from G4VSensitiveDetector
virtual G4int	GetCollectionID (G4int i)
Protected Attributes inherited from G4VSensitiveDetector
G4CollectionNameVector	collectionName
G4String	SensitiveDetectorName
G4String	thePathName
G4String	fullPathName
G4int	verboseLevel
G4bool	active
G4VReadOutGeometry *	ROgeometry
G4VSDFilter *	filter

Detailed Description

Definition at line 55 of file ML2SDWithVoxels.hh.

Constructor & Destructor Documentation

CML2SDWithVoxels()

CML2SDWithVoxels::CML2SDWithVoxels	(	G4String	name,
		G4int	saving_in_ROG_Voxels_every_events,
		G4int	seed,
		G4String	ROGOutFile,
		G4bool	bSaveROG,
		G4ThreeVector	centre,
		G4ThreeVector	halfSize,
		G4int	NumberOfVoxelsAlongX,
		G4int	NumberOfVoxelsAlongY,
		G4int	NumberOfVoxelsAlongZ
	)

Definition at line 49 of file ML2SDWithVoxels.cc.

: G4VSensitiveDetector(name),voxelsSum(0), voxelsSingle(0)
{
    saving_in_ROG_Voxels_every_events=voxSave;
    bSaveROG=bROG;
    bActive=true;
    nParticle=0;
    nParticleValatile=0;
    nTotalEvents=0;
    nSingleTotalEvents=0;
    density=1.;
    voxelVolume=0.;
    voxelMass=0.;
    nRecycling=1;

    G4String seedName;
    char a[10];
    sprintf(a,"%d", seed);
    seedName=(G4String)a;

    fullOutFileData=ROGOutFile+"_"+seedName+".txt";
    fullOutFileDataSingle="";
    if(bSaveROG)
    {
        centre=ctr;
        halfSize=hSiz;
        NumberOfVoxelsAlongX=NumVX;
        NumberOfVoxelsAlongY=NumVY;
        NumberOfVoxelsAlongZ=NumVZ;
        halfXVoxelDimensionX=halfSize.getX()/NumberOfVoxelsAlongX;
        halfXVoxelDimensionY=halfSize.getY()/NumberOfVoxelsAlongY;
        halfXVoxelDimensionZ=halfSize.getZ()/NumberOfVoxelsAlongZ;

        voxelVolume=halfXVoxelDimensionX*halfXVoxelDimensionY*halfXVoxelDimensionZ*8.;

// voxels to store and save the sum of the geometry configurations
        voxelsSum=new Svoxel**[NumberOfVoxelsAlongX];
        for (int ix=0; ix< NumberOfVoxelsAlongX; ix++)
        {
            voxelsSum[ix]=new Svoxel*[NumberOfVoxelsAlongY];
            for (int iy=0; iy< NumberOfVoxelsAlongY; iy++)
            {
                voxelsSum[ix][iy]=new Svoxel[NumberOfVoxelsAlongZ];
                for (int iz=0; iz< NumberOfVoxelsAlongZ; iz++)
                {
                    voxelsSum[ix][iy][iz].volumeId=-1;
                    voxelsSum[ix][iy][iz].depEnergy=0.;
                    voxelsSum[ix][iy][iz].depEnergy2=0.;
                    voxelsSum[ix][iy][iz].depEnergyNorm=0.;
                    voxelsSum[ix][iy][iz].depEnergyNormError=0.;
                    voxelsSum[ix][iy][iz].expDose=0.;
                    voxelsSum[ix][iy][iz].halfSize.set(halfXVoxelDimensionX, halfXVoxelDimensionY, halfXVoxelDimensionZ);
                    voxelsSum[ix][iy][iz].pos.set(2.*(ix)*halfXVoxelDimensionX  -halfSize.getX()+halfXVoxelDimensionX + centre.getX(), 
            2.*(iy)*halfXVoxelDimensionY  -halfSize.getY()+halfXVoxelDimensionY + centre.getY(), 
            2.*(iz)*halfXVoxelDimensionZ  -halfSize.getZ()+halfXVoxelDimensionZ + centre.getZ());
                    voxelsSum[ix][iy][iz].nEvents=0;
                }
            }
        }


// voxels to store and save the single geometry configuration
        voxelsSingle=new Svoxel**[NumberOfVoxelsAlongX];
        for (int ix=0; ix< NumberOfVoxelsAlongX; ix++)
        {
            voxelsSingle[ix]=new Svoxel*[NumberOfVoxelsAlongY];
            for (int iy=0; iy< NumberOfVoxelsAlongY; iy++)
            {
                voxelsSingle[ix][iy]=new Svoxel[NumberOfVoxelsAlongZ];
                for (int iz=0; iz< NumberOfVoxelsAlongZ; iz++)
                {
                    voxelsSingle[ix][iy][iz].volumeId=-1;
                    voxelsSingle[ix][iy][iz].depEnergy=0.;
                    voxelsSingle[ix][iy][iz].depEnergy2=0.;
                    voxelsSingle[ix][iy][iz].depEnergyNorm=0.;
                    voxelsSingle[ix][iy][iz].depEnergyNormError=0.;
                    voxelsSingle[ix][iy][iz].expDose=0.;
                    voxelsSingle[ix][iy][iz].halfSize.set(halfXVoxelDimensionX, halfXVoxelDimensionY, halfXVoxelDimensionZ);
                    voxelsSingle[ix][iy][iz].pos.set(2.*(ix)*halfXVoxelDimensionX  -halfSize.getX()+halfXVoxelDimensionX + centre.getX(), 
            2.*(iy)*halfXVoxelDimensionY  -halfSize.getY()+halfXVoxelDimensionY + centre.getY(), 
            2.*(iz)*halfXVoxelDimensionZ  -halfSize.getZ()+halfXVoxelDimensionZ + centre.getZ());
                    voxelsSingle[ix][iy][iz].nEvents=0;
                }
            }
        }
    }
}

Here is the call graph for this function:

~CML2SDWithVoxels()

CML2SDWithVoxels::~CML2SDWithVoxels

(

void

)

Definition at line 140 of file ML2SDWithVoxels.cc.

{
    if(bSaveROG)
    {
        delete [] voxelsSum;
        delete [] voxelsSingle;
    }
}

Member Function Documentation

EndOfEvent()

void CML2SDWithVoxels::EndOfEvent ( G4HCofThisEvent *  )

inlinevirtual

Reimplemented from G4VSensitiveDetector.

Definition at line 62 of file ML2SDWithVoxels.hh.

{}

getIdFromVolumeName()

G4int CML2SDWithVoxels::getIdFromVolumeName ( G4String  name )

private

Definition at line 210 of file ML2SDWithVoxels.cc.

{
    for (int i=0; i<(int)volumeNameIdLink.size(); i++)
    {
        if (volumeNameIdLink[i].volumeName==name)
        {
            return volumeNameIdLink[i].volumeId;
            break;
        }
    }
    return -1;
}

Here is the caller graph for this function:

getTotalNumberOfEvents()

G4int CML2SDWithVoxels::getTotalNumberOfEvents ( )

inline

Definition at line 63 of file ML2SDWithVoxels.hh.

{return nTotalEvents;}

Here is the caller graph for this function:

Initialize()

void CML2SDWithVoxels::Initialize ( G4HCofThisEvent *  )

inlinevirtual

Reimplemented from G4VSensitiveDetector.

Definition at line 61 of file ML2SDWithVoxels.hh.

{}

ProcessHits()

G4bool CML2SDWithVoxels::ProcessHits ( G4Step *  aStep, G4TouchableHistory *  ROHist  )

virtual

Implements G4VSensitiveDetector.

Definition at line 172 of file ML2SDWithVoxels.cc.

{
    
    if (bActive)
    {
        G4double energyDep = aStep->GetTotalEnergyDeposit();
        if (bSaveROG && energyDep>0.) 
        {
            G4int ix, iy, iz;
            G4String volumeName;

            ix=ROHist->GetReplicaNumber(2);
            iy=ROHist->GetReplicaNumber(0);
            iz=ROHist->GetReplicaNumber(1);

            density=aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetMaterial()->GetDensity();

            voxelMass=voxelVolume*density;
            energyDep/=voxelMass*nRecycling;
            voxelsSum[ix][iy][iz].volumeId=getIdFromVolumeName(aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetMaterial()->GetName());
            voxelsSum[ix][iy][iz].depEnergy+=energyDep;
            voxelsSum[ix][iy][iz].depEnergy2+=energyDep*energyDep;
            voxelsSum[ix][iy][iz].nEvents++;

            voxelsSingle[ix][iy][iz].volumeId=getIdFromVolumeName(aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetMaterial()->GetName());
            voxelsSingle[ix][iy][iz].depEnergy+=energyDep;
            voxelsSingle[ix][iy][iz].depEnergy2+=energyDep*energyDep;
            voxelsSingle[ix][iy][iz].nEvents++;
            nTotalEvents++;
            nSingleTotalEvents++;
            if (nTotalEvents%saving_in_ROG_Voxels_every_events==0 && nTotalEvents>0)
            {
                save();
            }
        }
    }
    return true;
}

Here is the call graph for this function:

resetVoxelsSingle()

void CML2SDWithVoxels::resetVoxelsSingle

(

)

Definition at line 148 of file ML2SDWithVoxels.cc.

{
    for (int ix=0; ix< NumberOfVoxelsAlongX; ix++)
    {
        for (int iy=0; iy< NumberOfVoxelsAlongY; iy++)
        {
            for (int iz=0; iz< NumberOfVoxelsAlongZ; iz++)
            {
                voxelsSingle[ix][iy][iz].volumeId=-1;
                voxelsSingle[ix][iy][iz].depEnergy=0.;
                voxelsSingle[ix][iy][iz].depEnergy2=0.;
                voxelsSingle[ix][iy][iz].depEnergyNorm=0.;
                voxelsSingle[ix][iy][iz].depEnergyNormError=0.;
                voxelsSingle[ix][iy][iz].expDose=0.;
                voxelsSingle[ix][iy][iz].halfSize.set(halfXVoxelDimensionX, halfXVoxelDimensionY, halfXVoxelDimensionZ);
                voxelsSingle[ix][iy][iz].pos.set(2.*(ix)*halfXVoxelDimensionX  -halfSize.getX()+halfXVoxelDimensionX + centre.getX(), 
        2.*(iy)*halfXVoxelDimensionY  -halfSize.getY()+halfXVoxelDimensionY + centre.getY(), 
        2.*(iz)*halfXVoxelDimensionZ  -halfSize.getZ()+halfXVoxelDimensionZ + centre.getZ());
                voxelsSingle[ix][iy][iz].nEvents=0;
            }
        }
    }
    nSingleTotalEvents=0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

save()

void CML2SDWithVoxels::save

(

)

Definition at line 222 of file ML2SDWithVoxels.cc.

{
std::cout<< "n. of events collected in the whole ROG phantom for all geometries: "<< nTotalEvents<< G4endl;
std::cout<< "n. of events collected in the whole ROG phantom for the current geometry: "<< nSingleTotalEvents<< G4endl;
    if (nTotalEvents>0)
    {saveData(fullOutFileData, voxelsSum);}
    if (nSingleTotalEvents>0)
    {saveData(fullOutFileDataSingle, voxelsSingle);}
}

Here is the call graph for this function:

Here is the caller graph for this function:

saveData()

void CML2SDWithVoxels::saveData ( G4String  Filename, Svoxel ***  voxels  )

private

Definition at line 231 of file ML2SDWithVoxels.cc.

{
    std::ofstream out;
    out.open(Filename, std::ios::out);
    out << "Sensitive Detector-Voxels. Total number of events, [mm]->centreX centreY centreZ HalfSizeX HalfSizeY HalfSizeZ minX maxX, minY maxY, minZ maxZ, Dx, Dy, Dz, nX, nY, nZ: \n";
    out <<nTotalEvents<<'\t';
    out <<centre.getX()/mm << '\t' << centre.getY()/mm<< '\t'<< centre.getZ()/mm<<'\t';
    out <<halfSize.getX()/mm << '\t' << halfSize.getY()/mm <<'\t'<< halfSize.getZ()/mm<<'\t';
    out <<(centre.getX()-halfSize.getX())/mm<<'\t'<<(centre.getX()+halfSize.getX())/mm<<'\t';
    out <<(centre.getY()-halfSize.getY())/mm<<'\t'<<(centre.getY()+halfSize.getY())/mm<<'\t';
    out <<(centre.getZ()-halfSize.getZ())/mm<<'\t'<<(centre.getZ()+halfSize.getZ())/mm<<'\t';
    out <<halfXVoxelDimensionX/mm<<'\t'<<halfXVoxelDimensionY/mm<<'\t'<<halfXVoxelDimensionZ/mm<<'\t';
    out <<NumberOfVoxelsAlongX <<'\t'<<NumberOfVoxelsAlongY <<'\t'<<NumberOfVoxelsAlongZ <<'\n';
    out << "Number of physical volumes: "<< volumeNameIdLink.size() << '\n';
    for (int i=0; i<(int)volumeNameIdLink.size(); i++)
    {
        out << volumeNameIdLink[i].volumeName  <<'\t'<< volumeNameIdLink[i].volumeId << G4endl;
    }


    out << "Phys Volume x [mm], y [mm], z [mm], ix, iy, iz, Dose [Gy], Dose2 [Gy^2], nEvents" << G4endl;

    for (int ix=0; ix< NumberOfVoxelsAlongX; ix++)
    {
        for (int iy=0; iy< NumberOfVoxelsAlongY; iy++)
        {
            for (int iz=0; iz< NumberOfVoxelsAlongZ; iz++)
            {
                if (voxels[ix][iy][iz].nEvents>0)
                {
                    out << voxels[ix][iy][iz].volumeId << '\t';
                    out << voxels[ix][iy][iz].pos.getX()/mm << '\t';
                    out << voxels[ix][iy][iz].pos.getY()/mm << '\t';
                    out << voxels[ix][iy][iz].pos.getZ()/mm << '\t';
                    out << ix << '\t';
                    out << iy << '\t';
                    out << iz << '\t';
                    out << voxels[ix][iy][iz].depEnergy/(joule/kg) << '\t';
                    out << voxels[ix][iy][iz].depEnergy2/((joule/kg)*(joule/kg)) << '\t';
                    out << voxels[ix][iy][iz].nEvents << G4endl;
                }
            }
        }
    }
    out.close();

}

Here is the call graph for this function:

Here is the caller graph for this function:

setActive()

void CML2SDWithVoxels::setActive ( G4bool  act )

inline

Definition at line 64 of file ML2SDWithVoxels.hh.

{bActive=act;}

Here is the call graph for this function:

setFullOutFileDataSingle()

void CML2SDWithVoxels::setFullOutFileDataSingle

(

G4String

val

)

Definition at line 278 of file ML2SDWithVoxels.cc.

{
    unsigned int ind = fullOutFileData.find(".txt");
    G4String onlyName=fullOutFileData.substr( 0, ind);
    if (val=="")
    {
        static unsigned int indGeom=0;
        char cT[5];
        sprintf(cT,"%d",indGeom);
        fullOutFileDataSingle=onlyName+"Single_"+G4String(cT)+".txt";
        indGeom++;
    }
    else
    {
        fullOutFileDataSingle=onlyName+val+".txt";
    }
}

Here is the caller graph for this function:

setRecycling()

void CML2SDWithVoxels::setRecycling ( int  recycling )

inline

Definition at line 66 of file ML2SDWithVoxels.hh.

{nRecycling=recycling;}

setVolumeNameIdLink()

void CML2SDWithVoxels::setVolumeNameIdLink ( std::vector< SvolumeNameId >  volNameIdLink )

inline

Definition at line 67 of file ML2SDWithVoxels.hh.

{volumeNameIdLink=volNameIdLink;}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

bActive

G4bool CML2SDWithVoxels::bActive

private

Definition at line 81 of file ML2SDWithVoxels.hh.

bSaveROG

G4bool CML2SDWithVoxels::bSaveROG

private

Definition at line 81 of file ML2SDWithVoxels.hh.

centre

G4ThreeVector CML2SDWithVoxels::centre

private

Definition at line 74 of file ML2SDWithVoxels.hh.

density

G4double CML2SDWithVoxels::density

private

Definition at line 82 of file ML2SDWithVoxels.hh.

fullOutFileData

G4String CML2SDWithVoxels::fullOutFileData

private

Definition at line 79 of file ML2SDWithVoxels.hh.

fullOutFileDataSingle

G4String CML2SDWithVoxels::fullOutFileDataSingle

private

Definition at line 79 of file ML2SDWithVoxels.hh.

halfSize

G4ThreeVector CML2SDWithVoxels::halfSize

private

Definition at line 74 of file ML2SDWithVoxels.hh.

halfXVoxelDimensionX

G4double CML2SDWithVoxels::halfXVoxelDimensionX

private

Definition at line 76 of file ML2SDWithVoxels.hh.

halfXVoxelDimensionY

G4double CML2SDWithVoxels::halfXVoxelDimensionY

private

Definition at line 76 of file ML2SDWithVoxels.hh.

halfXVoxelDimensionZ

G4double CML2SDWithVoxels::halfXVoxelDimensionZ

private

Definition at line 76 of file ML2SDWithVoxels.hh.

nParticle

G4int CML2SDWithVoxels::nParticle

private

Definition at line 80 of file ML2SDWithVoxels.hh.

nParticleValatile

G4int CML2SDWithVoxels::nParticleValatile

private

Definition at line 80 of file ML2SDWithVoxels.hh.

nRecycling

G4int CML2SDWithVoxels::nRecycling

private

Definition at line 80 of file ML2SDWithVoxels.hh.

nSingleTotalEvents

G4int CML2SDWithVoxels::nSingleTotalEvents

private

Definition at line 80 of file ML2SDWithVoxels.hh.

nTotalEvents

G4int CML2SDWithVoxels::nTotalEvents

private

Definition at line 80 of file ML2SDWithVoxels.hh.

NumberOfVoxelsAlongX

G4int CML2SDWithVoxels::NumberOfVoxelsAlongX

private

Definition at line 77 of file ML2SDWithVoxels.hh.

NumberOfVoxelsAlongY

G4int CML2SDWithVoxels::NumberOfVoxelsAlongY

private

Definition at line 77 of file ML2SDWithVoxels.hh.

NumberOfVoxelsAlongZ

G4int CML2SDWithVoxels::NumberOfVoxelsAlongZ

private

Definition at line 77 of file ML2SDWithVoxels.hh.

pos

G4ThreeVector CML2SDWithVoxels::pos

private

Definition at line 75 of file ML2SDWithVoxels.hh.

saving_in_ROG_Voxels_every_events

G4int CML2SDWithVoxels::saving_in_ROG_Voxels_every_events

private

Definition at line 80 of file ML2SDWithVoxels.hh.

volumeNameIdLink

std::vector<SvolumeNameId> CML2SDWithVoxels::volumeNameIdLink

private

Definition at line 83 of file ML2SDWithVoxels.hh.

voxelMass

G4double CML2SDWithVoxels::voxelMass

private

Definition at line 82 of file ML2SDWithVoxels.hh.

voxelsSingle

Svoxel *** CML2SDWithVoxels::voxelsSingle

private

Definition at line 78 of file ML2SDWithVoxels.hh.

voxelsSum

Svoxel*** CML2SDWithVoxels::voxelsSum

private

Definition at line 78 of file ML2SDWithVoxels.hh.

voxelVolume

G4double CML2SDWithVoxels::voxelVolume

private

Definition at line 82 of file ML2SDWithVoxels.hh.

---

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

• ML2SDWithVoxels.hh
• ML2SDWithVoxels.cc