Geant4.10
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions
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.

References IORTAnalysisManager::GetInstance(), SetAngleDiscoIORT0(), SetDetectorSize(), SetDetectorToPhantomPosition(), SetDiscoMaterialIORT(), SetDiscoMaterialIORT1(), SetDiscoXPositionIORT(), SetDiscoXPositionIORT1(), SetDiscoYPositionIORT(), SetDiscoZPositionIORT(), SetheightDiscoIORT(), SetheightDiscoIORT1(), SetinnerRadiusDiscoIORT(), SetinnerRadiusDiscoIORT1(), SetOuterRadiusDiscoIORT(), SetOuterRadiusDiscoIORT1(), SetPhantomMaterial(), SetPhantomPosition(), SetPhantomSize(), and UpdateGeometry().

72  : motherPhys(physicalTreatmentRoom), // pointer to WORLD volume
73  detectorSD(0), detectorROGeometry(0), matrix(0),
74  phantom(0), detector(0),
75  phantomLogicalVolume(0), detectorLogicalVolume(0),
76  phantomPhysicalVolume(0), detectorPhysicalVolume(0),
77  aRegion(0),
78 
79  solidDiscoIORT0(0),
80  logicDiscoIORT0(0),
81  physiDiscoIORT0(0),
82 
83  solidDiscoIORT(0),
84  logicDiscoIORT(0),
85  physiDiscoIORT(0),
86 
87  solidDiscoIORT1(0),
88  logicDiscoIORT1(0),
89  physiDiscoIORT1(0)
90 
91 {
92  IORTAnalysisManager::GetInstance();
93 
94  /* NOTE! that the IORTDetectorConstruction class
95  * does NOT inherit from G4VUserDetectorConstruction G4 class
96  * So the Construct() mandatory virtual method is inside another geometric class
97  * like the collimatorXXBeamLIne, ...
98  */
99 
100  // Messenger to change parameters of the phantom/detector geometry
101  detectorMessenger = new IORTDetectorMessenger(this);
102 
103  // Default detector voxels size
104  // 200 slabs along the beam direction (X)
105  sizeOfVoxelAlongX = 0.5 *CLHEP::mm; //
106  sizeOfVoxelAlongY = 0.5 *CLHEP::mm; //
107  sizeOfVoxelAlongZ = 0.5 *CLHEP::mm; //
108 
109  // Define here the material of the water phantom and of the detector
110  SetPhantomMaterial("G4_WATER");
111  // Construct geometry (messenger commands)
112  SetDetectorSize(7.*CLHEP::cm, 15.*CLHEP::cm, 15.*CLHEP::cm);
113  SetPhantomSize(20. *CLHEP::cm, 20. *CLHEP::cm, 20. *CLHEP::cm);
114  SetPhantomPosition(G4ThreeVector(4.5 *CLHEP::cm, 0. *CLHEP::cm, 0. *CLHEP::cm));
115  SetDetectorToPhantomPosition(G4ThreeVector(0. *CLHEP::cm, 2.5 *CLHEP::cm, 2.5 *CLHEP::cm));
116 
117  // Default protection disc geometry and materials
118  SetOuterRadiusDiscoIORT (40. *CLHEP::mm);
119  SetinnerRadiusDiscoIORT (0.*CLHEP::mm);
120  SetheightDiscoIORT (2.0*CLHEP::mm);
121  SetDiscoXPositionIORT (-11.0*CLHEP::mm);
122  SetDiscoYPositionIORT (0.0*CLHEP::mm);
123  SetDiscoZPositionIORT (0.0*CLHEP::mm);
124  SetDiscoMaterialIORT("G4_WATER");
125 
126  SetOuterRadiusDiscoIORT1 (40. *CLHEP::mm);
127  SetinnerRadiusDiscoIORT1 (0.*CLHEP::mm);
128  SetheightDiscoIORT1 (1.0*CLHEP::mm);
129  SetDiscoXPositionIORT1 (-8.0*CLHEP::mm);
130  SetDiscoMaterialIORT1("G4_WATER");
131 
132  SetAngleDiscoIORT0 (90.0 *CLHEP::deg);
133 
134  // Write virtual parameters to the real ones and check for consistency
135  UpdateGeometry();
136 }
IORTDetectorConstruction::SetDiscoZPositionIORT
void SetDiscoZPositionIORT(G4double zpos)
Definition: IORTDetectorConstruction.cc:587
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
IORTDetectorConstruction::SetOuterRadiusDiscoIORT1
void SetOuterRadiusDiscoIORT1(G4double outerr)
Definition: IORTDetectorConstruction.cc:594
IORTDetectorConstruction::SetDetectorToPhantomPosition
void SetDetectorToPhantomPosition(G4ThreeVector DetectorToPhantomPosition)
Definition: IORTDetectorConstruction.cc:549
IORTDetectorConstruction::SetPhantomSize
void SetPhantomSize(G4double sizeX, G4double sizeY, G4double sizeZ)
Definition: IORTDetectorConstruction.cc:520
IORTDetectorConstruction::SetDetectorSize
void SetDetectorSize(G4double sizeX, G4double sizeY, G4double sizeZ)
Definition: IORTDetectorConstruction.cc:528
IORTDetectorConstruction::SetheightDiscoIORT1
void SetheightDiscoIORT1(G4double height)
Definition: IORTDetectorConstruction.cc:606
IORTDetectorConstruction::SetinnerRadiusDiscoIORT1
void SetinnerRadiusDiscoIORT1(G4double innerr)
Definition: IORTDetectorConstruction.cc:600
IORTDetectorConstruction::SetDiscoMaterialIORT
G4bool SetDiscoMaterialIORT(G4String material)
Definition: IORTDetectorConstruction.cc:466
IORTDetectorConstruction::SetDiscoYPositionIORT
void SetDiscoYPositionIORT(G4double ypos)
Definition: IORTDetectorConstruction.cc:580
IORTDetectorConstruction::SetPhantomMaterial
G4bool SetPhantomMaterial(G4String material)
Definition: IORTDetectorConstruction.cc:438
IORTDetectorConstruction::SetDiscoXPositionIORT1
void SetDiscoXPositionIORT1(G4double xpos)
Definition: IORTDetectorConstruction.cc:612
IORTDetectorConstruction::SetAngleDiscoIORT0
void SetAngleDiscoIORT0(G4double phi0)
Definition: IORTDetectorConstruction.cc:618
IORTDetectorConstruction::SetPhantomPosition
void SetPhantomPosition(G4ThreeVector)
Definition: IORTDetectorConstruction.cc:543
IORTDetectorConstruction::SetinnerRadiusDiscoIORT
void SetinnerRadiusDiscoIORT(G4double innerr)
Definition: IORTDetectorConstruction.cc:561
IORTAnalysisManager::GetInstance
static IORTAnalysisManager * GetInstance()
Definition: IORTAnalysisManager.cc:81
IORTDetectorConstruction::SetheightDiscoIORT
void SetheightDiscoIORT(G4double height)
Definition: IORTDetectorConstruction.cc:567
IORTDetectorConstruction::SetOuterRadiusDiscoIORT
void SetOuterRadiusDiscoIORT(G4double outerr)
Definition: IORTDetectorConstruction.cc:555
IORTDetectorConstruction::SetDiscoMaterialIORT1
G4bool SetDiscoMaterialIORT1(G4String material)
Definition: IORTDetectorConstruction.cc:493
IORTDetectorConstruction::SetDiscoXPositionIORT
void SetDiscoXPositionIORT(G4double xpos)
Definition: IORTDetectorConstruction.cc:573
IORTDetectorConstruction::~IORTDetectorConstruction ( )

Definition at line 139 of file IORTDetectorConstruction.cc.

140 {
141  delete detectorROGeometry;
142  delete matrix;
143  delete detectorMessenger;
144 }

Member Function Documentation

void IORTDetectorConstruction::ConstructDisc ( )

Definition at line 255 of file IORTDetectorConstruction.cc.

References CLHEP::HepRotation::rotateY().

Referenced by UpdateGeometry().

256 {
257 // ---------------------------------------------------------------//
258  // 6.0 mm Protection Discs Volume //
259  // ---------------------------------------------------------------//
260  const G4double startAngleDiscoIORT0 = 0.*CLHEP::deg;
261  const G4double spanningAngleDiscoIORT0 = 360.*CLHEP::deg;
262 
263  //G4double phi0 = 180. *CLHEP::deg; // messenger
264 
265  // Matrix definition for a rotation (deg).
266  G4RotationMatrix rm0;
267  rm0.rotateY(AngleDiscoIORT0);
268 
269 
270  solidDiscoIORT0 = new G4Tubs("DiscoIORT0", innerRadiusDiscoIORT,
271  OuterRadiusDiscoIORT,
272  (heightDiscoIORT + heightDiscoIORT1),
273  startAngleDiscoIORT0,
274  spanningAngleDiscoIORT0);
275 
276  G4LogicalVolume* logDiscoIORT0 = new G4LogicalVolume(solidDiscoIORT0,
277  detectorMaterial, "DiscoIORT0Log", 0, 0, 0);
278 
279  physiDiscoIORT0 = new G4PVPlacement(G4Transform3D(rm0, G4ThreeVector((DiscoXPositionIORT + heightDiscoIORT1),DiscoYPositionIORT,DiscoZPositionIORT)),
280  "DiscoIORT0Phys", logDiscoIORT0, detectorPhysicalVolume, false, 0);
281 
282  white = new G4VisAttributes( G4Colour());
283  white -> SetVisibility(true);
284  // white -> SetForceSolid(true);
285  logDiscoIORT0 -> SetVisAttributes(white);
286 
287 // ---------------------------------------------------------------//
288  // 4.0 mm Aluminium Protection Disc //
289  // ---------------------------------------------------------------//
290  //G4bool isotopes = false;
291  // G4Material* leadNist = G4NistManager::Instance()->FindOrBuildMaterial("G4_Pb", isotopes);
292  // DiscoMaterialIORT = leadNist; // messenger
293  gray = new G4VisAttributes( G4Colour(0.5, 0.5, 0.5 ));
294  gray-> SetVisibility(true);
295  //gray -> SetForceWireframe(true);
296  //gray-> SetForceSolid(true);
297 
298  gray1 = new G4VisAttributes( G4Colour(0.7, 0.7, 0.7 ));
299  gray1-> SetVisibility(true);
300  //gray1 -> SetForceWireframe(true);
301  //gray1-> SetForceSolid(true);
302  // const G4double OuterRadiusDiscoIORT = 35. *CLHEP::mm; // messenger
303  // const G4double innerRadiusDiscoIORT = 0.*CLHEP::mm; // messenger
304  // const G4double heightDiscoIORT = 3.0*CLHEP::mm; // messenger
305  const G4double startAngleDiscoIORT = 0.*CLHEP::deg;
306  const G4double spanningAngleDiscoIORT = 360.*CLHEP::deg;
307  // const G4double DiscoXPositionIORT = -14.0*CLHEP::mm; // messenger
308 
309 //G4Material* DiscoMaterialIORT = G4NistManager::Instance()->FindOrBuildMaterial("G4_PLEXIGLASS", isotopes);// messenger
310 
311 
312  G4double phi = 0. *CLHEP::deg;
313 
314  // Matrix definition for a 90 deg rotation. Also used for other volumes
315  G4RotationMatrix rm;
316  rm.rotateY(phi);
317 
318 
319  solidDiscoIORT = new G4Tubs("DiscoIORT", innerRadiusDiscoIORT,
320  OuterRadiusDiscoIORT,
321  heightDiscoIORT,
322  startAngleDiscoIORT,
323  spanningAngleDiscoIORT);
324 
325  G4LogicalVolume* logDiscoIORT = new G4LogicalVolume(solidDiscoIORT,
326  DiscoMaterialIORT, "DiscoIORTLog", 0, 0, 0);
327 
328  physiDiscoIORT = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(0.,0.,(- heightDiscoIORT1))),
329  "DiscoIORTPhys", logDiscoIORT, physiDiscoIORT0, false, 0);
330 
331  logDiscoIORT -> SetVisAttributes(gray1);
332 
333 
334  // ---------------------------------------------------------------//
335  // 2.0 mm Lead Protection Disc //
336  // ---------------------------------------------------------------//
337 
338  // const G4double OuterRadiusDiscoIORT1 = 35. *CLHEP::mm;
339  // const G4double innerRadiusDiscoIORT1 = 0.*CLHEP::mm;
340  // const G4double heightDiscoIORT1 = 0.5*CLHEP::mm;
341  const G4double startAngleDiscoIORT1 = 0.*CLHEP::deg;
342  const G4double spanningAngleDiscoIORT1 = 360.*CLHEP::deg;
343 // const G4double DiscoXPositionIORT1 = -10.5*CLHEP::mm; messenger
344 // G4Material* DiscoMaterialIORT1 = G4NistManager::Instance()->FindOrBuildMaterial("G4_Cu", isotopes);// messenger
345 
346 
347 
348  solidDiscoIORT1 = new G4Tubs("DiscoIORT1", innerRadiusDiscoIORT1,
349  OuterRadiusDiscoIORT1,
350  heightDiscoIORT1,
351  startAngleDiscoIORT1,
352  spanningAngleDiscoIORT1);
353 
354  G4LogicalVolume* logDiscoIORT1 = new G4LogicalVolume(solidDiscoIORT1,
355  DiscoMaterialIORT1, "DiscoIORTLog1", 0, 0, 0);
356 
357  physiDiscoIORT1 = new G4PVPlacement(G4Transform3D(rm, G4ThreeVector(0.,0.,heightDiscoIORT)),
358  "DiscoIORTPhys1", logDiscoIORT1, physiDiscoIORT0, false, 0);
359  white = new G4VisAttributes( G4Colour());
360  white -> SetVisibility(true);
361  white -> SetForceSolid(true);
362  logDiscoIORT1 -> SetVisAttributes(gray);
363 
364 }
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4Tubs
Definition: G4Tubs.hh:84
CLHEP::HepRotation::rotateY
HepRotation & rotateY(double delta)
Definition: Rotation.cc:79
G4Transform3D
HepGeom::Transform3D G4Transform3D
Definition: G4Transform3D.hh:35
G4double
double G4double
Definition: G4Types.hh:76
void IORTDetectorConstruction::DeleteDisc ( )

Definition at line 705 of file IORTDetectorConstruction.cc.

References G4RunManager::GetRunManager().

706 {
707  delete solidDiscoIORT0;
708  delete logicDiscoIORT0;
709  delete physiDiscoIORT0;
710  delete solidDiscoIORT;
711  delete logicDiscoIORT;
712  delete physiDiscoIORT;
713  delete solidDiscoIORT1;
714  delete logicDiscoIORT1;
715  delete physiDiscoIORT1;
716  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
717  G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
718 }
G4RunManager::GetRunManager
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:74
G4LogicalVolume* IORTDetectorConstruction::GetDetectorLogicalVolume ( )
inline

Definition at line 169 of file IORTDetectorConstruction.hh.

169 { return detectorLogicalVolume;}
G4ThreeVector IORTDetectorConstruction::GetDetectorToPhantomPosition ( )
inline

Definition at line 84 of file IORTDetectorConstruction.hh.

References CLHEP::Hep3Vector::getX(), CLHEP::Hep3Vector::getY(), and CLHEP::Hep3Vector::getZ().

85 {
86  return G4ThreeVector(phantomSizeX/2 - detectorSizeX/2 + detectorPosition.getX(),
87  phantomSizeY/2 - detectorSizeY/2 + detectorPosition.getY(),
88  phantomSizeZ/2 - detectorSizeZ/2 + detectorPosition.getZ()
89  );
90 }
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
CLHEP::Hep3Vector::getY
double getY() const
CLHEP::Hep3Vector::getX
double getX() const
CLHEP::Hep3Vector::getZ
double getZ() const
G4ThreeVector IORTDetectorConstruction::GetDetectorToWorldPosition ( )
inline

Definition at line 78 of file IORTDetectorConstruction.hh.

Referenced by UpdateGeometry().

79  {
80  return phantomPosition + detectorPosition;
81  }
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.

References G4cout, G4endl, CLHEP::Hep3Vector::getX(), CLHEP::Hep3Vector::getY(), and CLHEP::Hep3Vector::getZ().

113 {
114 // Dimensions check... X Y and Z
115 // Firstly check what dimension we are modifying
116  {
117  if (detectorX > phantomX)
118  {
119  G4cout << "Error: Detector X dimension must be smaller or equal to the corrispondent of the phantom" << G4endl;
120  return false;
121  }
122  if ( (phantomX - detectorX) < detToPhantomPosition.getX())
123  {
124  G4cout << "Error: X dimension doesn't fit with detector to phantom relative position" << G4endl;
125  return false;
126  }
127  }
128 
129  {
130  if (detectorY > phantomY)
131  {
132  G4cout << "Error: Detector Y dimension must be smaller or equal to the corrispondent of the phantom" << G4endl;
133  return false;
134  }
135  if ( (phantomY - detectorY) < detToPhantomPosition.getY())
136  {
137  G4cout << "Error: Y dimension doesn't fit with detector to phantom relative position" << G4endl;
138  return false;
139  }
140  }
141 
142  {
143  if (detectorZ > phantomZ)
144  {
145  G4cout << "Error: Detector Z dimension must be smaller or equal to the corrispondent of the phantom" << G4endl;
146  return false;
147  }
148  if ( (phantomZ - detectorZ) < detToPhantomPosition.getZ())
149  {
150  G4cout << "Error: Z dimension doesn't fit with detector to phantom relative position" << G4endl;
151  return false;
152  }
153  }
154 
155  return true;
156 }
CLHEP::Hep3Vector::getY
double getY() const
CLHEP::Hep3Vector::getX
double getX() const
G4cout
G4GLOB_DLL std::ostream G4cout
CLHEP::Hep3Vector::getZ
double getZ() const
G4endl
#define G4endl
Definition: G4ios.hh:61
void IORTDetectorConstruction::PrintParameters ( )

Definition at line 721 of file IORTDetectorConstruction.cc.

References G4BestUnit, G4cout, G4endl, CLHEP::Hep3Vector::getX(), CLHEP::Hep3Vector::getY(), and CLHEP::Hep3Vector::getZ().

Referenced by UpdateGeometry().

722 {
723 
724  G4cout << "The (X,Y,Z) dimensions of the phantom are : (" <<
725  G4BestUnit( phantom -> GetXHalfLength()*2., "Length") << ',' <<
726  G4BestUnit( phantom -> GetYHalfLength()*2., "Length") << ',' <<
727  G4BestUnit( phantom -> GetZHalfLength()*2., "Length") << ')' << G4endl;
728 
729  G4cout << "The (X,Y,Z) dimensions of the detector are : (" <<
730  G4BestUnit( detector -> GetXHalfLength()*2., "Length") << ',' <<
731  G4BestUnit( detector -> GetYHalfLength()*2., "Length") << ',' <<
732  G4BestUnit( detector -> GetZHalfLength()*2., "Length") << ')' << G4endl;
733 
734  G4cout << "Displacement between Phantom and World is: ";
735  G4cout << "DX= "<< G4BestUnit(phantomPosition.getX(),"Length") <<
736  "DY= "<< G4BestUnit(phantomPosition.getY(),"Length") <<
737  "DZ= "<< G4BestUnit(phantomPosition.getZ(),"Length") << G4endl;
738 
739  G4cout << "The (X,Y,Z) sizes of the Voxels are: (" <<
740  G4BestUnit(sizeOfVoxelAlongX, "Length") << ',' <<
741  G4BestUnit(sizeOfVoxelAlongY, "Length") << ',' <<
742  G4BestUnit(sizeOfVoxelAlongZ, "Length") << ')' << G4endl;
743 
744  G4cout << "The number of Voxels along (X,Y,Z) is: (" <<
745  numberOfVoxelsAlongX << ',' <<
746  numberOfVoxelsAlongY <<',' <<
747  numberOfVoxelsAlongZ << ')' << G4endl;
748 
749 }
CLHEP::Hep3Vector::getY
double getY() const
G4BestUnit
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
Definition: G4SteppingVerbose.cc:53
CLHEP::Hep3Vector::getX
double getX() const
G4cout
G4GLOB_DLL std::ostream G4cout
CLHEP::Hep3Vector::getZ
double getZ() const
G4endl
#define G4endl
Definition: G4ios.hh:61
void IORTDetectorConstruction::SetAngleDiscoIORT0 ( G4double  phi0)

Definition at line 618 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

619 {
620 
621  AngleDiscoIORT0 = phi0;
622 }
void IORTDetectorConstruction::SetDetectorPosition ( )
inline

Definition at line 94 of file IORTDetectorConstruction.hh.

References CLHEP::Hep3Vector::getX(), CLHEP::Hep3Vector::getY(), CLHEP::Hep3Vector::getZ(), CLHEP::Hep3Vector::setX(), CLHEP::Hep3Vector::setY(), and CLHEP::Hep3Vector::setZ().

Referenced by UpdateGeometry().

95  {
96  // Adjust detector position
97  detectorPosition.setX(detectorToPhantomPosition.getX() - phantomSizeX/2 + detectorSizeX/2);
98  detectorPosition.setY(detectorToPhantomPosition.getY() - phantomSizeY/2 + detectorSizeY/2);
99  detectorPosition.setZ(detectorToPhantomPosition.getZ() - phantomSizeZ/2 + detectorSizeZ/2);
100 
101  //G4cout << "*************** DetectorToPhantomPosition " << detectorToPhantomPosition/cm << "\n";
102  //G4cout << "*************** DetectorPosition " << detectorPosition/cm << "\n";
103  }
CLHEP::Hep3Vector::getY
double getY() const
CLHEP::Hep3Vector::setY
void setY(double)
CLHEP::Hep3Vector::setZ
void setZ(double)
CLHEP::Hep3Vector::setX
void setX(double)
CLHEP::Hep3Vector::getX
double getX() const
CLHEP::Hep3Vector::getZ
double getZ() const
void IORTDetectorConstruction::SetDetectorSize ( G4double  sizeX,
G4double  sizeY,
G4double  sizeZ 
)

Definition at line 528 of file IORTDetectorConstruction.cc.

References SetVoxelSize().

Referenced by IORTDetectorConstruction().

529 {
530  if (sizeX > 0.) {detectorSizeX = sizeX;}
531  if (sizeY > 0.) {detectorSizeY = sizeY;}
532  if (sizeZ > 0.) {detectorSizeZ = sizeZ;}
533  SetVoxelSize(sizeOfVoxelAlongX, sizeOfVoxelAlongY, sizeOfVoxelAlongZ);
534 }
IORTDetectorConstruction::SetVoxelSize
void SetVoxelSize(G4double sizeX, G4double sizeY, G4double sizeZ)
Definition: IORTDetectorConstruction.cc:537
void IORTDetectorConstruction::SetDetectorToPhantomPosition ( G4ThreeVector  DetectorToPhantomPosition)

Definition at line 549 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

550 {
551  detectorToPhantomPosition = displ;
552 }
G4bool IORTDetectorConstruction::SetDiscoMaterialIORT ( G4String  material)

Definition at line 466 of file IORTDetectorConstruction.cc.

References G4cout, G4endl, G4RunManager::GetRunManager(), G4NistManager::Instance(), and EmPlot::SetMaterial().

Referenced by IORTDetectorConstruction().

467 {
468 
469  if (G4Material* dMat = G4NistManager::Instance()->FindOrBuildMaterial(material, false) )
470  {
471  DiscoMaterialIORT = dMat;
472 
473  if (logicDiscoIORT)
474  {
475  logicDiscoIORT -> SetMaterial(dMat);
476 
477  G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
478  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
479  G4cout << "The material of Protection disc 1 has been changed to " << material << G4endl;
480  }
481  }
482  else
483  {
484  G4cout << "WARNING: material \"" << material << "\" doesn't exist in NIST elements/materials"
485  " table [located in $G4INSTALL/source/materials/src/G4NistMaterialBuilder.cc]" << G4endl;
486  G4cout << "Use command \"/parameter/nist\" to see full materials list!" << G4endl;
487  return false;
488  }
489 
490  return true;
491 }
G4NistManager::Instance
static G4NistManager * Instance()
Definition: G4NistManager.cc:68
G4cout
G4GLOB_DLL std::ostream G4cout
EmPlot.SetMaterial
def SetMaterial
Definition: EmPlot.py:25
G4RunManager::GetRunManager
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:74
G4endl
#define G4endl
Definition: G4ios.hh:61
G4bool IORTDetectorConstruction::SetDiscoMaterialIORT1 ( G4String  material)

Definition at line 493 of file IORTDetectorConstruction.cc.

References G4cout, G4endl, G4RunManager::GetRunManager(), G4NistManager::Instance(), and EmPlot::SetMaterial().

Referenced by IORTDetectorConstruction().

494 {
495 
496  if (G4Material* d1Mat = G4NistManager::Instance()->FindOrBuildMaterial(material, false) )
497  {
498  DiscoMaterialIORT1 = d1Mat;
499 
500  if (logicDiscoIORT1)
501  {
502  logicDiscoIORT1 -> SetMaterial(d1Mat);
503 
504  G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
505  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
506  G4cout << "The material of Protection disc 2 has been changed to " << material << G4endl;
507  }
508  }
509  else
510  {
511  G4cout << "WARNING: material \"" << material << "\" doesn't exist in NIST elements/materials"
512  " table [located in $G4INSTALL/source/materials/src/G4NistMaterialBuilder.cc]" << G4endl;
513  G4cout << "Use command \"/parameter/nist\" to see full materials list!" << G4endl;
514  return false;
515  }
516 
517  return true;
518 }
G4NistManager::Instance
static G4NistManager * Instance()
Definition: G4NistManager.cc:68
G4cout
G4GLOB_DLL std::ostream G4cout
EmPlot.SetMaterial
def SetMaterial
Definition: EmPlot.py:25
G4RunManager::GetRunManager
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:74
G4endl
#define G4endl
Definition: G4ios.hh:61
void IORTDetectorConstruction::SetDiscoXPositionIORT ( G4double  xpos)

Definition at line 573 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

574 {
575 
576  DiscoXPositionIORT = xpos;
577 
578 }
void IORTDetectorConstruction::SetDiscoXPositionIORT1 ( G4double  xpos)

Definition at line 612 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

613 {
614 
615  DiscoXPositionIORT1 = xpos;
616 }
void IORTDetectorConstruction::SetDiscoYPositionIORT ( G4double  ypos)

Definition at line 580 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

581 {
582 
583  DiscoYPositionIORT = ypos;
584 
585 }
void IORTDetectorConstruction::SetDiscoZPositionIORT ( G4double  zpos)

Definition at line 587 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

588 {
589 
590  DiscoZPositionIORT = zpos;
591 
592 }
void IORTDetectorConstruction::SetheightDiscoIORT ( G4double  height)

Definition at line 567 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

568 {
569  if (height > 0.) {heightDiscoIORT = height;}
570 
571 }
void IORTDetectorConstruction::SetheightDiscoIORT1 ( G4double  height)

Definition at line 606 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

607 {
608  if (height > 0.) {heightDiscoIORT1 = height;}
609 
610 }
void IORTDetectorConstruction::SetinnerRadiusDiscoIORT ( G4double  innerr)

Definition at line 561 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

562 {
563  if (innerr > 0.) {innerRadiusDiscoIORT = innerr;}
564 
565 }
void IORTDetectorConstruction::SetinnerRadiusDiscoIORT1 ( G4double  innerr)

Definition at line 600 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

601 {
602  if (innerr > 0.) {innerRadiusDiscoIORT1 = innerr;}
603 
604 }
void IORTDetectorConstruction::SetOuterRadiusDiscoIORT ( G4double  outerr)

Definition at line 555 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

556 {
557  if (outerr > 0.) {OuterRadiusDiscoIORT = outerr;}
558 
559 }
void IORTDetectorConstruction::SetOuterRadiusDiscoIORT1 ( G4double  outerr)

Definition at line 594 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

595 {
596  if (outerr > 0.) {OuterRadiusDiscoIORT1 = outerr;}
597 
598 }
G4bool IORTDetectorConstruction::SetPhantomMaterial ( G4String  material)

Definition at line 438 of file IORTDetectorConstruction.cc.

References G4cout, G4endl, G4RunManager::GetRunManager(), G4NistManager::Instance(), and EmPlot::SetMaterial().

Referenced by IORTDetectorConstruction().

439 {
440 
441  if (G4Material* pMat = G4NistManager::Instance()->FindOrBuildMaterial(material, false) )
442  {
443  phantomMaterial = pMat;
444  detectorMaterial = pMat;
445  if (detectorLogicalVolume && phantomLogicalVolume)
446  {
447  detectorLogicalVolume -> SetMaterial(pMat);
448  phantomLogicalVolume -> SetMaterial(pMat);
449 
450  G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
451  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
452  G4cout << "The material of Phantom/Detector has been changed to " << material << G4endl;
453  }
454  }
455  else
456  {
457  G4cout << "WARNING: material \"" << material << "\" doesn't exist in NIST elements/materials"
458  " table [located in $G4INSTALL/source/materials/src/G4NistMaterialBuilder.cc]" << G4endl;
459  G4cout << "Use command \"/parameter/nist\" to see full materials list!" << G4endl;
460  return false;
461  }
462 
463  return true;
464 }
G4NistManager::Instance
static G4NistManager * Instance()
Definition: G4NistManager.cc:68
G4cout
G4GLOB_DLL std::ostream G4cout
EmPlot.SetMaterial
def SetMaterial
Definition: EmPlot.py:25
G4RunManager::GetRunManager
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:74
G4endl
#define G4endl
Definition: G4ios.hh:61
void IORTDetectorConstruction::SetPhantomPosition ( G4ThreeVector  pos)

Definition at line 543 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

544 {
545  phantomPosition = pos;
546 }
void IORTDetectorConstruction::SetPhantomSize ( G4double  sizeX,
G4double  sizeY,
G4double  sizeZ 
)

Definition at line 520 of file IORTDetectorConstruction.cc.

Referenced by IORTDetectorConstruction().

521 {
522  if (sizeX > 0.) phantomSizeX = sizeX;
523  if (sizeY > 0.) phantomSizeY = sizeY;
524  if (sizeZ > 0.) phantomSizeZ = sizeZ;
525 }
void IORTDetectorConstruction::SetVoxelSize ( G4double  sizeX,
G4double  sizeY,
G4double  sizeZ 
)

Definition at line 537 of file IORTDetectorConstruction.cc.

Referenced by SetDetectorSize().

538 {
539  if (sizeX > 0.) {sizeOfVoxelAlongX = sizeX;}
540  if (sizeY > 0.) {sizeOfVoxelAlongY = sizeY;}
541  if (sizeZ > 0.) {sizeOfVoxelAlongZ = sizeZ;}
542 }
void IORTDetectorConstruction::UpdateGeometry ( )

Definition at line 628 of file IORTDetectorConstruction.cc.

References ConstructDisc(), G4lrint(), GetDetectorToWorldPosition(), G4GeometryManager::GetInstance(), IORTMatrix::GetInstance(), G4RunManager::GetRunManager(), PrintParameters(), and SetDetectorPosition().

Referenced by IORTDetectorConstruction().

629 {
630  /*
631  * Check parameters consistency
632  */
633  ParametersCheck();
634 
635  G4GeometryManager::GetInstance() -> OpenGeometry();
636  if (phantom)
637  {
638  phantom -> SetXHalfLength(phantomSizeX/2);
639  phantom -> SetYHalfLength(phantomSizeY/2);
640  phantom -> SetZHalfLength(phantomSizeZ/2);
641  phantomPhysicalVolume -> SetTranslation(phantomPosition);
642  }
643  else ConstructPhantom();
644 
645  // Get the center of the detector
646  SetDetectorPosition();
647  if (detector)
648  {
649  detector -> SetXHalfLength(detectorSizeX/2);
650  detector -> SetYHalfLength(detectorSizeY/2);
651  detector -> SetZHalfLength(detectorSizeZ/2);
652  detectorPhysicalVolume -> SetTranslation(detectorPosition);
653  }
654  else ConstructDetector();
655 
656  // update disc function
657  delete solidDiscoIORT0;
658  delete logicDiscoIORT0;
659  delete physiDiscoIORT0;
660  delete solidDiscoIORT;
661  delete logicDiscoIORT;
662  delete physiDiscoIORT;
663  delete solidDiscoIORT1;
664  delete logicDiscoIORT1;
665  delete physiDiscoIORT1;
666  ConstructDisc();
667 
668 
669  // Round to nearest integer number of voxel
670  numberOfVoxelsAlongX = G4lrint(detectorSizeX / sizeOfVoxelAlongX);
671  sizeOfVoxelAlongX = ( detectorSizeX / numberOfVoxelsAlongX );
672 
673  numberOfVoxelsAlongY = G4lrint(detectorSizeY / sizeOfVoxelAlongY);
674  sizeOfVoxelAlongY = ( detectorSizeY / numberOfVoxelsAlongY );
675 
676  numberOfVoxelsAlongZ = G4lrint(detectorSizeZ / sizeOfVoxelAlongZ);
677  sizeOfVoxelAlongZ = ( detectorSizeZ / numberOfVoxelsAlongZ );
678 
679  //G4cout << "*************** DetectorToWorldPosition " << GetDetectorToWorldPosition()/cm << "\n";
680  ConstructSensitiveDetector(GetDetectorToWorldPosition());
681 
682  volumeOfVoxel = sizeOfVoxelAlongX * sizeOfVoxelAlongY * sizeOfVoxelAlongZ;
683  massOfVoxel = detectorMaterial -> GetDensity() * volumeOfVoxel;
684  // This will clear the existing matrix (together with all data inside it)!
685  matrix = IORTMatrix::GetInstance(numberOfVoxelsAlongX,
686  numberOfVoxelsAlongY,
687  numberOfVoxelsAlongZ,
688  massOfVoxel);
689 
690  // Initialize analysis
691 /*
692  IORTAnalysisManager* analysis = IORTAnalysisManager::GetInstance();
693 #ifdef G4ANALYSIS_USE_ROOT
694  analysis -> flush(); // Finalize the root file
695  analysis -> book();
696 #endif
697 */
698  // Inform the kernel about the new geometry
699  G4RunManager::GetRunManager() -> GeometryHasBeenModified();
700  G4RunManager::GetRunManager() -> PhysicsHasBeenModified();
701 
702  PrintParameters();
703 }
G4GeometryManager::GetInstance
static G4GeometryManager * GetInstance()
Definition: G4GeometryManager.cc:135
G4RunManager::GetRunManager
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:74
G4lrint
int G4lrint(double ad)
Definition: templates.hh:163
IORTMatrix::GetInstance
static IORTMatrix * GetInstance()
Definition: IORTMatrix.cc:61
The documentation for this class was generated from the following files:
Generated on Wed Apr 30 2014 15:57:25 for Geant4.10 by   doxygen 1.8.7