| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // Previous authors: G. Guerrieri, S. Guatelli
27 // Authors (since 2007): S. Guatelli, Universi
28 //
29 //
30 //
31 #include "G4MIRDPelvis.hh"
32
33 #include "globals.hh"
34 #include "G4SystemOfUnits.hh"
35 #include "G4SDManager.hh"
36 #include "G4VisAttributes.hh"
37 #include "G4HumanPhantomMaterial.hh"
38 #include "G4EllipticalTube.hh"
39 #include "G4RotationMatrix.hh"
40 #include "G4ThreeVector.hh"
41 #include "G4VPhysicalVolume.hh"
42 #include "G4PVPlacement.hh"
43 #include "G4SubtractionSolid.hh"
44 #include "G4Box.hh"
45 #include "G4VSolid.hh"
46 #include "G4LogicalVolume.hh"
47 #include "G4HumanPhantomColour.hh"
48
49 G4VPhysicalVolume* G4MIRDPelvis::Construct(con
50 const G4String& colourName, G4boo
51 {
52 auto* material = new G4HumanPhantomMaterial(
53
54 G4cout<<"Construct "<<volumeName<<" with mot
55
56 auto* skeleton = material -> GetMaterial("sk
57
58 delete material;
59 /*
60 G4double dx= 10.35 * cm;//12. *cm; // a2
61 G4double dy= 11.76 * cm;//12. * cm; // b2
62 G4double dz= 9.915 * cm; // z2/2
63
64 G4VSolid* outPelvis = new G4EllipticalTube
65
66 G4double dx_in = 9.75 * cm;//11.3 * cm; //
67 G4double dy_in = 11.07 *cm; //11.3* cm; //
68 G4double dz_in = 10. * cm;//11.0 *cm; // z
69
70 */
71 G4double dx= 12. *cm; // a2
72 G4double dy= 12. * cm; //b2
73 G4double dz= 11. * cm; // z2/2
74
75 G4VSolid* outPelvis = new G4EllipticalTube("
76
77 dx = 11.3 * cm; // a1
78 dy = 11.3* cm; // b1
79 dz = 12.0 *cm; // z2/2
80
81 G4VSolid* inPelvis = new G4EllipticalTube("I
82
83 G4double x = 28. * cm; // a2 * 2
84 G4double y = 28. * cm; //b2*2
85 G4double z = 24. *cm; // z2
86
87 auto* subPelvis = new G4Box("SubtrPelvis", x
88
89 auto* firstPelvis = new G4SubtractionSolid("
90 outPelvis,
91 inPelvis, nullptr, G4ThreeVec
92
93 auto* secondPelvis = new G4SubtractionSolid(
94 firstPelvis,
95 subPelvis, nullptr,
96 G4ThreeVector(0.0,
97 -14. * cm, 0.*cm));
98
99 auto* pelvis = new G4SubtractionSolid("Pelvi
100 nullptr,
101 G4ThreeVector(0.0,
102 22. * cm, -9. *cm));
103
104
105 auto* logicPelvis = new G4LogicalVolume(pelv
106 "logical" + volumeName, nullp
107
108
109 G4VPhysicalVolume* physPelvis = new G4PVPlac
110 // with respect to the trunk
111 "physicalPelvis",
112 logicPelvis,
113 mother,
114 false,
115 0, true);
116
117 // Visualization Attributes
118 // G4VisAttributes* PelvisVisAtt = new G4Vi
119
120 auto* colourPointer = new G4HumanPhantomColo
121 G4Colour colour = colourPointer -> GetColour
122 auto* PelvisVisAtt = new G4VisAttributes(col
123 PelvisVisAtt->SetForceSolid(wireFrame);
124 logicPelvis->SetVisAttributes(PelvisVisAtt);
125
126 G4cout << "Pelvis created !!!!!!" << G4endl;
127
128 // Testing Pelvis Volume
129 G4double PelvisVol = logicPelvis->GetSolid()
130 G4cout << "Volume of Pelvis = " << PelvisVol
131
132 // Testing Pelvis Material
133 G4String PelvisMat = logicPelvis->GetMateria
134 G4cout << "Material of Pelvis = " << PelvisM
135
136 // Testing Density
137 G4double PelvisDensity = logicPelvis->GetMat
138 G4cout << "Density of Material = " << Pelvis
139
140 // Testing Mass
141 G4double PelvisMass = (PelvisVol)*PelvisDens
142 G4cout << "Mass of Pelvis = " << PelvisMass/
143
144 return physPelvis;
145 }
146