DicomDetectorConstruction Class Reference

#include <DicomDetectorConstruction.hh>

Inheritance diagram for DicomDetectorConstruction:

Collaboration diagram for DicomDetectorConstruction:

Public Member Functions
	DicomDetectorConstruction ()
	~DicomDetectorConstruction ()
virtual G4VPhysicalVolume *	Construct ()
Public Member Functions inherited from G4VUserDetectorConstruction
	G4VUserDetectorConstruction ()
virtual	~G4VUserDetectorConstruction ()
virtual void	CloneSD ()
virtual void	CloneF ()
void	RegisterParallelWorld (G4VUserParallelWorld *)
G4int	ConstructParallelGeometries ()
void	ConstructParallelSD ()
G4int	GetNumberOfParallelWorld () const
G4VUserParallelWorld *	GetParallelWorld (G4int i) const

Protected Member Functions
void	InitialisationOfMaterials ()
void	ReadPhantomData ()
void	ReadPhantomDataFile (const G4String &fname)
void	MergeZSliceHeaders ()
G4Material *	BuildMaterialWithChangingDensity (const G4Material *origMate, float density, G4String newMateName)
void	ConstructPhantomContainer ()
virtual void	ConstructPhantom ()=0
void	SetScorer (G4LogicalVolume *voxel_logic)
virtual void	ConstructSDandField ()
Protected Member Functions inherited from G4VUserDetectorConstruction
void	SetSensitiveDetector (const G4String &logVolName, G4VSensitiveDetector *aSD, G4bool multi=false)
void	SetSensitiveDetector (G4LogicalVolume *logVol, G4VSensitiveDetector *aSD)

Protected Attributes
G4Material *	fAir
G4Box *	fWorld_solid
G4LogicalVolume *	fWorld_logic
G4VPhysicalVolume *	fWorld_phys
G4Box *	fContainer_solid
G4LogicalVolume *	fContainer_logic
G4VPhysicalVolume *	fContainer_phys
G4int	fNoFiles
std::vector< G4Material * >	fOriginalMaterials
std::vector< G4Material * >	fMaterials
size_t *	fMateIDs
std::map< G4int, G4double >	fDensityDiffs
std::vector< DicomPhantomZSliceHeader * >	fZSliceHeaders
DicomPhantomZSliceHeader *	fZSliceHeaderMerged
G4int	fNVoxelX
G4int	fNVoxelY
G4int	fNVoxelZ
G4double	fVoxelHalfDimX
G4double	fVoxelHalfDimY
G4double	fVoxelHalfDimZ
DicomPhantomZSliceMerged *	fMergedSlices
std::set< G4LogicalVolume * >	fScorers
G4bool	fConstructed

Detailed Description

Dicom detector construction

 - Start the building of the geometry
 - Initialisation of materials
 - Creation of the world
 - Reading of the DICOM data

History: 30.11.07 First version

Author

P. Arce

Definition at line 58 of file DicomDetectorConstruction.hh.

Constructor & Destructor Documentation

DicomDetectorConstruction()

DicomDetectorConstruction::DicomDetectorConstruction

(

)

Definition at line 53 of file DicomDetectorConstruction.cc.

 : G4VUserDetectorConstruction(),
   fAir(0),

   fWorld_solid(0),
   fWorld_logic(0),
   fWorld_phys(0),

   fContainer_solid(0),
   fContainer_logic(0),
   fContainer_phys(0),

   fNoFiles(0),
   fMateIDs(0),

   fZSliceHeaderMerged(0),

   fNVoxelX(0),
   fNVoxelY(0),
   fNVoxelZ(0),
   fVoxelHalfDimX(0),
   fVoxelHalfDimY(0),
   fVoxelHalfDimZ(0),

   fConstructed(false)
{

}

~DicomDetectorConstruction()

DicomDetectorConstruction::~DicomDetectorConstruction

(

)

Definition at line 83 of file DicomDetectorConstruction.cc.

{
}

Member Function Documentation

BuildMaterialWithChangingDensity()

G4Material * DicomDetectorConstruction::BuildMaterialWithChangingDensity ( const G4Material *  origMate, float  density, G4String  newMateName  )

protected

Definition at line 467 of file DicomDetectorConstruction.cc.

{
  //----- Copy original material, but with new density
  G4int nelem = origMate->GetNumberOfElements();
  G4Material* mate = new G4Material( newMateName, density*g/cm3, nelem,
                                     kStateUndefined, STP_Temperature );
  
  for( G4int ii = 0; ii < nelem; ii++ ){
    G4double frac = origMate->GetFractionVector()[ii];
    G4Element* elem = const_cast<G4Element*>(origMate->GetElement(ii));
    mate->AddElement( elem, frac );
  }
  
  return mate;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Construct()

G4VPhysicalVolume * DicomDetectorConstruction::Construct ( void  )

virtual

Implements G4VUserDetectorConstruction.

Reimplemented in DicomPartialDetectorConstruction.

Definition at line 88 of file DicomDetectorConstruction.cc.

{
  if(!fConstructed || fWorld_phys == 0) {
    fConstructed = true;
    InitialisationOfMaterials();

    //----- Build world
    G4double worldXDimension = 1.*m;
    G4double worldYDimension = 1.*m;
    G4double worldZDimension = 1.*m;

    fWorld_solid = new G4Box( "WorldSolid",
                             worldXDimension,
                             worldYDimension,
                             worldZDimension );

    fWorld_logic = new G4LogicalVolume( fWorld_solid,
                                       fAir,
                                       "WorldLogical",
                                       0, 0, 0 );

    fWorld_phys = new G4PVPlacement( 0,
                                    G4ThreeVector(0,0,0),
                                    "World",
                                    fWorld_logic,
                                    0,
                                    false,
                                    0 );

    ReadPhantomData();

    ConstructPhantomContainer();
    ConstructPhantom();
  }
    return fWorld_phys;
}

Here is the call graph for this function:

ConstructPhantom()

virtual void DicomDetectorConstruction::ConstructPhantom ( )

protectedpure virtual

Implemented in DicomPartialDetectorConstruction, DicomRegularDetectorConstruction, and DicomNestedParamDetectorConstruction.

Here is the caller graph for this function:

ConstructPhantomContainer()

void DicomDetectorConstruction::ConstructPhantomContainer ( )

protected

Definition at line 485 of file DicomDetectorConstruction.cc.

{
  //---- Extract number of voxels and voxel dimensions
  fNVoxelX = fZSliceHeaderMerged->GetNoVoxelX();
  fNVoxelY = fZSliceHeaderMerged->GetNoVoxelY();
  fNVoxelZ = fZSliceHeaderMerged->GetNoVoxelZ();
  
  fVoxelHalfDimX = fZSliceHeaderMerged->GetVoxelHalfX();
  fVoxelHalfDimY = fZSliceHeaderMerged->GetVoxelHalfY();
  fVoxelHalfDimZ = fZSliceHeaderMerged->GetVoxelHalfZ();
#ifdef G4VERBOSE
  G4cout << " fNVoxelX " << fNVoxelX << " fVoxelHalfDimX " << fVoxelHalfDimX 
         <<G4endl;
  G4cout << " fNVoxelY " << fNVoxelY << " fVoxelHalfDimY " << fVoxelHalfDimY 
         <<G4endl;
  G4cout << " fNVoxelZ " << fNVoxelZ << " fVoxelHalfDimZ " << fVoxelHalfDimZ 
         <<G4endl;
  G4cout << " totalPixels " << fNVoxelX*fNVoxelY*fNVoxelZ <<  G4endl;
#endif
  
  //----- Define the volume that contains all the voxels
  fContainer_solid = new G4Box("phantomContainer",fNVoxelX*fVoxelHalfDimX,
                               fNVoxelY*fVoxelHalfDimY,
                               fNVoxelZ*fVoxelHalfDimZ);
  fContainer_logic =
    new G4LogicalVolume( fContainer_solid,
   //the material is not important, it will be fully filled by the voxels
                         fMaterials[0],
                         "phantomContainer",
                         0, 0, 0 );
  //--- Place it on the world
  G4double fOffsetX = (fZSliceHeaderMerged->GetMaxX() + 
                       fZSliceHeaderMerged->GetMinX() ) /2.;
  G4double fOffsetY = (fZSliceHeaderMerged->GetMaxY() + 
                       fZSliceHeaderMerged->GetMinY() ) /2.;
  G4double fOffsetZ = (fZSliceHeaderMerged->GetMaxZ() + 
                       fZSliceHeaderMerged->GetMinZ() ) /2.;
  G4ThreeVector posCentreVoxels(fOffsetX,fOffsetY,fOffsetZ);
#ifdef G4VERBOSE
  G4cout << " placing voxel container volume at " << posCentreVoxels << G4endl;
#endif
  fContainer_phys =
    new G4PVPlacement(0,  // rotation
                      posCentreVoxels,
                      fContainer_logic,     // The logic volume
                      "phantomContainer",  // Name
                      fWorld_logic,  // Mother
                      false,           // No op. bool.
                      1);              // Copy number
  
  //fContainer_logic->SetVisAttributes(new G4VisAttributes(G4Colour(1.,0.,0.)));
}

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructSDandField()

void DicomDetectorConstruction::ConstructSDandField ( )

protectedvirtual

Reimplemented from G4VUserDetectorConstruction.

Definition at line 556 of file DicomDetectorConstruction.cc.

{
  
  G4cout << "\n\n\n\n\t CONSTRUCT SD AND FIELD \n\n\n" << G4endl;
  
  //G4SDManager* SDman = G4SDManager::GetSDMpointer();
  
  //SDman->SetVerboseLevel(1);
  
  //
  // Sensitive Detector Name
  G4String concreteSDname = "phantomSD";
  std::vector<G4String> scorer_names;
  scorer_names.push_back(concreteSDname);
  //------------------------
  // MultiFunctionalDetector
  //------------------------
  //
  // Define MultiFunctionalDetector with name.
  // declare MFDet as a MultiFunctionalDetector scorer
  G4MultiFunctionalDetector* MFDet = 
    new G4MultiFunctionalDetector(concreteSDname);
  //SDman->AddNewDetector( MFDet );                 // Register SD to SDManager
  //G4VPrimitiveScorer* dosedep = new G4PSDoseDeposit("DoseDeposit");
  G4VPrimitiveScorer* dosedep = 
    new G4PSDoseDeposit3D("DoseDeposit", fNVoxelX, fNVoxelY, fNVoxelZ);
  MFDet->RegisterPrimitive(dosedep);
  
  for(std::set<G4LogicalVolume*>::iterator ite = fScorers.begin(); 
      ite != fScorers.end(); ++ite) {
    SetSensitiveDetector(*ite, MFDet);
  }
  
  /*if(DicomRunAction::Instance()->GetDicomRun()) {
      DicomRunAction::Instance()->GetDicomRun()->ConstructMFD(scorer_names);
      }*/
 
}

Here is the call graph for this function:

InitialisationOfMaterials()

void DicomDetectorConstruction::InitialisationOfMaterials ( )

protected

Definition at line 126 of file DicomDetectorConstruction.cc.

{
    // Creating elements :
    G4double z, a, density;
    G4String name, symbol;

    G4Element* elC = new G4Element( name = "Carbon",
                                   symbol = "C",
                                   z = 6.0, a = 12.011 * g/mole );
    G4Element* elH = new G4Element( name = "Hydrogen",
                                   symbol = "H",
                                   z = 1.0, a = 1.008  * g/mole );
    G4Element* elN = new G4Element( name = "Nitrogen",
                                   symbol = "N",
                                   z = 7.0, a = 14.007 * g/mole );
    G4Element* elO = new G4Element( name = "Oxygen",
                                   symbol = "O",
                                   z = 8.0, a = 16.00  * g/mole );
    G4Element* elNa = new G4Element( name = "Sodium",
                                    symbol = "Na",
                                    z= 11.0, a = 22.98977* g/mole );
    G4Element* elS = new G4Element( name = "Sulfur",
                                   symbol = "S",
                                   z = 16.0,a = 32.065* g/mole );
    G4Element* elCl = new G4Element( name = "Chlorine",
                                    symbol = "P",
                                    z = 17.0, a = 35.453* g/mole );
    G4Element* elK = new G4Element( name = "Potassium",
                                   symbol = "P",
                                   z = 19.0, a = 30.0983* g/mole );
    G4Element* elP = new G4Element( name = "Phosphorus",
                                   symbol = "P",
                                   z = 30.0, a = 30.973976* g/mole );
    G4Element* elFe = new G4Element( name = "Iron",
                                    symbol = "Fe",
                                    z = 26, a = 56.845* g/mole );
    G4Element* elMg = new G4Element( name = "Magnesium",
                                    symbol = "Mg",
                                    z = 12.0, a = 24.3050* g/mole );
    G4Element* elCa = new G4Element( name="Calcium",
                                    symbol = "Ca",
                                    z = 20.0, a = 40.078* g/mole );

    // Creating Materials :
    G4int numberofElements;

    // Air
    fAir = new G4Material( "Air",
                          1.290*mg/cm3,
                          numberofElements = 2 );
    fAir->AddElement(elN, 0.7);
    fAir->AddElement(elO, 0.3);

    //  Lung Inhale
    G4Material* lunginhale = new G4Material( "LungInhale",
                                            density = 0.217*g/cm3,
                                            numberofElements = 9);
    lunginhale->AddElement(elH,0.103);
    lunginhale->AddElement(elC,0.105);
    lunginhale->AddElement(elN,0.031);
    lunginhale->AddElement(elO,0.749);
    lunginhale->AddElement(elNa,0.002);
    lunginhale->AddElement(elP,0.002);
    lunginhale->AddElement(elS,0.003);
    lunginhale->AddElement(elCl,0.002);
    lunginhale->AddElement(elK,0.003);

    // Lung exhale
    G4Material* lungexhale = new G4Material( "LungExhale",
                                            density = 0.508*g/cm3,
                                            numberofElements = 9 );
    lungexhale->AddElement(elH,0.103);
    lungexhale->AddElement(elC,0.105);
    lungexhale->AddElement(elN,0.031);
    lungexhale->AddElement(elO,0.749);
    lungexhale->AddElement(elNa,0.002);
    lungexhale->AddElement(elP,0.002);
    lungexhale->AddElement(elS,0.003);
    lungexhale->AddElement(elCl,0.002);
    lungexhale->AddElement(elK,0.003);

    // Adipose tissue
    G4Material* adiposeTissue = new G4Material( "AdiposeTissue",
                                               density = 0.967*g/cm3,
                                               numberofElements = 7);
    adiposeTissue->AddElement(elH,0.114);
    adiposeTissue->AddElement(elC,0.598);
    adiposeTissue->AddElement(elN,0.007);
    adiposeTissue->AddElement(elO,0.278);
    adiposeTissue->AddElement(elNa,0.001);
    adiposeTissue->AddElement(elS,0.001);
    adiposeTissue->AddElement(elCl,0.001);

    // Breast
    G4Material* breast = new G4Material( "Breast",
                                        density = 0.990*g/cm3,
                                        numberofElements = 8 );
    breast->AddElement(elH,0.109);
    breast->AddElement(elC,0.506);
    breast->AddElement(elN,0.023);
    breast->AddElement(elO,0.358);
    breast->AddElement(elNa,0.001);
    breast->AddElement(elP,0.001);
    breast->AddElement(elS,0.001);
    breast->AddElement(elCl,0.001);

    // Water
    G4Material* water = new G4Material( "Water",
                                       density = 1.0*g/cm3,
                                       numberofElements = 2 );
    water->AddElement(elH,0.112);
    water->AddElement(elO,0.888);

    // Muscle
    G4Material* muscle = new G4Material( "Muscle",
                                        density = 1.061*g/cm3,
                                        numberofElements = 9 );
    muscle->AddElement(elH,0.102);
    muscle->AddElement(elC,0.143);
    muscle->AddElement(elN,0.034);
    muscle->AddElement(elO,0.710);
    muscle->AddElement(elNa,0.001);
    muscle->AddElement(elP,0.002);
    muscle->AddElement(elS,0.003);
    muscle->AddElement(elCl,0.001);
    muscle->AddElement(elK,0.004);

    // Liver
    G4Material* liver = new G4Material( "Liver",
                                       density = 1.071*g/cm3,
                                       numberofElements = 9);
    liver->AddElement(elH,0.102);
    liver->AddElement(elC,0.139);
    liver->AddElement(elN,0.030);
    liver->AddElement(elO,0.716);
    liver->AddElement(elNa,0.002);
    liver->AddElement(elP,0.003);
    liver->AddElement(elS,0.003);
    liver->AddElement(elCl,0.002);
    liver->AddElement(elK,0.003);

    // Trabecular Bone
    G4Material* trabecularBone = new G4Material( "TrabecularBone",
                                                density = 1.159*g/cm3,
                                                numberofElements = 12 );
    trabecularBone->AddElement(elH,0.085);
    trabecularBone->AddElement(elC,0.404);
    trabecularBone->AddElement(elN,0.058);
    trabecularBone->AddElement(elO,0.367);
    trabecularBone->AddElement(elNa,0.001);
    trabecularBone->AddElement(elMg,0.001);
    trabecularBone->AddElement(elP,0.034);
    trabecularBone->AddElement(elS,0.002);
    trabecularBone->AddElement(elCl,0.002);
    trabecularBone->AddElement(elK,0.001);
    trabecularBone->AddElement(elCa,0.044);
    trabecularBone->AddElement(elFe,0.001);

    // Dense Bone
    G4Material* denseBone = new G4Material( "DenseBone",
                                           density = 1.575*g/cm3,
                                           numberofElements = 11 );
    denseBone->AddElement(elH,0.056);
    denseBone->AddElement(elC,0.235);
    denseBone->AddElement(elN,0.050);
    denseBone->AddElement(elO,0.434);
    denseBone->AddElement(elNa,0.001);
    denseBone->AddElement(elMg,0.001);
    denseBone->AddElement(elP,0.072);
    denseBone->AddElement(elS,0.003);
    denseBone->AddElement(elCl,0.001);
    denseBone->AddElement(elK,0.001);
    denseBone->AddElement(elCa,0.146);

    //----- Put the materials in a vector
    fOriginalMaterials.push_back(fAir); // rho = 0.00129
    fOriginalMaterials.push_back(lunginhale); // rho = 0.217
    fOriginalMaterials.push_back(lungexhale); // rho = 0.508
    fOriginalMaterials.push_back(adiposeTissue); // rho = 0.967
    fOriginalMaterials.push_back(breast ); // rho = 0.990
    fOriginalMaterials.push_back(water); // rho = 1.018
    fOriginalMaterials.push_back(muscle); // rho = 1.061
    fOriginalMaterials.push_back(liver); // rho = 1.071
    fOriginalMaterials.push_back(trabecularBone); // rho = 1.159
    fOriginalMaterials.push_back(denseBone); // rho = 1.575

}

Here is the call graph for this function:

Here is the caller graph for this function:

MergeZSliceHeaders()

void DicomDetectorConstruction::MergeZSliceHeaders ( )

protected

Definition at line 456 of file DicomDetectorConstruction.cc.

{
  //----- Images must have the same dimension ...
  fZSliceHeaderMerged = new DicomPhantomZSliceHeader( *fZSliceHeaders[0] );
  for( unsigned int ii = 1; ii < fZSliceHeaders.size(); ii++ ) {
    *fZSliceHeaderMerged += *fZSliceHeaders[ii];
  };
  
}

Here is the caller graph for this function:

ReadPhantomData()

void DicomDetectorConstruction::ReadPhantomData ( )

protected

Definition at line 315 of file DicomDetectorConstruction.cc.

{

    G4String dataFile = "Data.dat";
    std::ifstream finDF(dataFile.c_str());
    G4String fname;
    if(finDF.good() != 1 ) {
        G4String descript = "Problem reading data file: "+dataFile;
        G4Exception(" DicomDetectorConstruction::ReadPhantomData",
                    "",
                    FatalException,
                    descript);
    }

    G4int compression;
    finDF >> compression; // not used here

    finDF >> fNoFiles;
    for(G4int i = 0; i < fNoFiles; i++ ) {
        finDF >> fname;
        //--- Read one data file
        fname += ".g4dcm";
        ReadPhantomDataFile(fname);
    }

    //----- Merge data headers
    MergeZSliceHeaders();

    finDF.close();

}

Here is the call graph for this function:

Here is the caller graph for this function:

ReadPhantomDataFile()

void DicomDetectorConstruction::ReadPhantomDataFile ( const G4String &  fname )

protected

Definition at line 348 of file DicomDetectorConstruction.cc.

{
#ifdef G4VERBOSE
  G4cout << " DicomDetectorConstruction::ReadPhantomDataFile opening file " 
         << fname << G4endl;
#endif
  std::ifstream fin(fname.c_str(), std::ios_base::in);
  if( !fin.is_open() ) {
    G4Exception("DicomDetectorConstruction::ReadPhantomDataFile",
                "",
                FatalErrorInArgument,
                G4String("File not found " + fname ).c_str());
  }
  //----- Define density differences (maximum density difference to create
  // a new material)
  char* part = getenv( "DICOM_CHANGE_MATERIAL_DENSITY" );
  G4double densityDiff = -1.;
  if( part ) densityDiff = G4UIcommand::ConvertToDouble(part);
  if( densityDiff != -1. ) {
    for( unsigned int ii = 0; ii < fOriginalMaterials.size(); ii++ ){
      fDensityDiffs[ii] = densityDiff; //currently all materials with 
      // same difference
    }
  }else {
    if( fMaterials.size() == 0 ) { // do it only for first slice
      for( unsigned int ii = 0; ii < fOriginalMaterials.size(); ii++ ){
        fMaterials.push_back( fOriginalMaterials[ii] );
      }
    }
  }
  
  //----- Read data header
  DicomPhantomZSliceHeader* sliceHeader = new DicomPhantomZSliceHeader( fin );
  fZSliceHeaders.push_back( sliceHeader );
  
  //----- Read material indices
  G4int nVoxels = sliceHeader->GetNoVoxels();
  
  //--- If first slice, initiliaze fMateIDs
  if( fZSliceHeaders.size() == 1 ) {
    //fMateIDs = new unsigned int[fNoFiles*nVoxels];
    fMateIDs = new size_t[fNoFiles*nVoxels];
    
  }
  
  unsigned int mateID;
  // number of voxels from previously read slices
  G4int voxelCopyNo = (fZSliceHeaders.size()-1)*nVoxels;
  for( G4int ii = 0; ii < nVoxels; ii++, voxelCopyNo++ ){
    fin >> mateID;
    fMateIDs[voxelCopyNo] = mateID;
  }
  
  //----- Read material densities and build new materials if two voxels have
  //  same material but its density is in a different density interval 
  // (size of density intervals defined by densityDiff)
  G4double density;
  // number of voxels from previously read slices
  voxelCopyNo = (fZSliceHeaders.size()-1)*nVoxels;
  for( G4int ii = 0; ii < nVoxels; ii++, voxelCopyNo++ ){
    fin >> density;
    
    //-- Get material from list of original materials
    mateID = fMateIDs[voxelCopyNo];
    G4Material* mateOrig  = fOriginalMaterials[mateID];
    
    //-- Get density bin: middle point of the bin in which the current
    // density is included
    G4String newMateName = mateOrig->GetName();
    float densityBin = 0.;
    if( densityDiff != -1.) {
      densityBin = fDensityDiffs[mateID] * 
        (G4int(density/fDensityDiffs[mateID])+0.5);
      //-- Build the new material name
      newMateName += G4UIcommand::ConvertToString(densityBin);
    }
    
    //-- Look if a material with this name is already created
    //  (because a previous voxel was already in this density bin)
    unsigned int im;
    for( im = 0; im < fMaterials.size(); im++ ){
      if( fMaterials[im]->GetName() == newMateName ) {
        break;
      }
    }
    //-- If material is already created use index of this material
    if( im != fMaterials.size() ) {
      fMateIDs[voxelCopyNo] = im;
      //-- else, create the material
    } else {
      if( densityDiff != -1.) {
        fMaterials.push_back( BuildMaterialWithChangingDensity( mateOrig,
                                                  densityBin, newMateName ) );
        fMateIDs[voxelCopyNo] = fMaterials.size()-1;
      } else {
        G4cerr << " im " << im << " < " << fMaterials.size() << " name "
               << newMateName << G4endl;
        G4Exception("DicomDetectorConstruction::ReadPhantomDataFile",
                    "",
                    FatalErrorInArgument,
                    "Wrong index in material"); //it should never reach here
      }
    }
  }
  
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetScorer()

void DicomDetectorConstruction::SetScorer ( G4LogicalVolume *  voxel_logic )

protected

Definition at line 544 of file DicomDetectorConstruction.cc.

{
  
  G4cout << "\n\n\n\n\t SET SCORER : " << voxel_logic->GetName() 
         << " \n\n\n" << G4endl;
  
  fScorers.insert(voxel_logic);
  
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

fAir

G4Material* DicomDetectorConstruction::fAir

protected

Definition at line 96 of file DicomDetectorConstruction.hh.

fConstructed

G4bool DicomDetectorConstruction::fConstructed

protected

Definition at line 130 of file DicomDetectorConstruction.hh.

fContainer_logic

G4LogicalVolume* DicomDetectorConstruction::fContainer_logic

protected

Definition at line 104 of file DicomDetectorConstruction.hh.

fContainer_phys

G4VPhysicalVolume* DicomDetectorConstruction::fContainer_phys

protected

Definition at line 105 of file DicomDetectorConstruction.hh.

fContainer_solid

G4Box* DicomDetectorConstruction::fContainer_solid

protected

Definition at line 103 of file DicomDetectorConstruction.hh.

fDensityDiffs

std::map<G4int,G4double> DicomDetectorConstruction::fDensityDiffs

protected

Definition at line 115 of file DicomDetectorConstruction.hh.

fMateIDs

size_t* DicomDetectorConstruction::fMateIDs

protected

Definition at line 112 of file DicomDetectorConstruction.hh.

fMaterials

std::vector<G4Material*> DicomDetectorConstruction::fMaterials

protected

Definition at line 109 of file DicomDetectorConstruction.hh.

fMergedSlices

DicomPhantomZSliceMerged* DicomDetectorConstruction::fMergedSlices

protected

Definition at line 126 of file DicomDetectorConstruction.hh.

fNoFiles

G4int DicomDetectorConstruction::fNoFiles

protected

Definition at line 107 of file DicomDetectorConstruction.hh.

fNVoxelX

G4int DicomDetectorConstruction::fNVoxelX

protected

Definition at line 123 of file DicomDetectorConstruction.hh.

fNVoxelY

G4int DicomDetectorConstruction::fNVoxelY

protected

Definition at line 123 of file DicomDetectorConstruction.hh.

fNVoxelZ

G4int DicomDetectorConstruction::fNVoxelZ

protected

Definition at line 123 of file DicomDetectorConstruction.hh.

fOriginalMaterials

std::vector<G4Material*> DicomDetectorConstruction::fOriginalMaterials

protected

Definition at line 108 of file DicomDetectorConstruction.hh.

fScorers

std::set<G4LogicalVolume*> DicomDetectorConstruction::fScorers

protected

Definition at line 128 of file DicomDetectorConstruction.hh.

fVoxelHalfDimX

G4double DicomDetectorConstruction::fVoxelHalfDimX

protected

Definition at line 124 of file DicomDetectorConstruction.hh.

fVoxelHalfDimY

G4double DicomDetectorConstruction::fVoxelHalfDimY

protected

Definition at line 124 of file DicomDetectorConstruction.hh.

fVoxelHalfDimZ

G4double DicomDetectorConstruction::fVoxelHalfDimZ

protected

Definition at line 124 of file DicomDetectorConstruction.hh.

fWorld_logic

G4LogicalVolume* DicomDetectorConstruction::fWorld_logic

protected

Definition at line 100 of file DicomDetectorConstruction.hh.

fWorld_phys

G4VPhysicalVolume* DicomDetectorConstruction::fWorld_phys

protected

Definition at line 101 of file DicomDetectorConstruction.hh.

fWorld_solid

G4Box* DicomDetectorConstruction::fWorld_solid

protected

Definition at line 99 of file DicomDetectorConstruction.hh.

fZSliceHeaderMerged

DicomPhantomZSliceHeader* DicomDetectorConstruction::fZSliceHeaderMerged

protected

Definition at line 120 of file DicomDetectorConstruction.hh.

fZSliceHeaders

std::vector<DicomPhantomZSliceHeader*> DicomDetectorConstruction::fZSliceHeaders

protected

Definition at line 118 of file DicomDetectorConstruction.hh.

---

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

• DicomDetectorConstruction.hh
• DicomDetectorConstruction.cc