IORTDetectorConstruction Class Reference

#include <IORTDetectorConstruction.hh>

Public Member Functions
	IORTDetectorConstruction (G4VPhysicalVolume *)
	~IORTDetectorConstruction ()
G4ThreeVector	GetDetectorToWorldPosition ()
G4ThreeVector	GetDetectorToPhantomPosition ()
void	SetDetectorPosition ()
bool	IsInside (G4double detectorX, G4double detectorY, G4double detectorZ, G4double phantomX, G4double phantomY, G4double phantomZ, G4ThreeVector detToPhantomPosition)
G4bool	SetPhantomMaterial (G4String material)
void	SetVoxelSize (G4double sizeX, G4double sizeY, G4double sizeZ)
void	SetDetectorSize (G4double sizeX, G4double sizeY, G4double sizeZ)
void	SetPhantomSize (G4double sizeX, G4double sizeY, G4double sizeZ)
void	SetPhantomPosition (G4ThreeVector)
void	SetDetectorToPhantomPosition (G4ThreeVector DetectorToPhantomPosition)
void	UpdateGeometry ()
void	DeleteDisc ()
void	ConstructDisc ()
void	PrintParameters ()
G4LogicalVolume *	GetDetectorLogicalVolume ()
G4bool	SetDiscoMaterialIORT (G4String material)
void	SetOuterRadiusDiscoIORT (G4double outerr)
void	SetinnerRadiusDiscoIORT (G4double innerr)
void	SetheightDiscoIORT (G4double height)
void	SetDiscoXPositionIORT (G4double xpos)
void	SetDiscoYPositionIORT (G4double ypos)
void	SetDiscoZPositionIORT (G4double zpos)
G4bool	SetDiscoMaterialIORT1 (G4String material)
void	SetOuterRadiusDiscoIORT1 (G4double outerr)
void	SetinnerRadiusDiscoIORT1 (G4double innerr)
void	SetheightDiscoIORT1 (G4double height)
void	SetDiscoXPositionIORT1 (G4double xpos)
void	SetAngleDiscoIORT0 (G4double phi0)

Detailed Description

Definition at line 58 of file IORTDetectorConstruction.hh.

Constructor & Destructor Documentation

IORTDetectorConstruction::IORTDetectorConstruction

(

G4VPhysicalVolume *

physicalTreatmentRoom

)

Definition at line 71 of file IORTDetectorConstruction.cc.

  : motherPhys(physicalTreatmentRoom), // pointer to WORLD volume 
    detectorSD(0), detectorROGeometry(0), matrix(0), 
    phantom(0), detector(0),
    phantomLogicalVolume(0), detectorLogicalVolume(0), 
    phantomPhysicalVolume(0), detectorPhysicalVolume(0),
    aRegion(0),

    solidDiscoIORT0(0),
    logicDiscoIORT0(0),
    physiDiscoIORT0(0),

    solidDiscoIORT(0),
    logicDiscoIORT(0),
    physiDiscoIORT(0),

    solidDiscoIORT1(0),
    logicDiscoIORT1(0),
    physiDiscoIORT1(0)

{
  IORTAnalysisManager::GetInstance();

  /* NOTE! that the IORTDetectorConstruction class
   * does NOT inherit from G4VUserDetectorConstruction G4 class
   * So the Construct() mandatory virtual method is inside another geometric class
   * like the collimatorXXBeamLIne, ...
   */

  // Messenger to change parameters of the phantom/detector geometry
  detectorMessenger = new IORTDetectorMessenger(this);

  // Default detector voxels size
  // 200 slabs along the beam direction (X)
  sizeOfVoxelAlongX = 0.5 *CLHEP::mm; // 
  sizeOfVoxelAlongY = 0.5 *CLHEP::mm; //  
  sizeOfVoxelAlongZ = 0.5 *CLHEP::mm; // 

  // Define here the material of the water phantom and of the detector
  SetPhantomMaterial("G4_WATER"); 
  // Construct geometry (messenger commands)
  SetDetectorSize(7.*CLHEP::cm, 15.*CLHEP::cm, 15.*CLHEP::cm);    
  SetPhantomSize(20. *CLHEP::cm, 20. *CLHEP::cm, 20. *CLHEP::cm);   
  SetPhantomPosition(G4ThreeVector(4.5 *CLHEP::cm, 0. *CLHEP::cm, 0. *CLHEP::cm)); 
  SetDetectorToPhantomPosition(G4ThreeVector(0. *CLHEP::cm, 2.5 *CLHEP::cm, 2.5 *CLHEP::cm));  

  // Default protection disc geometry and materials
  SetOuterRadiusDiscoIORT (40. *CLHEP::mm);  
  SetinnerRadiusDiscoIORT (0.*CLHEP::mm);   
  SetheightDiscoIORT (2.0*CLHEP::mm);        
  SetDiscoXPositionIORT (-11.0*CLHEP::mm);
  SetDiscoYPositionIORT (0.0*CLHEP::mm);
  SetDiscoZPositionIORT (0.0*CLHEP::mm);
  SetDiscoMaterialIORT("G4_WATER");   

  SetOuterRadiusDiscoIORT1 (40. *CLHEP::mm);  
  SetinnerRadiusDiscoIORT1 (0.*CLHEP::mm);   
  SetheightDiscoIORT1 (1.0*CLHEP::mm);        
  SetDiscoXPositionIORT1 (-8.0*CLHEP::mm);
  SetDiscoMaterialIORT1("G4_WATER");

  SetAngleDiscoIORT0 (90.0 *CLHEP::deg);    

  // Write virtual parameters to the real ones and check for consistency      
  UpdateGeometry();
}

Here is the call graph for this function:

IORTDetectorConstruction::~IORTDetectorConstruction

(

)

Definition at line 139 of file IORTDetectorConstruction.cc.

{ 
    delete detectorROGeometry;  
    delete matrix;  
    delete detectorMessenger;
}

Member Function Documentation

void IORTDetectorConstruction::ConstructDisc

(

)

Definition at line 255 of file IORTDetectorConstruction.cc.

{
// ---------------------------------------------------------------//
  //                    6.0 mm Protection Discs Volume          //
    // ---------------------------------------------------------------//
  const G4double startAngleDiscoIORT0 = 0.*CLHEP::deg;
  const G4double spanningAngleDiscoIORT0 = 360.*CLHEP::deg;

  //G4double phi0 = 180. *CLHEP::deg; // messenger    

  // Matrix definition for a rotation (deg).       
  G4RotationMatrix rm0;               
  rm0.rotateY(AngleDiscoIORT0);

  
  solidDiscoIORT0 = new G4Tubs("DiscoIORT0", innerRadiusDiscoIORT, 
                    OuterRadiusDiscoIORT,
                    (heightDiscoIORT + heightDiscoIORT1), 
                    startAngleDiscoIORT0, 
                    spanningAngleDiscoIORT0);

  G4LogicalVolume* logDiscoIORT0 = new G4LogicalVolume(solidDiscoIORT0, 
                                  detectorMaterial, "DiscoIORT0Log", 0, 0, 0);

  physiDiscoIORT0 = new G4PVPlacement(G4Transform3D(rm0, G4ThreeVector((DiscoXPositionIORT + heightDiscoIORT1),DiscoYPositionIORT,DiscoZPositionIORT)),
                       "DiscoIORT0Phys", logDiscoIORT0, detectorPhysicalVolume, false, 0); 

  white = new G4VisAttributes( G4Colour());
  white -> SetVisibility(true);
  // white -> SetForceSolid(true);
  logDiscoIORT0 -> SetVisAttributes(white);

// ---------------------------------------------------------------//
  //                    4.0 mm Aluminium Protection Disc          //
    // ---------------------------------------------------------------//
  //G4bool isotopes = false;
 // G4Material* leadNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_Pb", isotopes);
 // DiscoMaterialIORT = leadNist;   // messenger
  gray = new G4VisAttributes( G4Colour(0.5, 0.5, 0.5 ));
  gray-> SetVisibility(true);
  //gray -> SetForceWireframe(true);
  //gray-> SetForceSolid(true);

  gray1 = new G4VisAttributes( G4Colour(0.7, 0.7, 0.7 ));
  gray1-> SetVisibility(true);
  //gray1 -> SetForceWireframe(true);
  //gray1-> SetForceSolid(true);
 // const G4double OuterRadiusDiscoIORT = 35. *CLHEP::mm;  // messenger
 // const G4double innerRadiusDiscoIORT = 0.*CLHEP::mm;   // messenger
 //  const G4double heightDiscoIORT = 3.0*CLHEP::mm;        // messenger
  const G4double startAngleDiscoIORT = 0.*CLHEP::deg;
  const G4double spanningAngleDiscoIORT = 360.*CLHEP::deg;
 // const G4double DiscoXPositionIORT = -14.0*CLHEP::mm; // messenger 

//G4Material* DiscoMaterialIORT = G4NistManager::Instance()->FindOrBuildMaterial("G4_PLEXIGLASS", isotopes);// messenger 
  
   
  G4double phi = 0. *CLHEP::deg;     

  // Matrix definition for a 90 deg rotation. Also used for other volumes       
  G4RotationMatrix rm;               
  rm.rotateY(phi);

  
  solidDiscoIORT = new G4Tubs("DiscoIORT", innerRadiusDiscoIORT, 
                    OuterRadiusDiscoIORT,
                    heightDiscoIORT, 
                    startAngleDiscoIORT, 
                    spanningAngleDiscoIORT);

  G4LogicalVolume* logDiscoIORT = new G4LogicalVolume(solidDiscoIORT, 
                                  DiscoMaterialIORT, "DiscoIORTLog", 0, 0, 0);

  physiDiscoIORT = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(0.,0.,(- heightDiscoIORT1))),
                       "DiscoIORTPhys", logDiscoIORT, physiDiscoIORT0, false, 0); 
  
  logDiscoIORT -> SetVisAttributes(gray1);
  
  
      // ---------------------------------------------------------------//
      //             2.0 mm Lead Protection Disc                          //
      // ---------------------------------------------------------------//
  
 // const G4double OuterRadiusDiscoIORT1 = 35. *CLHEP::mm;
 // const G4double innerRadiusDiscoIORT1 = 0.*CLHEP::mm;
 // const G4double heightDiscoIORT1 = 0.5*CLHEP::mm;
  const G4double startAngleDiscoIORT1 = 0.*CLHEP::deg;
  const G4double spanningAngleDiscoIORT1 = 360.*CLHEP::deg;
//  const G4double DiscoXPositionIORT1 = -10.5*CLHEP::mm; messenger
//  G4Material* DiscoMaterialIORT1 = G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu", isotopes);// messenger  
  
       
  
  solidDiscoIORT1 = new G4Tubs("DiscoIORT1", innerRadiusDiscoIORT1, 
                    OuterRadiusDiscoIORT1,
                    heightDiscoIORT1, 
                    startAngleDiscoIORT1, 
                    spanningAngleDiscoIORT1);

  G4LogicalVolume* logDiscoIORT1 = new G4LogicalVolume(solidDiscoIORT1, 
                                  DiscoMaterialIORT1, "DiscoIORTLog1", 0, 0, 0);

  physiDiscoIORT1 = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(0.,0.,heightDiscoIORT)),
                       "DiscoIORTPhys1", logDiscoIORT1, physiDiscoIORT0, false, 0); 
  white = new G4VisAttributes( G4Colour());
  white -> SetVisibility(true);
  white -> SetForceSolid(true);
  logDiscoIORT1 -> SetVisAttributes(gray);

}

Here is the call graph for this function:

Here is the caller graph for this function:

void IORTDetectorConstruction::DeleteDisc

(

)

Definition at line 697 of file IORTDetectorConstruction.cc.

{
    delete solidDiscoIORT0;
    delete logicDiscoIORT0;
    delete physiDiscoIORT0;
    delete solidDiscoIORT;
    delete logicDiscoIORT;
    delete physiDiscoIORT;
    delete solidDiscoIORT1;
    delete logicDiscoIORT1;
    delete physiDiscoIORT1;
    G4RunManager::GetRunManager() -> GeometryHasBeenModified();
    G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
}

Here is the call graph for this function:

G4LogicalVolume* IORTDetectorConstruction::GetDetectorLogicalVolume ( )

inline

Definition at line 169 of file IORTDetectorConstruction.hh.

{ return detectorLogicalVolume;}

G4ThreeVector IORTDetectorConstruction::GetDetectorToPhantomPosition ( )

inline

Definition at line 84 of file IORTDetectorConstruction.hh.

{
    return G4ThreeVector(phantomSizeX/2 - detectorSizeX/2 + detectorPosition.getX(),
                         phantomSizeY/2 - detectorSizeY/2 + detectorPosition.getY(),
                         phantomSizeZ/2 - detectorSizeZ/2 + detectorPosition.getZ()
                  );
}

Here is the call graph for this function:

G4ThreeVector IORTDetectorConstruction::GetDetectorToWorldPosition ( )

inline

Definition at line 78 of file IORTDetectorConstruction.hh.

  {
    return phantomPosition + detectorPosition;
  }

Here is the caller graph for this function:

bool IORTDetectorConstruction::IsInside ( G4double  detectorX, G4double  detectorY, G4double  detectorZ, G4double  phantomX, G4double  phantomY, G4double  phantomZ, G4ThreeVector  detToPhantomPosition  )

inline

Definition at line 106 of file IORTDetectorConstruction.hh.

{
// Dimensions check... X Y and Z
// Firstly check what dimension we are modifying
    {
        if (detectorX > phantomX) 
         {
            G4cout << "Error: Detector X dimension must be smaller or equal to the corrispondent of the phantom" << G4endl;
            return false;
         }
        if ( (phantomX - detectorX) < detToPhantomPosition.getX()) 
             {
            G4cout << "Error: X dimension doesn't fit with detector to phantom relative position" << G4endl;
            return false;
             }
    }

    {
        if (detectorY > phantomY) 
         {
            G4cout << "Error: Detector Y dimension must be smaller or equal to the corrispondent of the phantom" << G4endl;
            return false;
         }
        if ( (phantomY - detectorY) < detToPhantomPosition.getY()) 
         {
           G4cout << "Error: Y dimension doesn't fit with detector to phantom relative position" << G4endl;
           return false;
         }
    }            

    {
        if (detectorZ > phantomZ) 
         {
            G4cout << "Error: Detector Z dimension must be smaller or equal to the corrispondent of the phantom" << G4endl;
            return false;
         }
        if ( (phantomZ - detectorZ) < detToPhantomPosition.getZ()) 
         {
           G4cout << "Error: Z dimension doesn't fit with detector to phantom relative position" << G4endl;
           return false;
         }
    }

    return true;
}

Here is the call graph for this function:

void IORTDetectorConstruction::PrintParameters

(

)

Definition at line 713 of file IORTDetectorConstruction.cc.

{

    G4cout << "The (X,Y,Z) dimensions of the phantom are : (" << 
    G4BestUnit( phantom -> GetXHalfLength()*2., "Length") << ',' << 
    G4BestUnit( phantom -> GetYHalfLength()*2., "Length") << ',' << 
    G4BestUnit( phantom -> GetZHalfLength()*2., "Length") << ')' << G4endl; 
    
    G4cout << "The (X,Y,Z) dimensions of the detector are : (" << 
    G4BestUnit( detector -> GetXHalfLength()*2., "Length") << ',' << 
    G4BestUnit( detector -> GetYHalfLength()*2., "Length") << ',' << 
    G4BestUnit( detector -> GetZHalfLength()*2., "Length") << ')' << G4endl; 

    G4cout << "Displacement between Phantom and World is: "; 
    G4cout << "DX= "<< G4BestUnit(phantomPosition.getX(),"Length") << 
    "DY= "<< G4BestUnit(phantomPosition.getY(),"Length") << 
    "DZ= "<< G4BestUnit(phantomPosition.getZ(),"Length") << G4endl;

    G4cout << "The (X,Y,Z) sizes of the Voxels are: (" << 
    G4BestUnit(sizeOfVoxelAlongX, "Length")  << ',' << 
    G4BestUnit(sizeOfVoxelAlongY, "Length")  << ',' << 
    G4BestUnit(sizeOfVoxelAlongZ, "Length") << ')' << G4endl;

    G4cout << "The number of Voxels along (X,Y,Z) is: (" << 
    numberOfVoxelsAlongX  << ',' <<
    numberOfVoxelsAlongY  <<','  <<
    numberOfVoxelsAlongZ  << ')' << G4endl;

}

Here is the call graph for this function:

Here is the caller graph for this function:

void IORTDetectorConstruction::SetAngleDiscoIORT0

(

G4double

phi0

)

Definition at line 618 of file IORTDetectorConstruction.cc.

{
    
    AngleDiscoIORT0 = phi0;
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetDetectorPosition ( )

inline

Definition at line 94 of file IORTDetectorConstruction.hh.

  {
      // Adjust detector position
      detectorPosition.setX(detectorToPhantomPosition.getX() - phantomSizeX/2 + detectorSizeX/2);
      detectorPosition.setY(detectorToPhantomPosition.getY() - phantomSizeY/2 + detectorSizeY/2);
      detectorPosition.setZ(detectorToPhantomPosition.getZ() - phantomSizeZ/2 + detectorSizeZ/2);
     
    //G4cout << "*************** DetectorToPhantomPosition " << detectorToPhantomPosition/cm << "\n";
    //G4cout << "*************** DetectorPosition " << detectorPosition/cm << "\n";
  }

Here is the call graph for this function:

Here is the caller graph for this function:

void IORTDetectorConstruction::SetDetectorSize	(	G4double	sizeX,
		G4double	sizeY,
		G4double	sizeZ
	)

Definition at line 528 of file IORTDetectorConstruction.cc.

{
    if (sizeX > 0.) {detectorSizeX = sizeX;}
    if (sizeY > 0.) {detectorSizeY = sizeY;}
    if (sizeZ > 0.) {detectorSizeZ = sizeZ;}
    SetVoxelSize(sizeOfVoxelAlongX, sizeOfVoxelAlongY, sizeOfVoxelAlongZ);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void IORTDetectorConstruction::SetDetectorToPhantomPosition

(

G4ThreeVector

DetectorToPhantomPosition

)

Definition at line 549 of file IORTDetectorConstruction.cc.

{
    detectorToPhantomPosition = displ;
}

Here is the caller graph for this function:

G4bool IORTDetectorConstruction::SetDiscoMaterialIORT

(

G4String

material

)

Definition at line 466 of file IORTDetectorConstruction.cc.

{

    if (G4Material* dMat = G4NistManager::Instance()->FindOrBuildMaterial(material, false) )
    {
    DiscoMaterialIORT  = dMat;
    
    if (logicDiscoIORT) 
    {
        logicDiscoIORT -> SetMaterial(dMat); 
        
        G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
        G4RunManager::GetRunManager() -> GeometryHasBeenModified();
        G4cout << "The material of Protection disc 1 has been changed to " << material << G4endl;
    }
    }
    else
    {
    G4cout << "WARNING: material \"" << material << "\" doesn't exist in NIST elements/materials"
        " table [located in $G4INSTALL/source/materials/src/G4NistMaterialBuilder.cc]" << G4endl; 
    G4cout << "Use command \"/parameter/nist\" to see full materials list!" << G4endl; 
    return false;
    }

    return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool IORTDetectorConstruction::SetDiscoMaterialIORT1

(

G4String

material

)

Definition at line 493 of file IORTDetectorConstruction.cc.

{

    if (G4Material* d1Mat = G4NistManager::Instance()->FindOrBuildMaterial(material, false) )
    {
    DiscoMaterialIORT1  = d1Mat;
    
    if (logicDiscoIORT1) 
    {
        logicDiscoIORT1 -> SetMaterial(d1Mat); 
        
        G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
        G4RunManager::GetRunManager() -> GeometryHasBeenModified();
        G4cout << "The material of Protection disc 2 has been changed to " << material << G4endl;
    }
    }
    else
    {
    G4cout << "WARNING: material \"" << material << "\" doesn't exist in NIST elements/materials"
        " table [located in $G4INSTALL/source/materials/src/G4NistMaterialBuilder.cc]" << G4endl; 
    G4cout << "Use command \"/parameter/nist\" to see full materials list!" << G4endl; 
    return false;
    }

    return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void IORTDetectorConstruction::SetDiscoXPositionIORT

(

G4double

xpos

)

Definition at line 573 of file IORTDetectorConstruction.cc.

{
    
    DiscoXPositionIORT = xpos;
        
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetDiscoXPositionIORT1

(

G4double

xpos

)

Definition at line 612 of file IORTDetectorConstruction.cc.

{
    
    DiscoXPositionIORT1 = xpos;
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetDiscoYPositionIORT

(

G4double

ypos

)

Definition at line 580 of file IORTDetectorConstruction.cc.

{
    
    DiscoYPositionIORT = ypos;
        
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetDiscoZPositionIORT

(

G4double

zpos

)

Definition at line 587 of file IORTDetectorConstruction.cc.

{
    
    DiscoZPositionIORT = zpos;
        
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetheightDiscoIORT

(

G4double

height

)

Definition at line 567 of file IORTDetectorConstruction.cc.

{
    if (height > 0.) {heightDiscoIORT = height;}
    
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetheightDiscoIORT1

(

G4double

height

)

Definition at line 606 of file IORTDetectorConstruction.cc.

{
    if (height > 0.) {heightDiscoIORT1 = height;}
    
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetinnerRadiusDiscoIORT

(

G4double

innerr

)

Definition at line 561 of file IORTDetectorConstruction.cc.

{
    if (innerr > 0.) {innerRadiusDiscoIORT = innerr;}
    
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetinnerRadiusDiscoIORT1

(

G4double

innerr

)

Definition at line 600 of file IORTDetectorConstruction.cc.

{
    if (innerr > 0.) {innerRadiusDiscoIORT1 = innerr;}
    
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetOuterRadiusDiscoIORT

(

G4double

outerr

)

Definition at line 555 of file IORTDetectorConstruction.cc.

{
    if (outerr > 0.) {OuterRadiusDiscoIORT = outerr;}
    
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetOuterRadiusDiscoIORT1

(

G4double

outerr

)

Definition at line 594 of file IORTDetectorConstruction.cc.

{
    if (outerr > 0.) {OuterRadiusDiscoIORT1 = outerr;}
    
}

Here is the caller graph for this function:

G4bool IORTDetectorConstruction::SetPhantomMaterial

(

G4String

material

)

Definition at line 438 of file IORTDetectorConstruction.cc.

{

    if (G4Material* pMat = G4NistManager::Instance()->FindOrBuildMaterial(material, false) )
    {
    phantomMaterial  = pMat;
    detectorMaterial = pMat;
    if (detectorLogicalVolume && phantomLogicalVolume) 
    {
        detectorLogicalVolume -> SetMaterial(pMat); 
        phantomLogicalVolume ->  SetMaterial(pMat);

        G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
        G4RunManager::GetRunManager() -> GeometryHasBeenModified();
        G4cout << "The material of Phantom/Detector has been changed to " << material << G4endl;
    }
    }
    else
    {
    G4cout << "WARNING: material \"" << material << "\" doesn't exist in NIST elements/materials"
        " table [located in $G4INSTALL/source/materials/src/G4NistMaterialBuilder.cc]" << G4endl; 
    G4cout << "Use command \"/parameter/nist\" to see full materials list!" << G4endl; 
    return false;
    }

    return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void IORTDetectorConstruction::SetPhantomPosition

(

G4ThreeVector

pos

)

Definition at line 543 of file IORTDetectorConstruction.cc.

{
    phantomPosition = pos;
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetPhantomSize	(	G4double	sizeX,
		G4double	sizeY,
		G4double	sizeZ
	)

Definition at line 520 of file IORTDetectorConstruction.cc.

{
    if (sizeX > 0.) phantomSizeX = sizeX;
    if (sizeY > 0.) phantomSizeY = sizeY;
    if (sizeZ > 0.) phantomSizeZ = sizeZ;
}

Here is the caller graph for this function:

void IORTDetectorConstruction::SetVoxelSize	(	G4double	sizeX,
		G4double	sizeY,
		G4double	sizeZ
	)

Definition at line 537 of file IORTDetectorConstruction.cc.

{
    if (sizeX > 0.) {sizeOfVoxelAlongX = sizeX;}
    if (sizeY > 0.) {sizeOfVoxelAlongY = sizeY;}
    if (sizeZ > 0.) {sizeOfVoxelAlongZ = sizeZ;}
}

Here is the caller graph for this function:

void IORTDetectorConstruction::UpdateGeometry

(

)

Definition at line 628 of file IORTDetectorConstruction.cc.

{
    /* 
     * Check parameters consistency
     */
    ParametersCheck();

    G4GeometryManager::GetInstance() -> OpenGeometry();
    if (phantom)
    {
    phantom -> SetXHalfLength(phantomSizeX/2);
    phantom -> SetYHalfLength(phantomSizeY/2);
    phantom -> SetZHalfLength(phantomSizeZ/2);
    phantomPhysicalVolume -> SetTranslation(phantomPosition);
    }
    else   ConstructPhantom();

    // Get the center of the detector 
    SetDetectorPosition();
    if (detector)
    {
    detector -> SetXHalfLength(detectorSizeX/2);
    detector -> SetYHalfLength(detectorSizeY/2);
    detector -> SetZHalfLength(detectorSizeZ/2);
    detectorPhysicalVolume -> SetTranslation(detectorPosition);      
    }
    else    ConstructDetector();

    // update disc function
    delete solidDiscoIORT0;
    delete logicDiscoIORT0;
    delete physiDiscoIORT0;
    delete solidDiscoIORT;
    delete logicDiscoIORT;
    delete physiDiscoIORT;
    delete solidDiscoIORT1;
    delete logicDiscoIORT1;
    delete physiDiscoIORT1;
    ConstructDisc();
    

    // Round to nearest integer number of voxel 
    numberOfVoxelsAlongX = G4lrint(detectorSizeX / sizeOfVoxelAlongX);
    sizeOfVoxelAlongX = ( detectorSizeX / numberOfVoxelsAlongX );

    numberOfVoxelsAlongY = G4lrint(detectorSizeY / sizeOfVoxelAlongY);
    sizeOfVoxelAlongY = ( detectorSizeY / numberOfVoxelsAlongY );

    numberOfVoxelsAlongZ = G4lrint(detectorSizeZ / sizeOfVoxelAlongZ);
    sizeOfVoxelAlongZ = ( detectorSizeZ / numberOfVoxelsAlongZ );

    //G4cout << "*************** DetectorToWorldPosition " << GetDetectorToWorldPosition()/cm << "\n";
    ConstructSensitiveDetector(GetDetectorToWorldPosition());

    volumeOfVoxel = sizeOfVoxelAlongX * sizeOfVoxelAlongY * sizeOfVoxelAlongZ;
    massOfVoxel = detectorMaterial -> GetDensity() * volumeOfVoxel;
    //  This will clear the existing matrix (together with all data inside it)! 
    matrix = IORTMatrix::GetInstance(numberOfVoxelsAlongX, 
        numberOfVoxelsAlongY,
        numberOfVoxelsAlongZ,
        massOfVoxel);

    // Inform the kernel about the new geometry
    G4RunManager::GetRunManager() -> GeometryHasBeenModified();
    G4RunManager::GetRunManager() -> PhysicsHasBeenModified();

    PrintParameters();
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p02/examples/advanced/iort_therapy/include/IORTDetectorConstruction.hh
• source/geant4.10.03.p02/examples/advanced/iort_therapy/src/IORTDetectorConstruction.cc