Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/brachytherapy/src/BrachyDetectorConstructionI.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  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: S.Guatelli
 32 //
 33 ///    ****************************************
 34 //    *                                      *
 35 //    *    BrachyDetectorConstructionI.cc     *
 36 //    *                                      *
 37 //    ****************************************
 38 //
 39 //
 40 #include "BrachyDetectorConstructionI.hh"
 41 #include "globals.hh"
 42 #include "G4SystemOfUnits.hh"
 43 #include "G4CSGSolid.hh"
 44 #include "G4Sphere.hh"
 45 #include "G4MaterialPropertyVector.hh"
 46 #include "G4RunManager.hh"
 47 #include "G4Box.hh"
 48 #include "G4Tubs.hh"
 49 #include "G4LogicalVolume.hh"
 50 #include "G4ThreeVector.hh"
 51 #include "G4PVPlacement.hh"
 52 #include "G4Transform3D.hh"
 53 #include "G4RotationMatrix.hh"
 54 #include "G4TransportationManager.hh"
 55 #include "G4VisAttributes.hh"
 56 #include "G4Colour.hh"
 57 #include "G4NistManager.hh"
 58 
 59 BrachyDetectorConstructionI::BrachyDetectorConstructionI():
 60 fDefaultTub(nullptr), fCapsule(nullptr), fCapsuleTip(nullptr), fIodineCore(nullptr), fDefaultTubLog(nullptr),
 61 fCapsuleLog(nullptr), fCapsuleTipLog(nullptr), fIodineCoreLog(nullptr), fDefaultTubPhys(nullptr), 
 62 fCapsulePhys(nullptr),fCapsuleTipPhys1(nullptr),fCapsuleTipPhys2(nullptr), fIodineCorePhys(nullptr),
 63 fSimpleIodineVisAtt(nullptr), fSimpleCapsuleVisAtt(nullptr), fSimpleCapsuleTipVisAtt(nullptr)
 64 {}
 65 
 66 void BrachyDetectorConstructionI::ConstructIodine(G4VPhysicalVolume* mother)
 67 {
 68   // Model of the Bebig Isoseed I-125 brachy source
 69   G4NistManager* nist = G4NistManager::Instance();
 70   G4Material* titanium = nist -> FindOrBuildMaterial("G4_Ti");
 71   G4Material* iodine = nist -> FindOrBuildMaterial("G4_I");
 72   G4Material* air = nist -> FindOrBuildMaterial("G4_AIR");
 73   
 74   G4Colour  red     (1.0, 0.0, 0.0) ;
 75   G4Colour  magenta (1.0, 0.0, 1.0) ; 
 76   G4Colour  lblue   (0.0, 0.0, .75);
 77 
 78   // Air tub
 79   fDefaultTub = new G4Tubs("DefaultTub",0.*mm, 0.40*mm, 1.84*mm, 0.*deg, 360.*deg);
 80   fDefaultTubLog = new G4LogicalVolume(fDefaultTub,air,"DefaultTub_Log");
 81   fDefaultTubPhys = new G4PVPlacement(nullptr,
 82                                       G4ThreeVector(),
 83                                       "defaultTub_Phys",
 84                                       fDefaultTubLog,
 85                                       mother,
 86                                       false,
 87                                       0, true); 
 88   // Capsule main body ...
 89   fCapsule = new G4Tubs("fCapsule", 0.35*mm,0.40*mm,1.84*mm,0.*deg,360.*deg);
 90   fCapsuleLog = new G4LogicalVolume(fCapsule,titanium,"fCapsuleLog");
 91   fCapsulePhys = new G4PVPlacement(nullptr,
 92                                   G4ThreeVector(),
 93                                   "fCapsulePhys",
 94                                   fCapsuleLog,
 95                                   fDefaultTubPhys,
 96                                   false,
 97                                   0, true);
 98   // fCapsule tips
 99   fCapsuleTip = new G4Sphere("fCapsuleTip",
100                             0.*mm,
101                             0.40*mm,
102                             0.*deg,
103                             360.*deg,
104                             0.*deg,
105                             90.*deg);
106                             
107   fCapsuleTipLog = new G4LogicalVolume(fCapsuleTip,titanium,"fCapsuleTipLog");
108   fCapsuleTipPhys1 = new G4PVPlacement(nullptr,
109                                       G4ThreeVector(0.,0.,1.84*mm),
110                                       "IodinefCapsuleTipPhys1",
111                                       fCapsuleTipLog,
112                                       mother,
113                                       false,
114                                       0, true);
115 
116   auto rotateMatrix = new G4RotationMatrix();
117   rotateMatrix -> rotateX(180.0*deg);
118   fCapsuleTipPhys2 = new G4PVPlacement(rotateMatrix, 
119                                       G4ThreeVector(0,0,-1.84*mm),
120                                       "IodinefCapsuleTipPhys2",
121                                       fCapsuleTipLog,
122                                       mother,
123                                       false,
124                                       0, true);
125  
126   // Radiactive core ...
127   fIodineCore = new G4Tubs("ICore",0.085*mm,0.35*mm,1.75*mm,0.*deg,360.*deg);
128   fIodineCoreLog = new G4LogicalVolume(fIodineCore,iodine,"IodineCoreLog");
129   fIodineCorePhys = new G4PVPlacement(nullptr,
130                                      G4ThreeVector(0.,0.,0.),
131                                      "IodineCorePhys",
132                                      fIodineCoreLog,
133                                      fDefaultTubPhys,
134                                      false,
135                                      0, true);
136  
137   // Visual attributes ...
138   
139   fSimpleIodineVisAtt = new G4VisAttributes(magenta);
140   fSimpleIodineVisAtt -> SetVisibility(true);
141   fSimpleIodineVisAtt -> SetForceSolid(true);
142   fIodineCoreLog -> SetVisAttributes(fSimpleIodineVisAtt);
143 
144   fSimpleCapsuleVisAtt = new G4VisAttributes(red);
145   fSimpleCapsuleVisAtt -> SetVisibility(true);  
146   fSimpleCapsuleVisAtt -> SetForceWireframe(true);
147   fCapsuleLog -> SetVisAttributes( fSimpleCapsuleVisAtt);
148 
149   fSimpleCapsuleTipVisAtt = new G4VisAttributes(red);
150   fSimpleCapsuleTipVisAtt -> SetVisibility(true); 
151   fSimpleCapsuleTipVisAtt -> SetForceSolid(true);
152   fCapsuleTipLog -> SetVisAttributes(fSimpleCapsuleTipVisAtt);
153 }
154 
155 void BrachyDetectorConstructionI::CleanIodine()
156 {
157 delete fSimpleIodineVisAtt; fSimpleIodineVisAtt = nullptr;
158 delete fSimpleCapsuleVisAtt; fSimpleCapsuleVisAtt = nullptr;
159 delete fSimpleCapsuleTipVisAtt; fSimpleCapsuleTipVisAtt = nullptr;
160 delete fCapsuleTipPhys1; fCapsuleTipPhys1 = nullptr; 
161 delete fCapsuleTipPhys2; fCapsuleTipPhys2 = nullptr;
162 delete fIodineCorePhys; fIodineCorePhys = nullptr;
163 delete fCapsulePhys; fCapsulePhys = nullptr;
164 delete fDefaultTubPhys; fDefaultTubPhys = nullptr; 
165 delete fDefaultTubLog; fDefaultTubLog = nullptr;
166 delete fCapsuleLog; fCapsuleLog = nullptr;
167 delete fCapsuleTipLog; fCapsuleTipLog = nullptr;
168 delete fIodineCoreLog; fIodineCoreLog = nullptr;
169 delete fDefaultTub; fDefaultTub = nullptr;
170 delete fCapsule; fCapsule = nullptr;
171 delete fCapsuleTip; fCapsuleTip = nullptr;
172 delete fIodineCore; fIodineCore = nullptr;
173 
174 G4RunManager::GetRunManager() -> GeometryHasBeenModified();
175 }
176