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