Geant4  10.02.p02
BrachyDetectorConstructionLeipzig.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 // --------------------------------------------------------------
28 // GEANT 4 - Brachytherapy example
29 // --------------------------------------------------------------
30 //
31 // Code developed by:
32 // S.Guatelli
33 //
34 // *******************************************
35 // * *
36 // * BrachyDetectorConstructionLeipzig.cc *
37 // * *
38 // *******************************************
39 //
40 //
41 // $Id: BrachyDetectorConstructionLeipzig.cc 69765 2013-05-14 10:11:22Z gcosmo $
42 //
43 // Code by S. Guatelli
44 //
45 #include "globals.hh"
46 #include "G4SystemOfUnits.hh"
48 #include "G4CSGSolid.hh"
49 #include "G4Sphere.hh"
50 #include "G4RunManager.hh"
51 #include "G4Box.hh"
52 #include "G4Tubs.hh"
53 #include "G4LogicalVolume.hh"
54 #include "G4ThreeVector.hh"
55 #include "G4PVPlacement.hh"
56 #include "G4Transform3D.hh"
57 #include "G4RotationMatrix.hh"
59 #include "BrachyMaterial.hh"
60 #include "G4VisAttributes.hh"
61 #include "G4Colour.hh"
62 
63 // Leipzig Applicator ...
64 
66 capsule(0),capsuleTip(0), iridiumCore(0), applicator1(0), applicator2(0),
67 capsuleLog(0), capsuleTipLog(0), iridiumCoreLog(0), applicator1Log(0),
68 applicator2Log(0),capsulePhys(0), capsuleTipPhys(0),iridiumCorePhys(0),
69 applicator1Phys(0), applicator2Phys(0)
70 {
71  pMaterial = new BrachyMaterial();
72 }
73 
75 {
76  delete pMaterial;
77 }
78 
80 {
81  G4Colour red (1.0, 0.0, 0.0) ;
82  G4Colour lblue (0.0, 0.0, .75);
83 
84  G4Material* capsuleMat = pMaterial -> GetMat("Stainless steel");
85  G4Material* iridium = pMaterial -> GetMat("Iridium");
86  G4Material* tungsten = pMaterial -> GetMat("Tungsten");
87 
88  //Iridium source ...
89 
90  capsule = new G4Tubs("Capsule",0,0.55*mm,3.725*mm,0.*deg,360.*deg);
91  capsuleLog = new G4LogicalVolume(capsule,capsuleMat,"CapsuleLog");
92  capsulePhys = new G4PVPlacement(0, G4ThreeVector(0,0,-1.975*mm),"CapsulePhys",
93  capsuleLog,mother, //mother volume: phantom
94  false,0, true);
95 
96  // Capsule tip
97  capsuleTip = new G4Sphere("CapsuleTip",0.*mm,0.55*mm,0.*deg,360.*deg,0.*deg,90.*deg);
98  capsuleTipLog = new G4LogicalVolume(capsuleTip,capsuleMat,"CapsuleTipLog");
99  capsuleTipPhys = new G4PVPlacement(0,G4ThreeVector(0.,0.,1.75*mm),"CapsuleTipPhys",
100  capsuleTipLog,mother,false,0, true);
101  // Iridium core
102  iridiumCore = new G4Tubs("IrCore",0,0.30*mm,1.75*mm,0.*deg,360.*deg);
103  iridiumCoreLog = new G4LogicalVolume(iridiumCore, iridium, "IridiumCoreLog");
104  iridiumCorePhys = new G4PVPlacement(0,G4ThreeVector(0.,0.,1.975*mm),"IridiumCorePhys",
105  iridiumCoreLog,capsulePhys,false,0, true);
106 
107  //Leipzig Applicator is modelled with two different volumes
108  applicator1 = new G4Tubs("Appl1",5*mm,10.5*mm,12*mm,0.*deg,360.*deg);
109  applicator1Log = new G4LogicalVolume(applicator1,tungsten,"Appl1Log");
110  applicator1Phys = new G4PVPlacement(0,G4ThreeVector(0,0,4.0*mm),"Appl1Phys",applicator1Log,
111  mother,false,0, true);
112 
113  applicator2 = new G4Tubs("Appl2",0.55*mm,5.*mm,3.125*mm,0.*deg,360.*deg);
114  applicator2Log = new G4LogicalVolume(applicator2,tungsten,"Appl2");
115  applicator2Phys = new G4PVPlacement(0,G4ThreeVector(0,0,-4.875*mm),
116  "Appl2Phys",applicator2Log,mother,false,0, true);
117 
119  simpleCapsuleVisAtt -> SetVisibility(true);
120  simpleCapsuleVisAtt -> SetForceWireframe(true);
121  capsuleLog -> SetVisAttributes(simpleCapsuleVisAtt);
122 
124  simpleCapsuleTipVisAtt -> SetVisibility(true);
125  simpleCapsuleTipVisAtt -> SetForceSolid(true);
126  capsuleTipLog -> SetVisAttributes(simpleCapsuleTipVisAtt);
127  iridiumCoreLog -> SetVisAttributes(simpleCapsuleTipVisAtt);
128 
129  applicatorVisAtt = new G4VisAttributes(lblue);
130  applicatorVisAtt -> SetVisibility(true);
131  applicatorVisAtt -> SetForceWireframe(true);
132  applicator1Log -> SetVisAttributes(applicatorVisAtt);
133  applicator2Log -> SetVisAttributes(applicatorVisAtt);
134 }
136 {
140 delete applicator2Phys; applicator2Phys = 0;
141 delete applicator1Phys; applicator1Phys = 0;
142 delete iridiumCorePhys; iridiumCorePhys = 0;
143 delete capsuleTipPhys; capsuleTipPhys = 0;
144 delete capsulePhys; capsulePhys = 0;
145 delete applicator2Log; applicator2Log = 0;
146 delete applicator1Log; applicator1Log = 0;
147 delete iridiumCoreLog; iridiumCoreLog = 0;
148 delete capsuleTipLog; capsuleTipLog = 0;
149 delete capsuleLog; capsuleLog = 0;
150 }
Definition: test07.cc:36
CLHEP::Hep3Vector G4ThreeVector
Definition: G4Tubs.hh:85
static const double deg
Definition: G4SIunits.hh:151
static const double mm
Definition: G4SIunits.hh:114