| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // << 26 // Authors: S. Guatelli and M. G. Pia, INFN Genova, Italy
27 // Previous authors: G. Guerrieri, S. Guatelli <<
28 // Authors (since 2007): S. Guatelli, Universi <<
29 // 27 //
>> 28 // Based on code developed by the undergraduate student G. Guerrieri
>> 29 // Note: this is a preliminary beta-version of the code; an improved
>> 30 // version will be distributed in the next Geant4 public release, compliant
>> 31 // with the design in a forthcoming publication, and subject to a
>> 32 // design and code review.
30 // 33 //
31 #include "G4MIRDRibCage.hh" 34 #include "G4MIRDRibCage.hh"
32 35
>> 36 #include "G4Processor/GDMLProcessor.h"
33 #include "globals.hh" 37 #include "globals.hh"
34 #include "G4SystemOfUnits.hh" <<
35 #include "G4SDManager.hh" 38 #include "G4SDManager.hh"
36 #include "G4VisAttributes.hh" 39 #include "G4VisAttributes.hh"
37 #include "G4HumanPhantomMaterial.hh" <<
38 #include "G4SubtractionSolid.hh" <<
39 #include "G4EllipticalTube.hh" <<
40 #include "G4PVReplica.hh" <<
41 #include "G4Box.hh" <<
42 #include "G4PVPlacement.hh" <<
43 #include "G4HumanPhantomColour.hh" <<
44 G4MIRDRibCage::G4MIRDRibCage():fPhysRib1(nullp <<
45 fPhysRib4(nullptr), fPhysRib5(nu <<
46 fPhysRib7(nullptr), fPhysRib8(nu <<
47 fPhysRib10(nullptr), fPhysRib11( <<
48 {;} <<
49 <<
50 G4VPhysicalVolume* G4MIRDRibCage::Construct(co <<
51 const G4String& colourName, G4bo <<
52 <<
53 { <<
54 auto* material = new G4HumanPhantomMaterial( <<
55 <<
56 G4cout<<"Construct "<<volumeName<<" with mot <<
57 <<
58 auto* skeleton = material -> GetMaterial("sk <<
59 auto* soft = material -> GetMaterial("soft_t <<
60 <<
61 delete material; <<
62 <<
63 G4double dx= 17. *cm; // a2 <<
64 G4double dy= 9.8 * cm; //b2 <<
65 G4double thickness= 32.4 * cm; // z2/2 of ca <<
66 40
67 auto* outCage = new G4EllipticalTube("outCag << 41 G4MIRDRibCage::G4MIRDRibCage()
68 << 42 {
69 dx = 16.4 * cm; // a1 << 43 }
70 dy = 9.2 * cm; // b1 <<
71 G4double dz = 34. *cm; // z2/2 <<
72 44
73 auto* inCage = new G4EllipticalTube("inCage" << 45 G4MIRDRibCage::~G4MIRDRibCage()
>> 46 {
>> 47 sxp.Finalize();
>> 48 }
74 49
75 auto* cage = new G4SubtractionSolid("Cage", << 50 G4VPhysicalVolume* G4MIRDRibCage::ConstructRibCage(G4VPhysicalVolume* mother, G4String sex, G4bool sensitivity)
76 outCage, << 51 {
77 inCage, nullptr, G4ThreeVector << 52 // Initialize GDML Processor
>> 53 sxp.Initialize();
>> 54 config.SetURI( "gdmlData/"+sex+"/MIRDRibCage.gdml" );
>> 55 config.SetSetupName( "Default" );
>> 56 sxp.Configure( &config );
78 57
>> 58 // Run GDML Processor
>> 59 sxp.Run();
79 60
80 auto* logicRibCage = new G4LogicalVolume(cag <<
81 61
82 G4VPhysicalVolume* physRibCage = new G4PVPla << 62 G4LogicalVolume* logicRibCage = (G4LogicalVolume *)GDMLProcessor::GetInstance()->GetLogicalVolume("RibCageVolume");
83 // with respect to the trunk <<
84 "physicalRibCage", <<
85 logicRibCage, <<
86 mother, <<
87 false, <<
88 0, true); <<
89 <<
90 G4double xx = 17.*cm; <<
91 G4double yy = 9.8*cm; <<
92 G4double ribThickness = 1.4*cm; <<
93 auto* rib_out = new G4EllipticalTube("rib_ou <<
94 <<
95 xx = 16.5 *cm; <<
96 yy = 9.3 * cm; <<
97 G4double zz = 1.5 * cm; <<
98 auto* rib_in = new G4EllipticalTube("rib_in" <<
99 auto* rib = new G4SubtractionSolid("rib",rib <<
100 <<
101 auto* logicRib= new G4LogicalVolume(rib, ske <<
102 <<
103 fPhysRib1 = new G4PVPlacement(nullptr,G4Thre <<
104 // with respect to the trunk <<
105 "physicalRib", <<
106 logicRib, <<
107 physRibCage, <<
108 false, <<
109 0, true); <<
110 <<
111 fPhysRib2 = new G4PVPlacement(nullptr,G4Thre <<
112 // with respect to the trunk <<
113 "physicalRib", <<
114 logicRib, <<
115 physRibCage, <<
116 false, <<
117 0, true); <<
118 <<
119 fPhysRib3 = new G4PVPlacement(nullptr,G4Thre <<
120 // with respect to the trunk <<
121 "physicalRib", <<
122 logicRib, <<
123 physRibCage, <<
124 false, <<
125 0, true); <<
126 <<
127 fPhysRib4 = new G4PVPlacement(nullptr,G4Thre <<
128 // with respect to the trunk <<
129 "physicalRib", <<
130 logicRib, <<
131 physRibCage, <<
132 false, <<
133 0, true); <<
134 63
135 fPhysRib5 = new G4PVPlacement(nullptr,G4Thre << 64 G4ThreeVector position = (G4ThreeVector)*GDMLProcessor::GetInstance()->GetPosition("RibCagePos");
136 // with respect to the trunk << 65 G4RotationMatrix* rm = (G4RotationMatrix*)GDMLProcessor::GetInstance()->GetRotation("RibCageRot");
137 "physicalRib", <<
138 logicRib, <<
139 physRibCage, <<
140 false, <<
141 0, true); <<
142 <<
143 fPhysRib6 = new G4PVPlacement(nullptr,G4Thre <<
144 // with respect to the trunk <<
145 "physicalRib", <<
146 logicRib, <<
147 physRibCage, <<
148 false, <<
149 0, true); <<
150 <<
151 fPhysRib7 = new G4PVPlacement(nullptr,G4Thre <<
152 // with respect to the trunk <<
153 "physicalRib", <<
154 logicRib, <<
155 physRibCage, <<
156 false, <<
157 0, true); <<
158 <<
159 fPhysRib8 = new G4PVPlacement(nullptr,G4Thre <<
160 // with respect to the trunk <<
161 "physicalRib", <<
162 logicRib, <<
163 physRibCage, <<
164 false, <<
165 0, true); <<
166 <<
167 fPhysRib9 = new G4PVPlacement(nullptr,G4Thre <<
168 // with respect to the trunk <<
169 "physicalRib", <<
170 logicRib, <<
171 physRibCage, <<
172 false, <<
173 0, true); <<
174 <<
175 fPhysRib10 = new G4PVPlacement(nullptr, G4Th <<
176 // with respect to the trunk <<
177 "physicalRib", <<
178 logicRib, <<
179 physRibCage, <<
180 false, <<
181 0, true); <<
182 <<
183 fPhysRib11 = new G4PVPlacement(nullptr,G4Thr <<
184 // with respect to the trunk <<
185 "physicalRib", <<
186 logicRib, <<
187 physRibCage, <<
188 false, <<
189 0, true); <<
190 <<
191 fPhysRib12 = new G4PVPlacement(nullptr,G4Thr <<
192 // with respect to the trunk <<
193 "physicalRib", <<
194 logicRib, <<
195 physRibCage, <<
196 false, <<
197 0, true); <<
198 66
199 // Visualization Attributes << 67 // Define rotation and position here!
200 logicRibCage -> SetVisAttributes(G4VisAttrib << 68 G4VPhysicalVolume* physRibCage = new G4PVPlacement(rm,position,
201 << 69 "physicalRibCage",
202 //G4VisAttributes* RibCageVisAtt = new G4Vis << 70 logicRibCage,
>> 71 mother,
>> 72 false,
>> 73 0);
>> 74
>> 75 // Sensitive Body Part
>> 76 if (sensitivity==true)
>> 77 {
>> 78 G4SDManager* SDman = G4SDManager::GetSDMpointer();
>> 79 logicRibCage->SetSensitiveDetector( SDman->FindSensitiveDetector("BodyPartSD") );
>> 80 }
203 81
204 auto* colourPointer = new G4HumanPhantomColo << 82 // Visualization Attributes
205 G4Colour colour = colourPointer -> GetColour << 83 G4VisAttributes* RibCageVisAtt = new G4VisAttributes(G4Colour(1.0,1.0,0.0));
206 auto* RibCageVisAtt = new G4VisAttributes(co << 84 RibCageVisAtt->SetForceSolid(true);
207 RibCageVisAtt->SetForceSolid(wireFrame); << 85 logicRibCage->SetVisAttributes(RibCageVisAtt);
208 logicRib->SetVisAttributes(RibCageVisAtt); <<
209 86
210 G4cout << "RibCage created !!!!!!" << G4endl 87 G4cout << "RibCage created !!!!!!" << G4endl;
211 // Testing Pelvis Volume <<
212 G4double RibCageVol = logicRib->GetSolid()-> <<
213 G4cout << "Volume of RibCage = " << ((RibCag <<
214 <<
215 // Testing RibCage Material <<
216 G4String RibCageMat = logicRib->GetMaterial( <<
217 G4cout << "Material of RibCage = " << RibCag <<
218 88
219 // Testing Density <<
220 G4double RibCageDensity = logicRib->GetMater <<
221 G4cout << "Density of Material = " << RibCag <<
222 <<
223 // Testing Mass <<
224 G4double RibCageMass = (RibCageVol)* RibCage <<
225 G4cout << "Mass of RibCage = " << (RibCageMa <<
226 <<
227 return physRibCage; 89 return physRibCage;
228 } 90 }
229 91