#include <IORTDetectorConstruction.hh>

Collaboration diagram for IORTDetectorConstruction:

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)

Private Member Functions
void	ConstructPhantom ()

void	ConstructDetector ()

void	ConstructSensitiveDetector (G4ThreeVector positionToWORLD)

void	ParametersCheck ()

Private Attributes
IORTDetectorMessenger *	detectorMessenger

G4VisAttributes *	red

G4VPhysicalVolume *	motherPhys

IORTDetectorSD *	detectorSD

IORTDetectorROGeometry *	detectorROGeometry

IORTMatrix *	matrix

G4Box *	phantom

G4Box *	detector

G4LogicalVolume *	phantomLogicalVolume

G4LogicalVolume *	detectorLogicalVolume

G4VPhysicalVolume *	phantomPhysicalVolume

G4VPhysicalVolume *	detectorPhysicalVolume

G4double	phantomSizeX

G4double	phantomSizeY

G4double	phantomSizeZ

G4double	detectorSizeX

G4double	detectorSizeY

G4double	detectorSizeZ

G4ThreeVector	phantomPosition

G4ThreeVector	detectorPosition

G4ThreeVector	detectorToPhantomPosition

G4double	sizeOfVoxelAlongX

G4double	sizeOfVoxelAlongY

G4double	sizeOfVoxelAlongZ

G4int	numberOfVoxelsAlongX

G4int	numberOfVoxelsAlongY

G4int	numberOfVoxelsAlongZ

G4double	volumeOfVoxel

G4double	massOfVoxel

G4Material *	phantomMaterial

G4Material *	detectorMaterial

G4Region *	aRegion

G4Tubs *	solidDiscoIORT0

G4LogicalVolume *	logicDiscoIORT0

G4VPhysicalVolume *	physiDiscoIORT0

G4double	AngleDiscoIORT0

G4VisAttributes *	white

G4VisAttributes *	gray

G4VisAttributes *	gray1

G4double	innerRadiusDiscoIORT

G4double	OuterRadiusDiscoIORT

G4double	heightDiscoIORT

G4double	DiscoXPositionIORT

G4double	DiscoYPositionIORT

G4double	DiscoZPositionIORT

G4Tubs *	solidDiscoIORT

G4LogicalVolume *	logicDiscoIORT

G4VPhysicalVolume *	physiDiscoIORT

G4Material *	DiscoMaterialIORT

G4double	innerRadiusDiscoIORT1

G4double	OuterRadiusDiscoIORT1

G4double	heightDiscoIORT1

G4double	DiscoXPositionIORT1

G4Tubs *	solidDiscoIORT1

G4LogicalVolume *	logicDiscoIORT1

G4VPhysicalVolume *	physiDiscoIORT1

G4Material *	DiscoMaterialIORT1

Detailed Description

Definition at line 58 of file IORTDetectorConstruction.hh.

Constructor & Destructor Documentation

◆ IORTDetectorConstruction()

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::~IORTDetectorConstruction ( )

Definition at line 139 of file IORTDetectorConstruction.cc.

 { 
     delete detectorROGeometry;  
     delete matrix;  
     delete detectorMessenger;
 }

Member Function Documentation

◆ ConstructDetector()

void IORTDetectorConstruction::ConstructDetector ( )

private

Definition at line 204 of file IORTDetectorConstruction.cc.

 {
 
     // Definition of the solid volume of the Detector
     detector = new G4Box("Detector", 
              detectorSizeX/2, 
              detectorSizeY/2, 
              detectorSizeZ/2);
     
     // Definition of the logic volume of the Phantom
     detectorLogicalVolume = new G4LogicalVolume(detector,
                         detectorMaterial,
                         "DetectorLog",
                         0,0,0);
 // Definition of the physical volume of the Phantom 
     detectorPhysicalVolume = new G4PVPlacement(0, 
                            detectorPosition, // Setted by displacement 
                            "DetectorPhys", 
                            detectorLogicalVolume, 
                            phantomPhysicalVolume, 
                            false,0);
 
 // Visualisation attributes of the detector 
     //skyBlue = new G4VisAttributes( G4Colour(135/255. , 206/255. ,  235/255. ));
     G4VisAttributes * skyBlue1 = new G4VisAttributes( G4Colour(135/255. , 206/255. ,  235/255. ));
     //skyBlue1 -> SetForceWireframe(true);
     //skyBlue1 -> SetForceSolid(true);
     //skyBlue -> SetVisibility(true);
     //skyBlue -> SetForceSolid(true);
     //skyBlue -> SetForceWireframe(false);
     //detectorLogicalVolume -> SetVisAttributes(skyBlue);
     detectorLogicalVolume -> SetVisAttributes(skyBlue1);
     
    // detectorLogicalVolume -> SetVisAttributes(G4VisAttributes::Invisible);
 
 
   // **************
   // Cut per Region    
   // **************
   
   // A smaller cut is fixed in the phantom to calculate the energy deposit with the
   // required accuracy 
     if (!aRegion)
     {
     aRegion = new G4Region("DetectorLog");
     detectorLogicalVolume -> SetRegion(aRegion);
     aRegion -> AddRootLogicalVolume(detectorLogicalVolume);
     }
 
 }

Here is the caller graph for this function:

◆ ConstructDisc()

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:

◆ ConstructPhantom()

void IORTDetectorConstruction::ConstructPhantom ( )

private

Definition at line 151 of file IORTDetectorConstruction.cc.

 {
     // Definition of the solid volume of the Phantom
     phantom = new G4Box("Phantom", 
             phantomSizeX/2, 
             phantomSizeY/2, 
             phantomSizeZ/2);
     
 // Definition of the logical volume of the Phantom
     phantomLogicalVolume = new G4LogicalVolume(phantom, 
                          phantomMaterial, 
                          "phantomLog", 0, 0, 0);
   
     // Definition of the physics volume of the Phantom
     phantomPhysicalVolume = new G4PVPlacement(0,
                                         phantomPosition,
                         "phantomPhys",
                         phantomLogicalVolume,
                         motherPhys,
                         false,
                         0);
 
 // Visualisation attributes of the phantom
     red = new G4VisAttributes(G4Colour(255/255., 0/255. ,0/255.));
     red -> SetVisibility(true);
     //red -> SetForceSolid(true);
     //red -> SetForceWireframe(true);
     phantomLogicalVolume -> SetVisAttributes(red); 
     //phantomLogicalVolume -> SetVisAttributes(G4VisAttributes::Invisible);
 }

Here is the caller graph for this function:

◆ ConstructSensitiveDetector()

void IORTDetectorConstruction::ConstructSensitiveDetector ( G4ThreeVector positionToWORLD )

private

Definition at line 367 of file IORTDetectorConstruction.cc.

 {  
     // Install new Sensitive Detector and ROGeometry 
     delete detectorROGeometry; // this should be safe in C++ also if we have a NULL pointer
     //if (detectorSD) detectorSD->PrintAll();
     //delete detectorSD;
     // Sensitive Detector and ReadOut geometry definition
     G4SDManager* sensitiveDetectorManager = G4SDManager::GetSDMpointer();
 
     static G4String sensitiveDetectorName = "Detector"; 
     if (!detectorSD)
     {
         // The sensitive detector is instantiated
         detectorSD = new IORTDetectorSD(sensitiveDetectorName);
     }
     // The Read Out Geometry is instantiated
     static G4String ROGeometryName = "DetectorROGeometry";
     detectorROGeometry = new IORTDetectorROGeometry(ROGeometryName,
                                 detectorToWorldPosition,
                                 detectorSizeX/2,  // controllare che sia necessario /2
                                 detectorSizeY/2,  // CONFERMATO!!! ci vuole!!!
                                 detectorSizeZ/2,
                                 numberOfVoxelsAlongX,
                                 numberOfVoxelsAlongY,
                                 numberOfVoxelsAlongZ);
 
     G4cout << "Instantiating new Read Out Geometry \"" << ROGeometryName << "\""<< G4endl;
     // This will invoke Build() IORTDetectorROGeometry virtual method 
     detectorROGeometry -> BuildROGeometry();
     // Attach ROGeometry to SDetector
     detectorSD -> SetROgeometry(detectorROGeometry);
     //sensitiveDetectorManager -> Activate(sensitiveDetectorName, true);
     if (!sensitiveDetectorManager -> FindSensitiveDetector(sensitiveDetectorName, false))
     {
         G4cout << "Registering new DetectorSD \"" << sensitiveDetectorName << "\""<< G4endl;
         // Register user SD
         sensitiveDetectorManager -> AddNewDetector(detectorSD);
         // Attach SD to detector logical volume
         detectorLogicalVolume -> SetSensitiveDetector(detectorSD);
     }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ DeleteDisc()

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:

Here is the caller graph for this function:

◆ GetDetectorLogicalVolume()

G4LogicalVolume* IORTDetectorConstruction::GetDetectorLogicalVolume ( )

inline

Definition at line 169 of file IORTDetectorConstruction.hh.

169 { return detectorLogicalVolume;}

IORTDetectorConstruction::detectorLogicalVolume

G4LogicalVolume * detectorLogicalVolume

Definition: IORTDetectorConstruction.hh:200

Here is the call graph for this function:

◆ GetDetectorToPhantomPosition()

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:

◆ GetDetectorToWorldPosition()

G4ThreeVector IORTDetectorConstruction::GetDetectorToWorldPosition ( )

inline

Definition at line 78 of file IORTDetectorConstruction.hh.

   {
     return phantomPosition + detectorPosition;
   }

Here is the caller graph for this function:

◆ IsInside()

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:

Here is the caller graph for this function:

◆ ParametersCheck()

void IORTDetectorConstruction::ParametersCheck ( )

private

Definition at line 408 of file IORTDetectorConstruction.cc.

 {
     // Check phantom/detector sizes & relative position
     if (!IsInside(detectorSizeX, 
         detectorSizeY, 
         detectorSizeZ,
         phantomSizeX,
         phantomSizeY,
         phantomSizeZ,
         detectorToPhantomPosition
         ))
       G4Exception("IORTDetectorConstruction::ParametersCheck()", "IORT0001", FatalException, "Error: Detector is not fully inside Phantom!");
 
     // Check Detector sizes respect to the voxel ones
 
     if ( detectorSizeX < sizeOfVoxelAlongX) {
       G4Exception("IORTDetectorConstruction::ParametersCheck()", "IORT0002", FatalException, "Error: Detector X size must be bigger or equal than that of Voxel X");
     }
     if ( detectorSizeY < sizeOfVoxelAlongY) {
       G4Exception("IORTDetectorConstruction::ParametersCheck()", "IORT0003", FatalException, "Error: Detector X size must be bigger or equal than that of Voxel Y");    
     }
     if ( detectorSizeZ < sizeOfVoxelAlongZ) {
       G4Exception("IORTDetectorConstruction::ParametersCheck()", "IORT0004", FatalException, "Error: Detector X size must be bigger or equal than that of Voxel Z");
     }
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ PrintParameters()

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:

◆ SetAngleDiscoIORT0()

void IORTDetectorConstruction::SetAngleDiscoIORT0 ( G4double phi0 )

Definition at line 618 of file IORTDetectorConstruction.cc.

 {
     
     AngleDiscoIORT0 = phi0;
 }

Here is the caller graph for this function:

◆ SetDetectorPosition()

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:

◆ SetDetectorSize()

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:

◆ SetDetectorToPhantomPosition()

void IORTDetectorConstruction::SetDetectorToPhantomPosition ( G4ThreeVector DetectorToPhantomPosition )

Definition at line 549 of file IORTDetectorConstruction.cc.

 {
     detectorToPhantomPosition = displ;
 }

Here is the caller graph for this function:

◆ SetDiscoMaterialIORT()

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:

◆ SetDiscoMaterialIORT1()

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:

◆ SetDiscoXPositionIORT()

void IORTDetectorConstruction::SetDiscoXPositionIORT ( G4double xpos )

Definition at line 573 of file IORTDetectorConstruction.cc.

 {
     
     DiscoXPositionIORT = xpos;
         
 }

Here is the caller graph for this function:

◆ SetDiscoXPositionIORT1()

void IORTDetectorConstruction::SetDiscoXPositionIORT1 ( G4double xpos )

Definition at line 612 of file IORTDetectorConstruction.cc.

 {
     
     DiscoXPositionIORT1 = xpos;
 }

Here is the caller graph for this function:

◆ SetDiscoYPositionIORT()

void IORTDetectorConstruction::SetDiscoYPositionIORT ( G4double ypos )

Definition at line 580 of file IORTDetectorConstruction.cc.

 {
     
     DiscoYPositionIORT = ypos;
         
 }

Here is the caller graph for this function:

◆ SetDiscoZPositionIORT()

void IORTDetectorConstruction::SetDiscoZPositionIORT ( G4double zpos )

Definition at line 587 of file IORTDetectorConstruction.cc.

 {
     
     DiscoZPositionIORT = zpos;
         
 }

Here is the caller graph for this function:

◆ SetheightDiscoIORT()

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:

◆ SetheightDiscoIORT1()

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:

◆ SetinnerRadiusDiscoIORT()

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:

◆ SetinnerRadiusDiscoIORT1()

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:

◆ SetOuterRadiusDiscoIORT()

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:

◆ SetOuterRadiusDiscoIORT1()

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:

◆ SetPhantomMaterial()

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:

◆ SetPhantomPosition()

void IORTDetectorConstruction::SetPhantomPosition ( G4ThreeVector pos )

Definition at line 543 of file IORTDetectorConstruction.cc.

 {
     phantomPosition = pos;
 }

Here is the caller graph for this function:

◆ SetPhantomSize()

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:

◆ SetVoxelSize()

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:

◆ UpdateGeometry()

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();
 }