| Geant4 Cross Reference |
1 // ******************************************* 1
2 // * License and Disclaimer
3 // *
4 // * The Geant4 software is copyright of th
5 // * the Geant4 Collaboration. It is provided
6 // * conditions of the Geant4 Software License
7 // * LICENSE and available at http://cern.ch/
8 // * include a list of copyright holders.
9 // *
10 // * Neither the authors of this software syst
11 // * institutes,nor the agencies providing fin
12 // * work make any representation or warran
13 // * regarding this software system or assum
14 // * use. Please see the license in the file
15 // * for the full disclaimer and the limitatio
16 // *
17 // * This code implementation is the result
18 // * technical work of the GEANT4 collaboratio
19 // * By using, copying, modifying or distri
20 // * any work based on the software) you ag
21 // * use in resulting scientific publicati
22 // * acceptance of all terms of the Geant4 Sof
23 // *******************************************
24 //
25 //
26 // Author: HoangTRAN, 20/2/2019
27 #include "G4DNABoundingBox.hh"
28 #include <array>
29 #include "G4UnitsTable.hh"
30 using std::array;
31 using std::initializer_list;
32 //....oooOO0OOooo........oooOO0OOooo........oo
33
34 G4DNABoundingBox::G4DNABoundingBox(G4DNABoundi
35 {
36 fxhi = rhs.fxhi;
37 fxlo = rhs.fxlo;
38 fyhi = rhs.fyhi;
39 fylo = rhs.fylo;
40 fzhi = rhs.fzhi;
41 fzlo = rhs.fzlo;
42 }
43 ////....oooOO0OOooo........oooOO0OOooo........
44
45 G4DNABoundingBox::G4DNABoundingBox(const initi
46 {
47 fxhi = 0.;
48 fxlo = 0.;
49 fyhi = 0.;
50 fylo = 0.;
51 fzhi = 0.;
52 fzlo = 0.;
53 std::copy(l.begin(), l.end(), &fxhi);
54 }
55
56 //....oooOO0OOooo........oooOO0OOooo........oo
57
58 G4DNABoundingBox& G4DNABoundingBox::operator=(
59 {
60 fxhi = rhs.fxhi;
61 fxlo = rhs.fxlo;
62 fyhi = rhs.fyhi;
63 fylo = rhs.fylo;
64 fzhi = rhs.fzhi;
65 fzlo = rhs.fzlo;
66 return *this;
67 }
68
69 G4double G4DNABoundingBox::Volume() const
70 {
71 auto LengthY = this->Getyhi() - this->Getylo
72 auto LengthX = this->Getxhi() - this->Getxlo
73 auto LengthZ = this->Getzhi() - this->Getzlo
74 G4double V = LengthY * LengthX * LengthZ;
75 return V;
76 }
77
78 //....oooOO0OOooo........oooOO0OOooo........oo
79
80 void G4DNABoundingBox::resize(G4ThreeVector pi
81 {
82 for(size_t i = 0; i < 8; i++)
83 {
84 const G4ThreeVector& point = pics[i];
85 if(point.x() < fxlo)
86 {
87 fxlo = point.x();
88 }
89 if(point.x() > fxhi)
90 {
91 fxhi = point.x();
92 }
93
94 if(point.y() < fylo)
95 {
96 fylo = point.y();
97 }
98 if(point.y() > fyhi)
99 {
100 fyhi = point.y();
101 }
102
103 if(point.z() < fzlo)
104 {
105 fzlo = point.z();
106 }
107 if(point.z() > fzhi)
108 {
109 fzhi = point.z();
110 }
111 }
112 }
113
114 //--------------------------------------------
115
116 G4DNABoundingBox G4DNABoundingBox::translate(c
117 {
118 G4DNABoundingBox output;
119
120 output.fxhi = fxhi + trans.x();
121 output.fxlo = fxlo + trans.x();
122
123 output.fyhi = fyhi + trans.y();
124 output.fylo = fylo + trans.x();
125
126 output.fzhi = fzhi + trans.z();
127 output.fzlo = fzlo + trans.z();
128
129 return output;
130 }
131 //....oooOO0OOooo........oooOO0OOooo........oo
132
133 G4bool G4DNABoundingBox::contains(const G4DNAB
134 {
135 return fxlo <= other.fxlo && fxhi >= other.f
136 fyhi >= other.fyhi && fzlo <= other.f
137 }
138
139 //....oooOO0OOooo........oooOO0OOooo........oo
140
141 G4bool G4DNABoundingBox::contains(const G4Thre
142 {
143 return fxlo <= point.x() && fxhi >= point.x(
144 fyhi >= point.y() && fzlo <= point.z(
145 }
146
147 //....oooOO0OOooo........oooOO0OOooo........oo
148
149 G4bool G4DNABoundingBox::overlap(const G4DNABo
150 G4DNABounding
151 {
152 if(contains(other))
153 {
154 *out = other;
155 return true;
156 }
157 if(other.contains(*this))
158 {
159 *out = *this;
160 return true;
161 }
162
163 // Check if there is no intersection
164 if(fxhi < other.fxlo || fxlo > other.fxhi ||
165 fylo > other.fyhi || fzhi < other.fzlo ||
166 {
167 *out = invalid;
168 return false;
169 }
170
171 G4double upperX = std::min(fxhi, other.fxhi)
172 G4double upperY = std::min(fyhi, other.fyhi)
173 G4double upperZ = std::min(fzhi, other.fzhi)
174
175 G4double lowerX = std::max(fxlo, other.fxlo)
176 G4double lowerY = std::max(fylo, other.fylo)
177 G4double lowerZ = std::max(fzlo, other.fzlo)
178
179 *out = G4DNABoundingBox{ upperX, lowerX, upp
180 return true;
181 }
182
183 //....oooOO0OOooo........oooOO0OOooo........oo
184
185 G4bool G4DNABoundingBox::overlap(const G4Three
186 const G4doubl
187 {
188 G4double x = query.x() - this->middlePoint()
189 G4double y = query.y() - this->middlePoint()
190 G4double z = query.z() - this->middlePoint()
191
192 x = std::abs(x);
193 y = std::abs(y);
194 z = std::abs(z);
195
196 if((x > (Radius + this->halfSideLengthInX())
197 (y > (Radius + this->halfSideLengthInY())
198 (z > (Radius + this->halfSideLengthInZ())
199 {
200 return false; //
201 }
202 G4int num_less_extent = static_cast<int>(x <
203 static_cast<int>(y <
204 static_cast<int>(z <
205
206 if(num_less_extent > 1)
207 {
208 return true;
209 }
210
211 x = std::max(x - this->halfSideLengthInX(),
212 y = std::max(y - this->halfSideLengthInY(),
213 z = std::max(z - this->halfSideLengthInZ(),
214
215 G4double norm = std::sqrt(x * x + y * y + z
216
217 return (norm < Radius);
218 }
219 //....oooOO0OOooo........oooOO0OOooo........oo
220
221 G4bool G4DNABoundingBox::contains(const G4Thre
222 const G4doub
223 {
224 G4ThreeVector temp = query - this->middlePoi
225 G4double x = std::abs(temp.x()) + th
226 G4double y = std::abs(temp.y()) + th
227 G4double z = std::abs(temp.z()) + th
228 G4double norm = std::sqrt(x * x + y * y
229 return (norm < Radius);
230 }
231 //....oooOO0OOooo........oooOO0OOooo........oo
232
233 G4bool G4DNABoundingBox::contains(const G4Thre
234 const G4Thre
235 const G4doub
236 {
237 return (((query - Point).mag()) < Radius);
238 }
239
240 //....oooOO0OOooo........oooOO0OOooo........oo
241
242 array<G4DNABoundingBox, 8> G4DNABoundingBox::p
243 {
244 G4double xmid = (fxhi + fxlo) / 2.;
245 G4double ymid = (fyhi + fylo) / 2.;
246 G4double zmid = (fzhi + fzlo) / 2.;
247
248 std::array<G4DNABoundingBox, 8> ret{ {
249 G4DNABoundingBox{ xmid, fxlo, ymid, fylo,
250 fzlo }, // bottom left
251 G4DNABoundingBox{ fxhi, xmid, ymid, fylo,
252 fzlo }, // bottom right
253 G4DNABoundingBox{ xmid, fxlo, fyhi, ymid,
254 G4DNABoundingBox{ fxhi, xmid, fyhi, ymid,
255 fzlo }, // bottom right
256 G4DNABoundingBox{ xmid, fxlo, ymid, fylo,
257 G4DNABoundingBox{ fxhi, xmid, ymid, fylo,
258 G4DNABoundingBox{ xmid, fxlo, fyhi, ymid,
259 G4DNABoundingBox{ fxhi, xmid, fyhi, ymid,
260 } };
261 return ret;
262 }
263 //....oooOO0OOooo........oooOO0OOooo........oo
264
265 G4bool G4DNABoundingBox::operator==(const G4DN
266 {
267 return (fxhi == rhs.fxhi && fxlo == rhs.fxlo
268 fylo == rhs.fylo && fzhi == rhs.fzhi
269 (std::isnan(fxhi) && std::isnan(rhs.f
270 std::isnan(rhs.fxlo) && std::isnan(f
271 std::isnan(fylo) && std::isnan(rhs.f
272 std::isnan(rhs.fzhi) && std::isnan(f
273 }
274
275 G4bool G4DNABoundingBox::operator!=(const G4DN
276 {
277 return !operator==(rhs);
278 }
279 //....oooOO0OOooo........oooOO0OOooo........oo
280
281 std::ostream& operator<<(std::ostream& stream,
282 {
283 stream << "{" << G4BestUnit(rhs.fxhi, "Lengt
284 << G4BestUnit(rhs.fxlo, "Length") <<
285 << G4BestUnit(rhs.fyhi, "Length") <<
286 << G4BestUnit(rhs.fylo, "Length") <<
287 << G4BestUnit(rhs.fzhi, "Length") <<
288 << G4BestUnit(rhs.fzlo, "Length") <<
289 << "}";
290 return stream;
291 }
292