| 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 /// \file runAndEvent/RE02/src/RE02NestedPhant
27 /// \brief Implementation of the RE02NestedPha
28 //
29 //
30 //////////////////////////////////////////////
31 #include "RE02NestedPhantomParameterisation.hh
32
33 #include "G4Box.hh"
34 #include "G4LogicalVolume.hh"
35 #include "G4Material.hh"
36 #include "G4ThreeVector.hh"
37 #include "G4VPhysicalVolume.hh"
38 #include "G4VTouchable.hh"
39
40 //============================================
41 // (RE02NestedPhantomParameterisation)
42 //
43 // (Description)
44 // Class for nested parameterisation.
45 // This parameterisation handles material
46 //
47 // T.Aso Created. Nov.2007.
48 //
49 //////////////////////////////////////////////
50
51 //....oooOO0OOooo........oooOO0OOooo........oo
52 RE02NestedPhantomParameterisation ::RE02Nested
53 const G4ThreeVector& voxelSize, G4int nz, st
54 : G4VNestedParameterisation(),
55 fdX(voxelSize.x()),
56 fdY(voxelSize.y()),
57 fdZ(voxelSize.z()),
58 fNz(nz),
59 fMat(mat)
60 {
61 // Position of voxels.
62 // x and y positions are already defined in
63 // by using replicated volume. Here only we
64 // of voxles.
65 fpZ.clear();
66 G4double zp;
67 for (G4int iz = 0; iz < fNz; iz++) {
68 zp = (-fNz + 1 + 2 * iz) * fdZ;
69 fpZ.push_back(zp);
70 }
71 }
72
73 //....oooOO0OOooo........oooOO0OOooo........oo
74 RE02NestedPhantomParameterisation::~RE02Nested
75 {
76 fpZ.clear();
77 }
78
79 //....oooOO0OOooo........oooOO0OOooo........oo
80 //
81 // Material assignment to geometry.
82 //
83 G4Material* RE02NestedPhantomParameterisation
84
85
86 {
87 if (parentTouch == 0)
88 return fMat[0]; // protection for initial
89 // vis at idle state
90 // Copy number of voxels.
91 // Copy number of X and Y are obtained from
92 // Copy nymber of Z is the copy number of cu
93 G4int ix = parentTouch->GetReplicaNumber(0);
94 G4int iy = parentTouch->GetReplicaNumber(1);
95 G4int iz = copyNo;
96 // For demonstration purpose,a couple of mat
97 G4Material* mat = 0;
98 if (ix % 2 == 0 && iy % 2 == 0 && iz % 2 ==
99 mat = fMat[0];
100 else
101 mat = fMat[1];
102
103 return mat;
104 }
105
106 //....oooOO0OOooo........oooOO0OOooo........oo
107 //
108 // Number of Materials
109 // Material scanner is required for preparing
110 // stating simulation, so that G4 has to know
111 G4int RE02NestedPhantomParameterisation::GetNu
112 {
113 return fMat.size();
114 }
115
116 //....oooOO0OOooo........oooOO0OOooo........oo
117 //
118 // GetMaterial
119 // This is needed for material scanner and re
120 //
121 G4Material* RE02NestedPhantomParameterisation:
122 {
123 return fMat[i];
124 }
125
126 //....oooOO0OOooo........oooOO0OOooo........oo
127 //
128 // Transformation of voxels.
129 //
130 void RE02NestedPhantomParameterisation ::Compu
131
132 {
133 G4ThreeVector position(0., 0., fpZ[copyNo]);
134 physVol->SetTranslation(position);
135 }
136
137 //....oooOO0OOooo........oooOO0OOooo........oo
138 //
139 // Dimensions are always same in this RE02 exa
140 //
141 void RE02NestedPhantomParameterisation ::Compu
142
143 {
144 box.SetXHalfLength(fdX);
145 box.SetYHalfLength(fdY);
146 box.SetZHalfLength(fdZ);
147 }
148