Implements a G4VNestedParameterisation. More...

#include <Doxymodules_medical.h>

Inheritance diagram for DicomNestedPhantomParameterisation:

Public Types
typedef std::map< G4String, G4VisAttributes * >	ColourMap_t

Public Member Functions
	DicomNestedPhantomParameterisation (const G4ThreeVector &voxelSize, std::vector< G4Material * > &mat, G4int fnZ_=0, G4int fnY_=0, G4int fnX_=0, G4String colorFile=fDefaultColorFile)

	~DicomNestedPhantomParameterisation ()

G4Material *	ComputeMaterial (G4VPhysicalVolume currentVol, const G4int repNo, const G4VTouchable parentTouch)

G4int	GetNumberOfMaterials () const

G4Material *	GetMaterial (G4int idx) const

unsigned int	GetMaterialIndex (unsigned int copyNo) const

void	SetMaterialIndices (size_t *matInd)

void	SetNoVoxels (unsigned int nx, unsigned int ny, unsigned int nz)

void	ComputeTransformation (const G4int no, G4VPhysicalVolume *currentPV) const

void	ComputeDimensions (G4Box &, const G4int, const G4VPhysicalVolume *) const

const ColourMap_t &	GetColourMap () const

ColourMap_t &	GetColourMap ()

Static Public Attributes
static G4String	fDefaultColorFile

Private Member Functions
void	ComputeDimensions (G4Trd &, const G4int, const G4VPhysicalVolume *) const

void	ComputeDimensions (G4Trap &, const G4int, const G4VPhysicalVolume *) const

void	ComputeDimensions (G4Cons &, const G4int, const G4VPhysicalVolume *) const

void	ComputeDimensions (G4Sphere &, const G4int, const G4VPhysicalVolume *) const

void	ComputeDimensions (G4Ellipsoid &, const G4int, const G4VPhysicalVolume *) const

void	ComputeDimensions (G4Orb &, const G4int, const G4VPhysicalVolume *) const

void	ComputeDimensions (G4Torus &, const G4int, const G4VPhysicalVolume *) const

void	ComputeDimensions (G4Para &, const G4int, const G4VPhysicalVolume *) const

void	ComputeDimensions (G4Hype &, const G4int, const G4VPhysicalVolume *) const

void	ComputeDimensions (G4Tubs &, const G4int, const G4VPhysicalVolume *) const

void	ComputeDimensions (G4Polycone &, const G4int, const G4VPhysicalVolume *) const

void	ComputeDimensions (G4Polyhedra &, const G4int, const G4VPhysicalVolume *) const

void	ReadColourData (G4String)

Private Attributes
G4double	fdX

G4double	fdY

G4double	fdZ

G4int	fnX

G4int	fnY

G4int	fnZ

std::vector< G4Material * >	fMaterials

size_t *	fMaterialIndices

ColourMap_t	fColours

std::map< G4int, G4VisAttributes * >	fColours2

std::vector< G4double >	fpZ

Detailed Description

Implements a G4VNestedParameterisation.

Definition at line 24 of file Doxymodules_medical.h.

Member Typedef Documentation

◆ ColourMap_t

typedef std::map<G4String, G4VisAttributes*> DicomNestedPhantomParameterisation::ColourMap_t

Definition at line 68 of file DicomNestedPhantomParameterisation.hh.

Constructor & Destructor Documentation

◆ DicomNestedPhantomParameterisation()

DicomNestedPhantomParameterisation::DicomNestedPhantomParameterisation	(	const G4ThreeVector &	voxelSize,
		std::vector< G4Material * > &	mat,
		G4int	fnZ_ = `0`,
		G4int	fnY_ = `0`,
		G4int	fnX_ = `0`,
		G4String	colorFile = `fDefaultColorFile`
	)

Definition at line 50 of file DicomNestedPhantomParameterisation.cc.

:
  //G4VNestedParameterisation(),
  fdX(voxelSize.x()), fdY(voxelSize.y()), fdZ(voxelSize.z()),
  fnX(fnX_), fnY(fnY_), fnZ(fnZ_),
  fMaterials(mat),
  fMaterialIndices(0)
{
    ReadColourData(colorfile);
}

◆ ~DicomNestedPhantomParameterisation()

DicomNestedPhantomParameterisation::~DicomNestedPhantomParameterisation ( )

Definition at line 66 of file DicomNestedPhantomParameterisation.cc.

67{

68}

Member Function Documentation

◆ ComputeMaterial()

G4Material * DicomNestedPhantomParameterisation::ComputeMaterial	(	G4VPhysicalVolume *	currentVol,
		const G4int	repNo,
		const G4VTouchable *	parentTouch
	)

Definition at line 107 of file DicomNestedPhantomParameterisation.cc.

{
 
    // protection for initialization and vis at idle state
    //
    if(parentTouch == nullptr)
        return fMaterials[0];
 
    // Copy number of voxels.
    // Copy number of X and Y are obtained from replication number.
    // Copy nymber of Z is the copy number of current voxel.
    G4int ix = parentTouch->GetReplicaNumber(0);
    G4int iy = parentTouch->GetReplicaNumber(1);
 
    G4int copyID = ix + fnX*iy + fnX*fnY*iz;
 
    std::size_t matIndex = GetMaterialIndex(copyID);
    G4Material* mate = fMaterials[matIndex];
 
    if(G4VVisManager::GetConcreteInstance() && physVol)
    {
        G4String mateName = fMaterials.at(matIndex)->GetName();
        std::string::size_type iuu = mateName.find("__");
        if( iuu != std::string::npos )
            mateName = mateName.substr( 0, iuu );
 
        if(0 < fColours.count(mateName))
            physVol->GetLogicalVolume()->SetVisAttributes(
                fColours.find(mateName)->second);
        else
        {
            bool found = false;
            for(const auto& itr : fColours)
            {
                G4String mat_color = itr.first;
                auto len = mat_color.length();
                if(mateName.find(mat_color) == 0 &&
                   mateName.length() > len && mateName[len] == '_')
                {
                    physVol->GetLogicalVolume()->SetVisAttributes(
                        fColours.find(mat_color)->second);
                    found = true;
                }
                if(found)
                    break;
            }
            if(!found)
            {
 
                static uintmax_t n = 0;
                if(n++ < 100)
                    G4cout << "Unknown material name " << mateName
                           << " for index " << matIndex << G4endl;
                if(fColours2.find(matIndex) != fColours2.end())
                    physVol->GetLogicalVolume()->SetVisAttributes(
                        fColours2.find(matIndex)->second);
                else
                    physVol->GetLogicalVolume()->SetVisAttributes(
                        fColours.begin()->second);
            }
        }
        physVol -> GetLogicalVolume()->SetMaterial(mate);
    }
 
    return mate;
}

◆ GetNumberOfMaterials()

G4int DicomNestedPhantomParameterisation::GetNumberOfMaterials ( ) const

Definition at line 189 of file DicomNestedPhantomParameterisation.cc.

{
    return G4int(fMaterials.size());
}

◆ GetMaterial()

G4Material * DicomNestedPhantomParameterisation::GetMaterial ( G4int idx ) const

Definition at line 199 of file DicomNestedPhantomParameterisation.cc.

{
    return fMaterials[i];
}

◆ GetMaterialIndex()

unsigned int DicomNestedPhantomParameterisation::GetMaterialIndex ( unsigned int copyNo ) const

Definition at line 177 of file DicomNestedPhantomParameterisation.cc.

{
    //return *(fMaterialIndices+copyNo);
    return unsigned(fMaterialIndices[copyNo]);
}

◆ SetMaterialIndices()

void DicomNestedPhantomParameterisation::SetMaterialIndices ( size_t * matInd )

inline

Definition at line 94 of file DicomNestedPhantomParameterisation.hh.

94{ fMaterialIndices = matInd; }

◆ SetNoVoxels()

void DicomNestedPhantomParameterisation::SetNoVoxels	(	unsigned int	nx,
		unsigned int	ny,
		unsigned int	nz
	)

Definition at line 98 of file DicomNestedPhantomParameterisation.cc.

{
  fnX = nx;
  fnY = ny;
  fnZ = nz;
}

◆ ComputeTransformation()

void DicomNestedPhantomParameterisation::ComputeTransformation	(	const G4int	no,
		G4VPhysicalVolume *	currentPV
	)		const

Definition at line 207 of file DicomNestedPhantomParameterisation.cc.

{
    // Position of voxels.
    // x and y positions are already defined in DetectorConstruction by using
    // replicated volume. Here only we need to define is z positions of voxels.
    physVol->SetTranslation(G4ThreeVector(0.,0.,(2.*static_cast<double>(copyNo)
                                                +1.)*fdZ - fdZ*fnZ));
}

◆ ComputeDimensions() [1/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Box &	box,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

Definition at line 221 of file DicomNestedPhantomParameterisation.cc.

{
    box.SetXHalfLength(fdX);
    box.SetYHalfLength(fdY);
    box.SetZHalfLength(fdZ);
}

◆ GetColourMap() [1/2]

const ColourMap_t & DicomNestedPhantomParameterisation::GetColourMap ( ) const

inline

Definition at line 105 of file DicomNestedPhantomParameterisation.hh.

105{ return fColours; }

◆ GetColourMap() [2/2]

ColourMap_t & DicomNestedPhantomParameterisation::GetColourMap ( )

inline

Definition at line 106 of file DicomNestedPhantomParameterisation.hh.

106{ return fColours; }

◆ ComputeDimensions() [2/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Trd &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 110 of file DicomNestedPhantomParameterisation.hh.

111 {}

◆ ComputeDimensions() [3/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Trap &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 112 of file DicomNestedPhantomParameterisation.hh.

113 {}

◆ ComputeDimensions() [4/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Cons &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 114 of file DicomNestedPhantomParameterisation.hh.

115 {}

◆ ComputeDimensions() [5/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Sphere &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 116 of file DicomNestedPhantomParameterisation.hh.

117 {}

◆ ComputeDimensions() [6/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Ellipsoid &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 118 of file DicomNestedPhantomParameterisation.hh.

119 {}

◆ ComputeDimensions() [7/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Orb &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 120 of file DicomNestedPhantomParameterisation.hh.

121 {}

◆ ComputeDimensions() [8/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Torus &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 122 of file DicomNestedPhantomParameterisation.hh.

123 {}

◆ ComputeDimensions() [9/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Para &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 124 of file DicomNestedPhantomParameterisation.hh.

125 {}

◆ ComputeDimensions() [10/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Hype &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 126 of file DicomNestedPhantomParameterisation.hh.

127 {}

◆ ComputeDimensions() [11/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Tubs &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 128 of file DicomNestedPhantomParameterisation.hh.

129 {}

◆ ComputeDimensions() [12/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Polycone &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 130 of file DicomNestedPhantomParameterisation.hh.

131 {}

◆ ComputeDimensions() [13/13]

void DicomNestedPhantomParameterisation::ComputeDimensions	(	G4Polyhedra &	,
		const G4int	,
		const G4VPhysicalVolume *
	)		const

inlineprivate

Definition at line 132 of file DicomNestedPhantomParameterisation.hh.

133 {}

◆ ReadColourData()

void DicomNestedPhantomParameterisation::ReadColourData ( G4String colourFile )

private

Definition at line 71 of file DicomNestedPhantomParameterisation.cc.

{
    //----- Add a G4VisAttributes for materials not defined in file;
    G4VisAttributes* blankAtt = new G4VisAttributes;
    blankAtt->SetVisibility( FALSE );
    fColours["Default"] = blankAtt;
 
    std::ifstream fin(colourFile.c_str());
    G4int nMate;
    G4String mateName;
    G4double cred, cgreen, cblue, copacity;
    fin >> nMate;
    for( G4int ii = 0; ii < nMate; ii++ )
    {
        fin >> mateName;
        if(fin.eof())
            break;
        fin >> cred >> cgreen >> cblue >> copacity;
        G4Colour colour( cred, cgreen, cblue, copacity );
        G4VisAttributes* visAtt = new G4VisAttributes( colour );
        visAtt->SetVisibility(true);
        fColours[mateName] = visAtt;
        fColours2[ii] = new G4VisAttributes(*visAtt);
    }
}