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 | // $Id: G4GDMLWriteSolids.cc,v 1.59 2008/11/21 09:32:46 gcosmo Exp $ |
---|
28 | // GEANT4 tag $Name: geant4-09-02-ref-02 $ |
---|
29 | // |
---|
30 | // class G4GDMLWriteSolids Implementation |
---|
31 | // |
---|
32 | // Original author: Zoltan Torzsok, November 2007 |
---|
33 | // |
---|
34 | // -------------------------------------------------------------------- |
---|
35 | |
---|
36 | #include "G4GDMLWriteSolids.hh" |
---|
37 | |
---|
38 | void G4GDMLWriteSolids:: |
---|
39 | BooleanWrite(xercesc::DOMElement* solidsElement, |
---|
40 | const G4BooleanSolid* const boolean) |
---|
41 | { |
---|
42 | G4String tag("undefined"); |
---|
43 | if (dynamic_cast<const G4IntersectionSolid*>(boolean)) |
---|
44 | { tag = "intersection"; } else |
---|
45 | if (dynamic_cast<const G4SubtractionSolid*>(boolean)) |
---|
46 | { tag = "subtraction"; } else |
---|
47 | if (dynamic_cast<const G4UnionSolid*>(boolean)) |
---|
48 | { tag = "union"; } |
---|
49 | |
---|
50 | G4VSolid* firstPtr = const_cast<G4VSolid*>(boolean->GetConstituentSolid(0)); |
---|
51 | G4VSolid* secondPtr = const_cast<G4VSolid*>(boolean->GetConstituentSolid(1)); |
---|
52 | |
---|
53 | G4ThreeVector firstpos,firstrot,pos,rot; |
---|
54 | |
---|
55 | if (const G4DisplacedSolid* disp |
---|
56 | = dynamic_cast<const G4DisplacedSolid*>(firstPtr)) |
---|
57 | { |
---|
58 | firstpos = disp->GetObjectTranslation(); |
---|
59 | firstrot = GetAngles(disp->GetObjectRotation()); |
---|
60 | firstPtr = disp->GetConstituentMovedSolid(); |
---|
61 | } |
---|
62 | |
---|
63 | if (const G4DisplacedSolid* disp |
---|
64 | = dynamic_cast<const G4DisplacedSolid*>(secondPtr)) |
---|
65 | { |
---|
66 | pos = disp->GetObjectTranslation(); |
---|
67 | rot = GetAngles(disp->GetObjectRotation()); |
---|
68 | secondPtr = disp->GetConstituentMovedSolid(); |
---|
69 | } |
---|
70 | |
---|
71 | AddSolid(firstPtr); // At first add the constituent solids! |
---|
72 | AddSolid(secondPtr); |
---|
73 | |
---|
74 | const G4String& name = GenerateName(boolean->GetName(),boolean); |
---|
75 | const G4String& firstref = GenerateName(firstPtr->GetName(),firstPtr); |
---|
76 | const G4String& secondref = GenerateName(secondPtr->GetName(),secondPtr); |
---|
77 | |
---|
78 | xercesc::DOMElement* booleanElement = NewElement(tag); |
---|
79 | booleanElement->setAttributeNode(NewAttribute("name",name)); |
---|
80 | xercesc::DOMElement* firstElement = NewElement("first"); |
---|
81 | firstElement->setAttributeNode(NewAttribute("ref",firstref)); |
---|
82 | booleanElement->appendChild(firstElement); |
---|
83 | xercesc::DOMElement* secondElement = NewElement("second"); |
---|
84 | secondElement->setAttributeNode(NewAttribute("ref",secondref)); |
---|
85 | booleanElement->appendChild(secondElement); |
---|
86 | solidsElement->appendChild(booleanElement); |
---|
87 | // Add the boolean solid AFTER the constituent solids! |
---|
88 | |
---|
89 | if ( (std::fabs(pos.x()) > kLinearPrecision) |
---|
90 | || (std::fabs(pos.y()) > kLinearPrecision) |
---|
91 | || (std::fabs(pos.z()) > kLinearPrecision) ) |
---|
92 | { |
---|
93 | PositionWrite(booleanElement,name+"_pos",pos); |
---|
94 | } |
---|
95 | |
---|
96 | if ( (std::fabs(rot.x()) > kAngularPrecision) |
---|
97 | || (std::fabs(rot.y()) > kAngularPrecision) |
---|
98 | || (std::fabs(rot.z()) > kAngularPrecision) ) |
---|
99 | { |
---|
100 | RotationWrite(booleanElement,name+"_rot",rot); |
---|
101 | } |
---|
102 | |
---|
103 | if ( (std::fabs(firstpos.x()) > kLinearPrecision) |
---|
104 | || (std::fabs(firstpos.y()) > kLinearPrecision) |
---|
105 | || (std::fabs(firstpos.z()) > kLinearPrecision) ) |
---|
106 | { |
---|
107 | FirstpositionWrite(booleanElement,name+"_fpos",firstpos); |
---|
108 | } |
---|
109 | |
---|
110 | if ( (std::fabs(firstrot.x()) > kAngularPrecision) |
---|
111 | || (std::fabs(firstrot.y()) > kAngularPrecision) |
---|
112 | || (std::fabs(firstrot.z()) > kAngularPrecision) ) |
---|
113 | { |
---|
114 | FirstrotationWrite(booleanElement,name+"_frot",firstrot); |
---|
115 | } |
---|
116 | } |
---|
117 | |
---|
118 | void G4GDMLWriteSolids:: |
---|
119 | BoxWrite(xercesc::DOMElement* solidsElement, const G4Box* const box) |
---|
120 | { |
---|
121 | const G4String& name = GenerateName(box->GetName(),box); |
---|
122 | |
---|
123 | xercesc::DOMElement* boxElement = NewElement("box"); |
---|
124 | boxElement->setAttributeNode(NewAttribute("name",name)); |
---|
125 | boxElement->setAttributeNode(NewAttribute("x",2.0*box->GetXHalfLength()/mm)); |
---|
126 | boxElement->setAttributeNode(NewAttribute("y",2.0*box->GetYHalfLength()/mm)); |
---|
127 | boxElement->setAttributeNode(NewAttribute("z",2.0*box->GetZHalfLength()/mm)); |
---|
128 | boxElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
129 | solidsElement->appendChild(boxElement); |
---|
130 | } |
---|
131 | |
---|
132 | void G4GDMLWriteSolids:: |
---|
133 | ConeWrite(xercesc::DOMElement* solidsElement, const G4Cons* const cone) |
---|
134 | { |
---|
135 | const G4String& name = GenerateName(cone->GetName(),cone); |
---|
136 | |
---|
137 | xercesc::DOMElement* coneElement = NewElement("cone"); |
---|
138 | coneElement->setAttributeNode(NewAttribute("name",name)); |
---|
139 | coneElement-> |
---|
140 | setAttributeNode(NewAttribute("rmin1",cone->GetInnerRadiusMinusZ()/mm)); |
---|
141 | coneElement-> |
---|
142 | setAttributeNode(NewAttribute("rmax1",cone->GetOuterRadiusMinusZ()/mm)); |
---|
143 | coneElement-> |
---|
144 | setAttributeNode(NewAttribute("rmin2",cone->GetInnerRadiusPlusZ()/mm)); |
---|
145 | coneElement-> |
---|
146 | setAttributeNode(NewAttribute("rmax2",cone->GetOuterRadiusPlusZ()/mm)); |
---|
147 | coneElement-> |
---|
148 | setAttributeNode(NewAttribute("z",2.0*cone->GetZHalfLength()/mm)); |
---|
149 | coneElement-> |
---|
150 | setAttributeNode(NewAttribute("startphi",cone->GetStartPhiAngle()/degree)); |
---|
151 | coneElement-> |
---|
152 | setAttributeNode(NewAttribute("deltaphi",cone->GetDeltaPhiAngle()/degree)); |
---|
153 | coneElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
154 | coneElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
155 | solidsElement->appendChild(coneElement); |
---|
156 | } |
---|
157 | |
---|
158 | void G4GDMLWriteSolids:: |
---|
159 | ElconeWrite(xercesc::DOMElement* solidsElement, |
---|
160 | const G4EllipticalCone* const elcone) |
---|
161 | { |
---|
162 | const G4String& name = GenerateName(elcone->GetName(),elcone); |
---|
163 | |
---|
164 | xercesc::DOMElement* elconeElement = NewElement("elcone"); |
---|
165 | elconeElement->setAttributeNode(NewAttribute("name",name)); |
---|
166 | elconeElement->setAttributeNode(NewAttribute("dx",elcone->GetSemiAxisX()/mm)); |
---|
167 | elconeElement->setAttributeNode(NewAttribute("dy",elcone->GetSemiAxisY()/mm)); |
---|
168 | elconeElement->setAttributeNode(NewAttribute("zmax",elcone->GetZMax()/mm)); |
---|
169 | elconeElement->setAttributeNode(NewAttribute("zcut",elcone->GetZTopCut()/mm)); |
---|
170 | elconeElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
171 | solidsElement->appendChild(elconeElement); |
---|
172 | } |
---|
173 | |
---|
174 | void G4GDMLWriteSolids:: |
---|
175 | EllipsoidWrite(xercesc::DOMElement* solidsElement, |
---|
176 | const G4Ellipsoid* const ellipsoid) |
---|
177 | { |
---|
178 | const G4String& name = GenerateName(ellipsoid->GetName(),ellipsoid); |
---|
179 | |
---|
180 | xercesc::DOMElement* ellipsoidElement = NewElement("ellipsoid"); |
---|
181 | ellipsoidElement->setAttributeNode(NewAttribute("name",name)); |
---|
182 | ellipsoidElement-> |
---|
183 | setAttributeNode(NewAttribute("ax",ellipsoid->GetSemiAxisMax(0)/mm)); |
---|
184 | ellipsoidElement-> |
---|
185 | setAttributeNode(NewAttribute("by",ellipsoid->GetSemiAxisMax(1)/mm)); |
---|
186 | ellipsoidElement-> |
---|
187 | setAttributeNode(NewAttribute("cz",ellipsoid->GetSemiAxisMax(2)/mm)); |
---|
188 | ellipsoidElement-> |
---|
189 | setAttributeNode(NewAttribute("zcut1",ellipsoid->GetZBottomCut()/mm)); |
---|
190 | ellipsoidElement-> |
---|
191 | setAttributeNode(NewAttribute("zcut2",ellipsoid->GetZTopCut()/mm)); |
---|
192 | ellipsoidElement-> |
---|
193 | setAttributeNode(NewAttribute("lunit","mm")); |
---|
194 | solidsElement->appendChild(ellipsoidElement); |
---|
195 | } |
---|
196 | |
---|
197 | void G4GDMLWriteSolids:: |
---|
198 | EltubeWrite(xercesc::DOMElement* solidsElement, |
---|
199 | const G4EllipticalTube* const eltube) |
---|
200 | { |
---|
201 | const G4String& name = GenerateName(eltube->GetName(),eltube); |
---|
202 | |
---|
203 | xercesc::DOMElement* eltubeElement = NewElement("eltube"); |
---|
204 | eltubeElement->setAttributeNode(NewAttribute("name",name)); |
---|
205 | eltubeElement->setAttributeNode(NewAttribute("dx",eltube->GetDx()/mm)); |
---|
206 | eltubeElement->setAttributeNode(NewAttribute("dy",eltube->GetDy()/mm)); |
---|
207 | eltubeElement->setAttributeNode(NewAttribute("dz",eltube->GetDz()/mm)); |
---|
208 | eltubeElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
209 | solidsElement->appendChild(eltubeElement); |
---|
210 | } |
---|
211 | |
---|
212 | void G4GDMLWriteSolids:: |
---|
213 | XtruWrite(xercesc::DOMElement* solidsElement, |
---|
214 | const G4ExtrudedSolid* const xtru) |
---|
215 | { |
---|
216 | const G4String& name = GenerateName(xtru->GetName(),xtru); |
---|
217 | |
---|
218 | xercesc::DOMElement* xtruElement = NewElement("xtru"); |
---|
219 | xtruElement->setAttributeNode(NewAttribute("name",name)); |
---|
220 | xtruElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
221 | solidsElement->appendChild(xtruElement); |
---|
222 | |
---|
223 | const G4int NumVertex = xtru->GetNofVertices(); |
---|
224 | |
---|
225 | for (G4int i=0;i<NumVertex;i++) |
---|
226 | { |
---|
227 | xercesc::DOMElement* twoDimVertexElement = NewElement("twoDimVertex"); |
---|
228 | xtruElement->appendChild(twoDimVertexElement); |
---|
229 | |
---|
230 | const G4TwoVector& vertex = xtru->GetVertex(i); |
---|
231 | |
---|
232 | twoDimVertexElement->setAttributeNode(NewAttribute("x",vertex.x()/mm)); |
---|
233 | twoDimVertexElement->setAttributeNode(NewAttribute("y",vertex.y()/mm)); |
---|
234 | } |
---|
235 | |
---|
236 | const G4int NumSection = xtru->GetNofZSections(); |
---|
237 | |
---|
238 | for (G4int i=0;i<NumSection;i++) |
---|
239 | { |
---|
240 | xercesc::DOMElement* sectionElement = NewElement("section"); |
---|
241 | xtruElement->appendChild(sectionElement); |
---|
242 | |
---|
243 | const G4ExtrudedSolid::ZSection section = xtru->GetZSection(i); |
---|
244 | |
---|
245 | sectionElement->setAttributeNode(NewAttribute("zOrder",i)); |
---|
246 | sectionElement->setAttributeNode(NewAttribute("zPosition",section.fZ/mm)); |
---|
247 | sectionElement-> |
---|
248 | setAttributeNode(NewAttribute("xOffset",section.fOffset.x()/mm)); |
---|
249 | sectionElement-> |
---|
250 | setAttributeNode(NewAttribute("yOffset",section.fOffset.y()/mm)); |
---|
251 | sectionElement-> |
---|
252 | setAttributeNode(NewAttribute("scalingFactor",section.fScale)); |
---|
253 | } |
---|
254 | } |
---|
255 | |
---|
256 | void G4GDMLWriteSolids:: |
---|
257 | HypeWrite(xercesc::DOMElement* solidsElement, const G4Hype* const hype) |
---|
258 | { |
---|
259 | const G4String& name = GenerateName(hype->GetName(),hype); |
---|
260 | |
---|
261 | xercesc::DOMElement* hypeElement = NewElement("hype"); |
---|
262 | hypeElement->setAttributeNode(NewAttribute("name",name)); |
---|
263 | hypeElement->setAttributeNode(NewAttribute("rmin", |
---|
264 | hype->GetInnerRadius()/mm)); |
---|
265 | hypeElement->setAttributeNode(NewAttribute("rmax", |
---|
266 | hype->GetOuterRadius()/mm)); |
---|
267 | hypeElement->setAttributeNode(NewAttribute("inst", |
---|
268 | hype->GetInnerStereo()/degree)); |
---|
269 | hypeElement->setAttributeNode(NewAttribute("outst", |
---|
270 | hype->GetOuterStereo()/degree)); |
---|
271 | hypeElement->setAttributeNode(NewAttribute("z", |
---|
272 | 2.0*hype->GetZHalfLength()/mm)); |
---|
273 | hypeElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
274 | hypeElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
275 | solidsElement->appendChild(hypeElement); |
---|
276 | } |
---|
277 | |
---|
278 | void G4GDMLWriteSolids:: |
---|
279 | OrbWrite(xercesc::DOMElement* solidsElement, const G4Orb* const orb) |
---|
280 | { |
---|
281 | const G4String& name = GenerateName(orb->GetName(),orb); |
---|
282 | |
---|
283 | xercesc::DOMElement* orbElement = NewElement("orb"); |
---|
284 | orbElement->setAttributeNode(NewAttribute("name",name)); |
---|
285 | orbElement->setAttributeNode(NewAttribute("r",orb->GetRadius()/mm)); |
---|
286 | orbElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
287 | solidsElement->appendChild(orbElement); |
---|
288 | } |
---|
289 | |
---|
290 | void G4GDMLWriteSolids:: |
---|
291 | ParaWrite(xercesc::DOMElement* solidsElement, const G4Para* const para) |
---|
292 | { |
---|
293 | const G4String& name = GenerateName(para->GetName(),para); |
---|
294 | |
---|
295 | const G4ThreeVector simaxis = para->GetSymAxis(); |
---|
296 | const G4double alpha = std::atan(para->GetTanAlpha()); |
---|
297 | const G4double theta = std::acos(simaxis.z()); |
---|
298 | const G4double phi = (simaxis.z() != 1.0) |
---|
299 | ? (std::atan(simaxis.y()/simaxis.x())) : (0.0); |
---|
300 | |
---|
301 | xercesc::DOMElement* paraElement = NewElement("para"); |
---|
302 | paraElement->setAttributeNode(NewAttribute("name",name)); |
---|
303 | paraElement->setAttributeNode(NewAttribute("x", |
---|
304 | 2.0*para->GetXHalfLength()/mm)); |
---|
305 | paraElement->setAttributeNode(NewAttribute("y", |
---|
306 | 2.0*para->GetYHalfLength()/mm)); |
---|
307 | paraElement->setAttributeNode(NewAttribute("z", |
---|
308 | 2.0*para->GetZHalfLength()/mm)); |
---|
309 | paraElement->setAttributeNode(NewAttribute("alpha",alpha/degree)); |
---|
310 | paraElement->setAttributeNode(NewAttribute("theta",theta/degree)); |
---|
311 | paraElement->setAttributeNode(NewAttribute("phi",phi/degree)); |
---|
312 | paraElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
313 | paraElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
314 | solidsElement->appendChild(paraElement); |
---|
315 | } |
---|
316 | |
---|
317 | void G4GDMLWriteSolids:: |
---|
318 | ParaboloidWrite(xercesc::DOMElement* solidsElement, |
---|
319 | const G4Paraboloid* const paraboloid) |
---|
320 | { |
---|
321 | const G4String& name = GenerateName(paraboloid->GetName(),paraboloid); |
---|
322 | |
---|
323 | xercesc::DOMElement* paraboloidElement = NewElement("paraboloid"); |
---|
324 | paraboloidElement->setAttributeNode(NewAttribute("name",name)); |
---|
325 | paraboloidElement->setAttributeNode(NewAttribute("rlo", |
---|
326 | paraboloid->GetRadiusMinusZ()/mm)); |
---|
327 | paraboloidElement->setAttributeNode(NewAttribute("rhi", |
---|
328 | paraboloid->GetRadiusPlusZ()/mm)); |
---|
329 | paraboloidElement->setAttributeNode(NewAttribute("dz", |
---|
330 | paraboloid->GetZHalfLength()/mm)); |
---|
331 | paraboloidElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
332 | solidsElement->appendChild(paraboloidElement); |
---|
333 | } |
---|
334 | |
---|
335 | void G4GDMLWriteSolids:: |
---|
336 | PolyconeWrite(xercesc::DOMElement* solidsElement, |
---|
337 | const G4Polycone* const polycone) |
---|
338 | { |
---|
339 | const G4String& name = GenerateName(polycone->GetName(),polycone); |
---|
340 | |
---|
341 | xercesc::DOMElement* polyconeElement = NewElement("polycone"); |
---|
342 | polyconeElement->setAttributeNode(NewAttribute("name",name)); |
---|
343 | polyconeElement->setAttributeNode(NewAttribute("startphi", |
---|
344 | polycone->GetOriginalParameters()->Start_angle/degree)); |
---|
345 | polyconeElement->setAttributeNode(NewAttribute("deltaphi", |
---|
346 | polycone->GetOriginalParameters()->Opening_angle/degree)); |
---|
347 | polyconeElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
348 | polyconeElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
349 | solidsElement->appendChild(polyconeElement); |
---|
350 | |
---|
351 | const size_t num_zplanes = polycone->GetOriginalParameters()->Num_z_planes; |
---|
352 | const G4double* z_array = polycone->GetOriginalParameters()->Z_values; |
---|
353 | const G4double* rmin_array = polycone->GetOriginalParameters()->Rmin; |
---|
354 | const G4double* rmax_array = polycone->GetOriginalParameters()->Rmax; |
---|
355 | |
---|
356 | for (size_t i=0; i<num_zplanes; i++) |
---|
357 | { |
---|
358 | ZplaneWrite(polyconeElement,z_array[i],rmin_array[i],rmax_array[i]); |
---|
359 | } |
---|
360 | } |
---|
361 | |
---|
362 | void G4GDMLWriteSolids:: |
---|
363 | PolyhedraWrite(xercesc::DOMElement* solidsElement, |
---|
364 | const G4Polyhedra* const polyhedra) |
---|
365 | { |
---|
366 | const G4String& name = GenerateName(polyhedra->GetName(),polyhedra); |
---|
367 | |
---|
368 | xercesc::DOMElement* polyhedraElement = NewElement("polyhedra"); |
---|
369 | polyhedraElement->setAttributeNode(NewAttribute("name",name)); |
---|
370 | polyhedraElement->setAttributeNode(NewAttribute("startphi", |
---|
371 | polyhedra->GetOriginalParameters()->Start_angle/degree)); |
---|
372 | polyhedraElement->setAttributeNode(NewAttribute("deltaphi", |
---|
373 | polyhedra->GetOriginalParameters()->Opening_angle/degree)); |
---|
374 | polyhedraElement->setAttributeNode(NewAttribute("numsides", |
---|
375 | polyhedra->GetOriginalParameters()->numSide)); |
---|
376 | polyhedraElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
377 | polyhedraElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
378 | solidsElement->appendChild(polyhedraElement); |
---|
379 | |
---|
380 | const size_t num_zplanes = polyhedra->GetOriginalParameters()->Num_z_planes; |
---|
381 | const G4double* z_array = polyhedra->GetOriginalParameters()->Z_values; |
---|
382 | const G4double* rmin_array = polyhedra->GetOriginalParameters()->Rmin; |
---|
383 | const G4double* rmax_array = polyhedra->GetOriginalParameters()->Rmax; |
---|
384 | |
---|
385 | const G4double convertRad = |
---|
386 | std::cos(0.5*polyhedra->GetOriginalParameters()->Opening_angle |
---|
387 | / polyhedra->GetOriginalParameters()->numSide); |
---|
388 | |
---|
389 | for (size_t i=0;i<num_zplanes;i++) |
---|
390 | { |
---|
391 | ZplaneWrite(polyhedraElement,z_array[i], |
---|
392 | rmin_array[i]*convertRad, rmax_array[i]*convertRad); |
---|
393 | } |
---|
394 | } |
---|
395 | |
---|
396 | void G4GDMLWriteSolids:: |
---|
397 | SphereWrite(xercesc::DOMElement* solidsElement, const G4Sphere* const sphere) |
---|
398 | { |
---|
399 | const G4String& name = GenerateName(sphere->GetName(),sphere); |
---|
400 | |
---|
401 | xercesc::DOMElement* sphereElement = NewElement("sphere"); |
---|
402 | sphereElement->setAttributeNode(NewAttribute("name",name)); |
---|
403 | sphereElement->setAttributeNode(NewAttribute("rmin", |
---|
404 | sphere->GetInsideRadius()/mm)); |
---|
405 | sphereElement->setAttributeNode(NewAttribute("rmax", |
---|
406 | sphere->GetOuterRadius()/mm)); |
---|
407 | sphereElement->setAttributeNode(NewAttribute("startphi", |
---|
408 | sphere->GetStartPhiAngle()/degree)); |
---|
409 | sphereElement->setAttributeNode(NewAttribute("deltaphi", |
---|
410 | sphere->GetDeltaPhiAngle()/degree)); |
---|
411 | sphereElement->setAttributeNode(NewAttribute("starttheta", |
---|
412 | sphere->GetStartThetaAngle()/degree)); |
---|
413 | sphereElement->setAttributeNode(NewAttribute("deltatheta", |
---|
414 | sphere->GetDeltaThetaAngle()/degree)); |
---|
415 | sphereElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
416 | sphereElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
417 | solidsElement->appendChild(sphereElement); |
---|
418 | } |
---|
419 | |
---|
420 | void G4GDMLWriteSolids:: |
---|
421 | TessellatedWrite(xercesc::DOMElement* solidsElement, |
---|
422 | const G4TessellatedSolid* const tessellated) |
---|
423 | { |
---|
424 | const G4String& solid_name = tessellated->GetName(); |
---|
425 | const G4String& name = GenerateName(solid_name, tessellated); |
---|
426 | |
---|
427 | xercesc::DOMElement* tessellatedElement = NewElement("tessellated"); |
---|
428 | tessellatedElement->setAttributeNode(NewAttribute("name",name)); |
---|
429 | tessellatedElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
430 | tessellatedElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
431 | solidsElement->appendChild(tessellatedElement); |
---|
432 | |
---|
433 | std::map<G4ThreeVector, G4String> vertexMap; |
---|
434 | |
---|
435 | const size_t NumFacets = tessellated->GetNumberOfFacets(); |
---|
436 | size_t NumVertex = 0; |
---|
437 | |
---|
438 | for (size_t i=0;i<NumFacets;i++) |
---|
439 | { |
---|
440 | const G4VFacet* facet = tessellated->GetFacet(i); |
---|
441 | const size_t NumVertexPerFacet = facet->GetNumberOfVertices(); |
---|
442 | |
---|
443 | G4String FacetTag; |
---|
444 | |
---|
445 | if (NumVertexPerFacet==3) { FacetTag="triangular"; } else |
---|
446 | if (NumVertexPerFacet==4) { FacetTag="quadrangular"; } |
---|
447 | else |
---|
448 | { |
---|
449 | G4Exception("G4GDMLWriteSolids::TessellatedWrite()", "InvalidSetup", |
---|
450 | FatalException, "Facet should contain 3 or 4 vertices!"); |
---|
451 | } |
---|
452 | |
---|
453 | xercesc::DOMElement* facetElement = NewElement(FacetTag); |
---|
454 | tessellatedElement->appendChild(facetElement); |
---|
455 | |
---|
456 | for (size_t j=0; j<NumVertexPerFacet; j++) |
---|
457 | { |
---|
458 | std::stringstream name_stream; |
---|
459 | std::stringstream ref_stream; |
---|
460 | |
---|
461 | name_stream << "vertex" << (j+1); |
---|
462 | ref_stream << solid_name << "_v" << NumVertex; |
---|
463 | |
---|
464 | const G4String& name = name_stream.str(); // facet's tag variable |
---|
465 | G4String ref = ref_stream.str(); // vertex tag to be associated |
---|
466 | |
---|
467 | // Now search for the existance of the current vertex in the |
---|
468 | // map of cached vertices. If existing, do NOT store it as |
---|
469 | // position in the GDML file, so avoiding duplication; otherwise |
---|
470 | // cache it in the local map and add it as position in the |
---|
471 | // "define" section of the GDML file. |
---|
472 | |
---|
473 | const G4ThreeVector& vertex = facet->GetVertex(j); |
---|
474 | |
---|
475 | if(vertexMap.find(vertex) != vertexMap.end()) // Vertex is cached |
---|
476 | { |
---|
477 | ref = vertexMap[vertex]; // Set the proper tag for it |
---|
478 | } |
---|
479 | else // Vertex not found |
---|
480 | { |
---|
481 | vertexMap.insert(std::make_pair(vertex,ref)); // Cache vertex and ... |
---|
482 | AddPosition(ref, vertex); // ... add it to define section! |
---|
483 | NumVertex++; |
---|
484 | } |
---|
485 | |
---|
486 | // Now create association of the vertex with its facet |
---|
487 | // |
---|
488 | facetElement->setAttributeNode(NewAttribute(name,ref)); |
---|
489 | } |
---|
490 | } |
---|
491 | } |
---|
492 | |
---|
493 | void G4GDMLWriteSolids:: |
---|
494 | TetWrite(xercesc::DOMElement* solidsElement, const G4Tet* const tet) |
---|
495 | { |
---|
496 | const G4String& solid_name = tet->GetName(); |
---|
497 | const G4String& name = GenerateName(solid_name, tet); |
---|
498 | |
---|
499 | std::vector<G4ThreeVector> vertexList = tet->GetVertices(); |
---|
500 | |
---|
501 | xercesc::DOMElement* tetElement = NewElement("tet"); |
---|
502 | tetElement->setAttributeNode(NewAttribute("name",name)); |
---|
503 | tetElement->setAttributeNode(NewAttribute("vertex1",solid_name+"_v1")); |
---|
504 | tetElement->setAttributeNode(NewAttribute("vertex2",solid_name+"_v2")); |
---|
505 | tetElement->setAttributeNode(NewAttribute("vertex3",solid_name+"_v3")); |
---|
506 | tetElement->setAttributeNode(NewAttribute("vertex4",solid_name+"_v4")); |
---|
507 | tetElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
508 | solidsElement->appendChild(tetElement); |
---|
509 | |
---|
510 | AddPosition(solid_name+"_v1",vertexList[0]); |
---|
511 | AddPosition(solid_name+"_v2",vertexList[1]); |
---|
512 | AddPosition(solid_name+"_v3",vertexList[2]); |
---|
513 | AddPosition(solid_name+"_v4",vertexList[3]); |
---|
514 | } |
---|
515 | |
---|
516 | void G4GDMLWriteSolids:: |
---|
517 | TorusWrite(xercesc::DOMElement* solidsElement, const G4Torus* const torus) |
---|
518 | { |
---|
519 | const G4String& name = GenerateName(torus->GetName(),torus); |
---|
520 | |
---|
521 | xercesc::DOMElement* torusElement = NewElement("torus"); |
---|
522 | torusElement->setAttributeNode(NewAttribute("name",name)); |
---|
523 | torusElement->setAttributeNode(NewAttribute("rmin",torus->GetRmin()/mm)); |
---|
524 | torusElement->setAttributeNode(NewAttribute("rmax",torus->GetRmax()/mm)); |
---|
525 | torusElement->setAttributeNode(NewAttribute("rtor",torus->GetRtor()/mm)); |
---|
526 | torusElement-> |
---|
527 | setAttributeNode(NewAttribute("startphi",torus->GetSPhi()/degree)); |
---|
528 | torusElement-> |
---|
529 | setAttributeNode(NewAttribute("deltaphi",torus->GetDPhi()/degree)); |
---|
530 | torusElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
531 | torusElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
532 | solidsElement->appendChild(torusElement); |
---|
533 | } |
---|
534 | |
---|
535 | void G4GDMLWriteSolids:: |
---|
536 | TrapWrite(xercesc::DOMElement* solidsElement, const G4Trap* const trap) |
---|
537 | { |
---|
538 | const G4String& name = GenerateName(trap->GetName(),trap); |
---|
539 | |
---|
540 | const G4ThreeVector& simaxis = trap->GetSymAxis(); |
---|
541 | const G4double phi = (simaxis.z() != 1.0) |
---|
542 | ? (std::atan(simaxis.y()/simaxis.x())) : (0.0); |
---|
543 | const G4double theta = std::acos(simaxis.z()); |
---|
544 | const G4double alpha1 = std::atan(trap->GetTanAlpha1()); |
---|
545 | const G4double alpha2 = std::atan(trap->GetTanAlpha2()); |
---|
546 | |
---|
547 | xercesc::DOMElement* trapElement = NewElement("trap"); |
---|
548 | trapElement->setAttributeNode(NewAttribute("name",name)); |
---|
549 | trapElement->setAttributeNode(NewAttribute("z", |
---|
550 | 2.0*trap->GetZHalfLength()/mm)); |
---|
551 | trapElement->setAttributeNode(NewAttribute("theta",theta/degree)); |
---|
552 | trapElement->setAttributeNode(NewAttribute("phi",phi/degree)); |
---|
553 | trapElement->setAttributeNode(NewAttribute("y1", |
---|
554 | 2.0*trap->GetYHalfLength1()/mm)); |
---|
555 | trapElement->setAttributeNode(NewAttribute("x1", |
---|
556 | 2.0*trap->GetXHalfLength1()/mm)); |
---|
557 | trapElement->setAttributeNode(NewAttribute("x2", |
---|
558 | 2.0*trap->GetXHalfLength2()/mm)); |
---|
559 | trapElement->setAttributeNode(NewAttribute("alpha1",alpha1/degree)); |
---|
560 | trapElement->setAttributeNode(NewAttribute("y2", |
---|
561 | 2.0*trap->GetYHalfLength2()/mm)); |
---|
562 | trapElement->setAttributeNode(NewAttribute("x3", |
---|
563 | 2.0*trap->GetXHalfLength3()/mm)); |
---|
564 | trapElement->setAttributeNode(NewAttribute("x4", |
---|
565 | 2.0*trap->GetXHalfLength4()/mm)); |
---|
566 | trapElement->setAttributeNode(NewAttribute("alpha2",alpha2/degree)); |
---|
567 | trapElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
568 | trapElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
569 | solidsElement->appendChild(trapElement); |
---|
570 | } |
---|
571 | |
---|
572 | void G4GDMLWriteSolids:: |
---|
573 | TrdWrite(xercesc::DOMElement* solidsElement, const G4Trd* const trd) |
---|
574 | { |
---|
575 | const G4String& name = GenerateName(trd->GetName(),trd); |
---|
576 | |
---|
577 | xercesc::DOMElement* trdElement = NewElement("trd"); |
---|
578 | trdElement->setAttributeNode(NewAttribute("name",name)); |
---|
579 | trdElement->setAttributeNode(NewAttribute("x1", |
---|
580 | 2.0*trd->GetXHalfLength1()/mm)); |
---|
581 | trdElement->setAttributeNode(NewAttribute("x2", |
---|
582 | 2.0*trd->GetXHalfLength2()/mm)); |
---|
583 | trdElement->setAttributeNode(NewAttribute("y1", |
---|
584 | 2.0*trd->GetYHalfLength1()/mm)); |
---|
585 | trdElement->setAttributeNode(NewAttribute("y2", |
---|
586 | 2.0*trd->GetYHalfLength2()/mm)); |
---|
587 | trdElement->setAttributeNode(NewAttribute("z", |
---|
588 | 2.0*trd->GetZHalfLength()/mm)); |
---|
589 | trdElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
590 | solidsElement->appendChild(trdElement); |
---|
591 | } |
---|
592 | |
---|
593 | void G4GDMLWriteSolids:: |
---|
594 | TubeWrite(xercesc::DOMElement* solidsElement, const G4Tubs* const tube) |
---|
595 | { |
---|
596 | const G4String& name = GenerateName(tube->GetName(),tube); |
---|
597 | |
---|
598 | xercesc::DOMElement* tubeElement = NewElement("tube"); |
---|
599 | tubeElement->setAttributeNode(NewAttribute("name",name)); |
---|
600 | tubeElement->setAttributeNode(NewAttribute("rmin", |
---|
601 | tube->GetInnerRadius()/mm)); |
---|
602 | tubeElement->setAttributeNode(NewAttribute("rmax", |
---|
603 | tube->GetOuterRadius()/mm)); |
---|
604 | tubeElement->setAttributeNode(NewAttribute("z", |
---|
605 | 2.0*tube->GetZHalfLength()/mm)); |
---|
606 | tubeElement->setAttributeNode(NewAttribute("startphi", |
---|
607 | tube->GetStartPhiAngle()/degree)); |
---|
608 | tubeElement->setAttributeNode(NewAttribute("deltaphi", |
---|
609 | tube->GetDeltaPhiAngle()/degree)); |
---|
610 | tubeElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
611 | tubeElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
612 | solidsElement->appendChild(tubeElement); |
---|
613 | } |
---|
614 | |
---|
615 | void G4GDMLWriteSolids:: |
---|
616 | TwistedboxWrite(xercesc::DOMElement* solidsElement, |
---|
617 | const G4TwistedBox* const twistedbox) |
---|
618 | { |
---|
619 | const G4String& name = GenerateName(twistedbox->GetName(),twistedbox); |
---|
620 | |
---|
621 | xercesc::DOMElement* twistedboxElement = NewElement("twistedbox"); |
---|
622 | twistedboxElement->setAttributeNode(NewAttribute("name",name)); |
---|
623 | twistedboxElement->setAttributeNode(NewAttribute("x", |
---|
624 | 2.0*twistedbox->GetXHalfLength()/mm)); |
---|
625 | twistedboxElement->setAttributeNode(NewAttribute("y", |
---|
626 | 2.0*twistedbox->GetYHalfLength()/mm)); |
---|
627 | twistedboxElement->setAttributeNode(NewAttribute("z", |
---|
628 | 2.0*twistedbox->GetZHalfLength()/mm)); |
---|
629 | twistedboxElement->setAttributeNode(NewAttribute("PhiTwist", |
---|
630 | twistedbox->GetPhiTwist()/degree)); |
---|
631 | twistedboxElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
632 | twistedboxElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
633 | solidsElement->appendChild(twistedboxElement); |
---|
634 | } |
---|
635 | |
---|
636 | void G4GDMLWriteSolids:: |
---|
637 | TwistedtrapWrite(xercesc::DOMElement* solidsElement, |
---|
638 | const G4TwistedTrap* const twistedtrap) |
---|
639 | { |
---|
640 | const G4String& name = GenerateName(twistedtrap->GetName(),twistedtrap); |
---|
641 | |
---|
642 | xercesc::DOMElement* twistedtrapElement = NewElement("twistedtrap"); |
---|
643 | twistedtrapElement->setAttributeNode(NewAttribute("name",name)); |
---|
644 | twistedtrapElement->setAttributeNode(NewAttribute("y1", |
---|
645 | 2.0*twistedtrap->GetY1HalfLength()/mm)); |
---|
646 | twistedtrapElement->setAttributeNode(NewAttribute("x1", |
---|
647 | 2.0*twistedtrap->GetX1HalfLength()/mm)); |
---|
648 | twistedtrapElement->setAttributeNode(NewAttribute("x2", |
---|
649 | 2.0*twistedtrap->GetX2HalfLength()/mm)); |
---|
650 | twistedtrapElement->setAttributeNode(NewAttribute("y2", |
---|
651 | 2.0*twistedtrap->GetY2HalfLength()/mm)); |
---|
652 | twistedtrapElement->setAttributeNode(NewAttribute("x3", |
---|
653 | 2.0*twistedtrap->GetX3HalfLength()/mm)); |
---|
654 | twistedtrapElement->setAttributeNode(NewAttribute("x4", |
---|
655 | 2.0*twistedtrap->GetX4HalfLength()/mm)); |
---|
656 | twistedtrapElement->setAttributeNode(NewAttribute("z", |
---|
657 | 2.0*twistedtrap->GetZHalfLength()/mm)); |
---|
658 | twistedtrapElement->setAttributeNode(NewAttribute("Alph", |
---|
659 | twistedtrap->GetTiltAngleAlpha()/degree)); |
---|
660 | twistedtrapElement->setAttributeNode(NewAttribute("Theta", |
---|
661 | twistedtrap->GetPolarAngleTheta()/degree)); |
---|
662 | twistedtrapElement->setAttributeNode(NewAttribute("Phi", |
---|
663 | twistedtrap->GetAzimuthalAnglePhi()/degree)); |
---|
664 | twistedtrapElement->setAttributeNode(NewAttribute("PhiTwist", |
---|
665 | twistedtrap->GetPhiTwist()/degree)); |
---|
666 | twistedtrapElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
667 | twistedtrapElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
668 | |
---|
669 | solidsElement->appendChild(twistedtrapElement); |
---|
670 | } |
---|
671 | |
---|
672 | void G4GDMLWriteSolids:: |
---|
673 | TwistedtrdWrite(xercesc::DOMElement* solidsElement, |
---|
674 | const G4TwistedTrd* const twistedtrd) |
---|
675 | { |
---|
676 | const G4String& name = GenerateName(twistedtrd->GetName(),twistedtrd); |
---|
677 | |
---|
678 | xercesc::DOMElement* twistedtrdElement = NewElement("twistedtrd"); |
---|
679 | twistedtrdElement->setAttributeNode(NewAttribute("name",name)); |
---|
680 | twistedtrdElement->setAttributeNode(NewAttribute("x1", |
---|
681 | 2.0*twistedtrd->GetX1HalfLength()/mm)); |
---|
682 | twistedtrdElement->setAttributeNode(NewAttribute("x2", |
---|
683 | 2.0*twistedtrd->GetX2HalfLength()/mm)); |
---|
684 | twistedtrdElement->setAttributeNode(NewAttribute("y1", |
---|
685 | 2.0*twistedtrd->GetY1HalfLength()/mm)); |
---|
686 | twistedtrdElement->setAttributeNode(NewAttribute("y2", |
---|
687 | 2.0*twistedtrd->GetY2HalfLength()/mm)); |
---|
688 | twistedtrdElement->setAttributeNode(NewAttribute("z", |
---|
689 | 2.0*twistedtrd->GetZHalfLength()/mm)); |
---|
690 | twistedtrdElement->setAttributeNode(NewAttribute("PhiTwist", |
---|
691 | twistedtrd->GetPhiTwist()/degree)); |
---|
692 | twistedtrdElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
693 | twistedtrdElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
694 | solidsElement->appendChild(twistedtrdElement); |
---|
695 | } |
---|
696 | |
---|
697 | void G4GDMLWriteSolids:: |
---|
698 | TwistedtubsWrite(xercesc::DOMElement* solidsElement, |
---|
699 | const G4TwistedTubs* const twistedtubs) |
---|
700 | { |
---|
701 | const G4String& name = GenerateName(twistedtubs->GetName(),twistedtubs); |
---|
702 | |
---|
703 | xercesc::DOMElement* twistedtubsElement = NewElement("twistedtubs"); |
---|
704 | twistedtubsElement->setAttributeNode(NewAttribute("name",name)); |
---|
705 | twistedtubsElement->setAttributeNode(NewAttribute("twistedangle", |
---|
706 | twistedtubs->GetPhiTwist()/degree)); |
---|
707 | twistedtubsElement->setAttributeNode(NewAttribute("endinnerrad", |
---|
708 | twistedtubs->GetInnerRadius()/mm)); |
---|
709 | twistedtubsElement->setAttributeNode(NewAttribute("endouterrad", |
---|
710 | twistedtubs->GetOuterRadius()/mm)); |
---|
711 | twistedtubsElement->setAttributeNode(NewAttribute("zlen", |
---|
712 | 2.0*twistedtubs->GetZHalfLength()/mm)); |
---|
713 | twistedtubsElement->setAttributeNode(NewAttribute("phi", |
---|
714 | twistedtubs->GetDPhi()/degree)); |
---|
715 | twistedtubsElement->setAttributeNode(NewAttribute("aunit","deg")); |
---|
716 | twistedtubsElement->setAttributeNode(NewAttribute("lunit","mm")); |
---|
717 | solidsElement->appendChild(twistedtubsElement); |
---|
718 | } |
---|
719 | |
---|
720 | void G4GDMLWriteSolids:: |
---|
721 | ZplaneWrite(xercesc::DOMElement* element, const G4double& z, |
---|
722 | const G4double& rmin, const G4double& rmax) |
---|
723 | { |
---|
724 | xercesc::DOMElement* zplaneElement = NewElement("zplane"); |
---|
725 | zplaneElement->setAttributeNode(NewAttribute("z",z/mm)); |
---|
726 | zplaneElement->setAttributeNode(NewAttribute("rmin",rmin/mm)); |
---|
727 | zplaneElement->setAttributeNode(NewAttribute("rmax",rmax/mm)); |
---|
728 | element->appendChild(zplaneElement); |
---|
729 | } |
---|
730 | |
---|
731 | void G4GDMLWriteSolids::SolidsWrite(xercesc::DOMElement* gdmlElement) |
---|
732 | { |
---|
733 | G4cout << "G4GDML: Writing solids..." << G4endl; |
---|
734 | |
---|
735 | solidsElement = NewElement("solids"); |
---|
736 | gdmlElement->appendChild(solidsElement); |
---|
737 | |
---|
738 | solidList.clear(); |
---|
739 | } |
---|
740 | |
---|
741 | void G4GDMLWriteSolids::AddSolid(const G4VSolid* const solidPtr) |
---|
742 | { |
---|
743 | for (size_t i=0; i<solidList.size(); i++) // Check if solid is |
---|
744 | { // already in the list! |
---|
745 | if (solidList[i] == solidPtr) { return; } |
---|
746 | } |
---|
747 | |
---|
748 | solidList.push_back(solidPtr); |
---|
749 | |
---|
750 | if (const G4BooleanSolid* const booleanPtr |
---|
751 | = dynamic_cast<const G4BooleanSolid*>(solidPtr)) |
---|
752 | { BooleanWrite(solidsElement,booleanPtr); } else |
---|
753 | if (const G4Box* const boxPtr |
---|
754 | = dynamic_cast<const G4Box*>(solidPtr)) |
---|
755 | { BoxWrite(solidsElement,boxPtr); } else |
---|
756 | if (const G4Cons* const conePtr |
---|
757 | = dynamic_cast<const G4Cons*>(solidPtr)) |
---|
758 | { ConeWrite(solidsElement,conePtr); } else |
---|
759 | if (const G4EllipticalCone* const elconePtr |
---|
760 | = dynamic_cast<const G4EllipticalCone*>(solidPtr)) |
---|
761 | { ElconeWrite(solidsElement,elconePtr); } else |
---|
762 | if (const G4Ellipsoid* const ellipsoidPtr |
---|
763 | = dynamic_cast<const G4Ellipsoid*>(solidPtr)) |
---|
764 | { EllipsoidWrite(solidsElement,ellipsoidPtr); } else |
---|
765 | if (const G4EllipticalTube* const eltubePtr |
---|
766 | = dynamic_cast<const G4EllipticalTube*>(solidPtr)) |
---|
767 | { EltubeWrite(solidsElement,eltubePtr); } else |
---|
768 | if (const G4ExtrudedSolid* const xtruPtr |
---|
769 | = dynamic_cast<const G4ExtrudedSolid*>(solidPtr)) |
---|
770 | { XtruWrite(solidsElement,xtruPtr); } else |
---|
771 | if (const G4Hype* const hypePtr |
---|
772 | = dynamic_cast<const G4Hype*>(solidPtr)) |
---|
773 | { HypeWrite(solidsElement,hypePtr); } else |
---|
774 | if (const G4Orb* const orbPtr |
---|
775 | = dynamic_cast<const G4Orb*>(solidPtr)) |
---|
776 | { OrbWrite(solidsElement,orbPtr); } else |
---|
777 | if (const G4Para* const paraPtr |
---|
778 | = dynamic_cast<const G4Para*>(solidPtr)) |
---|
779 | { ParaWrite(solidsElement,paraPtr); } else |
---|
780 | if (const G4Paraboloid* const paraboloidPtr |
---|
781 | = dynamic_cast<const G4Paraboloid*>(solidPtr)) |
---|
782 | { ParaboloidWrite(solidsElement,paraboloidPtr); } else |
---|
783 | if (const G4Polycone* const polyconePtr |
---|
784 | = dynamic_cast<const G4Polycone*>(solidPtr)) |
---|
785 | { PolyconeWrite(solidsElement,polyconePtr); } else |
---|
786 | if (const G4Polyhedra* const polyhedraPtr |
---|
787 | = dynamic_cast<const G4Polyhedra*>(solidPtr)) |
---|
788 | { PolyhedraWrite(solidsElement,polyhedraPtr); } else |
---|
789 | if (const G4Sphere* const spherePtr |
---|
790 | = dynamic_cast<const G4Sphere*>(solidPtr)) |
---|
791 | { SphereWrite(solidsElement,spherePtr); } else |
---|
792 | if (const G4TessellatedSolid* const tessellatedPtr |
---|
793 | = dynamic_cast<const G4TessellatedSolid*>(solidPtr)) |
---|
794 | { TessellatedWrite(solidsElement,tessellatedPtr); } else |
---|
795 | if (const G4Tet* const tetPtr |
---|
796 | = dynamic_cast<const G4Tet*>(solidPtr)) |
---|
797 | { TetWrite(solidsElement,tetPtr); } else |
---|
798 | if (const G4Torus* const torusPtr |
---|
799 | = dynamic_cast<const G4Torus*>(solidPtr)) |
---|
800 | { TorusWrite(solidsElement,torusPtr); } else |
---|
801 | if (const G4Trap* const trapPtr |
---|
802 | = dynamic_cast<const G4Trap*>(solidPtr)) |
---|
803 | { TrapWrite(solidsElement,trapPtr); } else |
---|
804 | if (const G4Trd* const trdPtr |
---|
805 | = dynamic_cast<const G4Trd*>(solidPtr)) |
---|
806 | { TrdWrite(solidsElement,trdPtr); } else |
---|
807 | if (const G4Tubs* const tubePtr |
---|
808 | = dynamic_cast<const G4Tubs*>(solidPtr)) |
---|
809 | { TubeWrite(solidsElement,tubePtr); } else |
---|
810 | if (const G4TwistedBox* const twistedboxPtr |
---|
811 | = dynamic_cast<const G4TwistedBox*>(solidPtr)) |
---|
812 | { TwistedboxWrite(solidsElement,twistedboxPtr); } else |
---|
813 | if (const G4TwistedTrap* const twistedtrapPtr |
---|
814 | = dynamic_cast<const G4TwistedTrap*>(solidPtr)) |
---|
815 | { TwistedtrapWrite(solidsElement,twistedtrapPtr); } else |
---|
816 | if (const G4TwistedTrd* const twistedtrdPtr |
---|
817 | = dynamic_cast<const G4TwistedTrd*>(solidPtr)) |
---|
818 | { TwistedtrdWrite(solidsElement,twistedtrdPtr); } else |
---|
819 | if (const G4TwistedTubs* const twistedtubsPtr |
---|
820 | = dynamic_cast<const G4TwistedTubs*>(solidPtr)) |
---|
821 | { TwistedtubsWrite(solidsElement,twistedtubsPtr); } |
---|
822 | else |
---|
823 | { |
---|
824 | G4String error_msg = "Unknown solid: " + solidPtr->GetName(); |
---|
825 | G4Exception("G4GDMLWriteSolids::AddSolid()", "ReadError", |
---|
826 | FatalException, error_msg); |
---|
827 | } |
---|
828 | } |
---|